Tissue of origin determines cancer-associated CpG island promoter hypermethylation patterns

Background

Aberrant CpG island promoter DNA hypermethylation is frequently observed in cancer and is believed to contribute to tumor progression by silencing the expression of tumor suppressor genes. Previously, we observed that promoter hypermethylation in breast cancer reflects cell lineage rather than tumor progression and occurs at genes that are already repressed in a lineage-specific manner. To investigate the generality of our observation we analyzed the methylation profiles of 1,154 cancers from 7 different tissue types.

Results

We find that 1,009 genes are prone to hypermethylation in these 7 types of cancer. Nearly half of these genes varied in their susceptibility to hypermethylation between different cancer types. We show that the expression status of hypermethylation prone genes in the originator tissue determines their propensity to become hypermethylated in cancer; specifically, genes that are normally repressed in a tissue are prone to hypermethylation in cancers derived from that tissue. We also show that the promoter regions of hypermethylation-prone genes are depleted of repetitive elements and that DNA sequence around the same promoters is evolutionarily conserved. We propose that these two characteristics reflect tissue-specific gene promoter architecture regulating the expression of these hypermethylation prone genes in normal tissues.

Conclusions

As aberrantly hypermethylated genes are already repressed in pre-cancerous tissue, we suggest that their hypermethylation does not directly contribute to cancer development via silencing. Instead aberrant hypermethylation reflects developmental history and the perturbation of epigenetic mechanisms maintaining these repressed promoters in a hypomethylated state in normal cells.     

High Mutability of the Tumor Suppressor Genes RASSF1 and RBSP3 (CTDSPL) in Cancer

Background

Many different genetic alterations are observed in cancer cells. Individual cancer genes display point mutations such as base changes, insertions and deletions that initiate and promote cancer growth and spread. Somatic hypermutation is a powerful mechanism for generation of different mutations. It was shown previously that somatic hypermutability of proto-oncogenes can induce development of lymphomas.

Methodology/Principal Findings

We found an exceptionally high incidence of single-base mutations in the tumor suppressor genes RASSF1 and RBSP3 (CTDSPL) both located in 3p21.3 regions, LUCA and AP20 respectively. These regions contain clusters of tumor suppressor genes involved in multiple cancer types such as lung, kidney, breast, cervical, head and neck, nasopharyngeal, prostate and other carcinomas. Altogether in 144 sequenced RASSF1A clones (exons 1–2), 129 mutations were detected (mutation frequency, MF = 0.23 per 100 bp) and in 98 clones of exons 3–5 we found 146 mutations (MF = 0.29). In 85 sequenced RBSP3 clones, 89 mutations were found (MF = 0.10). The mutations were not cytidine-specific, as would be expected from alterations generated by AID/APOBEC family enzymes, and appeared de novo during cell proliferation. They diminished the ability of corresponding transgenes to suppress cell and tumor growth implying a loss of function. These high levels of somatic mutations were found both in cancer biopsies and cancer cell lines.

Conclusions/Significance

This is the first report of high frequencies of somatic mutations in RASSF1 and RBSP3 in different cancers suggesting it may underlay the mutator phenotype of cancer. Somatic hypermutations in tumor suppressor genes involved in major human malignancies offer a novel insight in cancer development, progression and spread.     

Tumor Endothelial Inflammation Predicts Clinical Outcome in Diverse Human Cancers

Background

Vascular endothelial cells contribute to the pathogenesis of numerous human diseases by actively regulating the stromal inflammatory response; however, little is known regarding the role of endothelial inflammation in the growth of human tumors and its influence on the prognosis of human cancers.

Methods

Using an experimental model of tumor necrosis factor-alpha (TNF-α)-mediated inflammation, we characterized inflammatory gene expression in immunopurified tumor-associated endothelial cells. These genes formed the basis of a multivariate molecular predictor of overall survival that was trained and validated in four types of human cancer.

Results

We report that expression of experimentally derived tumor endothelial genes distinguished pathologic tissue specimens from normal controls in several human diseases associated with chronic inflammation. We trained these genes in human cancer datasets and defined a six-gene inflammatory signature that predicted significantly reduced overall survival in breast cancer, colon cancer, lung cancer, and glioma. This endothelial-derived signature predicted outcome independently of, but cooperatively with, standard clinical and pathological prognostic factors. Consistent with these findings, conditioned culture media from human endothelial cells stimulated by pro-inflammatory cytokines accelerated the growth of human colon and breast tumors in immunodeficient mice as compared with conditioned media from untreated endothelial cells.

Conclusions

This study provides the first prognostic cancer gene signature derived from an experimental model of tumor-associated endothelial inflammation. These findings support the notion that activation of inflammatory pathways in non-malignant tumor-infiltrating endothelial cells contributes to tumor growth and progression in multiple human cancers. Importantly, these results identify endothelial-derived factors that could serve as potential targets for therapy in diverse human cancers.     

Functional Hypermethylation of ALDH1L1 , PLCL2 , and PPP2R3A in Colon Cancer

DNA hypermethylation and mutations are key mechanisms for the downregulation of tumor suppressor genes. NotI-microarrays allowed us to detect hypermethylation and/or deletions in 180 NotI sites associated with 188 genes of human chromosome 3, in 24 paired (tumor/normal) colon samples. The most frequent aberrations (in more than 20% of tumor samples) were detected in the promoter regions of 20 genes. Expression and promoter methylation of these genes were analyzed using the data for paired colon samples from The Cancer Genome Atlas project. Three genes—ALDH1L1, PLCL2, and PPP2R3A—revealed a more than two-fold average decrease in expression and a negative correlation between mRNA level and promoter hypermethylation. The expression of these three genes was then evaluated in 30 paired colon samples by quantitative PCR. Frequent (in more than 60% of cases) and significant (5–9-fold on average) mRNA level decrease was found for each of the genes in the tumor samples. The results indicate a suppressor role of the ALDH1L1, PLCL2, and PPP2R3A genes in colon cancer, as well as functional significance of hypermethylation in the downregulation of these genes.     

Examining the Relationship between Pre-Malignant Breast Lesions,Carcinogenesis and Tumor Evolution in the Mammary Epithelium Using an Agent-Based Model

Introduction

Breast cancer, the product of numerous rare mutational events that occur over an extended time period, presents numerous challenges to investigators interested in studying the transformation from normal breast epithelium to malignancy using traditional laboratory methods, particularly with respect to characterizing transitional and pre-malignant states. Dynamic computational modeling can provide insight into these pathophysiological dynamics, and as such we use a previously validated agent-based computational model of the mammary epithelium (the DEABM) to investigate the probabilistic mechanisms by which normal populations of ductal cells could transform into states replicating features of both pre-malignant breast lesions and a diverse set of breast cancer subtypes.

Methods

The DEABM consists of simulated cellular populations governed by algorithms based on accepted and previously published cellular mechanisms. Cells respond to hormones, undergo mitosis, apoptosis and cellular differentiation. Heritable mutations to 12 genes prominently implicated in breast cancer are acquired via a probabilistic mechanism. 3000 simulations of the 40-year period of menstrual cycling were run in wild-type (WT) and BRCA1-mutated groups. Simulations were analyzed by development of hyperplastic states, incidence of malignancy, hormone receptor and HER-2 status, frequency of mutation to particular genes, and whether mutations were early events in carcinogenesis.

Results

Cancer incidence in WT (2.6%) and BRCA1-mutated (45.9%) populations closely matched published epidemiologic rates. Hormone receptor expression profiles in both WT and BRCA groups also closely matched epidemiologic data. Hyperplastic populations carried more mutations than normal populations and mutations were similar to early mutations found in ER+ tumors (telomerase, E-cadherin, TGFB, RUNX3, p < .01). ER- tumors carried significantly more mutations and carried more early mutations in BRCA1, c-MYC and genes associated with epithelial-mesenchymal transition.

Conclusions

The DEABM generates diverse tumors that express tumor markers consistent with epidemiologic data. The DEABM also generates non-invasive, hyperplastic populations, analogous to atypia or ductal carcinoma in situ (DCIS), via mutations to genes known to be present in hyperplastic lesions and as early mutations in breast cancers. The results demonstrate that agent-based models are well-suited to studying tumor evolution through stages of carcinogenesis and have the potential to be used to develop prevention and treatment strategies.     

Expression Profiling of Mitochondrial Voltage-Dependent Anion Channel-1 Associated Genes Predicts Recurrence-Free Survival in Human Carcinomas

Background

Mitochondrial voltage-dependent anion channels (VDACs) play a key role in mitochondria-mediated apoptosis. Both in vivo and in vitro evidences indicate that VDACs are actively involved in tumor progression. Specifically, VDAC-1, one member of the VDAC family, was thought to be a potential anti-cancer therapeutic target. Our previous study demonstrated that the human gene VDAC1 (encoding the VDAC-1 isoform) was significantly up-regulated in lung tumor tissue compared with normal tissue. Also, we found a significant positive correlation between the gene expression of VDAC1 and histological grade in breast cancer. However, the prognostic power of VDAC1 and its associated genes in human cancers is largely unknown.

Methods

We systematically analyzed the expression pattern of VDAC1 and its interacting genes in breast, colon, liver, lung, pancreatic, and thyroid cancers. The genes differentially expressed between normal and tumor tissues in human carcinomas were identified.

Results

The expression level of VDAC1 was uniformly up-regulated in tumor tissue compared with normal tissue in breast, colon, liver, lung, pancreatic, and thyroid cancers. Forty-four VDAC1 interacting genes were identified as being commonly differentially expressed between normal and tumor tissues in human carcinomas. We designated VDAC1 and the 44 dysregulated interacting genes as the VDAC1 associated gene signature (VAG). We demonstrate that the VAG signature is a robust prognostic biomarker to predict recurrence-free survival in breast, colon, and lung cancers, and is independent of standard clinical and pathological prognostic factors.

Conclusions

VAG represents a promising prognostic biomarker in human cancers, which may enhance prediction accuracy in identifying patients at higher risk for recurrence. Future therapies aimed specifically at VDAC1 associated genes may lead to novel agents in the treatment of cancer.     

Activin Signaling in Microsatellite Stable Colon Cancers Is Disrupted by a Combination of Genetic and Epigenetic Mechanisms

Background

Activin receptor 2 (ACVR2) is commonly mutated in microsatellite unstable (MSI) colon cancers, leading to protein loss, signaling disruption, and larger tumors. Here, we examined activin signaling disruption in microsatellite stable (MSS) colon cancers.

Methods

Fifty-one population-based MSS colon cancers were assessed for ACVR1, ACVR2 and pSMAD2 protein. Consensus mutation-prone portions of ACVR2 were sequenced in primary cancers and all exons in colon cancer cell lines. Loss of heterozygosity (LOH) was evaluated for ACVR2 and ACVR1, and ACVR2 promoter methylation by methylation-specific PCR and bisulfite sequencing and chromosomal instability (CIN) phenotype via fluorescent LOH analysis of 3 duplicate markers. ACVR2 promoter methylation and ACVR2 expression were assessed in colon cancer cell lines via qPCR and IP-Western blots. Re-expression of ACVR2 after demethylation with 5-aza-2′-deoxycytidine (5-Aza) was determined. An additional 26 MSS colon cancers were assessed for ACVR2 loss and its mechanism, and ACVR2 loss in all tested cancers correlated with clinicopathological criteria.

Results

Of 51 MSS colon tumors, 7(14%) lost ACVR2, 2 (4%) ACVR1, and 5(10%) pSMAD2 expression. No somatic ACVR2 mutations were detected. Loss of ACVR2 expression was associated with LOH at ACVR2 (p<0.001) and ACVR2 promoter hypermethylation (p<0.05). ACVR2 LOH, but not promoter hypermethylation, correlated with CIN status. In colon cancer cell lines with fully methylated ACVR2 promoter, loss of ACVR2 mRNA and protein expression was restored with 5-Aza treatment. Loss of ACVR2 was associated with an increase in primary colon cancer volume (p<0.05).

Conclusions

Only a small percentage of MSS colon cancers lose expression of activin signaling members. ACVR2 loss occurs through LOH and ACVR2 promoter hypermethylation, revealing distinct mechanisms for ACVR2 inactivation in both MSI and MSS subtypes of colon cancer.     

Frequent Epigenetic Silencing of the Folate-Metabolising Gene Cystathionine-Beta-Synthase in Gastrointestinal Cancer

Background

Both gastric and colorectal cancers (CRC) are the most frequently occurring malignancies worldwide with the overall survival of these patients remains unsatisfied. Identification of tumor suppressor genes (TSG) silenced by promoter CpG methylation uncovers mechanisms of tumorigenesis and identifies new epigenetic biomarkers for early cancer detection and prognosis assessment. Cystathionine-beta-synthase (CBS) functions in the folate metabolism pathway, which is intricately linked to methylation of genomic DNA. Dysregulation of DNA methylation contributes substantially to cancer development.

Methodology/Principal Findings

To identify potential TSGs silenced by aberrant promoter methylation in CRC, we analyzed tumor and adjacent tissues from CRC cases using the Illumina Human Methylation45 BeadChip. We identified hypermethylation of the CBS gene in CRC samples, compared to adjacent tissues. Methylation and decreased mRNA expression of CBS were detected in most CRC cell lines by methylation-specific PCR and semiquantitative RT-PCR, as well as in gastric cancer. Treatment with 5-aza-2''-deoxycytidine and/or trichostatin A reversed methylation and restored CBS mRNA expression indicating a direct effect. Aberrant methylation was further detected in 31% of primary CRCs (29 of 96) and 55% of gastric tumors (11 of 20). In contrast, methylation was seldom found in normal tissues adjacent to the tumor. CBS methylation was associated with KRAS mutations in primary CRCs (P = 0.04, by χ²-test). However, no association was found between CBS methylation or KRAS mutations with cancer relapse/metastasis in Stage II CRC patients.

Conclusion

A novel finding from this study is that the folate metabolism enzyme CBS mRNA levels are frequently downregulated through CpG methylation of the CBS gene in gastric cancer and CRC, suggesting that CBS functions as a tumor suppressor gene. These findings warrant further study of CBS as an epigenetic biomarker for molecular diagnosis of gastrointestinal cancers.     

Transcriptional Regulation and Biological Functions of Selenium-Binding Protein 1 in Colorectal Cancer In Vitro and in Nude Mouse Xenografts

Background

It has been shown that selenium-binding protein 1 (SBP1) is significantly downregulated in different human cancers. Its regulation and function have not yet been established.

Methodology and Principal Findings

We show that the SBP1 promoter is hypermethylated in colon cancer tissues and human colon cancer cells. Treatment with 5′-Aza-2′-deoxycytidine leads to demethylation of the SBP1 promoter and to an increase of SBP1 promoter activity, rescues SBP1 mRNA and protein expression in human colon cancer cells. Additionally, overexpression of SBP1 sensitizes colon cancer cells to H₂O₂-induced apoptosis, inhibits cancer cell migration in vitro and inhibits tumor growth in nude mice.

Conclusion and Significance

These data demonstrate that SBP1 has tumor suppressor functions that are inhibited in colorectal cancer through epigenetic silencing.     

Prevalence of Germline Mutations in Genes Engaged in DNA Damage Repair by Homologous Recombination in Patients with Triple-Negative and Hereditary Non-Triple-Negative Breast Cancers

Purpose

This study sought to assess the prevalence of common germline mutations in several genes engaged in the repair of DNA double-strand break by homologous recombination in patients with triple-negative breast cancers and hereditary non-triple-negative breast cancers. Tumors deficient in this type of DNA damage repair are known to be especially sensitive to DNA cross-linking agents (e.g., platinum drugs) and to poly(ADP-ribose) polymerase (PARP) inhibitors.

Methods

Genetic testing was performed for 36 common germline mutations in genes engaged in the repair of DNA by homologous recombination, i.e., BRCA1, BRCA2, CHEK2, NBN, ATM, PALB2, BARD1, and RAD51D, in 202 consecutive patients with triple-negative breast cancers and hereditary non-triple-negative breast cancers.

Results

Thirty five (22.2%) of 158 patients in the triple-negative group carried mutations in genes involved in DNA repair by homologous recombination, while 10 (22.7%) of the 44 patients in the hereditary non-triple-negative group carried such mutations. Mutations in BRCA1 were most frequent in patients with triple-negative breast cancer (18.4%), and mutations in CHEK2 were most frequent in patients with hereditary non-triple-negative breast cancers (15.9%). In addition, in the triple-negative group, mutations in CHEK2, NBN, and ATM (3.8% combined) were found, while mutations in BRCA1, NBN, and PALB2 (6.8% combined) were identified in the hereditary non-triple-negative group.

Conclusions

Identifying mutations in genes engaged in DNA damage repair by homologous recombination other than BRCA1/2 can substantially increase the proportion of patients with triple-negative breast cancer and hereditary non-triple-negative breast cancer who may be eligible for therapy using PARP inhibitors and platinum drugs.     

Somatic mutation profiles of MSI and MSS colorectal cancer identified by whole exome next generation sequencing and bioinformatics analysis

Background

Colorectal cancer (CRC) is with approximately 1 million cases the third most common cancer worldwide. Extensive research is ongoing to decipher the underlying genetic patterns with the hope to improve early cancer diagnosis and treatment. In this direction, the recent progress in next generation sequencing technologies has revolutionized the field of cancer genomics. However, one caveat of these studies remains the large amount of genetic variations identified and their interpretation.

Methodology/Principal Findings

Here we present the first work on whole exome NGS of primary colon cancers. We performed 454 whole exome pyrosequencing of tumor as well as adjacent not affected normal colonic tissue from microsatellite stable (MSS) and microsatellite instable (MSI) colon cancer patients and identified more than 50,000 small nucleotide variations for each tissue. According to predictions based on MSS and MSI pathomechanisms we identified eight times more somatic non-synonymous variations in MSI cancers than in MSS and we were able to reproduce the result in four additional CRCs. Our bioinformatics filtering approach narrowed down the rate of most significant mutations to 359 for MSI and 45 for MSS CRCs with predicted altered protein functions. In both CRCs, MSI and MSS, we found somatic mutations in the intracellular kinase domain of bone morphogenetic protein receptor 1A, BMPR1A, a gene where so far germline mutations are associated with juvenile polyposis syndrome, and show that the mutations functionally impair the protein function.

Conclusions/Significance

We conclude that with deep sequencing of tumor exomes one may be able to predict the microsatellite status of CRC and in addition identify potentially clinically relevant mutations.     

Integrated analyses of DNA methylation and hydroxymethylation reveal tumor suppressive roles of ECM1, ATF5, and EOMES in human hepatocellular carcinoma

Background

Differences in 5-hydroxymethylcytosine, 5hmC, distributions may complicate previous observations of abnormal cytosine methylation statuses that are used for the identification of new tumor suppressor gene candidates that are relevant to human hepatocarcinogenesis. The simultaneous detection of 5-methylcytosine and 5-hydroxymethylcytosine is likely to stimulate the discovery of aberrantly methylated genes with increased accuracy in human hepatocellular carcinoma.

Results

Here, we performed ultra-performance liquid chromatography/tandem mass spectrometry and single-base high-throughput sequencing, Hydroxymethylation and Methylation Sensitive Tag sequencing, HMST-seq, to synchronously measure these two modifications in human hepatocellular carcinoma samples. After identification of differentially methylated and hydroxymethylated genes in human hepatocellular carcinoma, we integrate DNA copy-number alterations, as determined using array-based comparative genomic hybridization data, with gene expression to identify genes that are potentially silenced by promoter hypermethylation.

Conclusions

We report a high enrichment of genes with epigenetic aberrations in cancer signaling pathways. Six genes were selected as tumor suppressor gene candidates, among which, ECM1, ATF5 and EOMES are confirmed via siRNA experiments to have potential anti-cancer functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0533-9) contains supplementary material, which is available to authorized users.     

Frequent alteration of MLL3 frameshift mutations in microsatellite deficient colorectal cancer

Background

MLL3 is a histone 3- lysine 4 methyltransferase with tumor-suppressor properties that belongs to a family of chromatin regulator genes potentially altered in neoplasia. Mutations in MLL3 were found in a whole genome analysis of colorectal cancer but have not been confirmed by a separate study.

Methods and Results

We analyzed mutations of coding region and promoter methylation in MLL3 using 126 cases of colorectal cancer. We found two isoforms of MLL3 and DNA sequencing revealed frameshift and other mutations affecting both isoforms of MLL3 in colorectal cancer cells and 19 of 134 (14%) primary colorectal samples analyzed. Moreover, frameshift mutations were more common in cases with microsatellite instability (31%) both in CRC cell lines and primary tumors. The largest isoform of MLL3 is transcribed from a CpG island-associated promoter that has highly homology with a pseudo-gene on chromosome 22 (psiTPTE22). Using an assay which measured both loci simultaneously we found prominent age related methylation in normal colon (from 21% in individuals less than 25 years old to 56% in individuals older than 70, R = 0.88, p<0.001) and frequent hypermethylation (83%) in both CRC cell lines and primary tumors. We next studied the two loci separately and found that age and cancer related methylation was solely a property of the pseudogene CpG island and that the MLL3 loci was unmethylated.

Conclusions

We found that frameshift mutations of MLL3 in both CRC cells and primary tumor that were more common in cases with microsatellite instability. Moreover, we have shown CpG island-associated promoter of MLL3 gene has no DNA methylation in CRC cells but also primary tumor and normal colon, and this region has a highly homologous of pseudo gene (psiTPTE22) that was age relate DNA methylation.     

Differential sensitivity of ERBB2 kinase domain mutations towards lapatinib

Background

Overexpression of the ERBB2 kinase is observed in about one-third of breast cancer patients and the dual ERBB1/ERBB2 kinase inhibitor lapatinib was recently approved for the treatment of advanced ERBB2-positive breast cancer. Mutations in the ERBB2 receptor have recently been reported in breast cancer at diagnosis and also in gastric, colorectal and lung cancer. These mutations may have an impact on the clinical responses achieved with lapatinib in breast cancer and may also have a potential impact on the use of lapatinib in other solid cancers. However, the sensitivity of lapatinib towards clinically observed ERBB2 mutations is not known.

Methodology/Principal Findings

We cloned a panel of 8 clinically observed ERBB2 mutations, established stable cell lines and characterized their sensitivity towards lapatinib and alternative ERBB2 inhibitors. Both lapatinib-sensitive and lapatinib-resistant ERBB2 mutations were observed. Interestingly, we were able to generate lapatinib resistance mutations in wt-ERBB2 cells incubated with lapatinib for prolonged periods of time. This indicates that these resistance mutations may also cause secondary resistance in lapatinib-treated patients. Lapatinib-resistant ERBB2 mutations were found to be highly resistant towards AEE788 treatment but remained sensitive towards the dual irreversible inhibitors CL-387785 and WZ-4002.

Conclusions/Significance

Patients harbouring certain ERBB2 kinase domain mutations at diagnosis may not benefit from lapatinib treatment. Moreover, secondary lapatinib resistance may develop due to kinase domain mutations. Irreversible ERBB2 inhibitors may offer alternative treatment options for breast cancer and other solid tumor patients harbouring lapatinib resistance mutations. In addition, these inhibitors may be of interest in the scenario of secondary lapatinib resistance.