The down-regulation of the CYP2C19 gene is associated with aggressive tumor potential and the poorer recurrence-free survival of hepatocellular carcinoma

Arachidonic acid metabolism CYP450 pathway is deregulated in hepatocellular carcinoma and associated with microvascular invasion

Arachidonic acid metabolism plays a crucial role in the development and progression of inflammatory and metabolic liver diseases. However, its role in hepatocellular carcinoma (HCC) remains unclear. In this study, we investigated the expression of key genes involved in the arachidonic acid metabolism pathway in HCC using a combination of bioinformatics, proteomics and immunohistochemistry analyses. Through a comprehensive analysis of publicly available datasets, clinical HCC tissues, and tissue microarrays, we compared the expression of hepatic arachidonic acid metabolic genes. We observed significant downregulation of cytochrome P450 (CYP450) pathway genes at both the messenger RNA and protein levels in HCC tissues compared to normal liver tissues. Furthermore, we observed a strong correlation between the deregulation of the arachidonic acid metabolism CYP450 pathway and the pathological features and prognosis of HCC. Specifically, the expression of CYP2C8/9/18/19 was significantly correlated with pathological grade (r = -.484, p < .0001), vascular invasion (r = -.402, p < .0001), aspartate transaminase (r = -.246, p = .025), gamma-glutamyl transpeptidase (r = -.252, p = .022), alkaline phosphatase (r = -.342, p = .002), alpha-fetoprotein (r = -.311, p = .004) and carbohydrate antigen 19-9 (r = -.227, p = .047). Moreover, we discovered a significant association between CYP450 pathway activity and vascular invasion in HCC. Collectively, these data indicate that arachidonic acid CYP450 metabolic pathway deregulation is implicated in HCC progression and may be a potential predictive factor for early recurrence in patients with HCC.

Construction of HBV-HCC prognostic model and immune characteristics based on potential genes mining through protein interaction networks

OBJECTIVE

To search for human protein-coding genes related to hepatocellular carcinoma (HCC) in the context of hepatitis B virus (HBV) infection, and perform prognosis risk assessment.

METHODS

Genes related to HBV-HCC were selected through literature screening and protein-protein interaction (PPI) network database analysis. Prognosis potential genes (PPGs) were identified using Cox regression analysis. Patients were divided into high-risk and low-risk groups based on PPGs, and risk scores were calculated. Kaplan-Meier plots were used to analyze overall survival rates, and the results were predicted based on clinicopathological variables. Association analysis was also conducted with immune infiltration, immune therapy, and drug sensitivity. Experimental verification of the expression of PPGs was done in patient liver cancer tissue and normal liver tissue adjacent to tumors.

RESULTS

The use of a prognosis potential genes risk assessment model can reliably predict the prognosis risk of patients, demonstrating strong predictive ability. Kaplan-Meier analysis showed that the overall survival rate of the low-risk group was significantly higher than that of the high-risk group. There were significant differences between the two subgroups in terms of immune infiltration and IC50 association analysis. Experimental verification revealed that CYP2C19, FLNC, and HNRNPC were highly expressed in liver cancer tissue, while UBE3A was expressed at a lower level.

CONCLUSION

PPGs can be used to predict the prognosis risk of HBV-HCC patients and play an important role in the diagnosis and treatment of liver cancer. They also reveal their potential role in the tumor immune microenvironment, clinical-pathological characteristics, and prognosis.

Remnant liver function is associated with long-term survival in patients with hepatocellular carcinoma undergoing hepatectomy

It is important to assess the prognosis and intervene before and after surgery in patients with hepatocellular carcinoma. This study aims to elucidate the association of outcomes and residual liver function after hepatectomy. A total of 176 patients who underwent the initial resection for hepatocellular carcinoma between January 2011 and March 2021 at Jichi Medical University were included. Hepatic clearance of the remnant liver was measured using 99mTc-galactosyl serum albumin scintigraphy. The log-rank test was used to analyze survival using the Kaplan-Meier method. Hazard ratios (HR) and 95% confidence intervals (CI) for overall survival were calculated using Cox's proportional hazard model. In multivariate analysis, microvascular invasion, intraoperative blood loss, and hepatic clearance of the remnant liver were independently associated with overall survival. Hepatic clearance of the remnant liver was independently associated with recurrence free survival. This is the first report to show that lower residual liver function is associated with shorter survival in patients with hepatocellular carcinoma undergoing hepatectomy. Preoperative determination of remnant liver function may allow assessment of prognosis in patients planned to undergo resection of hepatocellular carcinoma. Preservation of liver functional reserve may be crucial for improved long-term outcomes after hepatectomy.

Exploring the potential of ToxCast™ data for mechanism-based prioritization of chemicals in regulatory context: Case study with priority existing chemicals (PECs) under K-REACH

Recent studies have highlighted the potential of ToxCast™ database to mechanism-based prioritization of chemicals. To explore the applicability of ToxCast data in the context of regulatory inventory chemicals, we screened 510 priority existing chemicals (PECs) regulated under the Act on the Registration and Evaluation of Chemical Substances (K-REACH) using ToxCast bioassays. In our analysis, a hit-call data matrix containing 298984 chemical-gene interactions was computed for 949 bioassays with the intended target genes, which enabled the identification of the putative toxicity mechanisms. Based on the reactivity to the chemicals, we analyzed 412 bioassays whose intended target gene families were cytochrome P450, oxidoreductase, transporter, nuclear receptor, steroid hormone, and DNA-binding. We also identified 141 chemicals based on their reactivity in the bioassays. These chemicals are mainly in consumer products including colorants, preservatives, air fresheners, and detergents. Our analysis revealed that in vitro bioactivities were involved in the relevant mechanisms inducing in vivo toxicity; however, this was not sufficient to predict more hazardous chemicals. Overall, the current results point to a potential and limitation in using ToxCast data for chemical prioritization in regulatory context in the absence of suitable in vivo data.

Comparative RNA-Sequencing Analysis Reveals High Complexity and Heterogeneity of Transcriptomic and Immune Profiles in Hepatocellular Carcinoma Tumors of Viral (HBV, HCV) and Non-Viral Etiology

Background and Objectives: Hepatocellular carcinoma (HCC), the most common type of primary liver cancer, is the leading cause of cancer-related mortality. It arises and progresses against fibrotic or cirrhotic backgrounds mainly due to infection with hepatitis viruses B (HBV) or C (HCV) or non-viral causes that lead to chronic inflammation and genomic changes. A better understanding of molecular and immune mechanisms in HCC subtypes is needed. Materials and Methods: To identify transcriptional changes in primary HCC tumors with or without hepatitis viral etiology, we analyzed the transcriptomes of 24 patients by next-generation sequencing. Results: We identified common and unique differentially expressed genes for each etiological tumor group and analyzed the expression of SLC, ATP binding cassette, cytochrome 450, cancer testis, and heat shock protein genes. Metascape functional enrichment analysis showed mainly upregulated cell-cycle pathways in HBV and HCV and upregulated cell response to stress in non-viral infection. GeneWalk analysis identified regulator, hub, and moonlighting genes and highlighted CCNB1, ACTN2, BRCA1, IGF1, CDK1, AURKA, AURKB, and TOP2A in the HCV group and HSF1, HSPA1A, HSP90AA1, HSPB1, HSPA5, PTK2, and AURKB in the group without viral infection as hub genes. Immune infiltrate analysis showed that T cell, cytotoxic, and natural killer cell markers were significantly more highly expressed in HCV than in non-viral tumors. Genes associated with monocyte activation had the highest expression levels in HBV, while high expression of genes involved in primary adaptive immune response and complement receptor activity characterized tumors without viral infection. Conclusions: Our comprehensive study underlines the high degree of complexity of immune profiles in the analyzed groups, which adds to the heterogeneous HCC genomic landscape. The biomarkers identified in each HCC group might serve as therapeutic targets.

Blocking hepatocarcinogenesis by a cytochrome P450 family member with female-preferential expression

OBJECTS

The incidence of hepatocellular carcinoma (HCC) shows an obvious male dominance in rodents and humans. We aimed to identify the key autosomal liver-specific sex-related genes and investigate their roles in hepatocarcinogenesis.

DESIGN

Two HCC cohorts (n=551) with available transcriptome and metabolome data were used. Class comparisons of omics data and ingenuity pathway analysis were performed to explore sex-related molecules and their associated functions. Functional assays were employed to investigate roles of the key candidates, including cellular assays, molecular assays and multiple orthotopic HCC mouse models.

RESULTS

A global comparison of multiple omics data revealed 861 sex-related molecules in non-tumour liver tissues between female and male HCC patients, which denoted a significant suppression of cancer-related diseases and functions in female liver than male. A member of cytochrome P450 family, CYP39A1, was one of the top liver-specific candidates with significantly higher levels in female vs male liver. In HCC tumours, CYP39A1 expression was dramatically reduced in over 90% HCC patients. Exogenous CYP39A1 significantly blocked tumour formation in both female and male mice and partially reduced the sex disparity of hepatocarcinogenesis. The HCC suppressor role of CYP39A1 did not rely on its known P450 enzyme activity but its C-terminal region, by which CYP39A1 impeded the transcriptional activation activity of c-Myc, leading to a significant inhibition of hepatocarcinogenesis.

CONCLUSIONS

The liver-specific CYP39A1 with female-preferential expression was a strong suppressor of HCC development. Strategies to up-regulate CYP39A1 might be promising methods for HCC treatment in both women and men in future.

Epigenetic priming in chronic liver disease impacts the transcriptional and genetic landscapes of hepatocellular carcinoma

Hepatocellular carcinomas (HCCs) usually arise from chronic liver disease (CLD). Precancerous cells in chronically inflamed environments may be 'epigenetically primed', sensitising them to oncogenic transformation. We investigated whether epigenetic priming in CLD may affect HCC outcomes by influencing the genomic and transcriptomic landscapes of HCC. Analysis of DNA methylation arrays from 10 paired CLD-HCC identified 339 shared dysregulated CpG sites and 18 shared differentially methylated regions compared with healthy livers. These regions were associated with dysregulated expression of genes with relevance in HCC, including ubiquitin D (UBD), cytochrome P450 family 2 subfamily C member 19 (CYP2C19) and O-6-methylguanine-DNA methyltransferase (MGMT). Methylation changes were recapitulated in an independent cohort of nine paired CLD-HCC. High CLD methylation score, defined using the 124 dysregulated CpGs in CLD and HCC in both cohorts, was associated with poor survival, increased somatic genetic alterations and TP53 mutations in two independent HCC cohorts. Oncogenic transcriptional and methylation dysregulation is evident in CLD and compounded in HCC. Epigenetic priming in CLD sculpts the transcriptional landscape of HCC and creates an environment favouring the acquisition of genetic alterations, suggesting that the extent of epigenetic priming in CLD could influence disease outcome.

Identification of Key Functional Modules and Immunomodulatory Regulators of Hepatocellular Carcinoma

Despite the advances in the treatment of hepatocellular carcinoma (HCC), the prognosis of HCC patients remains unsatisfactory due to postsurgical recurrence and treatment resistance. Therefore, it is important to reveal the mechanisms underlying HCC and identify potential therapeutic targets against HCC, which could facilitate the development of novel therapies. Based on 12 HCC samples and 12 paired paracancerous normal tissues, we identified differentially expressed mRNAs and lncRNAs using the "limma" package in R software. Moreover, we used the weighted gene coexpression network analysis (WGCNA) to analyze the expression data and screened hub genes. Furthermore, we performed pathway enrichment analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. In addition, the relative abundance of a given gene set was estimated by single-sample Gene Set Enrichment Analysis. We identified 687 differentially expressed mRNAs and 260 differentially expressed lncRNAs. A total of 6 modules were revealed by WGCNA, and MT1M and MT1E genes from the red module were identified as hub genes. Moreover, pathway analysis revealed the top 10 enriched KEGG pathways of upregulated or downregulated genes. Additionally, we also found that CD58 might act as an immune checkpoint gene in HCC via PD1/CTLA4 pathways and regulate the levels of tumor-infiltrating immune cells in HCC tissues, which might be an immunotherapeutic target in HCC. Our research identified key functional modules and immunomodulatory regulators for HCC, which might offer novel diagnostic biomarkers and/or therapeutic targets for cancer immunotherapy.

Global analysis of gene expression signature and diagnostic/prognostic biomarker identification of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common cancers in the world. The landscape of HCC's molecular alteration signature has been explored over the last few decades. Even so, more comprehensive research is still needed to improve understanding of tumorigenesis and progression of HCC, as well as to identify potential biomarkers for the malignancy. In this research, a comprehensive bioinformatics analysis was conducted based on the publicly available databases from both the Cancer Genome Atlas (TCGA) program and the gene expression omnibus (GEO) database. R/Bioconductor was used to analyze differentially expressed genes (DEGs) between HCC tumor and normal control (NC) samples, and then a protein-protein interaction (PPI) network of DEGs was established through the STRING platform. Finally, the application of specific candidate genes as diagnostic or prognostic biomarkers of HCC was explored and evaluated by ROC and survival analysis. A total of 310 DEGs were detected in the HCC tumor samples. Thirty-six hub DEGs in the PPI network and 10 candidates of the 36 genes showed significant alterations in tumor expression, including CDKN3, TOP2A, UBE2C, CDC20, PBK, ASPM, KIF20A, NCAPG, CCNB2, CYP3A4. The 10-gene signature had relatively significant effects when distinguishing tumors from normal samples (sensitivity >70%, specificity >70%, AUC >0.8, p < 0.001). Eight candidate genes were negatively correlated with the overall survival rate of the patients (p < 0.05) and were all up-regulated in HCC tumor samples. The age and gender factors had no significant impact on the overall survival rate of HCC patients (p > 0.05), and the TNM stage status factor had a significant negative prognosis correlation (p < 0.05). This research provides evidence for a better understanding of tumorigenesis and progression of HCC and helps to explore candidate targets for disease diagnosis and treatment.

Hepatocellular carcinoma: Gene expression profiling and regulation of xenobiotic-metabolizing cytochromes P450

Hepatocellular carcinoma (HCC) remains a highly prevalent and deadly disease, being among the top causes of cancer-related deaths worldwide. Despite the fact that the liver is the major site of biotransformation, studies on drug metabolizing enzymes in HCC are scarce. It is known that malignant transformation of hepatocytes leads to a significant alteration of their metabolic functions and overall deregulation of gene expression. Advanced stages of the disease are thus frequently associated with liver failure, and severe alteration of drug metabolism. However, the impact of dysregulation of metabolic enzymes on therapeutic efficacy and toxicity in HCC patients is largely unknown. Here we demonstrate a significant down-regulation in European Caucasian patients of cytochromes P450 (CYPs), the major xenobiotic-metabolizing enzymes, in HCC tumour samples as compared to their surrounding non-cancerous (reference) tissue. Moreover, we report for the first time the association of the unique CYP profiles with specific transcriptome changes, and interesting correlations with expression levels of nuclear receptors and with the histological grade of the tumours. Integrated analysis has suggested certain co-expression profiles of CYPs with lncRNAs that need to be further characterized. Patients with large tumours with down-regulated CYPs could be more vulnerable to drug toxicity; on the other hand, such tumours would eliminate drugs more slowly and should be more sensitive to pharmacotherapy (except in the case of pro-drugs where activation is necessary).

Identification of a robust functional subpathway signature for pancreatic ductal adenocarcinoma by comprehensive and integrated analyses

BACKGROUND

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy and its mortality continues to rise globally. Because of its high heterogeneity and complex molecular landscapes, published gene signatures have demonstrated low specificity and robustness. Functional signatures containing a group of genes involved in similar biological functions may display a more robust performance.

METHODS

The present study was designed to excavate potential functional signatures for PDAC by analyzing maximal number of datasets extracted from available databases with a recently developed method of FAIME (Functional Analysis of Individual Microarray Expression) in a comprehensive and integrated way.

RESULTS

Eleven PDAC datasets were extracted from GEO, ICGC and TCGA databases. By systemically analyzing these datasets, we identified a robust functional signature of subpathway (path:00982_1), which belongs to the drug metabolism-cytochrome P450 pathway. The signature has displayed a more powerful and robust capacity in predicting prognosis, drug response and chemotherapeutic efficacy for PDAC, particularly for the classical subtype, in comparison with published gene signatures and clinically used TNM staging system. This signature was verified by meta-analyses and validated in available cell line and clinical datasets with chemotherapeutic efficacy.

CONCLUSION

The present study has identified a novel functional PDAC signature, which has the potential to improve the current systems for predicting the prognosis and monitoring drug response, and to serve a linkage to therapeutic options for combating PDAC. However, the involvement of path:00982_1 subpathway in the metabolism of anti-PDAC chemotherapeutic drugs, particularly its biological interpretation, requires a further investigation. Video Abstract.

Molecular Pathological Differences in Global Gene Expression between Two Sustained Proliferative Lesions, Nodular Regenerative Hepatocellular Hyperplasia and Hepatocellular Adenoma, in Mice

Long-term exposure to piperonyl butoxide (PBO) induces multiple nodular masses along with hepatocellular tumors in the liver of mice. The histopathological features of the nodules led to our diagnosis of nodular regenerative hepatocellular hyperplasia (NRH). However, because of the lack of data on the biological characteristics of NRH, whether this lesion is truly nonneoplastic remains unknown. In this study, the molecular characteristics of NRH were compared with those of hepatocellular adenoma (HCA) by global gene expression analysis. Six-week-old male ICR mice were fed a diet containing 6,000 ppm PBO for 43 weeks to induce NRH and HCA development. Complementary DNA microarray analysis was performed using messenger RNA extracted from NRH and HCA frozen sections collected by laser microdissection. Hierarchical cluster analysis showed that all NRH samples clustered together but were separate from the HCA cluster. Pathway analysis revealed activation of the cell cycle and Delta-Notch signaling in both lesions, but the latter was more upregulated in HCA. Downregulation of cytochrome p450 enzymes was observed in NRH, but not in HCA. These results imply that NRH differs from HCA in terms of not only morphological but also molecular characteristics.