Predicting Cell Cycle Regulated Genes by Causal Interactions

The fundamental difference between classic and modern biology is that technological innovations allow to generate high-throughput data to get insights into molecular interactions on a genomic scale. These high-throughput data can be used to infer gene networks, e. g., the transcriptional regulatory or signaling network, representing a blue print of the current dynamical state of the cellular system. However, gene networks do not provide direct answers to biological questions, instead, they need to be analyzed to reveal functional information of molecular working mechanisms. In this paper we propose a new approach to analyze the transcriptional regulatory network of yeast to predict cell cycle regulated genes. The novelty of our approach is that, in contrast to all other approaches aiming to predict cell cycle regulated genes, we do not use time series data but base our analysis on the prior information of causal interactions among genes. The major purpose of the present paper is to predict cell cycle regulated genes in S. cerevisiae. Our analysis is based on the transcriptional regulatory network, representing causal interactions between genes, and a list of known periodic genes. No further data are used. Our approach utilizes the causal membership of genes and the hierarchical organization of the transcriptional regulatory network leading to two groups of periodic genes with a well defined direction of information flow. We predict genes as periodic if they appear on unique shortest paths connecting two periodic genes from different hierarchy levels. Our results demonstrate that a classical problem as the prediction of cell cycle regulated genes can be seen in a new light if the concept of a causal membership of a gene is applied consequently. This also shows that there is a wealth of information buried in the transcriptional regulatory network whose unraveling may require more elaborate concepts than it might seem at first.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

Identification of 5-Fluorouracil-inducible Target Genes Using cDNA Microarray Profiling

The fluoropyrimidine 5-Fluorouracil (5-FU) is widely used in the treatment of cancer. To identify novel downstream mediators of tumor cell response to 5-FU, we used DNA microarray technology to identify genes that are transcriptionally activated by 5-FU treatment in the MCF-7 breast cancer cell line. Of 2400 genes analyzed, 619 were up-regulated by >3-fold. Highly up-regulated genes (>6-fold) with signal intensities of >3000 were analyzed by Northern blot. Genes that were consistently found to be up-regulated were spermine/spermidine acetyl transferase (SSAT), annexin II, thymosin-beta-10, chaperonin-10, and MAT-8. Treatment of MCF-7 cells with the antifolate tomudex and DNA-damaging agent oxaliplatin also resulted in up-regulation of each of these targets. The 5-FU-induced activation of MAT-8, thymosin-beta-10, and chaperonin-10 was abrogated by inactivation of p53 in MCF-7 cells, whereas induction of SSAT and annexin II was significantly reduced in the absence of p53. Moreover, each of these genes contained more than one potential p53-binding site, suggesting that p53 may play an important regulatory role in 5-FU-induced expression of these genes. In addition, we found that basal expression levels of SSAT, annexin II, thymosin beta-10, and chaperonin-10 were increased (by approximately 2-3-fold), and MAT-8 expression dramatically increased (by approximately 10-fold) in a 5-FU-resistant colorectal cancer cell line (H630-R10) compared with the parental H630 cell line, suggesting these genes may be useful biomarkers of resistance. These results demonstrate the potential of DNA microarrays to identify novel genes involved in mediating the response of tumor cells to chemotherapy.

BRCA1 regulates the interferon gamma-mediated apoptotic response

BRCA1 is a tumor suppressor gene implicated in transcriptional regulation. We have generated cell lines with inducible expression of BRCA1 as a tool to identify downstream targets that may be important mediators of BRCA1 function. Oligonucleotide array-based expression profiling identified 11 previously described interferon regulated genes that were up-regulated following inducible expression of BRCA1. Northern blot analysis revealed that a subset of the identified targets including IRF-7, MxA, and ISG-54 were synergistically up-regulated by BRCA1 in the presence of interferon gamma (IFN-gamma) but not interferons alpha or beta. Importantly, IFN-gamma-mediated induction of IRF-7 and MxA was attenuated in the BRCA1 mutant cell line HCC1937, an effect that was rescued following reconstitution of exogenous wild type BRCA1 in these cells. Furthermore, reconstituted BRCA1 sensitized HCC1937 cells to IFN-gamma-induced apoptotic cell death. This study identifies BRCA1 as a component of the IFN-gamma-regulated signaling pathway and suggests that BRCA1 may play a role in the regulation of IFN-gamma-mediated apoptosis.

SGK1 activity in Na+ absorbing airway epithelial cells monitored by assaying NDRG1-Thr(346/356/366) phosphorylation

Studies of HeLa cells and serum- and glucocorticoid-regulated kinase 1 (SGK1) knockout mice identified threonine residues in the n-myc downstream-regulated gene 1 protein (NDRG1-Thr(346/356/366)) that are phosphorylated by SGK1 but not by related kinases (Murray et al., Biochem J 385:1-12, 2005). We have, therefore, monitored the phosphorylation of NDRG1-Thr(346/356/366) in order to explore the changes in SGK1 activity associated with the induction and regulation of the glucocorticoid-dependent Na+ conductance (G (Na)) in human airway epithelial cells. Transient expression of active (SGK1-S422D) and inactive (SGK1-K127A) SGK1 mutants confirmed that activating SGK1 stimulates NDRG1-Thr(346/356/366) phosphorylation. Although G (Na) is negligible in hormone-deprived cells, these cells displayed basal SGK1 activity that was sensitive to LY294002, an inhibitor of 3-phosphatidylinositol phosphate kinase (PI3K). Dexamethasone (0.2 mu M) acutely activated SGK1 and the peak of this response (2-3 h) coincided with the induction of G (Na), and both responses were PI3K-dependent. While these data suggest that SGK1 might mediate the rise in G (Na), transient expression of the inactive SGK1-K127A mutant did not affect the hormonal induction of G (Na) but did suppress the activation of SGK1. Dexamethasone-treated cells grown on permeable supports formed confluent epithelial sheets that generated short circuit current due to electrogenic Na+ absorption. Forskolin and insulin both stimulated this current and the response to insulin, but not forskolin, was LY294002-sensitive and associated with the activation of SGK1. While these data suggest that SGK1 is involved in the control of G (Na), its role may be minor, which could explain why sgk1 knockout has different effects upon different tissues.

Mechanisms underlying the metabolic actions of galegine that contribute to weight loss in mice

Background and purpose: Galegine and guanidine, originally isolated from Galega officinalis, led to the development of the biguanides. The weight-reducing effects of galegine have not previously been studied and the present investigation was undertaken to determine its mechanism(s) of action.

Experimental approach: Body weight and food intake were examined in mice. Glucose uptake and acetyl-CoA carboxylase activity were studied in 3T3-L1 adipocytes and L6 myotubes and AMP activated protein kinase (AMPK) activity was examined in cell lines. The gene expression of some enzymes involved in fat metabolism was examined in 3T3-L1 adipocytes.

Key results: Galegine administered in the diet reduced body weight in mice. Pair-feeding indicated that at least part of this effect was independent of reduced food intake. In 3T3-L1 adipocytes and L6 myotubes, galegine (50 µm-3 mm) stimulated glucose uptake. Galegine (1–300 µm) also reduced isoprenaline-mediated lipolysis in 3T3-L1 adipocytes and inhibited acetyl-CoA carboxylase activity in 3T3-L1 adipocytes and L6 myotubes. Galegine (500 µm) down-regulated genes concerned with fatty acid synthesis, including fatty acid synthase and its upstream regulator SREBP. Galegine (10 µm and above) produced a concentration-dependent activation of AMP activated protein kinase (AMPK) in H4IIE rat hepatoma, HEK293 human kidney cells, 3T3-L1 adipocytes and L6 myotubes.

Conclusions and implications: Activation of AMPK can explain many of the effects of galegine, including enhanced glucose uptake and inhibition of acetyl-CoA carboxylase. Inhibition of acetyl-CoA carboxylase both inhibits fatty acid synthesis and stimulates fatty acid oxidation, and this may to contribute to the in vivo effect of galegine on body weight.

Connective tissue growth factor antagonizes transforming growth factor-beta 1/Smad signalling in renal mesangial cells

The critical involvement of TGF-beta 1 (transforming growth factor-beta 1) in DN (diabetic nephropathy) is well established. However, the role of CTGF (connective tissue growth factor) in regulating the complex interplay of TGF-beta 1 signalling networks is poorly understood. The purpose of the present study was to investigate co-operative signalling between CTGF and TGF-beta 1 and its physiological significance. CTGF was determined to bind directly to the T beta RIII (TGF-beta type III receptor) and antagonize TGF-beta 1-induced Smad phosphorylation and transcriptional responses via its N-terminal half. Furthermore, TGF-beta 1 binding to its receptor was inhibited by CTGF. A consequent shift towards non-canonical TGF-beta 1 signalling and expression of a unique profile of differentially regulated genes was observed in CTGF/TGF-beta 1-treated mesangial cells. Decreased levels of Smad2/3 phosphorylation were evident in STZ (streptozotocin)-induced diabetic mice, concomitant with increased levels of CTGF Knockdown of T beta RIII restored TGF-beta 1-mediated Smad signalling and cell contractility, suggesting that T beta RIII is key for CTGF-mediated regulation of TGF-beta 1. Comparison of gene expression profiles from CTGF/TGF-beta 1-treated mesangial cells and human renal biopsy material with histological diagnosis of DN revealed significant correlation among gene clusters. In summary, mesangial cell responses to TGF-beta 1 are regulated by cross-talk with CTGF, emphasizing the potential utility of targeting CTGF in DN.

The androgen receptor induces a distinct transcriptional program in castration-resistant prostate cancer in man

The androgen receptor (AR) regulates prostate cell growth in man, and prostate cancer is the commonest cancer in men in the UK. We present a comprehensive analysis of AR binding sites in human prostate cancer tissues, including castrate-resistant prostate cancer (CRPC). We identified thousands of AR binding sites in CRPC tissue, most of which were not identified in PC cell lines. Many adjacent genes showed AR regulation in xenografts but not in cultured LNCaPs, demonstrating an in-vivo-restricted set of AR-regulated genes. Functional studies support a model of altered signaling in vivo that directs AR binding. We identified a 16 gene signature that outperformed a larger in-vitro-derived signature in clinical data sets, showing the importance of persistent AR signaling in CRPC.

Global identification of androgen response elements

Chromatin immunoprecipitation (ChIP) is an invaluable tool in the study of transcriptional regulation. ChIP methods require both a priori knowledge of the transcriptional regulators which are important for a given biological system and high-quality specific antibodies for these targets. The androgen receptor (AR) is known to play essential roles in male sexual development, in prostate cancer and in the function of many other AR-expressing cell types (e.g. neurons and myocytes). As a ligand-activated transcription factor the AR also represents an endogenous, inducible system to study transcriptional biology. Therefore, ChIP studies of the AR can make use of treatment contrast experiments to define its transcriptional targets. To date several studies have mapped AR binding sites using ChIP in combination with genome tiling microarrays (ChIP-chip) or direct sequencing (ChIP-seq), mainly in prostate cancer cell lines and with varying degrees of genomic coverage. These studies have provided new insights into the DNA sequences to which the AR can bind, identified AR cooperating transcription factors, mapped thousands of potential AR regulated genes and provided insights into the biological processes regulated by the AR. However, further ChIP studies will be required to fully characterise the dynamics of the AR-regulated transcriptional programme, to map the occupancy of different AR transcriptional complexes which result in different transcriptional output and to delineate the transcriptional networks downstream of the AR.

A chromatin-independent role of Polycomb-like 1 to stabilize p53 and promote cellular quiescence

Polycomb-like proteins 1-3 (PCL1-3) are substoichiometric components of the Polycomb-repressive complex 2 (PRC2) that are essential for association of the complex with chromatin. However, it remains unclear why three proteins with such apparent functional redundancy exist in mammals. Here we characterize their divergent roles in both positively and negatively regulating cellular proliferation. We show that while PCL2 and PCL3 are E2F-regulated genes expressed in proliferating cells, PCL1 is a p53 target gene predominantly expressed in quiescent cells. Ectopic expression of any PCL protein recruits PRC2 to repress the INK4A gene; however, only PCL2 and PCL3 confer an INK4A-dependent proliferative advantage. Remarkably, PCL1 has evolved a PRC2- and chromatin-independent function to negatively regulate proliferation. We show that PCL1 binds to and stabilizes p53 to induce cellular quiescence. Moreover, depletion of PCL1 phenocopies the defects in maintaining cellular quiescence associated with p53 loss. This newly evolved function is achieved by the binding of the PCL1 N-terminal PHD domain to the C-terminal domain of p53 through two unique serine residues, which were acquired during recent vertebrate evolution. This study illustrates the functional bifurcation of PCL proteins, which act in both a chromatin-dependent and a chromatin-independent manner to regulate the INK4A and p53 pathways.

The Early Effects of Rapid Androgen Deprivation on Human Prostate Cancer

The androgen receptor (AR) is the dominant growth factor in prostate cancer (PCa). Therefore, understanding how ARs regulate the human transcriptome is of paramount importance. The early effects of castration on human PCa have not previously been studied 27 patients medically castrated with degarelix 7 d before radical prostatectomy. We used mass spectrometry, immunohistochemistry, and gene expression array (validated by reverse transcription-polymerase chain reaction) to compare resected tumour with matched, controlled, untreated PCa tissue. All patients had levels of serum androgen, with reduced levels of intraprostatic androgen at prostatectomy. We observed differential expression of known androgen-regulated genes (TMPRSS2, KLK3, CAMKK2, FKBP5). We identified 749 genes downregulated and 908 genes upregulated following castration. AR regulation of α-methylacyl-CoA racemase expression and three other genes (FAM129A, RAB27A, and KIAA0101) was confirmed. Upregulation of oestrogen receptor 1 (ESR1) expression was observed in malignant epithelia and was associated with differential expression of ESR1-regulated genes and correlated with proliferation (Ki-67 expression).

PATIENT SUMMARY: This first-in-man study defines the rapid gene expression changes taking place in prostate cancer (PCa) following castration. Expression levels of the genes that the androgen receptor regulates are predictive of treatment outcome. Upregulation of oestrogen receptor 1 is a mechanism by which PCa cells may survive despite castration.

Transcriptional profiling of the oesophageal gland region of male worms of Schistosoma mansoni

The intestinal tract of schistosomes opens at the mouth and leads into the foregut or oesophageal region that is lined with syncytium continuous with the apical cytoplasm of the tegument. The oesophagus is surrounded by a specialised gland, the oesophageal gland. This gland releases materials into the lumen of the oesophagus and the region is thought to initiate the lysis of erythrocytes and neutralisation of immune effectors of the host. The oesophageal region is present in the early invasive schistosomulum, a stage potentially targetable by anti-schistosome vaccines. We used a 44k oligonucleotide microarray to identify highly up-regulated genes in microdissected frozen sections of the oesophageal gland of male worms of S. mansoni. We show that 122 genes were up-regulated 2-fold or higher in the oesophageal gland compared with a whole male worm tissue control. The enriched genes included several associated with lipid metabolism and transmembrane transport as well as some micro-exon genes. Since the oesophageal gland is important in the initiation of digestion and the fact that it develops early after invasion of the mammalian host, further study of selected highly up-regulated functionally important genes in this tissue may reveal new anti-schistosome intervention targets for schistosomiasis control.

CTCF modulates Estrogen Receptor function through specific chromatin and nuclear matrix interactions

Enhancer regions and transcription start sites of estrogen-target regulated genes are connected by means of Estrogen Receptor long-range chromatin interactions. Yet, the complete molecular mechanisms controlling the transcriptional output of engaged enhancers and subsequent activation of coding genes remain elusive. Here, we report that CTCF binding to enhancer RNAs is enriched when breast cancer cells are stimulated with estrogen. CTCF binding to enhancer regions results in modulation of estrogen-induced gene transcription by preventing Estrogen Receptor chromatin binding and by hindering the formation of additional enhancer-promoter ER looping. Furthermore, the depletion of CTCF facilitates the expression of target genes associated with cell division and increases the rate of breast cancer cell proliferation. We have also uncovered a genomic network connecting loci enriched in cell cycle regulator genes to nuclear lamina that mediates the CTCF function. The nuclear lamina and chromatin interactions are regulated by estrogen-ER. We have observed that the chromatin loops formed when cells are treated with estrogen establish contacts with the nuclear lamina. Once there, the portion of CTCF associated with the nuclear lamina interacts with enhancer regions, limiting the formation of ER loops and the induction of genes present in the loop. Collectively, our results reveal an important, unanticipated interplay between CTCF and nuclear lamina to control the transcription of ER target genes, which has great implications in the rate of growth of breast cancer cells.

Glucocorticoid receptor and Klf4 co-regulate anti-inflammatory genes in keratinocytes

The glucocorticoid (GC) receptor (GR) and Kruppel-like factor Klf4 are transcription factors that play major roles in skin homeostasis. However, whether these transcription factors cooperate in binding genomic regulatory regions in epidermal keratinocytes was not known. Here, we show that in dexamethasone-treated keratinocytes GR and Klf4 are recruited to genomic regions containing adjacent GR and KLF binding motifs to control transcription of the anti-inflammatory genes Tsc22d3 and Zfp36. GR- and Klf4 loss of function experiments showed total GR but partial Klf4 requirement for full gene induction in response to dexamethasone. In wild type keratinocytes induced to differentiate, GR and Klf4 protein expression increased concomitant with Tsc22d3 and Zfp36 up-regulation. In contrast, GR-deficient cells failed to differentiate or fully induce Klf4, Tsc22d3 and Zfp36 correlating with increased expression of the epithelium-specific Trp63, a known transcriptional repressor of Klf4. The identified transcriptional cooperation between GR and Klf4 may determine cell-type specific regulation and have implications for developing therapies for skin diseases.