SUMMARY Liver cancer has the second highest worldwide cancer mortality rate and has limited therapeutic options. We analyzed 363 hepatocellular carcinoma (HCC) cases by whole exome sequencing and DNA copy number analyses, and 196 HCC also by DNA methylation, RNA, miRNA, and proteomic expression. DNA sequencing and mutation analysis identified significantly mutated genes including LZTR1, EEF1A1, SF3B1, and SMARCA4. Significant alterations by mutation or down-regulation by hypermethylation in genes likely to result in HCC metabolic reprogramming (ALB, APOB, and CPS1) were observed. Integrative molecular HCC subtyping incorporating unsupervised clustering of five data platforms identified three subtypes, one of which was associated with poorer prognosis in three HCC cohorts. Integrated analyses enabled development of a p53 target gene expression signature correlating with poor survival. Potential therapeutic targets for which inhibitors exist include WNT signaling, MDM4, MET, VEGFA, MCL1, IDH1, TERT, and immune checkpoint proteins CTLA-4, PD-1, and PD-L1.

SUMMARY To investigate the cell-intrinsic aging mechanisms that erode the function of somatic stem cells during aging, we have conducted a comprehensive integrated genomic analysis of young and aged cells. We profiled the transcriptome, DNA methylome, and histone modifications of young and old murine hematopoietic stem cells (HSCs). Transcriptome analysis indicated reduced TGFβ signaling and perturbation of genes involved in HSC proliferation and differentiation. Aged HSCs exhibited broader H3K4me3 peaks across HSC identity and self-renewal genes, and showed increased DNA methylation at transcription factor binding sites associated with differentiation-promoting genes combined with a reduction at genes associated with HSC maintenance. Together these changes reinforce HSC self-renewal and diminish differentiation, paralleling phenotypic HSC aging behavior. Ribosomal biogenesis emerged as a particular target of aging, with increased transcription of ribosomal protein and RNA genes, and hypomethylation of rRNA genes. This dataset will serve as a reference for future epigenomic analysis of stem cell aging.

Alternative polyadenylation (APA) is a pervasive mechanism in the regulation of most human genes, and its implication in diseases including cancer is only beginning to be appreciated. Since conventional APA profiling has not been widely adopted, global cancer APA studies are very limited. Here we develop a novel bioinformatics algorithm (DaPars) for the de novo identification of dynamic APAs from standard RNA-seq. When applied to 358 TCGA Pan-Cancer tumor/normal pairs across 7 tumor types, DaPars reveals 1,346 genes with recurrent and tumor-specific APAs. Most APA genes (91%) have shorter 3′ UTRs in tumors that can avoid miRNA-mediated repression, including glutaminase (GLS), a key metabolic enzyme for tumor proliferation. Interestingly, selected APA events add strong prognostic power beyond common clinical and molecular variables, suggesting their potential as novel prognostic biomarkers. Finally, our results implicate CstF64, an essential polyadenylation factor, as a master regulator of 3′ UTR shortening across multiple tumor types.

The global shortening of mRNAs through alternative polyadenylation (APA) that occurs during enhanced cellular proliferation represents an important, yet poorly understood mechanism of regulated gene expression1,2. The 3′UTR truncation of growth promoting mRNA transcripts that relieves intrinsic microRNA- and AU-rich element-mediated repression has been observed to correlate with cellular transformation3; however, the importance to tumorigenicity of RNA 3′ end processing factors that potentially govern APA is unknown. Here, we have identified CFIm25 as a broad repressor of proximal poly(A) site usage that, when depleted, increases cell proliferation. Applying a regression model on standard RNA-seq data for novel APA events, we identified at least 1,450 genes with shortened 3′UTRs after CFIm25 knockdown, representing 11% of significantly expressed mRNA in HeLa cells. Dramatic increases in expression of several known oncogenes including Cyclin D1 are observed as a consequence of CFIm25 depletion. Importantly, we identified a subset of CFIm25-regulated APA genes with shortened 3′UTRs in glioblastoma (GBM) tumors that have reduced CFIm25 expression. Downregulation of CFIm25 expression in glioblastoma cells enhances their tumorigenic properties and increases tumor size while CFIm25 overexpression reduces these properties and inhibits tumor growth. These findings identify a pivotal role of the CFIm25 in governing APA and reveal a previously unknown connection between CFIm25 and glioblastoma tumorigenicity.

Tumor suppressors are mostly defined by inactivating mutations in tumors, yet little is known about their epigenetic features in normal cells. Through integrative analysis of 1,134 genome-wide epigenetic profiles, mutations from >8,200 tumor-normal pairs, and our experimental data from clinical samples, we discovered broad H3K4me3 (wider than 4 kb) as the first epigenetic signature for tumor suppressors in normal cells. Broad H3K4me3 is associated with increased transcription elongation and enhancer activity together leading to exceptionally high gene expression, and is distinct from other broad epigenetic features, such as super-enhancers. Broad H3K4me3 conserved across normal cells may represent pan-cancer tumor suppressors, such as P53 and PTEN, whereas cell-type-specific broad H3K4me3 may indicate cell-identity genes and cell-type-specific tumor suppressors. Furthermore, widespread shortening of broad H3K4me3 in cancers is associated with repression of tumor suppressors. Together, the broad H3K4me3 epigenetic signature provides mutation-independent information for the discovery and characterization of novel tumor suppressors.