Recurrence-associated gene signature in patients with stage I non-small-cell lung cancer.

Recurrent gene mutations and fusions in cancer patients are likely to be associated with cancer progression or recurrence by Vogelstein et al. (Science (80-) 340, 1546-1558 (2013)). In this study, we investigated gene mutations and fusions that recurrently occurred in early-stage cancer patients with stage I non-small-cell cancer (NSCLC). Targeted exome sequencing was performed to profile the variants and confirmed their fidelity at the gene and pathway levels through comparison with data for stage I lung cancer patients, which was obtained from The Cancer Genome Atlas (TCGA). Next, we identified prognostic gene mutations (ATR, ERBB3, KDR, and MUC6), fusions (GOPC-ROS1 and NTRK1-SH2D2A), and VEGF signaling pathway associated with cancer recurrence. To infer the functional implication of the recurrent variants in early-stage cancers, the extent of their selection pattern was investigated, and they were shown to be under positive selection, implying a selective advantage for cancer progression. Specifically, high selection scores were observed in the variants with significantly high risks for recurrence. Taken together, the results of this study enabled us to identify recurrent gene mutations and fusions in a stage I NSCLC cohort and to demonstrate positive selection, which had implications regarding cancer recurrence.

Longitudinal change of genetic variations in cetuximab-treated metastatic colorectal cancer.

Recurrent gene mutations and copy number alterations in cancer patients are presumably associated with resistance to targeted therapy. In the present study, we assessed the gene mutations and copy number alterations that recurrently occurred in cetuximab-treated patients with metastatic colorectal cancer (mCRC). Targeted next-generation sequencing was performed in the tumor samples obtained pre- and postcetuximab treatment to assess the variations that occurred during cetuximab treatment. Moreover, we identified the emergent gene mutations (CDK6, EPHA3, ERCC2, MYC, PCMTD1, PIK3CA, PRIM2, RICTOR, and ZNRF3) and copy number alterations (ARAF, BCL2, BRCA2, EGFR, MYC, and SMAD4) that were recurrently observed only in postprogression samples and not in pretreatment or posttreatment samples from patients revealing clinical response. Furthermore, to identify the feasible candidate variations implicated in treatment resistance, we examined the variants with clonal expansion during treatment and discovered PCBP1 as a variant associated with posttreatment progression. Various recurrent mutations were enriched in the TGF-beta signaling pathway. Collectively, we identified recurrent variations in mCRC samples exhibiting post-cetuximab progression. Additionally, future studies are required to evaluate the therapeutic potential of these variations.

Predicting clinical benefit of immunotherapy by antigenic or functional mutations affecting tumour immunogenicity.

Neoantigen burden is regarded as a fundamental determinant of response to immunotherapy. However, its predictive value remains in question because some tumours with high neoantigen load show resistance. Here, we investigate our patient cohort together with a public cohort by our algorithms for the modelling of peptide-MHC binding and inter-cohort genomic prediction of therapeutic resistance. We first attempt to predict MHC-binding peptides at high accuracy with convolutional neural networks. Our prediction outperforms previous methods in > 70% of test cases. We then develop a classifier that can predict resistance from functional mutations. The predictive genes are involved in immune response and EGFR signalling, whereas their mutation patterns reflect positive selection. When integrated with our neoantigen profiling, these anti-immunogenic mutations reveal higher predictive power than known resistance factors. Our results suggest that the clinical benefit of immunotherapy can be determined by neoantigens that induce immunity and functional mutations that facilitate immune evasion.

Network perturbation by recurrent regulatory variants in cancer.

Cancer driving genes have been identified as recurrently affected by variants that alter protein-coding sequences. However, a majority of cancer variants arise in noncoding regions, and some of them are thought to play a critical role through transcriptional perturbation. Here we identified putative transcriptional driver genes based on combinatorial variant recurrence in cis-regulatory regions. The identified genes showed high connectivity in the cancer type-specific transcription regulatory network, with high outdegree and many downstream genes, highlighting their causative role during tumorigenesis. In the protein interactome, the identified transcriptional drivers were not as highly connected as coding driver genes but appeared to form a network module centered on the coding drivers. The coding and regulatory variants associated via these interactions between the coding and transcriptional drivers showed exclusive and complementary occurrence patterns across tumor samples. Transcriptional cancer drivers may act through an extensive perturbation of the regulatory network and by altering protein network modules through interactions with coding driver genes.

Chromatin structure-based prediction of recurrent noncoding mutations in cancer.

Recurrence is a hallmark of cancer-driving mutations. Recurrent mutations can arise at the same site or affect the same gene at different sites. Here we identified a set of mutations arising in individual samples and altering different cis-regulatory elements that converge on a common gene via chromatin interactions. The mutations and genes identified in this fashion showed strong relevance to cancer, in contrast to noncoding mutations with site-specific recurrence only. We developed a prediction method that identifies potentially recurrent mutations on the basis of the features shared by mutations whose recurrence is observed in a given cohort. Our method was capable of accurately predicting recurrent mutations at the level of target genes but not mutations recurring at the same site. We experimentally validated predicted mutations in distal regulatory regions of the TERT gene. In conclusion, we propose a novel approach to discovering potential cancer-driving mutations in noncoding regions.

Global transcription network incorporating distal regulator binding reveals selective cooperation of cancer drivers and risk genes.

Global network modeling of distal regulatory interactions is essential in understanding the overall architecture of gene expression programs. Here, we developed a Bayesian probabilistic model and computational method for global causal network construction with breast cancer as a model. Whereas physical regulator binding was well supported by gene expression causality in general, distal elements in intragenic regions or loci distant from the target gene exhibited particularly strong functional effects. Modeling the action of long-range enhancers was critical in recovering true biological interactions with increased coverage and specificity overall and unraveling regulatory complexity underlying tumor subclasses and drug responses in particular. Transcriptional cancer drivers and risk genes were discovered based on the network analysis of somatic and genetic cancer-related DNA variants. Notably, we observed that the risk genes were functionally downstream of the cancer drivers and were selectively susceptible to network perturbation by tumorigenic changes in their upstream drivers. Furthermore, cancer risk alleles tended to increase the susceptibility of the transcription of their associated genes. These findings suggest that transcriptional cancer drivers selectively induce a combinatorial misregulation of downstream risk genes, and that genetic risk factors, mostly residing in distal regulatory regions, increase transcriptional susceptibility to upstream cancer-driving somatic changes.

Collagen triple helix repeat containing 1 (CTHRC1) acts via ERK-dependent induction of MMP9 to promote invasion of colorectal cancer cells.

Collagen triple helix repeat-containing 1 (CTHRC1) is known to be aberrantly upregulated in most human solid tumors, although the functional roles of CTHRC1 in colorectal cancer remain unclear. In this study, we investigated the occurrence of CTHRC1 upregulation and its role in vivo and in vitro. The expression profile and clinical importance of CTHRC1 were examined by reverse transcription-polymerase chain reaction and immunohistochemical analyses in normal and tumor patient samples. CTHRC1 was detectable in normal tissues, but also was highly expressed in tumor specimens. CTHRC1 upregulation was significantly associated with demethylation of the CTHRC1 promoter in colon cancer cell lines and tumor tissues. Clinicopathologic analyses showed that nodal status and expression of CTHRC1 (95% CI 0.999-3.984, p=0.05) were significant prognostic factors for disease-free survival. Promoter CpG methylation and hypermethylation status were measured by bisulfite sequencing and pyrosequencing analysis. Furthermore, we showed that overexpression of CTHRC1 in the SW480 and HT-29 cell lines increased invasiveness, an effect mediated by extracellular signal-regulated kinase (ERK)-dependent upregulation of matrix metalloproteinase 9 (MMP9). Consistent with this, we found that knockdown of CTHRC1 attenuated ERK activation and cancer cell invasivity. These results demonstrate that CTHRC1 expression is elevated in human colon cancer cell lines and clinical specimens, and promotes cancer cell invasivity through ERK-dependent induction of MMP9 expression. Our results further suggest that high levels of CTHRC1 expression are associated with poor clinical outcomes.

Genome-wide reorganization of histone H2AX toward particular fragile sites on cell activation.

γH2AX formation by phosphorylation of the histone variant H2AX is the key process in the repair of DNA lesions including those arising at fragile sites under replication stress. Here we demonstrate that H2AX is dynamically reorganized to preoccupy γH2AX hotspots on increased replication stress by activated cell proliferation and that H2AX is enriched in aphidicolin-induced replisome stalling sites in cycling cells. Interestingly, H2AX enrichment was particularly found in genomic regions that replicate in early S phase. High transcription activity, a hallmark of early replicating fragile sites, was a determinant of H2AX localization. Subtelomeric H2AX enrichment was also attributable to early replication and high gene density. In contrast, late replicating and infrequently transcribed regions, including common fragile sites and heterochromatin, lacked H2AX enrichment. In particular, heterochromatin was inaccessible to H2AX incorporation, maybe partly explaining the cause of mutation accumulation in cancer heterochromatin. Meanwhile, H2AX in actively dividing cells was intimately colocalized with INO80. INO80 silencing reduced H2AX levels, particularly at the INO80-enriched sites. Our findings suggest that active DNA replication is accompanied with the specific localization of H2AX and INO80 for efficient damage repair or replication-fork stabilization in actively transcribed regions.

Epigenetic regulation of microRNA-10b and targeting of oncogenic MAPRE1 in gastric cancer.

MicroRNAs act as negative regulators of gene expression, and the altered expression of microRNAs by epigenetic mechanisms is strongly implicated in carcinogenesis. Here we report that the microRNA-10b gene (miR-10b) was silenced in gastric cancer cells by promoter methylation. In this study, using a methylation array and bisulfate pyrosequencing analysis, we found that miR-10b promoter CpGs were heavily methylated in gastric cancers. Clinicopathologic data showed that miR-10b methylation increased with patient age and occurred significantly more frequently in intestinal-type (28/44, 64%) than in diffuse-type (22/56, 39%) gastric cancers (P = 0.016). In addition, miR-10b methylation was also associated with an increase in expression of the oncogene that encodes microtubule-associated protein, RP/EB family, member 1 (MAPRE1; P = 0.004), which was identified as a potential miR-10b target. After 5-aza-2'-deoxycytidine treatment of gastric cancer cells, miR-10b methylation was significantly decreased, and expression of miR-10b and HOXD4, which is 1 kb downstream of miR-10b, was greatly restored. Moreover, decreased MAPRE1 expression coincided with increased miR-10b expression, suggesting that miR-10b targets MAPRE1 transcription. We also found that transfection with precursor miR-10b into gastric cancer cells dramatically decreased MAPRE1 mRNA and protein, resulting in a significant decrease in colony formation and cell growth rates. Thus, we show a tumor-suppressive role for miR-10b in gastric carcinogenesis. miR-10b methylation may be a useful molecular biomarker for assessing the risk of gastric cancer development, and modulation of miR-10b may represent a therapeutic approach for treating gastric cancer.