Genetic variants in the EPCAM gene is associated with the prognosis of transarterial chemoembolization treated hepatocellular carcinoma with portal vein tumor thrombus

Useful genes for predicting the efficacy of transarterial chemoembolization in hepatocellular carcinoma

Transarterial chemoembolization (TACE) is generally used to treat patients with hepatocellular carcinoma (HCC), a common and deadly cancer; however, its efficacy varies according to factors such as tumor volume, stage, serum alpha-fetoprotein level, and chosen feeding artery. In addition, gene-related factors have been recently suggested to be involved in the regulation and prediction of TACE outcomes. Accordingly, genes could serve as effective biomarkers to select patients who can benefit from TACE. These gene-related factors can activate signaling pathways affecting cancer cell survival while regulating the epithelial-mesenchymal transition, angiogenesis, and the tumor microenvironment, all directly associated with tumor progression, thereby affecting TACE efficacy. Moreover, this disordered gene expression is associated with poor prognosis in patients with HCC, including TACE resistance, postoperative recurrence, and metastasis. To identify the exact relationship between various genes and TACE efficacy, this review summarizes the involvement of protein-coding and non-coding genes and single nucleotide polymorphisms in TACE efficacy for predicting the efficacy of TACE; the present findings may help improve the efficacy of TACE in clinical settings.

rs62139665 Polymorphism in the Promoter Region of EpCAM Is Associated With Hepatitis C Virus-Related Hepatocellular Carcinoma Risk in Egyptians

Hepatocellular carcinoma (HCC) is a universal health problem that is particularly alarming in Egypt. The major risk factor for HCC is hepatitis C virus (HCV) infection which is a main burden in Egypt. The epithelial cell adhesion molecule (EpCAM) is a stem cell marker involved in the tumorigenesis and progression of many malignancies, including HCC. We investigated the association of -935 C/G single nucleotide polymorphism in EpCAM promoter region (rs62139665) with HCC risk, EpCAM expression and overall survival in Egyptians. A total of 266 patients (128 HCV and 138 HCC cases) and 117 age- and sex-matched controls participated in this study. Genotyping, performed using allelic discrimination and confirmed by sequencing, revealed a significant association between EpCAM rs62139665 and HCC susceptibility, with higher GG genotype and G allele distribution in HCC patients than in non-HCC subjects. Such association was not detected in HCV patients compared to controls. EpCAM gene expression levels, determined in blood by RT-qPCR, and its serum protein expression levels, determined by ELISA, were significantly higher in GG relative to GC+CC genotype carriers in HCV and HCC patients in a recessive model. ROC analysis of EpCAM protein levels revealed significant discriminatory power between HCC patients and non-HCC subjects, with improved diagnostic accuracy when combining α-fetoprotein and EpCAM compared to that of α-fetoprotein alone. Altogether, EpCAM rs62139665 polymorphism is significantly associated with HCC and with EpCAM gene and protein expression levels in the Egyptian population. Moreover, serum EpCAM levels may hold promise for HCC diagnosis and for improving the diagnostic accuracy of α-fetoprotein.

Integrative Functional Genomics Implicated the Key T-/B-Cell Deficiency Regulator RAG1 in Transarterial Chemoembolization of Hepatocellular Carcinoma

Transarterial chemoembolization (TACE) has significantly prolonged overall survival (OS) of unresectable hepatocellular carcinoma (HCC) patients. Unfortunately, there are still a portion of patients without therapeutic responses to TACE. Although genome-wide association studies identified multiple HCC susceptibility SNPs, it is still largely unclear how genome-wide identified functional SNPs impacting gene expression contribute to the prognosis of TACE-treated HCC patients. In this study, we developed an integrative functional genomics methodology to identify gene expression-related SNPs significantly contributing to prognosis of TACE-treated HCC patients across the whole genome. Employing integration of data from expression quantitative trait locus (eQTLs) analyses of The Cancer Genome Atlas (TCGA) liver hepatocellular carcinoma (LIHC) as well as the 1000 Genomes project, we successfully annotated 60 gene expression-related SNPs which are associated with OS of the TCGA patients. After genotyping these 60 SNPs in our TACE cohort, we identified four SNPs (rs12574873, rs12513391, rs34597395, and rs35624901) which are significantly associated with OS of HCC patients treated with TACE. For instance, multivariate Cox proportional hazards model indicated that the rs35624901 Deletion.Deletion (Del.Del) genotype carriers had markedly prolonged OS and a 55% decreased death risk compared with individuals with the GG genotype after TACE therapy (p = 8.3 × 10^–5). In support of this, the rs35624901 Del.Del genotype is correlated to higher expression of RAG1, a key T-/B-cell deficiency regulator. Our findings reported the first evidence supporting the prognostic value of four eQTL SNPs in TACE-treated HCC patients. Importantly, our data implicated that antitumor immunity might contribute to TACE efficiency for unresectable HCC patients.

Gene Signature Associated with Upregulation of the Wnt/β-Catenin Signaling Pathway Predicts Tumor Response to Transarterial Embolization

PURPOSE

To identify gene mutations in tumors undergoing transarterial embolization and explore the relationship between gene mutations and tumor response to embolization.

MATERIALS AND METHODS

This was a retrospective review that included 17 patients with primary or metastatic liver tumors treated with embolization and had specimens analyzed for a 341-gene panel next-generation sequence assay. Pathologic conditions included hepatocellular, carcinoid, pancreatic neuroendocrine, melanoma, medullary thyroid, and liver acinar-cell carcinoma. Disease, procedure data, and tumor response data were collected. Dimensionality reduction was performed by using principal component analysis. A linear support vector machine was used to learn a prediction rule and identify the genes most predictive of objective tumor response (partial or complete) per modified Response Evaluation Criteria In Solid Tumors. Cross-validation was used to test the prediction on the holdout set. Permutation testing was used to determine statistical significance of prediction accuracy. Recursive feature elimination was used to identify the most predictive genes.

RESULTS

At 4 months after embolization, 9 tumors showed a response and 8 did not. Using the top two principal components, prediction accuracy of the gene mutation signature was 70% (±11%), which was statistically significant (P < .05). The most predictive genes were CTNNB1, MEN1, and NCOR1: three genes associated with the Wnt/β-catenin and hypoxia signaling pathways.

CONCLUSIONS

This study identifies gene mutations in tumors treated with transarterial embolization. A gene-mutation signature obtained from the mutation data suggests that upregulation of the Wnt/β-catenin signaling pathway may be associated with sensitivity to embolization.

Hepatic transcriptomic alterations for N,N-dimethyl-p-toluidine (DMPT) and p-toluidine after 5-day exposure in rats

N,N-dimethyl-p-toluidine (DMPT), an accelerant for methyl methacrylate monomers in medical devices, was a liver carcinogen in male and female F344/N rats and B6C3F1 mice in a 2-year oral exposure study. p-Toluidine, a structurally related chemical, was a liver carcinogen in mice but not in rats in an 18-month feed exposure study. In this current study, liver transcriptomic data were used to characterize mechanisms in DMPT and p-toluidine liver toxicity and for conducting benchmark dose (BMD) analysis. Male F344/N rats were exposed orally to DMPT or p-toluidine (0, 1, 6, 20, 60 or 120 mg/kg/day) for 5 days. The liver was examined for lesions and transcriptomic alterations. Both chemicals caused mild hepatic toxicity at 60 and 120 mg/kg and dose-related transcriptomic alterations in the liver. There were 511 liver transcripts differentially expressed for DMPT and 354 for p-toluidine at 120 mg/kg/day (false discovery rate threshold of 5 %). The liver transcriptomic alterations were characteristic of an anti-oxidative damage response (activation of the Nrf2 pathway) and hepatic toxicity. The top cellular processes in gene ontology (GO) categories altered in livers exposed to DMPT or p-toluidine were used for BMD calculations. The lower confidence bound benchmark doses for these chemicals were 2 mg/kg/day for DMPT and 7 mg/kg/day for p-toluidine. These studies show the promise of using 5-day target organ transcriptomic data to identify chemical-induced molecular changes that can serve as markers for preliminary toxicity risk assessment.

Alterations of telomere length and mtDNA copy number are associated with overall survival in hepatocellular carcinoma patients treated with transarterial chemoembolization

PURPOSE

Increasing evidence suggests that alterations in mitochondrial DNA (mtDNA) copy number (mtDNAcn) and relative telomere length (RTL) may be implicated in the tumorigenesis of several malignancies. Alterations of both RTL and mtDNAcn are generally accepted as independent biomarkers for predicting risk and prognosis in various cancers. The aim of this study was to evaluate the prognostic value of combining leukocyte RTL with mtDNAcn (RTL-mtDNAcn) in hepatocellular carcinoma (HCC).

METHODS

RTL and mtDNAcn in peripheral blood leukocytes (PBLs) were measured using a real-time PCR-based method in a total of 250 HCC patients treated with transcatheter arterial chemoembolization (TACE). We evaluated the associations between RTL and/or mtDNAcn and HCC overall survival using Kaplan-Meier curve analysis and Cox proportional hazards regression model.

RESULTS

We found that patients with longer leukocyte RTL or lower mtDNAcn had shorter overall survival time. The univariate analysis (HR 1.63, 95 % CI 1.23-2.17, P = 7.7 × 10(-4)) and multivariate analysis (HR 1.78, 95 % CI 1.31-2.42, P = 2.4 × 10(-4)) indicated that longer leukocyte RTL was significantly associated with poorer OS in HCC patients. Kaplan-Meier curve analysis showed that patients with longer RTL had shorter overall survival time than those with shorter RTL (log-rank P = 0.001). Patients with lower mtDNA copy number was significantly associated with poorer OS by Cox proportional hazards model using both univariate (HR 1.60, 95 % CI 1.21-2.13, P = 0.001) and multivariate analyses (HR 1.77, 95 % CI 1.30-2.41, P = 2.8 × 10(-4)). Kaplan-Meier curve analysis showed that patients with lower mtDNA content had significantly shorter overall survival time than those with higher mtDNA content (log-rank P = 0.001). Furthermore, combination of leukocyte RTL and mtDNAcn significantly improved the efficacy of predicting HCC prognosis. Patients with longer RTL and lower mtDNAcn exhibited a significantly poorer overall survival in both the univariate analysis (HR 2.21, 95 % CI 1.52-3.22, P = 3.5 × 10(-5)) and multivariate analysis (HR 2.60, 95 % CI 1.73-3.90, P = 4.3 × 10(-6)). The effect on patient prognosis was more evident in patients with longer RTL and lower mtDNAcn than in those with shorter RTL and lower mtDNA (HR 2.11, 95 % CI 1.34-3.32, P = 0.001) or in those with longer RTL and higher mtDNA (HR 2.10, 95 % CI 1.34-3.27, P = 0.001).

CONCLUSIONS

Our data suggest that combination of leukocyte RTL-mtDNAcn may be a potential efficient prognostic marker for HCC patients receiving the TACE treatment.

Polymorphisms in Genes of Tricarboxylic Acid Cycle Key Enzymes Are Associated with Early Recurrence of Hepatocellular Carcinoma

Alterations of activity and expression in tricarboxylic acid (TCA) cycle key enzymes have been indicated in several malignancies, including hepatocellular carcinoma (HCC). They play an important role in the progression of cancer. However, the impact of single nucleotide polymorphisms (SNPs) in genes encoding these key enzymes on the recurrence of HCC has not been investigated. In this study, we genotyped 17 SNPs in genes encoding TCA cycle key enzymes and analyzed their association with recurrence-free survival (RFS) in a cohort of 492 Chinese HCC patients by Cox proportional hazard model and survival tree analysis. We identified 7 SNPs in SDHC, SDHD, FH, and IDH2 genes to be significantly associated with the RFS of HCC patients. Moreover, all these SNPs were associated with the early recurrence (within 2 years after surgery) risk of diseases. Cumulative effect analysis showed that these SNPs exhibited a dose-dependent effect on the overall and early recurrence. Further stratified analysis suggested that number of risk genotypes modified the protective effect on HCC recurrence conferred by transcatheter arterial chemoembolization treatment. Finally, the survival tree analysis revealed that SNP rs10789859 in SDHD gene was the primary factor contributing to HCC recurrence in our population. To the best of our knowledge, we for the first time observed the association between SNPs in genes encoding TCA cycle key enzymes and HCC recurrence risk. Further observational and functional studies are needed to validate our findings and generalize its clinical usage.

Gene expression in hepatocellular carcinoma: pilot study of potential transarterial chemoembolization response biomarkers

PURPOSE

To perform a feasibility study to explore the relationship between hepatocellular carcinoma genetics and transarterial chemoembolization treatment response to identify potential biomarkers associated with enhanced treatment efficacy.

MATERIALS AND METHODS

In this single-institution study, pretreatment hepatocellular carcinoma biopsy specimens for tumors in 19 patients (14 men, five women; mean age, 59 y) treated with chemoembolization between 2007 and 2013 were analyzed for a panel of 60 chemotherapy-sensitivity, hypoxia, mitosis, and inflammatory genes with the QuantiGene Plex 2.0 mRNA detection assay. Demographic, disease, and procedure data and tumor response outcomes were collected. Quantitative mRNA levels were compared based on radiologic response between tumors exhibiting complete response (CR) versus partial response (PR).

RESULTS

The study sample included 19 biopsy specimens from tumors (mean size, 3.0 cm; grade 1, n = 6; grade 2, n = 9; grade 3, n = 4) in patients treated with a mean of two conventional chemoembolization sessions. Thirteen and six tumors exhibited CR and PR, respectively, at a mean of 116 days after treatment. Tumors with CR showed a significant increase in (P < .05) or trend toward (P < .1) greater (range, 1.49-3.50 fold) pretreatment chemotherapy-sensitivity and mitosis (ATF4, BAX, CCNE1, KIF11, NFX1, PPP3CA, SNX1, TOP2A, and TOP2B) gene mRNA expression compared with tumors with PR, in addition to lower CXCL10 levels (0.48-fold), and had significantly (P < .05) higher (1.65-fold) baseline VEGFA levels.

CONCLUSIONS

Genetic signatures may allow prechemoembolization stratification of tumor response probability, and gene analysis may therefore offer an opportunity to personalize locoregional therapy by enhancing treatment modality allocation. Further corroboration of identified markers and exploration of their respective predictive capacity thresholds is necessary.