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.

Stratifying ICIs-responsive tumor microenvironment in HCC: from parsing out immune-hypoxic crosstalk to clinically applicable MRI-radiomics models.

BACKGROUND: We aimed to redefine Immune checkpoint inhibitors (ICIs)-responsive "hot" TME and develop a corresponding stratification model to maximize ICIs-efficacy in Hepatocellular Carcinoma (HCC). METHODS: Hypoxic scores were designed, and the relevance to immunotherapy responses were validated in pan-cancers through single cell analysis. Multi-omics analysis using the hypoxic scores and immune infiltrate abundance was performed to redefine the ICIs-responsive TME subtype in HCC patients from TCGA (n = 363) and HCCDB database (n = 228). The immune hypoxic stress index (IHSI) was constructed to stratify the ICIs-responsive TME subtype, with exploring biological mechanism in vitro and in vivo. MRI-radiomics models were built for clinical applicability. RESULTS: The hypoxic scores were lower in the dominant cell-subclusters of responders in pan-cancers. The higher immune infiltrate-normoxic (HIN) subtype was redefined as the ICIs-responsive TME. Stratification of the HIN subtype using IHSI effectively identified ICIs-responders in Melanoma (n = 122) and urological cancer (n = 22). TRAF3IP3, the constituent gene of IHSI, was implicated in ICIs-relevant "immune-hypoxic" crosstalk by stimulating MAVS/IFN-I pathway under normoxic condition. MRI-radiomics models assessing TRAF3IP3 with HIF1A expression (AUC > 0.80) screened ICIs-Responders in HCC cohort (n = 75). CONCLUSION: The hypoxic-immune stratification redefined ICIs-responsive TME and provided MRI-Radiomics models for initial ICIs-responders screening, with IHSI facilitating further identification.

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.

Targeting HNRNPM Inhibits Cancer Stemness and Enhances Antitumor Immunity in Wnt-activated Hepatocellular Carcinoma.

BACKGROUND & AIMS: Cancer stemness and immune evasion are closely associated and play critical roles in tumor development and resistance to immunotherapy. However, little is known about the underlying molecular mechanisms that coordinate this association. METHODS: The expressions of heterogeneous nuclear ribonucleoprotein M (HNRNPM) in 240 hepatocellular carcinoma (HCC) samples, public databases, and liver development databases were analyzed. Chromatin immunoprecipitation assays were performed to explore the associations between stem-cell transcription factors and HNRNPM. HNRNPM-regulated alternative splicing (AS) and its binding motif were identified by RNA-seq and RIP-seq. HNRNPM-specific antisense oligonucleotides were developed to explore potential therapeutic targets in HCC. CD8+ T cells that were co-cultured with tumor cells were sorted by flow cytometry assays. RESULTS: We identified an elevated oncofetal splicing factor in HCC, HNRNPM, that unifies and regulates the positive association between cancer stemness and immune evasion. HNRNPM knockdown abolished HCC tumorigenesis and diminished cancer stem cell properties in vitro and in vivo. Mechanistically, HNRNPM regulated the AS of MBD2 by binding its flanking introns, whose isoforms played opposing roles. Although MBD2a and MBD2c competitively bound to CpG islands in the FZD3 promoter, MBD2a preferentially increased FZD3 expression and then activated the WNT/ß-catenin pathway. Interestingly, FZD3 and ß-catenin further provided additional regulation by targeting OCT4 and SOX2. We found that HNRNPM inhibition significantly promoted CD8+ T cell activation and that HNRNPM- antisense oligonucleotides effectively inhibited WNT/ß-catenin to enhance anti-programmed cell death protein-1 immunotherapy by promoting CD8+ T cell infiltration. CONCLUSIONS: HNRNPM has a tumor-intrinsic function in generating an immunosuppressive HCC environment through an AS-dependent mechanism and demonstrates proof of the concept of targeting HNRNPM in tailoring HCC immunotherapeutic approaches.

CircRPN2 Inhibits Aerobic Glycolysis and Metastasis in Hepatocellular Carcinoma.

Although circular RNAs (circRNA) are known to modulate tumor initiation and progression, their role in hepatocellular carcinoma (HCC) metastasis remains poorly understood. Here, three metastasis-associated circRNAs identified in a previous circRNA-sequencing study were screened and validated in two HCC cohorts. CircRPN2 was downregulated in highly metastatic HCC cell lines and HCC tissues with metastasis. Patients with HCC with lower circRPN2 levels displayed shorter overall survival and higher rates of cumulative recurrence. Mechanistic studies in vitro and in vivo revealed that circRPN2 binds to enolase 1 (ENO1) and accelerates its degradation to promote glycolytic reprogramming through the AKT/mTOR pathway, thereby inhibiting HCC metastasis. CircRPN2 also acted as a competing endogenous RNA for miR-183-5p, which increases forkhead box protein O1 (FOXO1) expression to suppress glucose metabolism and tumor progression. In clinical samples, circRPN2 expression negatively correlated with ENO1 and positively correlated with FOXO1, and expression of circRPN2, either alone or in combination with ENO1 and FOXO1, was a novel indicator of HCC prognosis. These data support a model wherein circRPN2 inhibits HCC aerobic glycolysis and metastasis via acceleration of ENO1 degradation and regulation of the miR-183-5p/FOXO1 axis, suggesting that circRPN2 represents a possible therapeutic target in HCC. SIGNIFICANCE: The circRNA circRPN2 is a potential prognostic biomarker and therapeutic target in hepatocellular carcinoma that suppresses aerobic glycolysis and metastasis.

Whole-genome sequencing reveals the evolutionary trajectory of HBV-related hepatocellular carcinoma early recurrence.

Patients with hepatocellular carcinoma (HCC) have poor long-term survival following curative resection because of the high rate of tumor early recurrence. Little is known about the trajectory of genomic evolution from primary to early-recurrent HCC. In this study, we performed whole-genome sequencing (WGS) on 40 pairs of primary and early-recurrent hepatitis B virus (HBV)-related HCC tumors from patients who received curative resection, and from four patients whose primary and recurrent tumor were extensively sampled. We identified two recurrence patterns: de novo recurrence (18/40), which developed genetically independently of the primary tumor and carried different HCC drivers, and ancestral recurrence (22/40), which was clonally related to the primary tumor and progressed more rapidly than de novo recurrence. We found that the recurrence location was predictive of the recurrence pattern: distant recurrence tended to display the de novo pattern, whereas local recurrence tended to display the ancestral pattern. We then uncovered the evolutionary trajectories based on the subclonal architecture, driver-gene mutations, and mutational processes observed in the primary and recurrent tumors. Multi-region WGS demonstrated spatiotemporal heterogeneity and polyclonal, monophyletic dissemination in HCC ancestral recurrence. In addition, we identified recurrence-specific mutations and copy-number gains in BCL9, leading to WNT/ß-catenin signaling activation and an immune-excluded tumor microenvironment, which suggests that BCL9 might serve as a new therapeutic target for recurrent HCC. Collectively, our results allow us to view with unprecedented clarity the genomic evolution during HBV-related HCC early recurrence, providing an important molecular foundation for enhanced understanding of HCC with implications for personalized therapy to improve patient survival.

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.

Associations among the mutational landscape, immune microenvironment, and prognosis in Chinese patients with hepatocellular carcinoma.

Recent studies suggested that the immune microenvironment and mutational landscape are associated with the response to immune-based therapy in several types of cancer. The roles of those factors in Chinese HCC remain largely unknown. In this study, we obtained 182 FFPE samples of HCC cohort that were previously subjected to NGS (49 WGS, 18 WES, and 115 targeted sequencing). We performed immunohistochemistry to detect CD3, CD4, CD8, CD57, Foxp3, CD68, CD66b, and PD-L1 expression in the samples. We identified diverse associations between the mutational landscape and the immune microenvironment in the HCC samples. High mutational burden and an aristolochic acid-dominated mutational signature were both correlated with elevated tumoral PD-L1 expression and CD3+ T-cell infiltration and high numbers of CD68+ TAMs and CD66b+ TANs. CD4+ and CD8+ T cells exhibited lower infiltration levels in tumors with mutations in AXIN1/CTNNB1 and in tumors with aflatoxin-dominant mutational signatures. Moreover, tumors with TP53 mutations had less CD8+ T-cell infiltration and more Foxp3+ Treg-cell infiltration than those without TP53 mutations. Kaplan-Meier survival analysis revealed that the presence of CD8+, Foxp3+, CD66b+, or CD68+ immune cells; tumoral PD-L1 expression alone; or the presence of CD8+ or Foxp3+ cells combined with TP53 mutation were predictive of recurrence and poor overall survival after curative resection. In conclusion, the association between the mutational landscape and the immune microenvironment warrants further analysis to determine its impact on patient outcomes to guide personalized immune-based therapy for Chinese patients with HCC.

Peritumoral plasmacytoid dendritic cells predict a poor prognosis for intrahepatic cholangiocarcinoma after curative resection.

BACKGROUND: Plasmacytoid dendritic cells (pDCs) are present in various primary and metastatic human neoplasms; however, their clinical significance in intrahepatic cholangiocarcinoma is not clear. METHODS: To evaluate pDCs' distributions in and around tumors as well as their potential function and predictive value for prognosis in patients undergoing curative resection, we performed immunohistochemistry to examine the expression of pDC marker BDCA2, and CD3, CD4, CD8 and Foxp3 in intratumoral and peritumoral tissues from 359 patients with intrahepatic cholangiocarcinoma and compared with prognostic and clinicopathologic factors. RESULTS: Results showed that patients with high numbers of BDCA2+ pDCs in peritumoral tissues were more likely to have elevated levels of carbohydrate antigen 19-9 and gamma-glutamyl transferase, larger and more tumors, advanced tumor-node-metastasis staging, more vascular/bile duct invasion, and lymphatic metastasis in association with greater chance of recurrence and shorter overall survival. Peritumoral tissues with larger numbers of pDCs also showed increased Foxp3+ regulatory T cell infiltration, both of which were found to be independent factors for predicting time to recurrence and overall survival. By contrast, patient outcomes were not associated with the presence of intratumoral pDCs. CONCLUSIONS: Peritumoral pDC infiltration may indicate an immune tolerogenic peritumor microenvironment and can be used to predict a poor prognosis for patients undergoing curative resection for intrahepatic cholangiocarcinoma.

Prognostic factors influencing outcome in adult liver transplantation using hypernatremic organ donation after brain death.

BACKGROUND: Hypernatremic donors was regarded as the expanded criteria donors in liver transplantation. The study was to investigate the effects of donor hypernatremia on the outcomes of liver transplantation and identify the prognostic factors possibly contributing to the poor outcomes. METHODS: Donor serum sodium levels before procurement were categorized as normal sodium (< 155 mmol/L), moderate high sodium (155-170 mmol/L), and severe high sodium (≥ 170 mmol/L). Furthermore, we subdivided the 142 hypernatremic donors (≥ 155 mmol/L) into two subgroups: subgroup A, the exposure time of liver grafts from hypernatremia to reperfusion was < 36 h; and subgroup B, the exposure time was ≥ 36 h. The outcomes included initial graft function, survival rates of grafts and recipients, graft loss and early events within the first year following liver transplantation. RESULTS: There were no significant differences in the 1-year survival rates of grafts and recipients, 1-year graft loss rates and early events among the normal, moderate high and severe high sodium groups. However, the overall survival rates of grafts and recipients in subgroup A were significantly higher than those in subgroup B. Cox model showed that the exposure time (HR = 1.117; 95% CI: 1.053-1.186; P < 0.001), cold ischemia time (HR = 1.015; 95% CI: 1.006-1.024; P = 0.001) and MELD (HR = 1.061; 95% CI: 1.003-1.121; P = 0.037) were the important prognostic factors contributing to the poor outcomes of recipients with hypernatremic donors. CONCLUSIONS: The level of donor sodium immediately before organ procurement does not have negative effects on the early outcomes following adult liver transplantation. For hypernatremia liver donors, minimization of the exposure time from hypernatremia to reperfusion is critical to prevent graft loss.

MACROD2 deficiency promotes hepatocellular carcinoma growth and metastasis by activating GSK-3ß/ß-catenin signaling.

Structural variations (SVs) influence the development and progression of multiple types of cancer. The genes affected by SVs in hepatocellular carcinoma (HCC) and their contribution to tumor growth and metastasis remain unknown. In this study, through whole-genome sequencing (WGS), we identified MACROD2 as the gene most frequently affected by SVs, which were associated with low MACROD2 expression levels. Low MACROD2 expression was predictive of tumor recurrence and poor overall survival. MACROD2 expression was decreased in HCC cell lines, especially those with high metastatic potential. MACROD2 knockdown in HCC cells markedly enhanced proliferation and invasiveness in vitro and tumor progression in vivo and promoted epithelial-mesenchymal transition (EMT). By contrast, MACROD2 overexpression reversed EMT and inhibited HCC growth and metastasis. Mechanistically, MACROD2 deficiency suppressed glycogen synthase kinase-3ß (GSK-3ß) activity and activated ß-catenin signaling, which mediated the effect of MACROD2 on HCC. In clinical HCC samples, decreased MACROD2 expression was correlated with the activation of GSK-3ß/ß-catenin signaling and the EMT phenotype. Overall, our results revealed that MACROD2 is frequently affected by SVs in HCC, and its deficiency promotes tumor growth and metastasis by activating GSK-3ß/ß-catenin signaling.

HNRNPAB-regulated lncRNA-ELF209 inhibits the malignancy of hepatocellular carcinoma.

Our previous study demonstrated that heterogeneous nuclear ribonucleoprotein AB (HNRNPAB) is a key gene that facilitates metastasis of hepatocellular carcinoma (HCC). However, the molecular mechanisms behind this relationship are not fully understood. In our study, we utilized long-noncoding RNA (lncRNA) microarrays to identify a HNRNPAB-regulated lncRNA named lnc-ELF209. Our findings from chromatin immunoprecipitation assays indicate that HNRNPAB represses lnc-ELF209 transcription by directly binding to its promoter region. We also analyzed clinical samples from HCC patients and cell lines with quantitative real-time polymerase chain reactions, RNA in situ hybridization and immunohistochemistry, and found that there is a negative relationship between HNRNPAB and lnc-ELF209 expression. Up/downregulation assays and rescue assays indicate that lnc-ELF209 inhibits cell migration, invasion and epithelial-mesenchymal transition regulated by HNRNPAB. This suggests a new regulatory mechanism for HNRNPAB-promoted HCC progression. RNA pull-down and LC-MS/MS were used to determine triosephosphate isomerase, heat shock protein 90-beta and vimentin may be involved in the tumor-suppressed function of lnc-ELF209. Furthermore, we found lnc-ELF209 could stabilize TPI protein expression. We also found that lnc-ELF209 overexpression in HCCLM3 cell resulted in a lower rate of lung metastatic, which suggested a less aggressive HCC phenotype. Collectively, these findings offer new insights into the regulatory mechanisms that underlie HNRNPAB cancer-promoting activities and demonstrate that lnc-ELF209 is a HNRNPAB-regulated lncRNA that may play an important role in the inhibition of HCC progression.

Circular RNA Sequencing Identifies CircASAP1 as a Key Regulator in Hepatocellular Carcinoma Metastasis.

BACKGROUND AND AIMS: There is growing evidence that single-stranded, circular RNA (circRNA) plays a key role in the development of certain cancers, including hepatocellular carcinoma (HCC). It is less clear, however, what role circRNA plays in HCC metastasis. APPROACH AND RESULTS: In this study, through circRNA sequencing, we identified a circRNA: circASAP1 (a circRNA derived from exons 2 and 3 of the ASAP1 gene, hsa_circ_0085616), which is associated with pulmonary metastasis after curative resection in patients with HCC. CircASAP1 was overexpressed in HCC cell lines with high metastatic potential and in metastatic HCCs. In vitro, circASAP1 promoted cell proliferation, colony formation, migration, and invasion, and in vivo, it enhanced tumor growth and pulmonary metastasis. Mechanism studies showed that circASAP1 acts as a competing endogenous RNA for microRNA 326 (miR-326) and microRNA 532-5p (miR-532-5p), both of which are tumor suppressors in HCC. We found that mitogen-activated protein kinase (MAPK) 1 and colony stimulating factor (CSF)-1 were direct common targets for microRNA 326 (miR-326) and microRNA 532-5p (miR-532-5p), which were regulated by circASAP1. CircASAP1 promotes HCC cell proliferation and invasion by regulating miR-326/miR-532-5p-MAPK1 signaling and, furthermore, mediates tumor-associated macrophage infiltration by regulating the miR-326/miR-532-5p-CSF-1 pathway. Clinical HCC samples exhibited a positive correlation between circASAP1 expression and levels of CSF-1, MAPK1, and CD68+ tumor-associated macrophages, all of which were predictive of patient outcomes. CONCLUSION: We identified circASAP1 as a key regulator of HCC metastasis that acts on miR-326/miR-532-5p-MAPK1/CSF-1 signaling and serves as a prognostic predictor in patients with HCC.

The miR-561-5p/CX3CL1 Signaling Axis Regulates Pulmonary Metastasis in Hepatocellular Carcinoma Involving CX3CR1+ Natural Killer Cells Infiltration.

Natural killer (NK) cell can inhibit tumor initiation and regulates metastatic dissemination, acting as key mediators of the innate immune response. Intrinsic factors modulating NK cells infiltration and its anticancer activity remain poorly characterized. We investigated the roles of dysregulation of micro(mi)RNAs and NK cells in progression of hepatocellular carcinoma (HCC). Methods: Small RNA sequencing were used to detect the miRNA profiles of tumor tissues from HCC patients with (n=14) or without (n=13) pulmonary metastasis and HCC cell lines with different pulmonary metastatic potentials. Chemokine expression profiling and bioinformatics were used to detect the downstream target of candidate target. In gain- and loss-of-function assays were used to investigate the role of miRNA in HCC progression. Different subsets of NK cells were isolated and used for chemotaxis and functional assays in vivo and in vitro. In situ hybridization and immunohistochemical analyses were performed to detect the expression of miRNA in tumor tissues from 242 HCC patients undergoing curative resection from 2010. Results: Three miRNAs (miR-137, miR-149-5p, and miR-561-5p) were identified to be associated with pulmonary metastasis in patients with HCC. miR-561-5p was most highly overexpressed in metastatic HCC tissues and high-metastatic-potential HCC cell lines. In gain- and loss-of-function assays in a murine model, miR-561-5p promoted tumor growth and spread to the lungs. Yet, miR-561-5p did not appear to affect cellular proliferation and migration in vitro. Bioinformatics and chemokine expression profiling identified chemokine (C-X3-C motif) ligand 1 (CX3CL1) as a potential target of miR-561-5p. Furthermore, miR-561-5p promoted tumorigenesis and metastasis via CX3CL1-dependent regulation of CX3CR1+ NK cell infiltration and function. CX3CR1+ NK cells demonstrated stronger in vivo anti-metastatic activity relative to CX3CR1- NK cells. CX3CL1 stimulated chemotactic migration and cytotoxicity in CX3CR1+ NK cells via STAT3 signaling. Blockade of CX3CL1, CX3CR1, or of pSTAT3 signaling pathways attenuated the antitumor responses. Clinical samples exhibited a negative correlation between miR-561-5p expression and levels of CX3CL1 and CX3CR1+ NK cells. High miR-561-5p abundance, low CX3CL1 levels, and low numbers of CX3CR1+ NK cells were associated with adverse prognosis. Conclusion: We delineated a miR-561-5p/CX3CL1/NK cell axis that drives HCC metastasis and demonstrated that CX3CR1+ NK cells serve as potent antitumor therapeutic effectors.

Genomic sequencing identifies WNK2 as a driver in hepatocellular carcinoma and a risk factor for early recurrence.

BACKGROUND & AIMS: Early recurrence of hepatocellular carcinoma (HCC) after curative resection is common. However, the association between genetic mechanisms and early HCC recurrence, especially in Chinese patients, remains largely unknown. METHODS: We performed whole-genome sequencing (49 cases), whole-exome sequencing (18 cases), and deep targeted sequencing (115 cases) on 182 primary HCC samples. Focusing on WNK2, we used Sanger sequencing and qPCR to evaluate all the coding exons and copy numbers of that gene in an additional 554 HCC samples. We also explored the functional effect and mechanism of WNK2 on tumor growth and metastasis. RESULTS: We identified 5 genes (WNK2, RUNX1T1, CTNNB1, TSC1, and TP53) harboring somatic mutations that correlated with early tumor recurrence after curative resection in 182 primary HCC samples. Focusing on WNK2, the overall somatic mutation and copy number loss occurred in 5.3% (39/736) and 27.2% (200/736), respectively, of the total 736 HCC samples. Both types of variation were associated with lower WNK2 protein levels, higher rates of early tumor recurrence, and shorter overall survival. Biofunctional investigations revealed a tumor-suppressor role of WNK2: its inactivation led to ERK1/2 signaling activation in HCC cells, tumor-associated macrophage infiltration, and tumor growth and metastasis. CONCLUSIONS: Our results delineate genomic events that characterize Chinese HCCs and identify WNK2 as a driver of early HCC recurrence after curative resection. LAY SUMMARY: We applied next-generation sequencing and conducted an in-depth genomic analysis of hepatocellular carcinomas from a Chinese patient cohort. The results delineate the genomic events that characterize hepatocellular carcinomas in Chinese patients and identify WNK2 as a driver associated with early tumor recurrence after curative resection.

Intratumoral plasmacytoid dendritic cells as a poor prognostic factor for hepatocellular carcinoma following curative resection.

Plasmacytoid dendritic cells (pDCs) are present in various primary and metastatic human neoplasms; however, their clinical significance in hepatocellular carcinoma (HCC) is unclear. In this study, we investigated the distribution, prognostic value, and potential function of pDCs in HCC patients undergoing curative resection. We performed immunohistochemical analyses of whole tumor sections from 224 patients to assess the expression of BDCA2, CD3, CD4, CD8, Foxp3, granzyme B, IL-17, and CD34. The findings were validated using tissue microarrays from another two independent cohorts totaling 841 HCC patients undergoing curative resection. Our results demonstrated that high numbers of BDCA2+ pDCs within tumors correlated with high alpha-fetoprotein levels, greater vascular invasion, advanced tumor-node-metastasis stage, shorter overall survival, and a higher recurrence rate. However, patient outcomes were not associated with pDCs in peritumoral stromal or nontumor tissues. Furthermore, an increase in intratumoral pDCs was associated with increased intratumoral infiltration of Foxp3+ regulatory T cells and IL-17-producing cells and correlated with tumor vascular density. Univariate and multivariate analyses revealed that the presence of intratumoral pDCs alone or in combination with regulatory T and/or IL-17-producing cells was an independent predictor of time to recurrence and overall survival. In conclusion, our study demonstrated that intratumoral infiltration by pDCs is a novel indicator for poor prognosis in patients with HCC, possibly through the induction of an immune tolerogenic and inflammatory tumor microenvironment comprising regulatory T and IL-17-producing cells. An assessment of the combination of these cells represents a superior predictor of patient outcome.

Development and validation of a new tumor-based gene signature predicting prognosis of HBV/HCV-included resected hepatocellular carcinoma patients.

BACKGROUND: Due to the phenotypic and molecular diversity of hepatocellular carcinomas (HCC), it is still a challenge to determine patients' prognosis. We aim to identify new prognostic markers for resected HCC patients. METHODS: 274 patients were retrospectively identified and samples collected from Zhongshan hospital, Fudan University. We analyzed the gene expression patterns of tumors and compared expression patterns with patient survival times. We identified a "9-gene signature" associated with survival by using the coefficient and regression formula of multivariate Cox model. This molecular signature was then validated in three patients cohorts from internal cohort (n = 69), TCGA (n = 369) and GEO dataset (n = 80). RESULTS: We identified 9-gene signature consisting of ZC2HC1A, MARCKSL1, PTGS1, CDKN2B, CLEC10A, PRDX3, PRKCH, MPEG1 and LMO2. The 9-gene signature was used, combined with clinical parameters, to fit a multivariable Cox model to the training cohort (concordance index, ci = 0.85), which was successfully validated (ci = 0.86 for internal cohort; ci = 0.78 for in silico cohort). The signature showed improved performance compared with clinical parameters alone (ci = 0.70). Furthermore, the signature predicted patient prognosis than previous gene signatures more accurately. It was also used to stratify early-stage, HBV or HCV-infected patients into low and high-risk groups, leading to significant differences in survival in training and validation (P < 0.001). CONCLUSIONS: The 9-gene signature, in which four were upregulated (ZC2HC1A, MARCKSL1, PTGS1, CDKN2B) and five (CLEC10A, PRDX3, PRKCH, MPEG1, LMO2) were downregulated in HCC with poor prognosis, stratified HCC patients into low and high risk group significantly in different clinical settings, including receiving adjuvant transarterial chemoembolization and especially in early stage disease. This new signature should be validated in prospective studies to stratify patients in clinical decisions.