A Multi-Omics Network of a Seven-Gene Prognostic Signature for Non-Small Cell Lung Cancer

MicroRNA-Based Discovery of Biomarkers, Therapeutic Targets, and Repositioning Drugs for Breast Cancer

Breast cancer treatment can be improved with biomarkers for early detection and individualized therapy. A set of 86 microRNAs (miRNAs) were identified to separate breast cancer tumors from normal breast tissues (n = 52) with an overall accuracy of 90.4%. Six miRNAs had concordant expression in both tumors and breast cancer patient blood samples compared with the normal control samples. Twelve miRNAs showed concordant expression in tumors vs. normal breast tissues and patient survival (n = 1093), with seven as potential tumor suppressors and five as potential oncomiRs. From experimentally validated target genes of these 86 miRNAs, pan-sensitive and pan-resistant genes with concordant mRNA and protein expression associated with in-vitro drug response to 19 NCCN-recommended breast cancer drugs were selected. Combined with in-vitro proliferation assays using CRISPR-Cas9/RNAi and patient survival analysis, MEK inhibitors PD19830 and BRD-K12244279, pilocarpine, and tremorine were discovered as potential new drug options for treating breast cancer. Multi-omics biomarkers of response to the discovered drugs were identified using human breast cancer cell lines. This study presented an artificial intelligence pipeline of miRNA-based discovery of biomarkers, therapeutic targets, and repositioning drugs that can be applied to many cancer types.

MicroRNA, mRNA, and Proteomics Biomarkers and Therapeutic Targets for Improving Lung Cancer Treatment Outcomes

The majority of lung cancer patients are diagnosed with metastatic disease. This study identified a set of 73 microRNAs (miRNAs) that classified lung cancer tumors from normal lung tissues with an overall accuracy of 96.3% in the training patient cohort (n = 109) and 91.7% in unsupervised classification and 92.3% in supervised classification in the validation set (n = 375). Based on association with patient survival (n = 1016), 10 miRNAs were identified as potential tumor suppressors (hsa-miR-144, hsa-miR-195, hsa-miR-223, hsa-miR-30a, hsa-miR-30b, hsa-miR-30d, hsa-miR-335, hsa-miR-363, hsa-miR-451, and hsa-miR-99a), and 4 were identified as potential oncogenes (hsa-miR-21, hsa-miR-31, hsa-miR-411, and hsa-miR-494) in lung cancer. Experimentally confirmed target genes were identified for the 73 diagnostic miRNAs, from which proliferation genes were selected from CRISPR-Cas9/RNA interference (RNAi) screening assays. Pansensitive and panresistant genes to 21 NCCN-recommended drugs with concordant mRNA and protein expression were identified. DGKE and WDR47 were found with significant associations with responses to both systemic therapies and radiotherapy in lung cancer. Based on our identified miRNA-regulated molecular machinery, an inhibitor of PDK1/Akt BX-912, an anthracycline antibiotic daunorubicin, and a multi-targeted protein kinase inhibitor midostaurin were discovered as potential repositioning drugs for treating lung cancer. These findings have implications for improving lung cancer diagnosis, optimizing treatment selection, and discovering new drug options for better patient outcomes.

Inferencing Bulk Tumor and Single-Cell Multi-Omics Regulatory Networks for Discovery of Biomarkers and Therapeutic Targets

There are insufficient accurate biomarkers and effective therapeutic targets in current cancer treatment. Multi-omics regulatory networks in patient bulk tumors and single cells can shed light on molecular disease mechanisms. Integration of multi-omics data with large-scale patient electronic medical records (EMRs) can lead to the discovery of biomarkers and therapeutic targets. In this review, multi-omics data harmonization methods were introduced, and common approaches to molecular network inference were summarized. Our Prediction Logic Boolean Implication Networks (PLBINs) have advantages over other methods in constructing genome-scale multi-omics networks in bulk tumors and single cells in terms of computational efficiency, scalability, and accuracy. Based on the constructed multi-modal regulatory networks, graph theory network centrality metrics can be used in the prioritization of candidates for discovering biomarkers and therapeutic targets. Our approach to integrating multi-omics profiles in a patient cohort with large-scale patient EMRs such as the SEER-Medicare cancer registry combined with extensive external validation can identify potential biomarkers applicable in large patient populations. These methodologies form a conceptually innovative framework to analyze various available information from research laboratories and healthcare systems, accelerating the discovery of biomarkers and therapeutic targets to ultimately improve cancer patient survival outcomes.

Multi-Omics Immune Interaction Networks in Lung Cancer Tumorigenesis, Proliferation, and Survival

There are currently no effective biomarkers for prognosis and optimal treatment selection to improve non-small cell lung cancer (NSCLC) survival outcomes. This study further validated a seven-gene panel for diagnosis and prognosis of NSCLC using RNA sequencing and proteomic profiles of patient tumors. Within the seven-gene panel, ZNF71 expression combined with dendritic cell activities defined NSCLC patient subgroups (n = 966) with distinct survival outcomes (p = 0.04, Kaplan-Meier analysis). ZNF71 expression was significantly associated with the activities of natural killer cells (p = 0.014) and natural killer T cells (p = 0.003) in NSCLC patient tumors (n = 1016) using Chi-squared tests. Overexpression of ZNF71 resulted in decreased expression of multiple components of the intracellular intrinsic and innate immune systems, including dsRNA and dsDNA sensors. Multi-omics networks of ZNF71 and the intracellular intrinsic and innate immune systems were computed as relevant to NSCLC tumorigenesis, proliferation, and survival using patient clinical information and in-vitro CRISPR-Cas9/RNAi screening data. From these networks, pan-sensitive and pan-resistant genes to 21 NCCN-recommended drugs for treating NSCLC were selected. Based on the gene associations with patient survival and in-vitro CRISPR-Cas9, RNAi, and drug screening data, MEK1/2 inhibitors PD-198306 and U-0126, VEGFR inhibitor ZM-306416, and IGF-1R inhibitor PQ-401 were discovered as potential targeted therapy that may also induce an immune response for treating NSCLC.

Hub Genes in Non-Small Cell Lung Cancer Regulatory Networks

There are currently no accurate biomarkers for optimal treatment selection in early-stage non-small cell lung cancer (NSCLC). Novel therapeutic targets are needed to improve NSCLC survival outcomes. This study systematically evaluated the association between genome-scale regulatory network centralities and NSCLC tumorigenesis, proliferation, and survival in early-stage NSCLC patients. Boolean implication networks were used to construct multimodal networks using patient DNA copy number variation, mRNA, and protein expression profiles. T statistics of differential gene/protein expression in tumors versus non-cancerous adjacent tissues, dependency scores in in vitro CRISPR-Cas9/RNA interference (RNAi) screening of human NSCLC cell lines, and hazard ratios in univariate Cox modeling of the Cancer Genome Atlas (TCGA) NSCLC patients were correlated with graph theory centrality metrics. Hub genes in multi-omics networks involving gene/protein expression were associated with oncogenic, proliferative potentials and poor patient survival outcomes (p < 0.05, Pearson's correlation). Immunotherapy targets PD1, PDL1, CTLA4, and CD27 were ranked as top hub genes within the 10th percentile in most constructed multi-omics networks. BUB3, DNM1L, EIF2S1, KPNB1, NMT1, PGAM1, and STRAP were discovered as important hub genes in NSCLC proliferation with oncogenic potential. These results support the importance of hub genes in NSCLC tumorigenesis, proliferation, and prognosis, with implications in prioritizing therapeutic targets to improve patient survival outcomes.

Downregulation of zinc finger protein 71 expression in oral squamous cell carcinoma tissues and its underlying molecular mechanism

BACKGROUND

Patients with oral squamous cell carcinoma (OSCC) have poor prognoses due to a limited understanding of the pathogenesis of OSCC. Zinc finger protein (ZNF) is the largest transcription factor family in the human genome and exert diverse and important functions. Nevertheless, the exact expression status and molecular mechanism of ZNF71 have not been described in OSCC. Therefore, this study aimed to identify the specific expression level of ZNF71 in OSCC tissues and to further interpret the potential molecular mechanism of ZNF71 in the pathogenesis of OSCC.

METHODS

In-house immunohistochemical staining of 116 OSCC samples and 29 non-OSCC samples was employed to detect the expression status of ZNF71 at the protein level of OSCC tissues. Single-cell RNA sequencing data from 7 OSCC samples was used to explore the expression landscape of ZNF71 in different cell types from OSCC tissues. High-throughput RNA sequencing data and gene chips data from 893 OSCC samples and 301 non-OSCC samples were utilized to identify the specific expression level of ZNF71 at the bulk mRNA level of OSCC tissues. Here, standardized mean difference (SMD) value was applied to calculate the expression differences between OSCC group and non-OSCC group. Multiple datasets were included; hence, the results were considered to be more reliable. Sensitivity analysis was conducted to evaluate the stability of the results. Enrichment analysis and immune infiltration analysis were used to explore the underlying molecular mechanism of ZNF71 in OSCC.

RESULTS

ZNF71 was significantly downregulated in OSCC tissues at the protein level (SMD = -1.96, 95 % confidence interval [95 % CI]: -2.43 to -1.50). ZNF71 was absent in various cell types from OSCC tissues including cancerous epithelial cells and tumor-infiltrating immune cells. ZNF71 was downregulated in OSCC tissues at the bulk mRNA level (SMD = -0.38, 95 % CI: -0.75 to -0.02). Enrichment analysis showed that positively and differentially co-expressed genes mainly concentrated on "herpes simplex virus 1 infection" and "regulation of plasma membrane bounded cell projection organization", and negatively and differentially co-expressed genes mainly participated in "cell cycle" and "DNA metabolic process". Moreover, the putative target genes of ZNF71 mainly participated in "cellular respiration" and "protein catabolic process". Finally, immune infiltration analysis revealed that ZNF71 expression was positively correlated with multiple immune cells including activated B cells, memory B cells, and natural killer (NK) cells, and negatively correlated with various immune cells, including CD56 bright NK cells, neutrophil, and immature dendritic cells.

CONCLUSION

The downregulation of ZNF71 may influence the initiation and promotion of OSCC by reducing immune infiltration, accelerating cell cycle progression, and affecting metabolic process, and this requires further research.

Construction of a redox-related gene signature for overall survival prediction and immune infiltration in non-small-cell lung cancer

Background: An imbalance in the redox homeostasis has been reported in multiple cancers and is associated with a poor prognosis of disease. However, the prognostic value of redox-related genes in non-small-cell lung cancer (NSCLC) remains unclear.

Methods: RNA sequencing data, DNA methylation data, mutation, and clinical data of NSCLC patients were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. Redox-related differentially expressed genes (DEGs) were used to construct the prognostic signature using least absolute shrinkage and selection operator (LASSO) regression analysis. Kaplan–Meier survival curve and receiver operator characteristic (ROC) curve analyses were applied to validate the accuracy of the gene signature. Nomogram and calibration plots of the nomogram were constructed to predict prognosis. Pathway analysis was performed using gene set enrichment analysis. The correlations of risk score with tumor stage, immune infiltration, DNA methylation, tumor mutation burden (TMB), and chemotherapy sensitivity were evaluated. The prognostic signature was validated using GSE31210, GSE26939, and GSE68465 datasets. Real-time polymerase chain reaction (PCR) was used to validate dysregulated genes in NSCLC.

Results: A prognostic signature was constructed using the LASSO regression analysis and was represented as a risk score. The high-risk group was significantly correlated with worse overall survival (OS) (p < 0.001). The area under the ROC curve (AUC) at the 5-year stage was 0.657. The risk score was precisely correlated with the tumor stage and was an independent prognostic factor for NSCLC. The constructed nomogram accurately predicted the OS of patients after 1-, 3-, and 5-year periods. DNA replication, cell cycle, and ECM receptor interaction were the main pathways enriched in the high-risk group. In addition, the high-risk score was correlated with higher TMB, lower methylation levels, increased infiltrating macrophages, activated memory CD4⁺ T cells, and a higher sensitivity to chemotherapy. The signature was validated in GSE31210, GSE26939, and GSE68465 datasets. Real-time PCR validated dysregulated mRNA expression levels in NSCLC.

Conclusions: A prognostic redox-related gene signature was successfully established in NSCLC, with potential applications in the clinical setting.

Single B Cell Gene Co-Expression Networks Implicated in Prognosis, Proliferation, and Therapeutic Responses in Non-Small Cell Lung Cancer Bulk Tumors

Simple Summary

This study presents novel insights on dysregulated B cell proliferation networks in non-small cell lung cancer (NSCLC). Within this network, a nine-gene signature demonstrated prognostic and predictive indications in more than 1400 NSCLC patients using their gene and protein expression profiles in bulk tumors. Furthermore, novel therapeutic candidates are identified to improve NSCLC treatment outcomes.

Abstract

In NSCLC, there is a pressing need for immunotherapy predictive biomarkers. The processes underlying B-cell dysfunction, as well as their prognostic importance in NSCLC, are unknown. Tumor-specific B-cell gene co-expression networks were constructed by comparing the Boolean implication modeling of single-cell RNA sequencing of NSCLC tumor B cells and normal B cells. Proliferation genes were selected from the networks using in vitro CRISPR-Cas9/RNA interfering (RNAi) screening data in more than 92 human NSCLC epithelial cell lines. The prognostic and predictive evaluation was performed using public NSCLC transcriptome and proteome profiles. A B cell proliferation and prognostic gene co-expression network was present only in normal lung B cells and missing in NSCLC tumor B cells. A nine-gene signature was identified from this B cell network that provided accurate prognostic stratification using bulk NSCLC tumor transcriptome (n = 1313) and proteome profiles (n = 103). Multiple genes (HLA-DRA, HLA-DRB1, OAS1, and CD74) differentially expressed in NSCLC B cells, peripheral blood lymphocytes, and tumor T cells had concordant prognostic indications at the mRNA and protein expression levels. The selected genes were associated with drug sensitivity/resistance to 10 commonly used NSCLC therapeutic regimens. Lestaurtinib was discovered as a potential repositioning drug for treating NSCLC.