Genome-wide DNA methylation profiling of recurrent and non-recurrent chordomas

Chordomas are an aggressive rare type of malignant bone tumors arising from the remnant of the notochord. Chordomas occur mainly in vertebral bones and account for 1–4% of malignant bone tumors. Management and treatment of chordomas are difficult as they are resistant to conventional chemotherapy; therefore, they are mainly treated with surgery and radiation therapy. In this study, we performed DNA methylation profiling of 26 chordomas and normal nucleus pulposus samples plus UCH-1 chordoma cell line using the Illumina Infinium HumanMethylation450 BeadChips. Combined bisulfite restriction analysis and bisulfite sequencing was used to confirm the methylation data. Gene expression was analyzed using RT-PCR before and after 5-aza-2’-deoxycytidine (5-azaDC) treatment of chordoma cell lines. Analysis of the HumanMethylation450 BeadChip data led to the identification of 8,819 loci (2.9%) that were significantly differentially methylated (>0.2 average β-value difference) between chordomas and nucleus pulposus samples (adjusted P < 0.05). Among these, 5,868 probes (66.5%) were hypomethylated, compared to 2,951 (33.5%) loci that were hypermethylated in chordomas compared to controls. From the 2,951 differentially hypermethylated probes, 33.3% were localized in the promoter region (982 probes) and, among these, 104 probes showed cancer-specific hypermethylation. Ingenuity Pathway Analysis indicates that the cancer-specific differentially methylated loci are involved in various networks including cancer disease, nervous system development and function, cell death and survival, cellular growth, cellular development, and proliferation. Furthermore, we identified a subset of probes that were differentially methylated between recurrent and non-recurrent chordomas. BeadChip methylation data was confirmed for these genes and gene expression was shown to be upregulated in methylated chordoma cell lines after treatment with 5-azaDC. Understanding epigenetic changes in chordomas may provide insights into chordoma tumorigenesis and development of epigenetic biomarkers.     

Recurrent chromosomal copy number alterations in sporadic chordomas

The molecular events in chordoma pathogenesis have not been fully delineated, particularly with respect to copy number changes. Understanding copy number alterations in chordoma may reveal critical disease mechanisms that could be exploited for tumor classification and therapy. We report the copy number analysis of 21 sporadic chordomas using array comparative genomic hybridization (CGH). Recurrent copy changes were further evaluated with immunohistochemistry, methylation specific PCR, and quantitative real-time PCR. Similar to previous findings, large copy number losses, involving chromosomes 1p, 3, 4, 9, 10, 13, 14, and 18, were more common than copy number gains. Loss of CDKN2A with or without loss of CDKN2B on 9p21.3 was observed in 16/20 (80%) unique cases of which six (30%) showed homozygous deletions ranging from 76 kilobases to 4.7 megabases. One copy loss of the 10q23.31 region which encodes PTEN was found in 16/20 (80%) cases. Loss of CDKN2A and PTEN expression in the majority of cases was not attributed to promoter methylation. Our sporadic chordoma cases did not show hotspot point mutations in some common cancer gene targets. Moreover, most of these sporadic tumors are not associated with T (brachyury) duplication or amplification. Deficiency of CDKN2A and PTEN expression, although shared across many other different types of tumors, likely represents a key aspect of chordoma pathogenesis. Sporadic chordomas may rely on mechanisms other than copy number gain if they indeed exploit T/brachyury for proliferation.     

Genotyping Cancer-Associated Genes in Chordoma Identifies Mutations in Oncogenes and Areas of Chromosomal Loss Involving CDKN2A,PTEN, and SMARCB1

The molecular mechanisms underlying chordoma pathogenesis are unknown. We therefore sought to identify novel mutations to better understand chordoma biology and to potentially identify therapeutic targets. Given the relatively high costs of whole genome sequencing, we performed a focused genetic analysis using matrix-assisted laser desorption/ionization-time of flight mass spectrometer (Sequenom iPLEX genotyping). We tested 865 hotspot mutations in 111 oncogenes and selected tumor suppressor genes (OncoMap v. 3.0) of 45 human chordoma tumor samples. Of the analyzed samples, seven were identified with at least one mutation. Six of these were from fresh frozen samples, and one was from a paraffin embedded sample. These observations were validated using an independent platform using homogeneous mass extend MALDI-TOF (Sequenom hME Genotyping). These genetic alterations include: ALK (A877S), CTNNB1 (T41A), NRAS (Q61R), PIK3CA (E545K), PTEN (R130), CDKN2A (R58*), and SMARCB1 (R40*). This study reports on the largest comprehensive mutational analysis of chordomas performed to date. To focus on mutations that have the greatest chance of clinical relevance, we tested only oncogenes and tumor suppressor genes that have been previously implicated in the tumorigenesis of more common malignancies. We identified rare genetic changes that may have functional significance to the underlying biology and potential therapeutics for chordomas. Mutations in CDKN2A and PTEN occurred in areas of chromosomal copy loss. When this data is paired with the studies showing 18 of 21 chordoma samples displaying copy loss at the locus for CDKN2A, 17 of 21 chordoma samples displaying copy loss at PTEN, and 3 of 4 chordoma samples displaying deletion at the SMARCB1 locus, we can infer that a loss of heterozygosity at these three loci may play a significant role in chordoma pathogenesis.     

Dnmt3a Protects Active Chromosome Domains against Cancer-Associated Hypomethylation

Changes in genomic DNA methylation patterns are generally assumed to play an important role in the etiology of human cancers. The Dnmt3a enzyme is required for the establishment of normal methylation patterns, and mutations in Dnmt3a have been described in leukemias. Deletion of Dnmt3a in a K-ras–dependent mouse lung cancer model has been shown to promote tumor progression, which suggested that the enzyme might suppress tumor development by stabilizing DNA methylation patterns. We have used whole-genome bisulfite sequencing to comprehensively characterize the methylomes from Dnmt3a wildtype and Dnmt3a-deficient mouse lung tumors. Our results show that profound global methylation changes can occur in K-ras–induced lung cancer. Dnmt3a wild-type tumors were characterized by large hypomethylated domains that correspond to nuclear lamina-associated domains. In contrast, Dnmt3a-deficient tumors showed a uniformly hypomethylated genome. Further data analysis revealed that Dnmt3a is required for efficient maintenance methylation of active chromosome domains and that Dnmt3a-deficient tumors show moderate levels of gene deregulation in these domains. In summary, our results uncover conserved features of cancer methylomes and define the role of Dnmt3a in maintaining DNA methylation patterns in cancer.     

DNA Methylation Profile Distinguishes Clear Cell Sarcoma of the Kidney from Other Pediatric Renal Tumors

A number of specific, distinct neoplastic entities occur in the pediatric kidney, including Wilms’ tumor, clear cell sarcoma of the kidney (CCSK), congenital mesoblastic nephroma (CMN), rhabdoid tumor of the kidney (RTK), and the Ewing’s sarcoma family of tumors (ESFT). By employing DNA methylation profiling using Illumina Infinium HumanMethylation27, we analyzed the epigenetic characteristics of the sarcomas including CCSK, RTK, and ESFT in comparison with those of the non-neoplastic kidney (NK), and these tumors exhibited distinct DNA methylation profiles in a tumor-type-specific manner. CCSK is the most frequently hypermethylated, but least frequently hypomethylated, at CpG sites among these sarcomas, and exhibited 490 hypermethylated and 46 hypomethylated CpG sites in compared with NK. We further validated the results by MassARRAY, and revealed that a combination of four genes was sufficient for the DNA methylation profile-based differentiation of these tumors by clustering analysis. Furthermore, THBS1 CpG sites were found to be specifically hypermethylated in CCSK and, thus, the DNA methylation status of these THBS1 sites alone was sufficient for the distinction of CCSK from other pediatric renal tumors, including Wilms’ tumor and CMN. Moreover, combined bisulfite restriction analysis could be applied for the detection of hypermethylation of a THBS1 CpG site. Besides the biological significance in the pathogenesis, the DNA methylation profile should be useful for the differential diagnosis of pediatric renal tumors.     

Identification of tumor suppressors and oncogenes from genomic and epigenetic features in ovarian cancer

The identification of genetic and epigenetic alterations from primary tumor cells has become a common method to identify genes critical to the development and progression of cancer. We seek to identify those genetic and epigenetic aberrations that have the most impact on gene function within the tumor. First, we perform a bioinformatic analysis of copy number variation (CNV) and DNA methylation covering the genetic landscape of ovarian cancer tumor cells. We separately examined CNV and DNA methylation for 42 primary serous ovarian cancer samples using MOMA-ROMA assays and 379 tumor samples analyzed by The Cancer Genome Atlas. We have identified 346 genes with significant deletions or amplifications among the tumor samples. Utilizing associated gene expression data we predict 156 genes with altered copy number and correlated changes in expression. Among these genes CCNE1, POP4, UQCRB, PHF20L1 and C19orf2 were identified within both data sets. We were specifically interested in copy number variation as our base genomic property in the prediction of tumor suppressors and oncogenes in the altered ovarian tumor. We therefore identify changes in DNA methylation and expression for all amplified and deleted genes. We statistically define tumor suppressor and oncogenic features for these modalities and perform a correlation analysis with expression. We predicted 611 potential oncogenes and tumor suppressors candidates by integrating these data types. Genes with a strong correlation for methylation dependent expression changes exhibited at varying copy number aberrations include CDCA8, ATAD2, CDKN2A, RAB25, AURKA, BOP1 and EIF2C3. We provide copy number variation and DNA methylation analysis for over 11,500 individual genes covering the genetic landscape of ovarian cancer tumors. We show the extent of genomic and epigenetic alterations for known tumor suppressors and oncogenes and also use these defined features to identify potential ovarian cancer gene candidates.     

Non-methylated Genomic Sites Coincidence Cloning (NGSCC): an approach to large scale analysis of hypomethylated CpG patterns at predetermined genomic loci

We have developed a new approach to the analysis of hypomethylated CpG patterns within predetermined, megabase long, genome regions. The approach, which we term Non-methylated Genomic Sites Coincidence Cloning (NGSCC), includes three main steps. First, total genomic DNA is digested with a methylation sensitive restriction endonuclease, such as Hpa II or Hha I. Then the fragments corresponding to the genomic area of interest are selected. To this end the fragmented genome DNA is hybridized with a mixture of clones (BACs, cosmids etc.) representing a given region and digested with the same restriction enzyme(s). A special version of the coincidence cloning procedure was developed to make this hybridization selection highly efficient and specific. Finally, fragments of the locus under study are mapped and sequenced. The technique proved to be efficient and specific. As a test, it was applied to the analysis of hypomethylated CpG patterns along the 1-Mb D19S208-COX7A1 (Chr 19q13.12) locus, on human chromosome 19, in normal testis and in seminoma tissues. Some differences in the distribution of hypomethylated CpGs between the two tissues were demonstrated. The methylation profiles in both tissues revealed a clear trend to clustering of non-methylated sites. We also analyzed the expression of genes located within hypomethylated clusters in both tissues. It was shown that, whereas the expression of some of the genes investigated was correlated with hypomethylation of the region, other genes were expressed regardless of their methylation status. NGSCC thus promises to be a useful approach for the analysis of the role of dynamic epigenetic factors in genome function.Communicated by G. P. Georgiev     

A Novel Chordoma Xenograft Allows In Vivo Drug Testing and Reveals the Importance of NF-κB Signaling in Chordoma Biology

Chordoma is a rare primary bone malignancy that arises in the skull base, spine and sacrum and originates from remnants of the notochord. These tumors are typically resistant to conventional chemotherapy, and to date there are no FDA-approved agents to treat chordoma. The lack of in vivo models of chordoma has impeded the development of new therapies for this tumor. Primary tumor from a sacral chordoma was xenografted into NOD/SCID/IL-2R γ-null mice. The xenograft is serially transplantable and was characterized by both gene expression analysis and whole genome SNP genotyping. The NIH Chemical Genomics Center performed high-throughput screening of 2,816 compounds using two established chordoma cell lines, U-CH1 and U-CH2B. The screen yielded several compounds that showed activity and two, sunitinib and bortezomib, were tested in the xenograft. Both agents slowed the growth of the xenograft tumor. Sensitivity to an inhibitor of IκB, as well as inhibition of an NF-κB gene expression signature demonstrated the importance of NF-κB signaling for chordoma growth. This serially transplantable chordoma xenograft is thus a practical model to study chordomas and perform in vivo preclinical drug testing.     

Long interspersed nuclear element-1 hypomethylation in cancer: biology and clinical applications

Epigenetic changes in long interspersed nuclear element-1s (LINE-1s or L1s) occur early during the process of carcinogenesis. A lower methylation level (hypomethylation) of LINE-1 is common in most cancers, and the methylation level is further decreased in more advanced cancers. Consequently, several previous studies have suggested the use of LINE-1 hypomethylation levels in cancer screening, risk assessment, tumor staging, and prognostic prediction. Epigenomic changes are complex, and global hypomethylation influences LINE-1s in a generalized fashion. However, the methylation levels of some loci are dependent on their locations. The consequences of LINE-1 hypomethylation are genomic instability and alteration of gene expression. There are several mechanisms that promote both of these consequences in cis. Therefore, the methylation levels of different sets of LINE-1s may represent certain phenotypes. Furthermore, the methylation levels of specific sets of LINE-1s may indicate carcinogenesis-dependent hypomethylation. LINE-1 methylation pattern analysis can classify LINE-1s into one of three classes based on the number of methylated CpG dinucleotides. These classes include hypermethylation, partial methylation, and hypomethylation. The number of partial and hypermethylated loci, but not hypomethylated LINE-1s, is different among normal cell types. Consequently, the number of hypomethylated loci is a more promising marker than methylation level in the detection of cancer DNA. Further genome-wide studies to measure the methylation level of each LINE-1 locus may improve PCR-based methylation analysis to allow for a more specific and sensitive detection of cancer DNA or for an analysis of certain cancer phenotypes. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s13148-011-0032-8) contains supplementary material, which is available to authorized users.     

Breast Tumors with Elevated Expression of 1q Candidate Genes Confer Poor Clinical Outcome and Sensitivity to Ras/PI3K Inhibition

Genomic aberrations are common in cancers and the long arm of chromosome 1 is known for its frequent amplifications in breast cancer. However, the key candidate genes of 1q, and their contribution in breast cancer pathogenesis remain unexplored. We have analyzed the gene expression profiles of 1635 breast tumor samples using meta-analysis based approach and identified clinically significant candidates from chromosome 1q. Seven candidate genes including exonuclease 1 (EXO1) are consistently over expressed in breast tumors, specifically in high grade and aggressive breast tumors with poor clinical outcome. We derived a EXO1 co-expression module from the mRNA profiles of breast tumors which comprises 1q candidate genes and their co-expressed genes. By integrative functional genomics investigation, we identified the involvement of EGFR, RAS, PI3K / AKT, MYC, E2F signaling in the regulation of these selected 1q genes in breast tumors and breast cancer cell lines. Expression of EXO1 module was found as indicative of elevated cell proliferation, genomic instability, activated RAS/AKT/MYC/E2F1 signaling pathways and loss of p53 activity in breast tumors. mRNA–drug connectivity analysis indicates inhibition of RAS/PI3K as a possible targeted therapeutic approach for the patients with activated EXO1 module in breast tumors. Thus, we identified seven 1q candidate genes strongly associated with the poor survival of breast cancer patients and identified the possibility of targeting them with EGFR/RAS/PI3K inhibitors.     

Genome-wide CpG island methylation and intergenic demethylation propensities vary among different tumor sites

The epigenetic landscape of cancer includes both focal hypermethylation and broader hypomethylation in a genome-wide manner. By means of a comprehensive genomic analysis on 6637 tissues of 21 tumor types, we here show that the degrees of overall methylation in CpG island (CGI) and demethylation in intergenic regions, defined as ‘backbone’, largely vary among different tumors. Depending on tumor type, both CGI methylation and backbone demethylation are often associated with clinical, epidemiological and biological features such as age, sex, smoking history, anatomic location, histological type and grade, stage, molecular subtype and biological pathways. We found connections between CGI methylation and hypermutability, microsatellite instability, IDH1 mutation, 19p gain and polycomb features, and backbone demethylation with chromosomal instability, NSD1 and TP53 mutations, 5q and 19p loss and long repressive domains. These broad epigenetic patterns add a new dimension to our understanding of tumor biology and its clinical implications.     

Estimating genomic instability mediated by Alu retroelements in breast cancer

Alu-PCR is a relatively simple technique that can be used to investigate genomic instability in cancer. This technique allows identification of the loss, gain or amplification of gene sequences based on the analysis of segments between two Alu elements coupled with quantitative and qualitative analyses of the profiles obtained from tumor samples, surgical margins and blood. In this work, we used Alu-PCR to identify gene alterations in ten patients with invasive ductal breast cancer. Several deletions and insertions were identified, indicating genomic instability in the tumor and adjacent normal tissue. Although not associated with specific genes, the alterations, which involved chromosomal bands 1p36.23, 1q41, 11q14.3, 13q14.2, occurred in areas of well-known genomic instability in breast and other types of cancer. These results indicate the potential usefulness of Alu-PCR in identifying altered gene sequences in breast cancer. However, caution is required in its application since the Alu primer can produce non-specific amplification.     

Human BRCA1-associated breast cancer: no increase in numerical chromosomal instability compared to sporadic tumors

BRCA1 is a major gatekeeper of genomic stability. Acting in multiple central processes like double-strand break repair, centrosome replication, and checkpoint control, BRCA1 participates in maintaining genomic integrity and protects the cell against genomic instability. Chromosomal instability (CIN) as part of genomic instability is an inherent characteristic of most solid tumors and is also involved in breast cancer development. In this study, we determined the extent of CIN in 32 breast cancer tumors of women with a BRCA1 germline mutation compared to 62 unselected breast cancers. We applied fluorescence in situ hybridization (FISH) with centromere-specific probes for the chromosomes 1, 7, 8, 10, 17, and X and locus-specific probes for 3q27 (BCL6), 5p15.2 (D5S23), 5q31 (EGR1), 10q23.3 (PTEN), and 14q32 (IGH@) on formalin-fixed paraffin-embedded tissue microarray sections. Our hypothesis of an increased level of CIN in BRCA1-associated breast cancer could not be confirmed by this approach. Surprisingly, we detected no significant difference in the extent of CIN in BRCA1-mutated versus sporadic tumors. The only exception was the CIN value for chromosome 1. Here, the extent of CIN was slightly higher in the group of sporadic tumors.     

Candidate Luminal B Breast Cancer Genes Identified by Genome,Gene Expression and DNA Methylation Profiling

Breast cancers (BCs) of the luminal B subtype are estrogen receptor-positive (ER+), highly proliferative, resistant to standard therapies and have a poor prognosis. To better understand this subtype we compared DNA copy number aberrations (CNAs), DNA promoter methylation, gene expression profiles, and somatic mutations in nine selected genes, in 32 luminal B tumors with those observed in 156 BCs of the other molecular subtypes. Frequent CNAs included 8p11-p12 and 11q13.1-q13.2 amplifications, 7q11.22-q34, 8q21.12-q24.23, 12p12.3-p13.1, 12q13.11-q24.11, 14q21.1-q23.1, 17q11.1-q25.1, 20q11.23-q13.33 gains and 6q14.1-q24.2, 9p21.3-p24,3, 9q21.2, 18p11.31-p11.32 losses. A total of 237 and 101 luminal B-specific candidate oncogenes and tumor suppressor genes (TSGs) presented a deregulated expression in relation with their CNAs, including 11 genes previously reported associated with endocrine resistance. Interestingly, 88% of the potential TSGs are located within chromosome arm 6q, and seven candidate oncogenes are potential therapeutic targets. A total of 100 candidate oncogenes were validated in a public series of 5,765 BCs and the overexpression of 67 of these was associated with poor survival in luminal tumors. Twenty-four genes presented a deregulated expression in relation with a high DNA methylation level. FOXO3, PIK3CA and TP53 were the most frequent mutated genes among the nine tested. In a meta-analysis of next-generation sequencing data in 875 BCs, KCNB2 mutations were associated with luminal B cases while candidate TSGs MDN1 (6q15) and UTRN (6q24), were mutated in this subtype. In conclusion, we have reported luminal B candidate genes that may play a role in the development and/or hormone resistance of this aggressive subtype.     

Integrated miRNA-mRNA Analysis Revealing the Potential Roles of miRNAs in Chordomas

Introduction

Emerging evidence suggests that microRNAs (miRNAs) are crucially involved in tumorigenesis and that paired expression profiles of miRNAs and mRNAs can be used to identify functional miRNA-target relationships with high precision. However, no studies have applied integrated analysis to miRNA and mRNA profiles in chordomas. The purpose of this study was to provide insights into the pathogenesis of chordomas by using this integrated analysis method.

Methods

Differentially expressed miRNAs and mRNAs of chordomas (n = 3) and notochord tissues (n = 3) were analyzed by using microarrays with hierarchical clustering analysis. Subsequently, the target genes of the differentially expressed miRNAs were predicted and overlapped with the differentially expressed mRNAs. Then, GO and pathway analyses were performed for the intersecting genes.

Results

The microarray analysis indicated that 33 miRNAs and 2,791 mRNAs were significantly dysregulated between the two groups. Among the 2,791 mRNAs, 911 overlapped with putative miRNA target genes. A pathway analysis showed that the MAPK pathway was consistently enriched in the chordoma tissue and that miR-149-3p, miR-663a, miR-1908, miR-2861 and miR-3185 likely play important roles in the regulation of MAPK pathways. Furthermore, the Notch signaling pathway and the loss of the calcification or ossification capacity of the notochord may also be involved in chordoma pathogenesis.

Conclusion

This study provides an integrated dataset of the miRNA and mRNA profiles in chordomas, and the results demonstrate that not only the MAPK pathway and its related miRNAs but also the Notch pathway may be involved in chordoma development. The occurrence of chordoma may be associated with dysfunctional calcification or ossification of the notochord.     

A survey of tandem repeat instabilities and associated gene expression changes in 35 colorectal cancers

Background

Colorectal cancer is a major contributor to cancer morbidity and mortality. Tandem repeat instability and its effect on cancer phenotypes remain so far poorly studied on a genome-wide scale.

Results

Here we analyze the genomes of 35 colorectal tumors and their matched normal (healthy) tissues for two types of tandem repeat instability, de-novo repeat gain or loss and repeat copy number variation. Specifically, we study for the first time genome-wide repeat instability in the promoters and exons of 18,439 genes, and examine the association of repeat instability with genome-scale gene expression levels. We find that tumors with a microsatellite instable (MSI) phenotype are enriched in genes with repeat instability, and that tumor genomes have significantly more genes with repeat instability compared to healthy tissues. Genes in tumor genomes with repeat instability in their promoters are significantly less expressed and show slightly higher levels of methylation. Genes in well-studied cancer-associated signaling pathways also contain significantly more unstable repeats in tumor genomes. Genes with such unstable repeats in the tumor-suppressor p53 pathway have lower expression levels, whereas genes with repeat instability in the MAPK and Wnt signaling pathways are expressed at higher levels, consistent with the oncogenic role they play in cancer.

Conclusions

Our results suggest that repeat instability in gene promoters and associated differential gene expression may play an important role in colorectal tumors, which is a first step towards the development of more effective molecular diagnostic approaches centered on repeat instability.

Electronic supplementary material

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