Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses.

The role of the serine/glycine metabolic pathway (SGP) has recently been demonstrated in tumors; however, the pathological relevance of the SGP in thyroid cancer remains unexplored. Here, we perform metabolomic profiling of 17 tumor-normal pairs; bulk transcriptomics of 263 normal thyroid, 348 papillary, and 21 undifferentiated thyroid cancer samples; and single-cell transcriptomes from 15 cases, showing the impact of mitochondrial one-carbon metabolism in thyroid tumors. High expression of serine hydroxymethyltransferase-2 (SHMT2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is associated with low thyroid differentiation scores and poor clinical features. A subpopulation of tumor cells with high mitochondrial one-carbon pathway activity is observed in the single-cell dataset. SHMT2 inhibition significantly compromises mitochondrial respiration and decreases cell proliferation and tumor size in vitro and in vivo. Collectively, our results highlight the importance of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer and suggest that SHMT2 is a potent therapeutic target.

Single Cell Analysis of Human Thyroid Reveals the Transcriptional Signatures of Aging.

The thyroid gland plays a critical role in the maintenance of whole-body metabolism. However, aging frequently impairs homeostatic maintenance by thyroid hormones due to increased prevalence of subclinical hypothyroidism associated with mitochondrial dysfunction, inflammation, and fibrosis. To understand the specific aging-related changes of endocrine function in thyroid epithelial cells, we performed single-cell RNA sequencing (RNA-seq) of 54 726 cells derived from pathologically normal thyroid tissues from 7 patients who underwent thyroidectomy. Thyroid endocrine epithelial cells were clustered into 5 distinct subpopulations, and a subset of cells was found to be particularly vulnerable with aging, showing functional deterioration associated with the expression of metallothionein (MT) and major histocompatibility complex class II genes. We further validated that increased expression of MT family genes are highly correlated with thyroid gland aging in bulk RNAseq datasets. This study provides evidence that aging induces specific transcriptomic changes across multiple cell populations in the human thyroid gland.

NgR1 is an NK cell inhibitory receptor that destabilizes the immunological synapse.

The formation of an immunological synapse (IS) is essential for natural killer (NK) cells to eliminate target cells. Despite an advanced understanding of the characteristics of the IS and its formation processes, the mechanisms that regulate its stability via the cytoskeleton are unclear. Here, we show that Nogo receptor 1 (NgR1) has an important function in modulating NK cell-mediated killing by destabilization of IS formation. NgR1 deficiency or blockade resulted in improved tumor control of NK cells by enhancing NK-to-target cell contact stability and regulating F-actin dynamics during IS formation. Patients with tumors expressing abundant NgR1 ligand had poor prognosis despite high levels of NK cell infiltration. Thus, our study identifies NgR1 as an immune checkpoint in IS formation and indicates a potential approach to improve the cytolytic function of NK cells in cancer immunotherapy.

Stepwise molecular mechanisms responsible for chemoresistance in bladder cancer cells.

Chemotherapy resistance is an obstacle to cancer therapy and is considered a major cause of recurrence. Thus, understanding the mechanisms of chemoresistance is critical to improving the prognosis of patients. Here, we have established a stepwise gemcitabine-resistant T24 bladder cancer cell line to understand the molecular mechanisms of chemoresistance within cancer cells. The characteristics of the stepwise chemoresistance cell line were divided into 4 phases (parental, early, intermediate, and late phases). These four phase cells showed increasingly aggressive phenotypes in vitro and in vivo experiments with increasing phases and revealed the molecular properties of the biological process from parent cells to phased gemcitabine-resistant cell line (GRC). Taken together, through the analysis of gene expression profile data, we have characterized gene set of each phase indicating the response to anticancer drug treatment. Specifically, we identified a multigene signature (23 genes including GATA3, APOBEC3G, NT5E, MYC, STC1, FOXD1, SMAD9) and developed a chemoresistance score consisting of that could predict eventual responsiveness to gemcitabine treatment. Our data will contribute to predicting chemoresistance and improving the prognosis of bladder cancer patients.

YAP1 activation is associated with the progression and response to immunotherapy of non-muscle invasive bladder cancer.

BACKGROUND: Despite the availability of several treatments for non-muscle-invasive bladder cancer (NMIBC), many patients are still not responsive to treatments, and the disease progresses. A new prognostic classifier can differentiate between treatment response and progression, and it could be used as a very important tool in patient decision-making regarding treatment options. In this study, we focused on the activation of Yes-associated protein 1 (YAP1), which is known to play a pivotal role in tumour progression and serves as a factor contributing to the mechanism of resistance to various relevant therapeutic agents. We further evaluated its potential as a novel prognostic agent. METHODS: We identified YAP1-associated gene signatures based on UC3-siYAP1 cells (n=8) and NMIBC cohort (n=460). Cross-validation was performed using 5 independent bladder cancer patient cohorts (n=1006). We also experimentally validated the changes of gene expression levels representing each subgroup. FINDINGS: The 976-gene signature based on YAP1-activation redefined three subgroups and had the benefits of Bacillus Calmette-Guérin (BCG) treatment in patients with NMIBC (hazard ratio 3.32, 95% CI 1.29-8.56, p = 0.01). The integrated analysis revealed that YAP1 activation was associated with the characterization of patients with high-risk NMIBC and the response to immunotherapy. INTERPRETATION: This study suggests that YAP1 activation has an important prognostic effect on bladder cancer progression and might be useful in the selection of immunotherapy. FUNDING: A funding list that contributed to this research can be found in the Acknowledgements section.

Transcriptomic Analysis of Papillary Thyroid Cancer: A Focus on Immune-Subtyping, Oncogenic Fusion, and Recurrence.

OBJECTIVES: Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer. METHODS: This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database. RESULTS: Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension). CONCLUSION: We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.

E2F1 Promotes Progression of Bladder Cancer by Modulating RAD54L Involved in Homologous Recombination Repair.

DNA repair defects are important factors in cancer development. High DNA repair activity can affect cancer progression and chemoresistance. DNA double-strand breaks in cancer cells caused by anticancer agents can be restored by non-homologous end joining (NHEJ) and homologous recombination repair (HRR). Our previous study has identified E2F1 as a key gene in bladder cancer progression. In this study, DNA repair genes related to E2F1 were analyzed, and RAD54L involved in HRR was identified. In gene expression analysis of bladder cancer patients, the survival of patients with high RAD54L expression was shorter with cancer progression than in patients with low RAD54L expression. This study also revealed that E2F1 directly binds to the promoter region of RAD54L and regulates the transcription of RAD54L related to the HRR pathway. This study also confirmed that DNA breaks are repaired by RAD54L induced by E2F1 in bladder cancer cells treated with MMC. In summary, RAD54L was identified as a new target directly regulated by E2F1. Our results suggest that, E2F1 and RAD54L could be used as diagnostic markers for bladder cancer progression and represent potential therapeutic targets.

Fanconi Anemia Pathway Activation by FOXM1 Is Critical to Bladder Cancer Recurrence and Anticancer Drug Resistance.

Although the 5-year survival rate of patients diagnosed with nonmuscle invasive bladder cancer (NMIBC) has reached 85%, more than 50% of patients suffer from frequent recurrences. To identify molecular targets associated with recurrence of NMIBC, we analyzed gene expression data and found that FOXM1 and FANCD2 were involved in recurrence. Therefore, we investigated how these genes were involved in the mechanism of recurrence and confirmed their usefulness as biomarkers. Investigation have shown that FOXM1 directly regulated the transcription of FANCD2, which is the key gene of the Fanconi anemia (FA) pathway. Depletion of FOXM1 resulted in DNA repair defects in the FA pathway and in decreased resistance to chemotherapy. Thus, the FANCD2-associated FA pathway activated by FOXM1 is an important mechanism involved in chemotherapy-related recurrence. In conclusion, FOXM1 and FANCD2 can be used as prognostic factors that are associated with high risk of recurrence and with anticancer drug resistance properties in NMIBC patients.

Transcriptional Profiling of Advanced Urothelial Cancer Predicts Prognosis and Response to Immunotherapy.

Recent investigations reported that some subtypes from the Lund or The Cancer Genome Atlas (TCGA) classifications were most responsive to PD-L1 inhibitor treatment. However, the association between previously reported subtypes and immune checkpoint inhibitor (ICI) therapy responsiveness has been insufficiently explored. Despite these contributions, the ability to predict the clinical applicability of immune checkpoint inhibitor therapy in patients remains a major challenge. Here, we aimed to re-classify distinct subtypes focusing on ICI responsiveness using gene expression profiling in the IMvigor 210 cohort (n = 298). Based on the hierarchical clustering analysis, we divided advanced urothelial cancer patients into three subgroups. To confirm a prognostic impact, we performed survival analysis and estimated the prognostic value in the IMvigor 210 and TCGA cohort. The activation of CD8+ T effector cells was common for patients of classes 2 and 3 in the TCGA and IMvigor 210 cohort. Survival analysis showed that patients of class 3 in the TCGA cohort had a poor prognosis, while patients of class 3 showed considerably prolonged survival in the IMvigor 210 cohort. One of the distinct characteristics of patients in class 3 is the inactivation of the TGFß and YAP/TAZ pathways and activation of the cell cycle and DNA replication and DNA damage (DDR). Based on our identified transcriptional patterns and the clinical outcomes of advanced urothelial cancer patients, we constructed a schematic summary. When comparing clinical and transcriptome data, patients with downregulation of the TGFß and YAP/TAZ pathways and upregulation of the cell cycle and DDR may be more responsive to ICI therapy.

Genomic Signature of the Standardized Uptake Value in 18F-Fluorodeoxyglucose Positron Emission Tomography in Breast Cancer.

The standardized uptake value (SUV), an indicator of the degree of glucose uptake in 18F-fluorodeoxyglucose positron emission tomography (FDG-PET), has been used for predicting the clinical behavior of malignant tumors. However, its characteristics have been insufficiently explored at the genomics level. Here, we aim to identify genomic signatures reflecting prognostic SUV characteristics in breast cancer (BRC). Through integrative genomic profiling of 3710 BRC patients, including 254 patients who underwent preoperative FDG-PET, we identified an SUV signature, which showed independent clinical utility for predicting BRC prognosis (hazard ratio [HR] 1.27, 95% confidence interval [CI] = 1.12 to 1.45, p = 2.23 × 10-4). The risk subgroups classified by the signature exhibited mutually exclusive mutation patterns of TP53 and PIK3CA and showed significantly different responsiveness to immunotherapy. Experimental assays revealed that a signaling axis defined by TP53-FOXM1 and its downstream effectors in glycolysis-gluconeogenesis, including LDHA, might be important mediators in the FDG-PET process. Our molecular characterizations support an understanding of glucose metabolism and poor prognosis in BRC with a high SUV, utilizable in clinical practice to assist other diagnostic tools.

Identification of an immunotherapy-responsive molecular subtype of bladder cancer.

BACKGROUND: Although various molecular subtypes of bladder cancer (BC) have been investigated, most of these studies have focused on muscle-invasive BC (MIBC). A few studies have investigated non-muscle-invasive BC (NMIBC) or NMIBC and MIBC together, but none has classified progressive NMIBC or immune checkpoint inhibitor (ICI)-based therapeutic responses in early-stage BC patients. METHODS: A total of 1,934 samples from seven patient cohorts were used. We performed unsupervised hierarchical clustering to stratify patients into distinct subgroups and constructed a classifier by applying SAM/PAM algorithms. We then investigated the association between molecular subtypes and immunotherapy responsiveness using various statistical methods. FINDINGS: We explored large-scale genomic datasets encompassing NMIBC and MIBC, redefining four distinct molecular subtypes, including a subgroup containing progressive NMIBC and MIBC with poor prognosis that would benefit from ICI treatment. This subgroup showed poor progression-free survival with the distinct features of high mutation load, activated cell cycle, and inhibited TGFß signalling. Importantly, we verified that BC patients with this subtype were significantly responsive to an anti-PD-L1 agent in the IMvigor210 cohort. INTERPRETATION: Our results reveal an immunotherapeutic option for ICI treatment of highly progressive NMIBC and MIBC with poor prognosis. FUNDING: This research was supported by the National Research Foundation of Korea grant funded by the Korean government, a grant from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute, funded by the Ministry of Health and Welfare, Republic of Korea, and a grant from the KRIBB Research Initiative Program.

Clinical and genomic landscape of gastric cancer with a mesenchymal phenotype.

Gastric cancer is a heterogeneous cancer, making treatment responses difficult to predict. Here we show that we identify two distinct molecular subtypes, mesenchymal phenotype (MP) and epithelial phenotype (EP), by analyzing genomic and proteomic data. Molecularly, MP subtype tumors show high genomic integrity characterized by low mutation rates and microsatellite stability, whereas EP subtype tumors show low genomic integrity. Clinically, the MP subtype is associated with markedly poor survival and resistance to standard chemotherapy, whereas the EP subtype is associated with better survival rates and sensitivity to chemotherapy. Integrative analysis shows that signaling pathways driving epithelial-to-mesenchymal transition and insulin-like growth factor 1 (IGF1)/IGF1 receptor (IGF1R) pathway are highly activated in MP subtype tumors. Importantly, MP subtype cancer cells are more sensitive to inhibition of IGF1/IGF1R pathway than EP subtype. Detailed characterization of these two subtypes could identify novel therapeutic targets and useful biomarkers for prognosis and therapy response.

A gene expression signature of FOXM1 predicts the prognosis of hepatocellular carcinoma.

FOXM1 (Forkhead box M1) is a key regulator of tumorigenesis. Previous studies demonstrated that FOXM1 overexpression was strongly correlated with poor prognosis in various cancers, including hepatocellular carcinoma (HCC). In this study, we examined an association between the gene expression signature of FOXM1 and HCC patient outcome. The co-expressed gene set of FOXM1, which is significantly associated with the prognosis of HCC patients, was identified by analyzing the gene expression profiles of 100 patients with HCC, and this gene set was validated in two independent HCC patient cohorts (n=573). In multivariate analysis, the co-expressed gene set of FOXM1 was the most significant prognostic factor for overall survival in patients with HCC (hazard ratio=1.706, 95% confidence intervals=1.176-2.475, P=0.005). We also analyzed different types of cancer, including pancreatic adenocarcinoma, lung adenocarcinoma, breast carcinoma and bladder urothelial carcinoma, to verify the association between the co-expressed gene set of FOXM1 and patient prognosis, and we found a consistent prognostic significance, regardless of tumor type. Finally, we identified a putative signaling pathway in which miR-34a acts as an upstream regulator of the FOXM1-MYC signaling network; this pathway may be ultimately responsible for the poor prognosis of HCC patients. The prognostic subgroups defined by the gene expression signature of FOXM1 could help predict high-risk patients and may guide selection of the best treatment strategy.

Bioinformatic identification of prognostic signature defined by copy number alteration and expression of CCNE1 in non-muscle invasive bladder cancer.

Non-muscle invasive bladder cancer (NMIBC) patients frequently fail to respond to treatment and experience disease progression because of their clinical and biological diversity. In this study, we identify a prognostic molecular signature for predicting the heterogeneity of NMIBC by using an integrative analysis of copy number and gene expression data. We analyzed the copy number and gene expression profiles of 404 patients with bladder cancer obtained from The Cancer Genome Atlas (TCGA) consortium. Of the 14 molecules with significant copy number alterations that were previously reported, 13 were significantly correlated with copy number and expression changes. Prognostic gene sets based on the 13 genes were developed, and their prognostic values were verified in three independent patient cohorts (n=501). Among them, a signature of CCNE1 and its coexpressed genes was significantly associated with disease progression and validated in the independent cohorts. The CCNE1 signature was an independent risk factor based on the result of a multivariate analysis (hazard ratio=6.849, 95% confidence interval=1.613-29.092, P=0.009). Finally, gene network and upstream regulator analyses revealed that NMIBC progression is potentially mediated by CCND1-CCNE1-SP1 pathways. The prognostic molecular signature defined by copy number and expression changes of CCNE1 suggests a novel diagnostic tool for predicting the likelihood of NMIBC progression.

Heterogeneous Niche Activity of Ex-Vivo Expanded MSCs as Factor for Variable Outcomes in Hematopoietic Recovery.

Ex-vivo expanded mesenchymal stromal cells (MSCs) are increasingly used for paracrine support of hematopoietic stem cell (HSC) regeneration, but inconsistent outcomes have hindered ongoing clinical trials. Here, we show that significant heterogeneity in the niche activity of MSCs is created during their culture in various serum-supplemented media. The MSCs cultured under stimulatory or non-stimulatory culture conditions exhibited differences in colony forming unit-fibroblast contents, expression levels of cross-talk molecules (Jagged-1 and CXCL-12) and their support for HSC self-renewal. Accordingly, the enhancing effects of MSCs on hematopoietic engraftment were only visible when HSCs were co-transplanted with MSCs under stimulatory conditions. Of note, these differences in MSCs and their effects on HSCs were readily reversed by switching the cultures, indicating that the difference in niche activity can be caused by distinct functional state, rather than by clonal heterogeneity. Supporting the findings, transcriptomic analysis showed distinct upstream signaling pathways such as inhibition of P53 and activation of ER-stress response gene ATF4 for MSCs under stimulatory conditions. Taken together, our study shows that the niche activity of MSCs can vary rapidly by the extrinsic cues during culture causing variable outcomes in hematopoietic recoveries, and point to the possibility that MSCs can be pre-screened for more predictable efficacy in various cell therapy trials.

A polymorphic minisatellite region of BORIS regulates gene expression and its rare variants correlate with lung cancer susceptibility.

Aberrant expression of BORIS/CTCFL (Brother of the Regulator of Imprinted Sites/CTCF-like protein) is reported in different malignancies. In this study, we characterized the entire promoter region of BORIS/CTCFL, including the CpG islands, to assess the relationship between BORIS expression and lung cancer. To simplify the construction of luciferase reporter cassettes with various-sized portions of the upstream region, genomic copies of BORIS were isolated using TAR cloning technology. We analyzed three promoter blocks: the GATA/CCAAT box, the CpG islands and the minisatellite region BORIS-MS2. Polymorphic minisatellite sequences were isolated from genomic DNA prepared from the blood of controls and cases. Of the three promoter blocks, the GATA/CCAAT box was determined to be a critical element of the core promoter, while the CpG islands and the BORIS-MS2 minisatellite region were found to act as regulators. Interestingly, the polymorphic minisatellite region BORIS-MS2 was identified as a negative regulator that repressed the expression levels of luciferase reporter cassettes less effectively in cancer cells compared with normal cells. We also examined the association between the size of BORIS-MS2 and lung cancer in a case-control study with 590 controls and 206 lung cancer cases. Rare alleles of BORIS-MS2 were associated with a statistically significantly increased risk of lung cancer (odds ratio, 2.04; 95% confidence interval, 1.02-4.08; and P=0.039). To conclude, our data provide information on the organization of the BORIS promoter region and gene regulation in normal and cancer cells. In addition, we propose that specific alleles of the BORIS-MS2 region could be used to identify the risk for lung cancer.

Inactivation of Hippo Pathway Is Significantly Associated with Poor Prognosis in Hepatocellular Carcinoma.

PURPOSE: The Hippo pathway is a tumor suppressor in the liver. However, the clinical significance of Hippo pathway inactivation in HCC is not clearly defined. We analyzed genomic data from human and mouse tissues to determine clinical relevance of Hippo pathway inactivation in HCC. EXPERIMENTAL DESIGN: We analyzed gene expression data from Mst1/2(-/-) and Sav1(-/-) mice and identified a 610-gene expression signature reflecting Hippo pathway inactivation in the liver [silence of Hippo (SOH) signature]. By integrating gene expression data from mouse models with those from human HCC tissues, we developed a prediction model that could identify HCC patients with an inactivated Hippo pathway and used it to test its significance in HCC patients, via univariate and multivariate Cox analyses. RESULTS: HCC patients (National Cancer Institute cohort, n = 113) with the SOH signature had a significantly poorer prognosis than those without the SOH signature [P < 0.001 for overall survival (OS)]. The significant association of the signature with poor prognosis was further validated in the Korean (n = 100, P = 0.006 for OS) and Fudan University cohorts (n = 242, P = 0.001 for OS). On multivariate analysis, the signature was an independent predictor of recurrence-free survival (HR, 1.6; 95% confidence interval, 1.12-2.28: P = 0.008). We also demonstrated significant concordance between the SOH HCC subtype and the hepatic stem cell HCC subtype that had been identified in a previous study (P < 0.001). CONCLUSIONS: Inactivation of the Hippo pathway in HCC is significantly associated with poor prognosis.

A Database of Gene Expression Profiles of Korean Cancer Genome.

Because there are clear molecular differences entailing different treatment effectiveness between Korean and non-Korean cancer patients, identifying distinct molecular characteristics of Korean cancers is profoundly important. Here, we report a web-based data repository, namely Korean Cancer Genome Database (KCGD), for searching gene signatures associated with Korean cancer patients. Currently, a total of 1,403 cancer genomics data were collected, processed and stored in our repository, an ever-growing database. We incorporated most widely used statistical survival analysis methods including the Cox proportional hazard model, log-rank test and Kaplan-Meier plot to provide instant significance estimation for searched molecules. As an initial repository with the aim of Korean-specific marker detection, KCGD would be a promising web application for users without bioinformatics expertise to identify significant factors associated with cancer in Korean.

Activation of EZH2 and SUZ12 Regulated by E2F1 Predicts the Disease Progression and Aggressive Characteristics of Bladder Cancer.

PURPOSE: Previous study identified E2F1 as a key mediator of non-muscle-invasive bladder cancer (NMIBC) progression. The aim of this study was to identify the E2F1-related genes associated with poor prognosis and aggressive characteristics of bladder cancer. EXPERIMENTAL DESIGN: Microarray analysis was performed to find E2F1-related genes associated with tumor progression and aggressiveness in the gene expression data from 165 primary patients with bladder cancer. The biologic activity of E2F1-related genes in tumor progression and aggressiveness was confirmed with experimental assays using bladder cancer cells and tumor xenograft assay. RESULTS: The expression of E2F1 was significantly associated with EZH2 and SUZ12. The overexpression of E2F1, EZH2, and SUZ12 enhanced cancer progression including cell colony formation, migration, and invasiveness. Knockdown of these genes reduced motility, blocked invasion, and decreased tumor size in vivo. E2F1 bound the proximal EZH2 and SUZ12 promoter to activate transcription, suggesting that E2F1 and its downstream effectors, EZH2 and SUZ12, could be important mediators for the cancer progression. In addition, we confirmed an association between these genes and aggressive characteristics. Interestingly, the treatment of anticancer drugs to the cells overexpressing E2F1, EZH2, and SUZ12 induced the expression of CD44, KLF4, OCT4, and ABCG2 known as cancer stem cell (CSC)-related genes. CONCLUSIONS: The link between E2F1, EZH2, and/or SUZ12 revealed that E2f1 directly regulates transcription of the EZH2 and SUZ12 genes. The signature of E2F1-EZH2-SUZ12 shows a predictive value for prognosis in bladder tumors and the E2F1-EZH2-SUZ12-driven transcriptional events may regulate the cancer aggressiveness and chemo-resistance, which may provide opportunity for development of new treatment modalities.

Genomic predictors for recurrence patterns of hepatocellular carcinoma: model derivation and validation.

BACKGROUND: Typically observed at 2 y after surgical resection, late recurrence is a major challenge in the management of hepatocellular carcinoma (HCC). We aimed to develop a genomic predictor that can identify patients at high risk for late recurrence and assess its clinical implications. METHODS AND FINDINGS: Systematic analysis of gene expression data from human liver undergoing hepatic injury and regeneration revealed a 233-gene signature that was significantly associated with late recurrence of HCC. Using this signature, we developed a prognostic predictor that can identify patients at high risk of late recurrence, and tested and validated the robustness of the predictor in patients (n = 396) who underwent surgery between 1990 and 2011 at four centers (210 recurrences during a median of 3.7 y of follow-up). In multivariate analysis, this signature was the strongest risk factor for late recurrence (hazard ratio, 2.2; 95% confidence interval, 1.3-3.7; p = 0.002). In contrast, our previously developed tumor-derived 65-gene risk score was significantly associated with early recurrence (p = 0.005) but not with late recurrence (p = 0.7). In multivariate analysis, the 65-gene risk score was the strongest risk factor for very early recurrence (<1 y after surgical resection) (hazard ratio, 1.7; 95% confidence interval, 1.1-2.6; p = 0.01). The potential significance of STAT3 activation in late recurrence was predicted by gene network analysis and validated later. We also developed and validated 4- and 20-gene predictors from the full 233-gene predictor. The main limitation of the study is that most of the patients in our study were hepatitis B virus-positive. Further investigations are needed to test our prediction models in patients with different etiologies of HCC, such as hepatitis C virus. CONCLUSIONS: Two independently developed predictors reflected well the differences between early and late recurrence of HCC at the molecular level and provided new biomarkers for risk stratification. Please see later in the article for the Editors' Summary.