Downregulation of CYB5D2 is associated with breast cancer progression

Identification of fatty acid metabolism-related clusters and immune infiltration features in hepatocellular carcinoma

Hepatocellular Carcinoma (HCC) is a type of liver cancer which is characterized by inflammation-associated tumor. The unique characteristics of tumor immune microenvironment in HCC contribute to hepatocarcinogenesis. It was also clarified that aberrant fatty acid metabolism (FAM) might accelerate tumor growth and metastasis of HCC. In this study, we aimed to identify fatty acid metabolism-related clusters and establish a novel prognostic risk model in HCC. Gene expression and corresponding clinical data were searched from the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) portal. From the TCGA database, by unsupervised clustering method, we determined three FAM clusters and two gene clusters with distinct clinicopathological and immune characteristics. Based on 79 prognostic genes identified from 190 differentially expressed genes (DEGs) among three FAM clusters, five prognostic DEGs (CCDC112, TRNP1, CFL1, CYB5D2, and SLC22A1) were determined to construct risk model by least absolute shrinkage and selection operator (LASSO) and multivariate cox regression analysis. Furthermore, the ICGC dataset was used to validate the model. In conclusion, the prognostic risk model constructed in this study exhibited excellent indicator performance of overall survival, clinical feature, and immune cell infiltration, which has the potential to be an effective biomarker for HCC immunotherapy.

Many or too many progesterone membrane receptors? Clinical implications

Several receptors for nongenomically initiated actions of progesterone (P4) exist, namely membrane-associated P4 receptors (MAPRs), membrane progestin receptors (mPRs), receptors for neurosteroids [GABA_A receptor (GABA_AR), NMDA receptor, sigma-1 and -2 receptors (S1R/S2R)], the classical genomic P4 receptor (PGR), and α/β hydrolase domain-containing protein 2 (ABHD2). Two drugs related to this field have been approved: brexanolone (Zulresso™) for the treatment of postpartum depression, and ganaxolone (Ztalmy™) for the treatment of CDKL5 deficiency disorder. Both are derivatives of P4 and target the GABA_AR. Several other indications are in clinical testing. CT1812 (Elayta™) is also being tested for the treatment of Alzheimer's disease (AD) in Phase 2 clinical trials, targeting the P4 receptor membrane component 1 (PGRMC1)/S2R complex. In this Review, we highlight emerging knowledge on the mechanisms of nongenomically initiated actions of P4 and its derivatives.

Integration of Bulk RNA Sequencing and Single-Cell RNA Sequencing to Reveal Uveal Melanoma Tumor Heterogeneity and Cells Related to Survival

Molecular classification based on transcriptional characteristics is often used to study tumor heterogeneity. Human cancer has different cell populations with distinct transcription in tumors, and their heterogeneity is the focus of tumor therapy. Our purpose was to explore the tumor heterogeneity of uveal melanoma (UM) through RNA sequencing (RNA-seq) and single-cell RNA sequencing (scRNA-seq). Based on the consensus clustering assays of the prognosis-related immune gene set, the immune subtype (IS) of UM and its corresponding immune characteristics were comprehensively analyzed. The heterogeneous cell groups and corresponding marker genes of UM were identified from GSE138433 using scRNA-seq analysis. Pseudotime trajectory analysis and SCENIC analysis were conducted to explore the trajectory of cell differentiation and the regulatory network of single-cell transcription factors (TFs). Based on 37 immune gene sets, UM was divided into two different immune subtypes (IS1 and IS2). The two kinds of ISs have different characteristics in prognosis, immune-related molecules, immune score, and immune cell infiltration. According to 11,988 cells of scRNA-seq data from six UM samples, 11 cell clusters and 10 cell types were identified. The subsets of C1, C4, C5, C8, and C9 were related to the prognosis of UM, and different TF–target gene regulatory networks were involved. These five cell subsets differentiated into 3 different states. Our results provided valuable information about the heterogeneity of UM tumors and the expression patterns of TFs in different cell types.

Construction and evaluation of a nomogram for predicting survival in patients with lung cancer

BACKGROUND

Lung cancer is a heterogeneous disease with a severe disease burden. Because the prognosis of patients with lung cancer varies, it is critical to identify effective biomarkers for prognosis prediction.

METHODS

A total of 2325 lung cancer patients were integrated into four independent sets (training set, validation set I, II and III) after removing batch effects in our study. We applied the microarray data algorithm to screen the differentially expressed genes in the training set. The most robust markers for prognosis were identified using the LASSO-Cox regression model, which was then used to create a Cox model and nomogram.

RESULTS

Through LASSO and multivariate Cox regression analysis, eight genes were identified as prognosis-associated hub genes, followed by the creation of prognosis-associated risk scores (PRS). The results of the Kaplan-Meier analysis in the three validation sets demonstrate the good predictive performance of PRS, with hazard ratios of 2.38 (95% confidence interval (CI), 1.61-3.53) in the validation set I, 1.35 (95% CI, 1.06-1.71) in the validation set II, and 2.71 (95% CI, 1.77-4.18) in the validation set III. Additionally, the PRS demonstrated superior survival prediction in subgroups by age, gender, p-stage, and histologic type (p < 0.0001). The complex model integrating PRS and clinical risk factors also have a good predictive performance for 3-year overall survival.

CONCLUSIONS

In this study, we developed a PRS signature to help predict the survival of lung cancer. By combining it with clinical risk factors, a nomogram was established to quantify the individual risk assessments.

Integrative Analysis of Epigenome and Transcriptome Data Reveals Aberrantly Methylated Promoters and Enhancers in Hepatocellular Carcinoma

DNA methylation is a key transcription regulator, whose aberration was ubiquitous and important in most cancers including hepatocellular carcinoma (HCC). Whole-genome bisulfite sequencing (WGBS) was conducted for comparison of DNA methylation in tumor and adjacent tissues from 33 HCC patients, accompanying RNA-seq to determine differentially methylated region-associated, differentially expressed genes (DMR-DEGs), which were independently replicated in the TCGA-LIHC cohort and experimentally validated via 5-aza-2-deoxycytidine (5-azadC) demethylation. A total of 9,867,700 CpG sites showed significantly differential methylation in HCC. Integrations of mRNA-seq, histone ChIP-seq, and WGBS data identified 611 high-confidence DMR-DEGs. Enrichment analysis demonstrated activation of multiple molecular pathways related to cell cycle and DNA repair, accompanying repression of several critical metabolism pathways such as tyrosine and monocarboxylic acid metabolism. In TCGA-LIHC, we replicated about 53% of identified DMR-DEGs and highlighted the prognostic significance of combinations of methylation and expression of nine DMR-DEGs, which were more efficient prognostic biomarkers than considering either type of data alone. Finally, we validated 22/23 (95.7%) DMR-DEGs in 5-azadC-treated LO2 and/or HepG2 cells. In conclusion, integration of epigenome and transcriptome data depicted activation of multiple pivotal cell cycle-related pathways and repression of several metabolic pathways triggered by aberrant DNA methylation of promoters and enhancers in HCC.

Analysis of Stop Codons within Prokaryotic Protein-Coding Genes Suggests Frequent Readthrough Events

Nonsense mutations turn a coding (sense) codon into an in-frame stop codon that is assumed to result in a truncated protein product. Thus, nonsense substitutions are the hallmark of pseudogenes and are used to identify them. Here we show that in-frame stop codons within bacterial protein-coding genes are widespread. Their evolutionary conservation suggests that many of them are not pseudogenes, since they maintain dN/dS values (ratios of substitution rates at non-synonymous and synonymous sites) significantly lower than 1 (this is a signature of purifying selection in protein-coding regions). We also found that double substitutions in codons—where an intermediate step is a nonsense substitution—show a higher rate of evolution compared to null models, indicating that a stop codon was introduced and then changed back to sense via positive selection. This further supports the notion that nonsense substitutions in bacteria are relatively common and do not necessarily cause pseudogenization. In-frame stop codons may be an important mechanism of regulation: Such codons are likely to cause a substantial decrease of protein expression levels.

SLFN5 influences proliferation and apoptosis by upregulating PTEN transcription via ZEB1 and inhibits the purine metabolic pathway in breast cancer

Human Schlafen-5 (SLFN5) is aberrantly involved in tumorigenesis in several types of cancer. However, its implications in breast cancer (BRCA) are unknown. Herein, we demonstrated that SLFN5 expression is negatively associated with the tumour growth of human BRCA using GEO database analysis and clinical sample immunostaining. Lentiviral overexpression of SLFN5 in BRCA cell lines suppressed tumourigenicity in nude mice. Knockdown and overexpression of SLFN5 in BRCA cell lines proved that SLFN5 can inhibit cell proliferation and colony formation and promote apoptosis by upregulating the transcription of a known cancer suppressor gene (the phosphatase and tensin homologue on chromosome 10, PTEN), resulting in molecular changes in the downstream AKT pathway and in proliferation/apoptosis. Lentiviral knockdown and overexpression of ZEB1 blocked the changes in the PTEN and AKT pathways and in the colony formation ability caused by SLFN5 knockdown and overexpression, respectively. Luciferase reporter assays demonstrated that ZEB1 can inhibit the PTEN promoter activity in MCF7 cells by binding to a motif in the PTEN promoter. Metabonomics analysis showed that SLFN5 influences many metabolic pathways and especially decreases purine metabolites, including inosine, xanthine, and hypoxanthine. In conclusion, our findings suggest that SLFN5 may be an important protective factor against BRCA, as it regulates PTEN transcription, the AKT pathway, and proliferation/apoptosis via ZEB1 mediation and inhibits the purine metabolic pathway. Thus, SLFN5 may be a potential therapeutic target for BRCA.

Construction and Validation of a 13-Gene Signature for Prognosis Prediction in Medulloblastoma

Background: Recent studies have identified several molecular subgroups of medulloblastoma associated with distinct clinical outcomes; however, no robust gene signature has been established for prognosis prediction. Our objective was to construct a robust gene signature-based model to predict the prognosis of patients with medulloblastoma.

Methods: Expression data of medulloblastomas were acquired from the Gene Expression Omnibus (GSE85217, n = 763; GSE37418, n = 76). To identify genes associated with overall survival (OS), we performed univariate survival analysis and least absolute shrinkage and selection operator (LASSO) Cox regression. A risk score model was constructed based on selected genes and was validated using multiple datasets. Differentially expressed genes (DEGs) between the risk groups were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and protein–protein interaction (PPI) analyses were performed. Network modules and hub genes were identified using Cytoscape. Furthermore, tumor microenvironment (TME) was evaluated using ESTIMATE algorithm. Tumor-infiltrating immune cells (TIICs) were inferred using CIBERSORTx.

Results: A 13-gene model was constructed and validated. Patients classified as high-risk group had significantly worse OS than those as low-risk group (Training set: p < 0.0001; Validation set 1: p < 0.0001; Validation set 2: p = 0.00052). The area under the curve (AUC) of the receiver operating characteristic (ROC) analysis indicated a good performance in predicting 1-, 3-, and 5-year OS in all datasets. Multivariate analysis integrating clinical factors demonstrated that the risk score was an independent predictor for the OS (validation set 1: p = 0.001, validation set 2: p = 0.004). We then identified 265 DEGs between risk groups and PPI analysis predicted modules that were highly related to central nervous system and embryonic development. The risk score was significantly correlated with programmed death-ligand 1 (PD-L1) expression (p < 0.001), as well as immune score (p = 0.035), stromal score (p = 0.010), and tumor purity (p = 0.010) in Group 4 medulloblastomas. Correlations between the 13-gene signature and the TIICs in Sonic hedgehog and Group 4 medulloblastomas were revealed.

Conclusion: Our study constructed and validated a robust 13-gene signature model estimating the prognosis of medulloblastoma patients. We also revealed genes and pathways that may be related to the development and prognosis of medulloblastoma, which might provide candidate targets for future investigation.