DECtp: Calling Differential Gene Expression Between Cancer and Normal Samples by Integrating Tumor Purity Information

Identifying the key hub genes linked with lung squamous cell carcinoma by examining the differentially expressed and survival genes

Lung Squamous Cell Carcinoma is characterised by significant alterations in RNA expression patterns, and a lack of early symptoms and diagnosis results in poor survival rates. Our study aimed to identify the hub genes involved in LUSC by differential expression analysis and their influence on overall survival rates in patients. Thus, identifying genes with the potential to serve as biomarkers and therapeutic targets. RNA sequence data for LUSC was obtained from TCGA and analysed using R Studio. Survival analysis was performed on DE genes. PPI network and hub gene analysis was performed on survival-relevant genes. Enrichment analysis was conducted on the PPI network to elucidate the functional roles of hub genes. Our analysis identified 2774 DEGs in LUSC patient datasets. Survival analysis revealed 511 genes with a significant impact on patient survival. Among these, 20 hub genes-FN1, ACTB, HGF, PDGFRB, PTEN, SNAI1, TGFBR1, ESR1, SERPINE1, THBS1, PDGFRA, VWF, BMP2, LEP, VTN, PXN, ABL1, ITGA3 and ANXA5-were found to have lower expression levels associated with better patient survival, whereas high expression of SOX2 correlated with longer survival. Enrichment analysis indicated that these hub genes are involved in critical cellular and cancer-related pathways. Our study has identified six key hub genes that are differentially expressed and exhibit significant influence over LUSC patient survival outcomes. Further, in vitro and in vivo studies must be conducted on the key genes for their utilisation as therapeutic targets and biomarkers in LUSC.

Automated Prognosis Marker Assessment in Breast Cancers Using BLEACH&STAIN Multiplexed Immunohistochemistry

Prognostic markers in routine clinical management of breast cancer are often assessed using RNA-based multi-gene panels that depend on fluctuating tumor purity. Multiplex fluorescence immunohistochemistry (mfIHC) holds the potential for an improved risk assessment. To enable automated prognosis marker detection (i.e., progesterone receptor [PR], estrogen receptor [ER], androgen receptor [AR], GATA3, TROP2, HER2, PD-L1, Ki67, TOP2A), a framework for automated breast cancer identification was developed and validated involving thirteen different artificial intelligence analysis steps and an algorithm for cell distance analysis using 11+1-marker-BLEACH&STAIN-mfIHC staining in 1404 invasive breast cancers of no special type (NST). The framework for automated breast cancer detection discriminated normal glands from malignant glands with an accuracy of 98.4%. This approach identified that five (PR, ER, AR, GATA3, PD-L1) of nine biomarkers were associated with prolonged overall survival (p ≤ 0.0095 each) and two of these (PR, AR) were found to be independent risk factors in multivariate analysis (p ≤ 0.0151 each). The combined assessment of PR-ER-AR-GATA3-PD-L1 as a five-marker prognosis score showed strong prognostic relevance (p < 0.0001) and was an independent risk factor in multivariate analysis (p = 0.0034). Automated breast cancer detection in combination with an artificial intelligence-based analysis of mfIHC enables a rapid and reliable analysis of multiple prognostic parameters. The strict limitation of the analysis to malignant cells excludes the impact of fluctuating tumor purity on assay precision.

IRF4 as a novel target involved in malignant transformation of oral submucous fibrosis into oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) in the context of oral submucous fibrosis (OSF) has a high incidence owing to undefined pathogenesis. Identifying key genes and exploring the underlying molecular mechanisms involved in the conversion of OSF into OSCC are in urgent need. Differentially expressed genes (DEGs) between OSCC and OSF were dug from GEO databases and a total of 170 DEGs were acquired. Functional association of DEGs were analyzed by GO and KEGG. Protein-protein interactions (PPIs) analysis was carried out and candidate biomarkers were identified by Gene co-expression analysis and Cox analyses. Hub genes were confirmed by qRT-PCR in tissues and cell lines, of which we found that IRF4 mRNA was successively up-regulated from Normal to OSF and then to OSCC and associated with immune infiltrating levels. In addition, Immunohistochemical (IHC) and Immunofluorescence (IF) assays were conducted to validate the consistent upregulation of IRF4 and the oncogene role of IRF4 in OSF and OSCC at translation level. IRF4 may be indicative biomarker in transformation of OSF into OSCC. High IRF4 expression contribute to increased immune infiltration of OSCC and may provide a novel diagnostic marker for OSCC patients translated from OSF.

Agent Repurposing for the Treatment of Advanced Stage Diffuse Large B-Cell Lymphoma Based on Gene Expression and Network Perturbation Analysis

Over 50% of diffuse large B-cell lymphoma (DLBCL) patients are diagnosed at an advanced stage. Although there are a few therapeutic strategies for DLBCL, most of them are more effective in limited-stage cancer patients. The prognosis of patients with advanced-stage DLBCL is usually poor with frequent recurrence and metastasis. In this study, we aimed to identify gene expression and network differences between limited- and advanced-stage DLBCL patients, with the goal of identifying potential agents that could be used to relieve the severity of DLBCL. Specifically, RNA sequencing data of DLBCL patients at different clinical stages were collected from the cancer genome atlas (TCGA). Differentially expressed genes were identified using DESeq2, and then, weighted gene correlation network analysis (WGCNA) and differential module analysis were performed to find variations between different stages. In addition, important genes were extracted by key driver analysis, and potential agents for DLBCL were identified according to gene-expression perturbations and the Crowd Extracted Expression of Differential Signatures (CREEDS) drug signature database. As a result, 20 up-regulated and 73 down-regulated genes were identified and 79 gene co-expression modules were found using WGCNA, among which, the thistle1 module was highly related to the clinical stage of DLBCL. KEGG pathway and GO enrichment analyses of genes in the thistle1 module indicated that DLBCL progression was mainly related to the NOD-like receptor signaling pathway, neutrophil activation, secretory granule membrane, and carboxylic acid binding. A total of 47 key drivers were identified through key driver analysis with 11 up-regulated key driver genes and 36 down-regulated key diver genes in advanced-stage DLBCL patients. Five genes (MMP1, RAB6C, ACCSL, RGS21 and MOCOS) appeared as hub genes, being closely related to the occurrence and development of DLBCL. Finally, both differentially expressed genes and key driver genes were subjected to CREEDS analysis, and 10 potential agents were predicted to have the potential for application in advanced-stage DLBCL patients. In conclusion, we propose a novel pipeline to utilize perturbed gene-expression signatures during DLBCL progression for identifying agents, and we successfully utilized this approach to generate a list of promising compounds.

Prediction of tumor purity from gene expression data using machine learning

MOTIVATION

Bulk tumor samples used for high-throughput molecular profiling are often an admixture of cancer cells and non-cancerous cells, which include immune and stromal cells. The mixed composition can confound the analysis and affect the biological interpretation of the results, and thus, accurate prediction of tumor purity is critical. Although several methods have been proposed to predict tumor purity using high-throughput molecular data, there has been no comprehensive study on machine learning-based methods for the estimation of tumor purity.

RESULTS

We applied various machine learning models to estimate tumor purity. Overall, the models predicted the tumor purity accurately and showed a high correlation with well-established gold standard methods. In addition, we identified a small group of genes and demonstrated that they could predict tumor purity well. Finally, we confirmed that these genes were mainly involved in the immune system.

AVAILABILITY

The machine learning models constructed for this study are available at https://github.com/BonilKoo/ML_purity.

A novel single-cell based method for breast cancer prognosis

Breast cancer prognosis is challenging due to the heterogeneity of the disease. Various computational methods using bulk RNA-seq data have been proposed for breast cancer prognosis. However, these methods suffer from limited performances or ambiguous biological relevance, as a result of the neglect of intra-tumor heterogeneity. Recently, single cell RNA-sequencing (scRNA-seq) has emerged for studying tumor heterogeneity at cellular levels. In this paper, we propose a novel method, scPrognosis, to improve breast cancer prognosis with scRNA-seq data. scPrognosis uses the scRNA-seq data of the biological process Epithelial-to-Mesenchymal Transition (EMT). It firstly infers the EMT pseudotime and a dynamic gene co-expression network, then uses an integrative model to select genes important in EMT based on their expression variation and differentiation in different stages of EMT, and their roles in the dynamic gene co-expression network. To validate and apply the selected signatures to breast cancer prognosis, we use them as the features to build a prediction model with bulk RNA-seq data. The experimental results show that scPrognosis outperforms other benchmark breast cancer prognosis methods that use bulk RNA-seq data. Moreover, the dynamic changes in the expression of the selected signature genes in EMT may provide clues to the link between EMT and clinical outcomes of breast cancer. scPrognosis will also be useful when applied to scRNA-seq datasets of different biological processes other than EMT.

Pan-cancer transcriptomic analysis dissects immune and proliferative functions of APOBEC3 cytidine deaminases

APOBEC3 cytidine deaminases are largely known for their innate immune protection from viral infections. Recently, members of the family have been associated with a distinct mutational activity in some cancer types. We report a pan-tissue, pan-cancer analysis of RNA-seq data specific to the APOBEC3 genes in 8,951 tumours, 786 cancer cell lines and 6,119 normal tissues. By deconvolution of levels of different cell types in tumour admixtures, we demonstrate that APOBEC3B (A3B), the primary candidate as a cancer mutagen, shows little association with immune cell types compared to its paralogues. We present a pipeline called RESPECTEx (REconstituting SPecific Cell-Type Expression) and use it to deconvolute cell-type specific expression levels in a given cohort of tumour samples. We functionally annotate APOBEC3 co-expressing genes, and create an interactive visualization tool which 'barcodes' the functional enrichment (http://fraternalilab.kcl.ac.uk/apobec-barcodes/). These analyses reveal that A3B expression correlates with cell cycle and DNA repair genes, whereas the other APOBEC3 members display specificity for immune processes and immune cell populations. We offer molecular insights into the functions of individual APOBEC3 proteins in antiviral and proliferative contexts, and demonstrate the diversification this family of enzymes displays at the transcriptomic level, despite their high similarity in protein sequences and structures.