Open reading frame dominance indicates protein-coding potential of RNAs.

Recent studies have identified numerous RNAs with both coding and noncoding functions. However, the sequence characteristics that determine this bifunctionality remain largely unknown. In the present study, we develop and test the open reading frame (ORF) dominance score, which we define as the fraction of the longest ORF in the sum of all putative ORF lengths. This score correlates with translation efficiency in coding transcripts and with translation of noncoding RNAs. In bacteria and archaea, coding and noncoding transcripts have narrow distributions of high and low ORF dominance, respectively, whereas those of eukaryotes show relatively broader ORF dominance distributions, with considerable overlap between coding and noncoding transcripts. The extent of overlap positively and negatively correlates with the mutation rate of genomes and the effective population size of species, respectively. Tissue-specific transcripts show higher ORF dominance than ubiquitously expressed transcripts, and the majority of tissue-specific transcripts are expressed in mature testes. These data suggest that the decrease in population size and the emergence of testes in eukaryotic organisms allowed for the evolution of potentially bifunctional RNAs.

Genomic landscape of chemical-induced lung tumors under Nrf2 different expression levels.

The transcription factor Nrf2 plays a crucial role in the anti-oxidative stress response, protection of DNA from injury and DNA repair mechanisms. Nrf2 activity reduces cancer initiation, but how Nrf2 affects whole-genome alterations upon carcinogenic stimulus remains unexplored. Although recent genome-wide analysis using next-generation sequencing revealed landscapes of nucleotide mutations and copy number alterations in various human cancers, genomic changes in murine cancer models have not been thoroughly examined. We elucidated the relationship between Nrf2 expression levels and whole exon mutation patterns using an ethyl-carbamate (urethane)-induced lung carcinogenesis model employing Nrf2-deficient and Keap1-kd mice, the latter of which express high levels of Nrf2. Exome analysis demonstrated that single nucleotide and trinucleotide mutation patterns and the Kras mutational signature differed significantly and were dependent on the expression level of Nrf2. The Nrf2-deficient tumors exhibited fewer copy number alterations relative to the Nrf2-wt and Keap1-kd tumors. The observed trend in genomic alterations likely prevented the Nrf2-deficient tumors from progressing into malignancy. For the first time, we present whole-exome sequencing results for chemically-induced lung tumors in the Nrf2 gain or loss of function mouse models. Our results demonstrate that different Nrf2 expression levels lead to distinct gene mutation patterns that underly different oncogenic mechanisms in each tumor genotype.

Prognostic biomarker study in patients with clinical stage I esophageal squamous cell carcinoma: JCOG0502-A1.

We undertook genomic analyses of Japanese patients with stage I esophageal squamous cell carcinoma (ESCC) to investigate the frequency of genomic alterations and the association with survival outcomes. Biomarker analysis was carried out for patients with clinical stage T1bN0M0 ESCC enrolled in JCOG0502 (UMIN000000551). Whole-exome sequencing (WES) was performed using DNA extracted from formalin-fixed, paraffin-embedded tissue of ESCC and normal tissue or blood sample. Single nucleotide variants (SNVs), insertions/deletions (indels), and copy number alterations (CNAs) were identified. We then evaluated the associations between each gene alteration with a frequency of 10% or more and progression-free survival (PFS) using a Cox regression model. We controlled for family-wise errors at 0.05 using the Bonferroni method. Among the 379 patients who were enrolled in JCOG0502, 127 patients were successfully analyzed using WES. The median patient age was 63 years (interquartile range, 57-67 years), and 78.0% of the patients ultimately underwent surgery. The 3-year PFS probability was 76.3%. We detected 20 genes with SNVs, indels, or amplifications with a frequency of 10% or more. Genomic alterations in FGF19 showed the strongest association with PFS with a borderline level of statistical significance of P = .00252 (Bonferroni-adjusted significance level is .0025). Genomic alterations in FGF4, MYEOV, CTTN, and ORAOV1 showed a marginal association with PFS (P < .05). These genomic alterations were all CNAs at chromosome 11q13.3. We have identified new genomic alterations associated with the poor efficacy of ESCC (T1bN0M0). These findings open avenues for the development of new potential treatments for patients with ESCC.

Single-cell DNA and RNA sequencing reveals the dynamics of intra-tumor heterogeneity in a colorectal cancer model.

BACKGROUND: Intra-tumor heterogeneity (ITH) encompasses cellular differences in tumors and is related to clinical outcomes such as drug resistance. However, little is known about the dynamics of ITH, owing to the lack of time-series analysis at the single-cell level. Mouse models that recapitulate cancer development are useful for controlled serial time sampling. RESULTS: We performed single-cell exome and transcriptome sequencing of 200 cells to investigate how ITH is generated in a mouse colorectal cancer model. In the model, a single normal intestinal cell is grown into organoids that mimic the intestinal crypt structure. Upon RNAi-mediated downregulation of a tumor suppressor gene APC, the transduced organoids were serially transplanted into mice to allow exposure to in vivo microenvironments, which play relevant roles in cancer development. The ITH of the transcriptome increased after the transplantation, while that of the exome decreased. Mutations generated during organoid culture did not greatly change at the bulk-cell level upon the transplantation. The RNA ITH increase was due to the emergence of new transcriptional subpopulations. In contrast to the initial cells expressing mesenchymal-marker genes, new subpopulations repressed these genes after the transplantation. Analyses of colorectal cancer data from The Cancer Genome Atlas revealed a high proportion of metastatic cases in human subjects with expression patterns similar to the new cell subpopulations in mouse. These results suggest that the birth of transcriptional subpopulations may be a key for adaptation to drastic micro-environmental changes when cancer cells have sufficient genetic alterations at later tumor stages. CONCLUSIONS: This study revealed an evolutionary dynamics of single-cell RNA and DNA heterogeneity in tumor progression, giving insights into the mesenchymal-epithelial transformation of tumor cells at metastasis in colorectal cancer.

tugHall: a simulator of cancer-cell evolution based on the hallmarks of cancer and tumor-related genes.

SUMMARY: The flood of recent cancer genomic data requires a coherent model that can sort out the findings to systematically explain clonal evolution and the resultant intra-tumor heterogeneity (ITH). Here, we present a new mathematical model designed to computationally simulate the evolution of cancer cells. The model connects the well-known hallmarks of cancer with the specific mutational states of tumor-related genes. The cell behavior phenotypes are stochastically determined, and the hallmarks probabilistically interfere with the phenotypic probabilities. In turn, the hallmark variables depend on the mutational states of tumor-related genes. Thus, our software can deepen our understanding of cancer-cell evolution and generation of ITH. AVAILABILITY AND IMPLEMENTATION: The open-source code is available in the repository https://github.com/nagornovys/Cancer_cell_evolution. CONTACT: mamkato@ncc.go.jp. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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

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

Comprehensive analysis of DNA adducts (DNA adductome analysis) in the liver of rats treated with 1,4-dioxane.

1,4-Dioxane is a genotoxic carcinogen, and its mutagenic properties were recently observed in the liver of guanine phosphoribosyl transferase (gpt) delta transgenic rats. However, the mechanisms of its genotoxicity remain unclear. We analyzed DNA adduct formation in rat livers following 1,4-dioxane treatment. After administering 1,4-dioxane in drinking water at doses of 0, 20, 200, and 5,000 ppm, liver adduct formation was analyzed by DNA adductome analysis. Adducts in treated rat livers were dose-dependently increased compared with those in the control group. Principal component analysis-discriminant analysis (PCA-DA) clearly revealed two clusters of DNA adducts, associated with 0 ppm and low-dose (20 ppm) 1,4-dioxane-treatment versus middle- and high-dose (200, 5,000 ppm)-treated rats. After confirming the intensity of each adduct, three adducts were screened as characteristic of 1,4-dioxane treatment. Two of the three candidates contained thymine or cytidine/uracil moieties. Another candidate was identified as 8-oxo-dG based on mass fragmentation together with high-resolution accurate-mass (HRAM) mass spectrometry data. Oxidative stress responses may partly explain the mechanisms of increased mutations in the liver of gpt delta rats following 1,4-dioxane treatment.

Frequent mutations of genes encoding vacuolar H+ -ATPase components in granular cell tumors.

Granular cell tumors (GCTs) are rare mesenchymal tumors that exhibit a characteristic morphology and a finely granular cytoplasm. The genetic alterations responsible for GCT tumorigenesis had been unknown until recently, when loss-of-function mutations of ATP6AP1 and ATP6AP2 were described. Thus, we performed whole-exome sequencing, RNA sequencing, and targeted sequencing of 51 GCT samples. From these genomic analyses, we identified mutations in genes encoding vacuolar H+ -ATPase (V-ATPase) components, including ATP6AP1 and ATP6AP2, in 33 (65%) GCTs. ATP6AP1 and ATP6AP2 mutations were found in 23 (45%) and 2 (4%) samples, respectively, and all were truncating or splice site mutations. In addition, seven other genes encoding V-ATPase components were also mutated, and three mutations in ATP6V0C occurred on the same amino acid (isoleucine 136). These V-ATPase component gene mutations were mutually exclusive, with one exception. These results suggest that V-ATPase function is impaired in GCTs not only by loss-of-function mutations of ATP6AP1 and ATP6AP2 but also through mutations of other subunits. Our findings provide additional support for the hypothesis that V-ATPase dysfunction promotes GCT tumorigenesis.

[Examination of Effectiveness Verification and Additional Items in Angiography and IVR of DRLs 2015].

Japanese Diagnostic Reference Levels (DRLs) were released as "Japan DRLs 2015" from Japan Network for Research and Information on Medical Exposure (J-RIME) in June 2015. In "Japan DRLs 2015", DRLs in angiography and interventional procedures are set at a fluoroscopic dose rate of 20 mGy/min at the interventional reference point using a phantom. In order to achieve optimization with DRLs, then it need to be revised regularly. Therefore, we (research group to examine the effect of Japan DRLs 2015 and the necessity of additional items in angiography and vascular interventions) examined the effects of "Japan DRLs 2015" on angiography and interventional procedures. And we also investigated for DRLs revision in the future. As a result, it turned out that it is important to create DRLs in medical procedures that can be effectively used in clinical settings.

[Evaluating Eye Lens Dose of Neurovascular and Cardiac Electrophysiology Interventional Physician].

PURPOSE: The International Commission on Radiological Protection recommended that interventional radiologies (IRs) have high radiation doses and that staff may also be exposed to high doses. In the present study, we measured the radiation exposure dose [3 mm dose equivalent, Hp (3) ] in the eye using an appropriate dosimeter placed next to the physician' s eye during neurovascular intervention procedure (Neuro-IR) and interventional cardiac electrophysiology procedure (EP-IR). METHOD: Physicians wore a direct eye dosemeter just lateral to the left eye and an additional direct eye dosemeter outside the radiation protective glasses close to their left eye. Additionally, a neck badge [0.07 mm dose equivalent, Hp (0.07) ] was worn outside the protective apron to the left of the neck, to compare the direct eye dosimeter estimated doses. The occupational eye lens dose was evaluated over a period of 6-month. RESULTS: The maximum Hp (3) of the Neuro-IR physician was estimated 5.1 mSv without the radiation protective glasses and 1.6 mSv with the radiation protective glasses. On the other hand, the maximum Hp (3) of the EP-IR physician was estimated 29 mSv without the radiation protective glasses and 15 mSv with the radiation protective glasses. CONCLUSION: Physicians eye lens dose [Hp (3) ] tended to be overestimated by the neck badge measurements [Hp (0.07)]. A correct evaluation of the lens dose [Hp (3) ] using the direct eye dosimeter is recommended. Although we found a positive correlation between Hp (0.07) and Hp (3), the value of R2 in the regression equation is low, we recommended that the eye lens dose estimated carefully from Hp (0.07).

[Nationwide Survey of Medical Radiation Exposure on Cardiovascular Examinations in Japan].

PURPOSE: It is very important to manage the radiation dose of cardiovascular interventional (CVI) procedures. Overseas, the diagnostic reference levels for cardiac interventional procedures were established with the air kerma at the patient entrance reference point (Ka,r) and the air kerma-area product (PKA). Although the Japan DRLs 2015 was established by the Japan Network for Research and Information on Medical Exposure (J-RIME), the Japan DRL for CVIs were established by fluoroscopic dose rates of 20 mGy/min at the patient entrance reference point with 20 cm thickness polymethyl methacrylate (PMMA) phantom. In the present our study, we performed a questionnaire survey of indicated values of angiographic parameters in CVI procedures. METHODS: A nationwide questionnaire was sent by post to 765 facilities. Question focused on angiographic technology, exposure parameters and radiation doses as the displayed dosimetric parameters on the angiographic machine. RESULTS: The recovery rate was 22.8% at 175 out of 765 facilities. In total 1728 cases of the coronary angiography (CAG), 1703 cases of the percutaneous coronary intervention (PCI), 962 cases of the radiofrequency catheter ablation (RFCA) and 377 cases of pediatric CVI. The 75th percentile value of Ka,r, PKA, fluoroscopy time (FT) and number of cine images (CI) for CAG, PCI, RFCA and pediatric CVI were 702, 2042, 644, and 159 mGy, respectively, 59.3, 152, 81.3, and 14.9 Gy・cm2, respectively, 10.2, 35.6, 61.1, and 35.6 min, respectively and 1503, 2672, 722, and 2378 images, respectively. Our investigation showed that the angiographic parameters were different in several CVI procedures. CONCLUSIONS: The displayed dosimetric parameters on the angiographic machine in CVI procedures showed different values. We should classify the dosimetric parameters for each procedure.

Feasibility and utility of a panel testing for 114 cancer-associated genes in a clinical setting: A hospital-based study.

Next-generation sequencing (NGS) of tumor tissue (ie, clinical sequencing) can guide clinical management by providing information about actionable gene aberrations that have diagnostic and therapeutic significance. Here, we undertook a hospital-based prospective study (TOP-GEAR project, 2nd stage) to investigate the feasibility and utility of NGS-based analysis of 114 cancer-associated genes (the NCC Oncopanel test). We examined 230 cases (comprising more than 30 tumor types) of advanced solid tumors, all of which were matched with nontumor samples. Gene profiling data were obtained for 187 cases (81.3%), 111 (59.4%) of which harbored actionable gene aberrations according to the Clinical Practice Guidelines for Next Generation Sequencing in Cancer Diagnosis and Treatment (Edition 1.0) issued by 3 major Japanese cancer-related societies. Twenty-five (13.3%) cases have since received molecular-targeted therapy according to their gene aberrations. These results indicate the utility of tumor-profiling multiplex gene panel testing in a clinical setting in Japan. This study is registered with UMIN Clinical Trials Registry (UMIN 000011141).

DNA Adductome Analysis Identifies N-Nitrosopiperidine Involved in the Etiology of Esophageal Cancer in Cixian, China.

Esophageal cancer is prevalent in Cixian, China, but the etiology of this disease remains largely unknown. Therefore, we explored this by conducting a DNA adductome analysis. Both tumorous and nontumorous tissues were collected from patients who underwent surgical procedures at Cixian Cancer Hospital and the Fourth Hospital of Hebei Medical University, which is in a low-incidence area. N2-(3,4,5,6-Tetrahydro-2H-pyran-2-yl)deoxyguanosine (THP-dG) was the major adduct detected in samples from esophageal cancer patients in Cixian. The precursor of THP-dG, N-nitrosopiperidine (NPIP), exhibited a strong mutagenic activity under metabolic activation in the Ames test and a significant dose-dependent increase in mutation frequency during an in vivo mutagenicity test with guanine phosphoribosyltransferase (gpt) delta rats. The NPIP-induced mutation was dominated by A:T to C:G transversions, followed by G:C to A:T and A:T to G:C transitions, in the liver and esophagus of animal samples. A similar mutational pattern was observed in the mutational signature of esophageal cancer patients that demonstrated weak correlation with THP-dG levels. These findings suggested that NPIP exposure is partly involved in the development of esophageal cancer in Cixian residents.

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.

Genomic and transcriptomic analysis of imatinib resistance in gastrointestinal stromal tumors.

Gastrointestinal stromal tumors represent the most common mesenchymal tumor of the digestive tract, driven by gain-of-function mutations in KIT. Despite its proven benefits, half of the patients treated with imatinib show disease progression within 2 years due to secondary resistance mutations in KIT. It remains unclear how the genomic and transcriptomic features change during the acquisition of imatinib resistance. Here, we performed exome sequencing and microarray transcription analysis for four imatinib-resistant cell lines and one cell line briefly exposed to imatinib. We also performed exome sequencing of clinical tumor samples. The cell line briefly exposed to imatinib exhibited few single-nucleotide variants and copy-number alterations, but showed marked upregulation of genes related to detoxification and downregulation of genes involved in cell cycle progression. Meanwhile, resistant cell lines harbored numerous genomic changes: amplified genes related to detoxification and deleted genes with cyclin-dependent kinase activity. Some variants in the resistant samples were traced back to the drug-sensitive samples, indicating the presence of ancestral subpopulations. The subpopulations carried variants associated with cell death. Pre-existing cancer cells with genetic alterations promoting apoptosis resistance may serve as a basis whereby cancer cells with critical mutations, such as secondary KIT mutations, can establish full imatinib resistance. © 2017 The Authors Genes, Chromosomes and Cancer Published by Wiley Periodicals, Inc.

Mismatch repair deficiency commonly precedes adenoma formation in Lynch Syndrome-Associated colorectal tumorigenesis.

Lynch syndrome is a cancer predisposition syndrome caused by germline mutations in mismatch repair (MMR) genes. MMR deficiency is a ubiquitous feature of Lynch syndrome-associated colorectal adenocarcinomas; however, it remains unclear when the MMR-deficient phenotype is acquired during tumorigenesis. To probe this issue, the present study examined genetic alterations and MMR statuses in Lynch syndrome-associated colorectal adenomas and adenocarcinomas, in comparison with sporadic adenomas. Among the Lynch syndrome-associated colorectal tumors, 68 of 86 adenomas (79%) and all adenocarcinomas were MMR-deficient, whereas all the sporadic adenomas were MMR-proficient, as determined by microsatellite instability testing and immunohistochemistry for MMR proteins. Sequencing analyses identified APC or CTNNB1 mutations in the majority of sporadic adenomas (58/84, 69%) and MMR-proficient Lynch syndrome-associated adenomas (13/18, 72%). However, MMR-deficient Lynch syndrome-associated adenomas had less APC or CTNNB1 mutations (25/68, 37%) and frequent frameshift RNF43 mutations involving mononucleotide repeats (45/68, 66%). Furthermore, frameshift mutations affecting repeat sequences constituted 14 of 26 APC mutations (54%) in MMR-deficient adenomas whereas these frameshift mutations were rare in MMR-proficient adenomas in patients with Lynch syndrome (1/12, 8%) and in sporadic adenomas (3/52, 6%). Lynch syndrome-associated adenocarcinomas exhibited mutation profiles similar to those of MMR-deficient adenomas. Considering that WNT pathway activation sufficiently drives colorectal adenoma formation, the distinct mutation profiles of WNT pathway genes in Lynch syndrome-associated adenomas suggest that MMR deficiency commonly precedes adenoma formation.

Genome-wide methylation profiles in primary intracranial germ cell tumors indicate a primordial germ cell origin for germinomas.

Intracranial germ cell tumors (iGCTs) are the second most common brain tumors among children under 14 in Japan. The World Health Organization classification recognizes several subtypes of iGCTs, which are conventionally subclassified into pure germinoma or non-germinomatous GCTs. Recent exhaustive genomic studies showed that mutations of the genes involved in the MAPK and/or PI3K pathways are common in iGCTs; however, the mechanisms of how different subtypes develop, often as a mixed-GCT, are unknown. To elucidate the pathogenesis of iGCTs, we investigated 61 GCTs of various subtypes by genome-wide DNA methylation profiling. We showed that pure germinomas are characterized by global low DNA methylation, a unique epigenetic feature making them distinct from all other iGCTs subtypes. The patterns of methylation strongly resemble that of primordial germ cells (PGC) at the migration phase, possibly indicating the cell of origin for these tumors. Unlike PGC, however, hypomethylation extends to long interspersed nuclear element retrotransposons. Histologically and epigenetically distinct microdissected components of mixed-GCTs shared identical somatic mutations in the MAPK or PI3K pathways, indicating that they developed from a common ancestral cell.

NCYM, a Cis-antisense gene of MYCN, encodes a de novo evolved protein that inhibits GSK3ß resulting in the stabilization of MYCN in human neuroblastomas.

The rearrangement of pre-existing genes has long been thought of as the major mode of new gene generation. Recently, de novo gene birth from non-genic DNA was found to be an alternative mechanism to generate novel protein-coding genes. However, its functional role in human disease remains largely unknown. Here we show that NCYM, a cis-antisense gene of the MYCN oncogene, initially thought to be a large non-coding RNA, encodes a de novo evolved protein regulating the pathogenesis of human cancers, particularly neuroblastoma. The NCYM gene is evolutionally conserved only in the taxonomic group containing humans and chimpanzees. In primary human neuroblastomas, NCYM is 100% co-amplified and co-expressed with MYCN, and NCYM mRNA expression is associated with poor clinical outcome. MYCN directly transactivates both NCYM and MYCN mRNA, whereas NCYM stabilizes MYCN protein by inhibiting the activity of GSK3ß, a kinase that promotes MYCN degradation. In contrast to MYCN transgenic mice, neuroblastomas in MYCN/NCYM double transgenic mice were frequently accompanied by distant metastases, behavior reminiscent of human neuroblastomas with MYCN amplification. The NCYM protein also interacts with GSK3ß, thereby stabilizing the MYCN protein in the tumors of the MYCN/NCYM double transgenic mice. Thus, these results suggest that GSK3ß inhibition by NCYM stabilizes the MYCN protein both in vitro and in vivo. Furthermore, the survival of MYCN transgenic mice bearing neuroblastoma was improved by treatment with NVP-BEZ235, a dual PI3K/mTOR inhibitor shown to destabilize MYCN via GSK3ß activation. In contrast, tumors caused in MYCN/NCYM double transgenic mice showed chemo-resistance to the drug. Collectively, our results show that NCYM is the first de novo evolved protein known to act as an oncopromoting factor in human cancer, and suggest that de novo evolved proteins may functionally characterize human disease.

Comprehensive screening of target molecules by next-generation sequencing in patients with malignant solid tumors: guiding entry into phase I clinical trials.

It is still controversial whether comprehensive genome screening of target molecules by next generation sequencing (NGS) is needed to increase clinical efficacy of investigational drugs or accelerate drug development, although several studies are being carried out. Therefore, we performed a prospective study to evaluate the feasibility of comprehensive gene screening in this setting. Our findings indicate that actionable alterations were identified in 45% of the analyzed patients, most frequently in those with breast cancer. Common actionable alterations were found in PIK3CA mutation, BRCA2 mutation, ERBB2 amplification, and CCND1 amplification. In total, 22% of the analyzed patients could be entered into phase I clinical trials, and 8% of them were treated with "matched" drugs. Among patients who received matched therapies, response and disease control rates were 33 and 78%, respectively. On the other hand, in the patients who received "non-matched" therapy, the objective response rate was 6%. We believe this data indicates that NGS-based molecular pre-screening is a potent platform for use before patient entry into phase I trials.