- Open Access
Meta-analysis of human cancer microarrays reveals GATA3 is integral to the estrogen receptor alpha pathway
Molecular Cancervolume 7, Article number: 49 (2008)
The transcription factor GATA3 has recently been shown to be necessary for mammary gland morphogenesis and luminal cell differentiation. There is also an increasing body of data linking GATA3 to the estrogen receptor α (ERα) pathway. Among these it was shown that GATA3 associates with the promoter of the ERα gene and ERα can reciprocally associate with the GATA3 gene. GATA3 has also been directly implicated in a differentiated phenotype in mouse models of mammary tumourigenesis. The purpose of our study was to compare coexpressed genes, by meta-analysis, of GATA3 and relate these to a similar analysis for ERα to determine the depth of overlap.
We have used a newly described method of meta-analysis of multiple cancer studies within the Oncomine database, focusing here predominantly upon breast cancer studies. We demonstrate that ERα and GATA3 reciprocally have the highest overlap with one another. Furthermore, we show that when both coexpression meta-analysis lists for ERα and GATA3 are compared there is a significant overlap between both and, like ERα, GATA3 coexpresses with ERα pathway partners such as pS2 (TFF1), TFF3, FOXA1, BCL2, ERBB4, XBP1, NRIP1, IL6ST, keratin 18(KRT18) and cyclin D1 (CCND1). Moreover, as these data are derived from human tumour samples this adds credence to previous cell-culture or murine based studies.
GATA3 is hypothesized to be integral to the ERα pathway given the following: (1) The large overlap of coexpressed genes as seen by meta-analysis, between GATA3 and ERα, (2) The highest coexpressing gene for GATA3 was ERα and vice-versa, (3) GATA3, like ERα, coexpresses with many well-known ERα pathway partners such as pS2.
While GATA3 has most intensively been studied in the immune system  GATA3 is also expressed in other biological environments such as the human endometrium epithelial cells, where levels are regulated in a cyclic manner . GATA3 levels are also considered a good prognostic biomarker in breast tumours. Specifically, in the luminal A subtype of breast cancer GATA3 has both a favorable prognostic outcome, and the highest ERα and GATA3 expression of all breast tumours . Consistent with this, basal-like tumours have the lowest GATA3 expression and the worst prognosis. GATA3 has also been shown in murine models to be essential to the development and maintenance of mammary luminal cells [4, 5]. There is also tentative data showing that different polymorphisms of the GATA3 gene may associate with differential susceptibility to breast cancer .
GATA3 levels have previously been correlated with expression of ERα  and both were shown to reciprocally regulate one another at the transcriptional level in a cell-culture based system in a cross-regulatory loop . Furthermore, in a meta-analysis of ERα 10 genes were proposed as classifier of ERα positive breast tumours, listing GATA3 as one of these . A study has also compared microarray experiments between estradiol-induced genes from MCF-7 cells, and transfected GATA3-induced genes from 293T cells to assess common upregulated genes .
In an elegant series of experiments utilizing MMTV-PyMT (polyoma middle T antigen) mice it was first shown that GATA3 expression was downregulated with the transition from adenoma to carcinoma in mammary tumours, and the expression was lost in lung metastases. Infection of the MMTV-PyMT carcinomas with GATA3 upregulated markers of differentiation and resulted in a dramatic 27-fold reduction in lung metasases . Further crossing of these mice with an inducible cre-WAP (whey acidic protein – specific to luminal mammary epithelial cells) driven knockout of GATA3, resulted in loss of markers of terminal differentiation, detachment from the basal membrane and apoptosis. This is consistent with the requirement of GATA3 in differentiated tumours.
As described in a recent study known pathway partners have been shown to yield a similar 'meta-analysis coexpression signature' i.e. having a significant overlap of coexpressed genes can link proteins to the same pathways . Thus performing independent meta-analyses for ERα and GATA3 (putative pathway partners), then comparing the results for overlapping genes would yield a highly significant number of genes if these transcription factors were in the same pathway. We report here not only that these meta-analyses have a high degree of overlap, but that genes identified are consistent with previous reports of the ERα pathway regulation. Additionally we show this correlation with previously identified ERα target genes by combining our meta-analysis data with both RT-PCR and genome-wide location analysis from other studies. These data not only confirm GATA3 as being a key player in the ERα pathway, but also give fresh insights into the pathway itself.
The following procedure was undertaken for independent meta-analyses of GATA3 or ERα: a co-expression gene search was performed within Oncomine . Twenty-one studies were chosen for analysis, most of which were breast cancer studies. The top 400 coexpressed genes were extracted and filtered to give one representative gene per study (removing duplicates and ESTs). These filtered genelists were then compared for repeating coexpressed genes over multiple studies. The frequency cutoff was 3 studies (14% of 21 studies). This generated a meta-analysis list for ERα or GATA3, which were then compared for overlap. As the overlap was high the stringency was increased to 4 studies (19%), the data of which is used for Table 1. Gene names were obtained using Genecards .
Reporter gene assays
MCF-7 Cells were grown in DMEM (minus phenol-red) with 10% charcoal-stripped FBS. SKBR3 were grown in DMEM with 10% FBS. MUC1 (-881 to +13) was cloned as a KpnI/XhoI fragment, and KRT18 (-2961 to +96) was cloned as a KpnI/BglII fragment. Both were generated by high-fidelity PCR from human genomic DNA (Roche), and were ligated into pGL4.20 (Promega). pS2 reporter has previously been described . Luciferase reporter gene assays were performed using standard protocols. Here 200–400 ng reporter were transfected with 200 ng pcDNA3 or pcDNA3-GATA3, and 3U/well of β-galactosidase protein (Sigma) as transfection efficiency control. Ten nM Tamoxifen (Sigma) was incubated for 14 h prior to cell assay.
Results and Discussion
Using the Oncomine™ integrated cancer profiling database GATA3 and ERα were searched for coexpressing genes . Coexpression data from 21 multi-array studies was extracted and analysed, separately for ERα and GATA3. While these studies varied in cancer-types, the overwhelming majority extracted for analysis were breast-cancer based [Additional file 1 and 2]. The frequency of coexpressing genes over the 21 studies was determined and the cutoff set to 3 studies or more (3 studies = 14% frequency overlap – [see Additional file 1 and 2]). Next, to ascertain the extent GATA3 may play a role in ERα pathways the frequency coexpression lists were compared for overlap. Interestingly, there was an extensive overlap between both GATA3 and ERα lists at the cutoff of 3 studies (Figure 1A). Increasing the cutoff to 4 or more studies (almost one-fifth of the studies) did not change the relative overlap with respect to total gene numbers, with 43% of the number of ERα coexpressed genes, and 56% of GATA3 coexpressed genes represented in the overlap (Figure 1B). The overlap data with the frequency cutoff of 4 studies is shown in Table 1.
Every technique has its caveats, and the limitation here is that we are assessing the common genes that are consistently coexpressed with ERα and GATA3 over many different human cancer studies. This implies that coexpressed genes are in the same pathways as GATA3 and ERα. However, the meta-analyses can only infer association within the same pathways, and pathway coexpression at the RNA level might not necessarily translate into protein level. Nevertheless, our data are strongly supported by previous knowledge of the ERα pathway.
A recent study identified 51 genes significantly upregulated in ERα positive breast tumours, using a real-time PCR based approach . Attesting to the stringency of the meta-analysis approach used here 32 of theses genes were found to overlap with the ERα coexpression list, while an identical number also overlapped with GATA3 (Table 2). This was reflected in a similar study comparing ERα over-expressed transcripts in both oligonucleotide microarray and SAGE platforms , where 27 genes common to the ERα pathway are represented here in our common ERα:GATA3 meta-analysis comparison [see Additional file 3]. These data not only acted as wide-ranging external validation for the individual meta-analyses, but also confirmed the extent of the involvement of GATA3 in ERα pathways.
Furthermore, when compared to a list of genome-wide promoters shown to be bound by ERα in MCF-7 cells  or on chromosomes 21 and 22 , 23 were identified in the ERα meta-analysis list, while 27 were identified within the GATA3 list (Table 3). This again supports both the validity of the meta-analysis technique used here, and the role of GATA3 in ERα pathways. It is also possible that the overlap would be even higher if the ERα genomic location analysis were performed on a pool of human ERα-positive breast tumour samples as opposed to a cell-culture model system. While not to detract from the power of a model system such as MCF-7 there are likely to be a great many differences between a homogeneous cell monolayer and a 3-dimensional cancer made up of a heterogeneous cell population.
Of the 10 classifier genes previously identified in a meta-analysis of ERα, the same 4 were identified in both meta-analyses of this study (ESR1, GATA3, FOXA1, SLC39A6) . Once again this adds credence to the high-quality data obtained in our current meta-analyses.
Implicating GATA3 in control of some of these gene products is a microarray experiment performed after overexpression of GATA3 in 293T cells . After expression of GATA3 elevated levels of TFF1, TFF3, KRT18, FOXA1, SLC9A3R1, TPD52, BCAS1 were observed, all of which we identified here for both GATA3 and ERα meta-analyses. While 293T are not breast cancer cells, it raises the question of how many more of our predicted pathway partners of GATA3 would be identified if the microarray were repeated in cells such as MCF-7 which also retain high ERα expression. In the example of SLC9A3R1 (NHERF1) which is a putative tumour suppressor, it was shown to increase growth of 2 breast cancer cell lines when knocked down by shRNA . If GATA3 does help to control expression of NHERF1 this might be one mechanism consistent with its role in the less-aggressive differentiated luminal A breast cancers. Another example is BCAS1 (NABC1) which is overexpressed in breast carcinomas but downregulated in colorectal tumours [22, 23]. Indeed, overexpression of NABC1 did not result in changes in cell-cycle or anchorage-dependent growth properties in NIH3T3 cells, implying it may not be intrinsically oncogenic .
As GATA3 is expressed in, and regulates, luminal epithelial cells and has also been shown to regulate the MUC1 gene it is no surprise that MUC1 is also mostly expressed in luminal breast epithelial cells as well as other glandular epithelia . MUC1, when abnormally expressed, leads to a loss of both cell-extracellular and cell-cell contacts. It has also been shown that MUC1 levels can be regulated by estrogen and ERα can bind putative binding sites derived from the MUC1 promoter in-vitro . Here we reveal that both GATA3 and ERα coexpress with MUC1 acting as further validation of the meta-analysis technique used here. Furthermore, transfected GATA3 can activate a MUC1 promoter reporter in MCF-7 cells, even in the presence of Tamoxifen i.e. independently to ERα activation. This activation could be repeated in the ERα-negative breast cancer cell line SKBR3 (Figure 2). The activation of ERα pathway genes was also observed with pS2 (TFF1) and KRT18 reporters (Figure 2). These data indicate that GATA3 can have its own impact on the ERα pathway and is not just acting indirectly via ERα.
It has also been postulated that, as the deletion of GATA3 in mammary primordia (by K14-Cre) resulted in an inability to form mammary placodes is similar to that of loss of LEF1, Msx1 and Msx2 these may all be intertwined in a transcriptional network [4, 27]. It is of interest that in our present study we observe MSX2 coexpression both with GATA3 and ERα, which helps to support this notion.
Using the meta-analysis data presented it is easy to build up transcriptional networks such as this and all of the data presented strongly supports (1) the quality of the meta-analysis results, (2) the concept that GATA3 is firmly entrenched within ERα pathways. Future in-depth analysis of the data presented may lead to novel aspects of ERα or GATA3 regulated pathways, and help to understand the etiology of ERα-positive breast cancers, and management of their outcomes.
While GATA3 has been identified previously in a meta-analysis of ERα only 10 genes were identified in total . Here we give an extensive list of coexpressed ERα genes and for the first time a reciprocal meta-analysis for GATA3 has been performed, and the results compared for overlap. This overlap was considerable, confirming the important role of GATA3 in the ERα pathway. The vital question raised is whether GATA3 is crucial to the ERα pathway only by regulation of ERα levels, or through further control of ERα-regulated genes in concert with ERα itself. The GATA3 overexpression microarray experiment in 293T cells, and our reporter gene assays certainly implies the latter . Genome-wide location analysis (ChIP-chip) of GATA3 in a well-established ERα system such as MCF-7 cells, as well as specific analysis of the ERα pathway in GATA3 conditional knockout mice will yield vital information regarding the extent that GATA3 is integral to the ERα pathway.
Ho IC, Pai SY: GATA-3 – not just for Th2 cells anymore. Cell Mol Immunol. 2007, 4: 15-29.
Inman D, Kawana K, Schust D, Lininger R, Young S: Cyclic regulation of T-Bet and GATA-3 in human endometrium. Reprod Sci. 2008, 15: 83-90. 10.1177/1933719107309690
Mehra R, Varambally S, Ding L, Shen R, Sabel MS, Ghosh D, Chinnaiyan AM, Kleer CG: Identification of GATA3 as a breast cancer prognostic marker by global gene expression meta-analysis. Cancer Res. 2005, 65: 11259-11264. 10.1158/0008-5472.CAN-05-2495
Asselin-Labat ML, Sutherland KD, Barker H, Thomas R, Shackleton M, Forrest NC, Hartley L, Robb L, Grosveld FG, Wees van der J, Lindeman GJ, Visvader JE: Gata-3 is an essential regulator of mammary-gland morphogenesis and luminal-cell differentiation. Nat Cell Biol. 2007, 9: 201-209. 10.1038/ncb1530
Kouros-Mehr H, Slorach EM, Sternlicht MD, Werb Z: GATA-3 maintains the differentiation of the luminal cell fate in the mammary gland. Cell. 2006, 127: 1041-1055. 10.1016/j.cell.2006.09.048
Garcia-Closas M, Troester MA, Qi Y, Langerod A, Yeager M, Lissowska J, Brinton L, Welch R, Peplonska B, Gerhard DS, Gram IT, Kristensen V, Borresen-Dale AL, Chanock S, Perou CM: Common genetic variation in GATA-binding protein 3 and differential susceptibility to breast cancer by estrogen receptor alpha tumor status. Cancer Epidemiol Biomarkers Prev. 2007, 16: 2269-2275. 10.1158/1055-9965.EPI-07-0449
van 't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a
Eeckhoute J, Keeton EK, Lupien M, Krum SA, Carroll JS, Brown M: Positive cross-regulatory loop ties GATA-3 to estrogen receptor alpha expression in breast cancer. Cancer Res. 2007, 67: 6477-6483. 10.1158/0008-5472.CAN-07-0746
Schneider J, Ruschhaupt M, Buness A, Asslaber M, Regitnig P, Zatloukal K, Schippinger W, Ploner F, Poustka A, Sultmann H: Identification and meta-analysis of a small gene expression signature for the diagnosis of estrogen receptor status in invasive ductal breast cancer. Int J Cancer. 2006, 119: 2974-2979. 10.1002/ijc.22234
Oh DS, Troester MA, Usary J, Hu Z, He X, Fan C, Wu J, Carey LA, Perou CM: Estrogen-regulated genes predict survival in hormone receptor-positive breast cancers. J Clin Oncol. 2006, 24: 1656-1664. 10.1200/JCO.2005.03.2755
Kouros-Mehr H, Bechis SK, Slorach EM, Littlepage LE, Egeblad M, Ewald AJ, Pai SY, Ho IC, Werb Z: GATA-3 links tumor differentiation and dissemination in a luminal breast cancer model. Cancer Cell. 2008, 13: 141-152. 10.1016/j.ccr.2008.01.011
Wilson BJ, Giguere V: Identification of novel pathway partners of p68 and p72 RNA helicases through Oncomine meta-analysis. BMC Genomics. 2007, 8: 419- 10.1186/1471-2164-8-419
Barry JB, Laganiere J, Giguere V: A single nucleotide in an estrogen-related receptor alpha site can dictate mode of binding and peroxisome proliferator-activated receptor gamma coactivator 1alpha activation of target promoters. Mol Endocrinol. 2006, 20: 302-310. 10.1210/me.2005-0313
Tozlu S, Girault I, Vacher S, Vendrell J, Andrieu C, Spyratos F, Cohen P, Lidereau R, Bieche I: Identification of novel genes that co-cluster with estrogen receptor alpha in breast tumor biopsy specimens, using a large-scale real-time reverse transcription-PCR approach. Endocr Relat Cancer. 2006, 13: 1109-1120. 10.1677/erc.1.01120
Abba MC, Hu Y, Sun H, Drake JA, Gaddis S, Baggerly K, Sahin A, Aldaz CM: Gene expression signature of estrogen receptor alpha status in breast cancer. BMC Genomics. 2005, 6: 37- 10.1186/1471-2164-6-37
Laganiere J, Deblois G, Lefebvre C, Bataille AR, Robert F, Giguere V: From the Cover: Location analysis of estrogen receptor alpha target promoters reveals that FOXA1 defines a domain of the estrogen response. Proc Natl Acad Sci USA. 2005, 102: 11651-11656. 10.1073/pnas.0505575102
Carroll JS, Liu XS, Brodsky AS, Li W, Meyer CA, Szary AJ, Eeckhoute J, Shao W, Hestermann EV, Geistlinger TR, Fox EA, Silver PA, Brown M: Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell. 2005, 122: 33-43. 10.1016/j.cell.2005.05.008
Usary J, Llaca V, Karaca G, Presswala S, Karaca M, He X, Langerod A, Karesen R, Oh DS, Dressler LG, Lonning PE, Strausberg RL, Chanock S, Borresen-Dale AL, Perou CM: Mutation of GATA3 in human breast tumors. Oncogene. 2004, 23: 7669-7678. 10.1038/sj.onc.1207966
Pan Y, Wang L, Dai JL: Suppression of breast cancer cell growth by Na+/H+ exchanger regulatory factor 1 (NHERF1). Breast Cancer Res. 2006, 8: R63- 10.1186/bcr1616
Collins C, Rommens JM, Kowbel D, Godfrey T, Tanner M, Hwang SI, Polikoff D, Nonet G, Cochran J, Myambo K, Jay KE, Froula J, Cloutier T, Kuo WL, Yaswen P, Dairkee S, Giovanola J, Hutchinson GB, Isola J, Kallioniemi OP, Palazzolo M, Martin C, Ericsson C, Pinkel D, Albertson D, Li WB, Gray JW: Positional cloning of ZNF217 and NABC1: genes amplified at 20q13.2 and overexpressed in breast carcinoma. Proc Natl Acad Sci USA. 1998, 95: 8703-8708. 10.1073/pnas.95.15.8703
Correa RG, de Carvalho AF, Pinheiro NA, Simpson AJ, de Souza SJ: NABC1 (BCAS1): alternative splicing and downregulation in colorectal tumors. Genomics. 2000, 65: 299-302. 10.1006/geno.2000.6172
Beardsley DI, Kowbel D, Lataxes TA, Mannino JM, Xin H, Kim WJ, Collins C, Brown KD: Characterization of the novel amplified in breast cancer-1 (NABC1) gene product. Exp Cell Res. 2003, 290: 402-413. 10.1016/S0014-4827(03)00353-7
Abba MC, Nunez MI, Colussi AG, Croce MV, Segal-Eiras A, Aldaz CM: GATA3 protein as a MUC1 transcriptional regulator in breast cancer cells. Breast Cancer Res. 2006, 8: R64- 10.1186/bcr1617
Zaretsky JZ, Barnea I, Aylon Y, Gorivodsky M, Wreschner DH, Keydar I: MUC1 gene overexpressed in breast cancer: structure and transcriptional activity of the MUC1 promoter and role of estrogen receptor alpha (ERalpha) in regulation of the MUC1 gene expression. Mol Cancer. 2006, 5: 57- 10.1186/1476-4598-5-57
Naylor MJ, Ormandy CJ: Gata-3 and mammary cell fate. Breast Cancer Res. 2007, 9: 302- 10.1186/bcr1661
We thank John Coligan, NIH, for pcDNA3-GATA3. Funding was provided by the Canadian Institute for Cancer Research (VG) and a McGill University Health Centre fellowship (BW).
BW conceived and designed the study, performed the meta-analyses, the reporter assays, and wrote the manuscript. VG critically reviewed the manuscript, and approved the final version.