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.

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.

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.

LINC01133 promotes hepatocellular carcinoma progression by sponging miR-199a-5p and activating annexin A2.

BACKGROUND: Long noncoding RNAs (lncRNAs) are functionally associated with cancer development and progression. Although gene copy number variation (CNV) is common in hepatocellular carcinoma (HCC), it is not known how CNV in lncRNAs affects HCC progression and recurrence. We aimed to identify a CNV-related lncRNA involved in HCC progression and recurrence and illustrate its underlying mechanisms and prognostic value. METHODS: We analyzed the whole genome sequencing (WGS) data of matched cancerous and noncancerous liver samples from 49 patients with HCC to identify lncRNAs with CNV. The results were validated in another cohort of 238 paired HCC and nontumor samples by TaqMan copy number assay. We preformed Kaplan-Meier analysis and log-rank test to identify lncRNA CNV with prognostic value. We conducted loss- and gain-of-function studies to explore the biological functions of LINC01133 in vitro and in vivo. The competing endogenous RNAs (ceRNAs) mechanism was clarified by microRNA sequencing (miR-seq), quantitative real-time PCR (qRT-PCR), western blot, and dual-luciferase reporter assays. We confirmed the binding mechanism between lncRNA and protein by RNA pull-down, RNA immunoprecipitation, qRT-PCR, and western blot analyses. RESULTS: Genomic copy numbers of LINC01133 were increased in HCC, which were positively related with the elevated expression of LINC01133. Increased copy number of LINC01133 predicted the poor prognosis in HCC patients. LINC01133 overexpression in HCC cells promoted proliferation and aggressive phenotypes in vitro, and facilitated tumor growth and lung metastasis in vivo, whereas LINC01133 knockdown had the opposite effects. LINC01133 sponged miR-199a-5p, resulting in enhanced expression of SNAI1, which induced epithelial-to-mesenchymal transition (EMT) in HCC cells. In addition, LINC01133 interacted with Annexin A2 (ANXA2) to activate the ANXA2/STAT3 signaling pathway. CONCLUSIONS: LINC01133 promotes HCC progression by sponging miR-199a-5p and interacting with ANXA2. LINC01133 CNV gain is predictive of poor prognosis in patients with HCC.

Decreased expression of serum semaphorin 3B is associated with poor prognosis of patients with hepatocellular carcinoma.

Semaphorin 3B (SEMA-3B), which belongs to the semaphorin family, has an important role in cell apoptosis and inhibition of angiogenesis. A previous study by our group revealed that SEMA-3B was downregulated in tumor tissues of patients with hepatocellular carcinoma (HCC) and exerts anti-motility and anti-invasion effects on tumor cells. However, the serum levels of SEMA-3B and their clinical significance have remained elusive; therefore, the aim of the present study was to monitor its expression in HCC and investigate its clinical significance. ELISA was used to determine the serum levels of SEMA-3B in 132 patients with HCC and 57 healthy individuals. The association between SEMA-3B and clinicopathological parameters was investigated. Serum SEMA-3B was indicated to be significantly decreased in patients with HCC as compared with that in the controls (P<0.05) and it was negatively associated with tumor size (P=0.039), encapsulation (P=0.002) and TNM stage (P=0.034). The prognosis of patients with low expression of SEMA-3B was poor. In conclusion, the results of the present study revealed that serum SEMA-3B is decreased in HCC and is negatively associated with prognosis; therefore, it may be used as a prognostic marker in HCC.