Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

Uncovering BRCA1-regulated signalling pathways by microarray-based expression profiling

The introduction of microarray technology to the scientific and medical communities has dramatically changed the way in which we now address basic biomedical questions. Expression profiling using microarrays facilitates an experimental approach where alterations in the transcript level of entire transcriptomes can be simultaneously assayed in response to defined stimuli. We have used microarray analysis to identify downstream transcriptional targets of the BRCA1 (Breast Cancer 1) tumour-suppressor gene as a means of defining its function. BRCA1 has been implicated in the predisposition to early onset breast and ovarian cancer and while its exact function remains to be defined, roles in DNA repair, cell-cycle control and transcriptional regulation have been implied. In the current study we have generated cell lines with tetracycline-regulated, inducible expression of BRCA1 as a tool to identify genes, which might represent important effectors of BRCA1 function. Oligonucleotide array-based expression profiling identified a number of genes that were upregulated at various times following inducible expression of BRCA1 including the DNA damage-responsive gene GADD45 (Growth Arrest after DNA Damage). Identified targets were confirmed by Northern blot analysis and their functional significance as BRCA1 targets examined.

Uncovering functionally relevant signaling pathways using microarray-based expression profiling

The introduction of microarray technology to the scientific and medical communities has fundamentally altered the way in which we now address basic biomedical questions. Microarrays technology facilitates a more complete and inclusive experimental approach where alterations in the transcript level of entire genomes can be simultaneously assayed in response to a variety of stimuli. Conceptually different approaches to the development of microarray technology have resulted in the generation of two different array formats: oligonucleotide arrays and cDNA arrays. The application of microarray and related technologies to identify specific targets of defined genes that have clearly been implicated in cancer progression requires a specific experimental approach. The objective of tiffs approach is to define changes in transcriptional profile that occur in response to modulating the expression level of the gene to be studied. The resulting altered expression profile can then be viewed as a blueprint by which that gene effects its cellular function. We have used oligonucleotide array-based expression profiling in collaboration with Affymetrix to identify downstream transcriptional targets of the BRCA1 tumor-suppressor gene as a means of defining its function. BRCA1 has been implicated in at least three functional pathways, namely, mediating the cellular response to DNA damage, as a cell cycle checkpoint protein and in the regulation of transcription. The physiological significance of these properties and their implications for the function of BRCA1 as a tumor-suppressor gene remain to be defined.

Expression of the NF-kappaB targets BCL2 and BIRC5/Survivin characterizes small B-cell and aggressive B-cell lymphomas, respectively

Nuclear factor kappa B (NF-kappaB) activation has been proposed as a cardinal feature of tumourigenesis, although the precise mechanism, frequency, relevance, and extent of NF-kappaB activation in lymphomas remain to be fully elucidated. In this study, expression profiling and tissue microarray studies of 209 and 323 non-Hodgkin's lymphomas (NHLs) respectively, including the most frequent sub-types of NHL, were employed to generate a hypothesis concerning the most common NF-kappaB targets in NHL. These analyses showed that NF-kappaB activation is a common phenomenon in NHL, resulting in the expression of distinct sets of NF-kappaB target genes, depending on the cell context. BCL2 and BIRC5/Survivin were identified as key NF-kappaB targets and their expression distinguished small and aggressive B-cell lymphomas, respectively. Interestingly, in the vast majority of B-cell lymphomas, the expression of these markers was mutually exclusive. A set of genes was identified whose expression correlates either with BIRC5/Survivin or with BCL2. BIRC5/Survivin expression, in contrast to BCL2, was associated with a signature of cell proliferation (overexpression of cell cycle control, DNA repair, and polymerase genes), which may contribute to the aggressive phenotype and poor prognosis of these lymphomas. Strikingly, mantle cell lymphoma and chronic lymphocytic leukaemia expressed highly elevated levels of BCL2 protein and mRNA, higher than that observed in reactive mantle zone cells or even in follicular lymphomas, where BCL2 expression is deregulated through the t(14;18) translocation. In parallel with this observation, BIRC5/Survivin expression was higher in Burkitt's lymphoma and diffuse large B-cell lymphoma than in non-tumoural germinal centre cells. In vitro studies confirmed that NF-kappaB activation contributes to the expression of both markers. In cell lines representing aggressive lymphomas, NF-kappaB inhibition resulted in a decrease in BIRC5/Survivin expression. Meanwhile, in chronic lymphocytic leukaemia (CLL)-derived lymphocytes, NF-kappaB inhibition resulted in a marked decrease in BCL2 expression.

Molecular pathology of prostate cancer:the key to identifying new biomarkers of disease

Microarray technology has recently accelerated the study of the molecular events involved in prostate cancer, offering the prospect of more precise prognosis and new therapeutic strategies. This review summarises current knowledge of the molecular pathology of prostate cancer. The expression and function of numerous genes have been shown to be altered in prostate cancer. Many of these genes are involved in cell cycle regulation, steroid hormone metabolism or regulation of gene expression. The mechanisms by which androgen independence arises are discussed, including cross-activation, gene amplification and point mutations of the androgen receptor. Analysis of changes in the levels of expression of large numbers of genes during prostate cancer progression have provided a better understanding of the basis of the disease, yielding new molecular markers, such as hepsin, with potential use in diagnosis and prognosis.

Wilms' tumor gene 1 (WT1) expression in childhood acute lymphoblastic leukemia:a wide range of WT1 expression levels, its impact on prognosis and minimal residual disease monitoring

Wilms' tumor gene 1 (WT1) is overexpressed in the majority (70-90%) of acute leukemias and has been identified as an independent adverse prognostic factor, a convenient minimal residual disease (MRD) marker and potential therapeutic target in acute leukemia. We examined WT1 expression patterns in childhood acute lymphoblastic leukemia (ALL), where its clinical implication remains unclear. Using a real-time quantitative PCR designed according to Europe Against Cancer Program recommendations, we evaluated WT1 expression in 125 consecutively enrolled patients with childhood ALL (106 BCP-ALL, 19 T-ALL) and compared it with physiologic WT1 expression in normal and regenerating bone marrow (BM). In childhood B-cell precursor (BCP)-ALL, we detected a wide range of WT1 levels (5 logs) with a median WT1 expression close to that of normal BM. WT1 expression in childhood T-ALL was significantly higher than in BCP-ALL (P<0.001). Patients with MLL-AF4 translocation showed high WT1 overexpression (P<0.01) compared to patients with other or no chromosomal aberrations. Older children (> or =10 years) expressed higher WT1 levels than children under 10 years of age (P<0.001), while there was no difference in WT1 expression in patients with peripheral blood leukocyte count (WBC) > or =50 x 10(9)/l and lower. Analysis of relapsed cases (14/125) indicated that an abnormal increase or decrease in WT1 expression was associated with a significantly increased risk of relapse (P=0.0006), and this prognostic impact of WT1 was independent of other main risk factors (P=0.0012). In summary, our study suggests that WT1 expression in childhood ALL is very variable and much lower than in AML or adult ALL. WT1, thus, will not be a useful marker for MRD detection in childhood ALL, however, it does represent a potential independent risk factor in childhood ALL. Interestingly, a proportion of childhood ALL patients express WT1 at levels below the normal physiological BM WT1 expression, and this reduced WT1 expression appears to be associated with a higher risk of relapse.

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.

A coactivator role of CARM1 in the dysregulation of β-catenin activity in colorectal cancer cell growth and gene expression

Aberrant activation of Wnt/β-catenin signaling, resulting in the expression of Wnt-regulated oncogenes, is recognized as a critical factor in the etiology of colorectal cancer. Occupancy of β-catenin at promoters of Wnt target genes drives transcription, but the mechanism of β-catenin action remains poorly understood. Here, we show that CARM1 (coactivator-associated arginine methyltransferase 1) interacts with β-catenin and positively modulates β-catenin-mediated gene expression. In colorectal cancer cells with constitutively high Wnt/β-catenin activity, depletion of CARM1 inhibits expression of endogenous Wnt/β-catenin target genes and suppresses clonal survival and anchorage-independent growth. We also identified a colorectal cancer cell line (RKO) with a low basal level of β-catenin, which is dramatically elevated by treatment with Wnt3a. Wnt3a also increased the expression of a subset of endogenous Wnt target genes, and CARM1 was required for the Wnt-induced expression of these target genes and the accompanying dimethylation of arginine 17 of histone H3. Depletion of β-catenin from RKO cells diminished the Wnt-induced occupancy of CARM1 on a Wnt target gene, indicating that CARM1 is recruited to Wnt target genes through its interaction with β-catenin and contributes to transcriptional activation by mediating events (including histone H3 methylation) that are downstream from the actions of β-catenin. Therefore, CARM1 is an important positive modulator of Wnt/β-catenin transcription and neoplastic transformation, and may thereby represent a novel target for therapeutic intervention in cancers involving aberrantly activated Wnt/β-catenin signaling.