High Expression of Three-Gene Signature Improves Prediction of Relapse-Free Survival in Estrogen Receptor-Positive and Node-Positive Breast Tumors

Abstract

The objective of the present study was to validate prognostic gene signature for estrogen receptor alpha-positive (ERα+) and lymph node (+) breast cancer for improved selection of patients for adjuvant therapy In our previous study, we identified a group of seven genes (GATA3, NTN4, SLC7A8, ENPP1, MLPH, LAMB2, and PLAT) that show elevated messenger RNA (mRNA) expression levels in ERα (+) breast cancer patient samples. The prognostic values of these genes were evaluated using gene expression data from three public data sets of breast cancer patients (n = 395). Analysis of ERα (+) breast cancer cohort (n = 195) showed high expression of GATA3, NTN4, and MLPH genes significantly associated with longer relapse-free survival (RFS). Next cohort of ERα (+) and node (+) samples (n = 109) revealed high mRNA expression of GATA3, SLC7A8, and MLPH significantly associated with longer RFS. Multivariate analysis of combined three-gene signature for ERα (+) cohort, and ERα (+) and node (+) cohorts showed better hazard ratio than individual genes. The validated three-gene signature sets for ERα (+) cohort, and ERα (+) and node (+) cohort may have potential clinical utility since they demonstrated predictive and prognostic ability in three independent public data sets.

Keywords

breast cancer gene signature estrogen receptor lymph node positive relapse-free survival prognosis biomarkers

Introduction

Breast cancer is a heterogeneous disease,¹ making it an ideal disease to study using microarrays since different expression patterns can be identified within distinct tumor groups. An increased understanding of the pathogenesis of breast cancer is imperative in the pursuit of innovative therapies for treatment and/or prognosis of patients. Gene expression studies based on microarray have been used extensively by cancer researchers to profile cancer subsets, predict patients’ outcome, and identify genes of clinical relevance.^2–5

Breast cancer is no longer a single disease, but it is a heterogeneous disease consisting of different subtypes on the molecular and histopathological levels with different prognostic and therapeutic outcomes.^6–9 Gene expression profiling has classified breast cancer into five biologically distinct intrinsic subtypes: luminal A, luminal B, human epidermal growth factor receptor 2 (HER2+), basal-like, and normal-like.^6–9 The luminal A subtype is estrogen receptor alpha positive (ERα+), progesterone receptor positive (PR+), and HER2 (–); luminal B subtype is ERα+, PR+, and HER+, and luminal B is associated with a relatively worse outcome. Both HER2 (+) and basal-like (ER-, PR-, and HER-) breast cancers have poor outcomes. Parker et al.⁹ developed an efficient classifier, called PAM50, to distinguish these five intrinsic subtypes using the expression of 50 “classifier genes.” In a more recent study, a large breast cancer patient cohort (n ~ 2000) was clustered into 10 molecularly defined subgroups with apparently distinct biology and disease-specific survival characteristics.¹⁰ In addition, different breast cancer subtypes have different treatment responses.^11,12

An important part of the diagnostic workup of all the breast cancer patients is the determination of the ERα status of the tumor. Clinically, an ERα (+) status is associated with improved prognosis, lower risk of relapse, and better overall survival,¹³ and that are key aspects for making decisions for endocrine therapy with antiestrogens. A major problem in clinical oncology is to distinguish the patients who are likely to present a relapse of the disease from those with a favorable prognosis. In recent years, it has been realized that apart from ERα, other factors are also important in deciding the therapeutic strategies of the patient. These include histological markers such as grade, tumor size, lymph node involvement, PR, and HER2 receptor status. Each of these has modest positive predictive value (30%-60%).^14–17 Moreover, the current histological classifications of breast cancer do not fully represent the diverse clinical outcome of the disease. Recent approaches for patient management, which utilize histological markers in conjunction with online statistical algorithms such as “Nottingham Prognostic Index” and “Adjuvant! Online”, fail to predict the course of the disease in a significant number of breast cancer patients.^18,19 Women with node (+), ERα (+), and HER2 (–) often receive adjuvant treatment with chemotherapy and hormonal therapy. Nevertheless, few patients eventually experience a recurrence. Thus, new tools are needed to allow improved definition of a risk of recurrence. If it were possible to predict cancer recurrence following standard therapy, these patients could be targeted for alternative treatment strategies.

Recently, we published gene expression profile of breast tumors and identified seven genes (GATA3, NTN4, SLC7A8, MLPH, ENPP1, LAMB2, and PLAT) that showed high messenger RNA (mRNA) expression levels in ERα (+) compared to ERα (–) breast tumors.²⁰ In the present study, we have delineated the association between mRNA expression of these aforementioned seven genes with clinicopathological parameters such as PR, HER2, tumor grade, and lymph node status. To demonstrate that these genes have significantly better prediction for RFS than standard prognostic factors, we investigated the clinical utility of the aforementioned seven genes on 340 patient samples from three public data sets as validation cohort and used Kaplan-Meier survival curve using univariate and multivariate analysis.

Materials and Methods

Patients and breast cancer tissues

Human cells and specimens were obtained for previous research,²⁰ deidentified, and reused for this study. This research was therefore exempt from the requirement for ethics committee approval under US regulations §46.101(b)(4). The research was conducted in accordance with the standards in the Declaration of Helsinki. Breast tumor samples were obtained from patients undergoing surgery after informed written consent (Apollo Hospital). The excised tumor specimens were immediately preserved in RNA-later (Life Technologies) and stored at 4 °C until shipment. All tumor samples utilized in this study were invasive ductal carcinoma (Supplementary Table 1). All the histopathologi-cal information used in the analysis was directly documented from the original pathology reports. Grading, tumor type, ERα status, PR status, and HER2 status had been routinely recorded at the Apollo Hospital.

Public microarray data sets

Due to limited availability of clinical information in our data set, we used independent data sets to assess the predictive ability of our gene signatures, which will give us additional confidence in clinical validity. The data sets with Gene Expression Omnibus (GEO) accession numbers GSE2740, GSE1992, and GSE2607 generated using Agilent microarray platform have been previously described.^21–23 A total of 141 patient samples in GSE2740, 152 patient samples in GSE1992, and 102 patient samples in GSE2607 were used. A total of 395 patient samples with histopathological information are available in Supplementary Table 2. Patient samples with missing survival information were omitted; hence, a total of 340 samples were used for further analysis (Table 1). The clinical data were extracted from the gene expression data files downloaded from GEO. The ERα status, nodal status, survival data, and gene expression data were used for Kaplan-Meier survival curve and Cox multivariate analysis.

Table 1

The range of expression for seven genes and a median value used to stratify patient samples from public data sets into the groups of high expression and low expression.

	GATA3	NTN4	SLC7A8	MLPH	ENPP1	LAMB2	PLAT
Cohort I: 195 estrogen receptor positive (ER +) pat ient samples
mRNA expression range	1.7 to-4.1	2.7 to-3.9	1.5 to-4.0	2.2 to -2.1	2.5 to -2.5	1.1 to-1.2	2.8 to-3.2
Cut-off value	0.0	0.0	0.0	0.0	0.0	0.0	0.0
Number of high expression samples	95	96	97	91	97	96	95
Number of low expression samples	100	99	98	104	98	99	100
Cohort II: 109 ER (+) and node (+) patient samples
mRNA expression range	1.4 to–3.0	1.0 to-3.3	2.3 to-1.4	2.8 to -2.7	1.6 to-2.0	1.2 to-1.1	2.8 to-3.2
Cut-off value	0.0	0.0	0.0	0.0	0.0	0.0	0.0
Number of high expression samples	22	22	22	22	22	22	22
Number of low expression samples	87	87	87	87	87	87	87
	Cohort III: 67 ER (+) samples with 3-gene combined expression			Cohort IV: 43 ER (+) Node (+) samples with 3-gene combined expression
	GATA3/NTN4 /MLPH			GATA3/SLC7A8/MLPH
Cut-off value	0.0			0.0
Number of high expression samples	30			20
Number of low expression samples	37			23

Table 2

The relation between mRNA levels and clinicopathological parameters.

		ESTROGEN RECEPTOR (ER)		PROGESTERONE RECEPTOR (PR)		HER2		HISTOLOGICAL GRADE			NUMBER OF LYMPH NODES POSITIVE			PHENOTYPE
		POSITIVE (40)	NEGATIVE (36)	POSITIVE (42)	NEGATIVE (34)	POSITIVE (22)	NEGATIVE (54)	I (8)	II (36)	III (26)	0(25)	1–3 (20)	>3 (22)	ER(+)	LUMINAL A SUBTYPE	LUMINAL B SUBTYPE	TRIPLE NEGATIVE (BASAL)
GATA3 mRNA	Median	5.7	0.33	5.0	0.79	2.2	1.43	11.2	2.3	0.42	0.76	0.4	2.7	6.5	8.26	4.49	0.22
	P-value		<0.0001		<0.0001		0.785		0.069			0.219				0.0039
NTN4 mRNA	Median	3.8	0.52	3.9	0.72	2.2	1.28	1.6	1.3	0.9	1.1	0.66	2.5	1.3	8	2.5	0.41
	P-value		<0.0001		0.0003		0.93		0.305			0.348				0.0127
SLC7A8 mRNA	Median	15.3	1.9	15.75	3.25	6.48	9.48	13.6	10.04	4.245	9.46	5.52	9.48	9.82	17	12.2	0.67
	P-value		<0.0001		<0.0001		0.428		0.54			0.454				0.014
MLPH mRNA	Median	6.85	0.3	6.75	0.405	4.385	2.66	21.7	3.3	1.36	2.02 5	1.15	5.46	6.15	6.88	5.16	0.13
	P-value		<0.0001		<0.0001		0.271		0.013			0.212				0.02
ENPP1 mRNA	Median	4.55	0.79	3.81	0.95	1.92	2.72	2.8	3.25	1.18	1.95	0.95	4.1	4.55	5.93	2.35	0.38
	P-value		<0.0001		<0.0002		0.838		0.253			0.527				0.0003
LAMB2 mRNA	Median	5.87	0.94	6.1	1.05	2.81	3.16	7.53	3	1.01	3	1.09	3.69	2.74	7.63	3.68	0.7
	P-value		0.0006		<0.0001		0.628		0.698			0.179				0.0014
PLAT mRNA	Median	1.62	0.59	1.62	0.6	1.28	1.05	2.41	1.7	0.67	1.15	0.78	2.26	0.69	2.96	1.04	0.21
	P-value		0.03		0.027		0.975		0.544			0.101				0.6759

Note: Significant P-values are in italics.

Statistical analysis

To visualize gene expression values using heat maps, the values for each probe were centered by subtracting the mean expression value across patients. No gene-specific scaling (standardization) was performed, and, thus, information about the relative signal strength between probes was retained. The color tone in the heat maps was calibrated so that saturated red and saturated green were reached at values equal to 3.5-fold the standard deviation of the expression values of the entire matrix. Red and green reflect high- and low-expression levels (log 2-transformed scale), respectively.

The gene expression data from 340 patient samples were dichotomized according to the median (cutoffvalue 0.0) of the complete cohort. The expression data higher than the median were grouped into the “high-expression” group, and the expression values lesser than the median were grouped into the “low-expression” group (Table 1). Mann-Whitney t-test was used to evaluate the difference between mRNA levels of genes and clinicopathological parameters. Survival distributions were estimated using the Kaplan-Meier method (univariate), and the significance of differences between survival rates was ascertained by the log-rank test using the GraphPad software. Candidate prognostic factors for RFS with a 0.05 significance level in univariate analysis were entered in a multivariate Cox model.²⁴ Multivariate analysis was evaluated by step-wise forward Cox's regression analysis. The Cox proportional hazard model was used to calculate the hazard ratio (HR) and 95% confidence interval (CI) in the analysis of relapse-free survival (RFS). RFS was measured from the date of diagnosis to relapse or censored at the last follow-up. A P-value less than 0.05 were considered statistically significant.

Results

Relationship between mRNA levels of seven genes and clinicopathological parameters

Initially, we sought to corroborate the findings of our recent study.²⁰ We correlated the mRNA expression of these seven genes with clinicopathological parameters. In our previous study,²⁰ we had performed reverse transcription quantitative polymerase chain reaction (RT-qPCR) analyses of 76 tumor specimens and showed that the expression of aforementioned seven genes was significantly associated with ERα (+) tumors (P < 0.01). Herein, we utilized the same RT-qPCR data and classified it based on PR, HER2, tumor grade, and lymph node status (Supplementary Table 1). We observed that increased expression of the aforementioned seven genes was significantly associated with PR (+) breast tumors (P < 0.05). In contrast, no such association was found between mRNA expressions of these genes with HER2 receptor status and lymph node status (Table 2). Interestingly, six out of the seven genes did not show any association with regard to tumor grade. The only mRNA expression level of MLPH (P = 0.013) was significantly higher in grade I than grade III tumor. Given that the increased expression of each of the seven genes was associated with not only ERα (+) status but also PR (+) status, we ascertained if there was any correlation between the expression of these genes with ERα (+) and PR (+) breast tumors. We observed that this latter group of patients expresses statistically significant higher mRNA levels of GATA3, NTN4, SLC7A8, MLPH, ENPP1, LAMB2, and PLAT compared to triple negative [ERα (–)/PR (–)/HER2 (–)] breast tumors (Table 2). Interestingly among endocrine-treated patients, the previous report²⁵ also showed that the presence of both ERα and PR was a stronger marker for the benefit of adjuvant endocrine therapy than ER alone. Accordingly, high expression of our gene signature could be used as a prediction marker for endocrine therapy.

Meta-analysis of seven genes deregulated in ERα (+) breast tumors

Next, we sought to correlate the mRNA expression of these seven genes with long-term survival data. Since we did not have an adequate number of clinical samples, we utilized published data sets to validate our predictive gene set. Further use of independent data sets to assess the predictive ability of our gene signatures will give us additional confidence in clinical validity. Accordingly, gene expression data were obtained from three independent public data sets (n = 395).^21–23 The 340 patient samples from public data sets were considered for further analysis based on the available survival information (Table 1; Supplementary Table 2). Out of 340 samples, there were 195 ERα (+) and 145 ERα (–) samples. The seven genes dysregulated in ERα (+) breast tumors (GATA3, NTN4, SLC7A8, MLPH, ENPP1, LAMB2, and PLAT) were utilized to cluster all the tumors and, subsequently, the data were visualized in the form of heat map (Fig. 1). Out of 195 ERα (+) and 145 ERα (–) samples, our seven genes could classify 181 ERα (+) samples and 135 ERα (–) samples correctly with 93% positive prediction. These proved prediction value of mRNA levels of these seven genes to classify ERα (+) and ERα (–) samples.

Figure 1

Dendrogram of 340 breast cancer samples from public data sets. Unsupervised, hierarchical, uncentered Pearson distance (co-relation) clustering was performed to classify the seven genes into homogeneous clusters.

High mRNA expressions of GATA3, NTN4, and MLPH are associated with longer RFS in ERα (+) breast tumors

The gene expression values from the public data sets were dichotomized according to the median of the complete cohort, and expression data higher than the median were grouped into the high-expression group, and the expression values lesser than the median were grouped into the low-expression group (Table 1). Univariate analysis on ERα (+) test data sets (n = 195; Supplementary Table 3) revealed that high mRNA expression levels of GATA3 (P= 0.0003), NTN4 (P = 0.0011), SLC7A8 (P = 0.012), and MLPH (P = 0.0054) were significantly associated with longer RFS. Cox multivariate analysis revealed that GATA3 (P = 0.0167), NTN4 (P = 0.0044), and MLPH (P = 0.0321) were independent prognostic markers and significantly associated with RFS (Fig. 2; Table 3).

Figure 2

Kaplan-Meier survival curve using high and low mRNA expression among ERα (+) breast tumors from public data sets (n = 195). Univariate analysis (n = 195) revealed that high mRNA expression levels of (A) GATA3 (P = 0.0003), (B) NTN4 (P = 0.0011), (C) SLC7A8 (P = 0.012), (D) MLPH (P < 0.0054), (E) ENPP1 (P = 0.3191), (F) LAMB2 (P = 0.3664), and (G) PLAT(P = 0.119) were significantly associated with longer relapse-free survival. A significance of difference between survival rates was ascertained by the log-rank test and P-value less than 0.05 was considered statistically significant.

Table 3

The univariate and multivariate analysis in relation to RFS among 195 ERα (+) breast cancer patient samples from public data sets.

GENES	UNIVARIATE			MULTIVARIATE
	HAZARD RATIO	95% CI	P-VALUE	HAZARD RATIO	95% CI	P-VALUE	LIKELIHOOD RATIO
GATA3	0.248	0.1360 to 0.5570	0.0003	0.338	0.139 to 0.822	0.0167	271.19
NTN4	0.274	0.1542 to 0.6257	0.0011	0.275	0.113 to	0.0044	261.97
SLC7A8	0.392	0.1998 to 0.8195	0.012	N/A	N/A	NS	N/A
MLPH	0.367	0.1796 to 0.7425	0.0054	0.418	0.188 to	0.0321	256.55
ENPP1	1.417	0.7054 to 2.918	0.3191	N/A	N/A	N/A	N/A
LAMB2	0.725	0.3599 to 1.458	0.3664	N/A	N/A	N/A	N/A
PLAT	0.570	0.2805 to 1.156	0.119	N/A	N/A	N/A	N/A

Note: Significant P-values are in italics.

Abbreviations: NS, non-significant; N/A, not available; CI, confidence interval; ER, estrogen receptor; RFS, relapse-free survival.

High mRNA expressions of GATA3, SLC7A8, and MLPH are associated with longer RFS in ERα (+) and node (+) breast tumors

Having studied the prognostic significance of the seven dysregulated genes in ERα (+) patients, we next sought to perform our analysis with the clinically important issue of the metastatic spread of the tumor. The determination of the extent of lymph node involvement in primary breast cancer is the single most important risk factor in disease outcome. Accordingly, we next investigated the correlation of these seven dysregulated genes with ERα (+) and node (+) cohort (n = 109; Supplementary Table 4). In univariate analysis, high mRNA expression of GATA3 (P = 0.003), SLC7A8 (P = 0.0045), MLPH (P = 0.0021), and PLAT (P = 0.0311) showed significantly longer RFS. In Cox multivariate analysis, high mRNA expression of GATA3 (P= 0.0009), SLC7A8 (P = 0.0431), and MLPH (P = 0.0170) showed significantly longer RFS. In contrast, high mRNA expression of ENPP1 (P= 0.0001) was significantly associated with the worst RFS (Table 4). Interestingly, in a subset analysis of ERα (+) and node (–-) cohort of patients (n = 84), none of the genes persisted in multivariate analysis (Table 5).

Table 4

Univariate and multivariate analysis in relation to RFS among 109 ERα (+) and node (+) breast cancer patient samples from public data sets.

GENES	UNIVARIATE			MULTIVARIATE
	HAZARD RATIO	95% CI	P VALUE	HAZARD RATIO	95% CI	P-VALUE	LIKELIHOOD RATIO
GATA3	12861	0.1490 to	0.003	0.148	0.047 to	0.0009	191.67
NTN4	15296	0.2480 to 1.134	0.1018	N/A	N/A	N/A	N/A
SLC7A8	13033	0.1522 to	0.0045	0.372	0.140 to	0.0431	180.58
MLPH	12785	0.1392 to	0.0021	0.293	0.106 to	0.0170	210.32
ENPP1	2.27	1.099 to 5.103	0.0278	5.853	2.318 to	0.0001	200.47
LAMB2	1584	0.2566 to 1.209	0.1389	N/A	N/A	N/A	N/A
PLAT	14233	0.2055 to	0.0311	N/A	N/A	NS	N/A

Note: Significant P-values are in italics.

Abbreviations: NS, non-significant; N/A, not available; CI, confidence interval; ER, estrogen receptor; RFS, relapse-free survival.

Table 5

Univariate and multivariate analysis in relation to RFS among 84 ERα (+) and node (–) breast cancer patient samples from public data sets.

GENES	UNIVARIATE			MULTIVARIATE
	HAZARD RATIO	95% CI	P- VALUE	HAZARD RATIO	95% CI	P-VALUE	LIKELIHOOD RATIO
GATA3	0.451	0.06059 to 3.606	0.4656	N/A
NTN4	0.2637	0.02665 to 1.949	0.1769	N/A
SLC7A8	0	0.007922 to 0.5021	0.009	NS
MLPH	N/A	N/A	0.172	N/A
ENPP1	0.2637	0.02665 to 1.949	0.1769	N/A
LAMB2	0	0.01232 to 0.7328	0.0239	NS
PLAT	4.5	0.7140 to 78.62	0.0931	N/A

Note: Significant P-values are in italics.

Abbreviations: NS, non-significant; N/A, not available; CI, confidence interval; RFS, relapse-free survival; ER, estrogen receptor.

Elevated expression of three-gene signature improves RFS

In our previous analyses, we had ascertained the predictive power of “individual” genes. We next sought to determine the “combined” predictive power of the three-gene signature. Accordingly, patients who expressed high mRNA levels of “all” three genes were studied against those who expressed low mRNA levels of all three genes with the cutoff value of zero. Initially, we looked into the 195 samples of ERα (+) cohort. In multivariate analysis, the three genes GATA3, NTN4, and MLPH showed independent prognostic value and were significantly associated with RFS. Hence, to determine the combined prediction value of all three genes, the 67 patient samples expressing similar expression of all three genes were selected (Supplementary Table 5). Patients having the high mRNA expression of all three genes GATA3, NTN4, and MLPH (n = 30) were compared with patients having the low mRNA expression of all three genes (n = 37). Multivariate analyses revealed that the three-gene signature showed better HR (0.056) and longer RFS than individual gene signature (Table 6; Fig. 3). Next, we looked into the 109 samples of ERα (+) and node (+) cohort. To determine the combined prediction value of all three genes GATA3, SLC7A8, and MLPH in this cohort, the 43 patient samples expressing similar expression of all three genes were selected (Supplementary Table 6). Patient samples showing high mRNA expression of all three genes GATA3, SLC7A8, and MLPH (n = 20) were compared with patients having low mRNA expression (n = 23) and observed that the combined three-gene signature showed better HR (0.057) than individual genes (Table 6; Fig. 3).

Table 6

Univariate and multivariate analysis of three-gene signature in relation to RFS among ERα cohort (π = 67), and ERα (+) and node (+) samples (n = 43) from public data sets.

SAMPLE CLASSIFICATION	UNIVARIATE				MULTIVARIATE
	HIGH OR LOW EXPRESSION OF ALL3-GENES	HAZARD RATIO	95% CI	P-VALUE	HAZARD RATIO	95% CI	P-VALUE	LIKELIHOOD RATIO
ERα (+) cohort	GATA3/NTN4/MLPH	0.1058	0.052 to 0.376	<0.0001	0.056	0.0074 to 0.4268	0.0053	115.19
ERα (+) and node (+)	GATA3/SLC7A8/MLPH	0.0581	0.052 to 0.367	<0.0001	0.057	0.0076 to 0.4320	0.0055	99.93

Note: Significant P-values are in italics.

Abbreviations: CI, confidence interval; ER, estrogen receptor.

Figure 3

(A) Kaplan-Meier survival curve using high mRNA expression of all three genes GATA3, NTN4, and MLPH (n = 30) was compared with patients having the low mRNA expression of all three genes (n = 37) among ERα (+) breast cancer samples from public data sets (n = 67); (B) Kaplan-Meier survival curve using high mRNA expression of all three genes GATA3, SLC7A8, and MLPH (n = 20) was compared with patients having the low mRNA expression (n = 23) among ERα (+) and node (+) breast tumors from public data sets (n = 43). A significance of difference between survival rates was ascertained by the log-rank test and P-value less than 0.05 was considered statistically significant.

Discussion

The current method of determination of ERα status by immunohistochemistry under clinical setup provides information about the expression pattern of ERα with no information on possibly disabled downstream ER pathway.^14,17 Thus, it is plausible that the status of the ER pathway is also clinically relevant and may explain variable response to endocrine therapy in ERα (+) patients. Hence, measurements of gene expression profiles that reflect the activity of ER pathway could provide an important insight in understanding the behavior of breast cancers. We reported the expression pattern for seven genes that were regulated by ERα and demonstrated the differential expression in various phenotypes of breast cancer pathology. We observed the highest expression of five genes in ERα (+) and PR (+) cohort of patients (luminal A subtype). Interestingly, among endocrine-treated patients, previous report²⁵ also showed that the presence of both ERα and PR was a stronger marker for the benefit of adjuvant endocrine therapy than ER alone.

Due to limited availability of clinical information in our data set, we used independent data sets^21–23 to assess the predictive ability of our gene signatures, which will give us additional confidence in clinical validity. Meta-analysis of seven genes showed that a gene expression is reliable and robust to determine ER expression in three independent data sets comprising 340 tumor samples. These analyses strengthen our recent findings and suggest that a seven-gene signature can stratify/classify ERα (+) and ERα (–) tumors in an independent data set. Of the seven genes analyzed in the ERα (+) tumors, the high mRNA expression of GATA3, NTN4, and MLPH emerged as an independent prognostic factor in multivariate analysis. In ERα (+) and node (+) tumors, the high mRNA levels of GATA3, SLC7A8, and MLPH emerged as an independent prognostic factor in multivariate analysis. In our findings, we did not find any literature describing the importance of SLC7A8 and MLPH in prognosis and overall disease-free survival of breast cancer patients.

Classic parameters, such as ER, PR, HER2 status, number of lymph nodes positive, and tumor size, have been integrated into software applications such as Adjuvant! Online¹⁹ to help doctors in calculating a risk of relapse and benefit from adjuvant therapy. However, uncertainty remains in many cases even with the use of this software. The well-established prognostic and/or therapeutic breast cancer markers are hormone receptors (ER and PR),²⁶ HER2,²⁷ Ki-67 antigen,²⁸ tumor protein p53,²⁹ carbohydrate 15-3 and carcinoembryonic antigens (CA 15-3 and CEA),^30,31 and breast cancer susceptibility genes (BRCA1 and BRCA2).³² Gene signatures can complement classic prognostic factors to obtain more accurate prognostic information. The 70-gene signature (MammaPrint; Agendia) and the 21-gene signature (OncoType; Genomic Health) are being used in selected patients with early ERα (+) disease to identify those women who will be cured even if they do not receive adjuvant chemotherapy.^4,33 These signatures have been extensively studied and are widely used in Europe and USA.^34–36 The National Cancer Comprehensive Network guidelines indicate that the 21-gene signature can be considered in women with tumors >0.5 cm, HER2-negative disease, and node-negative disease.³⁷ Limitations of the 21-gene and 70-gene signatures are intended to be used by women with node-negative breast cancer diagnosis. Many other gene signatures have been developed and have undergone validation. One of them is the breast cancer gene expression ratio test, which only measures the ratio of HOXB13 to IL17BR^38,39 A high mRNA expression ratio was associated with a high risk of recurrence in tamoxifen-treated patients. Recently, the accuracy of this test could be improved by including proliferation-associated genes of the molecular grade index,⁴⁰ which is an RT-qPCR assay consisting of five genes that are able to identify a subgroup of ER (+) patients with a worse outcome despite endocrine therapy. The Rotterdam 76-gene signature was created on the basis of predicting the development of metastatic disease within 5 years using an unselected patient cohort regarding age, tumor size, grade, and hormone receptor status.^41,42

The 5-year survival for patients with the node-negative disease is 82.8% compared with 73% for 1-3 positive nodes, 45.7% for 4-12 positive nodes, and 28.4% for ≥13 positive nodes.⁴³ These data demonstrate that the risk of recurrence is significant enough with lymph node-positive disease to warrant adjuvant systemic therapy since, generally, a future risk of distant recurrence of 20% or greater is regarded significant enough to consider the risks of therapy. Hence, it is important to stratify ERα (+) and node (+) patients into a low- and high-risk group for RFS. In this study's sample population, out of 109 ER (+) and node (+) samples, 67 samples showed either high or low expression of all GATA3/NTN4/MLPH genes, whereas 43 samples showed either high or low expression of one or two GATA3/SLC7A8/MLPH genes. So the limitation of this gene signature is that not all patient samples will exhibit either low or high gene expression of all three genes and has to be excluded from the prediction.

The GATA3 transcription factor has been studied intensively in the immune system but has most recently been shown to be important in the context of breast cancer and the ERα pathway. Recently, it was demonstrated that GATA3 is required for estradiol stimulation of cell-cycle progression of breast cancer cells.⁴⁴ Meta-analysis of four breast cancer microarray data sets revealed GATA3 as a promising novel prognostic biomarker in breast cancer with HR of 0.12 and P-value of 0.05.⁴⁵ As reported in the present study, combined three-gene signature showed better HR and improved RFS. The netrins are a family of secreted proteins that are highly conserved through evolution. Three netrins, netrinl (NTN1), netrin 3 (NTN3), and netrin 4 (NTN4), have been identified. Both the netrins and their receptors are widely expressed and have been implicated in a wide range of development processes including axon guidance, angiogenesis, and mammary gland development.^46,47 In an earlier study, NTN4expression was demonstrated to be an independent predictor of improved outcome with a HR of 0.2544 and P-value of 0.015.⁴⁸ Our study reports three-gene signature with better HR and improved RFS in ERα (+) tumors.

In summary, we show that the high expression of estrogen-responsive three-gene signature GATA3/NTN4/ MLPH is associated with good prognosis in ERα (+) patients and GATA3/SLC7A8/MLPH is associated with good prognosis in ERα (+) and node (+) patients and are promising novel prognostic biomarkers in breast cancer. Our gene sets may have potential clinical utility since they demonstrated predictive ability in three independent public data sets. Although our results are promising, it needs to be validated in large cohort breast cancer patients. Such signatures will also be the starting point for functional profiling of genes and proteins to understand the biological processes associated with disease formation and progression.

Author Contributions

Conceived and designed the experiments: AT, MP, AB, HR. Analyzed the data: AT. Wrote the first draft of the manuscript: AT. Contributed to the revising manuscript critically for important intellectual content: AT, AB, BM, MP. All authors made an intellectual input and contributed to writing the paper. All authors reviewed and approved of the final manuscript.

Footnotes

Acknowledgment

Debarshi Chakrabarti provided scientific input, discussions and deliberations during the time of the study.

Supplementary Materials

References

Russo

, Hu

Y.F.

, Yang

, Russo

I.H.

Developmental, cellular, and molecular basis of human breast cancer. JNatl Cancer Inst Monogr. 2000; (27): 17–37.

Lonning

P.E.

, Sorlie

, Perou

C.M.

, Brown

P.O.

, Botstein

, Borresen-Dale

A.L.

Microarrays in primary breast cancer - lessons from chemotherapy studies. Endocr Relat Cancer. 2001; 8(3): 259–63.

Sotiriou

, Neo

S.Y.

, McShane

L.M.

. Breast cancer classification and prognosis based on gene expression profiles from a population-based study. Proc Natl Acad Sci USA. 2003; 100(18): 10393–8.

van de Vijver

M.J.

, He

Y.D.

, van't Veer

L.J.

. A gene-expression signature as a predictor of survival in breast cancer. N End J Med. 2002: 347(25): 1999–2009.

Huang

, Cheng

S.H.

, Dressman

. Gene expression predictors of breast cancer outcomes. Lancet. 2003; 361(9369): 1590–6.

Perou

C.M.

, Sørlie

, Eisen

M.B.

. Molecular portraits of human breast tumours. Nature. 2000; 406(6797): 747–52.

Sørlie

, Perou

C.M.

, Tibshirani

. Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci USA. 2001; 98(19): 10869–74.

Sorlie

, Tibshirani

, Parker

. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA. 2003; 100(14): 8418–23.

Parker

J.S.

, Mullins

, Cheang

M.C.

. Supervised risk predictor of breast cancer based on intrinsic subtypes. J Clin Oncol. 2009; 27(8): 1160–7.

10.

Curtis

, Shah

S.P.

, Chin

S.F.

.; METABRIC Group. The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012; 486(7403): 346–52.

11.

Rouzier

, Perou

C.M.

, Symmans

W.F.

. Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Clin Cancer Res. 2005; 11(16): 5678–85.

12.

Rody

, Karn

, Solbach

. The erbB2+ cluster of the intrinsic gene set predicts tumor response of breast cancer patients receiving neoadjuvant chemotherapy with docetaxel, doxorubicin and cyclophosphamide within the GEPAR-TRIO trial. Breast. 2007; 16(3): 235–40.

13.

Ring

B.Z.

, Seitz

R.S.

, Beck

. Novel prognostic immunohistochemical biomarker panel for estrogen receptor-positive breast cancer. J Clin Oncol. 2006; 24(19): 3039–47.

14.

Bonneterre

, Thürlimann

, Robertson

J.F.

. Anastrozole versus tamoxifen as first-line therapy for advanced breast cancer in 668 postmenopausal women: results of the Tamoxifen or Arimidex Randomized Group Efficacy and Tolerability Study. J Clin Oncol. 2000; 18(22): 3748–57.

15.

Colozza

, Azambuja

, Cardoso

, Sotiriou

, Larsimont

, Piccart

M.J.

Proliferative markers as prognostic and predictive tools in early breast cancer: where are we now?

Ann Oncol. 2005; 16(11): 1723–39.

16.

Hayes

D.F.

Prognostic and predictive factors revisited. Breast. 2005; 14(6): 493–9.

17.

Mouridsen

, Gershanovich

, Sun

. Superior efficacy of letrozole versus tamoxifen as first-line therapy for postmenopausal women with advanced breast cancer: results of a phase III study of the International Letrozole Breast Cancer Group. J Clin Oncol. 2001; 19(10): 2596–606.

18.

D'Eredita

, Giardina

, Martellotta

, Natale

, Ferrarese

Prognostic factors in breast cancer: the predictive value of the Nottingham Prognostic Index in patients with a long-term follow-up that were treated in a single institution. Eur J Cancer. 2001; 37(5): 591–6.

19.

Olivotto

I.A.

, Bajdik

C.D.

, Ravdin

P.M.

. Population-based validation of the prognostic model ADJUVANT! for early breast cancer./ Clin Oncol. 2005; 23(12): 2716–25.

20.

Thakkar

A.D.

, Raj

, Chakrabarti

. Identification of gene expression signature in estrogen receptor positive breast carcinoma. Biomark Cancer. 2010; 2: 1–15.

21.

, Fan

, Oh

D.S.

. The molecular portraits of breast tumors are conserved across microarray platforms. BMC Genomics. 2006; 7: 96.

22.

D.S.

, Troester

M.A.

, Usary

. Estrogen-regulated genes predict survival in hormone receptor-positive breast cancers. J Clin Oncol. 2006; 24(11): 1656–64.

23.

Perreard

, Fan

, Quackenbush

J.F.

. Classification and risk stratification of invasive breast carcinomas using a real-time quantitative RT-PCR assay. Breast Cancer Res. 2006; 8(2): R23.

24.

Cox

Multinomial regression models based on continuation ratios. Stat Med. 1988; 7(3): 435–41.

25.

Bardou

V.J.

, Arpino

, Elledge

R.M.

, Osborne

C.K.

, Clark

G.M.

Progesterone receptor status significantly improves outcome prediction over estrogen receptor status alone for adjuvant endocrine therapy in two large breast cancer databases. J Clin Oncol. 2003; 21(10): 1973–9.

26.

Althuis

M.D.

, Fergenbaum

J.H.

, Garcia-Closas

, Brinton

L.A.

, Madigan

M.P.

, Sherman

M.E.

Etiology of hormone receptor-defined breast cancer: a systematic review of the literature. Cancer Epidemiol Biomarkers Prev. 2004; 13(10): 1558–68.

27.

Slamon

D.J.

, Clark

G.M.

, Wong

S.G.

, Levin

W.J.

, Ullrich

, McGuire

W.L.

Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987; 235(4785): 177–82.

28.

Vielh

, Chevillard

, Mosseri

, Donatini

, Magdelenat

Ki67 index and S-phase fraction in human breast carcinomas. Comparison and correlations with prognostic factors. Am J Clin Pathol. 1990; 94(6): 681–6.

29.

Dumay

, Feugeas

J.P.

, Wittmer

. Distinct tumor protein p53 mutants in breast cancer subgroups. Int J Cancer. 2013; 132(5): 1227–31.

30.

Ebeling

F.G.

, Stieber

, Untch

. Serum CEA and CA 15-3 as prognostic factors in primary breast cancer. Br J Cancer. 2002; 86(8): 1217–22.

31.

Vizcarra

, Lluch

, Cibrián

. Value of CA 15.3 in breast cancer and comparison with CEA and TPA: a study of specificity in disease-free follow-up patients and sensitivity in patients at diagnosis of the first metastasis. Breast Cancer Res Treat. 1996; 37(3): 209–16.

32.

Senst

, Llacuachaqui

, Lubinski

.; Hereditary Breast Cancer Study Group. Parental origin of mutation and the risk of breast cancer in a prospective study of women with a BRCA1 or BRCA2 mutation. Clin Genet. 2013; 84(1): 43–6.

33.

Paik

, Shak

, Tang

. A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004; 351(27): 2817–26.

34.

Albain

K.S.

, Barlow

W.E.

, Shak

.; Breast Cancer Intergroup of North America. Prognostic and predictive value of the 21-gene recurrence score assay in postmenopausal women with node-positive, oestrogen-receptor-positive breast cancer on chemotherapy: a retrospective analysis of a randomised trial. Lancet Oncol. 2010; 11(1): 55–65.

35.

Mook

, Schmidt

M.K.

, Viale

.; TRANSBIG Consortium. The 70-gene prognosis-signature predicts disease outcome in breast cancer patients with 1–3 positive lymph nodes in an independent validation study. Breast Cancer Res Treat. 2009; 116(2): 295–302.

36.

Mittempergher

, de Ronde

J.J.

, Nieuwland

. Gene expression profiles from formalin fixed paraffin embedded breast cancer tissue are largely comparable to fresh frozen matched tissue. PLoS One. 2011; 6(2): e17163.

37.

Theriault

R.L.

, Carlson

R.W.

, Allred

.; National Comprehensive Cancer Network. Breast cancer, version 3.2013: featured updates to the NCCN guidelines. J Natl Compr Canc Netw. 2013; 11(7): 753–60. quiz 761.

38.

X.J.

, Hilsenbeck

S.G.

, Wang

. The HOXB13: IL17BR expression index is a prognostic factor in early-stage breast cancer./ Clin Oncol. 2006; 24(28): 4611–9.

39.

X.J.

, Salunga

, Dahiya

. A five-gene molecular grade index and HOXB13: IL17BR are complementary prognostic factors in early stage breast cancer. Clin Cancer Res. 2008; 14(9): 2601–8.

40.

Wang

, Dahiya

, Provencher

. The prognostic biomarkers HOXB13, IL17BR, and CHDH are regulated by estrogen in breast cancer. Clin Cancer Res. 2007; 13(21): 6327–34.

41.

Wang

, Klijn

J.G.

, Zhang

. Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005; 365(9460): 671–9.

42.

Loi

, Haibe-Kains

, Desmedt

. Definition of clinically distinct molecular subtypes in estrogen receptor-positive breast carcinomas through genomic grade. J Clin Oncol. 2007; 25(10): 1239–46.

43.

Fisher

, Bauer

, Wickerham

D.L.

. Relation of number of positive axillary nodes to the prognosis of patients with primary breast cancer. An NSABP update. Cancer. 1983; 52(9): 1551–7.

44.

Eeckhoute

, Keeton

E.K.

, Lupien

, Krum

S.A.

, Carroll

J.S.

, Brown

Positive cross-regulatory loop ties GATA-3 to estrogen receptor alpha expression in breast cancer. Cancer Res. 2007; 67(13): 6477–83.

45.

Mehra

, Varambally

, Ding

. Identification of GATA3 as a breast cancer prognostic marker by global gene expression meta-analysis. Cancer Res. 2005; 65(24): 11259–64.

46.

Hinck

The versatile roles of “axon guidance” cues in tissue morphogenesis. Dev Cell. 2004; 7(6): 783–93.

47.

Yamamura

, Miyoshi

, Tamaki

. mRNA expression level of estrogen-inducible gene, alpha 1-antichymotrypsin, is a predictor of early tumor recurrence in patients with invasive breast cancers. Cancer Sci. 2004; 95(11): 887–92.

48.

Esseghir

, Kennedy

, Seedhar

. Identification of NTN4, TRA1, and STC2 as prognostic markers in breast cancer in a screen for signal sequence encoding proteins. Clin Cancer Res. 2007; 13(11): 3164–73.