Search

Background Single-cell RNA sequencing (scRNA-seq) methods have been advantageous for quantifying cell-to-cell variation by profiling the transcriptomes of individual cells. For scRNA-seq data, variability in gene expression reflects the degree of variation in gene expression from one cell to another. Analyses that focus on cell–cell variability therefore are useful for going beyond changes based on average expression and, instead, identifying genes with homogeneous expression versus those that vary widely from cell to cell. Results We present a novel statistical framework, scShapes , for identifying differential distributions in single-cell RNA-sequencing data using generalized linear models. Most approaches for differential gene expression detect shifts in the mean value. However, as single-cell data are driven by overdispersion and dropouts, moving beyond means and using distributions that can handle excess zeros is critical. scShapes quantifies gene-specific cell-to-cell variability by testing for differences in the expression distribution while flexibly adjusting for covariates if required. We demonstrate that scShapes identifies subtle variations that are independent of altered mean expression and detects biologically relevant genes that were not discovered through standard approaches. Conclusions This analysis also draws attention to genes that switch distribution shapes from a unimodal distribution to a zero-inflated distribution and raises open questions about the plausible biological mechanisms that may give rise to this, such as transcriptional bursting. Overall, the results from scShapes help to expand our understanding of the role that gene expression plays in the transcriptional regulation of a specific perturbation or cellular phenotype. Our framework scShapes is incorporated into a Bioconductor R package (https://www.bioconductor.org/packages/release/bioc/html/scShapes.html). [ABSTRACT FROM AUTHOR]

In search of redox mechanisms in breast cancer, we uncovered a striking role for glutathione peroxidase 2 (GPx2) in oncogenic signaling and patient survival. GPx2 loss stimulates malignant progression due to reactive oxygen species/hypoxia inducible factor-α (HIF1α)/ VEGFA (vascular endothelial growth factor A) signaling, causing poor perfusion and hypoxia, which were reversed by GPx2 reexpression or HIF1α inhibition. Ingenuity Pathway Analysis revealed a link between GPx2 loss, tumor angiogenesis, metabolic modulation, and HIF1α signaling. Single-cell RNA analysis and bioenergetic profiling revealed that GPx2 loss stimulated the Warburg effect in most tumor cell subpopulations, except for one cluster, which was capable of oxidative phosphorylation and glycolysis, as confirmed by coexpression of phosphorylated-AMPK and GLUT1. These findings underscore a unique role for redox signaling by GPx2 dysregulation in breast cancer, underlying tumor heterogeneity, leading to metabolic plasticity and malignant progression. [ABSTRACT FROM AUTHOR]