A comprehensive performance evaluation, comparison, and integration of computational methods for detecting and estimating cross-contamination of human samples in cancer next-generation sequencing analysis.

Cross-sample contamination is one of the major issues in next-generation sequencing (NGS)-based molecular assays. This type of contamination, even at very low levels, can significantly impact the results of an analysis, especially in the detection of somatic alterations in tumor samples. Several contamination identification tools have been developed and implemented as a crucial quality-control step in the routine NGS bioinformatic pipeline. However, no study has been published to comprehensively and systematically investigate, evaluate, and compare these computational methods in the cancer NGS analysis. In this study, we comprehensively investigated nine state-of-the-art computational methods for detecting cross-sample contamination. To explore their application in cancer NGS analysis, we further compared the performance of five representative tools by qualitative and quantitative analyses using in silico and simulated experimental NGS data. The results showed that Conpair achieved the best performance for identifying contamination and predicting the level of contamination in solid tumors NGS analysis. Moreover, based on Conpair, we developed a Python script, Contamination Source Predictor (ConSPr), to identify the source of contamination. We anticipate that this comprehensive survey and the proposed tool for predicting the source of contamination will assist researchers in selecting appropriate cross-contamination detection tools in cancer NGS analysis and inspire the development of computational methods for detecting sample cross-contamination and identifying its source in the future.

Characterization and mitigation of artifacts derived from NGS library preparation due to structure-specific sequences in the human genome.

BACKGROUND: Hybridization capture-based targeted next generation sequencing (NGS) is gaining importance in routine cancer clinical practice. DNA library preparation is a fundamental step to produce high-quality sequencing data. Numerous unexpected, low variant allele frequency calls were observed in libraries using sonication fragmentation and enzymatic fragmentation. In this study, we investigated the characteristics of the artifact reads induced by sonication and enzymatic fragmentation. We also developed a bioinformatic algorithm to filter these sequencing errors. RESULTS: We used pairwise comparisons of somatic single nucleotide variants (SNVs) and insertions and deletions (indels) of the same tumor DNA samples prepared using both ultrasonic and enzymatic fragmentation protocols. Our analysis revealed that the number of artifact variants was significantly greater in the samples generated using enzymatic fragmentation than using sonication. Most of the artifacts derived from the sonication-treated libraries were chimeric artifact reads containing both cis- and trans-inverted repeat sequences of the genomic DNA. In contrast, chimeric artifact reads of endonuclease-treated libraries contained palindromic sequences with mismatched bases. Based on these distinctive features, we proposed a mechanistic hypothesis model, PDSM (pairing of partial single strands derived from a similar molecule), by which these sequencing errors derive from ultrasonication and enzymatic fragmentation library preparation. We developed a bioinformatic algorithm to generate a custom mutation "blacklist" in the BED region to reduce errors in downstream analyses. CONCLUSIONS: We first proposed a mechanistic hypothesis model (PDSM) of sequencing errors caused by specific structures of inverted repeat sequences and palindromic sequences in the natural genome. This new hypothesis predicts the existence of chimeric reads that could not be explained by previous models, and provides a new direction for further improving NGS analysis accuracy. A bioinformatic algorithm, ArtifactsFinder, was developed and used to reduce the sequencing errors in libraries produced using sonication and enzymatic fragmentation.

Helicobacter pylori CagA-mediated ether lipid biosynthesis promotes ferroptosis susceptibility in gastric cancer.

Helicobacter pylori, particularly cytotoxin-associated gene A (CagA)-positive strains, plays a key role in the progression of gastric cancer (GC). Ferroptosis, associated with lethal lipid peroxidation, has emerged to play an important role in malignant and infectious diseases, but the role of CagA in ferroptosis in cancer cells has not been determined. Here, we report that CagA confers GC cells sensitivity to ferroptosis both in vitro and in vivo. Mechanistically, CagA promotes the synthesis of polyunsaturated ether phospholipids (PUFA-ePLs), which is mediated by increased expression of alkylglycerone phosphate synthase (AGPS) and 1-acylglycerol-3-phosphate O-acyltransferase 3 (AGPAT3), leading to susceptibility to ferroptosis. This susceptibility is mediated by activation of the MEK/ERK/SRF pathway. SRF is a crucial transcription factor that increases AGPS transcription by binding to the AGPS promoter region. Moreover, the results demonstrated that CagA-positive cells are more sensitive to apatinib than are CagA-negative cells, suggesting that detecting the H. pylori CagA status may aid patient stratification for treatment with apatinib.

Survival analysis of comprehensive treatment in Chinese patients with metastatic melanoma: A retrospective analysis.

BACKGROUND: Most of the current progression of immune checkpoint inhibitors for malignant melanoma is based on data from Caucasians in Western countries, but the benefit of Chinese patients is limited, mainly due to different pathological subtypes. The patients in western countries are mainly skin melanoma (about 90%), while the acral and mucosal types are dominant in China, accounting for 41.8% and 22.6% respectively. Acral and mucosal melanoma have lower response rates to immunotherapy and chemotherapy. OBJECTIVE: Whether immune checkpoint inhibitors can improve the survival of Chinese patients with malignant melanoma, therefore, we conducted a retrospective analysis. METHODS: We analyzed 53 patients with metastatic melanoma treated in our hospital to evaluate the efficacy and safety of PD-1 mAb in Chinese patients with metastatic melanoma, and performed univariate and multivariate analyses of prognostic factors that may affect overall survival (OS). RESULTS: In a study of 125 patients with advanced malignant melanoma, 53 patients participated, with a median follow-up of 16 months. Among these, 69.8% died, and 30.2% remained on treatment. Median progression-free survival (PFS) was 6 months, and median OS was 19 months. Patients treated with immune checkpoint inhibitors had improved outcomes, with a median PFS of 7 months and a median OS of 24 months. Patients with bone metastasis and aberrant Lactate dehydrogenase (LDH) post-treatment had worse prognoses, while immunotherapy was a protective factor. Subgroup analysis showed the benefits of immunotherapy across various patient characteristics. No unexpected toxicities were observed. CONCLUSION: The study highlights the efficacy of immune checkpoint inhibitors, particularly PD-1 mAb, in improving survival outcomes for Chinese patients with metastatic melanoma.

OncoCTMiner: streamlining precision oncology trial matching via molecular profile analysis.

By establishing omics sequencing of patient tumors as a crucial element in cancer treatment, the extensive implementation of precision oncology necessitates effective and prompt execution of clinical studies for approving molecular-targeted therapies. However, the substantial volume of patient sequencing data, combined with strict clinical trial criteria, increasingly complicates the process of matching patients to precision oncology studies. To streamline enrollment in these studies, we developed OncoCTMiner, an automated pre-screening platform for molecular cancer clinical trials. Through manual tagging of eligibility criteria for 2227 oncology trials, we identified key bio-concepts such as cancer types, genes, alterations, drugs, biomarkers and therapies. Utilizing this manually annotated corpus along with open-source biomedical natural language processing tools, we trained multiple named entity recognition models specifically designed for precision oncology trials. These models analyzed 460 952 clinical trials, revealing 8.15 million precision medicine concepts, 9.32 million entity-criteria-trial triplets and a comprehensive precision oncology eligibility criteria database. Most significantly, we developed a patient-trial matching system based on cancer patients' clinical and genetic profiles, which can seamlessly integrate with the omics data analysis platform. This system expedites the pre-screening process for potentially suitable precision oncology trials, offering patients swifter access to promising treatment options. Database URL  https://oncoctminer.chosenmedinfo.com.

A unified DNA- and RNA-based NGS strategy for the analysis of multiple types of variants at the dual nucleic acid level in solid tumors.

BACKGROUND: Targeted next-generation sequencing (NGS) is a powerful and suitable approach to comprehensively identify multiple types of variants in tumors. RNA-based NGS is increasingly playing an important role in precision oncology. Both parallel and sequential DNA- and RNA-based approaches are expensive, burdensome, and have long turnaround times, which can be impractical in clinical practice. A streamlined, unified DNA- and RNA-based NGS approach is urgently needed in clinical practice. METHODS: A DNA/RNA co-hybrid capture sequencing (DRCC-Seq) approach was designed to capture pre-capture DNA and RNA libraries in a single tube and convert them into one NGS library. The performance of the DRCC-Seq approach was evaluated by a panel of reference standards and clinical samples. RESULTS: The average depth, DNA data ratio, capture ratio, and target coverage 250 (×) of the DNA panel data had a negative correlation with an increase in the proportion of RNA probes. The SNVs, indels, fusions, and MSI status were not affected by the proportion of RNA probes, but the copy numbers of the target genes were higher than expected in the standard materials, and many unexpected gene amplifications were found using D:R (1:2) and D:R (1:4) probe panels. The optimal ratio of DNA and RNA probes in the combined probe panel was 1:1 using the DRCC-Seq approach. The DRCC-Seq approach was feasible and reliable for detecting multiple types of variants in reference standards and real-world clinical samples. CONCLUSIONS: The DRCC-Seq approach is more cost-effective, with a shorter turnaround time and lower labor requirements than either parallel or sequential targeted DNA NGS and RNA NGS. It is feasible to identify multiple genetic variations at the DNA and RNA levels simultaneously in clinical practice.

Chromosome constitution and genetic relationships of Morus spp. revealed by genomic in situ hybridization.

BACKGROUND: Mulberry (Morus spp.) is an economically important woody plant, which has been used for sericulture (silk farming) for thousands of years. The genetic background of mulberry is complex due to polyploidy and frequent hybridization events. RESULTS: Comparative genomic in situ hybridization (cGISH) and self-GISH were performed to illustrate the chromosome constitution and genetic relationships of 40 mulberry accessions belonging to 12 species and three varietas in the Morus genus and containing eight different ploidy levels. We identified six homozygous cGISH signal patterns and one heterozygous cGISH signal pattern using four genomic DNA probes. Using cGISH and self-GISH data, we defined five mulberry sections (Notabilis, Nigra, Wittiorum, and Cathayana, all contained only one species; and Alba, which contained seven closely related species and three varietas, was further divided into two subsections) and proposed the genetic relationships among them. Differential cGISH signal patterns detected in section Alba allowed us to refine the genetic relationships among the closely related members of this section. CONCLUSIONS: We propose that GISH is an efficient tool to investigate the chromosome constitution and genetic relationships in mulberry. The results obtained here can be used to guide outbreeding of heterozygous perennial crops like mulberry.

Targeting of SLC25A22 boosts the immunotherapeutic response in KRAS-mutant colorectal cancer.

KRAS is an important tumor intrinsic factor driving immune suppression in colorectal cancer (CRC). In this study, we demonstrate that SLC25A22 underlies mutant KRAS-induced immune suppression in CRC. In immunocompetent male mice and humanized male mice models, SLC25A22 knockout inhibits KRAS-mutant CRC tumor growth with reduced myeloid derived suppressor cells (MDSC) but increased CD8+ T-cells, implying the reversion of mutant KRAS-driven immunosuppression. Mechanistically, we find that SLC25A22 plays a central role in promoting asparagine, which binds and activates SRC phosphorylation. Asparagine-mediated SRC promotes ERK/ETS2 signaling, which drives CXCL1 transcription. Secreted CXCL1 functions as a chemoattractant for MDSC via CXCR2, leading to an immunosuppressive microenvironment. Targeting SLC25A22 or asparagine impairs KRAS-induced MDSC infiltration in CRC. Finally, we demonstrate that the targeting of SLC25A22 in combination with anti-PD1 therapy synergizes to inhibit MDSC and activate CD8+ T cells to suppress KRAS-mutant CRC growth in vivo. We thus identify a metabolic pathway that drives immunosuppression in KRAS-mutant CRC.

Detecting liquid remnants of solid tumors treated with curative intent: Circulating tumor DNA as a biomarker of minimal residual disease (Review).

Circulating tumor DNA (ctDNA) has emerged as a promising biomarker of minimal residual disease (MRD) in solid tumors. There is increasing evidence to suggest that the detection of ctDNA following curative­intent treatments has high potential in anticipating future relapse in various solid tumors. Multiple liquid biopsy technical approaches and commercial platforms, including tumor­informed and tumor­agnostic ctDNA assays, have been developed for ctDNA­based MRD detection in solid tumors. Accurate ctDNA­based MRD analysis remains a critical technical challenge due to the very low concentration of ctDNA in peripheral blood samples, particularly in cancer patients following a curative­intent surgery or treatment. The present review summarizes the current key technical approaches that can be used to analyze ctDNA in the surveillance of MRD in solid tumors and provides a brief update on current commercial assays or platforms available for ctDNA­based MRD detection. The available evidence to date supporting ctDNA as a biomarker for detection of MRD in various types of solid tumors is also reviewed. In addition, technical and biological variables and considerations in pre­analytical and analytical steps associated with ctDNA­based MRD detection are discussed.

Up-regulation of oxytocin receptors on peripheral sensory neurons mediates analgesia in chemotherapy-induced neuropathic pain.

BACKGROUND AND PURPOSE: Chemotherapy-induced neuropathic pain (CINP) currently has limited effective treatment. Although the roles of oxytocin (OXT) and the oxytocin receptor (OXTR) in central analgesia have been well documented, the expression and function of OXTR in the peripheral nervous system remain unclear. Here, we evaluated the peripheral antinociceptive profiles of OXTR in CINP. EXPERIMENTAL APPROACH: Paclitaxel (PTX) was used to establish CINP. Quantitative real-time polymerase chain reaction (qRT-PCR), in situ hybridization, and immunohistochemistry were used to observe OXTR expression in dorsal root ganglia (DRG). The antinociceptive effects of OXT were assessed by hot-plate and von Frey tests. Whole-cell patch clamp was performed to record sodium currents, excitability of DRG neurons, and excitatory synapse transmission. KEY RESULTS: Expression of OXTR in DRG neurons was enhanced significantly after PTX treatment. Activation of OXTR exhibited antinociceptive effects, by decreasing the hyperexcitability of DRG neurons in PTX-treated mice. Additionally, OXTR activation up-regulated the phosphorylation of protein kinase C (pPKC) and, in turn, impaired voltage-gated sodium currents, particularly the voltage-gated sodium channel 1.7 (NaV 1.7) current, that plays an indispensable role in PTX-induced neuropathic pain. OXT suppressed excitatory transmission in the spinal dorsal horn as well as excitatory inputs from primary afferents in PTX-treated mice. CONCLUSION AND IMPLICATIONS: The OXTR in small-sized DRG neurons is up-regulated in CINP and its activation relieved CINP by inhibiting the neural excitability by impairment of NaV 1.7 currents via pPKC. Our results suggest that OXTR on peripheral sensory neurons is a potential therapeutic target to relieve CINP.

OncoPubMiner: a platform for mining oncology publications.

Updated and expert-quality knowledge bases are fundamental to biomedical research. A knowledge base established with human participation and subject to multiple inspections is needed to support clinical decision making, especially in the growing field of precision oncology. The number of original publications in this field has risen dramatically with the advances in technology and the evolution of in-depth research. Consequently, the issue of how to gather and mine these articles accurately and efficiently now requires close consideration. In this study, we present OncoPubMiner (https://oncopubminer.chosenmedinfo.com), a free and powerful system that combines text mining, data structure customisation, publication search with online reading and project-centred and team-based data collection to form a one-stop 'keyword in-knowledge out' oncology publication mining platform. The platform was constructed by integrating all open-access abstracts from PubMed and full-text articles from PubMed Central, and it is updated daily. OncoPubMiner makes obtaining precision oncology knowledge from scientific articles straightforward and will assist researchers in efficiently developing structured knowledge base systems and bring us closer to achieving precision oncology goals.

The cholesterol uptake regulator PCSK9 promotes and is a therapeutic target in APC/KRAS-mutant colorectal cancer.

Therapeutic targeting of KRAS-mutant colorectal cancer (CRC) is an unmet need. Here, we show that Proprotein Convertase Subtilisin/Kexin type 9 (PSCK9) promotes APC/KRAS-mutant CRC and is a therapeutic target. Using CRC patient cohorts, isogenic cell lines and transgenic mice, we identify that de novo cholesterol biosynthesis is induced in APC/KRAS mutant CRC, accompanied by increased geranylgeranyl diphosphate (GGPP)─a metabolite necessary for KRAS activation. PCSK9 is the top up-regulated cholesterol-related gene. PCSK9 depletion represses APC/KRAS-mutant CRC cell growth in vitro and in vivo, whereas PCSK9 overexpression induces oncogenesis. Mechanistically, PCSK9 reduces cholesterol uptake but induces cholesterol de novo biosynthesis and GGPP accumulation. GGPP is a pivotal metabolite downstream of PCSK9 by activating KRAS/MEK/ERK signaling. PCSK9 inhibitors suppress growth of APC/KRAS-mutant CRC cells, organoids and xenografts, especially in combination with simvastatin. PCSK9 overexpression predicts poor survival of APC/KRAS-mutant CRC patients. Together, cholesterol homeostasis regulator PCSK9 promotes APC/KRAS-mutant CRC via GGPP-KRAS/MEK/ERK axis and is a therapeutic target.

METTL3 Inhibits Antitumor Immunity by Targeting m6A-BHLHE41-CXCL1/CXCR2 Axis to Promote Colorectal Cancer.

BACKGROUND & AIMS: N6-Methyladenosine (m6A) is the most prevalent RNA modification and recognized as an important epitranscriptomic mechanism in colorectal cancer (CRC). We aimed to exploit whether and how tumor-intrinsic m6A modification driven by methyltransferase like 3 (METTL3) can dictate the immune landscape of CRC. METHODS: Mettl3 knockout mice, CD34+ humanized mice, and different syngeneic mice models were used. Immune cell composition and cytokine level were analyzed by flow cytometry and Cytokine 23-Plex immunoassay, respectively. M6A sequencing and RNA sequencing were performed to identify downstream targets and pathways of METTL3. Human CRC specimens (n = 176) were used to evaluate correlation between METTL3 expression and myeloid-derived suppressor cell (MDSC) infiltration. RESULTS: We demonstrated that silencing of METTL3 in CRC cells reduced MDSC accumulation to sustain activation and proliferation of CD4+ and CD8+ T cells, and eventually suppressed CRC in ApcMin/+Mettl3+/- mice, CD34+ humanized mice, and syngeneic mice models. Mechanistically, METTL3 activated the m6A-BHLHE41-CXCL1 axis by analysis of m6A sequencing, RNA sequencing, and cytokine arrays. METTL3 promoted BHLHE41 expression in an m6A-dependent manner, which subsequently induced CXCL1 transcription to enhance MDSC migration in vitro. However, the effect was negligible on BHLHE41 depletion, CXCL1 protein or CXCR2 inhibitor SB265610 administration, inferring that METTL3 promotes MDSC migration via BHLHE41-CXCL1/CXCR2. Consistently, depletion of MDSCs by anti-Gr1 antibody or SB265610 blocked the tumor-promoting effect of METTL3 in vivo. Importantly, targeting METTL3 by METTL3-single guide RNA or specific inhibitor potentiated the effect of anti-programmed cell death protein 1 (anti-PD1) treatment. CONCLUSIONS: Our study identifies METTL3 as a potential therapeutic target for CRC immunotherapy whose inhibition reverses immune suppression through the m6A-BHLHE41-CXCL1 axis. METTL3 inhibition plus anti-PD1 treatment shows promising antitumor efficacy against CRC.

Uncovering 1058 Novel Human Enteric DNA Viruses Through Deep Long-Read Third-Generation Sequencing and Their Clinical Impact.

BACKGROUND & AIMS: Lack of viral reference genomes poses a challenge to virome study. We investigated human gut virome and its clinical implication by ultra-deep metagenomic sequencing. METHODS: We extracted sufficient viral DNA from human feces for ultra-deep PacBio sequencing (>10 µg) and Illumina sequencing (>1 µg). Upon de novo assembly and 6 stages of strict filtering, viral genomes were generated and validated in 3 cohorts of 2819 published fecal metagenomes. Diagnostic performance of assembled viruses for colorectal cancer were tested in a training cohort and 2 independent validation cohorts. Virus mapping ratio, evolutionary history, and virus status (lytic or temperate) were also examined. RESULTS: The mean amount of extracted viral DNA increased by 14-fold compared with previous protocols. We obtained PacBio long reads and Illumina short reads with 290-fold higher depth than previous studies. We assembled and validated 1178 contigs as complete viral genomes, of which 1058 were newly identified. Thirteen viral genomes (398-839 kb) that are longer than the largest bacteriophage found in humans (393 kb) were discovered. Phylogenetic tree was constructed based on Hidden Markov Models alignment scores of 4 conserved viral proteins. Incorporating our assembled genomes into the National Center for Biotechnology Information database improved the mapping ratio of published metagenomes ≤18 times. Lytic viruses (75.9% ± 12.2% of total) were predominantly present in our sample. A biomarker panel of 14 novel viruses could discriminate patients with colorectal cancer from controls with an area under the receiver operating characteristics curve of 0.87 in the training cohort, which was validated with areas under the receiver operating characteristics curve of 0.85 and 0.73 in 2 independent cohorts. CONCLUSIONS: We uncovered 1058 novel human gut viruses. These findings can contribute to clinical diagnosis, current viral reference genome, and future virome investigation.

MarkovHC: Markov hierarchical clustering for the topological structure of high-dimensional single-cell omics data with transition pathway and critical point detection.

Clustering cells and depicting the lineage relationship among cell subpopulations are fundamental tasks in single-cell omics studies. However, existing analytical methods face challenges in stratifying cells, tracking cellular trajectories, and identifying critical points of cell transitions. To overcome these, we proposed a novel Markov hierarchical clustering algorithm (MarkovHC), a topological clustering method that leverages the metastability of exponentially perturbed Markov chains for systematically reconstructing the cellular landscape. Briefly, MarkovHC starts with local connectivity and density derived from the input and outputs a hierarchical structure for the data. We firstly benchmarked MarkovHC on five simulated datasets and ten public single-cell datasets with known labels. Then, we used MarkovHC to investigate the multi-level architectures and transition processes during human embryo preimplantation development and gastric cancer procession. MarkovHC found heterogeneous cell states and sub-cell types in lineage-specific progenitor cells and revealed the most possible transition paths and critical points in the cellular processes. These results demonstrated MarkovHC's effectiveness in facilitating the stratification of cells, identification of cell populations, and characterization of cellular trajectories and critical points.

Comprehensive review and evaluation of computational methods for identifying FLT3-internal tandem duplication in acute myeloid leukaemia.

Internal tandem duplication (ITD) of FMS-like tyrosine kinase 3 (FLT3-ITD) constitutes an independent indicator of poor prognosis in acute myeloid leukaemia (AML). AML with FLT3-ITD usually presents with poor treatment outcomes, high recurrence rate and short overall survival. Currently, polymerase chain reaction and capillary electrophoresis are widely adopted for the clinical detection of FLT3-ITD, whereas the length and mutation frequency of ITD are evaluated using fragment analysis. With the development of sequencing technology and the high incidence of FLT3-ITD mutations, a multitude of bioinformatics tools and pipelines have been developed to detect FLT3-ITD using next-generation sequencing data. However, systematic comparison and evaluation of the methods or software have not been performed. In this study, we provided a comprehensive review of the principles, functionality and limitations of the existing methods for detecting FLT3-ITD. We further compared the qualitative and quantitative detection capabilities of six representative tools using simulated and biological data. Our results will provide practical guidance for researchers and clinicians to select the appropriate FLT3-ITD detection tools and highlight the direction of future developments in this field. Availability: A Docker image with several programs pre-installed is available at https://github.com/niu-lab/docker-flt3-itd to facilitate the application of FLT3-ITD detection tools.

CircVAPA exerts oncogenic property in non-small cell lung cancer by the miR-876-5p/WNT5A axis.

BACKGROUND: Non-small cell lung cancer (NSCLC) is one of the most fatal malignant tumors. Emerging studies have clarified the crucial roles of circular RNAs (circRNAs) in the tumorigenesis of cancers. CircVAPA was demonstrated to function in some human cancers. The present study aimed to investigate the role of circVAPA in NSCLC. METHODS: A quantitative real-time polymerase chain reaction was used to measure the expression of genes. Actinomycin D and RNase R were employed to examine the stability of circVAPA. Cell-counting kit-8, 5-ethynyl-2'-deoxyuridine, Transwell and sphere formation assays, and well as western blot analysis, were conducted to examine the changes of NSCLC cells in response to circVAPA knockdown. A luciferase reporter assay was conducted for the molecular mechanism. RESULTS: Our findings demonstrated high expression of circVAPA in tissues and cell lines of NSCLC. Knockdown of circVAPA had a suppressive effect on cell proliferation, migration, invasion and stemness, and also inhibited tumor growth in vivo. Mechanistically, circVAPA acted as a competing endogenous RNA to up-regulate WNT5A by sponging miR-876-5p. Moreover, circVAPA activated Wnt/ß-catenin signaling by up-regulation of WNT5A. Rescue assays showed that silencing of miR-876-5p or overexpression of WNT5A reversed the circVAPA knockdown-mediated inhibition on cellular processes in NSCLC. CONCLUSIONS: CircVAPA promotes aggressive phenotypes of NSCLC cells by the miR-876-5p/WNT5A axis activating Wnt/ß-catenin signaling.

Anti-oncogene PTPN13 inactivation by hepatitis B virus X protein counteracts IGF2BP1 to promote hepatocellular carcinoma progression.

Hepatitis B x protein (HBx) affects cellular protein expression and participates in the tumorigenesis and progression of hepatitis B virus (HBV)-related hepatocellular carcinoma (HCC). Metabolic reprogramming contributed to the HCC development, but its role in HBV-related HCC remains largely unclear. Tyrosine-protein phosphatase nonreceptor type 13 (PTPN13) is a significant regulator in tumor development, however, its specific role in hepatocarcinogenesis remains to be explored. Here, we found that decreased PTPN13 expression was associated with HBV/HBx. Patients with low PTPN13 expression showed a poor prognosis. Functional assays revealed that PTPN13 inhibited proliferation and tumorigenesis in vitro and in vivo. Further mechanistic studies indicated that HBx inhibited PTPN13 expression by upregulating the expression of DNMT3A and interacting with DNMT3A. Furthermore, we found that DNMT3A bound to the PTPN13 promoter (-343 to -313 bp) in an epigenetically controlled manner associated with elevated DNA methylation and then inhibited PTPN13 transcription. In addition, we identified IGF2BP1 as a novel PTPN13-interacting gene and demonstrated that PTPN13 influences c-Myc expression by directly and competitively binding to IGF2BP1 to decrease the intracellular concentration of functional IGF2BP1. Overexpressing PTPN13 promoted c-Myc mRNA degradation independent of the protein tyrosine phosphatase (PTP) activity of PTPN13. Importantly, we discovered that the PTPN13-IGF2BP1-c-Myc axis was important for cancer cell growth through promoting metabolic reprogramming. We verified the significant negative correlations between PTPN13 expression and c-Myc, PSPH, and SLC7A1 expression in clinical HCC tissue samples. In summary, our findings demonstrate that PTPN13 is a novel regulator of HBV-related hepatocarcinogenesis and may play an important role in HCC. PTPN13 may serve as a prognostic marker and therapeutic target in HBV-related HCC patients.

PRL-3 facilitates Hepatocellular Carcinoma progression by co-amplifying with and activating FAK.

Background: In addition to protein tyrosine kinases, accumulating evidence has shown that protein tyrosine phosphatases (PTPs) are suitable therapeutic targets in cancer. PRL-3 is a PTP member that has been well studied in many malignant tumours. The goal of the present study was to elucidate the role of PRL-3 in hepatocellular carcinoma (HCC), which remains largely unknown. Methods: Bioinformatic and immunohistochemical analyses were performed to analyse PRL-3 expression in HCC tissue samples and determine its clinical relevance. PRL-3 gene copy number variations were evaluated by bioinformatic analysis and quantitative-genomic polymerase chain reaction. The biological functions of PRL-3 were investigated in vivo and vitro. Gene microarray assays, RT-qPCR, western blotting and luciferase experiments were performed to identify the downstream effectors of PRL-3 that mediate its functions in HCC. Results: PRL-3 expression was upregulated in HCC samples from public databases and in cohort samples from our centre. High PRL-3 expression was associated with poor prognosis. Copy number gains and amplification of chromosome 8q24.3 in HCC were determined to be positively correlated with the PRL-3 overexpression. PRL-3 overexpression promoted HCC cell proliferation, migration and adhesion, while its loss had the opposite effects. Further study showed that focal adhesion kinase (FAK) was co-amplified and co-expressed with PRL-3 in HCC. Interestingly, PRL-3 also promoted the phosphorylation of FAK, which subsequently mediated the oncogenic functions of PRL-3 in HCC cells. Moreover, TGFB1 was identified as a downstream molecule of PRL-3. TGF-ß signalling was shown to mediate the PRL-3-induced activation of FAK. Furthermore, the p38 and PI3K/AKT pathways were observed to mediate the PRL-3-induced expression of TGFB1 and the subsequent activation of FAK, while the activation of FAK in turn stimulated activation of the p38 and PI3K/AKT pathways, forming a PRL-3-triggered AKT/p38/TGFB1/FAK positive feedback loop. Conclusion: Collectively, our findings indicate that the PTP PRL-3 plays a crucial role in the progression of HCC and provides an example of how co-amplified genes work together in HCC.

In Colorectal Cancer Cells With Mutant KRAS, SLC25A22-Mediated Glutaminolysis Reduces DNA Demethylation to Increase WNT Signaling, Stemness, and Drug Resistance.

BACKGROUND & AIMS: Mutant KRAS promotes glutaminolysis, a process that uses steps from the tricarboxylic cycle to convert glutamine to α-ketoglutarate and other molecules via glutaminase and SLC25A22. This results in inhibition of demethylases and epigenetic alterations in cells that increase proliferation and stem cell features. We investigated whether mutant KRAS-mediated glutaminolysis affects the epigenomes and activities of colorectal cancer (CRC) cells. METHODS: We created ApcminKrasG12D mice with intestine-specific knockout of SLC25A22 (ApcminKrasG12DSLC25A22fl/fl mice). Intestine tissues were collected and analyzed by histology, immunohistochemistry, and DNA methylation assays; organoids were derived and studied for stem cell features, along with organoids derived from 2 human colorectal tumor specimens. Colon epithelial cells (1CT) and CRC cells (DLD1, DKS8, HKE3, and HCT116) that expressed mutant KRAS, with or without knockdown of SLC25A22 or other proteins, were deprived of glutamine or glucose and assayed for proliferation, colony formation, glucose or glutamine consumption, and apoptosis; gene expression patterns were analyzed by RNA sequencing, proteins by immunoblots, and metabolites by liquid chromatography-mass spectrometry, with [U-13C5]-glutamine as a tracer. Cells and organoids with knocked down, knocked out, or overexpressed proteins were analyzed for DNA methylation at CpG sites using arrays. We performed immunohistochemical analyses of colorectal tumor samples from 130 patients in Hong Kong (57 with KRAS mutations) and Kaplan-Meier analyses of survival. We analyzed gene expression levels of colorectal tumor samples in The Cancer Genome Atlas. RESULTS: CRC cells that express activated KRAS required glutamine for survival, and rapidly incorporated it into the tricarboxylic cycle (glutaminolysis); this process required SLC25A22. Cells incubated with succinate and non-essential amino acids could proliferate under glutamine-free conditions. Mutant KRAS cells maintained a low ratio of α-ketoglutarate to succinate, resulting in reduced 5-hydroxymethylcytosine-a marker of DNA demethylation, and hypermethylation at CpG sites. Many of the hypermethylated genes were in the WNT signaling pathway and at the protocadherin gene cluster on chromosome 5q31. CRC cells without mutant KRAS, or with mutant KRAS and knockout of SLC25A22, expressed protocadherin genes (PCDHAC2, PCDHB7, PCDHB15, PCDHGA1, and PCDHGA6)-DNA was not methylated at these loci. Expression of the protocadherin genes reduced WNT signaling to ß-catenin and expression of the stem cell marker LGR5. ApcminKrasG12DSLC25A22fl/fl mice developed fewer colon tumors than ApcminKrasG12D mice (P < .01). Organoids from ApcminKrasG12DSLC25A22fl/fl mice had reduced expression of LGR5 and other markers of stemness compared with organoids derived from ApcminKrasG12D mice. Knockdown of SLC25A22 in human colorectal tumor organoids reduced clonogenicity. Knockdown of lysine demethylases, or succinate supplementation, restored expression of LGR5 to SLC25A22-knockout CRC cells. Knockout of SLC25A22 in CRC cells that express mutant KRAS increased their sensitivity to 5-fluorouacil. Level of SLC25A22 correlated with levels of LGR5, nuclear ß-catenin, and a stem cell-associated gene expression pattern in human colorectal tumors with mutations in KRAS and reduced survival times of patients. CONCLUSIONS: In CRC cells that express activated KRAS, SLC25A22 promotes accumulation of succinate, resulting in increased DNA methylation, activation of WNT signaling to ß-catenin, increased expression of LGR5, proliferation, stem cell features, and resistance to 5-fluorouacil. Strategies to disrupt this pathway might be developed for treatment of CRC.