Incorporating genomic, transcriptomic and clinical data: a prognostic and stem cell-like MYC and PRC imbalance in high-risk neuroblastoma.

BACKGROUND: Previous studies suggested that cancer cells possess traits reminiscent of the biological mechanisms ascribed to normal embryonic stem cells (ESCs) regulated by MYC and Polycomb repressive complex 2 (PRC2). Several poorly differentiated adult tumors showed preferentially high expression levels in targets of MYC, coincident with low expression levels in targets of PRC2. This paper will reveal this ESC-like cancer signature in high-risk neuroblastoma (HR-NB), the most common extracranial solid tumor in children. METHODS: We systematically assembled genomic variants, gene expression changes, priori knowledge of gene functions, and clinical outcomes to identify prognostic multigene signatures. First, we assigned a new, individualized prognostic index using the relative expressions between the poor- and good-outcome signature genes. We then characterized HR-NB aggressiveness beyond these prognostic multigene signatures through the imbalanced effects of MYC and PRC2 signaling. We further analyzed Retinoic acid (RA)-induced HR-NB cells to model tumor cell differentiation. Finally, we performed in vitro validation on ZFHX3, a cell differentiation marker silenced by PRC2, and compared cell morphology changes before and after blocking PRC2 in HR-NB cells. RESULTS: A significant concurrence existed between exons with verified variants and genes showing MYCN-dependent expression in HR-NB. From these biomarker candidates, we identified two novel prognostic gene-set pairs with multi-scale oncogenic defects. Intriguingly, MYC targets over-represented an unfavorable component of the identified prognostic signatures while PRC2 targets over-represented a favorable component. The cell cycle arrest and neuronal differentiation marker ZFHX3 was identified as one of PRC2-silenced tumor suppressor candidates. Blocking PRC2 reduced tumor cell growth and increased the mRNA expression levels of ZFHX3 in an early treatment stage. This hypothesis-driven systems bioinformatics work offered novel insights into the PRC2-mediated tumor cell growth and differentiation in neuroblastoma, which may exert oncogenic effects together with MYC regulation. CONCLUSION: Our results propose a prognostic effect of imbalanced MYC and PRC2 moderations in pediatric HR-NB for the first time. This study demonstrates an incorporation of genomic landscapes and transcriptomic profiles into the hypothesis-driven precision prognosis and biomarker discovery. The application of this approach to neuroblastoma, as well as other cancer more broadly, could contribute to reduced relapse and mortality rates in the long term.

A c-Myc-regulated stem cell-like signature in high-risk neuroblastoma: A systematic discovery (Target neuroblastoma ESC-like signature).

c-Myc dysregulation is hypothesized to account for the 'stemness' - self-renewal and pluripotency - shared between embryonic stem cells (ESCs) and adult aggressive tumours. High-risk neuroblastoma (HR-NB) is the most frequent, aggressive, extracranial solid tumour in childhood. Using HR-NB as a platform, we performed a network analysis of transcriptome data and presented a c-Myc subnetwork enriched for genes previously reported as ESC-like cancer signatures. A subsequent drug-gene interaction analysis identified a pharmacogenomic agent that preferentially interacted with this HR-NB-specific, ESC-like signature. This agent, Roniciclib (BAY 1000394), inhibited neuroblastoma cell growth and induced apoptosis in vitro. It also repressed the expression of the oncogene c-Myc and the neural ESC marker CDK2 in vitro, which was accompanied by altered expression of the c-Myc-targeted cell cycle regulators CCND1, CDKN1A and CDKN2D in a time-dependent manner. Further investigation into this HR-NB-specific ESC-like signature in 295 and 243 independent patients revealed and validated the general prognostic index of CDK2 and CDKN3 compared with CDKN2D and CDKN1B. These findings highlight the very potent therapeutic benefits of Roniciclib in HR-NB through the targeting of c-Myc-regulated, ESC-like tumorigenesis. This work provides a hypothesis-driven systems computational model that facilitates the translation of genomic and transcriptomic signatures to molecular mechanisms underlying high-risk tumours.

Correction: Foxf Genes Integrate Tbx5 and Hedgehog Pathways in the Second Heart Field for Cardiac Septation.

[This corrects the article DOI: 10.1371/journal.pgen.1004604.].

The Cardiac TBX5 Interactome Reveals a Chromatin Remodeling Network Essential for Cardiac Septation.

Human mutations in the cardiac transcription factor gene TBX5 cause congenital heart disease (CHD), although the underlying mechanism is unknown. We report characterization of the endogenous TBX5 cardiac interactome and demonstrate that TBX5, long considered a transcriptional activator, interacts biochemically and genetically with the nucleosome remodeling and deacetylase (NuRD) repressor complex. Incompatible gene programs are repressed by TBX5 in the developing heart. CHD mis-sense mutations that disrupt the TBX5-NuRD interaction cause depression of a subset of repressed genes. Furthermore, the TBX5-NuRD interaction is required for heart development. Phylogenetic analysis showed that the TBX5-NuRD interaction domain evolved during early diversification of vertebrates, simultaneous with the evolution of cardiac septation. Collectively, this work defines a TBX5-NuRD interaction essential to cardiac development and the evolution of the mammalian heart, and when altered may contribute to human CHD.

Identification of epigenetic modifications that contribute to pathogenesis in therapy-related AML: Effective integration of genome-wide histone modification with transcriptional profiles.

BACKGROUND: Therapy-related, secondary acute myeloid leukemia (t-AML) is an increasingly frequent complication of intensive chemotherapy. This malignancy is often characterized by abnormalities of chromosome 7, including large deletions or chromosomal loss. A variety of studies suggest that decreased expression of the EZH2 gene located at 7q36.1 is critical in disease pathogenesis. This histone methyltransferase has been implicated in transcriptional repression through modifying histone H3 on lysine 27 (H3k27). However, the critical target genes of EZH2 and their regulatory roles remain unclear. METHOD: To characterize the subset of EZH2 target genes that might contribute to t-AML pathogenesis, we developed a novel computational analysis to integrate tissue-specific histone modifications and genome-wide transcriptional regulation. Initial integrative analysis utilized a novel "seq2gene" strategy to explore largely the target genes of chromatin immuneprecipitation sequencing (ChIP-seq) enriched regions. By combining seq2gene with our Phenotype-Genotype-Network (PGNet) algorithm, we enriched genes with similar expression profiles and genomic or functional characteristics into "biomodules". RESULTS: Initial studies identified SEMA3A (semaphoring 3A) as a novel oncogenic candidate that is regulated by EZH2-silencing, using data derived from both normal and leukemic cell lines as well as murine cells deficient in EZH2. A microsatellite marker at the SEMA3A promoter has been associated with chemosensitivity and radiosensitivity. Notably, our subsequent studies in primary t-AML demonstrate an expected up-regulation of SEMA3A that is EZH2-modulated. Furthermore, we have identified three biomodules that are co-expressed with SEMA3A and up-regulated in t-AML, one of which consists of previously characterized EZH2-repressed gene targets. The other two biomodules include MAPK8 and TATA box targets. Together, our studies suggest an important role for EZH2 targets in t-AML pathogenesis that warrants further study. CONCLUSION: These developed computational algorithms and systems biology strategies will enhance the knowledge discovery and hypothesis-driven analysis of multiple next generation sequencing data, for t-AML and other complex diseases.

Seq2pathway: an R/Bioconductor package for pathway analysis of next-generation sequencing data.

UNLABELLED: Seq2pathway is an R/Python wrapper for pathway (or functional gene-set) analysis of genomic loci, adapted for advances in genome research. Seq2pathway associates the biological significance of genomic loci with their target transcripts and then summarizes the quantified values on the gene-level into pathway scores. It is designed to isolate systematic disturbances and common biological underpinnings from next-generation sequencing (NGS) data. Seq2pathway offers Bioconductor users enhanced capability in discovering collective pathway effects caused by both coding genes and cis-regulation of non-coding elements. AVAILABILITY AND IMPLEMENTATION: The package is freely available at http://www.bioconductor.org/packages/release/bioc/html/seq2pathway.html. CONTACT: xyang2@uchicago.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Systematic computation with functional gene-sets among leukemic and hematopoietic stem cells reveals a favorable prognostic signature for acute myeloid leukemia.

BACKGROUND: Genes that regulate stem cell function are suspected to exert adverse effects on prognosis in malignancy. However, diverse cancer stem cell signatures are difficult for physicians to interpret and apply clinically. To connect the transcriptome and stem cell biology, with potential clinical applications, we propose a novel computational "gene-to-function, snapshot-to-dynamics, and biology-to-clinic" framework to uncover core functional gene-sets signatures. This framework incorporates three function-centric gene-set analysis strategies: a meta-analysis of both microarray and RNA-seq data, novel dynamic network mechanism (DNM) identification, and a personalized prognostic indicator analysis. This work uses complex disease acute myeloid leukemia (AML) as a research platform. RESULTS: We introduced an adjustable "soft threshold" to a functional gene-set algorithm and found that two different analysis methods identified distinct gene-set signatures from the same samples. We identified a 30-gene cluster that characterizes leukemic stem cell (LSC)-depleted cells and a 25-gene cluster that characterizes LSC-enriched cells in parallel; both mark favorable-prognosis in AML. Genes within each signature significantly share common biological processes and/or molecular functions (empirical p = 6e-5 and 0.03 respectively). The 25-gene signature reflects the abnormal development of stem cells in AML, such as AURKA over-expression. We subsequently determined that the clinical relevance of both signatures is independent of known clinical risk classifications in 214 patients with cytogenetically normal AML. We successfully validated the prognosis of both signatures in two independent cohorts of 91 and 242 patients respectively (log-rank p < 0.0015 and 0.05; empirical p < 0.015 and 0.08). CONCLUSION: The proposed algorithms and computational framework will harness systems biology research because they efficiently translate gene-sets (rather than single genes) into biological discoveries about AML and other complex diseases.

Foxf genes integrate tbx5 and hedgehog pathways in the second heart field for cardiac septation.

The Second Heart Field (SHF) has been implicated in several forms of congenital heart disease (CHD), including atrioventricular septal defects (AVSDs). Identifying the SHF gene regulatory networks required for atrioventricular septation is therefore an essential goal for understanding the molecular basis of AVSDs. We defined a SHF Hedgehog-dependent gene regulatory network using whole genome transcriptional profiling and GLI-chromatin interaction studies. The Forkhead box transcription factors Foxf1a and Foxf2 were identified as SHF Hedgehog targets. Compound haploinsufficiency for Foxf1a and Foxf2 caused atrioventricular septal defects, demonstrating the biological relevance of this regulatory network. We identified a Foxf1a cis-regulatory element that bound the Hedgehog transcriptional regulators GLI1 and GLI3 and the T-box transcription factor TBX5 in vivo. GLI1 and TBX5 synergistically activated transcription from this cis-regulatory element in vitro. This enhancer drove reproducible expression in vivo in the posterior SHF, the only region where Gli1 and Tbx5 expression overlaps. Our findings implicate Foxf genes in atrioventricular septation, describe the molecular underpinnings of the genetic interaction between Hedgehog signaling and Tbx5, and establish a molecular model for the selection of the SHF gene regulatory network for cardiac septation.

Network models of genome-wide association studies uncover the topological centrality of protein interactions in complex diseases.

BACKGROUND: While genome-wide association studies (GWAS) of complex traits have revealed thousands of reproducible genetic associations to date, these loci collectively confer very little of the heritability of their respective diseases and, in general, have contributed little to our understanding the underlying disease biology. Physical protein interactions have been utilized to increase our understanding of human Mendelian disease loci but have yet to be fully exploited for complex traits. METHODS: We hypothesized that protein interaction modeling of GWAS findings could highlight important disease-associated loci and unveil the role of their network topology in the genetic architecture of diseases with complex inheritance. RESULTS: Network modeling of proteins associated with the intragenic single nucleotide polymorphisms of the National Human Genome Research Institute catalog of complex trait GWAS revealed that complex trait associated loci are more likely to be hub and bottleneck genes in available, albeit incomplete, networks (OR=1.59, Fisher's exact test p < 2.24 × 10(-12)). Network modeling also prioritized novel type 2 diabetes (T2D) genetic variations from the Finland-USA Investigation of Non-Insulin-Dependent Diabetes Mellitus Genetics and the Wellcome Trust GWAS data, and demonstrated the enrichment of hubs and bottlenecks in prioritized T2D GWAS genes. The potential biological relevance of the T2D hub and bottleneck genes was revealed by their increased number of first degree protein interactions with known T2D genes according to several independent sources (p<0.01, probability of being first interactors of known T2D genes). CONCLUSION: Virtually all common diseases are complex human traits, and thus the topological centrality in protein networks of complex trait genes has implications in genetics, personal genomics, and therapy.