Colorectal cancer intrinsic subtypes predict chemotherapy benefit, deficient mismatch repair and epithelial-to-mesenchymal transition.

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.

Induction of iron homeostasis genes during estrogen-induced uterine growth and differentiation.

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.

Meeting report: Validation of toxicogenomics-based test systems: ECVAM-ICCVAM/NICEATM considerations for regulatory use.

This is the report of the first workshop "Validation of Toxicogenomics-Based Test Systems" held 11-12 December 2003 in Ispra, Italy. The workshop was hosted by the European Centre for the Validation of Alternative Methods (ECVAM) and organized jointly by ECVAM, the U.S. Interagency Coordinating Committee on the Validation of Alternative Methods (ICCVAM), and the National Toxicology Program (NTP) Interagency Center for the Evaluation of Alternative Toxicological Methods (NICEATM). The primary aim of the workshop was for participants to discuss and define principles applicable to the validation of toxicogenomics platforms as well as validation of specific toxicologic test methods that incorporate toxicogenomics technologies. The workshop was viewed as an opportunity for initiating a dialogue between technologic experts, regulators, and the principal validation bodies and for identifying those factors to which the validation process would be applicable. It was felt that to do so now, as the technology is evolving and associated challenges are identified, would be a basis for the future validation of the technology when it reaches the appropriate stage. Because of the complexity of the issue, different aspects of the validation of toxicogenomics-based test methods were covered. The three focus areas include a) biologic validation of toxicogenomics-based test methods for regulatory decision making, b) technical and bioinformatics aspects related to validation, and c) validation issues as they relate to regulatory acceptance and use of toxicogenomics-based test methods. In this report we summarize the discussions and describe in detail the recommendations for future direction and priorities.

Identification of early molecular pathways affected by paraquat in rat lung.

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.

Gene ontology mapping as an unbiased method for identifying molecular pathways and processes affected by toxicant exposure: application to acute effects caused by the rodent non-genotoxic carcinogen diethylhexylphthalate.

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.

Mapping molecular responses to xenoestrogens through Gene Ontology and pathway analysis of toxicogenomic data.

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.

Confirmation of uterotrophic activity of 3-(4-methylbenzylidine)camphor in the immature rat.

In this study we found that the ultraviolet sunscreen component 3-(4-methylbenzylidine)camphor (4MBC) is uterotrophic in immature rats when administered by either subcutaneous injection or oral gavage. These data confirm earlier reports of uterotrophic activity for this agent when administered to immature rats in the diet or by whole-body immersion; however, they are in contrast to negative unpublished immature rat uterotrophic assay results. Data also indicate that 4MBC binds to isolated rat uterine estrogen receptors and shows activity in a human estrogen receptor yeast transactivation assay; however, we considered both of these effects equivocal. In this study, we confirmed the original observation that 4MBC was active as a mitogen to MCF-7 breast cancer cells. We evaluated and discounted the possibility that the estrogenic activity of 4MBC is related to its bulky camphor group, which is of similar molecular dimensions to that of the weak estrogen kepone. Uncertainty remains regarding the mechanism of the uterotrophic activity of 4MBC.

The need to decide if all estrogens are intrinsically similar.

We used gene expression profiling to investigate whether the molecular effects induced by estrogens of different provenance are intrinsically similar. In this article we show that the physiologic estrogen 17-beta-estradiol, the phytoestrogen genistein, and the synthetic estrogen diethylstilbestrol alter the expression of the same 179 genes in the intact immature mouse uterus under conditions where each chemical has produced an equivalent gravimetric and histologic uterotrophic effect, using the standard 3-day assay protocol. Data are also presented indicating the limitations associated with comparison of gene expression profiles for different chemicals at times before the uterotrophic effects are fully realized. We conclude that the case has yet to be made for regarding synthetic estrogens as presenting a unique human hazard compared with phytoestrogens and physiologic estrogens. Key words: diethylstilbestrol, estrogen, gene expression, genistein, microarray, phytoestrogen, toxicogenomics, uterus.

Phenotypic anchoring of gene expression changes during estrogen-induced uterine growth.

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.

Toxicogenetics: applications and opportunities.

The response to drugs and environmental chemicals varies with genotype. Some patients react well to drugs, while others may not benefit, or may even respond adversely. Individuals also experience different reactions to environmental agents, such as allergens. The sequencing of the human genome and the large-scale identification of genome polymorphisms have provided opportunities for understanding the genetic basis for individual differences in response to potential toxicants: an area of study that has come to be known as toxicogenetics. In this article, we discuss the potential applications and implications of this evolving branch of toxicology.

The role of chromatin in molecular mechanisms of toxicity.

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.

Downregulation of lactoferrin by PPARalpha ligands: role in perturbation of hepatocyte proliferation and apoptosis.

PPARalpha (peroxisome proliferator activated receptor alpha) is a transcription factor that mediates the rodent liver tumorigenic responses to peroxisome proliferators via regulation of genes that remain to be identified. Using microarray gene expression profiling of mRNA from wild type versus PPARalpha null mice, we detected a 3- to 7-fold downregulation of hepatic lactoferrin (LF) in response to the PP, diethylhexylphthalate (DEHP; 1150 mg/kg). Northern blot analyses confirmed a significant downregulation of LF mRNA by DEHP in wild type mouse liver. Since LF has been reported to repress tumor necrosis factor-alpha (TNF-alpha), LF downregulation by PPs may permit TNF-alpha levels to rise, enhancing hepatocyte survival and proliferation. To test this hypothesis, we asked if exogenous LF could prevent the perturbation of hepatocyte growth by PPs but not by TNF-alpha. In vitro, the PPs monoethylhexylphthalate (MEHP; 500 microM, the active metabolite of DEHP) and another PP, nafenopin (50 microM) or exogenous TNF-alpha (5000 U/ml) induced hepatocyte proliferation and suppressed apoptosis. LF (200 microM) blocked the growth but not the peroxisome proliferation response to PPs but could not block the growth response to TNF-alpha. Immunocytochemistry using specific antibodies to LF but also to transferrin (TF), a related gene previously shown to contain a PP response element (PPRE), demonstrated that both LF and TF are expressed in murine liver. Furthermore, both were downregulated by DEHP in both wild type and PPARalpha null mouse liver. These data suggest that the regulation of iron binding proteins by PPARalpha ligands plays a role in PP-mediated liver growth, but not in peroxisome proliferation.

Assessment of glycosylation-dependent cell adhesion molecule 1 as a correlate of allergen-stimulated lymph node activation.

Early changes in gene expression have been identified by cDNA microarray technology. Analysis of draining auricular lymph node tissue sampled at 48 h following exposure to the potent contact allergen 2,4-dinitrofluorobenzene (DNFB) provided examples of up- and down-regulated genes, including onzin and guanylate binding protein 2, and glycosylation-dependent cell adhesion molecule 1 (GlyCAM-1), respectively. Allergen-induced changes in these three genes were confirmed in dose-response and kinetic analyses using Northern blotting and/or reverse transcription-polymerase chain reaction techniques. The results confirmed that these genes are robust and relatively sensitive markers of early changes provoked in the lymph node by contact allergen. Upon further investigation, it was found that altered expression of the adhesion molecule GlyCAM-1 was not restricted to treatment with DNFB. Topical sensitization of mice to a chemically unrelated contact allergen, oxazolone, was also associated with a decrease in the expression of mRNA for GlyCAM-1. Supplementary experiments revealed that changes in expression of this gene are independent of the stimulation by chemical allergens of proliferative responses by draining lymph node cells. Taken together these data indicate that the expression of GlyCAM-1 is down-regulated rapidly following epicutaneous treatment of mice with chemical allergens, but that this reduction is associated primarily with changes in lymph node cell number, or some other aspect of lymph node activation, rather than proliferation.

Genomic analysis of stress response genes.

Mammalian cells respond to a wide range of external stimuli including growth factors, peptide hormones, cytokines, osmotic stress, heat shock, pharmacological agents and toxicants via multiple signalling pathways. Genome-wide transcript profiling simultaneously monitors the gene expression programs downstream of all signal transduction pathways and can identify novel molecular targets for stress-inducing signals. Our laboratory has combined transcript profiling of cytotoxic compounds with experimental systems in which signalling components are disrupted (e.g. small molecule protein kinase inhibitors) to reveal the contribution of specific signalling pathways to the transcriptional response to toxicant-induced stress. A complementary approach for elucidating the molecular mechanisms that regulate transcriptional responses to toxicants involves DNA sequence analysis of gene regulatory regions obtained via data mining of recently completed mammalian genome sequences. Together, these approaches reveal the molecular mechanisms used to finely tune alterations in gene expression, enabling cells to react in an appropriate manner to external stress-inducing stimuli.

Perturbation of epigenetic status by toxicants.

It is becoming increasingly apparent that toxicant-induced changes in epigenetic status, particularly DNA methylation patterns, may play a role in some mechanisms of toxicity. Here, we discuss briefly the evidence that alterations in DNA methylation accompany, and may even promote, carcinogenesis induced by non-genotoxic chemicals. We also address recent findings indicating that the availability of dietary methyl donors can modulate DNA methylation levels and precipitate adverse effects.

Subtypes of primary colorectal tumors correlate with response to targeted treatment in colorectal cell lines.

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.

Molecular basis for the recognition of phosphorylated and phosphoacetylated histone h3 by 14-3-3.

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.

A unified theory of gene expression.

The human genome has been called "the blueprint for life." This master plan is realized through the process of gene expression. Recent progress has revealed that many of the steps in the pathway from gene sequence to active protein are connected, suggesting a unified theory of gene expression.

Sex steroids, ANGELS and osteoporosis.

Osteoporosis is characterized by reduced bone density and strength. Bone mass peaks between age 30 and 40 and then declines. This can be accelerated by factors including menopause and insufficient dietary calcium. Hormone replacement therapy (HRT) is currently the standard treatment for osteoporosis. However, growing concern over potential side effects of HRT has driven a search for alternative therapies. A recent report 1 reveals a potential alternative to HRT: a gender-neutral synthetic steroid that increases bone mass and strength without affecting reproductive organs. This compound acts via a novel extranuclear sex steroid receptor signaling mechanism that has important implications for nuclear receptor biology and human health.

FACT facilitates transcription-dependent nucleosome alteration.

The FACT (facilitates chromatin transcription) complex is required for transcript elongation through nucleosomes by RNA polymerase II (Pol II) in vitro. Here, we show that FACT facilitates Pol II-driven transcription by destabilizing nucleosomal structure so that one histone H2A-H2B dimer is removed during enzyme passage. We also demonstrate that FACT possesses intrinsic histone chaperone activity and can deposit core histones onto DNA. Importantly, FACT activity requires both of its constituent subunits and is dependent on the highly acidic C terminus of its larger subunit, Spt16. These findings define the mechanism by which Pol II can transcribe through chromatin without disrupting its epigenetic status.