Qing-Rong Chen

The MYCN oncogene is a direct target of miR-34a.

Loss of 1p36 heterozygosity commonly occurs with MYCN amplification in neuroblastoma tumors, and both are associated with an aggressive phenotype. Database searches identified five microRNAs that map to the commonly deleted region of 1p36 and we hypothesized that the loss of one or more of these microRNAs contributes to the malignant phenotype of MYCN-amplified tumors. By bioinformatic analysis, we identified that three out of the five microRNAs target MYCN and of these miR-34a caused the most significant suppression of cell growth through increased apoptosis and decreased DNA synthesis in neuroblastoma cell lines with MYCN amplification. Quantitative RT–PCR showed that neuroblastoma tumors with 1p36 loss expressed lower level of miR-34a than those with normal copies of 1p36. Furthermore, we demonstrated that MYCN is a direct target of miR-34a. Finally, using a series of mRNA expression profiling experiments, we identified other potential direct targets of miR-34a, and pathway analysis demonstrated that miR-34a suppresses cell-cycle genes and induces several neural-related genes. This study demonstrates one important regulatory role of miR-34a in cell growth and MYCN suppression in neuroblastoma.

Identification of FGFR4-activating mutations in human rhabdomyosarcomas that promote metastasis in xenotransplanted models

Rhabdomyosarcoma (RMS) is a childhood cancer originating from skeletal muscle, and patient survival is poor in the presence of metastatic disease. Few determinants that regulate metastasis development have been identified. The receptor tyrosine kinase FGFR4 is highly expressed in RMS tissue, suggesting a role in tumorigenesis, although its functional importance has not been defined. Here, we report the identification of mutations in FGFR4 in human RMS tumors that lead to its activation and present evidence that it functions as an oncogene in RMS. Higher FGFR4 expression in RMS tumors was associated with advanced-stage cancer and poor survival, while FGFR4 knockdown in a human RMS cell line reduced tumor growth and experimental lung metastases when the cells were transplanted into mice. Moreover, 6 FGFR4 tyrosine kinase domain mutations were found among 7 of 94 (7.5%) primary human RMS tumors. The mutants K535 and E550 increased autophosphorylation, Stat3 signaling, tumor proliferation, and metastatic potential when expressed in a murine RMS cell line. These mutants also transformed NIH 3T3 cells and led to an enhanced metastatic phenotype. Finally, murine RMS cell lines expressing the K535 and E550 FGFR4 mutants were substantially more susceptible to apoptosis in the presence of a pharmacologic FGFR inhibitor than the control cell lines expressing the empty vector or wild-type FGFR4. Together, our results demonstrate that mutationally activated FGFR4 acts as an oncogene, and these are what we believe to be the first known mutations in a receptor tyrosine kinase in RMS. These findings support the potential therapeutic targeting of FGFR4 in RMS.

Credentialing Preclinical Pediatric Xenograft Models Using Gene Expression and Tissue Microarray Analysis

Human tumor xenografts have been used extensively for rapid screening of the efficacy of anticancer drugs for the past 35 years. The selection of appropriate xenograft models for drug testing has been largely empirical and has not incorporated a similarity to the tumor type of origin at the molecular level. This study is the first comprehensive analysis of the transcriptome of a large set of pediatric xenografts, which are currently used for preclinical drug testing. Suitable models representing the tumor type of origin were identified. It was found that the characteristic expression patterns of the primary tumors were maintained in the corresponding xenografts for the majority of samples. Because a prerequisite for developing rationally designed drugs is that the target is expressed at the protein level, we developed tissue arrays from these xenografts and corroborated that high mRNA levels yielded high protein levels for two tested genes. The web database and availability of tissue arrays will allow for the rapid confirmation of the expression of potential targets at both the mRNA and the protein level for molecularly targeted agents. The database will facilitate the identification of tumor markers predictive of response to tested agents as well as the discovery of new molecular targets.

OmicCircos: A Simple-to-Use R Package for the Circular Visualization of Multidimensional Omics Data

Summary OmicCircos is an R software package used to generate high-quality circular plots for visualizing genomic variations, including mutation patterns, copy number variations (CNVs), expression patterns, and methylation patterns. Such variations can be displayed as scatterplot, line, or text-label figures. Relationships among genomic features in different chromosome positions can be represented in the forms of polygons or curves. Utilizing the statistical and graphic functions in an R/Bioconductor environment, OmicCircos performs statistical analyses and displays results using cluster, boxplot, histogram, and heatmap formats. In addition, OmicCircos offers a number of unique capabilities, including independent track drawing for easy modification and integration, zoom functions, link-polygons, and position-independent heatmaps supporting detailed visualization. Availability and Implementation OmicCircos is available through Bioconductor at http://www.bioconductor.org/packages/devel/bioc/html/OmicCircos.html. An extensive vignette in the package describes installation, data formatting, and workflow procedures. The software is open source under the Artistic—2.0 license.

BRD4 is a histone acetyltransferase that evicts nucleosomes from chromatin

Bromodomain protein 4 (BRD4) is a chromatin-binding protein implicated in cancer and autoimmune diseases that functions as a scaffold for transcription factors at promoters and super-enhancers. Although chromatin decompaction and transcriptional activation of target genes are associated with BRD4 binding, the mechanisms involved are unknown. We report that BRD4 is a histone acetyltransferase (HAT) that acetylates histones H3 and H4 with a pattern distinct from those of other HATs. Both mouse and human BRD4 have intrinsic HAT activity. Importantly, BRD4 acetylates H3 K122, a residue critical for nucleosome stability, thus resulting in nucleosome eviction and chromatin decompaction. Nucleosome clearance by BRD4 occurs genome wide, including at its targets MYC, FOS and AURKB (Aurora B kinase), resulting in increased transcription. These findings suggest a model wherein BRD4 actively links chromatin structure and transcription: it mediates chromatin decompaction by acetylating and evicting nucleosomes at target genes, thereby activating transcription.

Altered expression of cell cycle genes distinguishes aggressive neuroblastoma

In this study, gene expression profiling was performed on 103 neuroblastoma (NB) tumors, stages 1–4 with and without MYCN amplification, using cDNA microarrays containing 42 578 elements. Using principal component analysis (PCA) to analyse the relationships among these samples, we confirm that the global patterns of gene expression reflect the phenotype of the tumors. To explore the biological processes that may contribute to increasing aggressive phenotype of the tumors, we utilized a statistical approach based on PCA. We identified a specific subset of the cell cycle and/or chromosome segregation genes that distinguish stage 4 NB tumors from all lower stage tumors, including stage 3. Furthermore, the control of the kinetochore assembly emerges from the Gene Ontology analysis as one of the key biological processes associated with an aggressive NB phenotype. Finally, we establish that these genes are further upregulated in the most aggressive MYCN-amplified tumors.

microRNA Profiling Identifies Cancer-Specific and Prognostic Signatures in Pediatric Malignancies

Purpose: microRNAs have been shown to be involved in different human cancers. We therefore have performed expression profiles on a panel of pediatric tumors to identify cancer-specific microRNAs. We also investigated if microRNAs are coregulated with their host gene. Experimental Design: We performed parallel microRNAs and mRNA expression profiling on 57 tumor xenografts and cell lines representing 10 different pediatric solid tumors using microarrays. For those microRNAs that map to their host mRNA, we calculated correlations between them. Results: We found that the majority of cancer types clustered together based on their global microRNA expression profiles by unsupervised hierarchical clustering. Fourteen microRNAs were significantly differentially expressed between rhabdomyosarcoma and neuroblastoma, and 8 of them were validated in independent patient tumor samples. Exploration of the expression of microRNAs in relationship with their host genes showed that the expression for 43 of 68 (63%) microRNAs located inside known coding genes was significantly correlated with that of their host genes. Among these 43 microRNAs, 5 of 7 microRNAs in the OncomiR-1 cluster correlated significantly with their host gene MIRHG1 (P < 0.01). In addition, high expression of MIRHG1 was significantly associated with high stage and MYCN amplification in neuroblastoma tumors, and the expression level of MIRHG1 could predict the outcome of neuroblastoma patients independently from the current neuroblastoma risk-stratification in two independent patient cohorts. Conclusion: Pediatric cancers express cancer-specific microRNAs. The high expression of the OncomiR-1 host gene MIRHG1 correlates with poor outcome for patients with neuroblastoma, indicating important oncogenic functions of this microRNA cluster in neuroblastoma biology. (Clin Cancer Res 2009;15(17):5560–8)

Global Genomic and Proteomic Analysis Identifies Biological Pathways Related to High-Risk Neuroblastoma

Neuroblastoma (NB) is a heterogeneous pediatric tumor. To better understand the biological pathways involved in the development of high-risk neuroblastoma, we performed parallel global protein and mRNA expression profiling on NB tumors of stage 4 MYCN-amplified (4+) and stage 1 MYCN-not-amplified (1-) using isotope-coded affinity tags (ICAT) and Affymetrix U133plus2 microarray, respectively. A total of 1461 proteins represented by 2 or more peptides were identified from the quantitative ICAT analysis, of which 433 and 130 proteins are up- or down-regulated, respectively, in 4+ tumor compared to the 1- tumor. Pathway analysis of the differentially expressed proteins showed the enrichment of glycolysis, DNA replication and cell cycle processes in the up-regulated proteins and cell adhesion, nervous system development and cell differentiation processes in the down-regulated proteins in 4+ tumor; suggesting a less mature neural and a more invasive phenotype of 4+ tumor. Myc targets and ribosomal proteins are overrepresented in the 4+ tumors as expected; functional gene sets reported to be enriched in neural and embryonic stem cells are significantly enriched in the 4+ tumor, indicating the existence of a stemness signature in MYCN-amplified stage 4 tumor. In addition, protein and mRNA expression are moderately correlated (r = 0.51, p < 0.0001), as approximately half of the up-regulated proteins in 4+ tumor have elevated mRNA level (n = 208), and one-third of down-regulated proteins have lower mRNA expression (n = 47). Further biological network analysis revealed that the differentially expressed proteins closely interact with other proteins of known networks; the important role of MYCN is confirmed and other transcription factors identified in the network may have potential roles in the biology of NB tumor. We used global genomic and proteomic analysis to identify biologically relevant proteins and pathways important to NB progression and development that may provide new insights into the biology of advanced neuroblastoma.

An integrated cross-platform prognosis study on neuroblastoma patients

There have been several reports about the potential for predicting prognosis of neuroblastoma patients using microarray gene expression profiling of the tumors. However these studies have revealed an apparent diversity in the identity of the genes in their predictive signatures. To test the contribution of the platform to this discrepancy we applied the z-scoring method to minimize the impact of platform and combine gene expression profiles of neuroblastoma (NB) tumors from two different platforms, cDNA and Affymetrix. A total of 12442 genes were common to both cDNA and Affymetrix arrays in our data set. Two-way ANOVA analysis was applied to the combined data set for assessing the relative effect of prognosis and platform on gene expression. We found that 26.6% (3307) of the genes had significant impact on survival. There was no significant impact of microarray platform on expression after application of z-scoring standardization procedure. Artificial neural network (ANN) analysis of the combined data set in a leave-one-out prediction strategy correctly predicted the outcome for 90% of the samples. Hierarchical clustering analysis using the top-ranked 160 genes showed the great separation of two clusters, and the majority of matched samples from the different platforms were clustered next to each other. The ANN classifier trained with our combined cross-platform data for these 160 genes could predict the prognosis of 102 independent test samples with 71% accuracy. Furthermore it correctly predicted the outcome for 85/102 (83%) NB patients through the leave-one-out cross-validation approach. Our study showed that gene expression studies performed in different platforms could be integrated for prognosis analysis after removing variation resulting from different platforms.

Dysregulation of Ezrin Phosphorylation Prevents Metastasis and Alters Cellular Metabolism in Osteosarcoma

Ezrin links the plasma membrane to the actin cytoskeleton where it plays a pivotal role in the metastatic progression of several human cancers; however, the precise mechanistic basis for its role remains unknown. Here, we define transitions between active (phosphorylated open) and inactive (dephosphorylated closed) forms of Ezrin that occur during metastatic progression in osteosarcoma. In our evaluation of these conformations we expressed C-terminal mutant forms of Ezrin that are open (phosphomimetic T567D) or closed (phosphodeficient T567A) and compared their biologic characteristics to full-length wild-type Ezrin in osteosarcoma cells. Unexpectedly, cells expressing open, active Ezrin could form neither primary orthotopic tumors nor lung metastases. In contrast, cells expressing closed, inactive Ezrin were also deficient in metastasis but were unaffected in their capacity for primary tumor growth. By imaging single metastatic cells in the lung, we found that cells expressing either open or closed Ezrin displayed increased levels of apoptosis early after their arrival in the lung. Gene expression analysis suggested dysregulation of genes that are functionally linked to carbohydrate and amino acid metabolism. In particular, cells expressing closed, inactive Ezrin exhibited reduced lactate production and basal or ATP-dependent oxygen consumption. Collectively, our results suggest that dynamic regulation of Ezrin phosphorylation at amino acid T567 that controls structural transitions of this protein plays a pivotal role in tumor progression and metastasis, possibly in part by altering cellular metabolism.