Influence of cell cycle on responses of MCF-7 cells to benzo[a]pyrene.

BACKGROUND: Benzo[a]pyrene (BaP) is a widespread environmental genotoxic carcinogen that damages DNA by forming adducts. This damage along with activation of the aryl hydrocarbon receptor (AHR) induces complex transcriptional responses in cells. To investigate whether human cells are more susceptible to BaP in a particular phase of the cell cycle, synchronised breast carcinoma MCF-7 cells were exposed to BaP. Cell cycle progression was analysed by flow cytometry, DNA adduct formation was assessed by 32P-postlabeling analysis, microarrays of 44K human genome-wide oligos and RT-PCR were used to detect gene expression (mRNA) changes and Western blotting was performed to determine the expression of some proteins, including cytochrome P450 (CYP) 1A1 and CYP1B1, which are involved in BaP metabolism. RESULTS: Following BaP exposure, cells evaded G1 arrest and accumulated in S-phase. Higher levels of DNA damage occurred in S- and G2/M- compared with G0/G1-enriched cultures. Genes that were found to have altered expression included those involved in xenobiotic metabolism, apoptosis, cell cycle regulation and DNA repair. Gene ontology and pathway analysis showed the involvement of various signalling pathways in response to BaP exposure, such as the Catenin/Wnt pathway in G1, the ERK pathway in G1 and S, the Nrf2 pathway in S and G2/M and the Akt pathway in G2/M. An important finding was that higher levels of DNA damage in S- and G2/M-enriched cultures correlated with higher levels of CYP1A1 and CYP1B1 mRNA and proteins. Moreover, exposure of synchronised MCF-7 cells to BaP-7,8-diol-9,10-epoxide (BPDE), the ultimate carcinogenic metabolite of BaP, did not result in significant changes in DNA adduct levels at different phases of the cell cycle. CONCLUSIONS: This study characterised the complex gene response to BaP in MCF-7 cells and revealed a strong correlation between the varying efficiency of BaP metabolism and DNA damage in different phases of the cell cycle. Our results suggest that growth kinetics within a target-cell population may be important determinants of susceptibility and response to a genotoxic agent.

Elevated NCOR1 disrupts PPARalpha/gamma signaling in prostate cancer and forms a targetable epigenetic lesion.

The loss of anti-proliferative responsiveness in prostate cancer cell lines toward ligands for vitamin D receptor, retinoic acid receptors/retinoid X receptors and peroxisome proliferator activated receptor (PPAR)alpha/gamma may entail underlying epigenetic events, as ligand insensitivity reflects significantly altered messenger RNA expression of corepressors and histone-modifying enzymes. Expression patterns were dependent on phases of the cell cycle and associated with repressed basal gene expression of vitamin D receptor and PPARalpha/gamma target genes, for example CDKN1A [encodes p21((waf1/cip1))]. Elevated nuclear corepressor 1 (NCOR1) and nuclear corepressor 2/silencing mediator of retinoic acid and thyroid hormone receptor protein levels were detected in prostate cancer cell lines compared with non-malignant counterparts. Knockdown of the corepressor NCOR1 significantly elevated basal expression of a cohort of target genes, including CDKN1A. Both chemical [histone deacetylases inhibitor (HDACi)] and NCOR1 knockdown targeting enhanced anti-proliferative sensitivity toward PPARalpha/gamma ligands in prostate cancer cell lines. Pursuing PPARalpha/gamma signaling, microarray approaches were undertaken to identify pathways and genes regulated uniquely by a combination of PPARalpha/gamma activation and HDAC inhibition. Again, HDACi and knockdown approaches demonstrated that elevated NCOR1 expression and activity distorted PPARalpha/gamma gene targets centered on, for example cell cycle control, including CDKN1A and TGFBRAP1. Quantitative real time polymerase chain reaction validation and chromatin immunoprecipitation assays both confirmed that elevated NCOR1 disrupted the ability of PPARalpha/gamma to regulate key target genes (CDKN1A and TGFBRAP1). Interrogation of these relationships in prostate cancer samples using principal component and partial correlation analyses established significant interdependent relationships between NCOR1-PPARalpha/gamma and representative target genes, independently of androgen receptor expression. Therefore, we conclude that elevated NCOR1 distorts the actions of PPARalpha/gamma selectively and generates a potential epigenetic lesion with diagnostic and prognostic significance.

Identification through microarray gene expression analysis of cellular responses to benzo(a)pyrene and its diol-epoxide that are dependent or independent of p53.

Human colon carcinoma cells (HCT116) differing in p53 status were exposed to benzo(a)pyrene (BaP) or anti-benzo(a)pyrene-trans-7,8-dihydrodiol-9,10-epoxide (BPDE) and their gene expression responses compared by complementary DNA microarray technology. Exposure of cells to BPDE for up to 24 h resulted in gene expression profiles more distinguishable by duration of exposure than by p53 status, although a subset of genes were identified that had significantly different expression in p53 wild-type (WT) cells relative to p53-null cells. Apoptotic signalling genes were up-regulated in p53-WT cells but not in p53-null cells and, consistent with this, reduced viability and caspase activity were also p53 dependent. BPDE modulated cell cycle and histone genes in both cell lines and, in agreement with this, both cell lines accumulated in S phase. In p53-WT cells, G(2) arrest was also evident, which was associated with accumulation of CDKN1A. Regardless of p53 status, exposure to BaP for up to 48 h had subtle effects on gene transcription and had no influence on cell viability or cell cycle. Interestingly, DNA adduct formation after BaP, but not BPDE, exposure was p53 dependent with 10-fold lower levels detected in p53-null cells. Other cell lines were investigated for BaP-DNA adduct formation and in these the effect of p53 knockdown was also to reduce adduct formation. Taken together, these results give further insight into the role of p53 in the response of human cells to BaP and BPDE and suggest that loss of this tumour suppressor can influence the metabolic activation of BaP.

Accurate prediction of BRCA1 and BRCA2 heterozygous genotype using expression profiling after induced DNA damage.

PURPOSE: In this study, the differential gene expression changes following radiation-induced DNA damage in healthy cells from BRCA1/BRCA1 mutation carriers have been compared with controls using high-density microarray technology. We aimed to establish if BRCA1/BRCA2 mutation carriers could be distinguished from noncarriers based on expression profiling of normal cells. EXPERIMENTAL DESIGN: Short-term primary fibroblast cultures were established from skin biopsies from 10 BRCA1 and 10 BRCA2 mutation carriers and 10 controls, all of whom had previously had breast cancer. The cells were subjected to 15 Gy ionizing irradiation to induce DNA damage. RNA was extracted from all cell cultures, preirradiation and at 1 hour postirradiation. For expression profiling, 15 K spotted cDNA microarrays manufactured by the Cancer Research UK DNA Microarray Facility were used. Statistical feature selection was used with a support vector machine (SVM) classifier to determine the best feature set for predicting BRCA1 or BRCA2 heterozygous genotype. To investigate prediction accuracy, a nonprobabilistic classifier (SVM) and a probabilistic Gaussian process classifier were used. RESULTS: In the task of distinguishing BRCA1 and BRCA2 mutation carriers from noncarriers and from each other following radiation-induced DNA damage, the SVM achieved 90%, and the Gaussian process classifier achieved 100% accuracy. This effect could not be achieved without irradiation. In addition, the SVM identified a set of BRCA genotype predictor genes. CONCLUSIONS: We conclude that after irradiation-induced DNA damage, BRCA1 and BRCA2 mutation carrier cells have a distinctive expression phenotype, and this may have a future role in predicting genotypes, with application to clinical detection and classification of mutations.

Time- and concentration-dependent changes in gene expression induced by benzo(a)pyrene in two human cell lines, MCF-7 and HepG2.

BACKGROUND: The multi-step process of carcinogenesis can be more fully understood by characterizing gene expression changes induced in cells by carcinogens. In this study, expression microarrays were used to monitor the activity of 18,224 cDNA clones in MCF-7 and HepG2 cells exposed to the carcinogen benzo(a)pyrene (BaP) or its non-carcinogenic isomer benzo(e)pyrene (BeP). Time and concentration gene expression effects of BaP exposure have been assessed and linked to other measures of cellular stress to aid in the identification of novel genes/pathways involved in the cellular response to genotoxic carcinogens. RESULTS: BaP (0.25-5.0 muM; 6-48 h exposure) modulated 202 clones in MCF-7 cells and 127 in HepG2 cells, including 27 that were altered in both. In contrast, BeP did not induce consistent gene expression changes at the same concentrations. Significant time- and concentration-dependent responses to BaP were seen in both cell lines. Expression changes observed in both cell lines included genes involved in xenobiotic metabolism (e.g., CYP1B1, NQO1, MGST1, AKR1C1, AKR1C3,CPM), cell cycle regulation (e.g., CDKN1A), apoptosis/anti-apoptosis (e.g., BAX, IER3), chromatin assembly (e.g., histone genes), and oxidative stress response (e.g., TXNRD1). RTqPCR was used to validate microarray data. Phenotypic anchoring of the expression data to DNA adduct levels detected by 32P-postlabelling, cell cycle data and p53 protein expression identified a number of genes that are linked to these biological outcomes, thereby strengthening the identification of target genes. The overall response to BaP consisted of up-regulation of tumour suppressor genes and down-regulation of oncogenes promoting cell cycle arrest and apoptosis. Anti-apoptotic signalling that may increase cell survival and promote tumourigenesis was also evident. CONCLUSION: This study has further characterised the gene expression response of human cells after genotoxic insult, induced after exposure to concentrations of BaP that result in minimal cytotoxicity. We have demonstrated that investigating the time and concentration effect of a carcinogen on gene expression related to other biological end-points gives greater insight into cellular responses to such compounds and strengthens the identification of target genes.

Genome-wide screening for complete genetic loss in prostate cancer by comparative hybridization onto cDNA microarrays.

We demonstrate that comparative genomic hybridization (CGH) onto cDNA microarrays may be used to carry out genome-wide screens for regions of genetic loss, including homozygous (complete) deletions that may represent the possible location of tumour suppressor genes in human cancer. Screening of the prostate cancer cell lines LNCaP, PC3 and DU145 allowed the mapping of specific regions where genome copy number appeared altered and led to the identification of two novel regions of complete loss at 17q21.31 (500 kb spanning STAT3) and at 10q23.1 (50-350 kb spanning SFTPA2) in the PC3 cell line.

Gene expression profiling after radiation-induced DNA damage is strongly predictive of BRCA1 mutation carrier status.

PURPOSE: The impact of the presence of a germ-line BRCA1 mutation on gene expression in normal breast fibroblasts after radiation-induced DNA damage has been investigated. EXPERIMENTAL DESIGN: High-density cDNA microarray technology was used to identify differential responses to DNA damage in fibroblasts from nine heterozygous BRCA1 mutation carriers compared with five control samples without personal or family history of any cancer. Fibroblast cultures were irradiated, and their expression profile was compared using intensity ratios of the cDNA microarrays representing 5603 IMAGE clones. RESULTS: Class comparison and class prediction analysis has shown that BRCA1 mutation carriers can be distinguished from controls with high probability (approximately 85%). Significance analysis of microarrays and the support vector machine classifier identified gene sets that discriminate the samples according to their mutation status. These include genes already known to interact with BRCA1 such as CDKN1B, ATR, and RAD51. CONCLUSIONS: The results of this initial study suggest that normal cells from heterozygous BRCA1 mutation carriers display a different gene expression profile from controls in response to DNA damage. Adaptations of this pilot result to other cell types could result in the development of a functional assay for BRCA1 mutation status.

Serum angiogenic profile of patients with glioblastoma identifies distinct tumor subtypes and shows that TIMP-1 is a prognostic factor.

Angiogenesis plays a key role in glioblastoma biology and antiangiogenic agents are under clinical investigation with promising results. However, the angiogenic profiles of patients with glioblastoma and their clinical significance are not well understood. Here we characterize the serum angiogenic profile of patients with glioblastoma, and examine the prognostic significance of individual angiogenic factors. Serum samples from 36 patients with glioblastoma were collected on admission and simultaneously assayed for 48 angiogenic factors using protein microarrays. The data were analyzed using hierarchical cluster analysis. Vessel morphology was assessed histologically after immunostaining for the pan-endothelial marker CD31. Tumor samples were also immunostained for tissue inhibitor of metalloproteinase-1 (TIMP-1). Cluster analysis of the serum angiogenic profiles revealed 2 distinct subtypes of glioblastoma. The 2 subtypes had markedly different tumor microvessel densities. A low serum level of TIMP-1 was associated with significantly longer survival independent of patient age, performance status, or treatment. The serum angiogenic profile in patients with glioblastoma mirrors tumor biology and has prognostic value. Our data suggest the serum TIMP-1 level as an independent predictor of survival.

Correlation of overall survival with gene expression profiles in a prospective study of resectable esophageal cancer.

PURPOSE: Preoperative chemotherapy has demonstrated a survival benefit for patients with potentially resectable esophageal cancer; however, currently it is not possible to predict the benefit of this treatment for an individual patient. This prospective study was designed to correlate gene expression profiles with clinical outcome in this setting. PATIENTS AND METHODS: Eligible patients were deemed to have resectable disease after staging by computed tomography, endoscopic ultrasound, and laparoscopy as indicated and following discussion at the multidisciplinary team meeting. All patients received neoadjuvant platinum and fluoropyrimidine-based chemotherapy; and clinical data were entered prospectively onto a study-specific database. Total RNA was isolated from pretreatment tumor biopsies obtained at baseline endoscopy and analyzed using a cDNA array consisting of 22,000 cDNA clones. RESULTS: Of the patients with adequate follow-up accrued between 2002 and 2005, 35 satisfied the quality control measures for the microarray profiling. Median follow-up was 938 days. Supervised hierarchical clustering of normalized data revealed 165 significantly differentially expressed genes based on overall survival (OS; P < .01) with 2 distinct clusters: a poor outcome group: N = 17 (1 year OS 46.2%) and a good outcome group: N = 18 (1 year OS 100%). Genes identified included those previously associated with esophageal cancer and, interestingly, a group of genes encoding proteins involved in the regulation of the TOLL receptor-signaling pathway. CONCLUSION: This initial study has highlighted groups of tumors with distinct gene expression profiles based on survival and warrants further validation in a larger cohort. This approach may further our understanding of individual tumor biology and thus facilitate the development of tailored treatment.

Interlaboratory and interplatform comparison of microarray gene expression analysis of HepG2 cells exposed to benzo(a)pyrene.

Microarray technology is being used increasingly to study gene expression of biological systems on a large scale. Both interlaboratory and interplatform differences are known to contribute to variability in microarray data. In this study we have investigated data from different platforms and laboratories on the transcriptomic profile of HepG2 cells exposed to benzo(a)pyrene (BaP). RNA samples generated in two different laboratories were analyzed using both Agilent oligonucleotide microarrays and Cancer Research UK (CR-UK) cDNA microarrays. Comparability of the expression profiles was assessed at various levels including correlation and overlap between the data, clustering of the data and affected biological processes. Overlap and correlation occurred, but it was not possible to deduce whether choice of platform or interlaboratory differences contributed more to the data variation. Principal component analysis (PCA) and hierarchical clustering of the expression profiles indicated that the data were most clearly defined by duration of exposure to BaP, suggesting that laboratory and platform variability does not mask the biological effects. Real-time quantitative PCR was used to validate the two array platforms and indicated that false negatives, rather than false positives, are obtained with both systems. All together these results suggest that data from similar biological experiments analyzed on different microarray platforms can be combined to give a more complete transcriptomic profile. Each platform gives a slight variation in the BaP-gene expression response and, although it cannot be stated which is more correct, combining the two data sets is more informative than considering them individually.

AHR- and DNA-damage-mediated gene expression responses induced by benzo(a)pyrene in human cell lines.

Carcinogens induce complex transcriptional responses in cells that may hold key mechanistic information. Benzo(a)pyrene (BaP) modulation of transcription may occur through the activation of the aryl hydrocarbon receptor (AHR) or through responses to DNA damage. To characterize further the expression profiles induced by BaP in HepG2 and MCF-7 cells obtained in our previous study, they were compared to those induced by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), which activates AHR but does not bind to DNA, and anti-benzo(a)pyrene- trans-7,8-dihydrodiol-9,10-epoxide (BPDE), which binds directly to DNA but does not activate AHR. A total of 22 genes had altered expression in MCF-7 cells after both BaP and TCDD exposure, and a total of 29 genes had altered expression in HepG2 cells. In both cell lines, xenobiotic metabolism was upregulated through induction of NQO1, MGST1, and CYP1B1. A total of 78 expression changes were induced by both BaP and BPDE in MCF-7 cells, and a total of 29 expression changes were induced by both BaP and BPDE in HepG2 cells. These genes were predominantly involved in cell cycle regulation, apoptosis, and DNA repair. BaP and BPDE caused the repression of histone genes in both cell lines, suggesting that regulation of these genes is an important component of the DNA damage response. Interestingly, overlap of the BPDE and TCDD gene expression profiles was also observed. Furthermore, some genes were modulated by BaP but not by TCDD or BPDE, including induction of CRY1 and MAK, which may represent novel signaling pathways that are independent of both AHR activation and DNA damage. Promoter analysis identified candidate genes for direct transcriptional regulation by either AHR or p53. These analyses have further dissected and characterized the complex cellular response to BaP.

Identification of amplified and expressed genes in breast cancer by comparative hybridization onto microarrays of randomly selected cDNA clones.

Microarray analysis using sets of known human genes provides a powerful platform for identifying candidate oncogenes involved in DNA amplification events but suffers from the disadvantage that information can be gained only on genes that have been preselected for inclusion on the array. To address this issue, we have performed comparative genome hybridization (CGH) and expression analyses on microarrays of clones, randomly selected from a cDNA library, prepared from a cancer containing the DNA amplicon under investigation. Application of this approach to the BT474 breast carcinoma cell line, which contains amplicons at 20q13, 17q11-21, and 17q22-23, identified 50 amplified and expressed genes, including genes from these regions previously proposed as candidate oncogenes. When considered together with data from microarray expression profiles and Northern analyses, we were able to propose five genes as new candidate oncogenes where amplification in breast cancer cell lines was consistently associated with higher levels of RNA expression. These included the HB01 histone acetyl transferase gene at 17q22-23 and the TRAP100 gene, which encodes a thyroid hormone receptor-associated protein coactivator, at 17q11-21. The results demonstrate the utility of this microarray-based CGH approach in hunting for candidate oncogenes within DNA amplicons.