A human-specific switch of alternatively spliced AFMID isoforms contributes to TP53 mutations and tumor recurrence in hepatocellular carcinoma

A Review of the Evidence for Tryptophan and the Kynurenine Pathway as a Regulator of Stem Cell Niches in Health and Disease

Stem cells are ubiquitously found in various tissues and organs in the body, and underpin the body's ability to repair itself following injury or disease initiation, though repair can sometimes be compromised. Understanding how stem cells are produced, and functional signaling systems between different niches is critical to understanding the potential use of stem cells in regenerative medicine. In this context, this review considers kynurenine pathway (KP) metabolism in multipotent adult progenitor cells, embryonic, haematopoietic, neural, cancer, cardiac and induced pluripotent stem cells, endothelial progenitor cells, and mesenchymal stromal cells. The KP is the major enzymatic pathway for sequentially catabolising the essential amino acid tryptophan (TRP), resulting in key metabolites including kynurenine, kynurenic acid, and quinolinic acid (QUIN). QUIN metabolism transitions into the adjoining de novo pathway for nicotinamide adenine dinucleotide (NAD) production, a critical cofactor in many fundamental cellular biochemical pathways. How stem cells uptake and utilise TRP varies between different species and stem cell types, because of their expression of transporters and responses to inflammatory cytokines. Several KP metabolites are physiologically active, with either beneficial or detrimental outcomes, and evidence of this is presented relating to several stem cell types, which is important as they may exert a significant impact on surrounding differentiated cells, particularly if they metabolise or secrete metabolites differently. Interferon-gamma (IFN-γ) in mesenchymal stromal cells, for instance, highly upregulates rate-limiting enzyme indoleamine-2,3-dioxygenase (IDO-1), initiating TRP depletion and production of metabolites including kynurenine/kynurenic acid, known agonists of the Aryl hydrocarbon receptor (AhR) transcription factor. AhR transcriptionally regulates an immunosuppressive phenotype, making them attractive for regenerative therapy. We also draw attention to important gaps in knowledge for future studies, which will underpin future application for stem cell-based cellular therapies or optimising drugs which can modulate the KP in innate stem cell populations, for disease treatment.

Decoding the role of aberrant RNA alternative splicing in hepatocellular carcinoma: a comprehensive review

During eukaryotic gene expression, alternative splicing of messenger RNA precursors is critical in increasing protein diversity and regulatory complexity. Multiple transcript isoforms could be produced by alternative splicing from a single gene; they could eventually be translated into protein isoforms with deleted, added, or altered domains or produce transcripts containing premature termination codons that could be targeted by nonsense-mediated mRNA decay. Alternative splicing can generate proteins with similar, different, or even opposite functions. Increasingly strong evidence indicates that abnormal RNA splicing is a prevalent and crucial occurrence in cellular differentiation, tissue advancement, and the development and progression of cancer. Aberrant alternative splicing could affect cancer cell activities such as growth, apoptosis, invasiveness, drug resistance, angiogenesis, and metabolism. This systematic review provides a comprehensive overview of the impact of abnormal RNA alternative splicing on the development and progression of hepatocellular carcinoma.

Adhesion G Protein-Coupled Receptor G2 Promotes Hepatocellular Carcinoma Progression and Serves as a Neutrophil-Related Prognostic Biomarker

Adhesion G protein-coupled receptor G2 (ADGRG2) is an orphan adhesion G protein-coupled receptor (GPCR), which performs a tumor-promoting role in certain cancers; however, it has not been systematically investigated in hepatocellular carcinoma (HCC). In the current study, we utilized multiple databases to analyze the expression and diagnostic and prognostic value of ADGRG2 in HCC and its correlation with immune infiltration and inflammatory factors. The function and upstream regulatory miRNA of ADGRG2 were validated through qPCR, Western blot, CCK8, wound healing, and dual luciferase assays. It turned out that ADGRG2 was significantly higher in HCC and had a poor survival rate, especially in AFP ≤ 400 ng/mL subgroups. Functional enrichment analysis suggested that ADGRG2 may be involved in cancer pathways and immune-related pathways. In vitro, siRNA-mediated ADGRG2 silencing could inhibit the proliferation and migration of Huh7 and HepG2 cells. There was a highly significant positive correlation between ADGRG2 and neutrophils. Moreover, NET-related genes were filtered and confirmed, such as ENO1 and S100A9. Meanwhile, the high expression of ADGRG2 was also accompanied by the highest number of inflammatory cytokines, chemokines, and chemokine receptors and good immunotherapy efficacy. Finally, AGDGR2 may be sensitive to two drugs (PIK-93 and NPK76-II-72-1) and can be targeted by miR-326. In conclusion, ADGRG2 may serve as a novel biomarker and drug target for HCC diagnosis, immunotherapy, and prognosis and was related to neutrophils and the inflammatory process of liver cancer development.

The kynurenine pathway presents multi-faceted metabolic vulnerabilities in cancer

The kynurenine pathway (KP) and associated catabolites play key roles in promoting tumor progression and modulating the host anti-tumor immune response. To date, considerable focus has been on the role of indoleamine 2,3-dioxygenase 1 (IDO1) and its catabolite, kynurenine (Kyn). However, increasing evidence has demonstrated that downstream KP enzymes and their associated metabolite products can also elicit tumor-microenvironment immune suppression. These advancements in our understanding of the tumor promotive role of the KP have led to the conception of novel therapeutic strategies to target the KP pathway for anti-cancer effects and reversal of immune escape. This review aims to 1) highlight the known biological functions of key enzymes in the KP, and 2) provide a comprehensive overview of existing and emerging therapies aimed at targeting discrete enzymes in the KP for anti-cancer treatment.

Upregulation of hsa_circ_0002003 promotes hepatocellular carcinoma progression

BACKGROUND

Circular RNAs (circRNAs), which are involved in various human malignancies, have emerged as promising biomarkers. The present study aimed to investigate unique expression profiles of circRNAs in hepatocellular carcinoma (HCC) and identify novel biomarkers associated with HCC development and progression.

METHODS

CircRNA expression profiles of HCC tissues were jointly analyzed to identify differentially expressed circRNAs. Overexpression plasmid and siRNA targeting candidate circRNAs were used in functional assays in vitro. CircRNA-miRNA interactions were predicted using miRNAs expressed in the miRNA-seq dataset GSE76903. To further screen downstream genes targeted by the miRNAs, survival analysis and qRT-PCR were conducted to evaluate their prognostic role in HCC and construct a ceRNA regulatory network.

RESULTS

Three significantly upregulated circRNAs, hsa_circ_0002003, hsa_circ_0002454, and hsa_circ_0001394, and one significantly downregulated circRNA, hsa_circ_0003239, were identified and validated by qRT-PCR. Our in vitro data indicated that upregulation of hsa_circ_0002003 accelerated cell growth and metastasis. Mechanistically, DTYMK, DAP3, and STMN1, which were targeted by hsa-miR-1343-3p, were significantly downregulated in HCC cells when hsa_circ_0002003 was silenced and were significantly correlated with poor prognosis in patients with HCC.

CONCLUSION

Hsa_circ_0002003 may play critical roles in HCC pathogenesis and serve as a potential prognostic biomarker for HCC. Targeting the hsa_circ_0002003/hsa-miR-1343-3p/STMN1 regulatory axis could be an effective therapeutic strategy in patients with HCC.

Research on Predicting the Occurrence of Hepatocellular Carcinoma Based on Notch Signal-Related Genes Using Machine Learning Algorithms

BACKGROUND/AIMS

Hepatocellular carcinoma, a highly malignant tumor, is difficult to diagnose, treat, and predict the prognosis. Notch signaling pathway can affect hepatocellular carcinoma. We aimed to predict the occurrence of hepatocellular carcinoma based on Notch signal-related genes using machine learning algorithms.

MATERIALS AND METHODS

We downloaded hepatocellular carcinoma data from the Cancer Genome Atlas and Gene Expression Omnibus databases and used machine learning methods to screen the hub Notch signal-related genes. Machine learning classification was used to construct a prediction model for the classification and diagnosis of hepatocellular carcinoma cancer. Bioinformatics methods were applied to explore the expression of these hub genes in the hepatocellular carcinoma tumor immune microenvironment.

RESULTS

We identified 4 hub genes, namely, LAMA4, POLA2, RAD51, and TYMS, which were used as the final variables, and found that AdaBoostClassifie was the best algorithm for the classification and diagnosis model of hepatocellular carcinoma. The area under curve, accuracy, sensitivity, specificity, positive predictive value, negative predictive value, and F1 score of this model in the training set were 0.976, 0.881, 0.877, 0.977, 0.996, 0.500, and 0.932; respectively. The area under curves were 0.934, 0.863, 0.881, 0.886, 0.981, 0.489, and 0.926. The area under curve in the external validation set was 0.934. Immune cell infiltration was related to the expression of 4 hub genes. Patients in the low-risk group of hepatocellular carcinoma were more likely to have an immune escape.

CONCLUSION

The Notch signaling pathway was closely related to the occurrence and development of hepatocellular carcinoma. The hepatocellular carcinoma classification and diagnosis model established based on this had a high degree of reliability and stability.

Comprehensive analysis of full-length transcripts reveals novel splicing abnormalities and oncogenic transcripts in liver cancer

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. Although long-reads sequencing technologies would enable the sequencing of full-length transcripts, the data analysis is difficult. In this study, we developed an analysis pipeline named SPLICE and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC) and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 46,663 transcripts from the protein-coding genes in the HCCs and the matched non-cancerous livers, of which 5,366 (11.5%) were novel. A comparison of expression levels identified 9,933 differentially expressed transcripts (DETs) in 4,744 genes. Interestingly, 746 genes with DETs, including the LINE1-MET transcript, were not found by a gene-level analysis. We also found that fusion transcripts of transposable elements and hepatitis B virus (HBV) were overexpressed in HCCs. In vitro experiments on DETs showed that LINE1-MET and HBV-human transposable elements promoted cell growth. Furthermore, fusion gene detection showed novel recurrent fusion events that were not detected in the short-reads. These results suggest the efficiency of full-length transcriptome studies and the importance of splicing variants in carcinogenesis.

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. In this study, we developed an analysis pipeline named SPLICE for long-read transcriptome sequencing and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC), and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 5,366 novel transcripts and 9,933 differentially expressed transcripts in 4,744 genes between HCCs and non-cancerous livers. An analysis of hepatitis B virus (HBV) transcripts showed that fusion transcripts of the HBV gene and human transposable elements were overexpressed in HBV-infected HCCs. We also identified fusion genes that were not found in the short-reads. These results suggest that long-reads sequencing technologies provide a fuller understanding of cancer transcripts and that our method contributes to the analysis of transcriptome sequences by such technologies.

The Kynurenine Pathway and Cancer: Why Keep It Simple When You Can Make It Complicated

Simple Summary

The kynurenine pathway has two main physiological roles: (i) it protects specific organs such as the eyes and placenta from strong immune reactions and (ii) it additionally generate in the liver and kidney a metabolite essential to all cells of human body. Abnormal activation of this pathway is recurrently observed in numerous cancer types. Its two functions are hijacked to promote tumor growth and cancer cell dissemination through multiple mechanisms. Clinical assays including administration of inhibitors of this pathway have not yet been successful. The complex regulation of this pathway is likely the reason behind this failure. In this review, we try to give an overview of the current knowledge about this pathway, to point out the next challenges, and to propose alternative therapeutic routes.

Abstract

The kynurenine pathway has been highlighted as a gatekeeper of immune-privileged sites through its ability to generate from tryptophan a set of immunosuppressive metabolic intermediates. It additionally constitutes an important source of cellular NAD⁺ for the organism. Hijacking of its immunosuppressive functions, as recurrently observed in multiple cancers, facilitates immune evasion and promotes tumor development. Based on these observations, researchers have focused on characterizing indoleamine 2,3-dioxygenase (IDO1), the main enzyme catalyzing the first and limiting step of the pathway, and on developing therapies targeting it. Unfortunately, clinical trials studying IDO1 inhibitors have thus far not met expectations, highlighting the need to unravel this complex signaling pathway further. Recent advances demonstrate that these metabolites additionally promote tumor growth, metastatic dissemination and chemoresistance by a combination of paracrine and autocrine effects. Production of NAD⁺ also contributes to cancer progression by providing cancer cells with enhanced plasticity, invasive properties and chemoresistance. A comprehensive survey of this complexity is challenging but necessary to achieve medical success.

The oncogenic PI3K-induced transcriptomic landscape reveals key functions in splicing and gene expression regulation

The PI3K pathway regulates proliferation, survival, and metabolism and is frequently activated across human cancers. A comprehensive elucidation of how this signaling pathway controls transcriptional and co-transcriptional processes could provide new insights into the key functions of PI3K signaling in cancer. Here, we undertook a transcriptomic approach to investigate genome-wide gene expression and transcription factor (TF) activity changes, as well as splicing and isoform usage dynamics, downstream of PI3K. These analyses uncovered widespread alternatively spliced (AS) isoforms linked to proliferation, metabolism, and splicing in PIK3CA mutant cells, which were reversed by inhibition of PI3Kα. Analysis of paired tumor biopsies from PIK3CA-mutated breast cancer patients undergoing treatment with PI3Kα inhibitors identified widespread splicing alterations that affect specific isoforms in common with the preclinical models, and these alterations, namely PTK2/FRNK and AFMID isoforms, were validated as functional drivers of cancer cell growth or migration. Mechanistically, isoform-specific splicing factors mediated PI3K-dependent RNA splicing. Treatment with splicing inhibitors rendered breast cancer cells more sensitive to the PI3Kα inhibitor alpelisib, resulting in greater growth inhibition than alpelisib alone. This study provides the first comprehensive analysis of widespread splicing alterations driven by oncogenic PI3K in breast cancer. The atlas of PI3K-mediated splicing programs establishes a key role for the PI3K pathway in regulating splicing, opening new avenues for exploiting PI3K signaling as a therapeutic vulnerability in breast cancer.

An Alternatively Spliced Variant of METTL3 Mediates Tumor Suppression in Hepatocellular Carcinoma

Many post-transcriptional mRNA processing steps play crucial roles in tumorigenesis and the progression of cancers, such as N6-methyladenosine (m⁶A) modification and alternative splicing. Upregulation of methyltransferase-like 3 (METTL3), the catalytic core of the m⁶A methyltransferase complex, increases m⁶A levels and results in significant effects on the progression of hepatocellular carcinoma (HCC). However, alternative splicing of METTL3 has not been fully investigated, and the functions of its splice variants remain unclear. Here, we analyzed both our and online transcriptomic data, obtaining 13 splice variants of METTL3 in addition to canonical full-length METTL3-A in HCC cell lines and tissues. Validated by RT–qPCR and Western blotting, we found that METTL3-D, one of the splice variants expressing a truncated METTL3 protein, exhibits higher levels than METTL3-A in normal human livers but lower levels than METTL3-A in HCC tumor tissues and cell lines. Further functional assays demonstrated that METTL3-D expression decreased cellular m⁶A modification, inhibited the proliferation, migration, and invasion of HCC cells, and was negatively associated with the malignancy of patient tumors, exhibiting functions opposite to those of full-length METTL3-A. This study demonstrates that the METTL3-D splice variant is a tumor suppressor that could potentially be used as a target for HCC therapy.

Comprehensive Analysis of Alternative Splicing in Gastric Cancer Identifies Epithelial-Mesenchymal Transition Subtypes Associated with Survival

Alternatively spliced RNA isoforms are a hallmark of tumors, but their nature, prevalence, and clinical implications in gastric cancer have not been comprehensively characterized. We systematically profiled the splicing landscape of 83 gastric tumors and matched normal mucosa, identifying and experimentally validating eight splicing events that can classify all gastric cancers into three subtypes: epithelial-splicing (EpiS), mesenchymal-splicing (MesS), and hybrid-splicing. These subtypes were associated with distinct molecular signatures and epithelial-mesenchymal transition markers. Subtype-specific splicing events were enriched in motifs for splicing factors RBM24 and ESRP1, which were upregulated in MesS and EpiS tumors, respectively. A simple classifier based only on RNA levels of RBM24 and ESRP1, which can be readily implemented in the clinic, was sufficient to distinguish gastric cancer subtypes and predict patient survival in multiple independent patient cohorts. Overall, this study provides insights into alternative splicing in gastric cancer and the potential clinical utility of splicing-based patient classification.

SIGNIFICANCE

This study presents a comprehensive analysis of alternative splicing in the context of patient classification, molecular mechanisms, and prognosis in gastric cancer.

Comprehensive analysis and establishment of a prediction model of alternative splicing events reveal the prognostic predictor and immune microenvironment signatures in triple negative breast cancer

Background

Triple-negative breast cancer (TNBC) is widely concerning because of high malignancy and poor prognosis. There is increasing evidence that alternative splicing (AS) plays an important role in the development of cancer and the formation of the tumour microenvironment. However, comprehensive analysis of AS signalling in TNBC is still lacking and urgently needed.

Methods

Transcriptome and clinical data of 169 TNBC tissues and 15 normal tissues were obtained and integrated from the cancer genome atlas (TCGA), and an overview of AS events was downloaded from the SpliceSeq database. Then, differential comparative analysis was performed to obtain cancer-associated AS events (CAAS). Metascape was used to perform parent gene enrichment analysis based on CAAS. Unsupervised cluster analysis was performed to analyse the characteristics of immune infiltration in the microenvironment. A splicing network was established based on the correlation between CAAS events and splicing factors (SFs). We then constructed prediction models and assessed the accuracy of these models by receiver operating characteristic (ROC) curve and Kaplan–Meier survival analyses. Furthermore, a nomogram was adopted to predict the individualized survival rate of TNBC patients.

Results

We identified 1194 cancer-associated AS events (CAAS) and evaluated the enrichment of 981 parent genes. The top 20 parent genes with significant differences were mostly related to cell adhesion, cell component connection and other pathways. Furthermore, immune-related pathways were also enriched. Unsupervised clustering analysis revealed the heterogeneity of the immune microenvironment in TNBC. The splicing network also suggested an obvious correlation between SFs expression and CAAS events in TNBC patients. Univariate and multivariate Cox regression analyses showed that the survival-related AS events were detected, including some significant participants in the carcinogenic process. A nomogram incorporating risk, AJCC and radiotherapy showed good calibration and moderate discrimination.

Conclusion

Our study revealed AS events related to tumorigenesis and the immune microenvironment, elaborated the potential correlation between SFs and CAAS, established a prognostic model based on survival-related AS events, and created a nomogram to better predict the individual survival rate of TNBC patients, which improved our understanding of the relationship between AS events and TNBC.

Genetic Alterations and Transcriptional Expression of m6A RNA Methylation Regulators Drive a Malignant Phenotype and Have Clinical Prognostic Impact in Hepatocellular Carcinoma

Background: N6-methyladenosine (m⁶A) RNA methylation, associated with cancer initiation and progression, is dynamically regulated by the m⁶A RNA regulators. However, its role in liver carcinogenesis is poorly understood.

Methods: Three hundred seventy-one hepatocellular carcinoma (HCC) patients from The Cancer Genome Atlas database with sequencing and copy number variations/mutations data were included. Survival analysis was performed using Cox regression model. We performed gene set enrichment analysis to explore the functions associated with different HCC groups. Finally, we used a machine-learning model on selected regulators for developing a risk signature (m⁶Ascore) The prognostic value of m⁶Ascore was finally validated in another two GEO datasets.

Results: We demonstrated that 11 m⁶A RNA regulators are significantly differentially expressed among 371 HCC patients stratified by clinicopathological features (P<0.001). We then identified two distinct HCC clusters by applying consensus clustering to m⁶A RNA regulators. Compared with the cluster2 subgroup, the cluster1 subgroup correlates with poorer prognosis (P < 0.001). Moreover, the cell cycle, splicesome and notch signaling pathway are significantly enriched in the cluster1 subgroup. We further derived m⁶Ascore, using four m⁶A regulators, predicting HCC prognosis well at three (AUC = 0.7) or 5 years (AUC=0.7) in validation. The prognostic value of m⁶Ascore also was validated successfully in two GEO datasets (P < 0.05). Finally, we discovered that mutations and copy number variations of m⁶A regulators, conferring worse survival, are strongly associated with TP53 mutations in HCC.

Conclusions: We find a significant relationship between the alterations and different expressions causing increased m⁶A level and worse survival, especially in TP53-mutated HCC patients. Genetic alterations of m⁶A genes might cooperate with TP53 and its regulator targets in the HCC pathogenesis. Our m⁶Ascore may be applied in the clinical trials for patient stratification in HCC.

PSI-Sigma: a comprehensive splicing-detection method for short-read and long-read RNA-seq analysis

MOTIVATION

Percent Spliced-In (PSI) values are commonly used to report alternative pre-mRNA splicing (AS) changes. Previous PSI-detection tools were limited to specific AS events and were evaluated by in silico RNA-seq data. We developed PSI-Sigma, which uses a new PSI index, and we employed actual (non-simulated) RNA-seq data from spliced synthetic genes (RNA Sequins) to benchmark its performance (i.e. precision, recall, false positive rate and correlation) in comparison with three leading tools (rMATS, SUPPA2 and Whippet).

RESULTS

PSI-Sigma outperformed these tools, especially in the case of AS events with multiple alternative exons and intron-retention events. We also briefly evaluated its performance in long-read RNA-seq analysis, by sequencing a mixture of human RNAs and RNA Sequins with nanopore long-read sequencers.

AVAILABILITY AND IMPLEMENTATION

PSI-Sigma is implemented is available at https://github.com/wososa/PSI-Sigma.

Alternative Splicing in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) accounts for the majority of primary liver cancer cases, with more than 850,000 new diagnoses per year globally. Recent trends in the United States have shown that liver cancer mortality has continued to increase in both men and women, while 5-year survival remains below 20%. Understanding key mechanisms that drive chronic liver disease progression to HCC can reveal new therapeutic targets and biomarkers for early detection of HCC. In that regard, many studies have underscored the importance of alternative splicing as a source of novel HCC prognostic markers and disease targets. Alternative splicing of pre-mRNA provides functional diversity to the genome, and endows cells with the ability to rapidly remodel the proteome. Genes that control fundamental processes, such as metabolism, cell proliferation, and apoptosis, are altered globally in HCC by alternative splicing. This review highlights the major splicing factors, RNA binding proteins, transcriptional targets, and signaling pathways that are of key relevance to HCC. We highlight primary research from the past 3-5 years involving functional interrogation of alternative splicing in rodent and human liver, using both large-scale transcriptomic and focused mechanistic approaches. Because this is a rapidly advancing field, we anticipate that it will be transformative for the future of basic liver biology, as well as HCC diagnosis and management.

Whole transcriptome analysis of chemically induced hepatocellular carcinoma using RNA-sequencing analysis

RNA‐sequencing (RNA‐Seq) is a useful method to explore the molecular events in cells and tissues at the transcriptional level. However, comprehensive transcriptome analysis of hepatocarcinogenesis and progression is lacking. In this study, we aimed to characterize a dimethylnitrosamine (DEN) and carbon tetrachloride (CCl₄; DEN+CCl₄)‐induced hepatocellular carcinoma (HCC) mouse model by RNA‐Seq. In total, 2033 genes were up‐regulated and 841 genes were down‐regulated after DEN and CCl₄ stimulation. The differentially expressed genes were highly enriched for the Gene Ontology terms oxoacid metabolic process, carboxylic acid metabolic process, and organic acid metabolic process. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that the top five significantly overrepresented pathways were metabolic pathways, chemical carcinogenesis, steroid hormone biosynthesis, retinol metabolism and metabolism of xenobiotics by cytochrome P450. Moreover, a protein‐protein interaction network analysis indicated that Rous sarcoma oncogene (Src) may play a key role in DEN+CCl₄‐induced HCC. These results provide a comprehensive overview of transcriptome events in DEN+CCl₄‐induced HCC.