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

Genetic association of tertiary lymphoid structure-related gene signatures with HCC based on Mendelian randomization and machine learning and construction of prognosis model

BACKGROUND

Tertiary lymphoid structures (TLS) are formed in numerous cancer types. However, their value and significance in hepatocellular carcinoma (HCC) is unclear.

METHODS

We performed differential genes expression analysis of TLS-related Genes (TLSG) based on The Cancer Genome Atlas (TCGA) database, and performed Mendelian randomization (MR) analysis using expression quantitative trait loci, and then took their intersecting genes. A TLSG prognostic signature (TLSGPS)-based risk score was constructed using Least Absolute Shrinkage and Selection Operator (LASSO), univariate and multivariate COX regression analysis, and survival analysis was then performed. We used the International Cancer Genome Consortium for outside validation. We also performed biological function, tumor mutational burden, immune infiltration, single-cell analysis, CeRNA and drug sensitivity analysis based on TLSGPS.

RESULTS

Three TLSGs (HM13, CSTB, CDCA7L) were identified to construct the TLSGPS, which showed good predictive ability and outperformed most prognostic signatures. MR suggested that HM13 (OR = 0.9997, 95 %CI: 0.9994-0.9999, P = 0.014) and CSTB (OR = 0.9997, 95 %CI: 0.9995-0.9999, P = 0.048) were negatively correlated with the risk of HCC onset, while CDCA7L (OR = 1.0004, 1.0001-1.0007, P = 0.0161) was the opposite. The differences in biological functions between the TLSGPS-based high-risk group (HRG) and low-risk group (LRG) involved cell proliferation, differentiation, and drug metabolism. HRG plus high mutations exhibited extremely poor survival. HRG had higher abundance of immune cell-oncogenic phenotypes, higher immune escape ability, and greater sensitivity to Afatinib, Dasatinib, and Gefitinib.

CONCLUSION

3 TLSGs identified by machine learning and MR can predict the onset, prognosis and clinical treatment of HCC patients, and had significant genetic association with HCC.

Sohlh2 Promotes the Progression of Hepatocellular Carcinoma via TGM2-Mediated Autophagy

Liver cancer is the third leading cause of cancer-related deaths worldwide, with hepatocellular carcinoma (HCC) accounting for 85% of liver cancer-related deaths. Autophagy controls HCC cell growth, invasion, metastasis, drug resistance, and stemness. Spermatogenesis and oogenesis basic helix-loop-helix transcription factor 2 (Sohlh2) can bind to the E-boxes in the promoter regions of target genes, which are involved in multiple neoplasms. In this study, Sohlh2 was highly expressed in HCC tissues and was related to poor prognosis. Moreover, Sohlh2 overexpression promoted the proliferation, migration, invasion, and metastasis of HCC cells in vivo and in vitro. However, Sohlh2 silencing inhibited proliferation, migration, invasion, and metastasis of HCC cells in vivo and in vitro. Mechanistically, Sohlh2 could bind to the promoter of TGM2 and enhance its transcriptional activity, thereby enhancing the autophagy of HCC cells. Furthermore, Sohlh2 protein levels were positively associated with TGM2 expression in HCC tissues. Taken together, these results demonstrate that Sohlh2 can promote HCC progression via TGM2-mediated autophagy, implying that Sohlh2 is a promising candidate for HCC treatment.

When less is more: The association between the expression of polymorphic CYPs and AFB1-induced HCC

BACKGROUND

An individual's genetic fingerprint is emerging as a pivotal predictor of numerous disease- and treatment-related factors. Single nucleotide polymorphisms (SNPs) in drug-metabolizing enzymes play key roles in an individual's exposure to a malignancy-associated risk, such as Aflatoxin B1 (AFB1)-induced hepatocellular carcinoma (HCC).

AIM

This study aimed at reviewing literature on the polymorphisms that exist in CYP enzymes and their possible link with susceptibility to AFB1-induced HCC.

MATERIALS & METHODS

A set of keywords associated with the study subject of interest was used to search the Google Scholar and the PubMed database. The last ten years' worth of research projects were included in the results filter. The research involved HCC patients and any connection between polymorphic forms of CYP enzymes and their susceptibility to AFB1-induced HCC, including older but significant data.

RESULTS

Variations in CYP1A2 and CYP3A4 were reported to impact the rate and magnitude of AFB1 bio-activation, thus influencing an individual's vulnerability to develop HCC. In HCC patients, the activity of CYP isoforms varies, where increased activity has been reported with CYP2C9, CYP2D6, and CYP2E1, while CYP1A2, CYP2C8, and CYP2C19 exhibit decreased activity. CYP2D6*10 frequency has been discovered to differ considerably in HCC patients. Rs2740574 (an upstream polymorphism in CYP3A4 as detected in CYP3A4*1B) and rs776746 (which affects CYP3A5 RNA splicing), both of which influence CYP3A expression, thus impacting the variability of AFB1-epoxide adducts in HCC patients.

DISCUSSION

CYP1A2 is the primary enzyme accountable for the formation of harmful AFBO globally. CYP3A4, CYP3A5, CYP3A7, CYP2B7, and CYP3A3 are also implicated in the bio-activation of AFB1 to mutagenic metabolites. It is thought that CYP3A4 is the protein that interacts with AFB1 metabolism the most.

CONCLUSION

Polymorphic variants of CYP enzymes have a functional impact on the susceptibility to AFB1-induced HCC. Outlining such variation and their implications may provide deeper insights into approaching HCC in a more personalized manner for guiding future risk-assessment, diagnosis, and treatment.

Mechanism of atorvastatin in treating hepatocellular carcinoma: a study based on network pharmacology, molecular docking, and bioinformatics analysis

Hepatocellular carcinoma (HCC) is a tumor with high morbidity and mortality. Current research suggests that statins may aid in its prevention and treatment, while studies on the associated mechanisms remain limited. Therefore, we aim to reveal the mechanism of atorvastatin treatment for HCC by using network pharmacology and bioinformatics methods. The databases SwissTargetPrediction, PharmMapper, and DrugBank were utilized to obtain targets of atorvastatin, while GSE169289, GSE135631, and GSE207435 were used to identify differentially expressed genes (DEGs) for HCC. The overlap between the two groups was used to identify atorvastatin's target for treating HCC. Following protein-protein interaction (PPI) analysis, hub genes were identified using Cytoscape software and LASSO analysis. The hub genes were further validated using data from The Cancer Genome Atlas (TCGA) and The Human Protein Atlas (HPA) databases. To evaluate the clinical significance of the hub genes, Kaplan-Meier (KM) survival analysis and Cox analysis were conducted. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) were performed to investigate potential mechanisms. Finally, molecular docking analysis was performed to validate the interaction between atorvastatin and the hub genes. A total of 1948 DEGs of HCC and 380 targets of atorvastatin were identified, respectively. After taking the intersection, 79 genes were identified as potential targets of atorvastatin for HCC treatment. After multiple screening methods, CYP2C9 was ultimately identified as the hub gene. Analysis of data from TCGA and HPA databases showed reduced expression of CYP2C9 in HCC tissues. KM and Cox analysis showed a favorable prognosis for HCC patients with high CYP2C9 expression. KEGG and GSEA indicated that metabolism of xenobiotics by cytochrome P450, and PPAR signaling pathway could be the potential mechanisms for atorvastatin in treating HCC. Molecular docking analysis revealed that atorvastatin binds to CYP2C9 with a binding energy of - 8.837, indicating highly stable binding. CYP2C9 is associated with the prognosis of HCC patients and could serve as a potential target for atorvastatin treatment in HCC.

Identification of ribosomal protein S21 as a potential prognostic and immunotherapy biomarker for hepatocellular carcinoma

BACKGROUND

Recent studies show that ribosomal protein S21 (RPS21) plays a role in the development and progression of various malignancies. However, the biological value of RPS21 in hepatocellular carcinoma (HCC) and its association with immunotherapy remain unknown.

METHODS

Here, we examined the differential expression of RPS21 between HCC and normal liver tissues, using the TCGA, ICGC and GEO databases, followed by verification by reverse-transcription quantitative polymerase chain reaction (RT-qPCR) in LO2, SMMC7721, HepG2, and MHCC-97H cell lines. Kaplan-Meier and Cox regression analyses were applied to investigate how RPS21 expression influenced overall survival, and a nomogram was established to predict prognosis among HCC patients. We further analyzed how RPS21 expression was related to tumor immune microenvironment, immunotherapy efficiency, and genomic alterations, and investigated potential underlying mechanisms.

RESULTS

RPS21 upregulation was observed in HCC tissues and cell lines, compared to normal controls. Survival analysis revealed that RPS21 overexpression was significantly associated with poor clinical outcomes (all p < 0.05). Functional enrichment analyses indicated that differentially expressed genes relative to RPS21 expression were mainly involved in tumor response, proliferation, and metabolism. Additionally, RPS21 expression was positively correlated with the infiltration of activated CD4+ T cells and tumor mutational burden (all p < 0.05). Moreover, RPS21 was co-expressed with immune-related genes and immune checkpoint genes. Analyses of drug sensitivity predict that HCC patients with low RPS21 expression were more sensitive to targeted immunotherapy.

CONCLUSIONS

The present results suggested that RPS21 might be a promising prognostic marker and a potential immunotherapy target for patients with HCC.

Gene Expression Reprogramming by Citrate Supplementation Reduces HepG2 Cell Migration and Invasion

Citrate, which is obtained from oxaloacetate and acetyl-CoA by citrate synthase in mitochondria, plays a key role in both normal and cancer cell metabolism. In this work, we investigated the effect of 10 mM extracellular citrate supplementation on HepG2 cells. Gene expression reprogramming was evaluated by whole transcriptome analysis using gene set enrichment analysis (GSEA). The transcriptomic data were validated through analyzing changes in the mRNA levels of selected genes by qRT-PCR. Citrate-treated cells exhibited the statistically significant dysregulation of 3551 genes; 851 genes were upregulated and 822 genes were downregulated. GSEA identified 40 pathways affected by differentially expressed mRNAs. The most affected biological processes were related to lipid and RNA metabolism. Several genes of the cytochrome P450 family were upregulated in treated cells compared to controls, including the CYP3A5 gene, a tumor suppressor in hepatocellular carcinoma (HCC) that plays an important protective role in HCC metastasis. The citrate-induced dysregulation of cytochromes could both improve the effectiveness of chemotherapeutics used in combination and reduce the aggressiveness of tumors by diminishing cell migration and invasion.

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.

The cytotoxic and antiproliferative properties of ruthenium nitrosyl complexes and their modulation effect on cytochrome P450 in the HepG2 cell line

Ruthenium nitrosyl complexes are actively investigated as antitumor agents. Evaluation of potential interactions between cytochromes P450 (CYPs) with new compounds is carried out regularly during early drug development. In this study we have investigated the cytotoxic and antiproliferative activities of ruthenium nitrosyl complexes with methyl/ethyl esters of nicotinic and isonicotinic acids and γ-picoline against 2D and 3D cultures of human hepatocellular carcinoma HepG2 and non-cancer human lung fibroblasts MRC-5, assessed their photoinduced activity at λrad = 445 nm, and also evaluated their modulating effect on CYP3A4, CYP2C9, and CYP2C19. The study of cytotoxic and antiproliferative activities against 2D and 3D cell models was performed using phenotypic-based high content screening (HCS). The expression of CYP3A4, CYP2C9, and CYP2C19 mRNAs and CYP3A4 protein was examined using target-based HCS. The results of CYP3A4 mRNA expression were confirmed by real-time reverse transcription-polymerase chain reaction (RT-PCR). The ruthenium nitrosyl complexes exhibited a dose-dependent cytotoxic effect against HepG2 and MRC-5 cells. The cytotoxic activity of complexes with ethyl isonicotinate (1) and nicotinate (3, 4) was significantly lower for MRC-5 than for HepG2, for a complex with methyl isonicotinate (2) it was higher for MRC-5 than for HepG2, for a complex with γ-picoline (5) it was comparable for both lines. The antiproliferative effect of complexes 2 and 5 was one order of magnitude higher for MRC-5; for complexes 1, 3, and 4 it was comparable for both lines. The cytotoxic activity of all compounds for 3D HepG2 was lower than for 2D HepG2, with the exception of 4. Photoactivation affected the activity of complex 1 only. Its cytotoxic activity decreased, while the antiproliferative activity increased. The ruthenium nitrosyl complexes 1-4 acted as inducers of CYP3A4 and CYP2C19, while the complex with γ-picoline (5) induced of CYP3A4. Among the studied ruthenium nitrosyl complexes, the most promising potential antitumor compound is the ruthenium compound with methyl nicotinate (4).

Synthesis, Structure, and Magnetic and Biological Properties of Copper(II) Complexes with 1,3,4-Thiadiazole Derivatives

New coordination compounds of copper(II) with 2,5-bis(ethylthio)-1,3,4-thiadiazole (L1) and 2,5-bis(pyridylmethylthio)-1,3,4-thiadiazole (L2) with compositions Cu(L1)2Br2, Cu(L1)(C2N3)2, Cu(L2)Cl2, and Cu(L2)Br2 were prepared. The complexes were identified and studied by CHN analysis, infrared (IR) spectroscopy, powder X-Ray diffraction (XRD), and static magnetic susceptibility. The crystal structures of Cu(II) complexes with L1 were determined. The structures of the coordination core of complexes Cu(L2)Cl2 and Cu(L2)Br2 were determined by Extended X-ray absorption fine structure (EXAFS) spectroscopy. Magnetization measurements have revealed various magnetic states in the studied complexes, ranging from an almost ideal paramagnet in Cu(L1)2Br2 to alternating-exchange antiferromagnetic chains in Cu(L1)(C2N3)2, where double dicyanamide bridges provide an unusually strong exchange interaction (J1/kB ≈ -23.5 K; J2/kB ≈ -20.2 K) between Cu(II) ions. The cytotoxic activity of copper(II) complexes with L2 was estimated on the human cell lines of breast adenocarcinoma (MCF-7) and hepatocellular carcinoma (HepG2).

Enhanced Cytotoxic Activity of 6-Mercaptopurine-Loaded Solid Lipid Nanoparticles in Hepatic Cancer Treatment

6-Mercaptopurine (6-MCP) is an antiproliferative purine analog used in acute lymphoblastic leukemia, non-Hodgkin lymphoma, and inflammatory bowel disease (Crohn's disease, ulcerative colitis). Although 6-MCP has the great therapeutic potential for cancer and immunosuppressant-related diseases, 6-MCP is not readily soluble in water, presents a high first-pass effect, short half-life (0.5-1.5 h), and implies a low bioavailability (16%). On the contrary, solid lipid nanoparticles (SLNs) are prepared from solid lipids at room temperature and body temperature. In this study, SLNs were prepared w/o/w double emulsion-solvent evaporation method using Precirol ATO5 as matrix lipid. In the emulsion stabilization, surfactant (Tween 80) and polymeric stabilizer (polyvinyl alcohol [PVA]) were used. Two group formulations using Tween 80 and PVA were compared in terms of particle size, polydispersity index, zeta potential encapsulation efficiency%, and process yield%. Differential calorimetric analysis and release properties were examined for optimum formulation, and release kinetics were calculated. According to studies, sustained release was obtained with SLNs by the Korsmayer-Peppas kinetic model. The in vitro cytotoxicity studies were performed on the hepatocarcinoma (HEP3G) cell line. According to the results, successful SLN formulations were produced, and PVA was found best stabilizer. Optimum formulation exhibited significantly higher cytotoxic effects on HEP3G than on pure 6-MCP. These results demonstrated that solid lipid nanodrug delivery systems have great potential for formulation of 6-MCP.

Identification of biomarkers of hepatocellular carcinoma gene prognosis based on the immune-related lncRNA signature of transcriptome data

BACKGROUND

Long non-coding RNAs (lncRNAs) are well established to have an important role in cancer. The goal of this research was to investigate the prognostic usefulness of putative immune-related lncRNAs in hepatocellular carcinoma (HCC).

METHODS

The developed lncRNA signature was validated using 343 HCC patients from The Cancer Genome Atlas (TCGA) and 81 samples from Gene Expression Omnibus (GEO). Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO) analysis were used to analyze immune-related lncRNAs for HCC prognosis. Patients in the low-risk group survived substantially longer than those in the high-risk group (P < 0.05). The discovered signal might be a useful prognostic factor for predicting patient survival. Overall survival predicted some clinical net improvements, according to the nomogram. Numerous enrichment approaches (including gene set enrichment analysis) were utilized to investigate the underlying mechanisms.

RESULTS

Drug metabolism, mTOR, and p53 signaling pathways were associated with high-risk groups. When the expression of lncRNA PRRT3-AS1 was silenced in HepG2 cells, the proliferation, migration, and invasion abilities of HepG2 cells were decreased, and apoptosis was enhanced. In the supernatant from HepG2 cells with PRRT3-AS1 knockdown, the anti-inflammatory factors IL-10 and TGF-1 were induced, whereas the pro-inflammatory factors IL-1β, TNF-α, and IL-6 were reduced (P < 0.05). After PRRT3-AS1 knockdown, the protein expression of CD24, THY1, LYN, CD47, and TRAF2 in HepG2 cells was attenuated (P < 0.05).

CONCLUSION

The discovery of five immune-related lncRNA signatures has significant therapeutic significance for predicting patient prognosis and directing personalized treatment for patients with HCC, which requires additional prospective confirmation.

Exploring the mechanisms underlying the therapeutic effect of the Radix Bupleuri-Rhizoma Cyperi herb pair on hepatocellular carcinoma using multilevel data integration and molecular docking

Traditional Chinese medicine (TCM) is a promising and effective treatment for cancer with minimal side effects through a multi-active ingredient multitarget network. Radix Bupleuri and Rhizoma Cyperi are listed as herbs dispersing stagnated liver Qi in China. They have been used clinically to treat liver diseases for many years and recent pharmacological studies have shown that they inhibit the proliferation of hepatocellular carcinoma (HCC). However, the pharmacological mechanisms, potential targets, and clinical value of the Radix Bupleuri-Rhizoma Cyperi herb pair (CXP) for suppressing HCC growth have not been fully elucidated. We identified 44 CXP targets involved in the treatment of HCC using the GEO dataset and HERB database. An analysis of the Traditional Chinese Medicine System Pharmacology Database (TCMSP) showed that CXP exerts synergistic effects through 4 active ingredients, including quercetin, stigmasterol, isorhamnetin, and kaempferol. GO and KEGG analyses revealed that CXP mainly regulates HCC progression through metabolic pathways, the p53 signaling pathway, and the cell cycle. Additionally, we applied The Cancer Genome Atlas (TCGA)-liver hepatocellular carcinoma (LIHC) database to perform the expression patterns, clinical features, and prognosis of 6 genes (CCNB1, CDK1, CDK4, MYC, CDKN2A, and CHEK1) in cell cycle pathways to reveal that CXP suppresses HCC clinical therapeutic value. Moreover, based on molecular docking, we further verified that CXP exerts its anti-HCC activity through the interaction of multiple active components with cell cycle-related genes. We systematically revealed the potential pharmacological mechanisms and targets of CXP in HCC using multilevel data integration and molecular docking strategies.

A Ferroptosis-Related LncRNA Signature Associated with Prognosis, Tumor Immune Environment, and Genome Instability in Hepatocellular Carcinoma

Background

Ferroptosis is an iron-dependent form of cell death. In this study, we identified ferroptosis-related long noncoding RNAs (FRlncRNAs) to investigate their association with hepatocellular carcinoma (HCC) in prognosis, tumor immune environment, and genome instability.

Methods

Transcriptome profile data of HCC were retrieved from a public database. FRlncRNAs were identified by co-expression analysis. Patients were randomly divided into training and test cohorts. Univariate Cox analysis and Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression were performed to construct a risk model. Patients were divided into high- and low-risk groups based on the risk model. AUC and C index were used to assess the risk model. Survival analysis, immune status, and genome instability were compared between the two groups.

Results

Sixteen FRlncRNAs were identified and used to construct an FRlncRNA signature for the risk model. The Kaplan-Meier analysis revealed that patients in the high-risk group had poorer overall survival than patients in the low-risk group. The area under curve of the risk model was 0.879, 0.809, and 0.757 in the training cohort and 0.635, 0.688, and 0.739 in the test cohort at 1, 2, and 3 years, respectively. The risk model was an independent prognostic predictor and showed excellent prediction of prognosis compared with clinicopathological features. For the differentially expressed ferroptosis-related genes, many enriched metabolic pathways were identified in the functional enrichment analysis. Immune cells such as CD8+ T cells, macrophages M1, natural killer cells, and B cells, which may be associated with antitumor immune responses, differed between the high- and low-risk groups. Genome instability based on the risk model was also explored. A total of 61 genes were differently mutated between the two risk groups, and among them, TP53, HECW2, TRIM66, MCTP2, and KIAA1551 had the most significant mutation frequency differences.

Conclusion

The FRlncRNA signature is closely related with overall survival, tumor immune environment, and genome instability in HCC.

In silico, in vitro screening of antioxidant and anticancer potentials of bioactive secondary metabolites from an endophytic fungus (Curvularia sp.) from Phyllanthus niruri L

The main objective of this research work is to discover novel and efficient phytochemical substances from endophytic fungus found in medicinal plants. Curvularia geniculata L. (C. geniculata L.), an endophytic fungus isolated from Phyllanthus niruri L. (P. niruri L.), was tested against hepatoma cell lines (HepG2) in order to screen their antioxidant and anticancer potentials. The profiling of phytochemicals from the fungal extract was characterized using gas chromatography-mass spectrometry (GC-MS), and molecular docking was done for the identified compounds against one of the potential receptors predominantly present in the hepatocellular carcinoma cell lines. Among the phytochemicals found, 2-methyl-7-phenylindole had the highest binding affinity (- 8.8 kcal mol^-1) for the epidermal growth factor receptor (EGFR). The stability of 2-methyl-7-phenylindole in the EGFR-binding pockets was tested using in silico molecular dynamics simulation. The fungal extract showed the highest antioxidant activity as measured by DPPH, ABTS radical scavenging, and FRAP assays. In vitro cytotoxicity assay of fungal extract demonstrated the concentration-dependent cytotoxicity against HepG2 cells after 24 h, and the IC₅₀ (50% cell death) value was estimated to be 62.23 μg mL^-1. Typical morphological changes such as condensation of nuclei and deformed membrane structures are indicative of ongoing apoptosis. The mitochondria of HepG2 cells were also targeted by the endophytic fungal extract, which resulted in substantial generation of reactive oxygen species (ROS) leading to the destruction of mitochondrial transmembrane potential integrity. These outcomes suggest that the ethyl acetate extract of C. geniculata L. has the potential to be an antioxidant agent and further to be exploited in developing potential anticancer agents.

Aidi injection altered the activity of CYP2D4, CYP1A2, CYP2C19, CYP3A2, CYP2E1 and CYP2C11 in normal and diethylnitrosamine-induced hepatocellular carcinoma in rats

ETHNOPHARMACOLOGICAL RELEVANCE

Aidi injection (ADI), a traditional chinese medicine preparation, is widely used in combination with chemotherapy for the treatment of various malignant tumors, such as hepatocellular carcinoma (HCC). Studies have shown that changes in cytochrome P450 (CYP450) activity in disease states would affect the metabolism of drugs in vivo, especially liver diseases. However, the changes of Aidi injection on the activities of CYP2D4, CYP1A2, CYP2C19, CYP3A2, CYP2E1 and CYP2C11 in normal and HCC states are still unknown.

AIM OF THE STUDY

The cocktail probe drugs method was used to investigate the effects of ADI on the activity of CYP2D4, CYP1A2, CYP2C19, CYP3A2, CYP2E1 and CYP2C11 in normal and HCC rats.

MATERIALS AND METHODS

The HCC rats was induced by diethylnitrosamine (DEN). Then, both normal and HCC rats were randomly divided into 2 groups (n = 6). They were given saline or ADI (10 mL/kg/d, i.p) for 2 weeks, respectively. On the fifteenth day, cocktail probe mixing solution, including metoprolol (10 mg/kg), caffeine (1.0 mg/kg), omeprazole (2.0 mg/kg), midazolam (2.0 mg/kg), chlorzoxazone (4.0 mg/kg) and tolbutamide (0.5 mg/kg), was injected into tail vein of all rats in each group. The blood sample was obtained at specified time. After the protein is precipitated, six probe drugs are analyzed by ultra performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS).

RESULTS

Compared with control group, the activity of CYP3A2 and CYP2E1 was significantly lower in the ADI group. Compared with the model group, the activities of CYP1A2, CYP3A2, CYP2E1, and CYP2C11 enzymes in the ADI model group were significantly reduced. Additionally, the activity of CYP2D4, CYP1A2, CYP2C19, CYP3A2, CYP2E1 and CYP2C11 enzymes in model group was significantly lower than control group.

CONCLUSIONS

ADI can inhibit a lot of CYP450 enzyme, so it may reduce the dosage of chemotherapeutic drugs to reach the required plasma concentration of chemotherapeutic drugs, which is of great significance for the combination of anti-tumor chemotherapeutic drugs and is worthy of further in-depth study and clinical attention.

Acyl-CoA Binding Domain Containing 4 Polymorphism rs4986172 and Expression Can Serve as Overall Survival Biomarkers for Hepatitis B Virus-Related Hepatocellular Carcinoma Patients After Hepatectomy

Background

The aim of our study was to evaluate the potential of expression and single nucleotide polymorphism of Acyl-CoA binding domain containing 4 (ACBD4) gene as prognosis biomarkers in patients with hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC) after hepatectomy.

Methods

HBV-related HCC patients from the First Affiliated Hospital of Guangxi Medical University and GSE14520 were included in the current study, as well as The Cancer Genome Atlas (TCGA) HCC verification cohort. Prognostic analysis and multiple functional enrichment analysis methods were used to evaluate the prognostic value and potential biological functions of the ACBD4 gene in HBV-related HCC.

Results

We found that ACBD4 gene is highly expressed in normal liver tissues and markedly down-regulated in HBV-related HCC tissues. ACBD4 gene was significantly related to overall survival (OS) of HCC in TCGA and GSE14520 cohorts, and patients with low ACBD4 expression were markedly related to poor OS. Rs4986172 was observed as an OS biomarker after hepatectomy in the Guangxi HBV-related HCC cohort. The OS of rs4986172 GG genotype was worse than that of HCC patients with A allele (AA and AG genotypes). Multifunctional enrichment analysis suggested that ACBD4 gene is closely related to the metabolic, peroxisome proliferator-activated receptor and cytochrome P450 pathway. Through connectivity map, we also identified eight compounds that may be used as targeted therapeutic agents for ACBD4 gene in HBV-related HCC; these compounds were scopoletin, alfaxalone, bephenium hydroxynaphthoate, apramycin, 4,5-dianilinophthalimide, DL-thiorphan, aminohippuric acid and quinidine. Immune microenvironment analysis revealed that there were significant differences in immune scores of HBV-related HCC tumor tissues with different ACBD4 expression levels.

Conclusion

Our study reveals that ACBD4 expression and rs4986172 can be serve as biomarkers of OS in HBV-related HCC patients after hepatectomy.

Single-Cell Sequencing of Hepatocellular Carcinoma Reveals Cell Interactions and Cell Heterogeneity in the Microenvironment

BACKGROUND

Hepatocellular carcinoma (HCC) is the main histological subtype of liver cancer, which has the characteristics of poor prognosis and high fatality rate. Single-cell sequencing can provide quantitative and unbiased characterization of cell heterogeneity by analyzing the molecular profile of the whole genome of thousands of single cells. Thus, the purpose of this study was to identify novel prognostic markers for HCC based on single-cell sequencing data.

METHODS

Single-cell sequencing of 21 HCC samples and 256 normal liver tissue samples in the GSE124395 dataset was collected from the Gene Expression Omnibus (GEO) database. The quality-controlled cells were grouped by unsupervised cluster analysis and identified the marker genes of each cell cluster. Hereafter, these cell clusters were annotated by singleR and CellMarker according to the expression patterns of the marker genes. Pseudotime analysis was performed to construct the trajectory of cell evolution and to define hub genes in the evolution process. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were used to explore the potential regulatory mechanism of hub genes in HCC. Next, the differential expression of hub genes and the correlation of the expression of these genes with patients' survival and diagnosis were investigated in The Cancer Genome Atlas (TCGA) database.

RESULTS

A total of 9 clusters corresponding to 9 cell types, including NKT cells, hepatocytes, endothelial cells, Kupffer cells, EPCAM+ cells, cancer cells, plasma cells (B cells), immature B cells, and myofibroblasts were identified. We screened 63 key genes related to cell differentiation through trajectory analysis, which were enriched in the process of coagulation. Ultimately, we identified 10 survival-related hub genes in the TCGA database, namely ALDOB, APOC3, APOH, CYP2E1, CYP3A4, GC, HRG, LINC01554, PDK4, and TXN.

CONCLUSION

In conclusion, ALDOB, APOC3, APOH, CYP2E1, CYP3A4, GC, HRG, LINC01554, PDK4, and TXN may serve as hub genes in the diagnosis and prognosis for HCC.

Identification and validation of ferroptosis-associated gene-based on immune score as prognosis markers for hepatocellular carcinoma patients

Background

Ferroptosis has been found to affect the prognosis and immunotherapy of hepatocellular carcinoma (HCC). However, the association between ferroptosis-related genes and infiltrating immune cells in tumor immune microenvironment (TIME) has not been fully elucidated. This study aimed at establishing a prediction model for the progression of HCC using ferroptosis-associated genes based on immune score.

Methods

Transcriptomic, mutation and clinicopathological information were downloaded from TCGA and International Cancer Genome Consortium (ICGC) for this study. Construction of the prediction model was done by Lasso regression analysis. Estimation of the clustering ability of the prediction model was done by t-distributed stochastic neighbor embedding (t-SNE) and principal component analysis (PCA) analyses. Assessment of the accuracy of the prediction model was done by receiver operating characteristic (ROC) and Kaplan-Meier curves.

Results

A prediction model was formulated utilizing three ferroptosis-related genes (G6PD, SAT1 and SLC1A5). The model independently predicted the overall survival (OS). Differentially expressed genes (DEGs) linked to based on Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) analyses immune-associated pathways and functions. Single-sample gene set enrichment analysis (ssGSEA) strategy further confirmed the model was related to immune-associated functions as well as immune cell infiltration.

Conclusions

The three ferroptosis-associated gene-based prediction model was good at predicting the OS outcomes of HCC, improve HCC prognostication and treatment in the clinic.

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.

Suspended solids-associated toxicity of hydraulic fracturing flowback and produced water on early life stages of zebrafish (Danio rerio)

Hydraulic fracturing flowback and produced water (HF-FPW), which contains polyaromatic hydrocarbons (PAHs) and numerous other potential contaminants, is a complex wastewater produced during the recovery of tight hydrocarbon resources. Previous studies on HF-FPW have demonstrated various toxicological responses of aquatic organisms as consequences of combined exposure to high salinity, dissolved organic compounds and particle/suspended solids-bound pollutants. Noteworthy is the lack of studies illustrating the potentially toxic effects of the FPW suspended solids (FPW-SS). In this study, we investigated the acute and sublethal toxicity of suspended solids filtered from six authentic FPW sample collected from two fracturing wells, using a sediment contact assay based on early-life stages of zebrafish (Danio rerio). PAHs profiles and acute toxicity tests provided initial information on the toxic potency of the six samples. Upon exposure to sediment mixture at two selected doses (1.6 and 3.1 mg/mL), results showed adverse effects in larval zebrafish, as revealed by increased Ethoxyresorufin-O-deethylase (EROD) activity. Transcriptional alterations were also observed in xenobiotic biotransformation (ahr, pxr, cyp1a, cyp1b1, cyp1c1, cyp1c2, cyp3a65, udpgt1a1, udpgt5g1), antioxidant response (sod1, sod2, gpx1a, gpx1b) and hormone receptor signaling (esr1, esr2a, cyp19a1a, vtg1) genes. The results demonstrated that even separated from the complex aqueous FPW mixture, FPW-SS can induce toxicological responses in aquatic organisms' early life stages. Since FPW-SS could sediment to the bottom of natural wetland acting as a continuous source of contaminants, the current findings imply the likelihood of long-term environmental risks of polluted sediments on aquatic ecosystems due to FPW spills.

Long Non-coding RNA TMEM220-AS1 Suppressed Hepatocellular Carcinoma by Regulating the miR-484/MAGI1 Axis as a Competing Endogenous RNA

Long non-coding RNAs (lncRNAs) have a considerable regulatory influence on multiple biological processes. Nevertheless, the role of TMEM220-AS1 in hepatocellular carcinoma (HCC) remains unclear. We used The Cancer Genome Atlas (TCGA) database to analyze the differentially expressed lncRNAs. qRT-PCR was used to verify the results for a large population. The in vitro effects of TMEM220-AS1 on HCC cells were determined using Cell Counting Kit-8 (CCK-8), 5-ethynyl-2'-deoxyuridine (EdU), flow cytometry, and Transwell assays in HCC cells. We used qRT-PCR and western blotting to identify the epithelial-mesenchymal transition (EMT). Moreover, we performed bioinformatics analysis, western blotting, dual luciferase reporter gene assay, RNA pull-down, and RNA binding protein immunoprecipitation (RIP) to investigate the underlying molecular mechanisms of TMEM220-AS1 function. Finally, the function of TMEM220-AS1 was verified in vivo. The results showed that TMEM220-AS1 was expressed at considerably low levels in HCC. It was demonstrated that malignant phenotypes and EMT of HCC cells were promoted by the knock down of TMEM220-AS1 both in vivo and in vitro. TMEM220-AS1, which was detected primarily in the cytoplasm, functioned as an miRNA sponge to bind miR-484 and promote the level of membrane-associated guanylate kinase, WW, and PDZ domain containing 1 (MAGI1), thereby curbing the malignant phenotypes of HCC cells. In conclusion, low levels of TMEM220-AS1 promote proliferation and metastasis through the miR-484/MAGI1 axis in HCC.

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.

Whole-transcriptome and proteome analyses identify key differentially expressed mRNAs, miRNAs, lncRNAs and circRNAs associated with HCC

Hepatocellular carcinoma (HCC) is the most common subtype of primary liver cancer and one of the leading causes of cancer-related death worldwide. To gain more insights into the transcriptomic landscape and molecular mechanism of HCC, we performed TMT-labelled tandem mass spectrometry (n = 4) and whole-transcriptome sequencing (n = 3) based on HCC tumour (T) and adjacent normal (N) tissues from seven HCC patients. To comprehensively evaluate the gene-regulatory circuits in HCC, differential expression and enrichment analyses were performed on the differentially expressed proteins (DEPs), genes (DEGs), miRNAs (555), lncRNAs (29) and circRNAs (895). A total of 977 proteins and 243 genes were found to be differentially expressed in HCC tumours compared with adjacent normal tissues. HCC data from The Cancer Genome Atlas were used to validate the results. Combined with the results above, 56 DEP-DEGs with common changes in relative quantity were identified. Functional pathway analysis showed that the DEP-DEGs were mainly enriched in the spliceosome and various metabolic processes. Bioinformatics analysis showed that hsa-miR-1266-5p, hsa-miR-128-1-5p, hsa-miR-139-5p, hsa-miR-34b-3p and hsa-miR-570-3p were involved in the regulation of the hub genes mentioned above. The crucial coexpression (lncRNA-mRNA, circRNA-mRNA) and competing endogenous RNA interaction axes showed the possible functions of the lncRNAs and circRNAs. We explored potential cancer biomarkers by combining proteomic and transcriptomic studies. Our study provides a valuable resource for understanding regulatory mechanisms at the RNA level and may ultimately further assist in the development of diagnostic and/or therapeutic targets for HCC.

Identification and validation of ADME genes as prognosis and therapy markers for hepatocellular carcinoma patients

Purpose: ADME genes are genes involved in drug absorption, distribution, metabolism, and excretion (ADME). Previous studies report that expression levels of ADME-related genes correlate with prognosis of hepatocellular carcinoma (HCC) patients. However, the role of ADME gene expression on HCC prognosis has not been fully explored. The present study sought to construct a prediction model using ADME-related genes for prognosis of HCC.

Methods: Transcriptome and clinical data were retrieved from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC), which were used as training and validation cohorts, respectively. A prediction model was constructed using univariate Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO) analysis. Patients were divided into high- and low-risk groups based on the median risk score. The predictive ability of the risk signature was estimated through bioinformatics analyses.

Results: Six ADME-related genes (CYP2C9, ABCB6, ABCC5, ADH4, DHRS13, and SLCO2A1) were used to construct the prediction model with a good predictive ability. Univariate and multivariate Cox regression analyses showed the risk signature was an independent predictor of overall survival (OS). A single-sample gene set enrichment analysis (ssGSEA) strategy showed a significant relationship between risk signature and immune status. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed differentially expressed genes (DEGs) in the high- and low-risk groups were enriched in biological process (BP) associated with metabolic and cell cycle pathways.

Conclusion: A prediction model was constructed using six ADME-related genes for prediction of HCC prognosis. This signature can be used to improve HCC diagnosis, treatment, and prognosis in clinical use.

Six Metabolism Related mRNAs Predict the Prognosis of Patients With Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is one of the most common aggressive solid malignant tumors and current research regards HCC as a type of metabolic disease. This study aims to establish a metabolism-related mRNA signature model for risk assessment and prognosis prediction in HCC patients. Methods: HCC data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC) and Gene Enrichment Analysis (GSEA) website. Least absolute shrinkage and selection operator (LASSO) was used to screen out the candidate mRNAs and calculate the risk coefficient to establish the prognosis model. A high-risk group and low-risk group were separated for further study depending on their median risk score. The reliability of the prediction was evaluated in the validation cohort and the whole cohort. Results: A total of 548 differential mRNAs were identified from HCC samples (n = 374) and normal controls (n = 50), 45 of which were correlated with prognosis. A total of 373 samples met the screening criteria and there were randomly divided into the training cohort (n = 186) and the validation cohort (n = 187). In the training cohort, six metabolism-related mRNAs were used to construct a prognostic model with a LASSO regression model. Based on the risk model, the overall survival rate of the high-risk cohort was significantly lower than that of the low-risk cohort. The results of a time-ROC curve proved that the risk score (AUC = 0.849) had a higher prognostic value than the pathological grade, clinical stage, age or gender. Conclusion: The model constructed by the six metabolism-related mRNAs has a significant value for survival prediction and can be applied to guide the evaluation of HCC and the designation of clinical therapy.

Parameterization of Microsomal and Cytosolic Scaling Factors: Methodological and Biological Considerations for Scalar Derivation and Validation

Mathematical models that can predict the kinetics of compounds have been increasingly adopted for drug development and risk assessment. Data for these models may be generated from in vitro experimental systems containing enzymes contributing to metabolic clearance, such as subcellular tissue fractions including microsomes and cytosol. Extrapolation from these systems is facilitated by common scaling factors, known as microsomal protein per gram (MPPG) and cytosolic protein per gram (CPPG). Historically, parameterization of MPPG and CPPG has employed the use of recovery factors, commonly benchmarked to cytochromes P450 which work well in some contexts, but could be problematic for other enzymes. Here, we propose absolute quantification of protein content and supplementary assays to evaluate microsomal/cytosolic purity that should be employed. Examples include calculation of microsomal latency by mannose-6-phosphatase activity and immunoblotting of subcellular fractions with fraction-specific markers. Further considerations include tissue source, as disease states can affect enzyme expression and activity, and the methodology used for scalar parameterization. Regional- and organ-specific expression of enzymes, in addition to differences in organ physiology, is another important consideration. Because most efforts have focused on the liver that is, for the most part, homogeneous, derived scalars may not capture the heterogeneity of other major tissues contributing to xenobiotic metabolism including the kidneys and small intestine. Better understanding of these scalars, and how to appropriately derive them from extrahepatic tissues can provide support to the inferences made with physiologically based pharmacokinetic modeling, increase its accuracy in characterizing in vivo drug pharmacokinetics, and improve confidence in go-no-go decisions for clinical trials.

Pharmacogenomics to Predict Tumor Therapy Response: A Focus on ATP-Binding Cassette Transporters and Cytochromes P450

Pharmacogenomics is an evolving tool of precision medicine. Recently, due to the introduction of next-generation sequencing and projects generating "Big Data", a plethora of new genetic variants in pharmacogenes have been discovered. Cancer resistance is a major complication often preventing successful anticancer treatments. Pharmacogenomics of both somatic mutations in tumor cells and germline variants may help optimize targeted treatments and improve the response to conventional oncological therapy. In addition, integrative approaches combining copy number variations and long noncoding RNA profiling with germline and somatic variations seem to be a promising approach as well. In pharmacology, expression and enzyme activity are traditionally the more studied aspects of ATP-binding cassette transporters and cytochromes P450. In this review, we briefly introduce the field of pharmacogenomics and the advancements driven by next-generation sequencing and outline the possible roles of genetic variation in the two large pharmacogene superfamilies. Although the evidence needs further substantiation, somatic and copy number variants as well as rare variants and common polymorphisms in these genes could all affect response to cancer therapy. Regulation by long noncoding RNAs has also been shown to play a role. However, in all these areas, more comprehensive studies on larger sets of patients are needed.