Postoperative recurrence detection using individualized circulating tumor DNA in upper tract urothelial carcinoma.

Biomarkers that could detect the postoperative recurrence of upper tract urothelial carcinoma (UTUC) have not been established. In this prospective study, we aim to evaluate the utility of individualized circulating tumor DNA (ctDNA) monitoring using digital PCR (dPCR) as a tumor recurrence biomarker for UTUC in the perioperative period. Twenty-three patients who underwent radical nephroureterectomy (RNU) were included. In each patient, whole exome sequencing by next-generation sequencing and TERT promoter sequencing of tumor DNA were carried out. Case-specific gene mutations were selected from sequencing analysis to examine ctDNA by dPCR analysis. We also prospectively collected plasma and urine ctDNA from each patient. The longitudinal variant allele frequencies of ctDNA during the perioperative period were plotted. Case-specific gene mutations were detected in 22 cases (96%) from ctDNA in the preoperative samples. Frequently detected genes were TERT (39%), FGFR3 (26%), TP53 (22%), and HRAS (13%). In all cases, we obtained plasma and urine samples for 241 time points and undertook individualized ctDNA monitoring for 2 years after RNU. Ten patients with intravesical recurrence had case-specific ctDNA detected in urine at the time of recurrence. The mean lead time of urinary ctDNA in intravesical recurrence was 60 days (range, 0-202 days). Two patients with distal metastasis had case-specific ctDNA in plasma at the time of metastasis. In UTUC, tumor-specific gene mutations can be monitored postoperatively as ctDNA in plasma and urine. Individualized ctDNA might be a minimally invasive biomarker for the early detection of postoperative recurrence.

Immuno-genomic analysis reveals eosinophilic feature and favorable prognosis of female non-smoking esophageal squamous cell carcinomas.

Most of esophageal squamous cell carcinoma (ESCC) develop in smoking males in Japan, but the genomic etiology and immunological characteristics of rare non-smoking female ECSS remain unclear. To elucidate the genomic and immunological features of ESCC in non-smoking females, we analyzed whole-genome or transcriptome sequencing data from 94 ESCCs, including 20 rare non-smoking female cases. In addition, 31,611 immune cells were extracted from four ESCC tissues and subject to single-cell RNA-seq. We compared their immuno-genomic and microbiome profiles between non-smoking female and smoking ESCCs. Non-smoking females showed much better prognosis. Whole-genome sequencing analysis showed no significant differences in driver genes or copy number alterations depending on smoking status. The mutational signatures specifically observed in non-smoking females ESCC could be attributed to aging. Immune profiling from RNA-seq revealed that ESCC in non-smoking females had high tumor microenvironment signatures and a high abundance of eosinophils with a favorable prognosis. Single-cell RNA-sequencing of intratumor immune cells revealed gender differences of eosinophils and their activation in female cases. ESCCs in non-smoking females have age-related mutational signatures and gender-specific tumor immune environment with eosinophils, which is likely to contribute to their favorable prognosis.

Mitochondrial DNA Variants at Low-Level Heteroplasmy and Decreased Copy Numbers in Chronic Kidney Disease (CKD) Tissues with Kidney Cancer.

The human mitochondrial genome (mtDNA) is a circular DNA molecule with a length of 16.6 kb, which contains a total of 37 genes. Somatic mtDNA mutations accumulate with age and environmental exposure, and some types of mtDNA variants may play a role in carcinogenesis. Recent studies observed mtDNA variants not only in kidney tumors but also in adjacent kidney tissues, and mtDNA dysfunction results in kidney injury, including chronic kidney disease (CKD). To investigate whether a relationship exists between heteroplasmic mtDNA variants and kidney function, we performed ultra-deep sequencing (30,000×) based on long-range PCR of DNA from 77 non-tumor kidney tissues of kidney cancer patients with CKD (stages G1 to G5). In total, this analysis detected 697 single-nucleotide variants (SNVs) and 504 indels as heteroplasmic (0.5% ≤ variant allele frequency (VAF) < 95%), and the total number of detected SNVs/indels did not differ between CKD stages. However, the number of deleterious low-level heteroplasmic variants (pathogenic missense, nonsense, frameshift and tRNA) significantly increased with CKD progression (p < 0.01). In addition, mtDNA copy numbers (mtDNA-CNs) decreased with CKD progression (p < 0.001). This study demonstrates that mtDNA damage, which affects mitochondrial genes, may be involved in reductions in mitochondrial mass and associated with CKD progression and kidney dysfunction.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma.

BACKGROUND: Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated. METHODS: To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours. FINDINGS: RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients. INTERPRETATION: These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology. FUNDING: This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Variant spectrum of von Hippel-Lindau disease and its genomic heterogeneity in Japan.

Von Hippel-Lindau (VHL) disease is an autosomal dominant, inherited syndrome with variants in the VHL gene, causing predisposition to multi-organ neoplasms with vessel abnormality. Germline variants in VHL can be detected in 80-90% of patients clinically diagnosed with VHL disease. Here, we summarize the results of genetic tests for 206 Japanese VHL families, and elucidate the molecular mechanisms of VHL disease, especially in variant-negative unsolved cases. Of the 206 families, genetic diagnosis was positive in 175 families (85%), including 134 families (65%) diagnosed by exon sequencing (15 novel variants) and 41 (20%) diagnosed by multiplex ligation-dependent probe amplification (MLPA) (one novel variant). The deleterious variants were significantly enriched in VHL disease Type 1. Interestingly, five synonymous or non-synonymous variants within exon 2 caused exon 2 skipping, which is the first report of exon 2 skipping caused by several missense variants. Whole genome and target deep sequencing analysis were performed for 22 unsolved cases with no variant identified and found three cases with VHL mosaicism (variant allele frequency: 2.5-22%), one with mobile element insertion in the VHL promoter region, and two with a pathogenic variant of BAP1 or SDHB. The variants associated with VHL disease are heterogeneous, and for more accuracy of the genetic diagnosis of VHL disease, comprehensive genome and DNA/RNA analyses are required to detect VHL mosaicism, complicated structure variants and other related gene variants.

Hereditary cancer variants and homologous recombination deficiency in biliary tract cancer.

BACKGROUND & AIMS: The heritability and actionability of variants in homologous recombination-related genes in biliary tract cancers (BTCs) are uncertain. Although associations between BTC and BRCA germline variants have been reported, homologous recombination deficiency has not been investigated in BTCs. METHODS: We sequenced germline variants in 27 cancer-predisposing genes in 1,292 BTC cases and 37,583 controls without a personal nor family history of cancer. We compared pathogenic germline variant frequencies between cases and controls and documented the demographic and clinical characteristics of carriers. In addition, whole-genome sequencing of 45 BTC tissues was performed to evaluate homologous recombination deficiency status. RESULTS: Targeted sequencing identified 5,018 germline variants, which were classified into 317 pathogenic, 3,611 variants of uncertain significance, and 1,090 benign variants. Seventy-one BTC cases (5.5%) had at least one pathogenic variant among 27 cancer-predisposing genes. Pathogenic germline variants enriched in BTCs were present in BRCA1, BRCA2, APC, and MSH6 (p <0.00185). PALB2 variants were marginally associated with BTC (p = 0.01). APC variants were predominantly found in ampulla of Vater carcinomas. Whole-genome sequencing demonstrated that three BTCs with pathogenic germline variants in BRCA2 and PALB2, accompanied by loss of heterozygosity, displayed homologous recombination deficiency. Conversely, pathogenic germline variants without a second hit or variants of other homologous recombination-related genes such as ATM and BRIP1 showed homologous recombination-proficient phenotypes. CONCLUSIONS: In this study, we describe the heritability and actionability of variants in homologous recombination-related genes, which could be used to guide screening and therapeutic strategies for BTCs. IMPACT AND IMPLICATIONS: We found that 5.5% of biliary tract cancers (BTCs) in a Japanese population possessed hereditary cancer-predisposing gene alterations, including in BRCA and genes associated with colorectal cancer. Two hits in homologous recombination-related genes were required to confer a homologous recombination-deficient phenotype. PARP inhibitors and DNA-damaging regimens may be effective strategies against BTCs exhibiting homologous recombination deficiency. Hence, in this study, genome-wide sequencing has revealed a potential new therapeutic strategy that could be applied to a subset of BTCs.

Genomic features of renal cell carcinoma developed during end-stage renal disease and dialysis.

Patients with end-stage renal disease (ESRD) or receiving dialysis have a much higher risk for renal cell carcinoma (RCC), but carcinogenic mechanisms and genomic features remain little explored and undefined. This study's goal was to identify the genomic features of ESRD RCC and characterize them for associations with tumor histology and dialysis exposure. In this study, we obtained 33 RCCs, with various histological subtypes, that developed in ESRD patients receiving dialysis and performed whole-genome sequencing and transcriptome analyses. Driver events, copy-number alteration (CNA) analysis and mutational signature profiling were performed using an analysis pipeline that integrated data from germline and somatic SNVs, Indels and structural variants as well as CNAs, while transcriptome data were analyzed for differentially expressed genes and through gene set enrichment analysis. ESRD related clear cell RCCs' driver genes and mutations mirrored those in sporadic ccRCCs. Longer dialysis periods significantly correlated with a rare mutational signature SBS23, whose etiology is unknown, and increased mitochondrial copy number. All acquired cystic disease (ACD)-RCCs, which developed specifically in ESRD patients, showed chromosome 16q amplification. Gene expression analysis suggests similarity between certain ACD-RCCs and papillary RCCs and in TCGA papillary RCCs with chromosome 16 gain identified enrichment for genes related to DNA repair, as well as pathways related to reactive oxygen species, oxidative phosphorylation and targets of Myc. This analysis suggests that ESRD or dialysis could induce types of cellular stress that impact some specific types of genomic damage leading to oncogenesis.

Proteo-genomic characterization of virus-associated liver cancers reveals potential subtypes and therapeutic targets.

Primary liver cancer is a heterogeneous disease in terms of its etiology, histology, and therapeutic response. Concurrent proteomic and genomic characterization of a large set of clinical liver cancer samples can help elucidate the molecular basis of heterogeneity and thus serve as a valuable resource for personalized liver cancer treatment. In this study, we perform proteomic profiling of ~300 proteins on 259 primary liver cancer tissues with reverse-phase protein arrays, mutational analysis using whole genome sequencing and transcriptional analysis with RNA-Seq. Patients are of Japanese ethnic background and mainly HBV or HCV positive, providing insight into this important liver cancer subtype. Unsupervised classification of tumors based on protein expression profiles reveal three proteomic subclasses R1, R2, and R3. The R1 subclass is immunologically hot and demonstrated a good prognosis. R2 contains advanced proliferative tumor with TP53 mutations, high expression of VEGF receptor 2 and the worst prognosis. R3 is enriched with CTNNB1 mutations and elevated mTOR signaling pathway activity. Twenty-two proteins, including CDK1 and CDKN2A, are identified as potential prognostic markers. The proteomic classification presented in this study can help guide therapeutic decision making for liver cancer treatment.

Immuno-genomic profiling of biopsy specimens predicts neoadjuvant chemotherapy response in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is one of the most aggressive cancers and is primarily treated with platinum-based neoadjuvant chemotherapy (NAC). Some ESCCs respond well to NAC. However, biomarkers to predict NAC sensitivity and their response mechanism in ESCC remain unclear. We perform whole-genome sequencing and RNA sequencing analysis of 141 ESCC biopsy specimens before NAC treatment to generate a machine-learning-based diagnostic model to predict NAC reactivity in ESCC and analyzed the association between immunogenomic features and NAC response. Neutrophil infiltration may play an important role in ESCC response to NAC. We also demonstrate that specific copy-number alterations and copy-number signatures in the ESCC genome are significantly associated with NAC response. The interactions between the tumor genome and immune features of ESCC are likely to be a good indicator of therapeutic capability and a therapeutic target for ESCC, and machine learning prediction for NAC response is useful.

Actionability evaluation of biliary tract cancer by genome transcriptome analysis and Asian cancer knowledgebase.

INTRODUCTION: Treatment options for biliary tract cancer (BTC) are very limited. It is necessary to investigate actionable genes and candidate drugs using a sophisticated knowledgebase (KB) and characterize BTCs immunologically for evaluating the actionability of molecular and immune therapies. MATERIALS AND METHODS: The genomic and transcriptome data of 219 patients with BTC who underwent surgery were analyzed. Actionable mutations and candidate drugs were annotated using the largest available KB of the Asian population (CancerSCAN®). Predictive biomarkers of immune checkpoint inhibitors were analyzed using DNA and RNA sequencing data. RESULTS: Twenty-two actionable genes and 43 candidate drugs were annotated in 74 patients (33.8%). The most frequent actionable genes were PTEN (7.3%), CDKN2A (6.8%), KRAS (6.4%). BRCA2, CDKN2A, and FGFR2 mutations were most frequently identified in case of intrahepatic cholangiocarcinoma. PTEN and CDKN2A mutations were associated with significantly shorter overall survival. PD-L1 and PD-1 expression was significantly higher in case of extrahepatic cholangiocarcinoma and T-cell-high expression. In total, 49.7% of cases were evaluated as having actionability for molecular therapy or immune checkpoint inhibitors. CONCLUSIONS: Identifying actionable genes and candidate drugs using the KB contribute to the development of therapeutic drugs and personalized treatment for BTC.

Molecular Classification and Tumor Microenvironment Characterization of Gallbladder Cancer by Comprehensive Genomic and Transcriptomic Analysis.

Gallbladder cancer (GBC), a rare but lethal disease, is often diagnosed at advanced stages. So far, molecular characterization of GBC is insufficient, and a comprehensive molecular portrait is warranted to uncover new targets and classify GBC. We performed a transcriptome analysis of both coding and non-coding RNAs from 36 GBC fresh-frozen samples. The results were integrated with those of comprehensive mutation profiling based on whole-genome or exome sequencing. The clustering analysis of RNA-seq data facilitated the classification of GBCs into two subclasses, characterized by high or low expression levels of TME (tumor microenvironment) genes. A correlation was observed between gene expression and pathological immunostaining. TME-rich tumors showed significantly poor prognosis and higher recurrence rate than TME-poor tumors. TME-rich tumors showed overexpression of genes involved in epithelial-to-mesenchymal transition (EMT) and inflammation or immune suppression, which was validated by immunostaining. One non-coding RNA, miR125B1, exhibited elevated expression in stroma-rich tumors, and miR125B1 knockout in GBC cell lines decreased its invasion ability and altered the EMT pathway. Mutation profiles revealed TP53 (47%) as the most commonly mutated gene, followed by ELF3 (13%) and ARID1A (11%). Mutations of ARID1A, ERBB3, and the genes related to the TGF-ß signaling pathway were enriched in TME-rich tumors. This comprehensive analysis demonstrated that TME, EMT, and TGF-ß pathway alterations are the main drivers of GBC and provides a new classification of GBCs that may be useful for therapeutic decision-making.

An Organoid Biobank of Neuroendocrine Neoplasms Enables Genotype-Phenotype Mapping.

Gastroenteropancreatic (GEP) neuroendocrine neoplasm (NEN) that consists of neuroendocrine tumor and neuroendocrine carcinoma (NEC) is a lethal but under-investigated disease owing to its rarity. To fill the scarcity of clinically relevant models of GEP-NEN, we here established 25 lines of NEN organoids and performed their comprehensive molecular characterization. GEP-NEN organoids recapitulated pathohistological and functional phenotypes of the original tumors. Whole-genome sequencing revealed frequent genetic alterations in TP53 and RB1 in GEP-NECs, and characteristic chromosome-wide loss of heterozygosity in GEP-NENs. Transcriptome analysis identified molecular subtypes that are distinguished by the expression of distinct transcription factors. GEP-NEN organoids gained independence from the stem cell niche irrespective of genetic mutations. Compound knockout of TP53 and RB1, together with overexpression of key transcription factors, conferred on the normal colonic epithelium phenotypes that are compatible with GEP-NEN biology. Altogether, our study not only provides genetic understanding of GEP-NEN, but also connects its genetics and biological phenotypes.

Whole genome sequencing analysis identifies recurrent structural alterations in esophageal squamous cell carcinoma.

Esophageal squamous cell carcinoma (ESCC) is the predominant type of esophageal cancer in the Asian region, including Japan. A previous study reported mutational landscape of Japanese ESCCs by using exome sequencing. However, somatic structural alterations were yet to be explored. To provide a comprehensive mutational landscape, we performed whole genome sequencing (WGS) analysis of biopsy specimens from 20 ESCC patients in a Japanese population. WGS analysis identified non-silent coding mutations of TP53, ZNF750 and FAT1 in ESCC. We detected six mutational signatures in ESCC, one of which showed significant association with smoking status. Recurrent structural variations, many of which were chromosomal deletions, affected genes such as LRP1B, TTC28, CSMD1, PDE4D, SDK1 and WWOX in 25%-30% of tumors. Somatic copy number amplifications at 11q13.3 (CCND1), 3q26.33 (TP63/SOX2), and 8p11.23 (FGFR1) and deletions at 9p21.3 (CDKN2A) were identified. Overall, these multi-dimensional view of genomic alterations improve the understanding of the ESCC development at molecular level and provides future prognosis and therapeutic implications for ESCC in Japan.

Genomic and Transcriptomic Profiling of Combined Hepatocellular and Intrahepatic Cholangiocarcinoma Reveals Distinct Molecular Subtypes.

We performed genomic and transcriptomic sequencing of 133 combined hepatocellular and intrahepatic cholangiocarcinoma (cHCC-ICC) cases, including separate, combined, and mixed subtypes. Integrative comparison of cHCC-ICC with hepatocellular carcinoma and intrahepatic cholangiocarcinoma revealed that combined and mixed type cHCC-ICCs are distinct subtypes with different clinical and molecular features. Integrating laser microdissection, cancer cell fraction analysis, and single nucleus sequencing, we revealed both mono- and multiclonal origins in the separate type cHCC-ICCs, whereas combined and mixed type cHCC-ICCs were all monoclonal origin. Notably, cHCC-ICCs showed significantly higher expression of Nestin, suggesting Nestin may serve as a biomarker for diagnosing cHCC-ICC. Our results provide important biological and clinical insights into cHCC-ICC.