Overlay analysis of the oligonucleotide array gene expression profiles and copy number abnormalities as determined by array comparative genomic hybridization in medulloblastomas

OTX2 sustains a bivalent-like state of OTX2-bound promoters in medulloblastoma by maintaining their H3K27me3 levels

Recent studies showed frequent mutations in histone H3 lysine 27 (H3K27) demethylases in medulloblastomas of Group 3 and Group 4, suggesting a role for H3K27 methylation in these tumors. Indeed, trimethylated H3K27 (H3K27me3) levels were shown to be higher in Group 3 and 4 tumors compared to WNT and SHH medulloblastomas, also in tumors without detectable mutations in demethylases. Here, we report that polycomb genes, required for H3K27 methylation, are consistently upregulated in Group 3 and 4 tumors. These tumors show high expression of the homeobox transcription factor OTX2. Silencing of OTX2 in D425 medulloblastoma cells resulted in downregulation of polycomb genes such as EZH2, EED, SUZ12 and RBBP4 and upregulation of H3K27 demethylases KDM6A, KDM6B, JARID2 and KDM7A. This was accompanied by decreased H3K27me3 and increased H3K27me1 levels in promoter regions. Strikingly, the decrease of H3K27me3 was most prominent in promoters that bind OTX2. OTX2-bound promoters showed high levels of the H3K4me3 and H3K9ac activation marks and intermediate levels of the H3K27me3 inactivation mark, reminiscent of a bivalent modification. After silencing of OTX2, H3K27me3 levels strongly dropped, but H3K4me3 and H3K9ac levels remained high. OTX2-bound bivalent genes showed high expression levels in D425, but the expression of most of these genes did not change after OTX2 silencing and loss of the H3K27me3 mark. Maintaining promoters in a bivalent state by sustaining H3K27 trimethylation therefore seems to be an important function of OTX2 in medulloblastoma, while other transcription factors might regulate the actual expression levels of these genes.

Chromosome 9q and 16q loss identified by genome-wide pooled-analysis are associated with tumor aggressiveness in patients with classic medulloblastoma

Medulloblastoma (MB) is one of the most aggressive pediatric brain tumor. We report genome-wide pooled-analysis of classic MB variant of patients over 3 years of age at diagnosis. We combined array comparative genomic hybridization (aCGH) results from experimental analysis (31 cases) with two public databases (55 cases) in a final evaluation of 86 MBs. The most common chromosome structural aberrations were gains of 17q (45.3%), 1q (22.1%), and losses of 8p (15.1%), 10q (19.8%), 17p (37.2%), and 16q (16.3%). Isochromome (17q) was observed in 29.1% MBs. A significant association between poor patients survival and losses of 9q (p < 0.0023), 10q (p < 0.012), and 16q (p < 0.036) was observed. Univariate analysis showed association of 9q loss (p < 0.008) and 16q loss (p = 0.05) with adverse overall survival (OS). Chromosome 6 monosomy was a protective event although statistically borderline (p = 0.066). After adjusting for confounding factors, a poor OS was found for patients whose tumor has 9q loss [hazard ratio (HR) = 3.97; p < 0.006) or 16q loss (HR = 2.41; p = 0.038). Our results highlight the importance of genomic studies in different MB histological variants and indicate a genotype-phenotype correlation.

Integrative genomic analysis of medulloblastoma identifies a molecular subgroup that drives poor clinical outcome

PURPOSE

Medulloblastomas are heterogeneous tumors that collectively represent the most common malignant brain tumor in children. To understand the molecular characteristics underlying their heterogeneity and to identify whether such characteristics represent risk factors for patients with this disease, we performed an integrated genomic analysis of a large series of primary tumors.

PATIENTS AND METHODS

We profiled the mRNA transcriptome of 194 medulloblastomas and performed high-density single nucleotide polymorphism array and miRNA analysis on 115 and 98 of these, respectively. Non-negative matrix factorization-based clustering of mRNA expression data was used to identify molecular subgroups of medulloblastoma; DNA copy number, miRNA profiles, and clinical outcomes were analyzed for each. We additionally validated our findings in three previously published independent medulloblastoma data sets.

RESULTS

Identified are six molecular subgroups of medulloblastoma, each with a unique combination of numerical and structural chromosomal aberrations that globally influence mRNA and miRNA expression. We reveal the relative contribution of each subgroup to clinical outcome as a whole and show that a previously unidentified molecular subgroup, characterized genetically by c-MYC copy number gains and transcriptionally by enrichment of photoreceptor pathways and increased miR-183∼96∼182 expression, is associated with significantly lower rates of event-free and overall survivals.

CONCLUSION

Our results detail the complex genomic heterogeneity of medulloblastomas and identify a previously unrecognized molecular subgroup with poor clinical outcome for which more effective therapeutic strategies should be developed.

Regulation of cell cycle genes and induction of senescence by overexpression of OTX2 in medulloblastoma cell lines

The transcription factor orthodenticle homeobox 2 (OTX2) has been implicated in the pathogenesis of medulloblastoma, as it is often highly expressed and sometimes amplified in these tumors. Little is known of the downstream pathways regulated by OTX2. We therefore generated MED8A and DAOY medulloblastoma cell lines with doxycycline-inducible OTX2 expression. In both cell lines, OTX2 inhibited proliferation and induced a senescence-like phenotype with senescence-associated β-galactosidase activity. Expression profiles of time series after OTX2 induction in MED8A showed early upregulation of cell cycle genes related to the G(2)-M phase, such as AURKA, CDC25C, and CCNG2. Paradoxically, G(1)-S phase genes such as MYC, CDK4, CDK6, CCND1, and CCND2 were strongly downregulated, in line with the observed G(1) arrest. ChIP-on-chip analyses of OTX2 binding to promoter regions in MED8A and DAOY showed a strong enrichment for binding to the G(2)-M genes, suggesting a direct activation. Their mRNA expression correlated with OTX2 expression in primary tumors, underscoring the in vivo relevance of this regulation. OTX2 induction activated the P53 pathway in MED8A, but not in DAOY, which carries a mutated P53 gene. In DAOY cells, senescence-associated secretory factors, such as interleukin-6 and insulin-like growth factor binding protein 7, were strongly upregulated after OTX2 induction. We hypothesize that the imbalance in cell cycle stimulation by OTX2 leads to cellular senescence either by activating the P53 pathway or through the induction of secretory factors. Our data indicate that OTX2 directly induces a series of cell cycle genes but requires cooperating genes for an oncogenic acceleration of the cell cycle.

Copy number and gene expression alterations in radiation-induced papillary thyroid carcinoma from chernobyl pediatric patients

BACKGROUND

Following exposure to radiation during the Chernobyl fallout tragedy, papillary thyroid carcinoma (PTC) increased significantly in individuals who were children at the time of the accident. We have used two high-throughput, whole genome platforms to analyze radiation-induced PTCs from pediatric patients from the Chernobyl region.

METHODS

We performed comparative genomic hybridization using Affymetrix 50K Mapping arrays and gene expression profiling on 10 pediatric post-Chernobyl PTCs obtained from patients living in the region. We performed an overlay analysis of these two data sets.

RESULTS

Many regions of copy number alterations (CNAs) were detected including novel regions that had never been associated with PTCs. Increases in copy numbers were consistently found on chromosomes 1p, 5p, 9q, 12q, 13q, 16p, 21q, and 22q. Deletions were observed less frequently and were mapped to 1q, 6q, 9q, 10q, 13q, 14q, 21q, and 22q. Gene expression analysis revealed that most of the altered genes were also perturbed in sporadic adult PTC; however, 141 gene expression changes were found to be unique to the post-Chernobyl tumors. The genes with the highest increases in expression that were novel to the pediatric post-Chernobyl tumors were TESC, PDZRN4, TRAa/TRDa, GABBR2, and CA12. The genes showing the largest expression decreases included PAPSS2, PDLIM3, BEXI, ANK2, SORBS2, and PPARGCIA. An overlay analysis of the gene expression and CNA profiles was then performed. This analysis identified genes showing both CNAs and concurrent gene expression alterations. Many of these are commonly seen in sporadic PTC such as SERPINA, COL8A, and PDX, while others were unique to the radiation-induced profiles including CAMK2N1, AK1, DHRS3, and PDE9A.

CONCLUSIONS

This type of analysis allows an assessment of gene expression changes that are associated with a physical mechanism. These genes and chromosomal regions are potential markers for radiation-induced PTC.

Combining chromosomal arm status and significantly aberrant genomic locations reveals new cancer subtypes

Many types of tumors exhibit characteristic chromosomal losses or gains, as well as local amplifications and deletions. Within any given tumor type, sample specific amplifications and deletions are also observed. Typically, a region that is aberrant in more tumors, or whose copy number change is stronger, would be considered as a more promising candidate to be biologically relevant to cancer. We sought for an intuitive method to define such aberrations and prioritize them. We define V, the "volume" associated with an aberration, as the product of three factors: (a) fraction of patients with the aberration, (b) the aberration's length and (c) its amplitude. Our algorithm compares the values of V derived from the real data to a null distribution obtained by permutations, and yields the statistical significance (p-value) of the measured value of V. We detected genetic locations that were significantly aberrant, and combine them with chromosomal arm status (gain/loss) to create a succinct fingerprint of the tumor genome. This genomic fingerprint is used to visualize the tumors, highlighting events that are co-occurring or mutually exclusive. We apply the method on three different public array CGH datasets of Medulloblastoma and Neuroblastoma, and demonstrate its ability to detect chromosomal regions that were known to be altered in the tested cancer types, as well as to suggest new genomic locations to be tested. We identified a potential new subtype of Medulloblastoma, which is analogous to Neuroblastoma type 1.

Ploidy status and copy number aberrations in primary glioblastomas defined by integrated analysis of allelic ratios, signal ratios and loss of heterozygosity using 500K SNP Mapping Arrays

Background

Genomic hybridization platforms, including BAC-CGH and genotyping arrays, have been used to estimate chromosome copy number (CN) in tumor samples by detecting the relative strength of genomic signal. The methods rely on the assumption that the predominant chromosomal background of the samples is diploid, an assumption that is frequently incorrect for tumor samples. In addition to generally greater resolution, an advantage of genotyping arrays over CGH arrays is the ability to detect signals from individual alleles, allowing estimation of loss-of-heterozygosity (LOH) and allelic ratios to enhance the interpretation of copy number alterations. Copy number events associated with LOH potentially have the same genetic consequences as deletions.

Results

We have utilized allelic ratios to detect patterns that are indicative of higher ploidy levels. An integrated analysis using allelic ratios, total signal and LOH indicates that many or most of the chromosomes from 24 glioblastoma tumors are in fact aneuploid. Some putative whole-chromosome losses actually represent trisomy, and many apparent sub-chromosomal losses are in fact relative losses against a triploid or tetraploid background.

Conclusion

These results suggest a re-interpretation of previous findings based only on total signal ratios. One interesting observation is that many single or multiple-copy deletions occur at common putative tumor suppressor sites subsequent to chromosomal duplication; these losses do not necessarily result in LOH, but nonetheless occur in conspicuous patterns. The 500 K Mapping array was also capable of detecting many sub-mega base losses and gains that were overlooked by CGH-BAC arrays, and was superior to CGH-BAC arrays in resolving regions of complex CN variation.

Biological background of pediatric medulloblastoma and ependymoma: a review from a translational research perspective

Survival rates of pediatric brain tumor patients have significantly improved over the years due to developments in diagnostic techniques, neurosurgery, chemotherapy, radiotherapy, and supportive care. However, brain tumors are still an important cause of cancer-related deaths in children. Prognosis is still highly dependent on clinical characteristics, such as the age of the patient, tumor type, stage, and localization, but increased knowledge about the genetic and biological features of these tumors is being obtained and might be useful to further improve outcome for these patients. It has become clear that the deregulation of signaling pathways essential in brain development, for example, sonic hedgehog (SHH), Wnt, and Notch pathways, plays an important role in pathogenesis and biological behavior, especially for medulloblastomas. More recently, data have become available about the cells of origin of brain tumors and the possible existence of brain tumor stem cells. Newly developed array-based techniques for studying gene expression, protein expression, copy number aberrations, and epigenetic events have led to the identification of other potentially important biological abnormalities in pediatric medulloblastomas and ependymomas.

Role of DEAD box 1 in retinoblastoma and neuroblastoma

Analysis of hereditary and nonhereditary retinoblastoma led to the formulation of the two-hit hypothesis of cancer in the early 1970s. The two-hit hypothesis was validated in the 1980s when both copies of the RB1 gene were shown to be mutated in hereditary and nonhereditary retinoblastoma. However, consistent genetic abnormalities other than RB1 mutations suggest that additional events may be required for the formation of these malignant tumors. For example, MYCN amplification has long been known to occur in both retinoblastoma and neuroblastoma tumors and is strongly associated with poor prognosis in neuroblastoma. The DEAD box gene, DEAD box 1 (DDX1), is often coamplified with MYCN in both these childhood tumors. Here, we examine possible roles for DDX1 overexpression in retinoblastoma and neuroblastoma.

Genome-wide screening of copy number alterations and LOH events in renal cell carcinomas and integration with gene expression profile

Background

Clear cell renal carcinoma (RCC) is the most common and invasive adult renal cancer. For the purpose of identifying RCC biomarkers, we investigated chromosomal regions and individual genes modulated in RCC pathology. We applied the dual strategy of assessing and integrating genomic and transcriptomic data, today considered the most effective approach for understanding genetic mechanisms of cancer and the most sensitive for identifying cancer-related genes.

Results

We performed the first integrated analysis of DNA and RNA profiles of RCC samples using Affymetrix technology. Using 100K SNP mapping arrays, we assembled a genome-wide map of DNA copy number alterations and LOH areas. We thus confirmed the typical genetic signature of RCC but also identified other amplified regions (e.g. on chr. 4, 11, 12), deleted regions (chr. 1, 9, 22) and LOH areas (chr. 1, 2, 9, 13). Simultaneously, using HG-U133 Plus 2.0 arrays, we identified differentially expressed genes (DEGs) in tumor vs. normal samples. Combining genomic and transcriptomic data, we identified 71 DEGs in aberrant chromosomal regions and observed, in amplified regions, a predominance of up-regulated genes (27 of 37 DEGs) and a trend to clustering. Functional annotation of these genes revealed some already implicated in RCC pathology and other cancers, as well as others that may be novel tumor biomarkers.

Conclusion

By combining genomic and transcriptomic profiles from a collection of RCC samples, we identified specific genomic regions with concordant alterations in DNA and RNA profiles and focused on regions with increased DNA copy number. Since the transcriptional modulation of up-regulated genes in amplified regions may be attributed to the genomic alterations characteristic of RCC, these genes may encode novel RCC biomarkers actively involved in tumor initiation and progression and useful in clinical applications.

Candidate glioblastoma development gene identification using concordance between copy number abnormalities and gene expression level changes

Copy number abnormalities (CNAs) in tumor cells are presumed to affect expression levels of genes located in region of abnormality. To investigate this relationship we have surveyed the losses, gains and amplifications in 30 glioblastomas using array comparative genome hybridization and compared these data with gene expression changes in the same tumors using the Affymetrix U133Plus2.0 oligonucleotide arrays. The two datasets were overlaid using our in-house overlay tool which highlights concordance between CNAs and expression level changes for the same tumors. In this survey we have highlighted genes frequently overexpressed in amplified regions on chromosomes 1, 4, 11, and 12 and have identified novel amplicons on these chromosomes. Deletions of specific regions on chromosomes 9, 10, 11, 14, and 15 have also been correlated with reduced gene expression in the regions of minimal overlap. In addition we describe a novel approach for comparing gene expression levels between tumors based on the presence or absence of chromosome CNAs. This genome wide screen provides an efficient and comprehensive survey of genes which potentially serve as the drivers for the CNAs in GBM.

Genome wide copy number abnormalities in pediatric medulloblastomas as assessed by array comparative genome hybridization

Array-based comparative genomic hybridization was used to characterize 22 medulloblastomas in order to precisely define genetic alterations in these malignant childhood brain tumors. The 17p(-)/17q(+) copy number abnormality (CNA), consistent with the formation of isochromosome 17q, was the most common event (8/22). Amplifications in this series included MYCL, MYCN and MYC previously implicated in medulloblastoma pathogenesis, as well as novel amplicons on chromosomes 2, 4, 11 and 12. Losses involving chromosomes 1, 2, 8, 10, 11, 16 and 19 and gains of chromosomes 4, 7, 8, 9 and 18 were seen in greater than 20% of tumors in this series. A homozygous deletion in 11p15 defines the minimal region of loss on this chromosome arm. In order to map the minimal regions involved in losses, gains and amplifications, we combined aCGH data from this series with that of two others obtained using the same RPCI BAC arrays. As a result of this combined analysis of 72 samples, we have defined specific regions on chromosomes 1, 8p, 10q, 11p and 16q which are frequently involved in CNAs in medulloblastomas. Using high density oligonucleotide expression arrays, candidate genes were identified within these consistently involved regions in a subset of the tumors.