Plasma-only circulating tumor DNA analysis detects minimal residual disease and predicts early relapse in hepatocellular carcinoma patients undergoing curative resection

Background: Minimal residual disease (MRD) is considered an essential factor leading to relapse within 2 years (early relapse) after radical surgery, which is challenging to be detected by conventional imaging. Circulating tumor DNA (ctDNA) provides a novel approach for detecting MRD and predicting clinical outcomes. Here, we tried to construct a fixed panel for plasma-only ctDNA NGS to enable tumor-uninformed MRD detection in hepatocellular carcinoma (HCC).

Methods: Here, we performed the followings: (i) profiling genomic alteration spectrum of ctDNA from the Chinese HCC cohort consisting of 493 individuals by NGS; (ii) screening of MRD monitoring genes; and (iii) performance evaluation of MRD monitoring genes in predicting early relapse in the ZJZS2020 cohort comprising 20 HCC patients who underwent curative resection.

Results: A total of 493 plasma samples from the Chinese HCC cohort were detected using a 381/733-gene NGS panel to characterize the mutational spectrum of ctDNA. Most patients (94.1%, 464/493) had at least one mutation in ctDNA. The variants fell most frequently in TP53 (45.1%), LRP1B (20.2%), TERT (20.2%), FAT1 (16.2%), and CTNNB1 (13.4%). By customized filtering strategy, 13 MRD monitoring genes were identified, and any plasma sample with one or more MRD monitoring gene mutations was considered MRD-positive. In the ZJZS2020 cohort, MRD positivity presented a sensitivity of 75% (6/8) and a specificity of 100% (6/6) in identifying early postoperative relapse. The Kaplan-Meier analysis revealed a significantly short relapse-free survival (RFS; median RFS, 4.2 months vs. NR, P=0.002) in the MRD-positive patients versus those with MRD negativity. Cox regression analyses revealed MRD positivity as an independent predictor of poor RFS (HR 13.00, 95% CI 2.60-69.00, P=0.002).

Conclusions: We successfully developed a 13-gene panel for plasma-only MRD detection, which was effective and convenient for predicting the risk of early postoperative relapse in HCC.

Introduction

Liver cancer is the second leading cause of cancer death, and its incidence is growing worldwide (1–4). Over a million people per year are expected to be liver cancer patients by 2025 (5). Hepatocellular carcinoma (HCC) is the main type of primary liver cancer, comprising approximately 90% of patients. Although hepatectomy is a widely accepted treatment option for HCC patients with good liver function, the relapse rate of up to 60%-70% within 5 years after surgery remains a severe problem (6–8).

Relapse in multiple solid tumors, including HCC, can be divided into early relapse (≤2 years) mainly caused by minimal residual disease (MRD) following resection, and late relapse (>2 years) caused by de novo tumors arising in a microenvironment predisposed to carcinogenesis (9–14). Early relapse has been reported to account for over 60% of all relapsed HCC events (9, 15–19). Consequently, there has been much interest in detecting and eliminating MRD to prevent relapse or for early treatment of degeneration. Current postoperative surveillance methods are not sensitive or specific enough to detect MRD, such as monitoring clinical symptoms, tumor markers, and routine imaging. As a result, improving current therapeutic strategies and preventing recurrence might be achieved by establishing more precise MRD detection approaches.

Circulating cell-free DNA (cfDNA) is extracellular nucleic acid fragments released into the bloodstream due to apoptosis and necrosis from both healthy and malignant cells. There has been strong evidence that postoperative tumor-derived cfDNA (ctDNA) detection is correlated with MRD and could identify patients at high risk of relapse (20–23). Currently, there are two available ctDNA detection strategies for monitoring MRD: tumor-informed and tumor-uninformed (also referred to as tumor-agnostic, tumor-naïve, or plasma-only) assay.

The tumor-informed approach relies on tumor tissue sequencing to identify tumor-derived alterations for the design of patient-specific targeted gene panels for ctDNA tracking, which has presented effectiveness for monitoring MRD in some solid tumors after curative-intent treatment, including colon, lung, and pancreas cancer (24–27). A potential limitation is that tumor heterogeneity in both space and time can affect its performance and might generate some false negative results (28, 29). Additionally, designing individualized next-generation sequencing (NGS) panels also lengthens the turnaround time versus the fixed NGS panel. The tumor-agnostic approach only requires plasma cfDNA sequencing with a fixed panel, which endows it with the advantages of noninvasiveness, convenience, cost-effectiveness, and rapid turnaround time. The method’s drawback is the lack of sensitivity, which can be improved by incorporating serial longitudinal surveillance samples and examining a variety of biomarkers, such as ctDNA mutation and ctDNA methylation. Given the attractive advantages and redeemable disadvantages, this approach has increasingly been investigated in multiple solid tumors. Excitingly, one recent study reported that plasma-only MRD detection presented favorable sensitivity and specificity for predicting recurrence in colorectal cancer patients undergoing curative-intent surgery, comparable to the tumor-informed approach (30). However, few studies have examined whether a plasma-only ctDNA assay can identify MRD in HCC with clinically meaningful specificity and sensitivity.

In this study, we first performed a plasma-only ctDNA assay integrating genomic signatures to identify gene candidates associated with MRD in two cohorts. We then determined whether these genes could reliably identify MRD and predict early relapse in postoperative patients with HCC.

Materials and methods

Study design and population

There are three major phases in this study: (i) profiling genomic alteration spectrum of HCC using ctDNA from the Chinese HCC cohort; (ii) screening of MRD monitoring genes via excluding the genes mutated at high frequency in ctDNA of the ZJ2020 cohort with patients having relapse-free survival (RFS) more than 2 years after radical resection from the comprehensive mutational spectrum of ctDNA in the Chinese HCC patient cohort; and (iii) performance evaluation of MRD monitoring genes in predicting early relapse in the ZJZS2020 cohort comprising HCC patients underwent curative surgical resection. The study design and the CONSORT participant flow diagram are summarized in Supplementary Figure 1.

The Chinese HCC patient cohort consisted of 493 patients from the Zhujiang Hospital of Southern Medical University and the First Affiliated Hospital, Sun Yat-Sen University, who underwent plasma-based NGS using a 381/733-gene panel (3D Medicines Inc. Shanghai, China) between January 6, 2017 and June 2, 2020. The gene lists of these two panels were attached in Supplementary Tables 1, 2. The Venn diagram revealed 294 shared genes between these two panels (Supplementary Figure 2).

The ZJ2020 cohort included 24 operable HCC patients who underwent radical hepatectomy at the Zhujiang Hospital of Southern Medical University from May 25, 2010 to October 23, 2020, and remained disease-free for greater than two years postoperatively. Of 24 patients, 10 received transcatheter arterial chemoembolization (TACE) after surgery. The blood samples (≥10 ml per case) were collected within 1-4 weeks after surgery or completion of TACE if they received TACE. The plasma ctDNA was analyzed using a 381/733-gene NGS panel.

The ZJZS2020 cohort consisted of 20 operable HCC patients who had curative resection at the Zhujiang Hospital of Southern Medical University from July 8, 2019 to December 28, 2020. The blood samples were collected about one week after surgery and submitted to NGS with a 381/733-gene panel.

In the two surgery cohorts, ZJ2020 and ZJZS2020, clinicopathological data were collected through review of the medical records using a standardized case report form, including age, sex, liver function, tumor stage, maximum tumor diameter, tumor number, etiology, alpha-foetoprotein (AFP) level, microvascular invasion (MVI), and portal vein tumor thrombus (PVTT).

This study was conducted following the principles of the Declaration of Helsinki and approved by the Ethics Committee of the Zhujiang Hospital of Southern Medical University (Approval number 2022-KY-150-01). Written informed consent was obtained from all participants in this study.

Blood sample processing and cfDNA isolation

Whole blood (2 × 10 ml) from every patient was drawn into cell-free DNA BCT tubes (Streck). Blood samples were centrifuged, and cfDNA was extracted and quantified using the QiAmp Circulating Nucleic Acid Kit (Qiagen) and the Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific), respectively, as per the manufacturer’s instructions.

Library preparation and targeted capture

cfDNA libraries were constructed and barcoded with unique molecular identifiers (UMI). The libraries were then PCR-amplified and purified for target enrichment. For targeted capture, indexed libraries were subjected to probe-based hybridization with the customized NGS panel.

DNA sequencing, data processing, and variant calling

The captured libraries were subjected to 100-bp paired-end sequencing on an Illumina NovaSeq 6000. The average effective sequencing depth for the 381-gene and 733-gene NGS panels was 5161× (range, 1840-10664) and 6337× (range, 2194-14250), respectively. Raw reads were mapped to the reference human genome hg19. An in-house developed software was applied to generate a duplex consensus sequence by incorporating a dual UMI at the end of the DNA fragments. An in-house loci-specific variant detection model based on a binary test was also used to improve specificity, especially for variants with low allele frequency. Following these steps, the variants were filtered by supporting read count, strand bias, base quality, and mapping quality. Variant calling was also optimized to identify variants at short tandem repeat regions. Single-nucleotide polymorphism (SNPs) and indels were annotated using ANNOVAR. Only missense, stopgain, frameshift, and non-frameshift indel mutations were retained.

Statistical analysis

Gene mutation analysis was performed by R package maftools (31). Wilcoxon test was adopted to estimate differences between continuous variables. The difference in proportions between groups was assessed using the Chi-square test. Kaplan-Meier analysis with a log-rank test was conducted to compare RFS between groups. Cox regression analyses were applied to identify variables associated with RFS. Statistical analysis was performed using R version 4.0.3. All results were considered statistically significant when the P value was less than 0.05.

Results

Mutational landscape in HCC by ctDNA profiling

A total of 493 plasma samples from Chinese patients with HCC were detected to characterize the mutational spectrum of ctDNA, wherein 285 plasma samples were tested with the 381-gene panel, and 208 were determined with the 733-gene panel. Most patients (94.1%, 464/493) had at least one mutation in ctDNA. The top 30 mutated genes are shown in Figure 1. The top five most frequently mutated genes were TP53 (45.1%), LRP1B (20.2%), TERT (20.2%), FAT1 (16.2%), and CTNNB1 (13.4%) (Figure 1A). Variant classification showed that the missense mutations were the most common (Figure 1B). SNPs were the most common variant types. Among the single-nucleotide variation (SNV) class, C > T was the most common base substitution (Figure 1C). The median variant number per sample was 4 (range, 1-33) (Figure 1D). The pathway analysis showed that the top affected pathways were TP53, TGF-β, NRF2, Cell_cycle, RTK-RAS, and PI3K which occurred in a substantial number of HCC plasma samples (Figure 1E). Besides, the analysis of mutation interactions identified one mutually exclusive and 28 co-occurring events in the top 25 highly mutated genes (Figure 2A). We also used the OncodriveCLUST algorithm to identify potential driver mutations in HCC. The top six driver genes ordered according to z-score were RBM10, CHEK2, KRAS, BARD1, EZH2, and IDH1 (Figure 2B).

FIGURE 1

Figure 1 Mutational landscape by ctDNA profiling in the Chinese HCC cohort. (A) The ratio of the top 30 genes in the number of mutations. (B) Frequency of different mutation classifications. (C) Frequency of SNV class. (D) Variant number per HCC samples. (E) The number of mutated genes contained in the number of mutated samples per pathway. HCC, hepatocellular carcinoma; SNV, single nucleotide variation.

FIGURE 2

Figure 2 Identification of mutation interactions and cancer driver genes in the Chinese HCC cohort. (A) The co-occurrence or exclusive associations between top 25 mutation genes. (B) Detecting cancer driver genes based on positional clustering in HCC. Each dot represents a gene, and the size of the dot represents the number of clusters (mentioned inside square brackets) within which a fraction (X-axis) of total variants is accumulated. HCC, hepatocellular carcinoma.

Screening of MRD monitoring genes

It has been well documented that early relapse after surgery is primarily due to MRD (9–12). Thus we subtracted the genes mutated at high frequency in ctDNA of the ZJ2020 cohort from the comprehensive mutational spectrum of ctDNA in the Chinese HCC cohort to screen for MRD monitoring genes. Baseline characteristics of the ZJ2020 cohort are listed in Supplementary Table 3. This cohort comprised 24 HCC patients with RFS more than 2 years after radical resection. As of May 8, 2019, the median postoperative follow-up was 35.5 months. Twenty-one patients remained on RFS, and three experienced recurrence at postoperative 34.7, 46.4, and 52.8 months, respectively. We identified ≥1 mutation in 23 of 24 (95.8%) patients. The most frequently mutated genes included FAT1, FAM135B, BRCA2, ERBB2, LRP1B, TSC2, AR, ARID1A, ATM, MSH6, NOTCH3, RB1, and SPEN, with each being mutated in at least 10% of patients (Figure 3A). These 14 genes were considered to be unrelated to early relapse and were filtered out from the 27 genes with a mutation frequency of > 5% in the Chinese HCC cohort (32). The remaining 13 genes were regarded as MRD monitoring genes. As shown in Figure 3B, these 13 genes were TP53, TERT, CTNNB1, APC, RMB10, NTRK3, NOTCH1, NOTCH2, NF1, CREBBP, GLI3, CDKN2A, and EZH2. Any plasma sample with one or more MRD monitoring gene mutations was considered MRD-positive and otherwise MRD-negative (33, 34).

FIGURE 3

Figure 3 Screening of MRD monitoring genes. (A) The top 30 frequently mutated genes in the ZJ2020 cohort with HCC patients having RFS more than 2 years after radical resection. (B) The 14 genes with a mutation frequency of >10% in the ZJ2020 cohort are filtered out from the 27 genes with a mutation frequency of >5% in the Chinese HCC cohort, and the remaining 13 genes are used as MRD monitoring genes. HCC, hepatocellular carcinoma; MRD, minimal residual disease; RFS, relapse-free survival.

The predictive value of MRD status for early relapse

We further evaluated the predictive value of MRD status for early postoperative relapse in the ZJZS2020 cohort comprising 20 HCC patients who underwent R0 resection. Table 1 presents the baseline characteristics of this cohort. Patients were predominantly male (n=17, 85%), and the median age was 57.5 years (range, 42-76). Most patients presented HBV infection (14/20, 70.0%), Child-Pugh A (19/20, 95.0%), and a single tumor (15/20, 75.0%). About half of the patients had cirrhosis (11, 55.0%) and MVI (12/20, 60.0%) and were at BCLC stage C (10/20, 50.0%). The AFP level of ≥ 400 ng/mL was observed in three patients (15.0%). We identified ≥ 1 mutation in 19 of 20 (95.0%) patients. Based on the definition of MRD status, 20 patients were divided into two subgroups, MRD-positive (6, 30%) and MRD-negative (14, 70%) subgroups. The specific mutational information of 6 MRD-positive patients is listed in Table 2. Patient baseline characteristics were well balanced between subgroups. With a median postoperative follow-up of 14.1 months (range, 5.5-22.1 months), six MRD-positive cases and two of 14 MRD-negative cases had experienced an early relapse. MRD status presented a sensitivity of 75% (6/8) and a specificity of 100% (6/6) in identifying early postoperative relapse. Kaplan-Meier plots indicated that the MRD-positive patients had a significantly poor RFS (mRFS, 4.2 vs. NR months, P=0.002) compared with the MRD-negative patients (Figure 4A). Univariable analysis revealed that MRD status was the only significant variable to predict PFS (HR=13.00, 95% CI 2.60-69.00, P=0.002) (Figure 4B).

TABLE 1

Table 1 Basic characteristics of HCC patients in the ZJZS2020 cohort.

TABLE 2

Table 2 The specific mutational information of six MRD-positive HCC patients.

FIGURE 4

Figure 4 The predictive value of MRD status for early relapse. (A) Kaplan-Meier curves comparing RFS between MRD-positive and -negative patients with HCC. (B) The univariate Cox regression analyses reveal that only MRD status is significantly associated with RFS. HCC, hepatocellular carcinoma; MRD, minimal residual disease; RFS, relapse-free survival.

ctDNA dynamic monitoring

Of 20 patients in the ZJZS2020 cohort, two patients (Patients A and B) underwent longitudinal assessment of MRD status (Figures 5A, B). Patient A is a 76-year-old male with HCV positivity, and his surgical pathological examination demonstrated a grade II-III HCC with MVI and no definite nerve invasion. AFP was restored to the normal level from 7.05ug/L before the operation. The blood collected seven days after surgery was subjected to ctDNA analysis using the 381-gene panel. NF1 mutation with a variant allele fraction (VAF) of 0.0050933 was observed in the plasma sample, suggesting that this patient was MRD-positive. This patient was subsequently followed up, and his AFP level was examined every 1-2 months. On postoperative month 8, the AFP level increased to 11.68 ug/L; ctDNA analysis presented the mutations in three MRD monitoring genes (NF1, TP53, and CTNNB1) with the maximal VAF of 0.0345581; and the imaging evaluation showed no sign of abnormality. On postoperative month 9.6, computed tomography (CT) scanning revealed multiple metastases of both lungs and mediastinal lymph node metastases, and AFP level was increased up to13.97 μg/L, suggesting disease progression. Meanwhile, ctDNA analysis identified mutations in four monitoring genes, including TERT, TP53, CTNNB1, and EZH2, and the maximal VAF of 0.232117 was obviously higher than that of the previous two tests.

FIGURE 5

Figure 5 ctDNA dynamic monitoring. The dynamical changes in the number of MRD monitoring gene variants, the level of AFP/PIVKA-II in plasma samples, and imaging results of HCC patients A (A) and B (B) during HCC progression. HCC, hepatocellular carcinoma; MRD, minimal residual disease; AFP, alpha-foetoprotein; PIVKA-II, a protein induced by vitamin K absence or antagonist II.

Patient B is a 57-year-old male with HBV-positive liver cirrhosis. The postoperative pathology revealed a grade II-III HCC with MVI and hepatic capsule invasion. The analysis of ctDNA from the plasma collected seven days after surgery showed no mutations in the MRD monitoring genes, and the patient was MRD-negative. CREBBP mutation with a VAF of 0.0050933 was detected by ctDNA analysis on postoperative month 6, suggesting that the patient’s MRD status converted from negative to positive, while the imaging evaluation showed no sign of abnormality. The AFP level remained normal after surgery. The PIVKA-11, a potential biomarker complementary to AFP for HCC diagnosis, remained stable in the first six months after surgery and then increased rapidly, reaching 1146.71 mAU/ml in the tenth postoperative month. At 11 months after surgery, the PET-CT imaging results showed changes in the upper part of the right femur and increased metabolism, which was more considered for metastasis and disease recurrence; the ctDNA results still represented MRD-positive. These results showed that serial ctDNA monitoring provides dynamic information on somatic variants with pathogenic effects in tumors and might be superior to single detection to evaluate disease progression.

Discussion

In this study, we profiled genomic alterations by targeted NGS of ctDNA from Chinese HCC patients and identified MRD monitoring genes by filtering out the genes unrelated to early relapse to define MRD status. Evaluation of the prognostic value of MRD status in an independent cohort showed that MRD positivity was an independent prognostic factor for poor RFS. Additionally, the data of two patients undergoing multiple ctDNA tests suggest that the serial longitudinal MRD status tracking may be superior to single testing in predicting postoperative relapse of HCC patients.

Genomic profiles of ctDNA have been shown to be closely associated with types and stages of malignant tumors, suggesting that ctDNA could contribute to cancer diagnosis and personalized treatment (35, 36). The genomic characterization of plasma ctDNA of HCC patients has been reported by several studies, but these studies focus primarily on limited genes or in a small subset. For instance, in Caucasian populations, IKEDA et al. (37) reported the spectrum of genomic mutations in ctDNA from 14 HCC patients by a 68-gene NGS panel, and all patients had somatic alterations, of which the most predominant somatic mutations were in TP53 (57%), followed by CTNNB1 (29%), PTEN (7%), CDKN2A (7%), ARID1A (7%), and MET (7%). In Chinese populations, Cai and colleagues (38) analyzed the mutation profiles of plasma ctDNA from three HCC patients by a 574-gene NGS panel, highlighting that ctDNA may be able to overcome tumor heterogeneity and monitor therapeutic efficacy in real time. Yan et al. (39) utilized a 354-gene NGS panel to profile the genomic landscape of ctDNA from 26 HCC patients and found a positive ctDNA detection rate of 96.2%, wherein TP53 (50.00%) was the most common mutant gene, followed by AXIN1 (11.54%), BCOR (11.54%), CTNNB1 (11.54%), FANCE (11.54%), FANCM (11.54%), and NCOR1 (11.54%). Here we characterized the genomic alteration spectrum of ctDNA from 493 Chinese HCC patients by the 381/733-gene NGS panel. To our knowledge, this is the most comprehensive targeted NGS analysis of ctDNA in Chinese HCC. Consistently, TP53 was the most commonly mutated gene (45%). It is worth noting that the detection rate of ctDNA in the above studies of others and ours is higher than 90%. The high ctDNA detection rate may be attributed to the high vascular nature of HCC or/and the decreased hepatic clearance of ctDNA, because the liver is the main organ responsible for cfDNA clearance (40). Meanwhile, the high ctDNA detection rate also suggested that NGS-based ctDNA detection is technically feasible for identifying genomic variations in patients with HCC. Furthermore, the mutation profiling of ctDNA showed similar mutation frequencies for multiple known HCC drivers, including TP53, TERT, CTNNB1, CDKN2A, ARID1A, RB1, and LRP1B (41–43), compared to those of tissue samples reported by Wang et al. (44), supporting biological plausibility and feasibility of plasma ctDNA testing among patients with HCC.

MRD is considered an important factor generating early postoperative relapse, and ctDNA detection based on the tumor-informed or tumor-agnostic approach has been proven an effective method for detecting MRD in various tumor types (24–27, 30, 45, 46). While the published research on ctDNA detection in HCC is mainly based on tumor-informed approaches. For example, Zhu et al. (20) applied a tumor-informed ctDNA detection approach in HCC patients by a 197-gene panel and identified that postoperative ctDNA positivity was significantly associated with tumor recurrence. Similarly, Cai et al. (14) identified that postoperative ctDNA could indicate MRD more accurately than the conventional protein biomarkers AFP and DCP and was an independent prognostic factor for both RFS and overall survival. These findings show that tumor-informed ctDNA testing is feasible for monitoring MRD, but its broad application is inevitably constrained by the inherent limitations of the strategy, such as the inability to detect emerging mutations, long and much labor-intensive processing, and the resulting high cost. Here we tried to construct a fixed panel to detect MRD. The MRD monitoring genes that made up this fixed panel were identified by removing the genes with a mutation frequency of ≥ 10% in ctDNA of the ZJ2020 cohort with patients having RFS over 2 years after radical resection from the genes with a mutation frequency of ≥ 5% in ctDNA of Chinese HCC cohort. A total of 13 genes were identified as MRD monitoring genes. Among these 13 genes, TP53, TERT, CTNNB1, APC, and CDKN2A were well-defined driver genes for HCC (41–43, 47). The remaining eight genes have also been experimentally verified to be associated with HCC. RBM10, a member of the RNA binding motif gene family, plays a regulatory role in alternative splicing. As a tumor suppressor, its downregulation has been linked to tumor progression, metastasis, and poor prognosis in multiple human cancers, including HCC (48–50). NTRK3, located on chromosome 15q25, has been wildly reported as a tumor suppressor implicated in the modulation of cell growth, invasion, and migration in a diverse array of tumors, including HCC (51–53). NOTCH1 and NOTCH2 are the key regulators of stem cell proliferation, differentiation, and apoptosis. Numerous studies have revealed that NOTCH1 and NOTCH2 are involved in the development of HCC, and their activation contributes to HCC cell growth and aggressiveness and poor overall survival of HCC patients (54–57). NF1, one of the largest human genes, is located on chromosome 17, band q11.2.6. Recently, Lu et al. (58) identified NF1 as a critical driver for lenvatinib resistance in HCC, whose loss reactivates the PI3K/AKT and MAPK/ERK signaling pathways. CREBBP, encoding an acetyltransferase, is one of the most frequently mutated genes in small cell lung cancer. In the last few years, some studies documented that CREBBP activation is associated with the early recurrence of HCC, and targeting CREBBP attenuates HCC progression (59, 60). GLI3, a member of the Hedgehog signaling pathway, is upregulated in a variety of tumor types (61). The positive correlation of Gli3 with tumor progression has been observed in HCC (62), pancreatic cancer (63), colon cancer (64), ovarian cancer (65), breast cancer (65), and bladder cancer (66). EZH2 is the catalytic subunit of polycomb repressive complex 2, which methylates histone H3 lysine 27, silencing a number of tumor-suppressor genes, including E-cadherin. Available data have supported EZH2 involved in metastatic spread and tumor angiogenesis. Several recent studies not only linked this role of EZH2 to HCC development but also bridged its pro-oncogenic function to the decrease of PD-L1 and shaping of tumor immunosuppressive microenvironment, suggesting that HCC patients with EZH2 mutations should be treated with immunotherapy carefully (67–69). MRD positivity was defined as ≥ one MRD monitoring gene mutation in this fixed panel.

We further evaluated the predictive value of MRD status for post-surgery relapse in an independent HCC cohort and showed that patients who were MRD-positive postoperatively had a significantly high risk of early relapse versus those with negative MRD. The prognostic value of MRD positivity has been well-established as an independent prognostic factor for poor RFS. MRD status’s sensitivity and specificity to predict early postoperative relapse were 75% and 100%, respectively. These results showed that this panel provides a cost-effective and feasible approach to monitoring MRD and predicting early postoperative relapse. However, our data also indicated that in some cases, MRD detection at a single postoperative time point was insufficient to predict the prognosis of HCC patients, which might result from inadequate ctDNA release from MRD right after curative resection. Serial MRD monitoring could overcome this limitation and provide a full picture of MRD dynamics during the course of HCC.

Our study has several limitations. First, this is a retrospective study, and the screening strategy for MRD monitoring genes is developed by ourselves, which is reasonable but may not be strict enough. Nevertheless, our results are still beneficial to promoting a standardized sequencing process and a rigorous strategy for identifying MRD monitoring genes. Second, the cohort for MRD performance evaluation is relatively small, and serial ctDNA sampling is in a few patients. Further prospective studies focusing on large serial samplings will be needed to evaluate the performance of this fixed panel in detail.

Conclusion

In summary, this study demonstrates the feasibility of MRD evaluation based on a tumor-agnostic liquid biopsy approach for HCC patients and reports promising data. The established 13-gene panel can reliably predict early relapse after radical hepatectomy, contributing to the personalized management of HCC patients.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author/s.

Ethics statement

The study was approved by the ethical institutional review board of Zhujiang Hospital and performed following the guidelines of the Declaration of Helsinki. Written informed consent was obtained from individual participants.

Author contributions

Conception and design: MP, YC, YX, TC, and KZ. Provision of study materials or patients: MP. Collection and assembly of data: YX, TC, KZ, YW, LC, GH, HL, CZ, and SF. Data analysis and interpretation: YX, JiaPC, KZ, TC, JinPC, XZ, CC, MH, and KZ. Drafting of the manuscript: JinPC, TC, YX, JiaPC, and KZ. Review and/or revision of the manuscript: All authors. Accountable for all aspects of the work: All authors. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by a grant from the National Natural Science Foundation of China (No.82072627) and the Guangdong Basic and Applied Basic Research Foundation of China (No. 2021B1515230011).

Acknowledgments

The authors would like to acknowledge the technical support provided by Hao Chen from 3D Medicines, Inc.

Conflict of interest

TC, JinPC, XZ, CC, and MH are employees of 3D Medicines Inc.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1119744/full#supplementary-material

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Villanueva A. Hepatocellular carcinoma. N Engl J Med (2019) 380(15):1450–62. doi: 10.1056/NEJMra1713263

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Llovet JM, Kelley RK, Villanueva A, Singal AG, Pikarsky E, Roayaie S, et al. Hepatocellular carcinoma. Nat Rev Dis Primers (2021) 7(1):6. doi: 10.1038/s41572-020-00240-3

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Llovet JM, Kelley RK, Villanueva A, Singal AG, Pikarsky E, Roayaie S, et alHepatocellular carcinoma. Nat Rev Dis Primers (2021) 7(1):7. doi: 10.1038/s41572-021-00245-6

PubMed Abstract | CrossRef Full Text | Google Scholar

5. W. H. Organization. International agency for research on cancer. (Geneva, Switzerland: World Health Organization) (2019).

Google Scholar

6. Sulpice L, Rayar M, Boucher E, Pracht M, Meunier B, Boudjema K. Treatment of recurrent intrahepatic cholangiocarcinoma. J Br Surg (2012) 99(12):1711–7. doi: 10.1002/bjs.8953

CrossRef Full Text | Google Scholar

7. Hyder O, Hatzaras I, Sotiropoulos GC, Paul A, Alexandrescu S, Marques H, et al. Recurrence after operative management of intrahepatic cholangiocarcinoma. Surgery (2013) 153(6):811–8. doi: 10.1016/j.surg.2012.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zhang X-F, Beal E, Bagante F, Chakedis J, Weiss M, Popescu I, et al. Early versus late recurrence of intrahepatic cholangiocarcinoma after resection with curative intent. J Br Surg (2018) 105(7):848–56. doi: 10.1002/bjs.10676

CrossRef Full Text | Google Scholar

9. Imamura H, Matsuyama Y, Tanaka E, Ohkubo T, Hasegawa K, Miyagawa S, et al. Risk factors contributing to early and late phase intrahepatic recurrence of hepatocellular carcinoma after hepatectomy. J Hepatol (2003) 38(2):200–7. doi: 10.1016/S0168-8278(02)00360-4

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Sun Y, Wu L, Zhong Y, Zhou K, Hou Y, Wang Z, et al. Single-cell landscape of the ecosystem in early-relapse hepatocellular carcinoma. Cell (2021) 184(2):404–421 e16. doi: 10.1016/j.cell.2020.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sherman M. Recurrence of hepatocellular carcinoma. N Engl J Med (2008) 359(19):2045–7. doi: 10.1056/NEJMe0807581

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Zheng J, Kuk D, Gonen M, Balachandran VP, Kingham TP, Allen PJ, et al. Actual 10-year survivors after resection of hepatocellular carcinoma. Ann Surg Oncol (2017) 24(5):1358–66. doi: 10.1245/s10434-016-5713-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kim JH, Sohn BH, Lee H-S, Kim S-B, Yoo JE, Park Y-Y, et al. Genomic predictors for recurrence patterns of hepatocellular carcinoma: model derivation and validation. PloS Med (2014) 11(12):e1001770. doi: 10.1371/journal.pmed.1001770

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Cai Z, Chen G, Zeng Y, Dong X, Li Z, Huang Y, et al. Comprehensive liquid profiling of circulating tumor DNA and protein biomarkers in long-term follow-up patients with hepatocellular CarcinomaComprehensive ctDNA profile integration to monitor HCC. Clin Cancer Res (2019) 25(17):5284–94. doi: 10.1158/1078-0432.CCR-18-3477

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Poon RTP, Fan ST, Ng IOL, Lo CM, Liu CL, Wong J. Different risk factors and prognosis for early and late intrahepatic recurrence after resection of hepatocellular carcinoma. Cancer (2000) 89(3):500–7. doi: 10.1002/1097-0142(20000801)89:3<500::AID-CNCR4>3.0.CO;2-O

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Wu J-C, Huang Y-H, Chau G-Y, Su C-W, Lai C-R, Lee P-C, et al. Risk factors for early and late recurrence in hepatitis b-related hepatocellular carcinoma. J Hepatol (2009) 51(5):890–7. doi: 10.1016/j.jhep.2009.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Sohn W, Paik Y-H, Kim JM, Kwon CH, Joh JW, Cho JY, et al. HBV DNA and HBsAg levels as risk predictors of early and late recurrence after curative resection of HBV-related hepatocellular carcinoma. Ann Surg Oncol (2014) 21(7):2429–35. doi: 10.1245/s10434-014-3621-x

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Portolani N, Coniglio A, Ghidoni S, Giovanelli M, Benetti A, Tiberio GAM, et al. Early and late recurrence after liver resection for hepatocellular carcinoma: prognostic and therapeutic implications. Ann Surg (2006) 243(2):229. doi: 10.1097/01.sla.0000197706.21803.a1

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Choi G-H, Kim D-H, Kang C-M, Kim K-S, Choi J-S, Lee W-J, et al. Prognostic factors and optimal treatment strategy for intrahepatic nodular recurrence after curative resection of hepatocellular carcinoma. Ann Surg Oncol (2008) 15(2):618–29. doi: 10.1245/s10434-007-9671-6

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Zhu GQ, Liu WR, Tang Z, Qu WF, Fang Y, Jiang XF, et al. Serial circulating tumor DNA to predict early recurrence in patients with hepatocellular carcinoma: a prospective study. Mol Oncol (2022) 16(2):549–61. doi: 10.1002/1878-0261.13105

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Garcia-Murillas I, Schiavon G, Weigelt B, Ng C, Hrebien S, Cutts RJ, et al. Mutation tracking in circulating tumor DNA predicts relapse in early breast cancer. Sci Trans Med (2015) 7(302):302ra133–302ra133. doi: 10.1126/scitranslmed.aab0021

CrossRef Full Text | Google Scholar

22. Tarazona N, Gimeno-Valiente F, Gambardella V, Zuñiga S, Rentero-Garrido P, Huerta M, et al. Targeted next-generation sequencing of circulating-tumor DNA for tracking minimal residual disease in localized colon cancer. Ann Oncol (2019) 30(11):1804–12. doi: 10.1093/annonc/mdz390

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Chen G, Peng J, Xiao Q, Wu H-X, Wu X, Wang F, et al. Postoperative circulating tumor DNA as markers of recurrence risk in stages II to III colorectal cancer. J Hematol Oncol (2021) 14(1):1–11. doi: 10.1186/s13045-021-01089-z

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Tie J, Wang Y, Tomasetti C, Li L, Springer S, Kinde I, et al. Circulating tumor DNA analysis detects minimal residual disease and predicts recurrence in patients with stage II colon cancer. Sci Trans Med (2016) 8(346):346ra92–346ra92. doi: 10.1126/scitranslmed.aaf6219

CrossRef Full Text | Google Scholar

25. Chaudhuri AA, Chabon JJ, Lovejoy AF, Newman AM, Stehr H, Azad TD, et al. Early detection of molecular residual disease in localized lung cancer by circulating tumor DNA ProfilingEarly detection of lung cancer MRD by ctDNA profiling. Cancer Discov (2017) 7(12):1394–403. doi: 10.1158/2159-8290.CD-17-0716

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Punnoose EA, Atwal S, Liu W, Raja R, Fine BM, Hughes BG, et al. Evaluation of circulating tumor cells and circulating tumor DNA in non–small cell lung cancer: Association with clinical endpoints in a phase II clinical trial of pertuzumab and ErlotinibCTCs and HER family–targeted therapies in NSCLC. Clin Cancer Res (2012) 18(8):2391–401. doi: 10.1158/1078-0432.CCR-11-3148

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Takai E, Yachida S. Circulating tumor DNA as a liquid biopsy target for detection of pancreatic cancer. World J Gastroenterol (2016) 22(38):8480. doi: 10.3748/wjg.v22.i38.8480

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Gao Q, Wang ZC, Duan M, Lin YH, Zhou XY, Worthley DL, et al. Cell culture system for analysis of genetic heterogeneity within hepatocellular carcinomas and response to pharmacologic agents. Gastroenterology (2017) 152(1):232–42 e4. doi: 10.1053/j.gastro.2016.09.008

PubMed Abstract | CrossRef Full Text | Google Scholar

29. von Felden J, Craig AJ, Garcia-Lezana T, Labgaa I, Haber PK, D’Avola D, et al. Mutations in circulating tumor DNA predict primary resistance to systemic therapies in advanced hepatocellular carcinoma. Oncogene (2021) 40(1):140–51. doi: 10.1038/s41388-020-01519-1

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Parikh AR, Van Seventer EE, Siravegna G, Hartwig AV, Jaimovich A, He Y, et al. Minimal residual disease detection using a plasma-only circulating tumor DNA assay in patients with colorectal CancerPlasma-only ctDNA-guided MRD detection in patients with CRC. Clin Cancer Res (2021) 27(20):5586–94. doi: 10.1158/1078-0432.CCR-21-0410

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Mayakonda A, Lin D-C, Assenov Y, Plass C, Koeffler HP. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Evans EJ Jr., DeGregori J. Cells with cancer-associated mutations overtake our tissues as we age. Aging Cancer (2021) 2(3):82–97. doi: 10.1002/aac2.12037

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Azad TD, Chaudhuri AA, Fang P, Qiao Y, Esfahani MS, Chabon JJ, et al. Circulating tumor DNA analysis for detection of minimal residual disease after chemoradiotherapy for localized esophageal cancer. Gastroenterology (2020) 158(3):494–505.e6. doi: 10.1053/j.gastro.2019.10.039

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Tie J, Wang Y, Cohen J, Li L, Hong W, Christie M, et al. Circulating tumor DNA dynamics and recurrence risk in patients undergoing curative intent resection of colorectal cancer liver metastases: A prospective cohort study. PloS Med (2021) 18(5):e1003620. doi: 10.1371/journal.pmed.1003620

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Schwaederle M, Chattopadhyay R, Kato S, Fanta PT, Banks KC, Choi IS, et al. Genomic alterations in circulating tumor DNA from diverse cancer patients identified by next-generation SequencingGenomic alterations analysis by NGS of ctDNA. Cancer Res (2017) 77(19):5419–27. doi: 10.1158/0008-5472.CAN-17-0885

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Diaz LA Jr., Bardelli A. Liquid biopsies: genotyping circulating tumor DNA. J Clin Oncol (2014) 32(6):579. doi: 10.1200/JCO.2012.45.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Ikeda S, Tsigelny IF, Skjevik Å A, Kono Y, Mendler M, Kuo A, et al. Next-generation sequencing of circulating tumor DNA reveals frequent alterations in advanced hepatocellular carcinoma. Oncologist (2018) 23(5):586–93. doi: 10.1634/theoncologist.2017-0479

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Cai ZX, Chen G, Zeng YY, Dong XQ, Lin MJ, Huang XH, et al. Circulating tumor DNA profiling reveals clonal evolution and real-time disease progression in advanced hepatocellular carcinoma. Int J Cancer (2017) 141(5):977–85. doi: 10.1002/ijc.30798

PubMed Abstract | CrossRef Full Text | Google Scholar

39. An Y, Guan Y, Xu Y, Han Y, Wu C, Bao C, et al. The diagnostic and prognostic usage of circulating tumor DNA in operable hepatocellular carcinoma. Am J Trans Res (2019) 11(10):6462.

Google Scholar

40. Elshimali YI, Khaddour H, Sarkissyan M, Wu Y, Vadgama JV. The clinical utilization of circulating cell free DNA (CCFDNA) in blood of cancer patients. Int J Mol Sci (2013) 14(9):18925–58. doi: 10.3390/ijms140918925

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Fujimoto A, Totoki Y, Abe T, Boroevich KA, Hosoda F, Nguyen HH, et al. Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators. Nat Genet (2012) 44(7):760–4. doi: 10.1038/ng.2291

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Bailey MH, Tokheim C, Porta-Pardo E, Sengupta S, Bertrand D, Weerasinghe A, et al. Comprehensive characterization of cancer driver genes and mutations. Cell (2018) 173(2):371–385. e18. doi: 10.1016/j.cell.2018.02.060

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Guichard C, Amaddeo G, Imbeaud S, Ladeiro Y, Pelletier L, Maad IB, et al. Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma. Nat Genet (2012) 44(6):694–8. doi: 10.1038/ng.2256

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Wang S, Shi H, Liu T, Li M, Zhou S, Qiu X, et al. Mutation profile and its correlation with clinicopathology in Chinese hepatocellular carcinoma patients. Hepatobil Surg Nutr (2021) 10(2):172–9. doi: 10.21037/hbsn.2019.09.17

CrossRef Full Text | Google Scholar

45. Janni W, Huober J, Braun T, Müller V, Fink A, de Gregorio A, et al. Multiomic, plasma-only circulating tumor DNA (ctDNA) assay identifies breast cancer patients with minimal residual disease (MRD) and predicts distant recurrence. Cancer Res (2022) 82(12_Supplement):3403–3. doi: 10.1158/1538-7445.AM2022-3403

CrossRef Full Text | Google Scholar

46. Dawson S-J, Tsui DW, Murtaza M, Biggs H, Rueda OM, Chin S-F, et al. Analysis of circulating tumor DNA to monitor metastatic breast cancer. N Engl J Med (2013) 368(13):1199–209. doi: 10.1056/NEJMoa1213261

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Zhang Z, Xu L, Sun C. Comprehensive characterization of cancer genes in hepatocellular carcinoma genomes. Oncol Lett (2018) 15(2):1503–10. doi: 10.3892/ol.2017.7521

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zhao Z, Li J, Shen F. Protective effect of the RNA-binding protein RBM10 in hepatocellular carcinoma. Eur Rev Med Pharmacol Sci (2020) 24(11):6005–13.

PubMed Abstract | Google Scholar

49. Jiang L, Zhang K, Hu C, Zhang H, Wang D, Zhai W. Effect of RBM10 on regulating the proliferation and metastasis activity of human hepatocellular carcinoma cells by affecting the stability of miR-21. Biocell (2022) 46(4):969. doi: 10.32604/biocell.2022.017765

CrossRef Full Text | Google Scholar

50. Pang S-J, Sun Z, Lu W-F, Si-Ma H, Lin Z-P, Shi Y, et al. Integrated bioinformatics analysis and validation of the prognostic value of RBM10 expression in hepatocellular carcinoma. Cancer Manage Res (2022) 14:969. doi: 10.2147/CMAR.S349884

CrossRef Full Text | Google Scholar

51. Xiong D, Sheng Y, Ding S, Chen J, Tan X, Zeng T, et al. LINC00052 regulates the expression of NTRK3 by miR-128 and miR-485-3p to strengthen HCC cells invasion and migration. Oncotarget (2016) 7(30):47593. doi: 10.18632/oncotarget.10250

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Luo Y, Kaz AM, Kanngurn S, Welsch P, Morris SM, Wang J, et al. NTRK3 is a potential tumor suppressor gene commonly inactivated by epigenetic mechanisms in colorectal cancer. PloS Genet (2013) 9(7):e1003552. doi: 10.1371/journal.pgen.1003552

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Zhang H-F, Li W, Han Y-D. LINC00261 suppresses cell proliferation, invasion and notch signaling pathway in hepatocellular carcinoma. Cancer Biomarkers (2018) 21(3):575–82. doi: 10.3233/CBM-170471

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Gao J, Song Z, Chen Y, Xia L, Wang J, Fan R, et al. Deregulated expression of notch receptors in human hepatocellular carcinoma. Digest Liver Dis (2008) 40(2):114–21. doi: 10.1016/j.dld.2007.08.001

CrossRef Full Text | Google Scholar

55. Huang Q, Li J, Zheng J, Wei A. The carcinogenic role of the notch signaling pathway in the development of hepatocellular carcinoma. J Cancer (2019) 10(6):1570. doi: 10.7150/jca.26847

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Hayashi Y, Osanai M, Lee G-H. NOTCH2 signaling confers immature morphology and aggressiveness in human hepatocellular carcinoma cells. Oncol Rep (2015) 34(4):1650–8. doi: 10.3892/or.2015.4171

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Zhou L, Zhang N, Song W, You N, Li Q, Sun W, et al. The significance of Notch1 compared with Notch3 in high metastasis and poor overall survival in hepatocellular carcinoma. PloS One (2013) 8(2):e57382. doi: 10.1371/journal.pone.0057382

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Lu Y, Shen H, Huang W, He S, Chen J, Zhang D, et al. Genome-scale CRISPR-Cas9 knockout screening in hepatocellular carcinoma with lenvatinib resistance. Cell Death Discov (2021) 7(1):1–12. doi: 10.1038/s41420-021-00747-y

CrossRef Full Text | Google Scholar

59. Cai L-Y, Chen S-J, Xiao S-H, Sun Q-J, Ding C-H, Zheng B-N, et al. Targeting p300/CBP attenuates hepatocellular carcinoma progression through epigenetic regulation of metabolism. Cancer Res (2021) 81(4):860–72. doi: 10.1158/0008-5472.CAN-20-1323

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Cui Y, Wu X, Lin C, Zhang X, Ye L, Ren L, et al. AKIP1 promotes early recurrence of hepatocellular carcinoma through activating the wnt/β-catenin/CBP signaling pathway. Oncogene (2019) 38(27):5516–29. doi: 10.1038/s41388-019-0807-5

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Matissek SJ, Elsawa SF. GLI3: A mediator of genetic diseases, development and cancer. Cell Commun Signal (2020) 18(1):54. doi: 10.1186/s12964-020-00540-x

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Hyun J, Wang S, Kim J, Rao KM, Park SY, Chung I, et al. MicroRNA-378 limits activation of hepatic stellate cells and liver fibrosis by suppressing Gli3 expression. Nat Commun (2016) 7(1):1–16. doi: 10.1038/ncomms10993

CrossRef Full Text | Google Scholar

63. Steg A, Amm HM, Novak Z, Frost AR, Johnson MR. Gli3 mediates cell survival and sensitivity to cyclopamine in pancreatic cancer. Cancer Biol Ther (2010) 10(9):893–902. doi: 10.4161/cbt.10.9.13252

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Trnski D, Sabol M, Gojević A, Martinić M, Ozretić P, Musani V, et al. GSK3β and Gli3 play a role in activation of hedgehog-gli pathway in human colon cancer–targeting GSK3β downregulates the signaling pathway and reduces cell proliferation. Biochim Biophys Acta (BBA)-Molecular Basis Dis (2015) 1852(12):2574–84. doi: 10.1016/j.bbadis.2015.09.005

CrossRef Full Text | Google Scholar

65. Lin M, Zhu H, Shen Q, Sun L-Z, Zhu X. GLI3 and androgen receptor are mutually dependent for their malignancy-promoting activity in ovarian and breast cancer cells. Cell Signalling (2022) 92:110278. doi: 10.1016/j.cellsig.2022.110278

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Li J, Qiu M, An Y, Huang J, Gong C. miR-7-5p acts as a tumor suppressor in bladder cancer by regulating the hedgehog pathway factor Gli3. Biochem Biophys Res Commun (2018) 503(3):2101–7. doi: 10.1016/j.bbrc.2018.07.166

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Sudo T, Utsunomiya T, Mimori K, Nagahara H, Ogawa K, Inoue H, et al. Clinicopathological significance of EZH2 mRNA expression in patients with hepatocellular carcinoma. Br J Cancer (2005) 92(9):1754–8. doi: 10.1038/sj.bjc.6602531

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Xiao G, Jin L-L, Liu C-Q, Wang Y-C, Meng Y-M, Zhou Z-G, et al. EZH2 negatively regulates PD-L1 expression in hepatocellular carcinoma. J Immunother Cancer (2019) 7(1):1–15. doi: 10.1186/s40425-019-0784-9

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Guo B, Tan X, Cen H. EZH2 is a negative prognostic biomarker associated with immunosuppression in hepatocellular carcinoma. PloS One (2020) 15(11):e0242191. doi: 10.1371/journal.pone.0242191

PubMed Abstract | CrossRef Full Text | Google Scholar