Atsushi Niida

Integrated molecular analysis of clear-cell renal cell carcinoma

Clear-cell renal cell carcinoma (ccRCC) is the most prevalent kidney cancer and its molecular pathogenesis is incompletely understood. Here we report an integrated molecular study of ccRCC in which ≥100 ccRCC cases were fully analyzed by whole-genome and/or whole-exome and RNA sequencing as well as by array-based gene expression, copy number and/or methylation analyses. We identified a full spectrum of genetic lesions and analyzed gene expression and DNA methylation signatures and determined their impact on tumor behavior. Defective VHL-mediated proteolysis was a common feature of ccRCC, which was caused not only by VHL inactivation but also by new hotspot TCEB1 mutations, which abolished Elongin C–VHL binding, leading to HIF accumulation. Other newly identified pathways and components recurrently mutated in ccRCC included PI3K-AKT-mTOR signaling, the KEAP1-NRF2-CUL3 apparatus, DNA methylation, p53-related pathways and mRNA processing. This integrated molecular analysis unmasked new correlations between DNA methylation, gene mutation and/or gene expression and copy number profiles, enabling the stratification of clinical risks for patients with ccRCC.

DKK1, a negative regulator of Wnt signaling, is a target of the β -catenin/TCF pathway

Wnt signaling plays an important role in embryonic development and tumorigenesis. These biological effects are exerted by activation of the β-catenin/TCF transcription complex and consequent regulation of a set of downstream genes. TCF-binding elements have been found in the promoter regions of many TCF target genes and characterized by a highly conserved consensus sequence. Utilizing this consensus sequence, we performed an in silico screening for new TCF target genes. Through computational screening and subsequent experimental analysis, we identified a novel TCF target gene, DKK1, which has been shown to be a potent inhibitor of Wnt signaling. Our finding suggests the existence of a novel feedback loop in Wnt signaling.

The Tumor-Suppressive miR-497-195 Cluster Targets Multiple Cell-Cycle Regulators in Hepatocellular Carcinoma

MicroRNAs (miRNAs) are key post-transcriptional regulators of gene expression and commonly deregulated in carcinogenesis. To explore functionally crucial tumor-suppressive (TS)-miRNAs in hepatocellular carcinoma (HCC), we performed integrative function- and expression-based screenings of TS-miRNAs in six HCC cell lines. The screenings identified seven miRNAs, which showed growth-suppressive activities through the overexpression of each miRNA and were endogenously downregulated in HCC cell lines. Further expression analyses using a large panel of HCC cell lines and primary tumors demonstrated four miRNAs, miR-101, -195, -378 and -497, as candidate TS-miRNAs frequently silenced in HCCs. Among them, two clustered miRNAs miR-195 and miR-497 showed significant growth-suppressive activity with induction of G1 arrest. Comprehensive exploration of their targets using Argonute2-immunoprecipitation-deep-sequencing (Ago2-IP-seq) and genome-wide expression profiling after their overexpression followed by pathway analysis, revealed a significant enrichment of cell cycle regulators. Among the candidates, we successfully identified CCNE1, CDC25A, CCND3, CDK4, and BTRC as direct targets for miR-497 and miR-195. Moreover, target genes frequently upregulated in HCC in a tumor-specific manner, such as CDK6, CCNE1, CDC25A and CDK4, showed an inverse correlation in the expression of miR-195 and miR-497, and their targets. These results suggest the molecular pathway regulating cell cycle progression to be integrally altered by downregulation of miR-195 and miR-497 expression, leading to the aberrant cell proliferation in hepatocarcinogenesis.

Human homolog of NOTUM, overexpressed in hepatocellular carcinoma, is regulated transcriptionally by β-catenin/TCF

The Drosophila Notum gene, which is regulated by the Wingless pathway, encodes a secreted hydrolase that modifies heparan sulfate proteoglycans. In comparative analysis of the gene expression profiles in primary human hepatocellular carcinomas (HCC) and normal organs, we observed that the human ortholog of Drosophila Notum was overexpressed markedly in a subset of HCC, but expressed rarely in adult normal tissues. Immunoblotting confirmed the overexpression of NOTUM protein in 12 of 40 primary HCC cases (30%). High levels of NOTUM protein were significantly associated with intracellular (nuclear or cytoplasmic) accumulation of β‐catenin protein: all 10 HCC with high intracellular β‐catenin also had high NOTUM expression, whereas only 2 of 30 cases (6.7%) without intracellular β‐catenin had high NOTUM expression (P < 0.00001). NOTUM expression in HepG2 cells was downregulated significantly by induction of a dominant‐negative mutant of TCF4, a β‐catenin partner. In vivo binding of the β‐catenin/TCF complex to the NOTUM promoter was demonstrated by chromatin immunoprecipitation in HepG2 and SW480 cells, where canonical Wnt signaling is activated constitutively. These findings provide evidence that NOTUM is a novel target of β‐catenin/TCF4 and is upregulated in Wnt/β‐catenin signaling‐activated HCC. (Cancer Sci 2008; 99: 1139–1146)

Integrated Multiregional Analysis Proposing a New Model of Colorectal Cancer Evolution.

Understanding intratumor heterogeneity is clinically important because it could cause therapeutic failure by fostering evolutionary adaptation. To this end, we profiled the genome and epigenome in multiple regions within each of nine colorectal tumors. Extensive intertumor heterogeneity is observed, from which we inferred the evolutionary history of the tumors. First, clonally shared alterations appeared, in which C>T transitions at CpG site and CpG island hypermethylation were relatively enriched. Correlation between mutation counts and patients’ ages suggests that the early-acquired alterations resulted from aging. In the late phase, a parental clone was branched into numerous subclones. Known driver alterations were observed frequently in the early-acquired alterations, but rarely in the late-acquired alterations. Consistently, our computational simulation of the branching evolution suggests that extensive intratumor heterogeneity could be generated by neutral evolution. Collectively, we propose a new model of colorectal cancer evolution, which is useful for understanding and confronting this heterogeneous disease.

A novel cell-cycle-indicator, mVenus-p27K-, identifies quiescent cells and visualizes G0-G1 transition

The quiescent (G0) phase of the cell cycle is the reversible phase from which the cells exit from the cell cycle. Due to the difficulty of defining the G0 phase, quiescent cells have not been well characterized. In this study, a fusion protein consisting of mVenus and a defective mutant of CDK inhibitor, p27 (p27K−) was shown to be able to identify and isolate a population of quiescent cells and to effectively visualize the G0 to G1 transition. By comparing the expression profiles of the G0 and G1 cells defined by mVenus-p27K−, we have identified molecular features of quiescent cells. Quiescence is also an important feature of many types of stem cells, and mVenus-p27K−-transgenic mice enabled the detection of the quiescent cells with muscle stem cell markers in muscle in vivo. The mVenus-p27K− probe could be useful in investigating stem cells as well as quiescent cells.

Gene set-based module discovery in the breast cancer transcriptome

BackgroundAlthough microarray-based studies have revealed global view of gene expression in cancer cells, we still have little knowledge about regulatory mechanisms underlying the transcriptome. Several computational methods applied to yeast data have recently succeeded in identifying expression modules, which is defined as co-expressed gene sets under common regulatory mechanisms. However, such module discovery methods are not applied cancer transcriptome data.ResultsIn order to decode oncogenic regulatory programs in cancer cells, we developed a novel module discovery method termed EEM by extending a previously reported module discovery method, and applied it to breast cancer expression data. Starting from seed gene sets prepared based on cis-regulatory elements, ChIP-chip data, and gene locus information, EEM identified 10 principal expression modules in breast cancer based on their expression coherence. Moreover, EEM depicted their activity profiles, which predict regulatory programs in each subtypes of breast tumors. For example, our analysis revealed that the expression module regulated by the Polycomb repressive complex 2 (PRC2) is downregulated in triple negative breast cancers, suggesting similarity of transcriptional programs between stem cells and aggressive breast cancer cells. We also found that the activity of the PRC2 expression module is negatively correlated to the expression of EZH2, a component of PRC2 which belongs to the E2F expression module. E2F-driven EZH2 overexpression may be responsible for the repression of the PRC2 expression modules in triple negative tumors. Furthermore, our network analysis predicts regulatory circuits in breast cancer cells.ConclusionThese results demonstrate that the gene set-based module discovery approach is a powerful tool to decode regulatory programs in cancer cells.

Integrative bioinformatics analysis of transcriptional regulatory programs in breast cancer cells.

BackgroundMicroarray technology has unveiled transcriptomic differences among tumors of various phenotypes, and, especially, brought great progress in molecular understanding of phenotypic diversity of breast tumors. However, compared with the massive knowledge about the transcriptome, we have surprisingly little knowledge about regulatory mechanisms underling transcriptomic diversity.ResultsTo gain insights into the transcriptional programs that drive tumor progression, we integrated regulatory sequence data and expression profiles of breast cancer into a Bayesian Network, and searched for cis-regulatory motifs statistically associated with given histological grades and prognosis. Our analysis found that motifs bound by ELK1, E2F, NRF1 and NFY are potential regulatory motifs that positively correlate with malignant progression of breast cancer.ConclusionThe results suggest that these 4 motifs are principal regulatory motifs driving malignant progression of breast cancer. Our method offers a more concise description about transcriptome diversity among breast tumors with different clinical phenotypes.

A Novel Network Profiling Analysis Reveals System Changes in Epithelial-Mesenchymal Transition

Patient-specific analysis of molecular networks is a promising strategy for making individual risk predictions and treatment decisions in cancer therapy. Although systems biology allows the gene network of a cell to be reconstructed from clinical gene expression data, traditional methods, such as Bayesian networks, only provide an averaged network for all samples. Therefore, these methods cannot reveal patient-specific differences in molecular networks during cancer progression. In this study, we developed a novel statistical method called NetworkProfiler, which infers patient-specific gene regulatory networks for a specific clinical characteristic, such as cancer progression, from gene expression data of cancer patients. We applied NetworkProfiler to microarray gene expression data from 762 cancer cell lines and extracted the system changes that were related to the epithelial-mesenchymal transition (EMT). Out of 1732 possible regulators of E-cadherin, a cell adhesion molecule that modulates the EMT, NetworkProfiler, identified 25 candidate regulators, of which about half have been experimentally verified in the literature. In addition, we used NetworkProfiler to predict EMT-dependent master regulators that enhanced cell adhesion, migration, invasion, and metastasis. In order to further evaluate the performance of NetworkProfiler, we selected Krueppel-like factor 5 (KLF5) from a list of the remaining candidate regulators of E-cadherin and conducted in vitro validation experiments. As a result, we found that knockdown of KLF5 by siRNA significantly decreased E-cadherin expression and induced morphological changes characteristic of EMT. In addition, in vitro experiments of a novel candidate EMT-related microRNA, miR-100, confirmed the involvement of miR-100 in several EMT-related aspects, which was consistent with the predictions obtained by NetworkProfiler.

The AURKA/TPX2 axis drives colon tumorigenesis cooperatively with MYC

BACKGROUND The MYC oncogene has long been established as a central driver in many types of human cancers including colorectal cancer. However, the realization of MYC-targeting therapies remains elusive; as a result, synthetic lethal therapeutic approaches are alternatively being explored. A synthetic lethal therapeutic approach aims to kill MYC-driven tumors by targeting a certain co-regulator on the MYC pathway. PATIENTS AND METHODS We analyzed copy number and expression profiles from 130 colorectal cancer tumors together with publicly available datasets to identify co-regulators on the MYC pathway. Candidates were functionally tested by in vitro assays using colorectal cancer and normal fibroblast cell lines. Additionally, survival analyses were carried out on another 159 colorectal cancer patients and public datasets. RESULTS Our in silico screening identified two MYC co-regulator candidates, AURKA and TPX2, which are interacting mitotic regulators located on chromosome 20q. We found the two candidates showed frequent co-amplification with the MYC locus while expression levels of MYC and the two genes were positively correlated with those of MYC downstream target genes across multiple cancer types. In vitro, the aberrant expression of MYC, AURKA and TPX2 resulted in more aggressive anchorage-independent growth in normal fibroblast cells. Furthermore, knockdown of AURKA or TPX2, or treatment with an AURKA-specific inhibitor effectively suppressed the proliferation of MYC-expressing colorectal cancer cells. Additionally, combined high expression of MYC, AURKA and TPX2 proved to be a poor prognostic indicator of colorectal cancer patient survival. CONCLUSIONS Through bioinformatic analyses and experiments, we proposed TPX2 and AURKA as novel co-regulators on the MYC pathway. Inhibiting the AURKA/TPX2 axis would be a novel synthetic lethal therapeutic approach for MYC-driven cancers.