Genomic evolution and the impact of SLIT2 mutation in relapsed intrahepatic cholangiocarcinoma.

BACKGROUND AND AIMS: Intrahepatic cholangiocarcinoma (ICC) is aggressive and has high rates of relapse, conferring poor long-term survival after curative resection. Little is known about the genomic evolution that occurs during ICC relapse. APPROACH AND RESULTS: We conducted whole-exome sequencing of 30 paired primary and relapsed tumors from 10 patients with ICC who received curative resection. We sought to identify frequently altered genes, infer tumor subclonal architectures, and track genomic evolution from primary to relapsed tumors. We examined functional effects and the mechanism of action of SLIT2, a gene specifically mutated in relapsed tumors, on tumor growth and metastasis and the tumor microenvironment. Our results indicated that relapsed ICCs were genetically derived from intrahepatic dissemination of primary tumors. However, they acquired additional mutations while maintaining most drivers, such as TP53 and IDH1. Multiregion sequencing suggested polyclonal seeding of ICC dissemination. Four of 10 relapsed ICCs acquired SLIT2 mutations that were not present in the corresponding primary tumors. Validation in an expanded sample revealed SLIT2 mutations in 2.3% (1/44) of primary ICCs and 29.5% (13/44) of relapsed ICCs. Biofunctional investigations revealed that inactivating mutation of SLIT2 resulted in activation of PI3K-Akt signaling in ICC cells, directly enhanced neutrophil chemotaxis, mediated tumor-associated neutrophil infiltration, and contributed to ICC growth and metastasis. CONCLUSIONS: We characterized genomic evolution during ICC relapse and identified SLIT2 as a driver of tumor dissemination and tumor-associated neutrophil infiltration.

Tumor-derived exosomes induce immunosuppressive macrophages to foster intrahepatic cholangiocarcinoma progression.

BACKGROUND AND AIMS: Macrophages are prominent components of solid tumors and exhibit distinct functions in different tumor microenvironments. Exosomes are emerging as necessary mediators of the cross-talk between tumor cells and the microenvironment. However, the underlying mechanisms of exosomes involving into crosstalk between tumor cells and macrophages during disease progression of intrahepatic cholangiocarcinoma (ICC) have not been yet fully realized. APPROACH AND RESULTS: We found that the macrophages of ICC tumor tissues up-regulated the expression levels of immunosuppressive molecule programmed death-ligand 1 (PD-L1). Increased PD-L1+ macrophages in tumor tissues effectively suppressed T-cell immunity and correlated with poor survival rates in patients with ICC. High-throughput RNA-sequencing analysis that was performed to identify differential levels of microRNAs (miRNAs) between exosomes derived from ICC cells and primary human intrahepatic biliary epithelial cells revealed that miR-183-5p was increased in ICC cell-derived exosomes. Exosomal miR-183-5p inhibited phosphatase and tensin homolog (PTEN) expression, to subsequently affect the elevations on both phosphorylated AKT and PD-L1 expression in macrophages. Furthermore, macrophages that treated with ICC cell-derived exosomes significantly suppressed T-cell immunity in vitro and contributed to the growth and progression of ICC in vivo, which were reversible through blockages on PD-L1 of these macrophages. Finally, clinical data showed that up-regulated levels of plasma exosomal miR-183-5p correlated with poor prognosis of patients with ICC after curative resection. CONCLUSIONS: Tumor-derived exosomal miR-183-5p up-regulates PD-L1-expressing macrophages to foster immune suppression and disease progression in ICC through the miR-183-5p/PTEN/AKT/PD-L1 pathway. Exosomal miR-183-5p is a potential predictive biomarker for ICC progression and a potential target for development of therapeutic strategies against immune tolerance feature of ICC.

Long noncoding RNA lnc-H2AFV-1 promotes cell growth by regulating aberrant m6A RNA modification in head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is the most common malignant tumor in the oral and maxillofacial regions, and long noncoding RNAs (lncRNAs) play crucial roles in the occurrence and progression of HNSCC. The lncRNA lnc-H2AFV-1 was found to be upregulated in HNSCC tissues; however, the function of lnc-H2AFV-1 in regulating HNSCC proliferation and the potential molecular mechanism is unclear. The present study evaluated the expression of lnc-H2AFV-1 in HNSCC tissues using quantitative real-time PCR (qPCR) and associated abundant lnc-H2AFV-1 expression with tumor size. Functionally, lnc-H2AFV-1 significantly promoted the proliferation of HNSCC cells in vitro and in vivo. Quantified N6-methyladenosine (m6A) RNA methylation and dot blot assays revealed that total m6A methylation in HNSCC cells was accompanied by lnc-H2AFV-1 expression. Western blotting showed that the expression of methyltransferase-like (METTL) 3 and METTL14 was consistent with that of lnc-H2AFV-1, whereas the expression of demethylase fat mass and obesity-associated (FTO) was contrary to that of lnc-H2AFV-1. Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and MeRIP-qPCR revealed that lnc-H2AFV-1 overexpression led to the elevated expression and maximal m6A methylation of intraflagellar transport (IFT) 80 in HNSCC. In addition, METTL3/14 knockdown decreased IFT80 expression. Thus, our findings suggested that lnc-H2AFV-1 might be a biomarker that alters m6A modification by regulating the m6A methylases METTL3/14 and FTO and then mediating the downstream target IFT80 to promote HNSCC progression.

WDR5 facilitates EMT and metastasis of CCA by increasing HIF-1α accumulation in Myc-dependent and independent pathways.

Cholangiocarcinoma (CCA) is a highly aggressive malignancy with extremely poor prognoses. The oncogenic role and prognostic value of c-Myc in CCA is not well elucidated. WD repeat domain 5 (WDR5) is a critical regulatory factor directly interacting with c-Myc to regulate c-Myc recruitment at chromosomal locations, but the interaction of WDR5 and c-Myc in CCA was uncovered. In our study, we detected WDR5 and c-Myc expression in all CCA types, including intrahepatic (iCCA), perihilar (pCCA), and distal (dCCA) CCA, and evaluated their prognostic significance. Consequently, we demonstrated that WDR5 was significantly correlated with poor prognosis of CCA and that WDR5 and c-Myc co-expression was a more sensitive prognostic factor. With in vitro and in vivo experiments and bioinformatics, we showed that WDR5 interacted with the Myc box IIIb (MBIIIb) motif of c-Myc and facilitated Myc-induced HIF1A transcription, thereby promoting the epithelial-mesenchymal transition (EMT), invasion, and metastasis of CCA. Moreover, WDR5 enhanced hypoxia-inducible factor 1 subunit α (HIF-1α) accumulation by binding with histone deacetylase 2 (HDAC2) and increasing histone 3 lysine 4 acetylation (H3K4ac) deacetylation of the prolyl hydroxylase domain protein 2 (PHD2) promoter, resulting in the attenuation of chromatin opening and PHD2 expression, and eventually leading to HIF-1α stabilization and accumulation. In conclusion, WDR5 facilitated EMT and metastasis of CCA by increasing HIF-1α accumulation in a Myc-dependent pathway to promote HIF-1α transcription and a Myc-independent pathway to stabilize HIF-1α.

Characterization of tumor microbiome and associations with prognosis in intrahepatic cholangiocarcinoma.

BACKGROUND: The tumor microbiome has been characterized in several malignancies; however, no previous studies have investigated its role in intrahepatic cholangiocarcinoma (ICC). Hence, we explored the tumor microbiome and its association with prognosis in ICC. METHODS: One hundred and twenty-one ICC tumor samples and 89 adjacent normal tissues were profiled by 16S rRNA sequencing. Microbial differences between tumor and adjacent nontumoral liver tissues were assessed. Tumor microbial composition was then evaluated to detect its association with prognosis. Finally, a risk score calculated by the tumor microbiota was accessed by the least absolute shrinkage and selector operator method (Lasso) to predict prognosis of ICC. RESULTS: The tumor microbiome displayed a greater diversity than that in adjacent nontumoral liver tissues. Tumor samples were characterized by a higher abundance of Firmicutes, Actinobacteria, Bacteroidetes, and Acidobacteriota. Higher tumor microbial α diversity was associated with lymph node metastasis and predicted shortened overall survival (OS) and recurrence-free survival (RFS). A total of 11 bacteria were selected to generate the risk score by Lasso. This score showed potential in predicting OS, and was an independent risk factor for OS. CONCLUSION: In conclusion, our study characterized the tumor microbiome and revealed its role in predicting prognosis in ICC.

Development and Validation of a Stromal-Immune Signature to Predict Prognosis in Intrahepatic Cholangiocarcinoma.

Background: Intrahepatic cholangiocarcinoma (ICC) is a highly desmoplastic tumor with poor prognosis even after curative resection. We investigated the associations between the composition of the ICC stroma and immune cell infiltration and aimed to develop a stromal-immune signature to predict prognosis in surgically treated ICC. Patients and methods: We recruited 359 ICC patients and performed immunohistochemistry to detect α-smooth muscle actin (α-SMA), CD3, CD4, CD8, Foxp3, CD68, and CD66b. Aniline was used to stain collagen deposition. Survival analyses were performed to detect prognostic values of these markers. Recursive partitioning for a discrete-time survival tree was applied to define a stromal-immune signature with distinct prognostic value. We delineated an integrated stromal-immune signature based on immune cell subpopulations and stromal composition to distinguish subgroups with different recurrence-free survival (RFS) and overall survival (OS) time. Results: We defined four major patterns of ICC stroma composition according to the distributions of α-SMA and collagen: dormant (α-SMAlow/collagenhigh), fibrogenic (α-SMAhigh/collagenhigh), inert (α-SMAlow/collagenlow), and fibrolytic (α-SMAhigh/collagenlow). The stroma types were characterized by distinct patterns of infiltration by immune cells. We divided patients into six classes. Class I, characterized by high CD8 expression and dormant stroma, displayed the longest RFS and OS, whereas Class VI, characterized by low CD8 expression and high CD66b expression, displayed the shortest RFS and OS. The integrated stromal-immune signature was consolidated in a validation cohort. Conclusion: We developed and validated a stromal-immune signature to predict prognosis in surgically treated ICC. These findings provide new insights into the stromal-immune response to ICC.

Transcriptome Mapping of the Internal N7-Methylguanosine Methylome in Messenger RNAs in Human Oral Squamous Cell Carcinoma.

BACKGROUND: Internal N7-methylguanosine (m7G) methylation in mammalian messenger RNAs (mRNAs) is essential in disease development. However, the status of internally m7G-modified mRNAs in oral squamous cell carcinoma (OSCC) remains poorly understood. METHODS: Methylated RNA immunoprecipitation sequencing (MeRIP-seq) was used to identify the m7G modification level of mRNAs and the expression of mRNAs between OSCC and normal tissues. These differentially methylated and expressed genes were subjected to Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and construction of protein-protein interaction (PPI) networks. Quantitative real-time PCR (qPCR) assay was performed to detect the expression of Baculoviral IAP Repeat Containing 3 (BIRC3) in vitro. The biological function of BIRC3 in OSCC was clarified using CCK-8, Transwell migration and Western blot assays. RESULTS: The m7G-mRNA profile showed 9514 unique m7G peaks within 7455 genes in OSCC tissues. In addition, the most conserved m7G motif within mRNAs in OSCC was GGARG (R = G/A). The identified m7G peaks were mainly distributed in the coding sequence region within mRNAs in OSCC. GO enrichment and KEGG pathway analyses showed that m7G-modified genes were closely related to cancer progression. m7G-modified hub genes were screened from the constructed PPI networks. Furthermore, BIRC3 with high m7G methylation showed high expression in OSCC cell lines, as confirmed by qPCR assay. Functionally, the knockdown of BIRC3 significantly inhibited the proliferation and migration ability of CAL-27 cells in vitro functional assays. In addition, the relative expression of E-cadherin expression was elevated, while Vimentin and N-cadherin protein expression was decreased in CAL-27 cells transfected with si-BIRC3. This study suggests that BIRC3 could promote OSCC proliferation and migration, which may be associated with involvement in epithelial-mesenchymal transition (EMT) progression. CONCLUSIONS: This paper constructed a transcriptome map of internal m7G in mRNAs, which provides potential research value to study the role of m7G methylation in OSCC.

Association of BRAF Variants With Disease Characteristics, Prognosis, and Targeted Therapy Response in Intrahepatic Cholangiocarcinoma.

Importance: BRAF variants are associated with tumor progression; however, the prevalence of BRAF variant subtypes and their association with disease characteristics, prognosis, and targeted therapy response in patients with intrahepatic cholangiocarcinoma (ICC) are largely unknown. Objective: To explore the association of BRAF variant subtypes with disease characteristics, prognosis, and targeted therapy response in patients with ICC. Design, Setting, and Participants: In this cohort study, 1175 patients who underwent curative resection for ICC from January 1, 2009, through December 31, 2017, were evaluated at a single hospital in China. Whole-exome sequencing, targeted sequencing, and Sanger sequencing were performed to identify BRAF variants. The Kaplan-Meier method and log-rank test were used to compare overall survival (OS) and disease-free survival (DFS). Univariate and multivariate analyses were performed using Cox proportional hazards regression. Associations between BRAF variants and targeted therapy response were tested in 6 BRAF-variant, patient-derived organoid lines and in 3 of the patient donors of those lines. Data were analyzed from June 1, 2021, to March 15, 2022. Interventions: Hepatectomy in patients with ICC. Main Outcomes and Measures: The association of BRAF variant subtypes with OS and DFS. Results: Of 1175 patients with ICC, the mean (SD) age was 59.4 (10.4) years and 701 (59.7%) were men. A total of 20 different subtypes of BRAF somatic variance affecting 49 patients (4.2%) were identified; V600E was the most frequent allele in this cohort, accounting for 27% of the identified BRAF variants, followed by K601E (14%), D594G (12%), and N581S (6%). Compared with patients with non-V600E BRAF variants, patients with BRAF V600E variants were more likely to have large tumor size (10 of 13 [77%] vs 12 of 36 [33%]; P = .007), multiple tumors (7 of 13 [54%] vs 8 of 36 [22%]; P = .04), and more vascular/bile duct invasion (7 of 13 [54%] vs 8 of 36 [22%]; P = .04). Multivariate analysis revealed that BRAF V600E variants, but not overall BRAF variants or non-V600E BRAF variants, were associated with poor OS (hazard ratio [HR], 1.87; 95% CI, 1.05-3.33; P = .03) and DFS (HR, 1.66; 95% CI, 1.03-2.97; P = .04). There were also broad differences among organoids with different BRAF variant subtypes in sensitivity to BRAF or MEK inhibitors. Conclusions and Relevance: The findings of this cohort study suggest that there are broad differences among organoids with different BRAF variant subtypes in sensitivity to BRAF or MEK inhibitors. Identifying and classifying BRAF variants may be able to help guide precise treatment for patients with ICC.

Association of KRAS Variant Subtypes With Survival and Recurrence in Patients With Surgically Treated Intrahepatic Cholangiocarcinoma.

Importance: KRAS variants are associated with tumor progression; however, the prevalence of KRAS variant subtypes and their association with survival and recurrence in patients with intrahepatic cholangiocarcinoma (ICC) after curative resection are largely unknown. Objective: To explore the prognostic association of KRAS variant subtypes with survival and recurrence in patients with ICC. Design, Setting, and Participants: In this cohort study, patients who underwent curative resection for ICC from January 2009 through December 2016 at a single hospital in China were recruited, and whole-exome sequencing, targeted sequencing, and Sanger sequencing were performed to identify KRAS variants. Kaplan-Meier and log-rank tests were used to compare overall survival (OS) and disease-free survival (DFS). Univariate and multivariate analyses were performed using the Cox proportional hazards regression model. Data were analyzed from April 2020 to January 2021. Interventions: Hepatectomy in patients with ICC. Main Outcomes and Measures: The association of KRAS variant subtypes with OS and DFS. Results: Of 1024 included patients with ICC, 621 (60.6%) were male, and the mean (SD) age was 59.2 (10.2) years. A total of 14 different subtypes of KRAS somatic variants affecting 127 patients (12.4%) were identified. G12D was the most frequent allele in this cohort, accounting for 55 of 127 identified KRAS variants (43.3%), followed by G12V (25 [19.7%]), G12C (9 [7.1%]), and G13D (8 [6.3%]). Compared with patients with wild-type KRAS, patients with variant KRAS were more likely to have high levels of carbohydrate antigen 19-9 (92 of 127 [72.4%] vs 546 of 897 [60.9%]; P = .01) and γ-glutamyltransferase (72 of 127 [56.7%] vs 420 of 897 [46.8%]; P = .04). Multivariable analysis revealed that G12 KRAS variants but not non-G12 KRAS variants were independently associated with worse OS (hazard ratio [HR], 1.69; 95% CI, 1.31-2.18; P < .001) and DFS (HR, 1.47; 95% CI, 1.16-1.88; P = .002). Among the patients with G12 KRAS variants, the G12V KRAS variant was the strongest prognostic determinant for the worst OS (HR, 3.05; 95% CI, 1.94-4.79; P < .001) and DFS (HR, 1.79; 95% CI, 1.13-2.85; P = .01). Conclusions and Relevance: In this cohort study, the distribution of KRAS variant subtypes was characterized in a large cohort of patients with ICC from China. The presence of G12 KRAS variants but not non-G12 KRAS variants was associated with worse survival and increased risk of recurrence. Patients with the G12V variant exhibited the worst outcomes in the whole cohort.

Tumor-associated neutrophils and macrophages interaction contributes to intrahepatic cholangiocarcinoma progression by activating STAT3.

BACKGROUND: Tumor-associated neutrophils (TANs) and macrophages (TAMs) can each influence cancer growth and metastasis, but their combined effects in intrahepatic cholangiocarcinoma (ICC) remain unclear. METHODS: We explored the distributions of TANs and TAMs in patient-derived ICC samples by multiplex immunofluorescent staining and tested their separate and combined effects on ICC in vitro and in vivo. We then investigated the mechanistic basis of the effects using PCR array, western blot analysis and ELISA experiments. Finally, we validated our results in a tissue microarray composed of primary tumor tissues from 359 patients with ICC. RESULTS: The spatial distributions of TANs and TAMs were correlated with each other in patient-derived ICC samples. Interaction between TANs and TAMs enhanced the proliferation and invasion abilities of ICC cells in vitro and tumor progression in a mouse xenograft model of ICC. TANs and TAMs produced higher levels of oncostatin M and interleukin-11, respectively, in co-culture than in monoculture. Both of those cytokines activated STAT3 signaling in ICC cells. Knockdown of STAT3 abolished the protumor effect of TANs and TAMs on ICC. In tumor samples from patients with ICC, increased TAN and TAM levels were correlated with elevated p-STAT3 expression. All three of those factors were independent predictors of patient outcomes. CONCLUSIONS: TANs and TAMs interact to promote ICC progression by activating STAT3.