Aristotelis Tsirigos

Dynamic changes in the human methylome during differentiation

DNA methylation is a critical epigenetic regulator in mammalian development. Here, we present a whole-genome comparative view of DNA methylation using bisulfite sequencing of three cultured cell types representing progressive stages of differentiation: human embryonic stem cells (hESCs), a fibroblastic differentiated derivative of the hESCs, and neonatal fibroblasts. As a reference, we compared our maps with a methylome map of a fully differentiated adult cell type, mature peripheral blood mononuclear cells (monocytes). We observed many notable common and cell-type-specific features among all cell types. Promoter hypomethylation (both CG and CA) and higher levels of gene body methylation were positively correlated with transcription in all cell types. Exons were more highly methylated than introns, and sharp transitions of methylation occurred at exon-intron boundaries, suggesting a role for differential methylation in transcript splicing. Developmental stage was reflected in both the level of global methylation and extent of non-CpG methylation, with hESC highest, fibroblasts intermediate, and monocytes lowest. Differentiation-associated differential methylation profiles were observed for developmentally regulated genes, including the HOX clusters, other homeobox transcription factors, and pluripotence-associated genes such as POU5F1, TCF3, and KLF4. Our results highlight the value of high-resolution methylation maps, in conjunction with other systems-level analyses, for investigation of previously undetectable developmental regulatory mechanisms.

Accurate phylogenetic classification of variable-length DNA fragments.

Metagenome studies have retrieved vast amounts of sequence data from a variety of environments leading to new discoveries and insights into the uncultured microbial world. Except for very simple communities, the encountered diversity has made fragment assembly and the subsequent analysis a challenging problem. A taxonomic characterization of metagenomic fragments is required for a deeper understanding of shotgun-sequenced microbial communities, but success has mostly been limited to sequences containing phylogenetic marker genes. Here we present PhyloPythia, a composition-based classifier that combines higher-level generic clades from a set of 340 completed genomes with sample-derived population models. Extensive analyses on synthetic and real metagenome data sets showed that PhyloPythia allows the accurate classification of most sequence fragments across all considered taxonomic ranks, even for unknown organisms. The method requires no more than 100 kb of training sequence for the creation of accurate models of sample-specific populations and can assign fragments ≥1 kb with high specificity.

Genetic inactivation of the polycomb repressive complex 2 in T cell acute lymphoblastic leukemia

T-cell acute lymphoblastic leukemia (T-ALL) is an immature hematopoietic malignancy driven mainly by oncogenic activation of NOTCH1 signaling1. In this study we report the presence of loss-of-function mutations and deletions of EZH2 and SUZ12 genes, encoding critical components of the Polycomb Repressive Complex 2 (PRC2) complex2,3, in 25% of T-ALLs. To further study the role of the PRC2 complex in T-ALL, we used NOTCH1-induced animal models of the disease, as well as human T-ALL samples, and combined locus-specific and global analysis of NOTCH1-driven epigenetic changes. These studies demonstrated that activation of NOTCH1 specifically induces loss of the repressive mark lysine-27 tri-methylation of histone 3 (H3K27me3)4 by antagonizing the activity of the Polycomb Repressive Complex 2 (PRC2) complex. These studies demonstrate a tumor suppressor role for the PRC2 complex in human leukemia and suggest a hitherto unrecognized dynamic interplay between oncogenic NOTCH1 and PRC2 function for the regulation of gene expression and cell transformation.

Contrasting roles of histone 3 lysine 27 demethylases in acute lymphoblastic leukaemia.

T-cell acute lymphoblastic leukaemia (T-ALL) is a haematological malignancy with a dismal overall prognosis, including a relapse rate of up to 25%, mainly because of the lack of non-cytotoxic targeted therapy options. Drugs that target the function of key epigenetic factors have been approved in the context of haematopoietic disorders, and mutations that affect chromatin modulators in a variety of leukaemias have recently been identified; however, ‘epigenetic’ drugs are not currently used for T-ALL treatment. Recently, we described that the polycomb repressive complex 2 (PRC2) has a tumour-suppressor role in T-ALL. Here we delineated the role of the histone 3 lysine 27 (H3K27) demethylases JMJD3 and UTX in T-ALL. We show that JMJD3 is essential for the initiation and maintenance of T-ALL, as it controls important oncogenic gene targets by modulating H3K27 methylation. By contrast, we found that UTX functions as a tumour suppressor and is frequently genetically inactivated in T-ALL. Moreover, we demonstrated that the small molecule inhibitor GSKJ4 (ref. 5) affects T-ALL growth, by targeting JMJD3 activity. These findings show that two proteins with a similar enzymatic function can have opposing roles in the context of the same disease, paving the way for treating haematopoietic malignancies with a new category of epigenetic inhibitors.

Ketones and lactate increase cancer cell “stemness”, driving recurrence, metastasis and poor clinical outcome in breast cancer: Achieving personalized medicine via metabolo-genomics

Previously, we showed that high-energy metabolites (lactate and ketones) “fuel” tumor growth and experimental metastasis in an in vivo xenograft model, most likely by driving oxidative mitochondrial metabolism in breast cancer cells. To mechanistically understand how these metabolites affect tumor cell behavior, here we used genome-wide transcriptional profiling. Briefly, human breast cancer cells (MCF7) were cultured with lactate or ketones, and then subjected to transcriptional analysis (exon-array). Interestingly, our results show that treatment with these high-energy metabolites increases the transcriptional expression of gene profiles normally associated with “stemness,” including genes upregulated in embryonic stem (ES) cells. Similarly, we observe that lactate and ketones promote the growth of bonafide ES cells, providing functional validation. The lactate- and ketone-induced “gene signatures” were able to predict poor clinical outcome (including recurrence and metastasis) in a cohort of human breast cancer patients. Taken together, our results are consistent with the idea that lactate and ketone utilization in cancer cells promotes the “cancer stem cell” phenotype, resulting in significant decreases in patient survival. One possible mechanism by which these high-energy metabolites might induce stemness is by increasing the pool of Acetyl-CoA, leading to increased histone acetylation, and elevated gene expression. Thus, our results mechanistically imply that clinical outcome in breast cancer could simply be determined by epigenetics and energy metabolism, rather than by the accumulation of specific “classical” gene mutations. We also suggest that high-risk cancer patients (identified by the lactate/ketone gene signatures) could be treated with new therapeutics that target oxidative mitochondrial metabolism, such as the anti-oxidant and “mitochondrial poison” metformin. Finally, we propose that this new approach to personalized cancer medicine be termed “Metabolo-Genomics,” which incorporates features of both 1) cell metabolism and 2) gene transcriptional profiling. Importantly, this powerful new approach directly links cancer cell metabolism with clinical outcome, and new therapeutic strategies for inhibiting the TCA cycle and mitochondrial oxidative phosphorylation in cancer cells.

Loss of stromal caveolin-1 leads to oxidative stress, mimics hypoxia and drives inflammation in the tumor microenvironment, conferring the "reverse Warburg effect": a transcriptional informatics analysis with validation.

Cav-1 (-/-) deficient stromal cells are a new genetic model for myofibroblasts and cancer-associated fibroblasts. Using an unbiased informatics analysis of the transcriptional profile of Cav-1 (-/-) deficient mesenchymal stromal cells, we have now identified many of the major signaling pathways that are activated by a loss of Cav-1, under conditions of metabolic restriction (with low glucose media). Our informatics analysis suggests that a loss of Cav-1 induces oxidative stress, which mimics a constitutive pseudo-hypoxic state, leading to 1) aerobic glycolysis and 2) inflammation in the tumor stromal microenvironment. This occurs via the activation of 2 major transcription factors, namely HIF (aerobic glycolysis) and NF-kB (inflammation) in Cav-1 (-/-) stromal fibroblastic cells. Experimentally, we show that Cav-1 deficient stromal cells may possess defective mitochondria, due to the over-production of nitric oxide (NO), resulting in the tyrosine nitration of the mitochondrial respiratory chain components (such as complex I). Elevated levels of nitro-tyrosine were observed both in Cav-1 (-/-) stromal cells, and via acute knock-down with siRNA targeting Cav-1. Finally, metabolic restriction with mitochondrial (complex I) and glycolysis inhibitors was synthetically lethal with a Cav-1 (-/-) deficiency in mice. As such, Cav-1 deficient mice show a dramatically reduced mitochondrial reserve capacity. Thus, a mitochondrial defect in Cav-1 deficient stromal cells could drive oxidative stress, leading to aerobic glycolysis, and inflammation, in the tumor microenvironment. These stromal alterations may underlie the molecular basis of the “Reverse Warburg Effect”, and could provide the key to targeted anti-cancer therapies using metabolic inhibitors. In direct support of these findings, the transcriptional profile of Cav-1 (-/-) stromal cells overlaps significantly with Alzheimer’s disease, which is characterized by oxidative stress, NO over-production (peroxynitrite formation), inflammation, hypoxia, and mitochondrial dysfunction. We conclude that Cav-1 (-/-) deficient mice are a new whole-body animal model for an activated lethal tumor micro-environment, i.e., “tumor stroma” without the tumor. Since Cav-1 (-/-) mice are also an established animal model for pro-fibrotic disease, our current results may have implications for understanding the pathogenesis of scleroderma (systemic sclerosis) and pulmonary fibrosis, which are also related to abnormal mesenchymal stem cell function.

Hyperactivation of oxidative mitochondrial metabolism in epithelial cancer cells in situ: Visualizing the therapeutic effects of metformin in tumor tissue

We have recently proposed a new mechanism for explaining energy transfer in cancer metabolism. In this scenario, cancer cells behave as metabolic parasites, by extracting nutrients from normal host cells, such as fibroblasts, via the secretion of hydrogen peroxide as the initial trigger. Oxidative stress in the tumor microenvironment then leads to autophagy-driven catabolism, mitochondrial dys-function, and aerobic glycolysis. This, in turn, produces high-energy nutrients (such as L-lactate, ketones, and glutamine) that drive the anabolic growth of tumor cells, via oxidative mitochondrial metabolism. A logical prediction of this new “parasitic” cancer model is that tumor-associated fibroblasts should show evidence of mitochondrial dys-function (mitophagy and aerobic glycolysis). In contrast, epithelial cancer cells should increase their oxidative mitochondrial capacity. To further test this hypothesis, here we subjected frozen sections from human breast tumors to a staining procedure that only detects functional mitochondria. This method detects the in situ enzymatic activity of cytochrome C oxidase (COX), also known as Complex IV. Remarkably, cancer cells show an over-abundance of COX activity, while adjacent stromal cells remain essentially negative. Adjacent normal ductal epithelial cells also show little or no COX activity, relative to epithelial cancer cells. Thus, oxidative mitochondrial activity is selectively amplified in cancer cells. Although COX activity staining has never been applied to cancer tissues, it could now be used routinely to distinguish cancer cells from normal cells, and to establish negative margins during cancer surgery. Similar results were obtained with NADH activity staining, which measures Complex I activity, and succinate dehydrogenase (SDH) activity staining, which measures Complex II activity. COX and NADH activities were blocked by electron transport inhibitors, such as Metformin. This has mechanistic and clinical implications for using Metformin as an anti-cancer drug, both for cancer therapy and chemo-prevention. We also immuno-stained human breast cancers for a series of well-established protein biomarkers of metabolism. More specifically, we now show that cancer-associated fibroblasts over-express markers of autophagy (cathepsin B), mitophagy (BNIP3L), and aerobic glycolysis (MCT4). Conversely, epithelial cancer cells show the over-expression of a mitochondrial membrane marker (TOMM20), as well as key components of Complex IV (MT-CO1) and Complex II (SDH-B). We also validated our observations using a bioinformatics approach with data from >2,000 breast cancer patients, which showed the transcriptional upregulation of mitochondrial oxidative phosphorylation (OXPHOS) in human breast tumors (p < 10-20), and a specific association with metastasis. Therefore, upregulation of OXPHOS in epithelial tumor cells is a common feature of human breast cancers. In summary, our data provide the first functional in vivo evidence that epithelial cancer cells perform enhanced mitochondrial oxidative phosphorylation, allowing them to produce high amounts of ATP. Thus, we believe that mitochondria are both the “powerhouse” and “Achilles’ heel” of cancer cells.

The autophagic tumor stroma model of cancer or "battery-operated tumor growth": A simple solution to the autophagy paradox.

The role of autophagy in tumorigenesis is controversial. Both autophagy inhibitors (chloroquine) and autophagy promoters (rapamycin) block tumorigenesis by unknown mechanism(s). This is called the “Autophagy Paradox”. We have recently reported a simple solution to this paradox. We demonstrated that epithelial cancer cells use oxidative stress to induce autophagy in the tumor microenvironment. As a consequence, the autophagic tumor stroma generates recycled nutrients that can then be used as chemical building blocks by anabolic epithelial cancer cells. This model results in a net energy transfer from the tumor stroma to epithelial cancer cells (an energy imbalance), thereby promoting tumor growth. This net energy transfer is both unilateral and vectorial, from the tumor stroma to the epithelial cancer cells, representing a true host-parasite relationship. We have termed this new paradigm “The Autophagic Tumor Stroma Model of Cancer Cell Metabolism” or “Battery-Operated Tumor Growth”. In this sense, autophagy in the tumor stroma serves as a “battery” to fuel tumor growth, progression, and metastasis, independently of angiogenesis. Using this model, the systemic induction of autophagy will prevent epithelial cancer cells from using recycled nutrients, while the systemic inhibiton of autophagy will prevent stromal cells from producing recycled nutrients—both effectively “starving” cancer cells. We discuss the idea that tumor cells could become resistant to the systemic induction of autophagy, by the up-regulation of natural endogenous autophagy inhibitors in cancer cells. Alternatively, tumor cells could also become resistant to the systemic induction of autophagy, by the genetic silencing/deletion of pro-autophagic molecules, such as Beclin1. If autophagy resistance develops in cancer cells, then the systemic inhibition of autophagy would provide a therapeutic solution to this type of drug resistance, as it would still target autophagy in the tumor stroma. As such, an anti-cancer therapy that combines the alternating use of both autophagy promoters and autophagy inhibitors would be expected to prevent the onset of drug resistance. We also discuss why anti-angiogenic therapy has been found to promote tumor recurrence, progression, and metastasis. More specifically, anti-angiogenic therapy would induce autophagy in the tumor stroma via the induction of stromal hypoxia, thereby converting a non-aggressive tumor type to a “lethal” aggressive tumor phenotype. Thus, uncoupling the metabolic parasitic relationship between cancer cells and an autophagic tumor stroma may hold great promise for anti-cancer therapy. Finally, we believe that autophagy in the tumor stroma is the local microscopic counterpart of systemic wasting (cancer-associated cachexia), which is associated with advanced and metastatic cancers. Cachexia in cancer patients is not due to decreased energy intake, but instead involves an increased basal metabolic rate and increased energy expenditures, resulting in a negative energy balance. Importantly, when tumors were surgically excised, this increased metabolic rate returned to normal levels. This view of cachexia, resulting in energy transfer to the tumor, is consistent with our hypothesis. So, cancer-associated cachexia may start locally as stromal autophagy, and then spread systemically. As such, stromal autophagy may be the requisite precursor of systemic cancer-associated cachexia.

Autophagy and senescence in cancer-associated fibroblasts metabolically supports tumor growth and metastasis via glycolysis and ketone production.

Senescent fibroblasts are known to promote tumor growth. However, the exact mechanism remains largely unknown. An important clue comes from recent studies linking autophagy with the onset of senescence. Thus, autophagy and senescence may be part of the same physiological process, known as the autophagy-senescence transition (AST). To test this hypothesis, human fibroblasts immortalized with telomerase (hTERT-BJ1) were stably transfected with autophagy genes (BNIP3, CTSB or ATG16L1). Their overexpression was sufficient to induce a constitutive autophagic phenotype, with features of mitophagy, mitochondrial dysfunction and a shift toward aerobic glycolysis, resulting in L-lactate and ketone body production. Autophagic fibroblasts also showed features of senescence, with increased p21(WAF1/CIP1), a CDK inhibitor, cellular hypertrophy and increased β-galactosidase activity. Thus, we genetically validated the existence of the autophagy-senescence transition. Importantly, autophagic-senescent fibroblasts promoted tumor growth and metastasis, when co-injected with human breast cancer cells, independently of angiogenesis. Autophagic-senescent fibroblasts stimulated mitochondrial metabolism in adjacent cancer cells, when the two cell types were co-cultured, as visualized by MitoTracker staining. In particular, autophagic ATG16L1 fibroblasts, which produced large amounts of ketone bodies (3-hydroxy-butyrate), had the strongest effects and promoted metastasis by up to 11-fold. Conversely, expression of ATG16L1 in epithelial cancer cells inhibited tumor growth, indicating that the effects of autophagy are compartment-specific. Thus, autophagic-senescent fibroblasts metabolically promote tumor growth and metastasis, by paracrine production of high-energy mitochondrial fuels. Our current studies provide genetic support for the importance of “two-compartment tumor metabolism” in driving tumor growth and metastasis via a simple energy transfer mechanism. Finally, β-galactosidase, a known lysosomal enzyme and biomarker of senescence, was localized to the tumor stroma in human breast cancer tissues, providing in vivo support for our hypothesis. Bioinformatic analysis of genome-wide transcriptional profiles from tumor stroma, isolated from human breast cancers, also validated the onset of an autophagy-senescence transition. Taken together, these studies establish a new functional link between host aging, autophagy, the tumor microenvironment and cancer metabolism.

Glycolytic cancer associated fibroblasts promote breast cancer tumor growth, without a measurable increase in angiogenesis: Evidence for stromal-epithelial metabolic coupling

Previously, we proposed a new model for understanding the Warburg effect in tumorigenesis and metastasis. In this model, the stromal fibroblasts would undergo aerobic glycolysis (a.k.a., the Warburg effect)—producing and secreting increased pyruvate/lactate that could then be used by adjacent epithelial cancer cells as “fuel” for the mitochondrial TCA cycle, oxidative phosphorylation, and ATP production. To test this model more directly, here we used a matched set of metabolically well-characterized immortalized fibroblasts that differ in a single gene. CL3 fibroblasts show a shift towards oxidative metabolism, and have an increased mitochondrial mass. In contrast, CL4 fibroblasts show a shift towards aerobic glycolysis, and have a reduced mitochondrial mass. We validated these differences in CL3 and CL4 fibroblasts by performing an unbiased proteomics analysis, showing the functional upregulation of 4 glycolytic enzymes, namely ENO1, ALDOA, LDHA and TPI1, in CL4 fibroblasts. Many of the proteins that were upregulated in CL4 fibroblasts, as seen by unbiased proteomics, were also transcriptionally upregulated in the stroma of human breast cancers, especially in the patients that were prone to metastasis. Importantly, when CL4 fibroblasts were co-injected with human breast cancer cells (MDA-MB-231) in a xenograft model, tumor growth was dramatically enhanced. CL4 fibroblasts induced a >4-fold increase in tumor mass, and a near 8-fold increase in tumor volume, without any measurable increases in tumor angiogenesis. In parallel, CL3 and CL4 fibroblasts both failed to form tumors when they were injected alone, without epithelial cancer cells. Mechanistically, under co-culture conditions, CL4 glycolytic fibroblasts increased mitochondrial activity in adjacent breast cancer cells (relative to CL3 cells), consistent with the “Reverse Warburg Effect”. Notably, Western blot analysis of CL4 fibroblasts revealed a significant reduction in caveolin-1 (Cav-1) protein levels. In human breast cancer patients, a loss of stromal Cav-1 is associated with an increased risk of early tumor recurrence, metastasis, tamoxifen-resistance, and poor clinical outcome. Thus, loss of stromal Cav-1 may be an effective marker for predicting the “Reverse Warburg Effect” in the stroma of human breast cancer patients. As such, CL4 fibroblasts are a new attractive model for mimicking the “glycolytic phenotype” of cancer-associated fibroblasts. Nutrients derived from glycolytic cancer associated fibroblasts could provide an escape mechanism to confer drug-resistance during anti-angiogenic therapy, by effectively reducing the dependence of cancer cells on a vascular blood supply.