Long Noncoding RNA RP11-278A23.1, a Potential Modulator of p53 Tumor Suppression, Contributes to Colorectal Cancer Progression.

Recently, many studies revealed that long noncoding RNAs (lncRNAs) play important roles in cancers. To identify lncRNAs contributing to colorectal cancers, we screened lncRNAs through expression and survival analyses in datasets from The Cancer Genome Atlas (TCGA). The screen revealed that RP11-278A23.1 expression is significantly increased in colorectal cancer tissues compared with normal tissues and that high RP11-278A23.1 expression correlates with poor prognosis. The knockdown of RP11-278A23.1 inhibited the growth of and promoted apoptosis in colorectal cancer cells. Next, to comprehensively examine differentially expressed genes after RP11-278A23.1 knockdown, RNA sequencing was performed in HCT116 cells. The expression of p21, a p53 target gene, was significantly upregulated, and the expression of several p53 target proapoptotic genes was also altered. RP11-278A23.1 knockdown increased p53 expression at the translational level but not at the transcriptional level. Interestingly, RP11-278A23.1 knockdown also altered the expression of these proapoptotic genes in DLD1 cells with mutated p53 and in p53-knockout HCT116 cells. These results suggest that RP11-278A23.1 modifies the expression of these apoptosis-related genes in p53-dependent and p53-independent manners. In summary, lncRNA RP11-278A23.1 contributes to colorectal cancer progression by promoting cell growth and inhibiting apoptosis, suggesting that this lncRNA may be a useful therapeutic target.

The Clinical Validity of Urinary Pellet DNA Monitoring for the Diagnosis of Recurrent Bladder Cancer.

The aim of this study was to evaluate the clinical validity of monitoring urine pellet DNA (upDNA) of bladder cancer (BC) by digital PCR (dPCR) as a biomarker for early recurrence prediction, treatment efficacy evaluation, and no-recurrence corroboration. Tumor panel sequencing was first performed to select patient-unique somatic mutations to monitor both upDNA and circulating tumor DNA (ctDNA) by dPCR. For longitudinal monitoring using upDNA as well as plasma ctDNA, an average of 7.2 (range, 2 to 12) time points per case were performed with the dPCR assay for 32 previously treated and untreated patients with BC. Clinical recurrence based on imaging and urine cytology was compared using upDNA variant allele frequency (VAF) dynamics. A continuous increasing trend of upDNA VAF ≥1% was considered to indicate molecular recurrence. Most (30/32; 93.8%) cases showed at least one traceable somatic mutation. In 5 of 7 cases (71.4%) with clinical recurrence, upDNA VAF >1% was detected 7 to 15 months earlier than the imaging diagnosis. The upDNA VAF remained high after initial treatment for locally recurrent cases. The clinical validity of upDNA monitoring was confirmed with the observation that 26 of 30 cases (86.7%) were traceable. Local recurrences were not indicated by ctDNA alone. The results support the clinical validity of upDNA monitoring in the management of recurrent BC.

The frequency and pathogenicity of BRCA1 and BRCA2 variants in the general Japanese population.

Hereditary breast and ovarian cancer syndrome (HBOC) resulting from pathogenic variants of BRCA1 or BRCA2 is the most common and well-documented hereditary tumor. Although founder variants have been identified in population-based surveys in various countries, the types of variants are not uniform across races and regions. Recently, the Tohoku Medical Megabank Organization (ToMMo) released whole-genome sequence data including approximately 54,000 individuals from the general population of the Tohoku area in Japan. We analyzed these data and comprehensively identified the prevalence of BRCA1/2 pathogenic and truncating variants. We believe that an accurate understanding of the unique distribution and characteristics of pathogenic BRCA1/2 variants in Japan through this analysis will enable better surveillance and intervention for HBOC patients, not only in Japan but also worldwide.

Effects of temozolomide on tumor mutation burden and microsatellite instability in melanoma cells.

The efficacy of combination therapy with an immune checkpoint inhibitor (ICI) and cytotoxic chemotherapeutic agents has been investigated in cancer, including melanoma. Before ICIs were introduced, dacarbazine or temozolomide (TMZ) were used to treat melanoma. Several studies using glioma or colorectal cancer cells showed that TMZ can increase the tumor mutation burden (TMB) and induce mismatch repair (MMR) deficiency associated with microsatellite instability (MSI). These could increase immunoreactivity to an ICI, but this has not been evaluated in melanoma cells. We investigated the effects of TMZ on MSI status and TMB in melanoma cells. To evaluate the TMB, we performed whole-exome sequencing using genomic DNA from the human melanoma cell lines Mel18, A375, WM266-4, G361, and TXM18 before and after TMZ treatment. Polymerase chain reaction amplification of five mononucleotide repeat markers, BAT25, BAT26, NR21, NR24, and MONO27, was performed, and we analyzed changes in the MSI status. In all cell lines, the TMB was increased after TMZ treatment (the change amount of TMB with ≤ 5% variant allele frequency [VAF] was 18.0-38.3 mutations per megabase) even in the condition without obvious cytological damage. MSI after TMZ treatment was not observed in any cells. TMZ increased TMB but did not change MSI status in melanoma cells.

LMNA::NTRK1 Fusion-positive Leiomyosarcoma: Discrepancy between DNA-based Comprehensive Genomic Profiling and RNA Sequencing.

A 26-year-old man presented with a tumor in the left soleus muscle. The tumor was diagnosed as a locally advanced leiomyosarcoma. The patient was treated with irradiation followed by wide resection. One year after surgery, the patient presented with multiple lung metastases. Despite aggressive sequential chemotherapy, systemic metastatic tumors continued to develop. To explore therapeutic options for the patient, we performed DNA-based CGP with FoundationOne® CDx (F1). F1 identified anout-of-strand rearrangement of the NOS1AP::NTRK1 gene, which has not been previously reported. In contrast, RNA sequencing revealed an in-frame LMNA::NTRK1 gene, which is an oncogenic fusion gene.

TM4SF1-AS1 inhibits apoptosis by promoting stress granule formation in cancer cells.

Long noncoding RNAs (lncRNAs) play pivotal roles in tumor development. To identify dysregulated lncRNAs in gastric cancer (GC), we analyzed genome-wide trimethylation of histone H3 lysine 4 (H3K4me3) to screen for transcriptionally active lncRNA genes in the non-tumorous gastric mucosa of patients with GC and healthy individuals. We found that H3K4me3 at TM4SF1-AS1 was specifically upregulated in GC patients and that the expression of TM4SF1-AS1 was significantly elevated in primary and cultured GC cells. TM4SF1-AS1 contributes to GC cell growth in vitro and in vivo, and its oncogenic function is mediated, at least in part, through interactions with purine-rich element-binding protein α (Pur-α) and Y-box binding protein 1 (YB-1). TM4SF1-AS1 also activates interferon signaling in GC cells, which is dependent on Pur-α and RIG-I. Chromatin isolation by RNA purification (ChIRP)-mass spectrometry demonstrated that TM4SF1-AS1 was associated with several stress granule (SG)-related proteins, including G3BP2, RACK1, and DDX3. Notably, TM4SF1-AS1 promoted SG formation and inhibited apoptosis in GC cells by sequestering RACK1, an activator of the stress-responsive MAPK pathway, within SGs. TM4SF1-AS1-induced SG formation and apoptosis inhibition are dependent on Pur-α and YB-1. These findings suggested that TM4SF1-AS1 contributes to tumorigenesis by enhancing SG-mediated stress adaptation.

MYEOV overexpression induced by demethylation of its promoter contributes to pancreatic cancer progression via activation of the folate cycle/c-Myc/mTORC1 pathway.

BACKGROUND: While molecular targeted drugs and other therapies are being developed for many tumors, pancreatic cancer is still considered to be the malignant tumor with the worst prognosis. We started this study to identify prognostic genes and therapeutic targets of pancreatic cancer. METHODS: To comprehensively identify prognostic genes in pancreatic cancer, we investigated the correlation between gene expression and cancer-specific prognosis using transcriptome and clinical information datasets from The Cancer Genome Atlas (TCGA). In addition, we examined the effects of the suppression of candidate prognostic genes in pancreatic cancer cell lines. RESULT: We found that patients with high expression levels of MYEOV, a primate-specific gene with unknown function, had significantly shorter disease-specific survival times than those with low expression levels. Cox proportional hazards analysis revealed that high expression of MYEOV was significantly associated with poor survival and was an independent prognostic factor for disease-specific survival in pancreatic cancer patients. Analysis of multiple cancer samples revealed that the MYEOV promoter region is methylated in noncancer tissues but is demethylated in tumors, causing MYEOV overexpression in tumors. Notably, the knockdown of MYEOV suppressed the expression of MTHFD2 and other folate metabolism-related enzyme genes required for the synthesis of amino acids and nucleic acids and also restored the expression of c-Myc and mTORC1 repressors. CONCLUSION: There is a significant correlation between elevated MYEOV expression and poor disease-specific survival in pancreatic cancer patients. MYEOV enhances the activation of several oncogenic pathways, resulting in the induction of pancreatic cancer cell proliferation. Overall, MYEOV acts as an oncogene in pancreatic cancer. Furthermore, MYEOV may be a prognostic biomarker and serve as an 'actionable' therapeutic target for pancreatic cancers.

Multiomics identifies the link between intratumor steatosis and the exhausted tumor immune microenvironment in hepatocellular carcinoma.

BACKGROUND AND AIMS: Immunotherapy has become the standard-of-care treatment for hepatocellular carcinoma (HCC), but its efficacy remains limited. To identify immunotherapy-susceptible HCC, we profiled the molecular abnormalities and tumor immune microenvironment (TIME) of rapidly increasing nonviral HCC. APPROACHES AND RESULTS: We performed RNA-seq of tumor tissues in 113 patients with nonviral HCC and cancer genome sequencing of 69 genes with recurrent genetic alterations reported in HCC. Unsupervised hierarchical clustering classified nonviral HCCs into three molecular classes (Class I, II, III), which stratified patient prognosis. Class I, with the poorest prognosis, was associated with TP53 mutations, whereas class III, with the best prognosis, was associated with cadherin-associated protein beta 1 (CTNNB1) mutations. Thirty-eight percent of nonviral HCC was defined as an immune class characterized by a high frequency of intratumoral steatosis and a low frequency of CTNNB1 mutations. Steatotic HCC, which accounts for 23% of nonviral HCC cases, presented an immune-enriched but immune-exhausted TIME characterized by T cell exhaustion, M2 macrophage and cancer-associated fibroblast (CAF) infiltration, high PD-L1 expression, and TGF-ß signaling activation. Spatial transcriptome analysis suggested that M2 macrophages and CAFs may be in close proximity to exhausted CD8+ T cells in steatotic HCC. An in vitro study showed that palmitic acid-induced lipid accumulation in HCC cells upregulated PD-L1 expression and promoted immunosuppressive phenotypes of cocultured macrophages and fibroblasts. Patients with steatotic HCC, confirmed by chemical-shift MR imaging, had significantly longer PFS with combined immunotherapy using anti-PD-L1 and anti-VEGF antibodies. CONCLUSIONS: Multiomics stratified nonviral HCCs according to prognosis or TIME. We identified the link between intratumoral steatosis and immune-exhausted immunotherapy-susceptible TIME.

Individualized circulating tumor DNA monitoring in head and neck squamous cell carcinoma.

There is no useful biomarker to evaluate treatment response and early relapse in head and neck squamous cell carcinoma (HNSCC). Circulating tumor DNA (ctDNA) is a promising biomarker for detecting minimal residual diseases and monitoring treatment effect. We investigated whether individualized ctDNA analysis could help monitor treatment response and relapse in HNSCC. Mutation analysis of tumor and peripheral blood mononuclear cell (PBMC) DNAs of 26 patients with HNSCC was performed using a custom squamous cell carcinoma (SCC) panel. The identified individualized mutated genes were defined as ctDNA candidates. We investigated whether frequent ctDNA monitoring via digital PCR (dPCR) is clinically valid for HNSCC patients. TP53 was the most frequently mutated gene and was detected in 14 of 24 cases (58.2%), wherein two cases were excluded owing to the absence of tumor-specific mutations in the SCC panel. Six cases were excluded because of undesignable and unusable primer-probes for dPCR. Longitudinal ctDNA was monitored in a total of 18 cases. In seven cases, ctDNA tested positive again or did not test negative, and all seven cases relapsed after initial curative treatment. In 11 cases, after initial curative treatment, ctDNA remained negative and patients were alive without recurrence. Patients who remained negative for ctDNA during follow-up after initial curative treatment (n = 11) had a significantly better prognosis than those who reverted to ctDNA positivity (n = 7; p < 0.0001; log-rank test). Individualized ctDNA monitoring using SCC panel and dPCR might be a novel and promising biomarker for HNSCC.

DLEU1 promotes oral squamous cell carcinoma progression by activating interferon-stimulated genes.

Long noncoding RNAs (lncRNAs) are deeply involved in cancer development. We previously reported that DLEU1 (deleted in lymphocytic leukemia 1) is one of the lncRNAs overexpressed in oral squamous cell carcinoma (OSCC) cells, where it exhibits oncogenic activity. In the present study, we further clarified the molecular function of DLEU1 in the pathogenesis of OSCC. Chromatin immunoprecipitation-sequencing (ChIP-seq) analysis revealed that DLEU1 knockdown induced significant changes in the levels of histone H3 lysine 4 trimethylation (H3K4me3) and H3K27 acetylation (H3K27ac) in OSCC cells. Notably, DLEU1 knockdown suppressed levels of H3K4me3/ H3K27ac and expression of a number of interferon-stimulated genes (ISGs), including IFIT1, IFI6 and OAS1, while ectopic DLEU1 expression activated these genes. Western blot analysis and reporter assays suggested that DLEU1 upregulates ISGs through activation of JAK-STAT signaling in OSCC cells. Moreover, IFITM1, one of the ISGs induced by DLUE1, was frequently overexpressed in primary OSCC tumors, and its knockdown inhibited OSCC cell proliferation, migration and invasion. These findings suggest that DLEU1 exerts its oncogenic effects, at least in part, through activation of a series ISGs in OSCC cells.

Dead or alive? Pitfall of survival analysis with TCGA datasets.

We often encounter situations in which data from the TCGA that have been analyzed in papers we read or reviewed cannot be reproduced, even when TCGA datasets are used, especially in survival analyses. Therefore, we attempted to confirm the data source for TCGA survival analysis and found that several websites used to analyze the survival data of TCGA datasets inappropriately handle the survival data, causing differences in statistical analyses. This causes the misinterpretation of results because figures of survival analysis results in several papers are sometimes exactly as generated by these sites, and the results depend on only the tools provided by these sites. We would like to make this situation widely known and raise the problem for scientific soundness.

Early dynamics of circulating tumor DNA predict chemotherapy responses for patients with esophageal cancer.

We investigated whether early circulating tumor DNA (ctDNA) changes, measured using digital PCR (dPCR), can predict later chemotherapy responses in esophageal squamous cell cancer (ESCC). We compared the dynamics of ctDNA and tumor volumes during chemotherapy in 42 ESCC. The accuracy of predictions of later chemotherapy responses was evaluated by the ratio of the variant allele frequency of ctDNA (post-/pre-ctDNA) and the total tumor volume (post-/pre-volume) before and after an initial chemotherapy cycle using a receiver-operating characteristic curve analysis. Total positive and negative objective responses (ORs) were defined as either >50 or ≤50% reductions, respectively, in the total tumor volume at the end of first-line chemotherapy. Mutation screening of 43 tumors from 42 patients revealed 96 mutations. The pretreatment dPCR-ctDNA data were informative in 38 patients, using 70 selected mutations (1-3 per patient). The areas under the curve (AUCs) for the post-/pre-volume and post-/pre-ctDNA levels used in predicting the total OR were 0.85 and 0.88, respectively. The optimal cutoff value of post-/pre-ctDNA was 0.13. In 20 patients with post-/pre-volume ≥50%, the total OR could be predicted by the post-/pre-ctDNA with high accuracy; the AUC by post-/pre-ctDNA was higher than that by post-/pre-volume (0.85 versus 0.76, respectively). Patients with low post-/pre-ctDNA (n = 18) had a significantly better overall survival rate than those with high post-/pre-ctDNA (n = 20; P = 0.03). Early ctDNA changes after an initial cycle of chemotherapy predict later responses to treatment with high accuracy in ESCC patients.

Generation of functional liver organoids on combining hepatocytes and cholangiocytes with hepatobiliary connections ex vivo.

In the liver, the bile canaliculi of hepatocytes are connected to intrahepatic bile ducts lined with cholangiocytes, which remove cytotoxic bile from the liver tissue. Although liver organoids have been reported, it is not clear whether the functional connection between hepatocytes and cholangiocytes is recapitulated in those organoids. Here, we report the generation of a hepatobiliary tubular organoid (HBTO) using mouse hepatocyte progenitors and cholangiocytes. Hepatocytes form the bile canalicular network and secrete metabolites into the canaliculi, which are then transported into the biliary tubular structure. Hepatocytes in HBTO acquire and maintain metabolic functions including albumin secretion and cytochrome P450 activities, over the long term. In this study, we establish functional liver tissue incorporating a bile drainage system ex vivo. HBTO enable us to reproduce the transport of hepatocyte metabolites in liver tissue, and to investigate the way in which the two types of epithelial cells establish functional connections.

Genetic analyses of a secondary poroma and trichoblastoma in a HRAS-mutated sebaceous nevus.

A sebaceous nevus is a congenital skin hamartoma caused by postzygotic HRAS or KRAS mosaic mutations. With age, affected individuals may develop secondary tumors within a sebaceous nevus. RAS mutations are harbored from the onset of sebaceous nevus, and further mutations can be expected to be required in order to explain the initiation of secondary tumors. However, genetic analyses of the secondary tumors have not been conducted. Herein, we describe the rare coexistence of a poroma and a trichoblastoma arising in a sebaceous nevus. This is the first report of an investigation of multiple genes in a secondary tumor in an SN. First, HRAS c.37G>C, which is the common mutation in sebaceous nevus, was detected in all three lesions (sebaceous nevus, poroma, and trichoblastoma). Next, to elucidate the potential second-hit mutations in the secondary poroma and trichoblastoma, we applied a panel sequencing for skin cancers that was newly developed in our institution. Our comparison of the mutational profile of 95 skin cancer-related genes in each of the three lesions newly revealed TP53 p.R158P in the poroma and NOTCH2 p.G329S in the trichoblastoma. TP53 p.R158P has been determined as a pathogenic mutation in other tumors, and NOTCH2 p.G329S was a novel mutation. We identified two novel mutations that may have contributed to the pathogenesis of the secondary tumor's development. The roles of the mutations remain unclear.

Mutation Analysis of Radioresistant Early-Stage Cervical Cancer.

Radiotherapy is a definitive treatment for early-stage cervical cancer; however, a subset of this disease recurs locally, necessitating establishment of predictive biomarkers and treatment strategies. To address this issue, we performed gene panel-based sequencing of 18 stage IB cervical cancers treated with definitive radiotherapy, including two cases of local recurrence, followed by in vitro and in silico analyses. Simultaneous mutations in KRAS and SMAD4 (KRASmt/SMAD4mt) were detected only in a local recurrence case, indicating potential association of this mutation signature with radioresistance. In isogenic cell-based experiments, a combination of activating KRAS mutation and SMAD4 deficiency led to X-ray resistance, whereas either of these factors alone did not. Analysis of genomic data from 55,308 cancers showed a significant trend toward co-occurrence of mutations in KRAS and SMAD4. Gene Set Enrichment Analysis of the Cancer Cell Line Encyclopedia dataset suggested upregulation of the pathways involved in epithelial mesenchymal transition and inflammatory responses in KRASmt/SMAD4mt cancer cells. Notably, irradiation with therapeutic carbon ions led to robust killing of X-ray-resistant KRASmt/SMAD4mt cancer cells. These data indicate that the KRASmt/SMAD4mt signature is a potential predictor of radioresistance, and that carbon ion radiotherapy is a potential option to treat early-stage cervical cancers with the KRASmt/SMAD4mt signature.

Whole-exome sequencing of 79 xenografts as a potential approach for the identification of genetic variants associated with sensitivity to cytotoxic anticancer drugs.

Chemotherapy response remains unpredictable in most patients with cancer. In this study, we performed whole-exome sequencing of 79 cancer xenografts derived from human cancer tissues to identify genetic predictors of chemosensitivity to nine cytotoxic anticancer drugs. Xenografts were harvested from 12 organs with cancer and implanted into nude mice. The mice were exposed to one of nine cytotoxic anticancer drugs (5-fluorouracil, nimustine, adriamycin, cyclophosphamide, cisplatin, mitomycin C, methotrexate, vincristine, and vinblastine) to assess the correlation between chemosensitivity response and variant allele frequency. We found 162 candidate variants that were possibly associated with chemosensitivity to one or more of the nine anticancer drugs (P < 0.01). In a subgroup analysis of breast and gastric cancer xenografts, 78 and 67 variants, respectively, were possibly associated with chemosensitivity. This approach may help to contribute to the development of personalized treatments that may allow for the prescription of optimal chemotherapy regimens among patients with cancer.