Biclustering reveals breast cancer tumour subgroups with common clinical features and improves prediction of disease recurrence

BACKGROUND: Many studies have revealed correlations between breast tumour phenotypes, variations in gene expression, and patient survival outcomes. The molecular heterogeneity between breast tumours revealed by these studies has allowed prediction of prognosis and has underpinned stratified therapy, where groups of patients with particular tumour types receive specific treatments. The molecular tests used to predict prognosis and stratify treatment usually utilise fixed sets of genomic biomarkers, with the same biomarker sets being used to test all patients. In this paper we suggest that instead of fixed sets of genomic biomarkers, it may be more effective to use a stratified biomarker approach, where optimal biomarker sets are automatically chosen for particular patient groups, analogous to the choice of optimal treatments for groups of similar patients in stratified therapy. We illustrate the effectiveness of a biclustering approach to select optimal gene sets for determining the prognosis of specific strata of patients, based on potentially overlapping, non-discrete molecular characteristics of tumours. RESULTS: Biclustering identified tightly co-expressed gene sets in the tumours of restricted subgroups of breast cancer patients. The co-expressed genes in these biclusters were significantly enriched for particular biological annotations and gene regulatory modules associated with breast cancer biology. Tumours identified within the same bicluster were more likely to present with similar clinical features. Bicluster membership combined with clinical information could predict patient prognosis in conditional inference tree and ridge regression class prediction models. CONCLUSIONS: The increasing clinical use of genomic profiling demands identification of more effective methods to segregate patients into prognostic and treatment groups. We have shown that biclustering can be used to select optimal gene sets for determining the prognosis of specific strata of patients

MMP1 bimodal expression and differential response to inflammatory mediators is linked to promoter polymorphisms.

BACKGROUND: Identifying the functional importance of the millions of single nucleotide polymorphisms (SNPs) in the human genome is a difficult challenge. Therefore, a reverse strategy, which identifies functionally important SNPs by virtue of the bimodal abundance across the human population of the SNP-related mRNAs will be useful. Those mRNA transcripts that are expressed at two distinct abundances in proportion to SNP allele frequency may warrant further study. Matrix metalloproteinase 1 (MMP1) is important in both normal development and in numerous pathologies. Although much research has been conducted to investigate the expression of MMP1 in many different cell types and conditions, the regulation of its expression is still not fully understood. RESULTS: In this study, we used a novel but straightforward method based on agglomerative hierarchical clustering to identify bimodally expressed transcripts in human umbilical vein endothelial cell (HUVEC) microarray data from 15 individuals. We found that MMP1 mRNA abundance was bimodally distributed in un-treated HUVECs and showed a bimodal response to inflammatory mediator treatment. RT-PCR and MMP1 activity assays confirmed the bimodal regulation and DNA sequencing of 69 individuals identified an MMP1 gene promoter polymorphism that segregated precisely with the MMP1 bimodal expression. Chromatin immunoprecipitation (ChIP) experiments indicated that the transcription factors (TFs) ETS1, ETS2 and GATA3, bind to the MMP1 promoter in the region of this polymorphism and may contribute to the bimodal expression. CONCLUSIONS: We describe a simple method to identify putative bimodally expressed RNAs from transcriptome data that is effective yet easy for non-statisticians to understand and use. This method identified bimodal endothelial cell expression of MMP1, which appears to be biologically significant with implications for inflammatory disease. (271 Words).RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are

Gene expression profiling of breast tumours from New Zealand patients

AIMS: New Zealand has one of the highest rates of breast cancer incidence in the world. We investigated the gene expression profiles of breast tumours from New Zealand patients, compared them to gene expression profiles of international breast cancer cohorts and identified any associations between altered gene expression and the clinicopathological features of the tumours. METHODS: Affymetrix microarrays were used to measure the gene expression profiles of 106 breast tumours from New Zealand patients. Gene expression data from six international breast cancer cohorts were collated, and all the gene expression data were analysed using standard bioinformatic and statistical tools. RESULTS: Gene expression profiles associated with tumour ER and ERBB2 status, molecular subtype and selected gene expression signatures within the New Zealand cohort were consistent with those found in international cohorts. Significant differences in clinicopathological features such as tumour grade, tumour size and lymph node status were also observed between the New Zealand and international cohorts. CONCLUSIONS: Gene expression profiles, which are a sensitive indicator of tumour biology, showed no clear di¬fference between breast tumours from New Zealand patients and those from non-New Zealand patients. This suggests that other factors may contribute to the high and increasing breast cancer incidence in New Zealand compared to international populations

Positive regulation of c-Myc by cohesin is direct, and evolutionarily conserved

AbstractContact between sister chromatids from S phase to anaphase depends on cohesin, a large multi-subunit protein complex. Mutations in sister chromatid cohesion proteins underlie the human developmental condition, Cornelia de Lange syndrome. Roles for cohesin in regulating gene expression, sometimes in combination with CCCTC-binding factor (CTCF), have emerged. We analyzed zebrafish embryos null for cohesin subunit rad21 using microarrays to determine global effects of cohesin on gene expression during embryogenesis. This identified Rad21-associated gene networks that included myca (zebrafish c-myc), p53 and mdm2. In zebrafish, cohesin binds to the transcription start sites of p53 and mdm2, and depletion of either Rad21 or CTCF increased their transcription. In contrast, myca expression was strongly downregulated upon loss of Rad21 while depletion of CTCF had little effect. Depletion of Rad21 or the cohesin-loading factor Nipped-B in Drosophila cells also reduced expression of myc and Myc target genes. Cohesin bound the transcription start site plus an upstream predicted CTCF binding site at zebrafish myca. Binding and positive regulation of the c-Myc gene by cohesin is conserved through evolution, indicating that this regulation is likely to be direct. The exact mechanism of regulation is unknown, but local changes in histone modification associated with transcription repression at the myca gene were observed in rad21 mutants

PDX1 DNA methylation distinguishes two subtypes of pancreatic neuroendocrine neoplasms with a different prognosis

DNA methylation is a crucial epigenetic mechanism for gene expression regulation and cell differentiation. Furthermore, it was found to play a major role in multiple pathological processes, including cancer. In pancreatic neuroendocrine neoplasms (PNENs), epigenetic deregulation is also considered to be of significance, as the most frequently mutated genes have an important function in epigenetic regulation. However, the exact changes in DNA methylation between PNENs and the endocrine cells of the pancreas, their likely cell-of-origin, remain largely unknown. Recently, two subtypes of PNENs have been described which were linked to cell-of-origin and have a different prognosis. A difference in the expression of the transcription factor PDX1 was one of the key molecular differences. In this study, we performed an exploratory genome-wide DNA methylation analysis using Infinium Methylation EPIC arrays (Illumina) on 26 PNENs and pancreatic islets of five healthy donors. In addition, the methylation profile of the PDX1 region was used to perform subtyping in a global cohort of 83 PNEN, 2 healthy alpha cell and 3 healthy beta cell samples. In our exploratory analysis, we identified 26,759 differentially methylated CpGs and 79 differentially methylated regions. The gene set enrichment analysis highlighted several interesting pathways targeted by altered DNA methylation, including MAPK, platelet-related and immune system-related pathways. Using the PDX1 methylation in 83 PNEN, 2 healthy alpha cell and 3 healthy beta cell samples, two subtypes were identified, subtypes A and B, which were similar to alpha and beta cells, respectively. These subtypes had different clinicopathological characteristics, a different pattern of chromosomal alterations and a different prognosis, with subtype A having a significantly worse prognosis compared with subtype B (HR 0.22 [95% CI: 0.051–0.95], p = 0.043). Hence, this study demonstrates that several cancer-related pathways are differently methylated between PNENs and normal islet cells. In addition, we validated the use of the PDX1 methylation status for the subtyping of PNENs and its prognostic importance

A Gene Expression Signature of Invasive Potential in Metastatic Melanoma Cells

BACKGROUND: We are investigating the molecular basis of melanoma by defining genomic characteristics that correlate with tumour phenotype in a novel panel of metastatic melanoma cell lines. The aim of this study is to identify new prognostic markers and therapeutic targets that might aid clinical cancer diagnosis and management. PRINCIPAL FINDINGS: Global transcript profiling identified a signature featuring decreased expression of developmental and lineage specification genes including MITF, EDNRB, DCT, and TYR, and increased expression of genes involved in interaction with the extracellular environment, such as PLAUR, VCAN, and HIF1a. Migration assays showed that the gene signature correlated with the invasive potential of the cell lines, and external validation by using publicly available data indicated that tumours with the invasive gene signature were less melanocytic and may be more aggressive. The invasion signature could be detected in both primary and metastatic tumours suggesting that gene expression conferring increased invasive potential in melanoma may occur independently of tumour stage. CONCLUSIONS: Our data supports the hypothesis that differential developmental gene expression may drive invasive potential in metastatic melanoma, and that melanoma heterogeneity may be explained by the differing capacity of melanoma cells to both withstand decreased expression of lineage specification genes and to respond to the tumour microenvironment. The invasion signature may provide new possibilities for predicting which primary tumours are more likely to metastasize, and which metastatic tumours might show a more aggressive clinical course

Gene network inference and visualization tools for biologists: application to new human transcriptome datasets.

Gene regulatory networks inferred from RNA abundance data have generated significant interest, but despite this, gene network approaches are used infrequently and often require input from bioinformaticians. We have assembled a suite of tools for analysing regulatory networks, and we illustrate their use with microarray datasets generated in human endothelial cells. We infer a range of regulatory networks, and based on this analysis discuss the strengths and limitations of network inference from RNA abundance data. We welcome contact from researchers interested in using our inference and visualization tools to answer biological questions

A Bayesian Search for Transcriptional Motifs

Identifying transcription factor (TF) binding sites (TFBSs) is an important step towards understanding transcriptional regulation. A common approach is to use gaplessly aligned, experimentally supported TFBSs for a particular TF, and algorithmically search for more occurrences of the same TFBSs. The largest publicly available databases of TF binding specificities contain models which are represented as position weight matrices (PWM). There are other methods using more sophisticated representations, but these have more limited databases, or aren't publicly available. Therefore, this paper focuses on methods that search using one PWM per TF. An algorithm, MATCHTM, for identifying TFBSs corresponding to a particular PWM is available, but is not based on a rigorous statistical model of TF binding, making it difficult to interpret or adjust the parameters and output of the algorithm. Furthermore, there is no public description of the algorithm sufficient to exactly reproduce it. Another algorithm, MAST, computes a p-value for the presence of a TFBS using true probabilities of finding each base at each offset from that position. We developed a statistical model, BaSeTraM, for the binding of TFs to TFBSs, taking into account random variation in the base present at each position within a TFBS. Treating the counts in the matrices and the sequences of sites as random variables, we combine this TFBS composition model with a background model to obtain a Bayesian classifier. We implemented our classifier in a package (SBaSeTraM). We tested SBaSeTraM against a MATCHTM implementation by searching all probes used in an experimental Saccharomyces cerevisiae TF binding dataset, and comparing our predictions to the data. We found no statistically significant differences in sensitivity between the algorithms (at fixed selectivity), indicating that SBaSeTraM's performance is at least comparable to the leading currently available algorithm. Our software is freely available at: http://wiki.github.com/A1kmm/sbasetram/building-the-tools