Common genetic variation in alcohol-related hepatocellular carcinoma: a case-control genome-wide association study.

BACKGROUND: Hepatocellular carcinoma is a frequent consequence of alcohol-related liver disease, with variable incidence among heavy drinkers. We did a genome-wide association study (GWAS) to identify common genetic variants for alcohol-related hepatocellular carcinoma. METHODS: We conducted a two-stage case-control GWAS in a discovery cohort of 2107 unrelated European patients with alcohol-related liver disease aged 20-92 years recruited between Oct 22, 1993, and March 12, 2017. Cases were patients with alcohol-related hepatocellular carcinoma diagnosed by imaging or histology. Controls were patients with alcohol-related liver disease without hepatocellular carcinoma. We used an additive logistic regression model adjusted for the first ten principal components to assess genetic variants associated with alcohol-related hepatocellular carcinoma. We did another analysis with adjustment for age, sex, and liver fibrosis. New candidate associations (p<1 × 10-6) and variants previously associated with alcohol-related hepatocellular carcinoma were evaluated in a validation cohort of 1933 patients with alcohol-related liver disease aged 29-92 years recruited between July 21, 1995, and May 2, 2019. We did a meta-analysis of the two case-control cohorts. FINDINGS: The discovery cohort included 775 cases and 1332 controls. Of 7 962 325 variants assessed, we identified WNT3A-WNT9A (rs708113; p=1·11 × 10-8) and found support for previously reported regions associated with alcohol-related hepatocellular carcinoma risk at TM6SF2 (rs58542926; p=6·02 × 10-10), PNPLA3 (rs738409; p=9·29 × 10-7), and HSD17B13 (rs72613567; p=2·49 × 10-4). The validation cohort included 874 cases and 1059 controls and three variants were replicated: WNT3A-WNT9A (rs708113; p=1·17 × 10-3), TM6SF2 (rs58542926; p=4·06 × 10-5), and PNPLA3 (rs738409; p=1·17 × 10-4). All three variants reached GWAS significance in the meta-analysis: WNT3A-WNT9A (odds ratio 0·73, 95% CI 0·66-0·81; p=3·93 × 10-10), TM6SF2 (1·77, 1·52-2·07; p=3·84×10-13), PNPLA3 (1·34, 1·22-1·47; p=7·30 × 10-10). Adjustment for clinical covariates yielded similar results. We observed an additive effect of at-risk alleles on alcohol-related hepatocellular carcinoma. WNT3A-WNT9A rs708113 was not associated with liver fibrosis. INTERPRETATION: WNT3A-WNT9A is a susceptibility locus for alcohol-related hepatocellular carcinoma, suggesting an early role of the Wnt-ß-catenin pathway in alcohol-related hepatocellular carcinoma carcinogenesis. FUNDING: Ligue Nationale contre le Cancer, Bpifrance, INSERM, AFEF, CARPEM, Labex OncoImmunology, and Agence Nationale de la Recherche.

Bi-allelic hydroxymethylbilane synthase inactivation defines a homogenous clinico-molecular subtype of hepatocellular carcinoma.

BACKGROUND & AIMS: Acute intermittent porphyria (AIP), caused by heterozygous germline mutations of the heme synthesis pathway enzyme HMBS (hydroxymethylbilane synthase), confers a high risk of hepatocellular carcinoma (HCC) development. Yet, the role of HMBS in liver tumorigenesis remains unclear. METHODS: Herein, we explore HMBS alterations in a large series of 758 HCC cases, including 4 patients with AIP. We quantify the impact of HMBS mutations on heme biosynthesis pathway intermediates and we investigate the molecular and clinical features of HMBS-mutated tumors. RESULTS: We identify recurrent bi-allelic HMBS inactivation, both in patients with AIP acquiring a second somatic HMBS mutation and in sporadic HCC with 2 somatic hits. HMBS alterations are enriched in truncating mutations, in particular in splice regions, leading to abnormal transcript structures. Bi-allelic HMBS inactivation results in a massive accumulation of its toxic substrate porphobilinogen and synergizes with CTNNB1-activating mutations, leading to the development of well-differentiated tumors with a transcriptomic signature of Wnt/ß-catenin pathway activation and a DNA methylation signature related to ageing. HMBS-inactivated HCC mostly affects females, in the absence of fibrosis and classical HCC risk factors. CONCLUSIONS: These data identify HMBS as a tumor suppressor gene whose bi-allelic inactivation defines a homogenous clinical and molecular HCC subtype. LAY SUMMARY: Heme (the precursor to hemoglobin, which plays a key role in oxygen transport around the body) synthesis occurs in the liver and involves several enzymes including hydroxymethylbilane synthase (HMBS). HMBS mutations cause acute intermittent porphyria, a disease caused by the accumulation of toxic porphyrin precursors. Herein, we show that HMBS inactivation is also involved in the development of liver cancers with distinct clinical and molecular characteristics.

DNA Methylation Signatures Reveal the Diversity of Processes Remodeling Hepatocellular Carcinoma Methylomes.

BACKGROUND AND AIMS: DNA methylation patterns are highly rearranged in HCCs. However, diverse sources of variation are intermingled in cancer methylomes, precluding the precise characterization of underlying molecular mechanisms. We developed a computational framework (methylation signature analysis with independent component analysis [MethICA]) leveraging independent component analysis to disentangle the diverse processes contributing to DNA methylation changes in tumors. APPROACH AND RESULTS: Applied to a collection of 738 HCCs, MethICA unraveled 13 stable methylation components preferentially active in specific chromatin states, sequence contexts, and replication timings. These included signatures of general processes associated with sex and age but also signatures related to specific driver events and molecular subgroups. Catenin beta 1 mutations were major modulators of methylation patterns in HCC, characterized by a targeted hypomethylation of transcription factor 7-bound enhancers in the vicinity of Wnt target genes as well as a widespread hypomethylation of late-replicated partially methylated domains. By contrast, demethylation of early replicated highly methylated domains was a signature of replication stress, leading to an extensive hypomethylator phenotype in cyclin-activated HCC. Inactivating mutations of the chromatin remodeler AT-rich interactive domain-containing protein 1A were associated with epigenetic silencing of differentiation-promoting transcriptional networks, also detectable in cirrhotic liver. Finally, a hypermethylation signature targeting polycomb-repressed chromatin domains was identified in the G1 molecular subgroup with progenitor features. CONCLUSIONS: This study elucidates the diversity of processes remodeling HCC methylomes and reveals the epigenetic and transcriptional impact of driver alterations.

Telomere length is key to hepatocellular carcinoma diversity and telomerase addiction is an actionable therapeutic target.

BACKGROUND & AIMS: Telomerase activation is the earliest event in hepatocellular carcinoma (HCC) development. Thus, we aimed to elucidate the role of telomere length maintenance during liver carcinogenesis. METHODS: Telomere length was measured in the tumor and non-tumor liver tissues of 1,502 patients (978 with HCC) and integrated with TERT alterations and expression, as well as clinical and molecular (analyzed by genome, exome, targeted and/or RNA-sequencing) features of HCC. The preclinical efficacy of anti-TERT antisense oligonucleotides (ASO) was assessed in vitro in 26 cell lines and in vivo in a xenograft mouse model. RESULTS: Aging, liver fibrosis, male sex and excessive alcohol consumption were independent determinants of liver telomere attrition. HCC that developed in livers with long telomeres frequently had wild-type TERT with progenitor features and BAP1 mutations. In contrast, HCC that developed on livers with short telomeres were enriched in the non-proliferative HCC class and frequently had somatic TERT promoter mutations. In HCCs, telomere length is stabilized in a narrow biological range around 5.7 kb, similar to non-tumor livers, by various mechanisms that activate TERT expression. Long telomeres are characteristic of very aggressive HCCs, associated with the G3 transcriptomic subclass, TP53 alterations and poor prognosis. In HCC cell lines, TERT silencing with ASO was efficient in highly proliferative and poorly differentiated cells. Treatment for 3 to 16 weeks induced cell proliferation arrest in 12 cell lines through telomere shortening, DNA damage and activation of apoptosis. The therapeutic effect was also obtained in a xenograft mouse model. CONCLUSIONS: Telomere maintenance in HCC carcinogenesis is diverse, and is associated with tumor progression and aggressiveness. The efficacy of anti-TERT ASO treatment in cell lines revealed the oncogenic addiction to TERT in HCC, providing a preclinical rationale for anti-TERT ASO treatment in HCC clinical trials. LAY SUMMARY: Telomeres are repeated DNA sequences that protect chromosomes and naturally shorten in most adult cells because of the inactivation of the TERT gene, coding for the telomerase enzyme. Here we show that telomere attrition in the liver, modulated by aging, sex, fibrosis and alcohol, associates with specific clinical and molecular features of hepatocellular carcinoma, the most frequent primary liver cancer. We also show that liver cancer is dependent on TERT reactivation and telomere maintenance, which could be targeted through a novel therapeutic approach called antisense oligonucleotides.

Clinical Impact of Genomic Diversity From Early to Advanced Hepatocellular Carcinoma.

To date, genomic analyses of hepatocellular carcinoma (HCC) have been limited to early stages obtained from liver resection. We aim to describe the genomic profiling of HCC from early to advanced stages. We analyzed 801 HCC from 720 patients (410 resections, 137 transplantations, 122 percutaneous ablations, and 52 noncurative) for 190 gene expressions and for 31 gene mutations. Forty-one advanced HCC and 156 whole exome of Barcelona Clinic Liver Cancer (BCLC) 0/A were analyzed by whole-exome sequencing. Genomic profiling was correlated with tumor stages, clinical features, and survival. Our cohort included patients classified in BCLC stage 0 (9.4%), A (59.5%), B (16.2%), and C (14.9%). Among the overall 801 HCC, the most frequently mutated genes were telomerase reverse transcriptase (TERT) (58.1%), catenin beta 1 (CTNNB1) (30.7%), tumor protein 53 (TP53; 18.7%), AT-rich interaction domain 1A (ARID1A) (13%), albumin (11.4%), apolipoprotein B (APOB) (9.4%), and AXIN1 (9.2%). Advanced-stage HCC (BCLC B/C) showed higher frequencies of splicing factor 3b subunit 1 (SF3B1) (P = 0.0003), TP53 (P = 0.0006), and RB Transcriptional Corepressor 1 mutations (P = 0.03). G1-G6 transcriptomic classification and the molecular prognostic 5-gene score showed different distributions according to the stage of the disease and the type of treatment with an enrichment of G3 (P < 0.0001), poor prognostic score (P < 0.0001), and increased proliferation and dedifferentiation at the transcriptomic level in advanced HCC. The 5-gene score predicted survival in patients treated by resection (P < 0.0001) and ablation (P = 0.01) and in advanced HCC (P = 0.04). Twenty-two percent of advanced HCC harbored potentially druggable genetic alterations, and MET amplification was associated with complete tumor response in patients with advanced HCC treated by a specific MET inhibitor. Conclusion: Genomic analysis across the different stages of HCC revealed the mechanisms of tumor progression and helped to identify biomarkers of response to targeted therapies.

Adeno-associated virus in the liver: natural history and consequences in tumour development.

OBJECTIVE: Adeno-associated virus (AAV) is a defective mono-stranded DNA virus, endemic in human population (35%-80%). Recurrent clonal AAV2 insertions are associated with the pathogenesis of rare human hepatocellular carcinoma (HCC) developed on normal liver. This study aimed to characterise the natural history of AAV infection in the liver and its consequence in tumour development. DESIGN: Viral DNA was quantified in tumour and non-tumour liver tissues of 1461 patients. Presence of episomal form and viral mRNA expression were analysed using a DNAse/TaqMan-based assay and quantitative RT-PCR. In silico analyses using viral capture data explored viral variants and new clonal insertions. RESULTS: AAV DNA was detected in 21% of the patients, including 8% of the tumour tissues, equally distributed in two major viral subtypes: one similar to AAV2, the other hybrid between AAV2 and AAV13 sequences. Episomal viral forms were found in 4% of the non-tumour tissues, frequently associated with viral RNA expression and human herpesvirus type 6, the candidate natural AAV helper virus. In 30 HCC, clonal AAV insertions were recurrently identified in CCNA2, CCNE1, TERT, TNFSF10, KMT2B and GLI1/INHBE. AAV insertion triggered oncogenic overexpression through multiple mechanisms that differ according to the localisation of the integration site. CONCLUSION: We provided an integrated analysis of the wild-type AAV infection in the liver with the identification of viral genotypes, molecular forms, helper virus relationship and viral integrations. Clonal AAV insertions were positive selected during HCC development on non-cirrhotic liver challenging the notion of AAV as a non-pathogenic virus.

Recurrent chromosomal rearrangements of ROS1, FRK and IL6 activating JAK/STAT pathway in inflammatory hepatocellular adenomas.

BACKGROUND: Inflammatory hepatocellular adenomas (IHCAs) are benign liver tumours characterised by an activation of the janus kinase (JAK)/signal transducers and activators of transcription (STAT) pathway caused by oncogenic activating mutations. However, a subset of IHCA lacks of identified mutation explaining the inflammatory phenotype. METHODS: 657 hepatocellular adenomas developed in 504 patients were analysed for gene expression of 17 genes and for mutations in seven genes by sequencing. 22 non-mutated IHCAs were analysed by whole-exome and/or RNA sequencing. RESULTS: We identified 296 IHCA (45%), 81% of them were mutated in either IL6ST (61%), FRK (8%), STAT3 (5%), GNAS (3%) or JAK1 (2%). Among non-mutated IHCA, RNA sequencing identified recurrent chromosome rearrangement involving ROS1, FRK or IL6 genes. ROS1 fusions were identified in 8 IHCA, involving C-terminal part of genes highly expressed in the liver (PLG, RBP4, APOB) fused with exon 33-35 to 43 of ROS1 including the tyrosine kinase domain. In two cases a truncated ROS1 transcript from exon 36 to 43 was identified. ROS1 rearrangements were validated by fluorescence in situ hybridisation (FISH) and led to ROS1 overexpression. Among the 5 IHCA with FRK rearrangements, 5 different partners were identified (MIA3, MIA2, LMO7, PLEKHA5, SEC16B) fused to a common region in FRK that included exon 3-8. No overexpression of FRK transcript was detected but the predicted chimeric proteins lacked the auto-inhibitory SH2-SH3 domains. In two IHCA, we identified truncated 3'UTR of IL6 associated with overexpression of the transcript. CONCLUSION: Recurrent chromosomal alterations involving ROS1, FRK or IL6 genes lead to activation of the JAK/STAT pathway in IHCAs.

BAP1 mutations define a homogeneous subgroup of hepatocellular carcinoma with fibrolamellar-like features and activated PKA.

BACKGROUND & AIMS: DNAJB1-PRKACA fusion is a specific driver event in fibrolamellar carcinoma (FLC), a rare subtype of hepatocellular carcinoma (HCC) that occurs in adolescents and young adults. In older patients, molecular determinants of HCC with mixed histological features of HCC and FLC (mixed-FLC/HCC) remain to be discovered. METHODS: A series of 151 liver tumors including 126 HCC, 15 FLC, and 10 mixed-FLC/HCC were analyzed by RNAseq and whole-genome- or whole-exome sequencing. Western blots were performed to validate genomic discoveries. Results were validated using the TCGA database. RESULTS: Most of the mixed-FLC/HCC RNAseq clustered in a robust subgroup of 17 tumors, which all had mutations or translocations inactivating BAP1, the gene encoding BRCA1-associated protein-1. Like FLC, BAP1-HCC were significantly enriched in females, patients with a lack of chronic liver disease, and fibrotic tumors compared to non-BAP1 HCC. However, patients were older and had a poorer prognosis than those with FLC. BAP1 tumors were immune hot, showed progenitor features and did not show DNAJB1-PRKACA fusion, while almost none of these tumors had mutations in CTNNB1, TP53 and TERT promoter. In contrast, 80% of the BAP1 tumors showed a chromosome gain of PRKACA at 19p13, combined with a loss of PRKAR2A (coding for the inhibitory regulatory subunit of PKA) at 3p21, leading to a high PRKACA/PRKAR2A ratio at the mRNA and protein levels. CONCLUSION: We have characterized a subgroup of BAP1-driven HCC with fibrolamellar-like features and a dysregulation of the PKA pathway, which could be at the root of the clinical and histological similarities between BAP1 tumors and DNAJB1-PRKACA FLCs. LAY SUMMARY: Herein, we have defined a homogeneous subgroup of hepatocellular carcinomas in which the BAP1 gene is inactivated. This leads to the development of cancers with features similar to those of fibrolamellar carcinoma. These tumors more frequently develop in females without chronic liver disease or cirrhosis. The presence of PKA activation and T cell infiltrates suggest that these tumors could be treated with PKA inhibitors or immunomodulators.

Analysis of Liver Cancer Cell Lines Identifies Agents With Likely Efficacy Against Hepatocellular Carcinoma and Markers of Response.

BACKGROUND AND AIMS: Hepatocellular carcinomas (HCCs) are heterogeneous aggressive tumors with low rates of response to treatment at advanced stages. We screened a large panel of liver cancer cell lines (LCCLs) to identify agents that might be effective against HCC and markers of therapeutic response. METHODS: We performed whole-exome RNA and microRNA sequencing and quantification of 126 proteins in 34 LCCLs. We screened 31 anticancer agents for their ability to decrease cell viability. We compared genetic, RNA, and protein profiles of LCCLs with those of primary HCC samples and searched for markers of response. RESULTS: The protein, RNA and mutational signatures of the LCCLs were similar to those of the proliferation class of HCC, which is the most aggressive tumor type. Cell lines with alterations in genes encoding members of the Ras-MAPK signaling pathway and that required fibroblast growth factor (FGF)19 signaling via FGF receptor 4 for survival were more sensitive to trametinib than to FGF receptor 4 inhibitors. Amplification of FGF19 resulted in increased activity of FGF19 only in tumor cells that kept a gene expression pattern of hepatocyte differentiation. We identified single agents and combinations of agents that reduced viability of cells with features of the progenitor subclass of HCC. LCCLs with inactivating mutations in TSC1 and TSC2 were sensitive to the mammalian target of rapamycin inhibitor rapamycin, and cells with inactivating mutations in TP53 were sensitive to the Aurora kinase A inhibitor alisertib. Amplification of MET was associated with hypersensitivity to cabozantinib and the combination of sorafenib and inhibitors of MAP kinase 1 and MAP kinase2 had a synergistic antiproliferative effect. CONCLUSION: LCCLs can be screened for drugs and agents that might be effective for treatment of HCC. We identified genetic alterations and gene expression patterns associated with response to these agents. This information might be used to select patients for clinical trials.

A 17-Beta-Hydroxysteroid Dehydrogenase 13 Variant Protects From Hepatocellular Carcinoma Development in Alcoholic Liver Disease.

Recently, a loss of function variant (rs72613567) in 17-beta-hydroxysteroid dehydrogenase 13 (HSD17B13) has been identified as protective of nonalcoholic (NAFLD) and alcoholic liver disease (ALD). However, the role of this single-nucleotide polymorphism (SNP) in the development of hepatocellular carcinoma (HCC) is currently unknown. A total of 3,315 European patients with HCC (n = 1,109) or without HCC, but with chronic liver disease (CLD; n = 2,206), from four centers were analyzed either by whole-exome sequencing (WES; exploratory cohort, 285 HCC) or genotyped for HSD17B13 rs72613567 (validation cohort, 824 HCC and all CLD). We included a control group of 33,337 healthy European individuals from the Exome Aggregation Consortium. We compared distribution of genotype using the chi-square test and logistic regression. In the exploratory cohort analyzed by WES, frequency of the TA allele of HSD17B13 rs72613567 was significantly decreased in HCC patients compared to healthy controls (P = 1.52 × 10-06 ). In the validation cohort, frequency of TA allele carriers was also decreased in patients with CLD and without HCC (39%) compared to healthy individuals (47%; P < 0.0001). The protective effect of the TA allele of HSD17B13 rs72613567 was identified in patients with ALD (odds ratio [OR] = 0.73; 95% confidence interval [CI], 0.65-0.82; P < 0.0001), NAFLD (OR = 0.64; 95% CI, 0.49-0.83; P = 0.0007), and hepatitis C (OR = 0.71; 95% CI, 0.60-0.85; P = 0.0002). In patients with ALD, the proportion of TA allele carriers with HCC was significantly lower (32%) than in CLD patients without HCC (40%), even after adjustment for age, sex, and fibrosis (OR = 0.64; 95% CI, 0.46-0.87; P = 0.005). Conclusion: The HSD17B13 rs72613567 loss of function variant is protective of HCC development in patients with ALD.

Palimpsest: an R package for studying mutational and structural variant signatures along clonal evolution in cancer.

Summary: Cancer genomes are altered by various mutational processes and, like palimpsests, bear the signatures of these different processes. The Palimpsest R package provides a complete workflow for the characterization and visualization of mutational signatures and their evolution along tumor development. The package covers a wide range of functions for extracting both base substitution and structural variant signatures, inferring the clonality of each alteration and analyzing the evolution of mutational processes between early clonal and late subclonal events. Palimpsest also estimates the probability of each mutation being due to each process to predict the mechanisms at the origin of driver events. Palimpsest is an easy-to-use toolset for reconstructing the natural history of a tumor using whole exome or whole genome sequencing data. Availability and implementation: Palimpsest is freely available at www.github.com/FunGEST/Palimpsest. Supplementary information: Supplementary data are available at Bioinformatics online.

High Impact: The Role of Promiscuous Binding Sites in Polypharmacology.

The literature focuses on drug promiscuity, which is a drug's ability to bind to several targets, because it plays an essential role in polypharmacology. However, little work has been completed regarding binding site promiscuity, even though its properties are now recognized among the key factors that impact drug promiscuity. Here, we quantified and characterized the promiscuity of druggable binding sites from protein-ligand complexes in the high quality Mother Of All Databases while using statistical methods. Most of the sites (80%) exhibited promiscuity, irrespective of the protein class. Nearly half were highly promiscuous and able to interact with various types of ligands. The corresponding pockets were rather large and hydrophobic, with high sulfur atom and aliphatic residue frequencies, but few side chain atoms. Consequently, their interacting ligands can be large, rigid, and weakly hydrophilic. The selective sites that interacted with one ligand type presented less favorable pocket properties for establishing ligand contacts. Thus, their ligands were highly adaptable, small, and hydrophilic. In the dataset, the promiscuity of the site rather than the drug mainly explains the multiple interactions between the drug and target, as most ligand types are dedicated to one site. This underlines the essential contribution of binding site promiscuity to drug promiscuity between different protein classes.

Deep dive into the immune response against murine mesothelioma permits design of novel anti-mesothelioma therapeutics.

Given the need to improve the efficacy of standard-of-care immunotherapy (anti-CTLA-4 + anti-PD-1) in human malignant pleural mesothelioma (hMPM), we thoroughly characterized the immunobiology of the AB12 murine mesothelioma (MM) model, aiming to increase its accuracy in predicting the response of hMPM to immunotherapy and in designing novel anti-hMPM treatments. Specifically, we used immunologic, transcriptomic and survival analyses, to synchronize the MM tumor growth phases and immune evolution with the histo-molecular and immunological characteristics of hMPM while also determining the anti-MM efficacy of standard-of-care anti-hMPM immunotherapy as a benchmark that novel therapeutics should meet. We report that early-, intermediate- and advanced- AB12 tumors are characterized by a bell-shaped anti-tumor response that peaks in intermediate tumors and decays in advanced tumors. We further show that intermediate- and advanced- tumors match with immune active ("hot") and immune inactive ("cold") hMPM respectively, and that they respond to immunotherapy in a manner that corresponds well with its performance in real-life settings. Finally, we show that in advanced tumors, addition of cisplatin to anti CTLA-4 + anti PD-1 can extend mice survival and invigorate the decaying anti-tumor response. Therefore, we highlight this triple combination as a worthy candidate to improve clinical outcomes in hMPM.

Structure, Dynamics, and Impact of Replication Stress-Induced Structural Variants in Hepatocellular Carcinoma.

Oncogene activation leads to replication stress and promotes genomic instability. Here we combine optical mapping and whole-genome sequencing (WGS) to explore in depth the nature of structural variants (SV) induced by replication stress in cyclin-activated hepatocellular carcinomas (CCN-HCC). In addition to classical tandem duplications, CCN-HCC displayed frequent intra-chromosomal and interchromosomal templated insertion cycles (TIC), likely resulting from template switching events. Template switching preferentially involves active topologically associated domains that are proximal to one another within the 3D genome. Template sizes depend on the type of cyclin activation and are coordinated within each TIC. Replication stress induced continuous accumulation of SVs during CCN-HCC progression, fostering the acquisition of new driver alterations and large-scale copy-number changes at TIC borders. Together, this analysis sheds light on the mechanisms, dynamics, and consequences of SV accumulation in tumors with oncogene-induced replication stress. SIGNIFICANCE: Optical mapping and whole-genome sequencing integration unravels a unique signature of replication stress-induced structural variants that drive genomic evolution and the acquisition of driver events in CCN-HCC.

Integrated Genomic Analysis Identifies Driver Genes and Cisplatin-Resistant Progenitor Phenotype in Pediatric Liver Cancer.

Pediatric liver cancers (PLC) comprise diverse diseases affecting infants, children, and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of patients with hepatoblastoma, all before three years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between "hepatocytic," "liver progenitor," and "mesenchymal" molecular subgroups of hepatoblastoma. We showed that during chemotherapy, "liver progenitor" cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies. SIGNIFICANCE: PLCs are deadly when they resist chemotherapy, with limited alternative treatment options. Using a multiomics approach, we identified PLC driver genes and the cellular phenotype at the origin of cisplatin resistance. We validated new treatments targeting these molecular features in cell lines and xenografts.This article is highlighted in the In This Issue feature, p. 2355.

Multi-site tumor sampling highlights molecular intra-tumor heterogeneity in malignant pleural mesothelioma.

BACKGROUND: Malignant pleural mesothelioma (MPM) is a heterogeneous cancer. Better knowledge of molecular and cellular intra-tumor heterogeneity throughout the thoracic cavity is required to develop efficient therapies. This study focuses on molecular intra-tumor heterogeneity using the largest series to date in MPM and is the first to report on the multi-omics profiling of a substantial series of multi-site tumor samples. METHODS: Intra-tumor heterogeneity was investigated in 16 patients from whom biopsies were taken at distinct anatomical sites. The paired biopsies collected from apex, side wall, costo-diaphragmatic, or highest metabolic sites as well as 5 derived cell lines were screened using targeted sequencing. Whole exome sequencing, RNA sequencing, and DNA methylation were performed on a subset of the cohort for deep characterization. Molecular classification, recently defined histo-molecular gradients, and cell populations of the tumor microenvironment were assessed. RESULTS: Sequencing analysis identified heterogeneous variants notably in NF2, a key tumor suppressor gene of mesothelial carcinogenesis. Subclonal tumor populations were shared among paired biopsies, suggesting a polyclonal dissemination of the tumor. Transcriptome analysis highlighted dysregulation of cell adhesion and extracellular matrix pathways, linked to changes in histo-molecular gradient proportions between anatomic sites. Methylome analysis revealed the contribution of epigenetic mechanisms in two patients. Finally, significant changes in the expression of immune mediators and genes related to immunological synapse, as well as differential infiltration of immune populations in the tumor environment, were observed and led to a switch from a hot to a cold immune profile in three patients. CONCLUSIONS: This comprehensive analysis reveals patient-dependent spatial intra-tumor heterogeneity at the genetic, transcriptomic, and epigenetic levels and in the immune landscape of the tumor microenvironment. Results support the need for multi-sampling for the implementation of molecular-based precision medicine.