Transcriptional regulation of bidirectional gene pairs by 17-β-estradiol in MCF-7 breast cancer cells

Braz J Med Biol Res, February 2011, Volume 44(2) 112-122

Transcriptional regulation of bidirectional gene pairs by 17β-estradiol in MCF-7 breast cancer cells

S.A.B. Garcia and M.A. Nagai

Disciplina de Oncologia, Departamento de Radiologia, Faculdade de Medicina, Universidade de São Paulo, São Paulo, SP, Brasil

^{Correspondence and Footnotes}

Abstract

Using cDNA microarray analysis, we previously identified a set of differentially expressed genes in primary breast tumors based on the status of estrogen and progesterone receptors. In the present study, we performed an integrated computer-assisted and manual search of potential estrogen response element (ERE) binding sites in the promoter region of these genes to characterize their potential to be regulated by estrogen receptors (ER). Publicly available databases were used to annotate the position of these genes in the genome and to extract a 5’flanking region 2 kb upstream to 2 kb downstream of the transcription start site for transcription binding site analysis. The search for EREs and other binding sites was performed using several publicly available programs. Overall, approximately 40% of the genes analyzed were potentially able to be regulated by estrogen via ER. In addition, 17% of these genes are located very close to other genes organized in a head-to-head orientation with less than 1.0 kb between their transcript units, sharing a bidirectional promoter, and could be classified as bidirectional gene pairs. Using quantitative real-time PCR, we further investigated the effects of 17β-estradiol and antiestrogens on the expression of the bidirectional gene pairs in MCF-7 breast cancer cells. Our results showed that some of these gene pairs, such as TXNDC9/EIF5B, GALNS/TRAPPC2L, and SERINC1/PKIB, are modulated by 17β-estradiol via ER in MCF-7 breast cancer cells. Here, we also characterize the promoter region of potential ER-regulated genes and provide new information on the transcriptional regulation of bidirectional gene pairs.

Key words: Breast cancer; Estrogen receptor; Gene expression; Bidirectional promoters

Introduction

Estrogens are crucial for the development and maintenance of reproductive organs and play important roles in several target tissues, including bone and the cardiovascular and nervous systems. In the mammary gland, the estrogens mediate key physiological processes that are essential for normal growth, differentiation, and survival. A large body of evidence shows that estrogens, especially 17β-estradiol (E₂), also play a pivotal role in promoting mammary carcinogenesis (1). Most of the complex effects of estrogens on cell proliferation, differentiation, and survival are mediated by the estrogen receptors ERα and ERβ via the transcriptional regulation of ER target genes. In the classical mechanism of action, the estrogen receptor complex binds with high affinity to estrogen response elements (EREs) in the promoter region of target genes, interacts with the basal transcription machinery, and regulates gene expression. The consensus ERE, which was first identified in the promoter region of the Xenopus vitellogenin gene, is a perfect palindromic DNA motif composed of two inverted sequences of 5 bp separated by three nucleotides: 5’-GGTCAnnnTGACC-3’ (2). Computational and experimental studies have revealed that most promoter regions in estrogen target genes in the human genome do not contain the consensus ERE, but contain degenerate palindromic sequences showing one or more substitutions, ERE half-sites, direct repeats of half-palindromes, and other transcription factor binding sites, such as AP1, SP1, NF-κB, C/EBP, FoxA1, and octamers (2-5). ER also interacts with EREs containing variations in base substitution, although some nucleotide substitutions can be more harmful than others (6). Some symmetrical substitutions on each side of the palindrome may more dramatically affect the interaction with the ER than a unilateral replacement (6). ERE half-site sequences associated with SP-1 sites, direct repeats of ERE half-sites, and EREs with different spacing between the palindromic sequences also display estrogen responses (7-9).

Genome-wide analysis has indicated that 10% of human genes are organized in a head-to-head orientation and located on opposite DNA strands in close proximity to the transcription start sites (TSS), separated from them by less than 1 kb (10,11). These genes are classified as bidirectional gene pairs and share an intervening sequence containing regulatory elements called bidirectional promoters (12,13). Bidirectional gene pairs are evolutionarily conserved, suggesting that the close proximity might have been maintained for some functional reason (14). One possibility is that genes sharing promoters might be coordinately regulated in different tissues; however, the mechanisms underlying the transcriptional regulation of bidirectional promoters are still poorly understood (15,16). Much evidence indicates that bidirectional gene pairs might play a role in cancer; DNA repair genes that are involved in tumorigenesis are frequently found to be organized as bidirectional gene pairs (10,11). The expression profiles of ovarian and breast cancer show an enrichment of bidirectional gene pairs that include DNA repair genes, such as BRCA1, BRCA2, CKEK1, and FANC family members (17). Furthermore, the silencing of bidirectional gene pairs by DNA methylation has been reported in various human cancer cell lines and primary breast tumors (18-20).

In the present study, we performed an integrated computational and manual search to identify potential ER binding sites in the promoter region of 83 genes previously identified as being differentially expressed in primary breast tumors based on the presence of estrogen and progesterone receptors (21). We identified a subgroup of genes arranged in a head-to-head orientation and separated by a small intergenic region that could be classified as bidirectional promoters. Selected candidate ER target genes were further examined for their potential to be modulated by 17β-estradiol in MCF-7 breast cancer cells.

Material and Methods

ERE binding site search

The NCBI (http://www.ncbi.nlm.nih.gov), USCS Genome Bioinformatics (http://genome.uscs.edu) and ENSEMBL (http://www.ensembl.org) databases were used to annotate the position of 83 genes in the genome and to extract sequences 2 kb upstream to 2 kb downstream of the TSS for binding site analysis. The search for EREs and other binding sites was performed using several publicly available programs (Transfac 6.0, http://www.gene-regulation.com, and Dragon Genome Explorer, http://research.i2.a-star.edusg/promoter/). A list of ERE sequences was prepared from the reports extracted from PubMed and used for a manual search (6,22-24). The present study was approved by the Institutional Ethics Committee.

Cell lines and culture conditions

Human breast cancer cells, MCF-7 and MDA-MB-231, were obtained from the American Type Culture Collection (ATCC). Cells were cultured at 37°C in phenol red-free RPMI supplemented with 10% FCS and 1% penicillin-streptomycin (Gibco, USA) in an atmosphere of 5% CO₂/95% air. Before the treatments, cells were washed twice with PBS to remove residual serum and were grown in phenol red-free media containing 5% charcoal-stripped FCS for 48 h. After hormone deprivation, the cells were incubated with 10 nM 17β-estradiol for 2, 6, and 24 h or treated with 1 µM 4-hydroxytamoxifen or 1 µM ICI 182.780 (Faslodex) alone or in combination with 17β-estradiol for 24 h. Cells maintained in phenol red-free media and supplemented with 5% FBS were used as controls.

RNA extraction

After the treatments, the cells were washed twice with PBS and harvested. Total RNA was extracted from cultured cells using the guanidine isothyocyanate method (25). RNA samples were analyzed using 1% agarose gel electrophoresis and ethidium bromide staining. The quality of the RNA samples was estimated by the ratio of 28S:18S (2.0 and higher). All RNA samples were treated with DNaseI for 1 h at 37°C to eliminate genomic DNA contamination.

Quantitative real-time PCR (qPCR)

qPCR was performed using the GeneAmp 7500 Sequence Detection System (Applied Biosystems, USA). cDNA was generated using the High Capacity cDNA Archive kit (Applied Biosystems), and each cDNA sample was analyzed in duplicate. PCR was carried out in a total volume of 50 µL according to manufacturer instructions for SYBR Green PCR Core reagent (PE Applied Biosystems). Primer pairs were designed using the Primer3 software (http://frodo.wi.mit.edu/primer3/). The primers used for qPCR analysis were as follows: DIDO1, 5’-AGAAAGCTGCAAATCC AAGTG-3’ and 5’-GGCGTGTTTGAGGATACAGTC-3’; C20orf11, 5’-TTTGACAGTCCCGAGGAGTC-3’ and 5’-GTGTTGACTCGCGATTTTCAT-3’; TXNDC9, 5’-GTG AAAATGTGGTTTGCCATT-3’ and 5’-TGCTTTTTCCACATT CAGCTT-3’; EIF5B, 5’-AGATTCGAGGCCTCCTGTTAC-3’ and 5’-GGGTAAACCAGCCAATGTTTT-3’; GALNS, 5’-CAGGATCACCTCGGTCGT-3’ and 5’-GAATGGATTCTG GAGGTGTCA-3’; TRAPPC2L, 5’-GTCACCAACTCCAA GGTGAAG-3’ and 5’-TTGTAGAAGGGGTTGCACATC-3’; SERINC1, 5’-GTGATGGATCACTGGAGGATG-3’ and 5’-TCACGAGAGGGTTCATACCTG-3’; PKIB, 5’-CGG AATGCCTTACCAGACATC-3’ and 5’-TTCTTCATTTTG AGGCTTTTCC-3’; GTPBP2, 5’-CCTACACACCACCAT CTTTGG-3’ and 5’-ATGAAGAAGGGCACTTTCAGG-3’; MAD2L1BP, 5’-GCAGAGGAGATGCTGAAGAAG-3’ and 5’-GGCTGAGGACACTCTCCAGTT-3’; GAPDH, 5’-CCTCC AAAATCAAGTGGGGCG-3’ and 5’-GGGGCAGAGATG ATGACCCTT-3’.

The relative gene expression was normalized using GAPDH expression as an internal control. MCF-7 cells maintained in phenol red-free medium and supplemented with 5% FCS served as calibrator samples for hormone-treated cells. The results are reported as n-fold differences in the expression of the target gene relative to the expression of GAPDH and a calibrator sample. The relative expression was calculated by 2^-DDCT, where CT is the fluorescence threshold value, DCT = CT of the target gene - CT of the reference gene, and DDCT = DCT of the tumor sample - DCT of the calibrator sample.

Statistical analysis

Statistical analyses were performed by one-way ANOVA using the SPSS 10.0 software (Statistical Package for the Social Sciences; SPSS Inc., USA), and P values were calculated using the Bonferroni test and considered to be significant at P < 0.05.

Results

In a previous study using cDNA microarray, we identified differentially expressed transcripts between ER/progesterone receptor (PR)-positive and ER/PR-negative primary breast tumors, providing a new set of potential target genes for the ER signaling pathway. In the present study, in order to identify DNA binding sites that mediate the estrogen response, we performed a computer-assisted and manual search for putative ERE binding sites in the promoter region extracted from 83 of these potential ER target genes. The manual search was performed using a list of putative ERE binding sites based on data available in the literature. A summary of the results obtained from this analysis is shown in Figure 1. Overall, 43% (36/83) of the candidate genes showed a good potential for being regulated by estrogen via ER. The promoter regions of these genes contain perfect or imperfect palindromic sequences, SP1, AP1, and/or FOXA1 binding sites.

Interestingly, during the annotation of the 83 potential ER target genes, we identified 17 genes located in close proximity to other genes and organized in a head-to-head orientation, possibly sharing the same promoter region. Using criteria of an intergenic distance of <1000 bp spanning the TSS of the divergent gene pair to define bidirectional promoters (11), we identified 14 (16.8%) bidirectional gene pairs (Table 1). The size of the intergenic sequences ranged from 18 to 973 bp, and all of these bidirectional promoters contain CpG-rich sequences. Figure 2 shows a schematic representation of three of the bidirectional gene pairs identified. We also investigated whether the identified bidirectional gene pairs maintained their divergent head-to-head organization in other species. Although limited, our analysis showed that the majority of the bidirectional gene pairs found in this study is evolutionarily conserved, as shown in Figure 3 for the GALNS/TRAPPC2L divergent gene pair.

Our computer-assisted search for EREs revealed that the promoter regions of the majority of the bidirectional gene pairs showed good potential to be regulated by estrogen. As shown in Table 1, all 14 bidirectional promoters analyzed contain putative ERE binding sites (at least half EREs) and other binding sites for other transcription factors, such as SP1 and AP1 or FOXA1 binding sites, indicating that they have the potential to be regulated by ER. The putative EREs were also found to be conserved in different species (Figure 3). In addition, approximately 80% of these bidirectional promoters contain GABP binding sites (CCGGAA(GTG)). We further evaluated the expression of five of these bidirectional gene pairs in response to 17β-estradiol in MCF-7 breast cancer cells. The gene pairs GALNS/TRAPPC2L, TXNDC9/EIF5B, and SERINC1/PKIB were modulated by estradiol (Figure 4). However, the expression of the gene pairs GTPBP2/MAD2L1BP and DIDO1/C20orf11 was not affected in the MCF-7 cells exposed to 17β-estradiol or to the antiestrogens tamoxifen and ICI 182,780 (data not shown). The transcripts of GALNS and TRAPPC2L showed increased expression in MCF-7 cells maintained in stripped serum (deprived of steroids and growth factors) compared to the cells grown in normal serum (control cells) and were negatively modulated in the presence of 10 nM 17β-estradiol for 24 h (Figure 4A). A similar reduction in these transcripts was observed in MCF-7 cells treated with 1 μM tamoxifen or 1 μM tamoxifen plus 100 nM 17β-estradiol compared to the expression in MCF-7 cells maintained in stripped serum. However, in MCF-7 cells treated with 1 μM ICI 182,780 alone the GALNS/TRAPPC2L transcripts showed increased expression compared to control cells and remained elevated in cells treated with 1 μM ICI 182,780 plus 100 nM 17β-estradiol. The transcripts of the bidirectional gene pair TXNDC9 and EIF5B were positively regulated by 10 nM 17β-estradiol after 6 h of treatment (Figure 4B). In contrast to the gene pair GALNS/TRAPPC2L, the TXNDC9 and EIF5B transcripts were not significantly altered in MCF-7 cells grown in stripped serum compared to control cells, but both transcripts were significantly induced (>2-fold) in MCF-7 cells exposed to 10 nM 17β-estradiol for 6 h. After 24 h of 17β-estradiol treatment, the expression of these transcripts was similar to that observed in control cells. However, the transcripts of the TXNDC9 and EIF5B gene pair were not affected by treatment with the antiestrogens tamoxifen and ICI 182,780 alone or in combination with 17β-estradiol for 24 h compared to cells growing in normal or stripped serum. Interestingly, the transcripts of the bidirectional gene pair SERINC1/PKIB were divergently modulated by 17β-estradiol (Figure 4C). Expression of the PKIB transcripts was reduced in MCF-7 cells maintained in stripped serum and was substantially stimulated by 10 nM 17β-estradiol after 6 h of treatment, whereas the SERINC1 transcripts increased in MCF-7 cells grown in stripped serum and expression was significantly inhibited in the presence of 17β-estradiol. In addition, cells treated with 1 μM tamoxifen showed a reduction of protein kinase inhibitor beta (PKIB) levels, which was reversed in the presence of 100 nM 17β-estradiol. Treatment with ICI 182,780 dramatically reduced the expression of PKIB but did not affect the expression of SERINC1 transcripts compared to cells grown in stripped serum.

To determine whether 17β-estradiol or antiestrogens are able to modulate the transcripts of GALNS/TRAPPC2L, TXNDC9/EIF5B, and SERINC1/PKIB in ER-negative breast cancer cells, MDA-MD-231 cells was also treated with 17β-estradiol and/or antiestrogens for 6 and 24 h. Estrogen or antiestrogens were unable to regulate the transcripts of these bidirectional gene pairs (data not shown) in MDA-MB-231 cells.

Figure 1.
Summary of the results of the estrogen response element (ERE) search for 83 potential estrogen receptor target genes. The NCBI (http://www.ncbi.nlm.nih.gov) and USCU Genome Bioinformatics (http://www.genome.uscs.edu) databases were used to annotate the position of the 83 differentially expressed genes in the genome and to extract sequences 2 kb upstream and 2 kb downstream of the start sites for binding site analyses. The search for EREs and other binding sites was performed using several publicly available programs (Transfac 6.0, http://www.gene-regulation.com, and Dragon Genome Explorer, http://research.i2.a-star.edusg/promoter/). ¹Palindromic ERE consensus derived from vitellogenin A2 (GGTCAnnnTGACC) separated by 1 to 7 bp. ²Imperfect palindromic ERE carrying a single substitution. ³Imperfect palindromic ERE carrying two substitutions. ⁴Direct repeats separated by 3 bp.

Figure 2.
Schematic presentation of the bidirectional gene pairs GALNS/TRAPPC2L, TXNDC9/EIF5B, and SERINC1/PKIB. The transcriptional units of the bidirectional gene pairs are organized in a head-to-head orientation on chromosomes 16q24.3 (GALNS/TRAPPC2L), 2q11.2 (TXNDC9/EIF5B), and 16q24.3 (SERINC1/PKIB) (http://genome.uscs.edu). The transcript start site (TSS) for each gene is indicated by black arrows, which show the transcriptional direction. Open boxes indicate the first exons of each gene. Black boxes indicate intergenic sequences for each gene pair. The boxes with vertical bars represent CpG-rich sequences.

Figure 3.
Evolutionarily conserved estrogen receptor binding elements (EREs) in the potential promoter region of the bidirectional gene pair GALNS/TRAPPC2L across vertebrates. A, Schematic presentation of the bidirectional gene pair GALNS/TRAPPC2L in humans, chimpanzees, dogs, rats, and chickens. The black arrows and numbers 1 to 5 indicate the location of conserved EREs relative to the position of the transcript start site in the TRAPPC2L gene. B, 1 to 5 ERE sequences; EREs 2 (C) and 5 (D) and the surrounding sequences are shown in the alignment. The half-EREs are in bold. Asterisks denote highly conserved bases.

Figure 4.
Effects of 17β-estradiol (E₂) and the anti-estrogens tamoxifen (T) and ICI 182,780 (ICI), on the transcriptional regulation of the bidirectional gene pairs GALNS/TRAPPC2L (A), TXNDC9/EIF5B (B), and SERINC1/PKIB (C) in MCF-7 breast cancer cells. Synchronized MCF-7 cells were cultured in: NS, medium supplemented with 5% fetal calf serum (FCS); SS, medium supplemented with 5% charcoal-stripped FCS; E₂ - 2 h, exposed to E₂ (10 nM) for 2 h; E₂ - 6 h, exposed to 10 nM E₂ for 6 h; E₂ - 24 h, exposed to 10 nM E₂ for 24 h; T - 24 h, exposed to 1 μM tamoxifen for 24 h; T+E₂ - 24 h, exposed to 1 μM tamoxifen plus 100 nM E₂ for 24 h; ICI - 24 h, exposed to ICI 182,780 for 24 h, or ICI+E₂ - 24 h, exposed to ICI 182,780 plus E₂ (100 nM) for 24 h. Data are reported as means ± SD of three experiments. *P < 0.05 compared to NS and/or SS (one-way ANOVA followed by the Bonferroni post-test).

Table 1.
List of the bidirectional gene pairs identified.

Discussion

In the present study, we employed an integrated computer-assisted and manual search to identify EREs in the promoter region of a set of candidate ER target genes previously identified as being differentially expressed in primary breast tumors based on their estrogen and progesterone status (21). Using this approach, we characterized the promoter region of 83 genes for their potential to be regulated by ER. Although only a small percentage of the candidate ER target genes showed the presence of an ERE consensus sequence in their promoter region, more than 40% of these genes were found to have putative ERE, SP1, AP1, and/or FOXA1 binding sites in their promoter regions, indicating the potential to be regulated by estrogen. The maximum ER transcriptional transactivation is thought to be displayed by the interaction of a homo- or heterodimer of ER with the ERE consensus sequence (CAGGTCAnnnTGACCTG). However, the ERE consensus sequence is rarely found in the promoter regions of human genes (2). The promoter regions of estrogen-responsive genes usually contain more than one imperfect palindromic sequence separated by one, two, three, or more interstitial bases, direct repeats of ERE half-sites, and widespread ERE half-sites (2). Half-EREs have been shown to act synergistically, but the degree of synergy between EREs is modulated by the physical distance between EREs and by the distance between the ERE and the TATA-box (26). In addition, perfect or imperfect ERE half-sites adjacent to AP1 or SP1 binding sites in the promoter region of known estrogen-regulated genes, such as HSP27 and PR, were also found to be effective in conferring estrogen responsiveness (27,28).

During the process of gene annotation, we identified a subgroup of divergent gene pairs organized in a head-to-head orientation that could be classified as bidirectional gene pairs. Surprisingly, the prevalence of bidirectional gene pairs found in the set of candidate ER target genes analyzed was higher than the predicted frequency of 9.4-11% found in the human genome by computational genome-wide studies (11,29). An incremental prevalence of gene pairs organized in a head-to-head orientation has also been observed among the genes that encode proteins that function in DNA repair (10,17).

Bidirectional gene pairs show evolutionary conservation, which suggests a selective pressure to facilitate the expression of functionally related genes or those that need to be co-regulated by the same stimulus or in a cell-specific manner (10,11,29). In agreement with the literature, the majority of the bidirectional gene pairs identified in the present study were shown to be conserved during evolution, keeping the same organization in vertebrates and indicating that these genes organized in a head-to-head orientation play important roles (29,30). Our data also showed that some of the putative ERE binding sites were conserved in the bidirectional gene pairs in other species. Different degrees of conservation of ERE or ERE-SP1 sites have been reported in the literature (4,9,31). In addition, the majority of the bidirectional promoters analyzed also contain GABP binding sites, which were previously shown to be over-represented in bidirectional promoters (32,33). However, whether or not they are all likely to be functionally conserved is still being questioned.

The bidirectional head-to-head orientation suggests that the potential promoter regions of these gene pairs share cis-elements, such as ERE and SP1 sites, that could be transcriptionally regulated by estrogen. To validate our results, we evaluated the effects of estrogen and antiestrogens on the transcriptional regulation of five bidirectional gene pairs in the breast cancer cell line MCF-7. To our surprise, the expression of two of the bidirectional gene pairs, GTPBP2/MAD2L1BP and DIDO1/C20orf11, was not modulated in MCF-7 cells exposed to estrogen and antiestrogens. Since the activity of transcription factors in transcriptional transactivation depends on the state of the chromatin (34,35), we can speculate that, in MCF-7 cells, the E₂-ER complex has no accessibility to the binding sites in the promoter regions of these genes, possibly due to the state of chromatin organization. Further studies using different cell lines and agents that affect chromatin organization are required. On the other hand, distinct patterns of gene regulation by estrogen and antiestrogens were observed for the bidirectional gene pairs GALNS/TRAPPC2L, TXNDC9/EIF5B, and SERINC1/PKIB in MCF-7 cells. In contrast, in ER-negative breast cancer cells, MDA-MB-231, 17β-estradiol or antiestrogens (tamoxifen or ICI) were unable to modulate the mRNA expression of these bidirectional gene pairs (data not shown), indicating that the presence of ER is essential for the transcriptional transactivation of these bidirectional promoters by estrogen. In MCF-7 cells, two of the gene pairs analyzed displayed a positive correlation, with both gene transcripts being up-regulated (TXNDC9/EIFB5) or down-regulated (GALNS/TRAPPC2L) by estrogen. The differences observed in the effects of the antiestrogens tamoxifen and ICI 182,780 on the transcriptional regulation of these bidirectional gene pairs can be attributed to their mechanisms of action. Tamoxifen is a selective estrogen modulator that acts by inhibiting estrogen binding to its main target estrogen receptor; however, it also acts as an estrogen agonist in certain tissues or in specific regulation of ER-responsive genes (36). Compared to tamoxifen, ICI 182,780 is considered to be a pure antagonist that blocks ER-mediated transactivation down-regulating ER protein levels. Additionally, the gene pair SERINC1/PKIB showed a negative correlation with the up-regulation of PKIB transcripts and down-regulation of SERINC1 transcripts when MCF-7 cells were exposed to estrogen. The SERINC1/PKIB gene pairs may have distinct elements in the bidirectional promoter and further studies using reporter assays are needed to better define the critical sequences involved in this divergent regulation. On the other hand, in this study we evaluated the presence of EREs in the proximal promoter region only; however, the presence of EREs located in distal regions or of other cis-regulatory elements regulated by other factors may be necessary to modulate the divergent expression of this bidirectional gene pair.

Although genome-wide analysis showed that 45-68% of the gene pairs organized in a head-to-head orientation are co-expressed (16,29), the mechanisms involved in the transcriptional regulation of these bidirectional gene pairs appear to be complex and still poorly understood (13,29,37,38). Diverse patterns of gene regulation have been reported for different bidirectional gene pairs, including co-regulation, and divergent or asymmetric regulation (20,30,39,40).

In the present study, we characterized the promoter region of potential ER-regulated genes and provided new information on the transcriptional regulation of gene pairs organized in a head-to-head orientation. Further clinical and experimental studies are required to determine whether these genes play a role in breast cancer. The identification of ER target genes may improve our understanding of the role played by estrogens in breast cancer, permitting a better tailored hormone dependence and improved patient response to hormonal therapy.

References

1. Medina D. Mammary developmental fate and breast cancer risk. Endocr Relat Cancer 2005; 12: 483-495.

2. Gruber CJ, Gruber DM, Gruber IM, Wieser F, Huber JC. Anatomy of the estrogen response element. Trends Endocrinol Metab 2004; 15: 73-78.

3. Kushner PJ, Agard DA, Greene GL, Scanlan TS, Shiau AK, Uht RM, et al. Estrogen receptor pathways to AP-1. J Steroid Biochem Mol Biol 2000; 74: 311-317.

4. Carroll JS, Liu XS, Brodsky AS, Li W, Meyer CA, Szary AJ, et al. Chromosome-wide mapping of estrogen receptor binding reveals long-range regulation requiring the forkhead protein FoxA1. Cell 2005; 122: 33-43.

5. Carroll JS, Brown M. Estrogen receptor target gene: an evolving concept. Mol Endocrinol 2006; 20: 1707-1714.

6. Bourdeau V, Deschenes J, Metivier R, Nagai Y, Nguyen D, Bretschneider N, et al. Genome-wide identification of high-affinity estrogen response elements in human and mouse. Mol Endocrinol 2004; 18: 1411-1427.

7. Aumais JP, Lee HS, DeGannes C, Horsford J, White JH. Function of directly repeated half-sites as response elements for steroid hormone receptors. J Biol Chem 1996; 271: 12568-12577.

8. Anderson I, Gorski J. Estrogen receptor alpha interaction with estrogen response element half-sites from the rat prolactin gene. Biochemistry 2000; 39: 3842-3847.

9. O’Lone R, Frith MC, Karlsson EK, Hansen U. Genomic targets of nuclear estrogen receptors. Mol Endocrinol 2004; 18: 1859-1875.

10. Adachi N, Lieber MR. Bidirectional gene organization: a common architectural feature of the human genome. Cell 2002; 109: 807-809.

11. Trinklein ND, Aldred SF, Hartman SJ, Schroeder DI, Otillar RP, Myers RM. An abundance of bidirectional promoters in the human genome. Genome Res 2004; 14: 62-66.

12. Takai D, Jones PA. Origins of bidirectional promoters: computational analyses of intergenic distance in the human genome. Mol Biol Evol 2004; 21: 463-467.

13. Yang MQ, Elnitski LL. Prediction-based approaches to characterize bidirectional promoters in the mammalian genome. BMC Genomics 2008; 9 (Suppl 1): S2.

14. Yang MQ, Taylor J, Elnitski L. Comparative analyses of bidirectional promoters in vertebrates. BMC Bioinformatics 2008; 9 (Suppl 6): S9.

15. Zanotto E, Shah ZH, Jacobs HT. The bidirectional promoter of two genes for the mitochondrial translational apparatus in mouse is regulated by an array of CCAAT boxes interacting with the transcription factor NF-Y. Nucleic Acids Res 2007; 35: 664-677.

16. Yang MQ, Elnitski LL. Diversity of core promoter elements comprising human bidirectional promoters. BMC Genomics 2008; 9 (Suppl 2): S3.

17. Yang MQ, Koehly LM, Elnitski LL. Comprehensive annotation of bidirectional promoters identifies co-regulation among breast and ovarian cancer genes. PLoS Comput Biol 2007; 3: e72.

18. Shu J, Jelinek J, Chang H, Shen L, Qin T, Chung W, et al. Silencing of bidirectional promoters by DNA methylation in tumorigenesis. Cancer Res 2006; 66: 5077-5084.

19. Lin JC, Jeong S, Liang G, Takai D, Fatemi M, Tsai YC, et al. Role of nucleosomal occupancy in the epigenetic silencing of the MLH1 CpG island. Cancer Cell 2007; 12: 432-444.

20. De Souza Santos E, De Bessa SA, Netto MM, Nagai MA. Silencing of LRRC49 and THAP10 genes by bidirectional promoter hypermethylation is a frequent event in breast cancer. Int J Oncol 2008; 33: 25-31.

21. Nagai MA, Da Ros N, Neto MM, de Faria Junior SR, Brentani MM, Hirata R Jr, et al. Gene expression profiles in breast tumors regarding the presence or absence of estrogen and progesterone receptors. Int J Cancer 2004; 111: 892-899.

22. Klinge CM. Estrogen receptor interaction with estrogen response elements. Nucleic Acids Res 2001; 29: 2905-2919.

23. Bajic VB, Tan SL, Chong A, Tang S, Strom A, Gustafsson JA, et al. Dragon ERE Finder version 2: A tool for accurate detection and analysis of estrogen response elements in vertebrate genomes. Nucleic Acids Res 2003; 31: 3605-3607.

24. Vega VB, Lin CY, Lai KS, Kong SL, Xie M, Su X, et al. Multiplatform genome-wide identification and modeling of functional human estrogen receptor binding sites. Genome Biol 2006; 7: R82.

25. Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987; 162: 156-159.

26. Sathya G, Li W, Klinge CM, Anolik JH, Hilf R, Bambara RA. Effects of multiple estrogen responsive elements, their spacing, and location on estrogen response of reporter genes. Mol Endocrinol 1997; 11: 1994-2003.

27. Porter W, Saville B, Hoivik D, Safe S. Functional synergy between the transcription factor Sp1 and the estrogen receptor. Mol Endocrinol 1997; 11: 1569-1580.

28. Petz LN, Ziegler YS, Schultz JR, Kim H, Kemper JK, Nardulli AM. Differential regulation of the human progesterone receptor gene through an estrogen response element half site and Sp1 sites. J Steroid Biochem Mol Biol 2004; 88: 113-122.

29. Li YY, Yu H, Guo ZM, Guo TQ, Tu K, Li YX. Systematic analysis of head-to-head gene organization: evolutionary conservation and potential biological relevance. PLoS Comput Biol 2006; 2: e74.

30. Nakagawa H, Tategu M, Yamauchi R, Sasaki K, Sekimachi S, Yoshida K. Transcriptional regulation of an evolutionary conserved intergenic region of CDT2-INTS7. PLoS One 2008; 3: e1484.

31. Pan YF, Wansa KD, Liu MH, Zhao B, Hong SZ, Tan PY, et al. Regulation of estrogen receptor-mediated long range transcription via evolutionarily conserved distal response elements. J Biol Chem 2008; 283: 32977-32988.

32. Collins PJ, Kobayashi Y, Nguyen L, Trinklein ND, Myers RM. The ets-related transcription factor GABP directs bidirectional transcription. PLoS Genet 2007; 3: e208.

33. Lin JM, Collins PJ, Trinklein ND, Fu Y, Xi H, Myers RM, et al. Transcription factor binding and modified histones in human bidirectional promoters. Genome Res 2007; 17: 818-827.

34. Biddie SC, John S, Hager GL. Genome-wide mechanisms of nuclear receptor action. Trends Endocrinol Metab 2010; 21: 3-9.

35. Hager GL, McNally JG, Misteli T. Transcription dynamics. Mol Cell 2009; 35: 741-753.

36. Lewis JS, Jordan VC. Selective estrogen receptor modulators (SERMs): mechanisms of anticarcinogenesis and drug resistance. Mutat Res 2005; 591: 247-263.

37. Chen PY, Chang WS, Chou RH, Lai YK, Lin SC, Chi CY, et al. Two non-homologous brain diseases-related genes, SERPINI1 and PDCD10, are tightly linked by an asymmetric bidirectional promoter in an evolutionarily conserved manner. BMC Mol Biol 2007; 8: 2.

38. Tsai HK, Su CP, Lu MY, Shih CH, Wang D. Co-expression of adjacent genes in yeast cannot be simply attributed to shared regulatory system. BMC Genomics 2007; 8: 352.

39. Bellizzi D, Dato S, Cavalcante P, Covello G, Di Cianni F, Passarino G, et al. Characterization of a bidirectional promoter shared between two human genes related to aging: SIRT3 and PSMD13. Genomics 2007; 89: 143-150.

40. Voth H, Oberthuer A, Simon T, Kahlert Y, Berthold F, Fischer M. Co-regulated expression of HAND2 and DEIN by a bidirectional promoter with asymmetrical activity in neuroblastoma. BMC Mol Biol 2009; 10: 28.

Acknowledgments

We thank Emerson de Souza Santos for technical assistance. Research supported by FAPESP (#06/01026-0) and, in part, by CNPq (#304949/2006).

Address for correspondence: M.A. Nagai, Disciplina de Oncologia, Departamento de Radiologia, FM, USP, Av. Dr. Arnaldo, 455, 4º andar, 01246-903 São Paulo, SP, Brasil. Fax: [phone]. E-mail: [e-mail]

Received September 10, 2010. Accepted December 8, 2010. Available online December 24, 2010. Publishec February 7, 2011 .

The Brazilian Journal of Medical and Biological Research is partially financed by

Correspondence and Footnotes

Publication Dates

• Publication in this collection
  04 Feb 2011
• Date of issue
  Feb 2011

History

• Received
  10 Sept 2010
• Accepted
  08 Dec 2010