Jae-Geun Yoon

The SOX2 response program in glioblastoma multiforme: an integrated ChIP-seq, expression microarray, and microRNA analysis

BackgroundSOX2 is a key gene implicated in maintaining the stemness of embryonic and adult stem cells. SOX2 appears to re-activate in several human cancers including glioblastoma multiforme (GBM), however, the detailed response program of SOX2 in GBM has not yet been defined.ResultsWe show that knockdown of the SOX2 gene in LN229 GBM cells reduces cell proliferation and colony formation. We then comprehensively characterize the SOX2 response program by an integrated analysis using several advanced genomic technologies including ChIP-seq, microarray profiling, and microRNA sequencing. Using ChIP-seq technology, we identified 4883 SOX2 binding regions in the GBM cancer genome. SOX2 binding regions contain the consensus sequence wwTGnwTw that occurred 3931 instances in 2312 SOX2 binding regions. Microarray analysis identified 489 genes whose expression altered in response to SOX2 knockdown. Interesting findings include that SOX2 regulates the expression of SOX family proteins SOX1 and SOX18, and that SOX2 down regulates BEX1 (brain expressed X-linked 1) and BEX2 (brain expressed X-linked 2), two genes with tumor suppressor activity in GBM. Using next generation sequencing, we identified 105 precursor microRNAs (corresponding to 95 mature miRNAs) regulated by SOX2, including down regulation of miR-143, -145, -253-5p and miR-452. We also show that miR-145 and SOX2 form a double negative feedback loop in GBM cells, potentially creating a bistable system in GBM cells.ConclusionsWe present an integrated dataset of ChIP-seq, expression microarrays and microRNA sequencing representing the SOX2 response program in LN229 GBM cells. The insights gained from our integrated analysis further our understanding of the potential actions of SOX2 in carcinogenesis and serves as a useful resource for the research community.

A CD133-related gene expression signature identifies an aggressive glioblastoma subtype with excessive mutations

Cancer cells are heterogeneous and, it has been proposed, fall into at least two classes: the tumor-initiating cancer stem cells (CSC) and the more differentiated tumor cells. The transmembrane protein CD133 has been widely used to isolate putative CSC populations in several cancer types, but its validity as a CSC marker and hence its clinical ramifications remain controversial. Here, we conducted transcriptomic profiling of sorted CD133+ and CD133− cells from human glioblastoma multiforme (GBM) and, by subtractive analysis, established a CD133 gene expression signature composed of 214 differentially expressed genes. Extensive computational comparisons with a compendium of published gene expression profiles reveal that the CD133 gene signature transcriptionally resembles human ES cells and in vitro cultured GBM stem cells, and this signature successfully distinguishes GBM from lower-grade gliomas. More importantly, the CD133 gene signature identifies an aggressive subtype of GBM seen in younger patients with shorter survival who bear excessive genomic mutations as surveyed through the Cancer Genome Atlas Network GBM mutation spectrum. Furthermore, the CD133 gene signature distinguishes higher-grade breast and bladder cancers from their lower-grade counterparts. Our systematic analysis provides molecular and genetic support for the stem cell-like nature of CD133+ cells and an objective means for evaluating cancer aggressiveness.

Comprehensive Analysis of MGMT Promoter Methylation: Correlation with MGMT Expression and Clinical Response in GBM

O6-methylguanine DNA-methyltransferase (MGMT) promoter methylation has been identified as a potential prognostic marker for glioblastoma patients. The relationship between the exact site of promoter methylation and its effect on gene silencing, and the patients subsequent response to therapy, is still being defined. The aim of this study was to comprehensively characterize cytosine-guanine (CpG) dinucleotide methylation across the entire MGMT promoter and to correlate individual CpG site methylation patterns to mRNA expression, protein expression, and progression-free survival. To best identify the specific MGMT promoter region most predictive of gene silencing and response to therapy, we determined the methylation status of all 97 CpG sites in the MGMT promoter in tumor samples from 70 GBM patients using quantitative bisulfite sequencing. We next identified the CpG site specific and regional methylation patterns most predictive of gene silencing and improved progression-free survival. Using this data, we propose a new classification scheme utilizing methylation data from across the entire promoter and show that an analysis based on this approach, which we call 3R classification, is predictive of progression-free survival (HR  = 5.23, 95% CI [2.089–13.097], p<0.0001). To adapt this approach to the clinical setting, we used a methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) test based on the 3R classification and show that this test is both feasible in the clinical setting and predictive of progression free survival (HR  = 3.076, 95% CI [1.301–7.27], p = 0.007). We discuss the potential advantages of a test based on this promoter-wide analysis and compare it to the commonly used methylation-specific PCR test. Further prospective validation of these two methods in a large independent patient cohort will be needed to confirm the added value of promoter wide analysis of MGMT methylation in the clinical setting.

Massively Parallel Signature Sequencing and Bioinformatics Analysis Identifies Up-Regulation of TGFBI and SOX4 in Human Glioblastoma

Background A comprehensive network-based understanding of molecular pathways abnormally altered in glioblastoma multiforme (GBM) is essential for developing effective therapeutic approaches for this deadly disease. Methodology/Principal Findings Applying a next generation sequencing technology, massively parallel signature sequencing (MPSS), we identified a total of 4535 genes that are differentially expressed between normal brain and GBM tissue. The expression changes of three up-regulated genes, CHI3L1, CHI3L2, and FOXM1, and two down-regulated genes, neurogranin and L1CAM, were confirmed by quantitative PCR. Pathway analysis revealed that TGF- β pathway related genes were significantly up-regulated in GBM tumor samples. An integrative pathway analysis of the TGF β signaling network identified two alternative TGF−β signaling pathways mediated by SOX4 (sex determining region Y-box 4) and TGFBI (Transforming growth factor beta induced). Quantitative RT-PCR and immunohistochemistry staining demonstrated that SOX4 and TGFBI expression is elevated in GBM tissues compared with normal brain tissues at both the RNA and protein levels. In vitro functional studies confirmed that TGFBI and SOX4 expression is increased by TGF- β stimulation and decreased by a specific inhibitor of TGF- β receptor 1 kinase. Conclusions/Significance Our MPSS database for GBM and normal brain tissues provides a useful resource for the scientific community. The identification of non-SMAD mediated TGF−β signaling pathways acting through SOX4 and TGFBI (GENE ID:7045) in GBM indicates that these alternative pathways should be considered, in addition to the canonical SMAD mediated pathway, in the development of new therapeutic strategies targeting TGF−β signaling in GBM. Finally, the construction of an extended TGF- β signaling network with overlaid gene expression changes between GBM and normal brain extends our understanding of the biology of GBM.

Dickkopf-1 is an epigenetically silenced candidate tumor suppressor gene in medulloblastoma

Medulloblastoma is a heterogeneous pediatric brain tumor with significant therapy-related morbidity, its five-year survival rates ranging from 30% to 70%. Improvement in diagnosis and therapy requires better understanding of medulloblastoma pathology. We used whole-genome microarray analysis to identify putative tumor suppressor genes silenced by epigenetic mechanisms in medulloblastoma. This analysis yielded 714 up-regulated genes in immortalized medulloblastoma cell line D283 on treatment with histone deacetylase (HDAC) inhibitor trichostatin A (TSA). Dickkopf-1 (DKK1), a Wnt antagonist, was found to be up-regulated on HDAC inhibition. We examined DKK1 expression in primary medulloblastoma cells and patient samples by reverse transcriptase PCR and found it to be significantly down-regulated relative to normal cerebellum. Transfection of a DKK1 gene construct into D283 cell lines suppressed medulloblastoma tumor growth in colony focus assays by 60% (P < 0.001). In addition, adenoviral vector-mediated expression of DKK1 in medulloblastoma cells increased apoptosis fourfold (P < 0.001). These data reveal that inappropriate histone modifications might deregulate DKK1 expression in medulloblastoma tumorigenesis and block its tumor-suppressive activity.

Exploration of the gene fusion landscape of glioblastoma using transcriptome sequencing and copy number data

BackgroundRNA-seq has spurred important gene fusion discoveries in a number of different cancers, including lung, prostate, breast, brain, thyroid and bladder carcinomas. Gene fusion discovery can potentially lead to the development of novel treatments that target the underlying genetic abnormalities.ResultsIn this study, we provide comprehensive view of gene fusion landscape in 185 glioblastoma multiforme patients from two independent cohorts. Fusions occur in approximately 30-50% of GBM patient samples. In the Ivy Center cohort of 24 patients, 33% of samples harbored fusions that were validated by qPCR and Sanger sequencing. We were able to identify high-confidence gene fusions from RNA-seq data in 53% of the samples in a TCGA cohort of 161 patients. We identified 13 cases (8%) with fusions retaining a tyrosine kinase domain in the TCGA cohort and one case in the Ivy Center cohort. Ours is the first study to describe recurrent fusions involving non-coding genes. Genomic locations 7p11 and 12q14-15 harbor majority of the fusions. Fusions on 7p11 are formed in focally amplified EGFR locus whereas 12q14-15 fusions are formed by complex genomic rearrangements. All the fusions detected in this study can be further visualized and analyzed using our website: http://ivygap.swedish.org/fusions.ConclusionsOur study highlights the prevalence of gene fusions as one of the major genomic abnormalities in GBM. The majority of the fusions are private fusions, and a minority of these recur with low frequency. A small subset of patients with fusions of receptor tyrosine kinases can benefit from existing FDA approved drugs and drugs available in various clinical trials. Due to the low frequency and rarity of clinically relevant fusions, RNA-seq of GBM patient samples will be a vital tool for the identification of patient-specific fusions that can drive personalized therapy.

Next Generation Sequencing Analysis of miRNAs: MiR-127-3p Inhibits Glioblastoma Proliferation and Activates TGF-β Signaling by Targeting SKI

Glioblastoma (GBM) proliferation is a multistep process during which the expression levels of many genes that control cell proliferation, cell death, and genetic stability are altered. MicroRNAs (miRNAs) are emerging as important modulators of cellular signaling, including cell proliferation in cancer. In this study, using next generation sequencing analysis of miRNAs, we found that miR-127-3p was downregulated in GBM tissues compared with normal brain tissues; we validated this result by RT-PCR. We further showed that DNA demethylation and histone deacetylase inhibition resulted in downregulation of miR-127-3p. We demonstrated that miR-127-3p overexpression inhibited GBM cell growth by inducing G1-phase arrest both in vitro and in vivo. We showed that miR-127-3p targeted SKI (v-ski sarcoma viral oncogene homolog [avian]), RGMA (RGM domain family, member A), ZWINT (ZW10 interactor, kinetochore protein), SERPINB9 (serpin peptidase inhibitor, clade B [ovalbumin], member 9), and SFRP1 (secreted frizzled-related protein 1). Finally, we found that miR-127-3p suppressed GBM cell growth by inhibiting tumor-promoting SKI and activating the tumor suppression effect of transforming growth factor-β (TGF-β) signaling. This study showed, for the first time, that miR-127-3p and its targeted gene SKI, play important roles in GBM and may serve as potential targets for GBM therapy.

Transcriptional profiling of lymphoblast lines from subjects with panic disorder.

In attempts to isolate genetic vulnerability factors for panic disorder (PD), a number of investigators have used genome‐wide linkage or association analyses. But these attempts have been only modestly successful which suggests that alternative approaches may be needed to define the biology of PD. Therefore, using recently developed genome‐wide gene expression profiling, we explored whether transcriptional signatures associated with PD are present in lymphoblast cell line. The expression of 2,469 transcripts in lymphoblast cell lines from 16 subjects was arithmetically increased in every line and significantly increased overall and 354 transcripts was arithmetically decreased in every cell line and significantly decreased overall as compared to those lymphoblast lines from 17 subjects without a history of behavioral illness. Further sex specific analyses showed that in those 10 lines derived from female probands, the expression of a further 67 transcripts was arithmetically increased in every line and significantly increased overall and a further 332 transcripts was arithmetically decreased in every cell line and significantly decreased. Conversely, in cell lines from the six male probands, the expression of an additional 212 was arithmetically increased in every line and significantly increased overall and a further 332 transcripts was arithmetically decreased in every cell line. We conclude that lymphoblast cell lines derived from subjects with PD have significant, partially sex dependent changes in gene transcription. Further studies are necessary to correlate these changes in these hemopoetically derived cells with those changes postulated to occur in the CNS in association with PD.

Epigenetic regulation of wnt pathway antagonists in human glioblastoma multiforme.

Epigenetic inactivation of tumor suppressor genes is common in human cancer. Using a large-scale whole-genome approach in an earlier study, the authors identified epigenetically silenced genes with potential tumor suppressor function in glioblastoma (GBM). Three genes identified in this analysis-DKK1, SFRP1, and WIF1-are potent inhibitors of the Wnt signal transduction pathway. Here, the authors confirm decreased expression of these genes in GBM tumor tissue samples relative to nontumor brain tissue samples using real-time PCR. They then show that expression of all 3 genes is restored in T98 GBM cells by treatment with the histone deacetylase inhibitor Trichostatin A (TSA), but only DKK1 expression is restored by treatment with the demethylating agent 5-azacytidine. Bisulfite sequencing did not reveal significant methylation in the promoter region of DKK1, whereas histone acetylation and chromatin accessibility increased significantly for all 3 genes after TSA treatment. Ectopic expression of DKK1 significantly reduces colony formation and increases chemotherapy-induced apoptosis in T98 cells. Ectopic expression of the canonical Wnt pathway inhibitors WIF1 and SFRP1 shows a relative lack of response. Chronic Wnt3a stimulation only partially reverses growth suppression after DKK1 reexpression, whereas a specific inhibitor of the JNK pathway significantly reverses the effect of DKK1 reexpression on colony formation and apoptosis in T98 cells. These results support a potential growth-suppressive function for epigenetically silenced DKK1 in GBM and suggest that DKK1 restoration could modulate Wnt signaling through both canonical and noncanonical pathways.

Transcriptional profiling of subjects from the Iowa adoption studies.

Transcriptional profiling has been used to identify gene expression patterns indicative of general medical illnesses such as atherosclerosis. However, whether these methods can identify common psychiatric disorders has not been established. To answer this question with respect to nicotine use, we used genome‐wide expression profiling lymphoblast cell lines from six actively smoking Iowa Adoption Studies (IAS) subjects and nine “clean” control subjects, followed by real‐time PCR (RT‐PCR) of gene expression patterns in lymphoblast derived RNA from 94 subjects in the IAS. As compared to those from controls without a history of smoking (n = 9), the expression levels of 579 of 29,098 genes were significantly up‐regulated and expression levels of 584 of 29,098 genes were significantly down‐regulated in lymphoblast lines from currently smoking subjects (n = 6). RT‐PCR confirmation of four select RNA levels confirmed the validity of the overall profile and revealed highly significant relationships between the expression of some of these transcripts and (1) major depression, (2) antisocial personality, (3) nicotine dependence, and (4) cannabis dependence. We conclude that the use of expression profiling may contribute significant insights into the biology of complex behavioral disorders.