Search

Abstract Background Colorectal cancer (CRC) is a heterogeneous and biologically poorly understood disease. To tailor CRC treatment, it is essential to first model this heterogeneity by defining subtypes of patients with homogeneous biological and clinical characteristics and second match these subtypes to cell lines for which extensive pharmacological data is available, thus linking targeted therapies to patients most likely to respond to treatment. Methods We applied a new unsupervised, iterative approach to stratify CRC tumor samples into subtypes based on genome-wide mRNA expression data. By applying this stratification to several CRC cell line panels and integrating pharmacological response data, we generated hypotheses regarding the targeted treatment of different subtypes. Results In agreement with earlier studies, the two dominant CRC subtypes are highly correlated with a gene expression signature of epithelial-mesenchymal-transition (EMT). Notably, further dividing these two subtypes using iNMF (iterative Non-negative Matrix Factorization) revealed five subtypes that exhibit activation of specific signaling pathways, and show significant differences in clinical and molecular characteristics. Importantly, we were able to validate the stratification on independent, published datasets comprising over 1600 samples. Application of this stratification to four CRC cell line panels comprising 74 different cell lines, showed that the tumor subtypes are well represented in available CRC cell line panels. Pharmacological response data for targeted inhibitors of SRC, WNT, GSK3b, aurora kinase, PI3 kinase, and mTOR, showed significant differences in sensitivity across cell lines assigned to different subtypes. Importantly, some of these differences in sensitivity were in concordance with high expression of the targets or activation of the corresponding pathways in primary tumor samples of the same subtype. Conclusions The stratification presented here is robust, captures important features of CRC, and offers valuable insight into functional differences between CRC subtypes. By matching the identified subtypes to cell line panels that have been pharmacologically characterized, it opens up new possibilities for the development and application of targeted therapies for defined CRC patient sub-populations.

In most colorectal cancer (CRC) patients, outcome cannot be predicted because tumors with similar clinicopathological features can have differences in disease progression and treatment response. Therefore, a better understanding of the CRC biology is required to identify those patients who will benefit from chemotherapy and to find a more tailored therapy plan for other patients. Based on unsupervised classification of whole genome data from 188 stages I-IV CRC patients, a molecular classification was developed that consist of at least three major intrinsic subtypes (A-, B- and C-type). The subtypes were validated in 543 stages II and III patients and were associated with prognosis and benefit from chemotherapy. The heterogeneity of the intrinsic subtypes is largely based on three biological hallmarks of the tumor: epithelial-to-mesenchymal transition, deficiency in mismatch repair genes that result in high mutation frequency associated with microsatellite instability and cellular proliferation. A-type tumors, observed in 22% of the patients, have the best prognosis, have frequent BRAF mutations and a deficient DNA mismatch repair system. C-type patients (16%) have the worst outcome, a mesenchymal gene expression phenotype and show no benefit from adjuvant chemotherapy treatment. Both A-type and B-type tumors have a more proliferative and epithelial phenotype and B-types benefit from adjuvant chemotherapy. B-type tumors (62%) show a low overall mutation frequency consistent with the absence of DNA mismatch repair deficiency. Classification based on molecular subtypes made it possible to expand and improve CRC classification beyond standard molecular and immunohistochemical assessment and might help in the future to guide treatment in CRC patients., (© 2013 UICC.)

We have used global gene expression profiling, combined with pathway analysis tools, to identify in rats the molecular events associated with paraquat toxicity in the lung. Early (2, 8 and 18h) gene expression changes induced following intraperitoneal (i.p.) exposure to paraquat were measured in the caudal lobe of lungs using Affymetrix rat genome GeneChips (31,042 probe sets). A single high dose of paraquat dichloride (20mg/kg) was used that has been shown previously to cause in rats extensive lung fibrosis after 10 days. Hierarchical clustering of 543 paraquat-responsive genes (false discovery rate<0.05) revealed that under these conditions of exposure paraquat induces a staged transcriptional response in the rat lung that precedes the appearance of lung damage. We report here that many of the transcriptional responses to paraquat were rapid (being maximal at 2h post-dose), and that the predominant molecular functions and biological processes associated with these genes include membrane transport, oxidative stress, lung development, epithelial cell differentiation and transforming growth factor beta (TGF-beta) signalling. These data provide novel insights into the molecular pathways that lead to toxicity after exposure of the rat lung to paraquat.

We have previously used genome-wide transcript profiling to investigate the relationships between changes in gene expression and physiological alterations during the response of the immature mouse uterus to estrogens. Here we describe the identification of a functionally inter-related group of estrogen-responsive genes associated with iron homeostasis, including the iron-binding protein lactotransferrin, the ferroxidase ceruloplasmin, the iron delivery protein lipocalin 2 and the iron-exporter ferroportin. Quantitative real-time PCR revealed that the expression of these genes increases with time during the uterotrophic response, reaching maximal levels in the post-proliferative phase (between 48 and 72 h). In contrast, the heme biosynthesis genes aminolevulinic acid synthase 1 and 2 were maximally induced by estrogen at 2 and 4 h, respectively, prior to increased cell proliferation. Together, these data reveal that estrogen induces the temporally coordinated expression of iron homeostasis genes in the mouse uterus, and suggest an important role for iron metabolism during sex steroid hormone-induced uterine cell growth and differentiation.

Phosphorylation of histone H3 is implicated in transcriptional activation and chromosome condensation, but its immediate molecular function has remained obscure. By affinity chromatography of nuclear extracts against modified H3 tail peptides, we identified 14-3-3 isoforms as proteins that bind these tails in a strictly phosphorylation-dependent manner. Acetylation of lysines 9 and 14 does not impede 14-3-3 binding to serine 10-phosphorylated H3 tails. In vivo, 14-3-3 is inducibly recruited to c-fos and c-jun nucleosomes upon gene activation, concomitant with H3 phosphoacetylation. We have determined the structures of 14-3-3zeta complexed with serine 10-phosphorylated or phosphoacetylated H3 peptides. These reveal a distinct mode of 14-3-3/phosphopeptide binding and provide a structural understanding for the lack of effect of acetylation at lysines 9 and 14 on this interaction. 14-3-3 isoforms thus represent a class of proteins that mediate the effect of histone phosphorylation at inducible genes.

The recent sequencing of mammalian genomes has driven the development of genomic technologies, including microarray-based gene expression profiling, that allow simultaneous measurement of the expression levels of thousands of genes. Gene expression profiling applied to toxicology (toxicogenomics) has the potential to reveal, holistically, the molecular pathways and cellular processes that mediate the adverse responses to a toxicant. However, the initial output of a toxicogenomics experiment consists of a list of genes whose expression is altered upon toxicant exposure. In order to interpret these data in a biological context, new bioinformatic methods must be developed to place gene expression changes in the context of the underlying pathways and processes affected. One emerging approach is the application of Gene Ontology (GO) mapping and pathway analysis to gene expression profiling data. The utility of this in mechanistic toxicology will be illustrated using examples in which GO mapping of toxicogenomic data has provided novel insights into the molecular mechanisms induced by exposure to xenoestrogens.

Toxicogenomics has the potential to reveal the molecular pathways and cellular processes that mediate the adverse responses to a toxicant. However, the initial output of a toxicogenomic experiment often consists of large lists of genes whose expression is altered after toxicant exposure. To interpret gene expression changes in the context of underlying biological pathways and processes, new bioinformatics methods must be developed. We have used global gene expression profiling combined with an evaluation of Gene Ontology (GO) and pathway mapping tools as unbiased methods for identifying the molecular pathways and processes affected upon toxicant exposure. We chose to use the acute effects caused by the non-genotoxic carcinogen and peroxisome proliferator (PP) diethylhexylphthalate (DEHP) in the mouse liver as a model system. Consistent with what is known about the mode of action of DEHP, our GO analysis of transcript profiling data revealed a striking overrepresentation of genes associated with the peroxisomal cellular component, together with genes involved in carboxylic acid and lipid metabolism. Furthermore we reveal gene expression changes associated with additional biological functions, including complement activation, hemostasis, the endoplasmic reticulum overload response, and circadian rhythm. Together, these data reveal potential new pathways of PP action and shed new light on the mechanisms by which non-genotoxic carcinogens control hepatocyte hypertrophy and proliferation. We demonstrate that GO mapping can identify, in an unbiased manner, both known and novel DEHP-induced molecular changes in the mouse liver and is therefore a powerful approach for elucidating modes of toxicity based on toxicogenomic data.

A major challenge in the emerging field of toxicogenomics is to define the relationships between chemically induced changes in gene expression and alterations in conventional toxicologic parameters such as clinical chemistry and histopathology. We have explored these relationships in detail using the rodent uterotrophic assay as a model system. Gene expression levels, uterine weights, and histologic parameters were analyzed 1, 2, 4, 8, 24, 48, and 72 hr after exposure to the reference physiologic estrogen 17 beta-estradiol (E2). A multistep analysis method, involving unsupervised hierarchical clustering followed by supervised gene ontology-driven clustering, was used to define the transcriptional program associated with E2-induced uterine growth and to identify groups of genes that may drive specific histologic changes in the uterus. This revealed that uterine growth and maturation are preceded and accompanied by a complex, multistage molecular program. The program begins with the induction of genes involved in transcriptional regulation and signal transduction and is followed, sequentially, by the regulation of genes involved in protein biosynthesis, cell proliferation, and epithelial cell differentiation. Furthermore, we have identified genes with common molecular functions that may drive fluid uptake, coordinated cell division, and remodeling of luminal epithelial cells. These data define the mechanism by which an estrogen induces organ growth and tissue maturation, and demonstrate that comparison of temporal changes in gene expression and conventional toxicology end points can facilitate the phenotypic anchoring of toxicogenomic data.

Eukaryotic cells store their genetic information in the form of a highly organized nucleoprotein complex termed chromatin. The high degree of compaction of DNA within chromatin places severe constraints on proteins that require access to the DNA template to facilitate gene transcription, DNA replication, and DNA repair. As a consequence, eukaryotic cells have developed sophisticated mechanisms to allow chromatin to be rapidly decompacted locally for access by DNA-binding proteins. Once thought to play only a structural role, it now appears that chromatin plays a key regulatory role by marshalling access to the DNA template. We have reviewed the role played by chromatin in the cellular response to physiological and toxicological stimuli and described how changes in chromatin structure may in the future be used as markers of toxicity. We also review the evidence that chromatin itself is the direct target of certain toxicants and that toxicant-induced perturbations in chromatin structure may precipitate adverse effects.