TSGene 2.0: an updated literature-based knowledgebase for tumor suppressor genes

Comprehensive characterization of somatic mutations associated with chimeric RNAs in human cancers

Chimeric RNAs, which can arise from gene recombination at the DNA level or non-canonical splicing events at the RNA level, have been identified as important roles in human tumors. Dysregulated gene expression caused by somatic mutations and altered splicing patterns of oncogenes or tumor suppressor genes can contribute to the development of tumors. Therefore, investigating the formation mechanism of chimeric RNAs via somatic mutations is critical for understanding tumor pathogenesis. This project is the first to propose studying the association between somatic single nucleotide variants and chimeric RNAs, identifying around 2900 somatic SNVs affecting chimeric RNAs in pan-cancer level. The somatic SNVs on chimeric RNAs were commonly observed in various types of tumor tissues, providing a valuable resource for future study. Additionally, these SNVs show distinct tumor specificity, and those with high frequency had a significant impact on the survival time of patients with tumors. Further research revealed that somatic SNVs associated with chimeric RNA (chiR-SNVs) were typically found within 10 nt of the junction site of chimeric RNAs and had a particularly significant effect on chimeric RNAs from different chromosomes. The enrichment analysis revealed that chiR-SNVs were significantly overrepresented in oncogenes and genes related to RNA binding proteins involved in RNA splicing, which could imply that chiR-SNVs may disrupt the process of RNA splicing and induce the occurrence of chimeric RNAs. This study sheds light on the potential molecular interaction mechanism between somatic SNVs and chimeric RNAs, which opens up a new avenue for researching disease pathway and tumorigenesis development.

Identifying and evaluating a disulfidptosis-related gene signature to predict prognosis in colorectal adenocarcinoma patients

Disulfidptosis, a regulated form of cell death, has been recently reported in cancers characterized by high SLC7A11 expression, including invasive breast carcinoma, lung adenocarcinoma, and hepatocellular carcinoma. However, its role in colon adenocarcinoma (COAD) has been infrequently discussed. In this study, we developed and validated a prognostic model based on 20 disulfidptosis-related genes (DRGs) using LASSO and Cox regression analyses. The robustness and practicality of this model were assessed via a nomogram. Subsequent correlation and enrichment analysis revealed a relationship between the risk score, several critical cancer-related biological processes, immune cell infiltration, and the expression of oncogenes and cell senescence-related genes. POU4F1, a significant component of our model, might function as an oncogene due to its upregulation in COAD tumors and its positive correlation with oncogene expression. In vitro assays demonstrated that POU4F1 knockdown noticeably decreased cell proliferation and migration but increased cell senescence in COAD cells. We further investigated the regulatory role of the DRG in disulfidptosis by culturing cells in a glucose-deprived medium. In summary, our research revealed and confirmed a DRG-based risk prediction model for COAD patients and verified the role of POU4F1 in promoting cell proliferation, migration, and disulfidptosis.

The Neo-Open Reading Frame Peptides That Comprise the Tumor Framome Are a Rich Source of Neoantigens for Cancer Immunotherapy

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate the whole-genome and long-read transcript sequencing of cancers to identify the collection of neo-open reading frame peptides (NOP) expressed in tumors. We termed this collection of NOPs the tumor framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe a class of hidden NOPs that derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of noncoding regions of the genome downstream of a rearrangement breakpoint, i.e., where no gene annotation or evidence for transcription exists. The entire collection of NOPs represents a vast number of possible neoantigens particularly in tumors with many structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T cells specific for hidden NOPs in peripheral blood from a patient with lung cancer. This work highlights NOPs as a major source of possible neoantigens for personalized cancer immunotherapy and provides a rationale for analyzing the complete cancer genome and transcriptome as a basis for the detection of NOPs.

Oncogenes and tumor suppressor genes: functions and roles in cancers

Cancer, being the most formidable ailment, has had a profound impact on the human health. The disease is primarily associated with genetic mutations that impact oncogenes and tumor suppressor genes (TSGs). Recently, growing evidence have shown that X-linked TSGs have specific role in cancer progression and metastasis as well. Interestingly, our genome harbors around substantial portion of genes that function as tumor suppressors, and the X chromosome alone harbors a considerable number of TSGs. The scenario becomes even more compelling as X-linked TSGs are adaptive to key epigenetic processes such as X chromosome inactivation. Therefore, delineating the new paradigm related to X-linked TSGs, for instance, their crosstalk with autosome and involvement in cancer initiation, progression, and metastasis becomes utmost importance. Considering this, herein, we present a comprehensive discussion of X-linked TSG dysregulation in various cancers as a consequence of genetic variations and epigenetic alterations. In addition, the dynamic role of X-linked TSGs in sex chromosome-autosome crosstalk in cancer genome remodeling is being explored thoroughly. Besides, the functional roles of ncRNAs, role of X-linked TSG in immunomodulation and in gender-based cancer disparities has also been highlighted. Overall, the focal idea of the present article is to recapitulate the findings on X-linked TSG regulation in the cancer landscape and to redefine their role toward improving cancer treatment strategies.

Identifying and assessing a prognostic model based on disulfidptosis-related genes: implications for immune microenvironment and tumor biology in lung adenocarcinoma

Introduction

Lung cancer, with the highest global mortality rate among cancers, presents a grim prognosis, often diagnosed at an advanced stage in nearly 70% of cases. Recent research has unveiled a novel mechanism of cell death termed disulfidptosis, which is facilitated by glucose scarcity and the protein SLC7A11.

Methods

Utilizing the least absolute shrinkage and selection operator (LASSO) regression analysis combined with Cox regression analysis, we constructed a prognostic model focusing on disulfidptosis-related genes. Nomograms, correlation analyses, and enrichment analyses were employed to assess the significance of this model. Among the genes incorporated into the model, CHRNA5 was selected for further investigation regarding its role in LUAD cells. Biological functions of CHRNA5 were assessed using EdU, transwell, and CCK-8 assays.

Results

The efficacy of the model was validated through internal testing and an external validation set, with further evaluation of its robustness and clinical applicability using a nomogram. Subsequent correlation analyses revealed associations between the risk score and infiltration of various cancer types, as well as oncogene expression. Enrichment analysis also identified associations between the risk score and pivotal biological processes and KEGG pathways. Our findings underscore the significant impact of CHRNA5 on LUAD cell proliferation, migration, and disulfidptosis.

Conclusion

This study successfully developed and validated a robust prognostic model centered on disulfidptosis-related genes, providing a foundation for predicting prognosis in LUAD patients.

Oncogenic composite mutations can be predicted by co-mutations and their chromosomal location

Genetic heterogeneity in tumors can show a remarkable selectivity when two or more independent genetic events occur in the same gene. This phenomenon, called composite mutation, points toward a selective pressure, which frequently causes therapy resistance to mutation-specific drugs. Since composite mutations have been described to occur in sub-clonal populations, they are not always captured through biopsy sampling. Here, we provide a proof of concept to predict composite mutations to anticipate which patients might be at risk for sub-clonally driven therapy resistance. We found that composite mutations occur in 5% of cancer patients, mostly affecting the PIK3CA, EGFR, BRAF, and KRAS genes, which are common precision medicine targets. Furthermore, we found a strong and significant relationship between the frequencies of composite mutations with commonly co-occurring mutations in a non-composite context. We also found that co-mutations are significantly enriched on the same chromosome. These observations were independently confirmed using cell line data. Finally, we show the feasibility of predicting compositive mutations based on their co-mutations (AUC 0.62, 0.81, 0.82, and 0.91 for EGFR, PIK3CA, KRAS, and BRAF, respectively). This prediction model could help to stratify patients who are at risk of developing therapy resistance-causing mutations.

Network-based elucidation of colon cancer drug resistance mechanisms by phosphoproteomic time-series analysis

Aberrant signaling pathway activity is a hallmark of tumorigenesis and progression, which has guided targeted inhibitor design for over 30 years. Yet, adaptive resistance mechanisms, induced by rapid, context-specific signaling network rewiring, continue to challenge therapeutic efficacy. Leveraging progress in proteomic technologies and network-based methodologies, we introduce Virtual Enrichment-based Signaling Protein-activity Analysis (VESPA)-an algorithm designed to elucidate mechanisms of cell response and adaptation to drug perturbations-and use it to analyze 7-point phosphoproteomic time series from colorectal cancer cells treated with clinically-relevant inhibitors and control media. Interrogating tumor-specific enzyme/substrate interactions accurately infers kinase and phosphatase activity, based on their substrate phosphorylation state, effectively accounting for signal crosstalk and sparse phosphoproteome coverage. The analysis elucidates time-dependent signaling pathway response to each drug perturbation and, more importantly, cell adaptive response and rewiring, experimentally confirmed by CRISPR knock-out assays, suggesting broad applicability to cancer and other diseases.

Integrative analysis in head and neck cancer reveals distinct role of miRNome and methylome as tumour epigenetic drivers

Head and neck cancer is the sixth most common malignancy worldwide, with the relatively low 5-year survival rate, mainly because it is diagnosed at a late stage. Infection with HPV is a well known aetiology, which affects the nature of these cancers and patients' survival. Besides, it is considered that the main driving force for this type of cancer could be epigenetics. In this study we aimed to find potential epigenetic biomarkers, by integrating miRNome, methylome, and transcriptome analyses. From the fresh head and neck cancer tissue samples, we chose a group for miRNome, methylome and transcriptome profiling, in comparison to adequate control samples. Bioinformatics analyses are performed in R v4.2.2. Count normalisation and group differential expression for mRNA and the previously obtained miRNA count data was performed with DESeq2 v1.36. Gene set enrichment analysis was performed and visualised using gProfiler2 v0.2.1 Identification of miRNA targets was performed by querying in miRTarBase using multiMiR v1.18.0. Annotation of CpG sites merging into islands was obtained from RnBeads.hg19 v1.28.0. package. For the integrative analysis we performed kmeans clustering using stats v4.2.2 package, using 8-12 clusters and nstart 100. We found that transcriptome analysis divides samples into cancers and controls clusters, with no relation to HPV status or cancer anatomical location. Differentially expressed genes (n = 2781) were predominantly associated with signalling pathways of tumour progression. We identified a cluster of genes under the control of the transcription factor E2F that are significantly underexpressed in cancer tissue, as well as T cell immunity genes and genes related to regulation of transcription. Among overexpressed genes in tumours we found those that belong to cell cycle regulation and vasculature. A small number of genes were found significantly differentially expressed in HPV-positive versus HPV-negative tumours (for example NEFH, ZFR2, TAF7L, ZNF541, and TYMS). In this comprehensive study on an overlapping set of samples where the integration of miRNome, methylome and transcriptome analysis were performed for head and neck cancer, we demonstrated that the majority of genes were associated exclusively with miRNome or methylome and, to a lesser extent, under the control of both epigenetic mechanisms.

CancerMHL: the database of integrating key DNA methylation, histone modifications and lncRNAs in cancer

The discovery of key epigenetic modifications in cancer is of great significance for the study of disease biomarkers. Through the mining of epigenetic modification data relevant to cancer, some researches on epigenetic modifications are accumulating. In order to make it easier to integrate the effects of key epigenetic modifications on the related cancers, we established CancerMHL (http://www.positionprediction.cn/), which provide key DNA methylation, histone modifications and lncRNAs as well as the effect of these key epigenetic modifications on gene expression in several cancers. To facilitate data retrieval, CancerMHL offers flexible query options and filters, allowing users to access specific key epigenetic modifications according to their own needs. In addition, based on the epigenetic modification data, three online prediction tools had been offered in CancerMHL for users. CancerMHL will be a useful resource platform for further exploring novel and potential biomarkers and therapeutic targets in cancer. Database URL: http://www.positionprediction.cn/.

Polyploidy Promotes Hypertranscription, Apoptosis Resistance, and Ciliogenesis in Cancer Cells and Mesenchymal Stem Cells of Various Origins: Comparative Transcriptome In Silico Study

Mesenchymal stem cells (MSC) attract an increasing amount of attention due to their unique therapeutic properties. Yet, MSC can undergo undesirable genetic and epigenetic changes during their propagation in vitro. In this study, we investigated whether polyploidy can compromise MSC oncological safety and therapeutic properties. For this purpose, we compared the impact of polyploidy on the transcriptome of cancer cells and MSC of various origins (bone marrow, placenta, and heart). First, we identified genes that are consistently ploidy-induced or ploidy-repressed through all comparisons. Then, we selected the master regulators using the protein interaction enrichment analysis (PIEA). The obtained ploidy-related gene signatures were verified using the data gained from polyploid and diploid populations of early cardiomyocytes (CARD) originating from iPSC. The multistep bioinformatic analysis applied to the cancer cells, MSC, and CARD indicated that polyploidy plays a pivotal role in driving the cell into hypertranscription. It was evident from the upregulation of gene modules implicated in housekeeping functions, stemness, unicellularity, DNA repair, and chromatin opening by means of histone acetylation operating via DNA damage associated with the NUA4/TIP60 complex. These features were complemented by the activation of the pathways implicated in centrosome maintenance and ciliogenesis and by the impairment of the pathways related to apoptosis, the circadian clock, and immunity. Overall, our findings suggest that, although polyploidy does not induce oncologic transformation of MSC, it might compromise their therapeutic properties because of global epigenetic changes and alterations in fundamental biological processes. The obtained results can contribute to the development and implementation of approaches enhancing the therapeutic properties of MSC by removing polyploid cells from the cell population.

Spatial transcriptomics reveals a low extent of transcriptionally active hepatitis B virus integration in patients with HBsAg loss

OBJECTIVE

Hepatitis B virus (HBV) can integrate into the chromosomes of infected hepatocytes, contributing to the production of hepatitis B surface antigen (HBsAg) and to hepatocarcinogenesis. In this study, we aimed to explore whether transcriptionally active HBV integration events spread throughout the liver tissue in different phases of chronic HBV infection, especially in patients with HBsAg loss.

DESIGN

We constructed high-resolution spatial transcriptomes of liver biopsies containing 13 059 tissue spots from 18 patients with chronic HBV infection to analyse the occurrence and relative distribution of transcriptionally active viral integration events. Immunohistochemistry was performed to evaluate the expression of HBsAg and HBV core antigen. Intrahepatic covalently closed circular DNA (cccDNA) levels were quantified by real-time qPCR.

RESULTS

Spatial transcriptome sequencing identified the presence of 13 154 virus-host chimeric reads in 7.86% (1026 of 13 059) of liver tissue spots in all patients, including three patients with HBsAg loss. These HBV integration sites were randomly distributed on chromosomes and can localise in host genes involved in hepatocarcinogenesis, such as ALB, CLU and APOB. Patients who were receiving or had received antiviral treatment had a significantly lower percentage of viral integration-containing spots and significantly fewer chimeric reads than treatment-naïve patients. Intrahepatic cccDNA levels correlated well with viral integration events.

CONCLUSION

Transcriptionally active HBV integration occurred in chronically HBV-infected patients at different phases, including in patients with HBsAg loss. Antiviral treatment was associated with a decreased number and extent of transcriptionally active viral integrations, implying that early treatment intervention may further reduce the number of viral integration events.

Histone 3 lysine 9 acetylation-specific reprogramming regulates esophageal squamous cell carcinoma progression and metastasis

Dysregulation of histone acetylation is widely implicated in tumorigenesis, yet its specific roles in the progression and metastasis of esophageal squamous cell carcinoma (ESCC) remain unclear. Here, we profiled the genome-wide landscapes of H3K9ac for paired adjacent normal (Nor), primary ESCC (EC) and metastatic lymph node (LNC) esophageal tissues from three ESCC patients. Compared to H3K27ac, we identified a distinct epigenetic reprogramming specific to H3K9ac in EC and LNC samples relative to Nor samples. This H3K9ac-related reprogramming contributed to the transcriptomic aberration of targeting genes, which were functionally associated with tumorigenesis and metastasis. Notably, genes with gained H3K9ac signals in both primary and metastatic lymph node samples (common-gained gene) were significantly enriched in oncogenes. Single-cell RNA-seq analysis further revealed that the corresponding top 15 common-gained genes preferred to be enriched in mesenchymal cells with high metastatic potential. Additionally, in vitro experiment demonstrated that the removal of H3K9ac from the common-gained gene MSI1 significantly downregulated its transcription, resulting in deficiencies in ESCC cell proliferation and migration. Together, our findings revealed the distinct characteristics of H3K9ac in esophageal squamous cell carcinogenesis and metastasis, and highlighted the potential therapeutic avenue for intervening ESCC through epigenetic modulation via H3K9ac.

Transcriptional and Translational Regulation of Differentially Expressed Genes in Yucatan Miniswine Brain Tissues following Traumatic Brain Injury

Traumatic brain injury (TBI) is a leading cause of morbidity, disability, and mortality worldwide. Motor and cognitive deficits and emotional disturbances are long-term consequences of TBI. A lack of effective treatment for TBI-induced neural damage, functional impairments, and cognitive deficits makes it challenging in the recovery following TBI. One of the reasons may be the lack of knowledge underlying the complex pathophysiology of TBI and the regulatory factors involved in the cellular and molecular mechanisms of inflammation, neural regeneration, and injury repair. These mechanisms involve a change in the expression of various proteins encoded by genes whose expression is regulated by transcription factors (TFs) at the transcriptional level and microRNA (miRs) at the mRNA level. In this pilot study, we performed the RNA sequencing of injured tissues and non-injured tissues from the brain of Yucatan miniswine and analyzed the sequencing data for differentially expressed genes (DEGs) and the TFs and miRs regulating the expression of DEGs using in-silico analysis. We also compared the effect of the electromagnetic field (EMF) applied to the injured miniswine on the expression profile of various DEGs. The results of this pilot study revealed a few DEGs that were significantly upregulated in the injured brain tissue and the EMF stimulation showed effect on their expression profile.

CSSLdb: Discovery of cancer-specific synthetic lethal interactions based on machine learning and statistic inference

Synthetic lethality (SL) occurs when the inactivation of two genes results in cell death while the inactivation of either gene alone is non-lethal. SL-based therapy has become a promising anti-cancer treatment option with the increasing researches and applications in clinical practice, while the specific therapeutic opportunities for various cancers have not yet been comprehensively investigated. Herein, we described a computational approach based on machine learning and statistical inference to discover the cancer-specific synthetic lethal interactions. First, Random Forest and One-Class SVM were used to perform cancer unbiased prediction of synthetic lethality. Then, two strategies, including mutual exclusivity and differential expression, were used to screen cancer-specific synthetic lethal interactions, resulting in 14,582 SL gene pairs in 33 cancer types. Finally, we developed a freely available database of CSSLdb (Cancer Specific Synthetic Lethality Database, http://www.tmliang.cn/CSSL/) to present cancer-specific synthetic lethal genetic interactions, which would enrich the relevant research and contribute to underlying therapy strategies based on synthetic lethality.

Constructing and Evaluating a Mitophagy-Related Gene Prognostic Model: Implications for Immune Landscape and Tumor Biology in Lung Adenocarcinoma

Mitophagy, a conserved cellular mechanism, is crucial for cellular homeostasis through the selective clearance of impaired mitochondria. Its emerging role in cancer development has sparked interest, particularly in lung adenocarcinoma (LUAD). Our study aimed to construct a risk model based on mitophagy-related genes (MRGs) to predict survival outcomes, immune response, and chemotherapy sensitivity in LUAD patients. We mined the GeneCards database to identify MRGs and applied LASSO/Cox regression to formulate a prognostic model. Validation was performed using two independent Gene Expression Omnibus (GEO) cohorts. Patients were divided into high- and low-risk categories according to the median risk score. The high-risk group demonstrated significantly reduced survival. Multivariate Cox analysis confirmed the risk score as an independent predictor of prognosis, and a corresponding nomogram was developed to facilitate clinical assessments. Intriguingly, the risk score correlated with immune infiltration levels, oncogenic expression profiles, and sensitivity to anticancer agents. Enrichment analyses linked the risk score with key oncological pathways and biological processes. Within the model, MTERF3 emerged as a critical regulator of lung cancer progression. Functional studies indicated that the MTERF3 knockdown suppressed the lung cancer cell proliferation and migration, enhanced mitophagy, and increased the mitochondrial superoxide production. Our novel prognostic model, grounded in MRGs, promises to refine therapeutic strategies and prognostication in lung cancer management.

Pan-cancer analyses identify MIR210HG overexpression, epigenetic regulation and oncogenic role in human tumors and its interaction with the tumor microenvironment

BACKGROUND

Molecular entities showing dysregulation in multiple cancers may hold great biomarker or therapeutic potential. There is accumulating evidence that highlights the dysregulation of a long non-coding RNA, MIR210HG, in various cancers and its oncogenic role. However, a comprehensive analysis of MIR210HG expression pattern, molecular mechanisms, diagnostic or prognostic significance or evaluation of its interaction with tumor microenvironment across various cancers remains unstudied.

METHODS

A systematic pan-cancer analysis was done using multiple public databases and bioinformatic tools to study the molecular role and clinical significance of MIR210HG. We have analyzed expression patterns, genome alteration, transcriptional and epigenetic regulation, correlation with patient survival, immune infiltrates, co-expressed genes, interacting proteins, and pathways associated with MIR210HG.

RESULTS

The Pan cancer expression analysis of MIR210HG through various tumor datasets demonstrated that MIR210HG is significantly upregulated in most cancers and increased with the tumor stage in a subset of them. Furthermore, prognostic analysis revealed high MIR210HG expression is associated with poor overall and disease-free survival in specific cancer types. Genetic alteration analysis showed minimal alterations in the MIR210HG locus, indicating that overexpression in cancers is not due to gene amplification. The exploration of SNPs on MIR210HG suggested possible structural changes that may affect its interactions with the miRNAs. The correlation of MIR210HG with promoter methylation was found to be significantly negative in nature in majority of cancers depicting the possible epigenetic regulation of expression of MIR210HG. Additionally, MIR210HG showed negative correlations with immune cells and thus may have strong impact on the tumor microenvironment. Functional analysis indicates its association with hypoxia, angiogenesis, metastasis, and DNA damage repair processes. MIR210HG was found to interact with several proteins and potentially regulate chromatin modifications and transcriptional regulation.

CONCLUSIONS

A first pan-can cancer analysis of MIR210HG highlights its transcriptional and epigenetic deregulation and oncogenic role in the majority of cancers, its correlation with tumor microenvironment factors such as hypoxia and immune infiltration, and its potential as a prognostic biomarker and therapeutic target in several cancers.

Bioinformatics Prediction for Network-Based Integrative Multi-Omics Expression Data Analysis in Hirschsprung Disease

Hirschsprung's disease (HSCR) is a rare developmental disorder in which enteric ganglia are missing along a portion of the intestine. HSCR has a complex inheritance, with RET as the major disease-causing gene. However, the pathogenesis of HSCR is still not completely understood. Therefore, we applied a computational approach based on multi-omics network characterization and clustering analysis for HSCR-related gene/miRNA identification and biomarker discovery. Protein-protein interaction (PPI) and miRNA-target interaction (MTI) networks were analyzed by DPClusO and BiClusO, respectively, and finally, the biomarker potential of miRNAs was computationally screened by miRNA-BD. In this study, a total of 55 significant gene-disease modules were identified, allowing us to propose 178 new HSCR candidate genes and two biological pathways. Moreover, we identified 12 key miRNAs with biomarker potential among 137 predicted HSCR-associated miRNAs. Functional analysis of new candidates showed that enrichment terms related to gene ontology (GO) and pathways were associated with HSCR. In conclusion, this approach has allowed us to decipher new clues of the etiopathogenesis of HSCR, although molecular experiments are further needed for clinical validations.

Development of an invasion score based on metastasis-related pathway activity profiles for identifying invasive molecular subtypes of lung adenocarcinoma

The invasive capacity of lung adenocarcinoma (LUAD) is an important factor influencing patients' metastatic status and survival outcomes. However, there is still a lack of suitable biomarkers to evaluate tumor invasiveness. LUAD molecular subtypes were identified by unsupervised consistent clustering of LUAD. The differences in prognosis, tumor microenvironment (TME), and mutation were assessed among different subtypes. After that, the invasion-related gene score (IRGS) was constructed by genetic differential analysis, WGCNA analysis, and LASSO analysis, then we evaluated the relationship between IRGS and invasive characteristics, TME, and prognosis. The predictive ability of the IRGS was verified by in vitro experiments. Next, the "oncoPredict" R package and CMap were used to assess the potential value of IRGS in drug therapy. The results showed that LUAD was clustered into two molecular subtypes. And the C1 subtype exhibited a worse prognosis, higher stemness enrichment activity, less immune infiltration, and higher mutation frequency. Subsequently, IRGS developed based on molecular subtypes demonstrated a strong association with malignant characteristics such as invasive features, higher stemness scores, less immune infiltration, and worse survival. In vitro experiments showed that the higher IRGS LUAD cell had a stronger invasive capacity than the lower IRGS LUAD cell. Predictive analysis based on the "oncoPredict" R package showed that the high IRGS group was more sensitive to docetaxel, erlotinib, paclitaxel, and gefitinib. Among them, in vitro experiments verified the greater killing effect of paclitaxel on high IRGS cell lines. In addition, CMap showed that purvalanol-a, angiogenesis-inhibitor, and masitinib have potential therapeutic effects in the high IRGS group. In summary we identified and analyzed the molecular subtypes associated with the invasiveness of LUAD and developed IRGS that can efficiently predict the prognosis and invasive ability of the tumor. IRGS may be able to facilitate the precision treatment of LUAD to some extent.

EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

Shikonin reactivates TSGs GADD45B and PPP3CC to block NSCLC cell proliferation and migration through JNK/P38/MAPK signaling pathways

BACKGROUND

Shikonin, a natural naphthoquinone compound extracted from the Chinese traditional herbal medicine "Lithospermum erythrorhizon", possesses antitumor activity against various cancer types. Tumor-suppressor genes (TSGs) negatively regulate cell growth, proliferation, and differentiation, thereby inhibiting tumor formation. However, the molecular mechanism of action of shikonin on TSGs in non-small-cell lung cancer (NSCLC) remains unclear.

METHODS

The inhibitory effect of shikonin on the proliferation and migration abilities of lung cancer cells were measured by Cell Counting Kit 8 (CCK8) and wound healing assays. The alteration of genes by shikonin treatment was detected by mRNA high-throughput sequencing and further confirmed by qPCR and western blotting experiments. The dominant functions of the upregulated genes were analyzed by GO and KEGG profiling.

RESULTS

Shikonin inhibited the proliferation and migration of A549 and H1299 NSCLC cells in a dose-dependent manner. mRNA high-throughput sequencing revealed a total of 1794 upregulated genes in shikonin-treated NSCLC cells. Moreover, bioinformatic analysis of GO and KEGG profiling revealed that the up-regulated genes were mostly involved in the JNK/P38/MAPK signaling pathway, among which the expression of GADD45B and PPP3CC was significantly enhanced. Finally, we confirmed that GADD45B and PPP3CC were indeed upregulated in JNK/P38/MAPK pathway.

CONCLUSIONS

Taken together, these results suggested that shikonin might affect the expression of GADD45B and PPP3CC through the JNK/P38/MAPK pathway, therefore exerting an inhibitory effect on the proliferation and migration of cancer cells. To our knowledge, this is the first study reporting the role of shikonin in upregulating TSGs to activate the JNK/P38/MAPK signaling pathways in NSCLC.

AEG-1 as a Novel Therapeutic Target in Colon Cancer: A Study from Silencing AEG-1 in BALB/c Mice to Large Data Analysis

BACKGROUND

Astrocyte elevated gene-1 (AEG-1) is overexpressed in various malignancies. Exostosin-1 (EXT-1), a tumor suppressor, is an intermediate for malignant tumors. Understanding the mechanism behind the interaction between AEG-1 and EXT-1 may provide insights into colon cancer metastasis.

METHODS

AOM/DSS was used to induce tumor in BALB/c mice. Using an in vivo-jetPEI transfection reagent, transient transfection of AEG-1 and EXT-1 siRNAs were achieved. Histological scoring, immunohistochemical staining, and gene expression studies were performed from excised tissues. Data from the Cancer Genomic Atlas and GEO databases were obtained to identify the expression status of AEG-1 and itsassociation with the survival.

RESULTS

In BALB/c mice, the AOM+DSS treated mice developed necrotic, inflammatory and dysplastic changes in the colon with definite clinical symptoms such as loss of goblet cells, colon shortening, and collagen deposition. Administration of AEG-1 siRNA resulted in a substantial decrease in the disease activity index. Mice treated with EXT-1 siRNA showed diffusely reduced goblet cells. In vivo investigations revealed that PTCH-1 activity was influenced by upstream gene AEG-1, which in turn may affect EXT-1 activity. Data from The Cancer Genomic Atlas and GEO databases confirmed the upregulation of AEG-1 and downregulation of EXT-1 in cancer patients.

CONCLUSIONS

This study revealed that AEG-1 silencing might alter EXT-1 expression indirectly through PTCH-1, influencing cell-ECM interactions, and decreasing dysplastic changes, proliferation and invasion.

ENLIGHTENMENT: A Scalable Annotated Database of Genomics and NGS-Based Nucleotide Level Profiles

The revolution in sequencing technologies has enabled human genomes to be sequenced at a very low cost and time leading to exponential growth in the availability of whole-genome sequences. However, the complete understanding of our genome and its association with cancer is a far way to go. Researchers are striving hard to detect new variants and find their association with diseases, which further gives rise to the need for aggregation of this Big Data into a common standard scalable platform. In this work, a database named Enlightenment has been implemented which makes the availability of genomic data integrated from eight public databases, and DNA sequencing profiles of H. sapiens in a single platform. Annotated results with respect to cancer specific biomarkers, pharmacogenetic biomarkers and its association with variability in drug response, and DNA profiles along with novel copy number variants are computed and stored, which are accessible through a web interface. In order to overcome the challenge of storage and processing of NGS technology-based whole-genome DNA sequences, Enlightenment has been extended and deployed to a flexible and horizontally scalable database HBase, which is distributed over a hadoop cluster, which would enable the integration of other omics data into the database for enlightening the path towards eradication of cancer.

Shedding light on the dark genome: Insights into the genetic, CRISPR-based, and pharmacological dependencies of human cancers and disease aggressiveness

Investigating the human genome is vital for identifying risk factors and devising effective therapies to combat genetic disorders and cancer. Despite the extensive knowledge of the "light genome", the poorly understood "dark genome" remains understudied. In this study, we integrated data from 20,412 protein-coding genes in Pharos and 8,395 patient-derived tumours from The Cancer Genome Atlas (TCGA) to examine the genetic and pharmacological dependencies in human cancers and their treatment implications. We discovered that dark genes exhibited high mutation rates in certain cancers, similar to light genes. By combining the drug response profiles of cancer cells with cell fitness post-CRISPR-mediated gene knockout, we identified the crucial vulnerabilities associated with both dark and light genes. Our analysis also revealed that tumours harbouring dark gene mutations displayed worse overall and disease-free survival rates than those without such mutations. Furthermore, dark gene expression levels significantly influenced patient survival outcomes. Our findings demonstrated a similar distribution of genetic and pharmacological dependencies across the light and dark genomes, suggesting that targeting the dark genome holds promise for cancer treatment. This study underscores the need for ongoing research on the dark genome to better comprehend the underlying mechanisms of cancer and develop more effective therapies.

Transcriptome reprogramming through alternative splicing triggered by apigenin drives cell death in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is characterized by its aggressiveness and resistance to cancer-specific transcriptome alterations. Alternative splicing (AS) is a major contributor to the diversification of cancer-specific transcriptomes. The TNBC transcriptome landscape is characterized by aberrantly spliced isoforms that promote tumor growth and resistance, underscoring the need to identify approaches that reprogram AS circuitry towards transcriptomes, favoring a delay in tumorigenesis or responsiveness to therapy. We have previously shown that flavonoid apigenin is associated with splicing factors, including heterogeneous nuclear ribonucleoprotein A2 (hnRNPA2). Here, we showed that apigenin reprograms TNBC-associated AS transcriptome-wide. The AS events affected by apigenin were statistically enriched in hnRNPA2 substrates. Comparative transcriptomic analyses of human TNBC tumors and non-tumor tissues showed that apigenin can switch cancer-associated alternative spliced isoforms (ASI) to those found in non-tumor tissues. Apigenin preferentially affects the splicing of anti-apoptotic and proliferation factors, which are uniquely observed in cancer cells, but not in non-tumor cells. Apigenin switches cancer-associated aberrant ASI in vivo in TNBC xenograft mice by diminishing proliferation and increasing pro-apoptotic ASI. In accordance with these findings, apigenin increased apoptosis and reduced tumor proliferation, thereby halting TNBC growth in vivo. Our results revealed that apigenin reprograms transcriptome-wide TNBC-specific AS, thereby inducing apoptosis and hindering tumor growth. These findings underscore the impactful effects of nutraceuticals in altering cancer transcriptomes, offering new options to influence outcomes in TNBC treatments.

EZH2 Cooperates with BRD4-NUT to Drive NUT Carcinoma Growth by Silencing Key Tumor Suppressor Genes

NUT carcinoma is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of progrowth genes. BET bromodomain inhibitors (BETi) are a promising treatment for NUT carcinoma that can impede BRD4-NUT's ability to activate genes, but the efficacy of BETi as monotherapy is limited. Here, we demonstrated that enhancer of zeste homolog 2 (EZH2), which silences genes through establishment of repressive chromatin, is a dependency in NUT carcinoma. Inhibition of EZH2 with the clinical compound tazemetostat potently blocked growth of NUT carcinoma cells. Epigenetic and transcriptomic analysis revealed that tazemetostat reversed the EZH2-specific H3K27me3 silencing mark and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4-NUT-regulated genes. Indeed, H3K27me3 and H3K27ac domains were found to be mutually exclusive in NUT carcinoma cells. CDKN2A was identified as the only gene among all tazemetostat-derepressed genes to confer resistance to tazemetostat in a CRISPR-Cas9 screen. Combined inhibition of EZH2 and BET synergized to downregulate cell proliferation genes, resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In preclinical models, combined tazemetostat and BETi synergistically blocked tumor growth and prolonged survival of NUT carcinoma-xenografted mice, with complete remission without relapse in one cohort. Identification of EZH2 as a dependency in NUT carcinoma substantiates the reliance of NUT carcinoma tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary, chromatin regulatory pathways to maintain NUT carcinoma growth.

SIGNIFICANCE

Repression of tumor suppressor genes, including CDKN2A, by EZH2 provides a mechanistic rationale for combining EZH2 and BET inhibitors for the clinical treatment of NUT carcinoma. See related commentary by Kazansky and Kentsis, p. 3827.

Parvimonas micra, an oral pathobiont associated with colorectal cancer, epigenetically reprograms human colonocytes

Recently, an intestinal dysbiotic microbiota with enrichment in oral cavity bacteria has been described in colorectal cancer (CRC) patients. Here, we characterize and investigate one of these oral pathobionts, the Gram-positive anaerobic coccus Parvimonas micra. We identified two phylotypes (A and B) exhibiting different phenotypes and adhesion capabilities. We observed a strong association of phylotype A with CRC, with its higher abundance in feces and in tumoral tissue compared with the normal homologous colonic mucosa, which was associated with a distinct methylation status of patients. By developing an in vitro hypoxic co-culture system of human primary colonic cells with anaerobic bacteria, we show that P. micra phylotype A alters the DNA methylation profile promoters of key tumor-suppressor genes, oncogenes, and genes involved in epithelial-mesenchymal transition. In colonic mucosa of CRC patients carrying P. micra phylotype A, we found similar DNA methylation alterations, together with significant enrichment of differentially expressed genes in pathways involved in inflammation, cell adhesion, and regulation of actin cytoskeleton, providing evidence of P. micra's possible role in the carcinogenic process.

Identifying the driver miRNAs with somatic copy number alterations driving dysregulated ceRNA networks in cancers

BACKGROUND

MicroRNAs (miRNAs) play critical roles in cancer initiation and progression, which were critical components to maintain the dynamic balance of competing endogenous RNA (ceRNA) networks. Somatic copy number alterations (SCNAs) in the cancer genome could disturb the transcriptome level of miRNA to deregulate this balance. However, the driving effects of SCNAs of miRNAs were insufficiently understood.

METHODS

In this study, we proposed a method to dissect the functional roles of miRNAs under different copy number states and identify driver miRNAs by integrating miRNA SCNAs profile, miRNA-target relationships and expression profiles of miRNA, mRNA and lncRNA.

RESULTS

Applying our method to 813 TCGA breast cancer (BRCA) samples, we identified 29 driver miRNAs whose SCNAs significantly and concordantly regulated their own expression levels and further inversely dysregulated expression levels of their targets or disturbed the miRNA-target networks they directly involved. Based on miRNA-target networks, we further constructed dynamic ceRNA networks driven by driver SCNAs of miRNAs and identified three different patterns of SCNA interference in the miRNA-mediated dynamic ceRNA networks. Survival analysis of driver miRNAs showed that high-level amplifications of four driver miRNAs (including has-miR-30d-3p, has-mir-30b-5p, has-miR-30d-5p and has-miR-151a-3p) in 8q24 characterized a new BRCA subtype with poor prognosis and contributed to the dysfunction of cancer-associated hallmarks in a complementary way. The SCNAs of driver miRNAs across different cancer types contributed to the cancer development by dysregulating different components of the same cancer hallmarks, suggesting the cancer specificity of driver miRNA.

CONCLUSIONS

These results demonstrate the efficacy of our method in identifying driver miRNAs and elucidating their functional roles driven by endogenous SCNAs, which is useful for interpreting cancer genomes and pathogenic mechanisms.

Spatial proteo-transcriptomic profiling reveals the molecular landscape of borderline ovarian tumors and their invasive progression

Serous borderline tumors (SBT) are epithelial neoplastic lesions of the ovaries that commonly have a good prognosis. In 10-15% of cases, however, SBT will recur as low-grade serous cancer (LGSC), which is deeply invasive and responds poorly to current standard chemotherapy1,2,3. While genetic alterations suggest a common origin, the transition from SBT to LGSC remains poorly understood4. Here, we integrate spatial proteomics5 with spatial transcriptomics to elucidate the evolution from SBT to LGSC and its corresponding metastasis at the molecular level in both the stroma and the tumor. We show that the transition of SBT to LGSC occurs in the epithelial compartment through an intermediary stage with micropapillary features (SBT-MP), which involves a gradual increase in MAPK signaling. A distinct subset of proteins and transcripts was associated with the transition to invasive tumor growth, including the neuronal splicing factor NOVA2, which was limited to expression in LGSC and its corresponding metastasis. An integrative pathway analysis exposed aberrant molecular signaling of tumor cells supported by alterations in angiogenesis and inflammation in the tumor microenvironment. Integration of spatial transcriptomics and proteomics followed by knockdown of the most altered genes or pharmaceutical inhibition of the most relevant targets confirmed their functional significance in regulating key features of invasiveness. Combining cell-type resolved spatial proteomics and transcriptomics allowed us to elucidate the sequence of tumorigenesis from SBT to LGSC. The approach presented here is a blueprint to systematically elucidate mechanisms of tumorigenesis and find novel treatment strategies.

Identifying the personalized driver gene sets maximally contributing to abnormality of transcriptome phenotype in glioblastoma multiforme individuals

High heterogeneity in genome and phenotype of cancer populations made it difficult to apply population-based common driver genes to the diagnosis and treatment of cancer individuals. Characterizing and identifying the personalized driver mechanism for glioblastoma multiforme (GBM) individuals were pivotal for the realization of precision medicine. We proposed an integrative method to identify the personalized driver gene sets by integrating the profiles of gene expression and genetic alterations in cancer individuals. This method coupled genetic algorithm and random walk to identify the optimal gene sets that could explain abnormality of transcriptome phenotype to the maximum extent. The personalized driver gene sets were identified for 99 GBM individuals using our method. We found that genomic alterations in between one and seven driver genes could maximally and cumulatively explain the dysfunction of cancer hallmarks across GBM individuals. The driver gene sets were distinct even in GBM individuals with significantly similar transcriptomic phenotypes. Our method identified MCM4 with rare genetic alterations as previously unknown oncogenic genes, the high expression of which were significantly associated with poor GBM prognosis. The functional experiments confirmed that knockdown of MCM4 could significantly inhibit proliferation, invasion, migration, and clone formation of the GBM cell lines U251 and U118MG, and overexpression of MCM4 significantly promoted the proliferation, invasion, migration, and clone formation of the GBM cell line U87MG. Our method could dissect the personalized driver genetic alteration sets that are pivotal for developing targeted therapy strategies and precision medicine. Our method could be extended to identify key drivers from other levels and could be applied to more cancer types.

Landscape of enhancer disruption and functional screen in melanoma cells

BACKGROUND

The high mutation rate throughout the entire melanoma genome presents a major challenge in stratifying true driver events from the background mutations. Numerous recurrent non-coding alterations, such as those in enhancers, can shape tumor evolution, thereby emphasizing the importance in systematically deciphering enhancer disruptions in melanoma.

RESULTS

Here, we leveraged 297 melanoma whole-genome sequencing samples to prioritize highly recurrent regions. By performing a genome-scale CRISPR interference (CRISPRi) screen on highly recurrent region-associated enhancers in melanoma cells, we identified 66 significant hits which could have tumor-suppressive roles. These functional enhancers show unique mutational patterns independent of classical significantly mutated genes in melanoma. Target gene analysis for the essential enhancers reveal many known and hidden mechanisms underlying melanoma growth. Utilizing extensive functional validation experiments, we demonstrate that a super enhancer element could modulate melanoma cell proliferation by targeting MEF2A, and another distal enhancer is able to sustain PTEN tumor-suppressive potential via long-range interactions.

CONCLUSIONS

Our study establishes a catalogue of crucial enhancers and their target genes in melanoma growth and progression, and illuminates the identification of novel mechanisms of dysregulation for melanoma driver genes and new therapeutic targeting strategies.

Comprehensive single cell aging atlas of mammary tissues reveals shared epigenomic and transcriptomic signatures of aging and cancer

Aging is the greatest risk factor for breast cancer; however, how age-related cellular and molecular events impact cancer initiation is unknown. We investigate how aging rewires transcriptomic and epigenomic programs of mouse mammary glands at single cell resolution, yielding a comprehensive resource for aging and cancer biology. Aged epithelial cells exhibit epigenetic and transcriptional changes in metabolic, pro-inflammatory, or cancer-associated genes. Aged stromal cells downregulate fibroblast marker genes and upregulate markers of senescence and cancer-associated fibroblasts. Among immune cells, distinct T cell subsets (Gzmk+, memory CD4+, γδ) and M2-like macrophages expand with age. Spatial transcriptomics reveal co-localization of aged immune and epithelial cells in situ. Lastly, transcriptional signatures of aging mammary cells are found in human breast tumors, suggesting mechanistic links between aging and cancer. Together, these data uncover that epithelial, immune, and stromal cells shift in proportions and cell identity, potentially impacting cell plasticity, aged microenvironment, and neoplasia risk.

Synthetic viability induces resistance to immune checkpoint inhibitors in cancer cells

BACKGROUND

Immune checkpoint inhibitors (ICI) have revolutionized the treatment for multiple cancers. However, most of patients encounter resistance. Synthetic viability (SV) between genes could induce resistance. In this study, we established SV signature to predict the efficacy of ICI treatment for melanoma.

METHODS

We collected features and predicted SV gene pairs by random forest classifier. This work prioritized SV gene pairs based on CRISPR/Cas9 screens. SV gene pairs signature were constructed to predict the response to ICI for melanoma patients.

RESULTS

This study predicted robust SV gene pairs based on 14 features. Filtered by CRISPR/Cas9 screens, we identified 1,861 SV gene pairs, which were also related with prognosis across multiple cancer types. Next, we constructed the six SV pairs signature to predict resistance to ICI for melanoma patients. This study applied the six SV pairs signature to divide melanoma patients into high-risk and low-risk. High-risk melanoma patients were associated with worse response after ICI treatment. Immune landscape analysis revealed that high-risk melanoma patients had lower natural killer cells and CD8+ T cells infiltration.

CONCLUSIONS

In summary, the 14 features classifier accurately predicted robust SV gene pairs for cancer. The six SV pairs signature could predict resistance to ICI.

Functional Profiling of Soft Tissue Sarcoma Using Mechanistic Models

Soft tissue sarcoma is an umbrella term for a group of rare cancers that are difficult to treat. In addition to surgery, neoadjuvant chemotherapy has shown the potential to downstage tumors and prevent micrometastases. However, finding effective therapeutic targets remains a research challenge. Here, a previously developed computational approach called mechanistic models of signaling pathways has been employed to unravel the impact of observed changes at the gene expression level on the ultimate functional behavior of cells. In the context of such a mechanistic model, RNA-Seq counts sourced from the Recount3 resource, from The Cancer Genome Atlas (TCGA) Sarcoma project, and non-diseased sarcomagenic tissues from the Genotype-Tissue Expression (GTEx) project were utilized to investigate signal transduction activity through signaling pathways. This approach provides a precise view of the relationship between sarcoma patient survival and the signaling landscape in tumors and their environment. Despite the distinct regulatory alterations observed in each sarcoma subtype, this study identified 13 signaling circuits, or elementary sub-pathways triggering specific cell functions, present across all subtypes, belonging to eight signaling pathways, which served as predictors for patient survival. Additionally, nine signaling circuits from five signaling pathways that highlighted the modifications tumor samples underwent in comparison to normal tissues were found. These results describe the protective role of the immune system, suggesting an anti-tumorigenic effect in the tumor microenvironment, in the process of tumor cell detachment and migration, or the dysregulation of ion homeostasis. Also, the analysis of signaling circuit intermediary proteins suggests multiple strategies for therapy.

Evolutionary signatures of human cancers revealed via genomic analysis of over 35,000 patients

Recurring sequences of genomic alterations occurring across patients can highlight repeated evolutionary processes with significant implications for predicting cancer progression. Leveraging the ever-increasing availability of cancer omics data, here we unveil cancer's evolutionary signatures tied to distinct disease outcomes, representing "favored trajectories" of acquisition of driver mutations detected in patients with similar prognosis. We present a framework named ASCETIC (Agony-baSed Cancer EvoluTion InferenCe) to extract such signatures from sequencing experiments generated by different technologies such as bulk and single-cell sequencing data. We apply ASCETIC to (i) single-cell data from 146 myeloid malignancy patients and bulk sequencing from 366 acute myeloid leukemia patients, (ii) multi-region sequencing from 100 early-stage lung cancer patients, (iii) exome/genome data from 10,000+ Pan-Cancer Atlas samples, and (iv) targeted sequencing from 25,000+ MSK-MET metastatic patients, revealing subtype-specific single-nucleotide variant signatures associated with distinct prognostic clusters. Validations on several datasets underscore the robustness and generalizability of the extracted signatures.

Long-Read Sequencing Reveals Rapid Evolution of Immunity- and Cancer-Related Genes in Bats

Bats are exceptional among mammals for their powered flight, extended lifespans, and robust immune systems and therefore have been of particular interest in comparative genomics. Using the Oxford Nanopore Technologies long-read platform, we sequenced the genomes of two bat species with key phylogenetic positions, the Jamaican fruit bat (Artibeus jamaicensis) and the Mesoamerican mustached bat (Pteronotus mesoamericanus), and carried out a comprehensive comparative genomic analysis with a diverse collection of bats and other mammals. The high-quality, long-read genome assemblies revealed a contraction of interferon (IFN)-α at the immunity-related type I IFN locus in bats, resulting in a shift in relative IFN-ω and IFN-α copy numbers. Contradicting previous hypotheses of constitutive expression of IFN-α being a feature of the bat immune system, three bat species lost all IFN-α genes. This shift to IFN-ω could contribute to the increased viral tolerance that has made bats a common reservoir for viruses that can be transmitted to humans. Antiviral genes stimulated by type I IFNs also showed evidence of rapid evolution, including a lineage-specific duplication of IFN-induced transmembrane genes and positive selection in IFIT2. In addition, 33 tumor suppressors and 6 DNA-repair genes showed signs of positive selection, perhaps contributing to increased longevity and reduced cancer rates in bats. The robust immune systems of bats rely on both bat-wide and lineage-specific evolution in the immune gene repertoire, suggesting diverse immune strategies. Our study provides new genomic resources for bats and sheds new light on the extraordinary molecular evolution in this critically important group of mammals.

Methylation-related genes involved in renal carcinoma progression

Renal carcinomas are a group of malignant tumors often originating in the cells lining the small tubes in the kidney responsible for filtering waste from the blood and urine production. Kidney tumors arise from the uncontrolled growth of cells in the kidneys and are responsible for a large share of global cancer-related morbidity and mortality. Understanding the molecular mechanisms driving renal carcinoma progression results crucial for the development of targeted therapies leading to an improvement of patient outcomes. Epigenetic mechanisms such as DNA methylation are known factors underlying the development of several cancer types. There is solid experimental evidence of relevant biological functions modulated by methylation-related genes, associated with the progression of different carcinomas. Those mechanisms can often be associated to different epigenetic marks, such as DNA methylation sites or chromatin conformation patterns. Currently, there is no definitive method to establish clear relations between genetic and epigenetic factors that influence the progression of cancer. Here, we developed a data-driven method to find methylation-related genes, so we could find relevant bonds between gene co-expression and methylation-wide-genome regulation patterns able to drive biological processes during the progression of clear cell renal carcinoma (ccRC). With this approach, we found out genes such as ITK oncogene that appear hypomethylated during all four stages of ccRC progression and are strongly involved in immune response functions. Also, we found out relevant tumor suppressor genes such as RAB25 hypermethylated, thus potentially avoiding repressed functions in the AKT signaling pathway during the evolution of ccRC. Our results have relevant implications to further understand some epigenetic-genetic-affected roles underlying the progression of renal cancer.

EZH2 synergizes with BRD4-NUT to drive NUT carcinoma growth through silencing of key tumor suppressor genes

NUT carcinoma (NC) is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of pro-growth genes. BET bromodomain inhibitors (BETi) impede BRD4-NUT's ability to activate genes and are thus a promising treatment but limited as monotherapy. The role of gene repression in NC is unknown. Here, we demonstrate that EZH2, which silences genes through establishment of repressive chromatin, is a dependency in NC. Inhibition of EZH2 with the clinical compound tazemetostat (taz) potently blocked growth of NC cells. Epigenetic and transcriptomic analysis revealed that taz reversed the EZH2-specific H3K27me3 silencing mark, and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4- NUT-regulated genes. CDKN2A was identified as the only gene amongst all taz-derepressed genes to confer resistance to taz in a CRISPR-Cas9 screen. Combined EZH2 inhibition and BET inhibition synergized to downregulate cell proliferation genes resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In pre-clinical models, combined taz and BETi synergistically blocked growth and prolonged survival of NC-xenografted mice, with all mice cured in one cohort.

STATEMENT OF SIGNIFICANCE

Identification of EZH2 as a dependency in NC substantiates the reliance of NC tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary chromatin regulatory pathways to maintain NC growth. In particular, repression of CDKN2A expression by EZH2 provides a mechanistic rationale for combining EZH2i with BETi for the clinical treatment of NC.

Mutational landscape of cancer-driver genes across human cancers

The genetic mutations that contribute to the transformation of healthy cells into cancerous cells have been the subject of extensive research. The molecular aberrations that lead to cancer development are often characterised by gain-of-function or loss-of-function mutations in a variety of oncogenes and tumour suppressor genes. In this study, we investigate the genomic sequences of 20,331 primary tumours representing 41 distinct human cancer types to identify and catalogue the driver mutations present in 727 known cancer genes. Our findings reveal significant variations in the frequency of cancer gene mutations across different cancer types and highlight the frequent involvement of tumour suppressor genes (94%), oncogenes (93%), transcription factors (72%), kinases (64%), cell surface receptors (63%), and phosphatases (22%), in cancer. Additionally, our analysis reveals that cancer gene mutations are predominantly co-occurring rather than exclusive in all types of cancer. Notably, we discover that patients with tumours displaying different combinations of gene mutation patterns tend to exhibit variable survival outcomes. These findings provide new insights into the genetic landscape of cancer and bring us closer to a comprehensive understanding of the underlying mechanisms driving the development of various forms of cancer.

Synthetic Reader-Actuators Targeted to Polycomb-Silenced Genes Block Triple-Negative Breast Cancer Proliferation and Invasion

Scientists have used pharmacological inhibitors of polycomb proteins to restore the expression of tumor suppressor genes and stop cancer proliferation and invasion. A major limitation of this approach is that key transcriptional activators, such as TP53 and BAF SWI/SNF, are often mutated in cancer. Poor clinical results for polycomb-targeting therapies in solid cancers, including triple-negative breast cancer (TNBC), could discourage the further development of epigenetic monotherapies. Here, we performed epigenome actuation with a synthetic reader-actuator (SRA) that binds trimethylated histone H3 lysine 27 in polycomb chromatin and modulates core transcriptional activators. In SRA-expressing TNBC BT-549 cells, 122 genes become upregulated ≥2-fold, including the genes involved in cell death, cell cycle arrest, and migration inhibition. The SRA-expressing spheroids showed reduced size in Matrigel and loss of invasion. Therefore, targeting Mediator-recruiting regulators to silenced chromatin can activate tumor suppressors and stimulate anti-cancer phenotypes, and further development of robust gene regulators might benefit TNBC patients.

Identification of signature genes and drug candidates for primary plasma cell leukemia: An integrated system biology approach

BACKGROUND

Plasma cell leukemia (PCL) is one of the rare cancer which is characterized by the uncontrolled proliferation of plasma cells in peripheral blood and bone marrow. The aggressive behavior of the disease and high mortality rate among PCL patients makes it a thirst area to be explored.

METHODS

The dataset for PCL was obtained from the GEO database and was analyzed using GEO2R for differentially expressed genes. Further, the functional enrichment analysis was carried out for DEGs using DAVID. The protein-protein interactions (PPI) for DEGs were obtained using STRING 11.5 and were analyzed in Cytoscape 3.7.2. to obtain the key hub genes. These key hub genes were investigated for their interaction with suitable drug candidates using DGIdb, DrugMAP, and Schrodinger's version 2022-1.

RESULTS

Out of the total of 104 DEGs, 39 genes were up-regulated whereas 65 genes were down-regulated. A total of 11 biological processes, 2 cellular components, and 5 molecular functions were enriched along with the 7 KEGG pathways for the DEGs. Further, a total of 11 hub genes were obtained from the PPI of DEGs of which TP53, MAPK1, SOCS1, MBD3, and YES1 were the key hub genes. Oxaliplatin, mitoxantrone, and ponatinib were found to have the highest binding affinity towards the p53, MAPK1, and YES1 proteins respectively.

CONCLUSION

TP53, MAPK1, SOCS1, MBD3, and YES1 are the signature hub genes that might be responsible for the aggressive prognosis of PCL leading to poor survival rate. However, p53, MAPK1, and YES1 can be targeted with oxaliplatin, mitoxantrone, and ponatinib.

Register-based and genetic studies of Prader-Willi syndrome show a high frequency of gonadal tumors and a possible mechanism for tumorigenesis through imprinting relaxation

Prader-Willi syndrome (PWS) is a rare disease caused by a lack of expression of inherited imprinted genes in the paternally derived Prader-Willi critical region on chromosome 15q11.2-q13. It is characterized by poor feeding and hypotonia in infancy, intellectual disability, behavioral abnormalities, dysmorphic features, short stature, obesity, and hypogonadism. PWS is not a known cancer predisposition syndrome, but previous investigations regarding the prevalence of cancer in these patients suggest an increased risk of developing specific cancer types such as myeloid leukemia and testicular cancer. We present the results from a Swedish national population-based cohort study of 360 individuals with PWS and 18,000 matched comparisons. The overall frequency of cancer was not increased in our PWS cohort, but we found a high frequency of pediatric cancers. We also performed whole-genome sequencing of blood- and tumor-derived DNAs from a unilateral dysgerminoma in a 13-year-old girl with PWS who also developed bilateral ovarian sex cord tumors with annular tubules. In germline analysis, there were no additional findings apart from the 15q11.2-q13 deletion of the paternal allele, while a pathogenic activating KIT mutation was identified in the tumor. Additionally, methylation-specific multiplex ligation-dependent probe amplification revealed reduced methylation at the PWS locus in the dysgerminoma but not in the blood. In conclusion, our register-based study suggests an increased risk of cancer at a young age, especially testicular and ovarian tumors. We found no evidence of a general increase in cancer risk in patients with PWS. However, given our limited observational time, further studies with longer follow-up times are needed to clarify the lifetime cancer risk in PWS. We have also described the second case of locus-specific loss-of-imprinting in a germ cell tumor in PWS, suggesting a possible mechanism of carcinogenesis.

IGFBP2 secretion by mammary adipocytes limits breast cancer invasion

The progression of noninvasive ductal carcinoma in situ to invasive ductal carcinoma for patients with breast cancer results in a significantly poorer prognosis and is the precursor to metastatic disease. In this work, we have identified insulin-like growth factor-binding protein 2 (IGFBP2) as a potent adipocrine factor secreted by healthy breast adipocytes that acts as a barrier against invasive progression. In line with this role, adipocytes differentiated from patient-derived stromal cells were found to secrete IGFBP2, which significantly inhibited breast cancer invasion. This occurred through binding and sequestration of cancer-derived IGF-II. Moreover, depletion of IGF-II in invading cancer cells using small interfering RNAs or an IGF-II-neutralizing antibody ablated breast cancer invasion, highlighting the importance of IGF-II autocrine signaling for breast cancer invasive progression. Given the abundance of adipocytes in the healthy breast, this work exposes the important role they play in suppressing cancer progression and may help expound upon the link between increased mammary density and poorer prognosis.

Pan-cancer analysis of whole-genome doubling and its association with patient prognosis

BACKGROUND

Whole-genome doubling (WGD) is a common mutation in cancer. Various studies have suggested that WGD is associated with a poor prognosis in cancer. However, the detailed association between WGD occurrence and prognosis remains unclear. In this study, we aimed to elucidate the mechanism by which WGD affects prognosis using sequencing data from the Pan-Cancer Analysis of Whole Genomes (PCAWG) and The Cancer Genome Atlas.

METHODS

Whole-genome sequencing data of 23 cancer types were downloaded from PCAWG project. We defined the WGD event in each sample using the WGD status annotated using PCAWG. We used MutationTimeR to predict the relative timings of mutations and loss of heterozygosity (LOH) in WGD, thus evaluating their association with WGD. We also analyzed the association between WGD-associated factors and patient prognosis.

RESULTS

WGD was associated with several factors, e.g., length of LOH regions. Survival analysis using WGD-associated factors revealed that longer LOH regions and LOH in chr17 were associated with poor prognosis in samples with WGD (WGD samples) and samples without WGD (nWGD samples). In addition to these two factors, nWGD samples showed that the number of mutations in tumor suppressor genes was associated with prognosis. Moreover, we explored the genes associated with prognosis in both samples separately.

CONCLUSION

The prognosis-related factors in WGD samples differed significantly compared with those in nWGD samples. This study emphasizes the need for different treatment strategies for WGD and nWGD samples.

Gene expression networks involved in multiple cellular programs coexist in individual hepatocellular cancer cells

The gene expression networks of a single cell can be used to reveal cell type- and condition-specific patterns that account for cell states, cell identity, and its responses to environmental changes. We applied single cell sequencing datasets to define mRNA patterns and visualized potential cellular capacities among hepatocellular cancer cells. The expressing numbers and levels of genes were highly heterogenous among the cancer cells. The cellular characteristics were dependent strongly on the expressing numbers and levels of genes, especially oncogenes and anti-oncogenes, in an individual cancer cell. The transcriptional activations of oncogenes and anti-oncogenes were strongly linked to inherent multiple cellular programs, some of which oppose and contend against other processes, in a cancer cell. The gene expression networks of multiple cellular programs proliferation, differentiation, apoptosis, autophagy, epithelial-mesenchymal transition, ATP production, and neurogenesis coexisted in an individual cancer cell. The findings give rise a hypothesis that a cancer cell expresses balanced combinations of genes and undergoes a given biological process by rapidly transmuting gene expressing networks.

Identifying the Common Cell-Free DNA Biomarkers across Seven Major Cancer Types

Blood-based detection of circulating cell-free DNA (cfDNA) is a non-invasive and easily accessible method for early cancer detection. Despite the extensive utility of cfDNA, there are still many challenges to developing clinical biomarkers. For example, cfDNA with genetic alterations often composes a small portion of the DNA circulating in plasma, which can be confounded by cfDNA contributed by normal cells. Therefore, filtering out the potential false-positive cfDNA mutations from healthy populations will be important for cancer-based biomarkers. Additionally, many low-frequency genetic alterations are easily overlooked in a small number of cfDNA-based cancer tests. We hypothesize that the combination of diverse types of cancer studies on cfDNA will provide us with a new perspective on the identification of low-frequency genetic variants across cancer types for promoting early diagnosis. By building a standardized computational pipeline for 1358 cfDNA samples across seven cancer types, we prioritized 129 shard genetic variants in the major cancer types. Further functional analysis of the 129 variants found that they are mainly enriched in ribosome pathways such as cotranslational protein targeting the membrane, some of which are tumour suppressors, oncogenes, and genes related to cancer initiation. In summary, our integrative analysis revealed the important roles of ribosome proteins as common biomarkers in early cancer diagnosis.

Detection of germline variants with pathogenic potential in 48 patients with familial colorectal cancer by using whole exome sequencing

BACKGROUND

Hereditary genetic mutations causing predisposition to colorectal cancer are accountable for approximately 30% of all colorectal cancer cases. However, only a small fraction of these are high penetrant mutations occurring in DNA mismatch repair genes, causing one of several types of familial colorectal cancer (CRC) syndromes. Most of the mutations are low-penetrant variants, contributing to an increased risk of familial colorectal cancer, and they are often found in additional genes and pathways not previously associated with CRC. The aim of this study was to identify such variants, both high-penetrant and low-penetrant ones.

METHODS

We performed whole exome sequencing on constitutional DNA extracted from blood of 48 patients suspected of familial colorectal cancer and used multiple in silico prediction tools and available literature-based evidence to detect and investigate genetic variants.

RESULTS

We identified several causative and some potentially causative germline variants in genes known for their association with colorectal cancer. In addition, we identified several variants in genes not typically included in relevant gene panels for colorectal cancer, including CFTR, PABPC1 and TYRO3, which may be associated with an increased risk for cancer.

CONCLUSIONS

Identification of variants in additional genes that potentially can be associated with familial colorectal cancer indicates a larger genetic spectrum of this disease, not limited only to mismatch repair genes. Usage of multiple in silico tools based on different methods and combined through a consensus approach increases the sensitivity of predictions and narrows down a large list of variants to the ones that are most likely to be significant.

The Integrative Studies on the Functional A-to-I RNA Editing Events in Human Cancers

Adenosine-to-inosine (A-to-I) RNA editing, constituting nearly 90% of all RNA editing events in humans, has been reported to contribute to the tumorigenesis in diverse cancers. However, the comprehensive map for functional A-to-I RNA editing events in cancers is still insufficient. To fill this gap, we systematically and intensively analyzed multiple tumorigenic mechanisms of A-to-I RNA editing events in samples across 33 cancer types from The Cancer Genome Atlas. For individual candidate among ∼ 1,500,000 quantified RNA editing events, we performed diverse types of downstream functional annotations. Finally, we identified 24,236 potentially functional A-to-I RNA editing events, including the cases in APOL1, IGFBP3, GRIA2, BLCAP, and miR-589-3p. These events might play crucial roles in the scenarios of tumorigenesis, due to their tumor-related editing frequencies or probable effects on altered expression profiles, protein functions, splicing patterns, and microRNA regulations of tumor genes. Our functional A-to-I RNA editing events (https://ccsm.uth.edu/CAeditome/) will help better understand the cancer pathology from the A-to-I RNA editing aspect.

Possible involvement of silent mutations in cancer pathogenesis and evolution

Recent studies have shown that some silent mutations can be harmful to various processes. In this study, we performed a comprehensive in silico analysis to elucidate the effects of silent mutations on cancer pathogenesis using exome sequencing data derived from the Cancer Genome Atlas. We focused on the codon optimality scores of silent mutations, which were defined as the difference between the optimality of synonymous codons, calculated using the codon usage table. The relationship between cancer evolution and silent mutations showed that the codon optimality score of the mutations that occurred later in carcinogenesis was significantly higher than of those that occurred earlier. In addition, mutations with higher scores were enriched in genes involved in the cell cycle and cell division, while those with lower scores were enriched in genes involved in apoptosis and cellular senescence. Our results demonstrate that some silent mutations can be involved in cancer pathogenesis.

Decoding the lncRNAome Across Diverse Cellular Stresses Reveals Core p53-effector Pan-cancer Suppressive lncRNAs

Despite long non-coding RNAs (lncRNAs) emerging as key contributors to malignancies, their transcriptional regulation, tissue-type expression under different conditions, and functions remain largely unknown. Developing a combined computational and experimental framework, which integrates pan-cancer RNAi/CRISPR screens, and genomic, epigenetic, and expression profiles (including single-cell RNA sequencing), we report across multiple cancers, core p53-transcriptionally regulated lncRNAs, which were thought to be primarily cell/tissue-specific. These lncRNAs were consistently directly transactivated by p53 with different cellular stresses in multiple cell types and associated with pan-cancer cell survival/growth suppression and patient survival. Our prediction results were verified through independent validation datasets, our own patient cohort, and cancer cell experiments. Moreover, a top predicted p53-effector tumor-suppressive lncRNA (we termed PTSL) inhibited cell proliferation and colony formation by modulating the G₂ regulatory network, causing G₂ cell-cycle arrest. Therefore, our results elucidated previously unreported, high-confidence core p53-targeted lncRNAs that suppress tumorigenesis across cell types and stresses.

Significance

Identification of pan-cancer suppressive lncRNAs transcriptionally regulated by p53 across different cellular stresses by integrating multilayered high-throughput molecular profiles. This study provides critical new insights into the p53 tumor suppressor by revealing the lncRNAs in the p53 cell-cycle regulatory network and their impact on cancer cell growth and patient survival.

Systematic integration of protein-affecting mutations, gene fusions, and copy number alterations into a comprehensive somatic mutational profile

Somatic mutations occur as random genetic changes in genes through protein-affecting mutations (PAMs), gene fusions, or copy number alterations (CNAs). Mutations of different types can have a similar phenotypic effect (i.e., allelic heterogeneity) and should be integrated into a unified gene mutation profile. We developed OncoMerge to fill this niche of integrating somatic mutations to capture allelic heterogeneity, assign a function to mutations, and overcome known obstacles in cancer genetics. Application of OncoMerge to TCGA Pan-Cancer Atlas increased detection of somatically mutated genes and improved the prediction of the somatic mutation role as either activating or loss of function. Using integrated somatic mutation matrices increased the power to infer gene regulatory networks and uncovered the enrichment of switch-like feedback motifs and delay-inducing feedforward loops. These studies demonstrate that OncoMerge efficiently integrates PAMs, fusions, and CNAs and strengthens downstream analyses linking somatic mutations to cancer phenotypes.