Low expression of moonlight gene ALAD is correlated with poor prognosis in hepatocellular carcinoma

Construction of a prognostic model for breast cancer based on moonlighting genes

Breast cancer (BRCA) is a highly heterogeneous disease, with significant differences in prognosis among patients. Existing biomarkers and prognostic models have limited ability to predict BRCA prognosis. Moonlighting genes regulate tumor progression and are associated with cancer prognosis. This study aimed to construct a moonlighting gene-based prognostic model for BRCA. We obtained differentially expressed genes (DEGs) in BRCA from The Cancer Genome Atlas and intersected them with moonlighting genes from MoonProt to acquire differential moonlighting genes. GO and KEGG results showed main enrichment of these genes in the response of BRCA cells to environmental stimuli and pentose phosphate pathway. Based on moonlighting genes, we conducted drug prediction and validated results through cellular experiments. After ABCB1 knockdown, viability and proliferation of BRCA cells were significantly enhanced. Based on differential moonlighting genes, BRCA was divided into three subgroups, among which cluster2 had the highest survival rate and immunophenoscore and relatively low tumor mutation burden. TP53 had the highest mutation frequency in cluster2 and cluster3, while PIK3CA had a higher mutation frequency in cluster1, with the majority being missense mutations. Subsequently, we established an 11-gene prognostic model in the training set based on DEGs among subgroups using univariate Cox regression, LASSO regression, and multivariable Cox regression analyses. Model prognostic performance was verified in GEO, METABRIC and ICGC validation sets. In summary, this study obtained three BRCA moonlighting gene-related subtypes and constructed an 11-gene prognostic model. The 11-gene BRCA prognostic model has good predictive performance, guiding BRCA prognosis for clinical doctors.

Systematic proteome-wide Mendelian randomization using the human plasma proteome to identify therapeutic targets for lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the predominant histological subtype of lung cancer and the leading cause of cancer-related mortality. Identifying effective drug targets is crucial for advancing LUAD treatment strategies.

METHODS

This study employed proteome-wide Mendelian randomization (MR) and colocalization analyses. We collected data on 1394 plasma proteins from a protein quantitative trait loci (pQTL) study involving 4907 individuals. Genetic associations with LUAD were derived from the Transdisciplinary Research in Cancer of the Lung (TRICL) study, including 11,245 cases and 54,619 controls. We integrated pQTL and LUAD genome-wide association studies (GWASs) data to identify candidate proteins. MR utilizes single nucleotide polymorphisms (SNPs) as genetic instruments to estimate the causal effect of exposure on outcome, while Bayesian colocalization analysis determines the probability of shared causal genetic variants between traits. Our study applied these methods to assess causality between plasma proteins and LUAD. Furthermore, we employed a two-step MR to quantify the proportion of risk factors mediated by proteins on LUAD. Finally, protein-protein interaction (PPI) analysis elucidated potential links between proteins and current LUAD medications.

RESULTS

We identified nine plasma proteins significantly associated with LUAD. Increased levels of ALAD, FLT1, ICAM5, and VWC2 exhibited protective effects, with odds ratios of 0.79 (95% CI 0.72-0.87), 0.39 (95% CI 0.28-0.55), 0.91 (95% CI 0.72-0.87), and 0.85 (95% CI 0.79-0.92), respectively. Conversely, MDGA2 (OR, 1.13; 95% CI 1.08-1.19), NTM (OR, 1.12; 95% CI 1.09-1.16), PMM2 (OR, 1.35; 95% CI 1.18-1.53), RNASET2 (OR, 1.15; 95% CI 1.08-1.21), and TFPI (OR, 4.58; 95% CI 3.02-6.94) increased LUAD risk. Notably, none of the nine proteins showed evidence of reverse causality. Bayesian colocalization indicated that RNASET2, TFPI, and VWC2 shared the same variant with LUAD. Furthermore, NTM and FLT1 demonstrated interactions with targets of current LUAD medications. Additionally, FLT1 and TFPI are currently under evaluation as therapeutic targets, while NTM, RNASET2, and VWC2 are potentially druggable. These findings shed light on LUAD pathogenesis, highlighting the tumor-promoting effects of RNASET2, TFPI, and NTM, along with the protective effects of VWC2 and FLT1, providing a significant biological foundation for future LUAD therapeutic targets.

CONCLUSIONS

Our proteome-wide MR analysis highlighted RNASET2, TFPI, VWC2, NTM, and FLT1 as potential drug targets for further clinical investigation in LUAD. However, the specific mechanisms by which these proteins influence LUAD remain elusive. Targeting these proteins in drug development holds the potential for successful clinical trials, providing a pathway to prioritize and reduce costs in LUAD therapeutics.

Novel cell subtypes of SPP1 + S100P+, MS4A1-SPP1 + S100P+ were key subpopulations in intrahepatic cholangiocarcinoma

BACKGROUND

In this study, we integrated single-cell RNA sequencing (scRNA-seq) data to investigate cell heterogeneity and utilized MSigDB and CIBERSORTx to explore the pathways of major cell types and the relationships between different cell subtypes. Subsequently, we explored the correlation of cell subtypes with survival and used Gene Set Enrichment Analysis (GSEA) analyses to assess the pathways associated with the infiltration of specific cell subtypes. Finally, multiplex immunohistochemistry in tissue microarray cohort were performed to validate differences in protein level and their correlation with survival.

RESULTS

iCCA presented a unique immune ecosystem, with increased proportions of Epi (epithelial)-SPP1-2, Epi-S100P-1, Epi-DN (double negative for SPP1 and S100P expression)-1, Epi-DN-2, Epi-DP (double positive for SPP1 and S100P expression)-1, Plasma B-3, Plasma B-2, B-HSPA1A-1, B-HSPA1A-2 cells, and decreased proportions of B-MS4A1. High level of Epi-DN-2, Epi-SPP1-1, Epi-SPP1-2, B-MS4A1, and low level of Epi-DB-1, Epi-S100P-1, and Epi-S100P-2 was significantly associated with longer overall survival (OS), and high level of B-MS4A1_Low_Epi-DN-2_Low was associated with the shortest OS. Moreover, the results of MsigDB and GSEA suggest that bile acid metabolism is a crucial process in iCCA. Finally, we found that S100P+, SPP1+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ were highly expressed, whereas MS4A1 was lowly expressed in iCCA, and patients with high level of S100P+, SPP1 + S100P+, and MS4A1-SPP1 + S100P+ exhibited shorter survival.

CONCLUSIONS

We identified the cell heterogeneity of iCCA, found that iCCA is a unique immune ecosystem with many cell subtypes, and showed that the novel cell subtypes of SPP1 + S100P+ and MS4A1-SPP1 + S100P+ were key subpopulations in iCCA.

Network analysis of hepatocellular carcinoma liquid biopsies augmented by single-cell sequencing data

Liquid biopsy, the analysis of body fluids, represents a promising approach for disease diagnosis and prognosis with minimal intervention. Sequencing cell-free RNA derived from liquid biopsies has been very promising for the diagnosis of several diseases. Cancer research, in particular, has emerged as a prominent candidate since early diagnosis has been shown to be a critical determinant of disease prognosis. Although high-throughput analysis of liquid biopsies has uncovered many differentially expressed genes in the context of cancer, the functional connection between these genes is not investigated in depth. An important approach to remedy this issue is the construction of gene networks which describes the correlation patterns between different genes, thereby allowing to infer their functional organization. In this study, we aimed at characterizing extracellular transcriptome gene networks of hepatocellular carcinoma patients compared to healthy controls. Our analysis revealed a number of genes previously associated with hepatocellular carcinoma and uncovered their association network in the blood. Our study thus demonstrates the feasibility of performing gene co-expression network analysis from cell-free RNA data and its utility in studying hepatocellular carcinoma. Furthermore, we augmented cell-free RNA network analysis with single-cell RNA sequencing data which enables the contextualization of the identified network modules with cell-type specific transcriptomes from the liver.