Dissecting the expression landscape of cytochromes P450 in hepatocellular carcinoma: towards novel molecular biomarkers

Spatial transcriptomics unravels palmitoylation and zonation-dependent gene regulation by AEG-1 in mouse liver

Obesity-induced metabolic dysfunction-associated steatohepatitis (MASH) leads to hepatocellular carcinoma (HCC). Astrocyte-elevated gene-1/Metadherin (AEG-1/MTDH) plays a key role in promoting MASH and HCC. AEG-1 is palmitoylated at residue cysteine 75 (Cys75) and a knock-in mouse representing mutated Cys75 to serine (AEG-1-C75S) showed activation of MASH- and HCC-promoting gene signature when compared to wild-type littermates (AEG-1-WT). The liver consists of three zones, periportal, mid-lobular, and pericentral, and zone-specific dysregulated gene expression impairs metabolic homeostasis in the liver, contributing to MASH and HCC. Here, to elucidate how palmitoylation influences AEG-1-mediated gene regulation in regard to hepatic zonation, we performed spatial transcriptomics (ST) in the livers of AEG-1-WT and AEG-1-C75S littermates. ST identified six different clusters in livers and using zone- and cell-type-specific markers we attributed specific zones and cell types to specific clusters. Ingenuity Pathway Analysis (IPA) of differentially expressed genes in each cluster unraveled activation of pro-inflammatory and MASH- and HCC-promoting pathways, mainly in periportal and pericentral hepatocytes, in AEG-1-C75S liver compared to AEG-1-WT. Interestingly, in AEG-1-C75S liver, the mid-lobular zone exhibited widespread inhibition of xenobiotic metabolism pathways and inhibition of PXR/RXR and LXR/RXR activation, versus AEG-1-WT. In conclusion, AEG-1-C75S mutant exhibited zone-specific differential gene expression, which might contribute to metabolic dysfunction and dysregulated drug metabolism leading to MASH and HCC.

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.

Distinctive Metabolism-Associated Gene Clusters That Are Also Prognostic in Intrahepatic Cholangiocarcinoma and Hepatocellular Carcinoma

Objective

To offer new prognostic evaluations by exploring potentially distinctive genetic features of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (ICC).

Methods

There were 12 samples for gene expression profiling processes in this study. These included three HCC lesion samples and their matched adjacent nontumor liver tissues obtained from patients with HCC, as well as three ICC samples and their controls collected similarly. In addition to the expression matrix generated on our own, profiles of other cohorts from The Cancer Genome Atlas (TCGA) program and the Gene Expression Omnibus (GEO) were also employed in later bioinformatical analyses. Differential analyses, functional analyses, protein interaction network analyses, and gene set variation analyses were used to identify key genes. To establish the prognostic models, univariate/multivariate Cox analyses and subsequent stepwise regression were applied, with the Akaike information criterion evaluating the goodness of fitness.

Results

The top three pathways enriched in HCC were all metabolism-related; they were fatty acid degradation, retinol metabolism, and arachidonic acid metabolism. In ICC, on the other hand, additional pathways related to fat digestion and absorption and cholesterol metabolism were identified. Consistent characteristics of such a metabolic landscape were observed across different cohorts. A prognostic risk score model for calculating HCC risk was constructed, consisting of ADH4, ADH6, CYP2C9, CYP4F2, and RDH16. This signature predicts the 3-year survival with an AUC area of 0.708 (95%CI = 0.644 to 0.772). For calculating the risk of ICC, a prognostic risk score model was built upon the expression levels of CYP26A1, NAT2, and UGT2B10. This signature predicts the 3-year survival with an AUC area of 0.806 (95% CI = 0.664 to 0.947).

Conclusion

HCC and ICC share commonly abrupted pathways associated with the metabolism of fatty acids, retinol, arachidonic acids, and drugs, indicating similarities in their pathogenesis as primary liver cancers. On the flip side, these two types of cancer possess distinctive promising biomarkers for predicting overall survival or potential targeted therapies.

Role of miR-653 and miR-29c in downregulation of CYP1A2 expression in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is a major contributor to the worldwide cancer burden. Recent studies on HCC have demonstrated dramatic alterations in expression of several cytochrome P450 (CYP) family members that play a crucial role in biotransformation of many drugs and other xenobiotics; however, the mechanisms responsible for their deregulation remain unclear.

METHODS

We investigated a potential involvement of miRNAs in downregulation of expression of CYPs observed in HCC tumors. We compared miRNA expression profiles (TaqMan Array Human MicroRNA v3.0 TLDA qPCR) between HCC human patient tumors with strong (CYP-) and weak/no (CYP+) downregulation of drug-metabolizing CYPs. The role of significantly deregulated miRNAs in modulation of expression of the CYPs and associated xenobiotic receptors was then investigated in human liver HepaRG cells transfected with relevant miRNA mimics or inhibitors.

RESULTS

We identified five differentially expressed miRNAs in CYP- versus CYP+ tumors, namely miR-29c, miR-125b1, miR-505, miR-653 and miR-675. The two most-upregulated miRNAs found in CYP- tumor samples, miR-29c and miR-653, were found to act as efficient suppressors of CYP1A2 or AHR expression.

CONCLUSIONS

Our results revealed a novel role of miR-653 and miR-29c in regulation of expresion of CYPs involved in crucial biotransformation processes in liver, which are often deregulated during liver cancer progression.

Deregulation of signaling pathways controlling cell survival and proliferation in cancer cells alters induction of cytochrome P450 family 1 enzymes

Cytochrome P450 family 1 (CYP1) enzymes contribute both to metabolism of xenobiotics and to the control of endogenous levels of ligands of the aryl hydrocarbon receptor (AhR). Their activities, similar to other CYPs, can be altered in tumor tissues. Here, we examined a possible role of proliferative/survival pathways signaling, which is often deregulated in tumor cells, and possible links with p300 histone acetyltransferase (a transcriptional co-activator) in the control of CYP1 expression, focusing particularly on CYP1A1. Using cell models derived from human liver, we observed that the induction of CYP1A1 expression, as well as other CYP1 enzymes, was reduced in exponentially growing cells, as compared with their non-dividing counterparts. The siRNA-mediated inhibition of proliferation/pro-survival signaling pathway effectors (such as β-catenin and/or Hippo pathway effectors YAP/TAZ) increased the AhR ligand-induced CYP1A1 mRNA levels in liver HepaRG cells, and/or in colon carcinoma HCT-116 cells. The activation of proliferative Wnt/β-catenin signaling in HCT-116 cells reduced both the induction of CYP1 enzymes and the binding of p300 to the promoter of CYP1A1 or CYP1B1 genes. These results seem to indicate that aberrant proliferative signaling in tumor cells could suppress induction of CYP1A1 (or other CYP1 enzymes) via competition for p300 binding. This mechanism could be involved in modulation of the metabolism of both endogenous and exogenous substrates of CYP1A1 (and other CYP1 enzymes), with possible further consequences for alterations of the AhR signaling in tumor cells, or additional functional roles of CYP1 enzymes.

Prediction of hepatocellular carcinoma risk in patients with chronic liver disease from dynamic modular networks

Background

Discovering potential predictive risks in the super precarcinomatous phase of hepatocellular carcinoma (HCC) without any clinical manifestations is impossible under normal paradigm but critical to control this complex disease.

Methods

In this study, we utilized a proposed sequential allosteric modules (AMs)-based approach and quantitatively calculated the topological structural variations of these AMs.

Results

We found the total of 13 oncogenic allosteric modules (OAMs) among chronic hepatitis B (CHB), cirrhosis and HCC network used SimiNEF. We obtained the 11 highly correlated gene pairs involving 15 genes (r > 0.8, P < 0.001) from the 12 OAMs (the out-of-bag (OOB) classification error rate < 0.5) partial consistent with those in independent clinical microarray data, then a three-gene set (cyp1a2-cyp2c19-il6) was optimized to distinguish HCC from non-tumor liver tissues using random forests with an average area under the curve (AUC) of 0.973. Furthermore, we found significant inhibitory effect on the tumor growth of Bel-7402, Hep 3B and Huh7 cell lines in zebrafish treated with the compounds affected those three genes.

Conclusions

These findings indicated that the sequential AMs-based approach could detect HCC risk in the patients with chronic liver disease and might be applied to any time-dependent risk of cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-021-02791-9.

Hepatic transcriptome and DNA methylation patterns following perinatal and chronic BPS exposure in male mice

Background

Bisphenol S (BPS) is a common bisphenol A (BPA) substitute, since BPA is virtually banned worldwide. However, BPS and BPA have both endocrine disrupting properties. Their effects appear mostly in adulthood following perinatal exposures. The objective of the present study was to investigate the impact of perinatal and chronic exposure to BPS at the low dose of 1.5 μg/kg body weight/day on the transcriptome and methylome of the liver in 23 weeks-old C57BL6/J male mice.

Results

This multi-omic study highlights a major impact of BPS on gene expression (374 significant deregulated genes) and Gene Set Enrichment Analysis show an enrichment focused on several biological pathways related to metabolic liver regulation. BPS exposure also induces a hypomethylation in 58.5% of the differentially methylated regions (DMR). Systematic connections were not found between gene expression and methylation profile excepted for 18 genes, including 4 genes involved in lipid metabolism pathways (Fasn, Hmgcr, Elovl6, Lpin1), which were downregulated and featured differentially methylated CpGs in their exons or introns.

Conclusions

This descriptive study shows an impact of BPS on biological pathways mainly related to an integrative disruption of metabolism (energy metabolism, detoxification, protein and steroid metabolism) and, like most high-throughput studies, contributes to the identification of potential exposure biomarkers.

Graphical abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07294-3.

Epigenetic Regulation of Differentially Expressed Drug-Metabolizing Enzymes in Cancer

Drug metabolism is a biotransformation process of drugs, catalyzed by drug-metabolizing enzymes (DMEs), including phase I DMEs and phase II DMEs. The aberrant expression of DMEs occurs in the different stages of cancer. It can contribute to the development of cancer and lead to individual variations in drug response by affecting the metabolic process of carcinogen and anticancer drugs. Apart from genetic polymorphisms, which we know the most about, current evidence indicates that epigenetic regulation is also central to the expression of DMEs. This review summarizes differentially expressed DMEs in cancer and related epigenetic changes, including DNA methylation, histone modification, and noncoding RNAs. Exploring the epigenetic regulation of differentially expressed DMEs can provide a basis for implementing individualized and rationalized medication. Meanwhile, it can promote the development of new biomarkers and targets for the diagnosis, treatment, and prognosis of cancer. SIGNIFICANCE STATEMENT: This review summarizes the aberrant expression of DMEs in cancer and the related epigenetic regulation of differentially expressed DMEs. Exploring the epigenetic regulatory mechanism of DMEs in cancer can help us to understand the role of DMEs in cancer progression and chemoresistance. Also, it provides a basis for developing new biomarkers and targets for the diagnosis, treatment, and prognosis of cancer.

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).