Visualizing Single-Cell RNA-seq Data with Semisupervised Principal Component Analysis

Identifying collective variables (CVs) for chemical reactions is essential to reduce the 3N-dimensional energy landscape into lower dimensional basins and barriers of interest. However, in condensed phase processes, the nonmeaningful motions of bulk solvent often overpower the ability of dimensionality reduction methods to identify correlated motions that underpin collective variables. Yet solvent can play important indirect or direct roles in reactivity, and much can be lost through treatments that remove or dampen solvent motion. This has been amply demonstrated within principal component analysis (PCA), although less is known about the behavior of nonlinear dimensionality reduction methods, e.g., uniform manifold approximation and projection (UMAP), that have become recently utilized. The latter presents an interesting alternative to linear methods though often at the expense of interpretability. This work presents distance-attenuated projection methods of atomic coordinates that facilitate the application of both PCA and UMAP to identify collective variables in the presence of explicit solvent and further the specific identity of solvent molecules that participate in chemical reactions. The performance of both methods is examined in detail for two reactions where the explicit solvent plays very different roles within the collective variables. When applied to raw molecular dynamics data in solution, both PCA and UMAP representations are dominated by bulk solvent motions. On the other hand, when applied to data preprocessed by our attenuated projection methods, both PCA and UMAP identify the appropriate collective variables (though varying sensitivity is observed due to the presence of explicit solvent that results from the projection method). Importantly, this approach allows identification of specific solvent molecules that are relevant to the CVs and their importance.

BackgroundThis study aimed to use single-cell RNA-seq (scRNA-seq) to discover marker genes in endothelial cells (ECs) and construct a prognostic model for glioblastoma multiforme (GBM) patients in combination with traditional high-throughput RNA sequencing (bulk RNA-seq).MethodsBulk RNA-seq data was downloaded from The Cancer Genome Atlas (TCGA) and The China Glioma Genome Atlas (CGGA) databases. 10x scRNA-seq data for GBM were obtained from the Gene Expression Omnibus (GEO) database. The uniform manifold approximation and projection (UMAP) were used for downscaling and cluster identification. Key modules and differentially expressed genes (DEGs) were identified by weighted gene correlation network analysis (WGCNA). A non-negative matrix decomposition (NMF) algorithm was used to identify the different subtypes based on DEGs, and multivariate cox regression analysis to model the prognosis. Finally, differences in mutational landscape, immune cell abundance, immune checkpoint inhibitors (ICIs)-associated genes, immunotherapy effects, and enriched pathways were investigated between different risk groups.ResultsThe analysis of scRNA-seq data from eight samples revealed 13 clusters and four cell types. After applying Fisher’s exact test, ECs were identified as the most important cell type. The NMF algorithm identified two clusters with different prognostic and immunological features based on DEGs. We finally built a prognostic model based on the expression levels of four key genes. Higher risk scores were significantly associated with poorer survival outcomes, low mutation rates in IDH genes, and upregulation of immune checkpoints such as PD-L1 and CD276.ConclusionWe built and validated a 4-gene signature for GBM using 10 scRNA-seq and bulk RNA-seq data in this work.

Numerical modelling increasingly generates massive, high-dimensional spatio-temporal datasets. Exploring such datasets relies on effective visualization. This study presents a generic workflow to (i) project high-dimensional spatio-temporal data on a two-dimensional (2D) plane accurately (ii) compare dimensionality reduction techniques (DRTs) in terms of resolution and computational efficiency (iii) represent 2D projection spatially using a 2D perceptually uniform background color map. Machine learning (ML) based DRTs for data visualization i.e., principal component analysis (PCA), generative topographic mapping (GTM), t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP) are compared in terms of accuracy, resolution and computational efficiency to handle massive datasets. The accuracy of visualization is evaluated using a quality metric based on a co-ranking framework. The workflow is applied to an output of an Australian Water Resource Assessment (AWRA) model for Tasmania, Australia. The dataset consists of daily time series of nine components of the water balance at a 5 km grid cell resolution for the year 2017. The case study shows that PCA allows rapid visualization of global data structures, while t-SNE and UMAP allows more accurate representation of local trends. Furthermore, UMAP is computationally more efficient than t-SNE and least affected by the outliers compared to GTM.

We aimed to establish a rapid and precise method for identifying and quantifying multi-species vegetable oil (corn oil, olive oil (OLO), soybean oil, and sunflower oil (SUO)) adulterations in camellia oil (CAO), using soft ionization by chemical reaction in transfer–high-resolution mass spectrometry (SICRIT-HRMS) and machine learning methods. The results showed that SICRIT-HRMS could effectively characterize the volatile profiles of pure and adulterated CAO samples, including binary, ternary, quaternary, and quinary adulteration systems. The low m/z region (especially 100–300) exhibited importance to oil classification in multiple feature-selection methods. For qualitative detection, binary classification models based on convolutional neural networks (CNN), Random Forest (RF), and gradient boosting trees (GBT) algorithms showed high accuracies (98.70–100.00%) for identifying CAO adulteration under no dimensionality reduction (NON), principal component analysis (PCA), and uniform manifold approximation and projection (UMAP) strategies. The RF algorithm exhibited relatively high accuracy (96.25–99.45%) in multiclass classification. Moreover, the five models, including CNN, RF, support vector machines (SVM), logistic regression (LR), and GBT, exhibited different performances in distinguishing pure and adulterated CAO. Among 1093 blind oil samples, under NON, PCA, and UMAP: 10, 5, and 67 samples were misclassified by CNN model; 6, 7, and 41 samples were misclassified by RF model; 8, 9, and 82 samples were misclassified by SVM model; 17, 18, and 78 samples were misclassified by LR model; 7, 9, and 43 samples were misclassified by GBT model. For quantitative prediction, the PCA-CNN model performed optimally in predicting adulteration levels in CAO, especially with respect to OLO and SUO, exhibiting a high coefficient of determination for calibration (, 0.9664–0.9974) and coefficient of determination for prediction (, 0.9599–0.9963) values, low root mean square error of calibration (RMSEC, 0.9–5.3%) and root mean square error of prediction (RMSEP, 1.1–5.8%) values, and RPD (5.0–16.3) values greater than 3.0. These results indicate that SICRIT-HRMS combined with machine learning can rapidly and accurately identify and quantify multi-species vegetable oil adulterations in CAO, which provides a reference for developing non-targeted and high-throughput detection methods in edible oil authenticity.

This study used an odor sensing system with a 16-channel electrochemical sensor array to measure beef odors, aiming to distinguish odors under different storage days and processing temperatures for quality monitoring. Six storage days ranged from purchase (D0) to eight days (D8), with three temperature conditions: no heat (RT), boiling (100 °C), and frying (180 °C). Gas chromatography–mass spectrometry (GC-MS) analysis showed that odorants in the beef varied under different conditions. Compounds like acetoin and 1-hexanol changed significantly with the storage days, while pyrazines and furans were more detectable at higher temperatures. The odor sensing system data were visualized using principal component analysis (PCA) and uniform manifold approximation and projection (UMAP). PCA and unsupervised UMAP clustered beef odors by storage days but struggled with the processing temperatures. Supervised UMAP accurately clustered different temperatures and dates. Machine learning analysis using six classifiers, including support vector machine, achieved 57% accuracy for PCA-reduced data, while unsupervised UMAP reached 49.1% accuracy. Supervised UMAP significantly enhanced the classification accuracy, achieving over 99.5% with the dimensionality reduced to three or above. Results suggest that the odor sensing system can sufficiently enhance non-destructive beef quality and safety monitoring. This research advances electronic nose applications and explores data downscaling techniques, providing valuable insights for future studies.

Cluster-Based Insights into Cardiovascular and Autonomic Responses to Head-Up Tilt in Hypertension.

Crustacean aquaculture is expected to be a major source of fishery commodities in the near future. Hemocytes are key players of the immune system in shrimps; however, their classification, maturation, and differentiation are still under debate. To date, only discrete and inconsistent information on the classification of shrimp hemocytes has been reported, showing that the morphological characteristics are not sufficient to resolve their actual roles. Our present study using single-cell RNA sequencing revealed six types of hemocytes of Marsupenaeus japonicus based on their transcriptional profiles. We identified markers of each subpopulation and predicted the differentiation pathways involved in their maturation. We also predicted cell growth factors that might play crucial roles in hemocyte differentiation. Different immune roles among these subpopulations were suggested from the analysis of differentially expressed immune-related genes. These results provide a unified classification of shrimp hemocytes, which improves the understanding of its immune system.

In this paper, we study phase transitions of the q-state Potts model through a number of unsupervised machine learning techniques, namely Principal Component Analysis (PCA), k-means clustering, Uniform Manifold Approximation and Projection (UMAP), and Topological Data Analysis (TDA). Even though in all cases we are able to retrieve the correct critical temperatures $$T_\textrm{c}(q)$$ , for $$q=3,4$$ and 5, results show that non-linear methods as UMAP and TDA are less dependent on finite-size effects. This study may be considered as a benchmark for the use of different unsupervised machine learning algorithms in the investigation of phase transitions.

Comparison of dimension reduction techniques applied to the analysis of airborne radionuclide activity concentration

Advances in single-cell technologies have enabled high-resolution dissection of tissue composition. Several tools for dimensionality reduction are available to analyze the large number of parameters generated in single-cell studies. Recently, a nonlinear dimensionality-reduction technique, uniform manifold approximation and projection (UMAP), was developed for the analysis of any type of high-dimensional data. Here we apply it to biological data, using three well-characterized mass cytometry and single-cell RNA sequencing datasets. Comparing the performance of UMAP with five other tools, we find that UMAP provides the fastest run times, highest reproducibility and the most meaningful organization of cell clusters. The work highlights the use of UMAP for improved visualization and interpretation of single-cell data.

USP39 is involved in mRNA splicing and stress responses. This study analyzed USP39 expression and its clinical relevance in sepsis, aiming to provide a novel treatment target for sepsis. The RNA-seq datasets were used for analyzing the expression profile and clinical significance of USP39 in both sepsis and control groups. Functional pathway enrichment was performed using GSEA and GSVA with the clusterProfiler package (3.14.3). Immune infiltration was evaluated via single-sample GSEA (ssGSEA). The prognostic value of USP39 was assessed utilizing Cox regression, and its diagnostic performance was determined via receiver operating characteristic (ROC) curve analysis employing the survival (3.5.7) and pROC (1.18.5) packages, respectively. For single-cell RNA-seq (scRNA-seq) data (GSE175453), we performed quality control, normalization, principal component analysis (PCA), batch correction, clustering, and Uniform Manifold Approximation and Projection (UMAP) visualization employing the Seurat 4.4.0 package. THP-1 cells were treated with lipopolysaccharide (LPS) to mimic septic injury. USP39 was overexpressed via pcDNA3.1 transfection. Cellular assays were conducted to measure cell viability, apoptosis, and inflammatory cytokines (IL-1β, IL-6, TNF-α). Our study found that USP39 was significantly downregulated in sepsis patients and was associated with age, diabetes, and Intensive Care Unit (ICU)-acquired infections. Functional enrichment analysis revealed that high USP39 expression was linked to metabolic processes, primary immunodeficiency, and spliceosome, and immune response and inflammation pathway. In addition, USP39 also exhibited a high diagnostic value and was identified as an independent prognostic marker in sepsis. ssGSEA revealed higher infiltration of CD8+ T cells, activated B cells, and CD4+ T cells in the high USP39 expression group, while Th17 and NK cells were more abundant in the low USP39 expression group. According to the results of scRNA-seq analysis, USP39 exhibited significant differential expressions in monocytes and T cells, with low expression observed in sepsis. Furthermore, in vitro experiments confirmed that overexpression of USP39 significantly enhanced cell viability and suppressed apoptosis in LPS-induced THP-1 cells while also inhibiting the expression of inflammatory cytokines. The significant downregulation of USP39 was correlated with both clinical features and immune infiltration in sepsis, highlighting its role in immune regulation and its potential as a prognostic biomarker.

Ischemic stroke (IS) is a leading cause of adult disability worldwide. The inflammatory processes involved are complex, making it challenging to fully understand the pathological mechanisms of IS. Phagocytosis plays an important role in eliminating neurotoxic or damaged neurons resulting from inflammatory responses. This study employed bioinformatics methods to analyze single-cell RNA sequencing (scRNA-seq) data to investigate the cell types and molecular biological processes involved in IS. scRNA-seq data for IS were obtained from the Gene Expression Omnibus (GEO). Following sample screening and reprocessing, 5,582 single cells were identified from healthy controls and patients with IS. Uniform manifold approximation and projection (UMAP) was utilized to further explore the cellular composition in IS. Functional enrichment analysis of differentially expressed genes was conducted to identify transcriptional regulators, whereas cell developmental trajectories were predicted to uncover potential cell fate decisions. iTALK was employed to identify potential ligand-receptor axes within the cell-type immune microenvironment of IS. Based on scRNA-seq data analysis, we identified four cell types and their associated subclusters, along with genes exhibiting significant differential expression within these subclusters. Phagocytosis was significantly enriched in cell types linked to IS, while the differentiation trajectories of subpopulations in IS was different. Additionally, multiple receptor-ligand axes were identified, indicating diverse interactions within the immune microenvironment of IS. This study demonstrated that phagocytosis in IS cell types critically influences disease progression. It also predicted the trajectories of infarct cells. These findings provide valuable insights into the molecular and cellular mechanisms underlying IS and highlight potential pathways for therapeutic intervention.

This study compares the effectiveness of various dimensionality reduction (DR) techniques-including Principal Component Analysis (PCA), Uniform Manifold Approximation and Projection (UMAP), Multidimensional Scaling (MDS), and t-Distributed Stochastic Neighbor Embedding (t-SNE)-combined with a Support Vector Machine (SVM) classifier for a diagnostic task. The results showed that applying DR significantly enhanced performance over the baseline SVM (Accuracy: 0.68). Specifically, the fusion of t-SNE and an Autoencoder achieved the highest accuracy at 94 %, outperforming all other configurations. This demonstrates that leveraging advanced non-linear DR and feature learning methods provides a superior and highly accurate framework for classification.

Dimensionality reduction is often used to visualize complex expression profiling data. Here, we use the Uniform Manifold Approximation and Projection (UMAP) method on published transcript profiles of 1484 single gene deletions of Saccharomyces cerevisiae. Proximity in low-dimensional UMAP space identifies groups of genes that correspond to protein complexes and pathways, and finds novel protein interactions, even within well-characterized complexes. This approach is more sensitive than previous methods and should be broadly useful as additional transcriptome datasets become available for other organisms.

Transcriptomic Signature and Functional Abnormalities of Bone Marrow Mesenchymal Stromal Cells Mediate Disease Progression of Myelodysplastic/Myeloproliferative Neoplasms