Visual analysis of mass cytometry data by hierarchical stochastic neighbour embedding reveals rare cell types.
Ontology highlight
ABSTRACT:
SUBMITTER: van Unen V
PROVIDER: S-EPMC5700955 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Technological advances in mass spectrometry imaging (MSI) have contributed to growing interest in 3D MSI. However, the large size of 3D MSI data sets has made their efficient analysis and visualization and the identification of informative molecular patterns computationally challenging. Hierarchical stochastic neighbor embedding (HSNE), a nonlinear dimensionality reduction technique that aims at finding hierarchical and multiscale representations of large data sets, is a recent development that enables the analysis of millions of data points, with manageable time and memory complexities. We demonstrate that HSNE can be used to analyze large 3D MSI data sets at full mass spectral and spatial resolution. To benchmark the technique as well as demonstrate its broad applicability, we have analyzed a number of publicly available 3D MSI data sets, recorded from various biological systems and spanning different mass-spectrometry ionization techniques. We demonstrate that HSNE is able to rapidly identify regions of interest within these large high-dimensionality data sets as well as aid the identification of molecular ions that characterize these regions of interest; furthermore, through clearly separating measurement artifacts, the HSNE analysis exhibits a degree of robustness to measurement batch effects, spatially correlated noise, and mass spectral misalignment.
| S-EPMC5838640 | biostudies-literature
Project description:Flow cytometry is the prototypical assay for multi-parameter single cell analysis, and is essential in vaccine and biomarker research for the enumeration of antigen-specific lymphocytes that are often found in extremely low frequencies (0.1% or less). Standard analysis of flow cytometry data relies on visual identification of cell subsets by experts, a process that is subjective and often difficult to reproduce. An alternative and more objective approach is the use of statistical models to identify cell subsets of interest in an automated fashion. Two specific challenges for automated analysis are to detect extremely low frequency event subsets without biasing the estimate by pre-processing enrichment, and the ability to align cell subsets across multiple data samples for comparative analysis. In this manuscript, we develop hierarchical modeling extensions to the Dirichlet Process Gaussian Mixture Model (DPGMM) approach we have previously described for cell subset identification, and show that the hierarchical DPGMM (HDPGMM) naturally generates an aligned data model that captures both commonalities and variations across multiple samples. HDPGMM also increases the sensitivity to extremely low frequency events by sharing information across multiple samples analyzed simultaneously. We validate the accuracy and reproducibility of HDPGMM estimates of antigen-specific T cells on clinically relevant reference peripheral blood mononuclear cell (PBMC) samples with known frequencies of antigen-specific T cells. These cell samples take advantage of retrovirally TCR-transduced T cells spiked into autologous PBMC samples to give a defined number of antigen-specific T cells detectable by HLA-peptide multimer binding. We provide open source software that can take advantage of both multiple processors and GPU-acceleration to perform the numerically-demanding computations. We show that hierarchical modeling is a useful probabilistic approach that can provide a consistent labeling of cell subsets and increase the sensitivity of rare event detection in the context of quantifying antigen-specific immune responses.
| S-EPMC3708855 | biostudies-literature
Project description:Immunotherapy has emerged as a robust clinical strategy for cancer treatment. PD1/PD-L1 inhibitors have been used as second-line therapy for urothelial carcinoma due to the high tumor mutational burden. Despite the efficacy of the treatment is significant, the response rate is still poor. The tumor immune microenvironment plays a key role in the regulation of immunotherapeutic efficacy. However, a comprehensive understanding of the intricate microenvironment in clinical samples remains unclear. To obtain detailed systematic tumor immune profile, we performed an in-depth immunoassay on 12 human urothelial carcinoma tissue samples and 14 paratumor tissue samples using mass cytometry. Among the large number of cells assayed, we identified 71 T-cell phenotypes, 30 tumor-associated macrophage phenotypes. T cell marker expression profiles showed that almost all T cells in the tumor tissue were in a state of exhaustion. CD38 expression on tumor-associated macrophages (TAMs) was significantly higher than PDL1, and CD38+ TAMs were closely associated with immunosuppression. CD38 may be a more suitable target for immunotherapy in urothelial carcinoma compared to PD1/PDL1. This single-cell analysis of clinical samples expands our insights into the immune microenvironment of urothelial carcinoma and reveals potential biomarkers and targets for immunotherapy development.
| S-EPMC9210814 | biostudies-literature
Project description:Multidimensional single-cell analysis requires approaches to visualize complex data in intuitive 2D graphs. In this regard, t-distributed stochastic neighboring embedding (tSNE) is the most popular algorithm for single-cell RNA sequencing and cytometry by time-of-flight (CyTOF), but its application to polychromatic flow cytometry, including the recently developed 30-parameter platform, is still under investigation. We identified differential distribution of background values between samples, generated by either background calculation or spreading error (SE), as a major source of variability in polychromatic flow cytometry data representation by tSNE, ultimately resulting in the identification of erroneous heterogeneity among cell populations. Biexponential transformation of raw data and limiting SE during panel development dramatically improved data visualization. These aspects must be taken into consideration when using computational approaches as discovery tools in large sets of samples from independent experiments or immunomonitoring in clinical trials.
| S-EPMC6175173 | biostudies-literature
Project description:We present a discrete model of stochastic excitability by a low-dimensional set of delayed integral equations governing the probability in the rest state, the excited state, and the refractory state. The process is a random walk with discrete states and nonexponential waiting time distributions, which lead to the incorporation of memory kernels in the integral equations. We extend the equations of a single unit to the system of equations for an ensemble of globally coupled oscillators, derive the mean field equations, and investigate bifurcations of steady states. Conditions of destabilization are found, which imply oscillations of the mean fields in the stochastic ensemble. The relation between the mean field equations and the paradigmatic Kuramoto model is shown.
| S-EPMC2652551 | biostudies-other
Project description:Understanding the development and function of an organ requires the characterization of all of its cell types. Traditional methods for visualizing and isolating sub-populations of cells are based on mRNA or protein expression of only few known marker genes. The unequivocal identification of a specific marker gene, however, poses a major challenge, particularly if this cell type is rare. Identifying rare cell types, such as stem cells, short-lived progenitors, cancer stem cells, or circulating tumor cells is crucial to acquire a better understanding of normal or diseased tissue biology. To address this challenge we sequenced the transcriptome of hundreds of randomly selected cells from mouse intestinal organoids, cultured self-organizing epithelial structures that contain all cell lineages of the mammalian intestine. Organoid buds, like intestinal crypts, harbor stem cells that continuously differentiate into a variety of cell types, occurring at widely different abundances. Since available computational methods can only resolve more abundant cell types, we developed RaceID, an algorithm for rare cell type identification in complex populations of single cells. We demonstrate that this algorithm can resolve cell types represented by only a single cell in a population of randomly sampled organoid cells. We use this algorithm to identify Reg4 as a novel marker for enteroendocrine cells, a rare population of hormone producing intestinal cells. Next, we use Reg4 expression to enrich for these rare cells and investigate the heterogeneity within this population. Reassuringly, RaceID confirmed the existence of known enteroendocrine lineages, and moreover, discovered novel subtypes, which we subsequently validated in vivo. Having validated RaceID by this proof-of-principle experiment we then apply the algorithm to ex vivo isolated LGR5 positive cells and their direct progeny and demonstrate homogeneity of the stem cell pool. We envision broad applicability of our method for discovering rare cell types and the corresponding marker genes in healthy and diseased organs.
2015-08-19 | GSE62270 | GEO
Project description:The spatial relationships between body parts are a rich source of information for person perception, with even simple pairs of parts providing highly valuable information. Computation of these relationships would benefit from a hierarchical representation, where body parts are represented individually. We hypothesized that the human visual system makes use of such representations. To test this hypothesis, we used adaptation to determine whether observers were sensitive to changes in the length of one body part relative to another. Observers viewed forearm/upper arm pairs where the forearm had been either lengthened or shortened, judging the perceived length of the forearm. Observers then adapted to a variety of different stimuli (e.g., arms, objects, etc.) in different orientations and visual field locations. We found that following adaptation to distorted limbs, but not non-limb objects, observers experienced a shift in perceived forearm length. Furthermore, this effect partially transferred across different orientations and visual field locations. Taken together, these results suggest the effect arises in high level mechanisms specialized for specific body parts, providing evidence for a representation of bodies based on parts and their relationships.
| S-EPMC8142707 | biostudies-literature
Project description:BackgroundThe t-distributed Stochastic Neighbor Embedding (t-SNE) algorithm has emerged as one of the leading methods for visualising high-dimensional (HD) data in a wide variety of fields, especially for revealing cluster structure in HD single-cell transcriptomics data. However, t-SNE often fails to correctly represent hierarchical relationships between clusters and creates spurious patterns in the embedding. In this work we generalised t-SNE using shape-aware graph distances to mitigate some of the limitations of the t-SNE. Although many methods have been recently proposed to circumvent the shortcomings of t-SNE, notably Uniform manifold approximation (UMAP) and Potential of heat diffusion for affinity-based transition embedding (PHATE), we see a clear advantage of the proposed graph-based method.ResultsThe superior performance of the proposed method is first demonstrated on simulated data, where a significant improvement compared to t-SNE, UMAP and PHATE, based on quantitative validation indices, is observed when visualising imbalanced, nonlinear, continuous and hierarchically structured data. Thereafter the ability of the proposed method compared to the competing methods to create faithfully low-dimensional embeddings is shown on two real-world data sets, the single-cell transcriptomics data and the MNIST image data. In addition, the only hyper-parameter of the method can be automatically chosen in a data-driven way, which is consistently optimal across all test cases in this study.ConclusionsIn this work we show that the proposed shape-aware stochastic neighbor embedding method creates low-dimensional visualisations that robustly and accurately reveal key structures of high-dimensional data.
| S-EPMC9660178 | biostudies-literature
Project description:Matching 16S rRNA gene sequencing data to a metabolic reference database is a meaningful way to predict the metabolic function of bacteria and archaea, bringing greater insight to the working of the microbial community. However, some operational taxonomy units (OTUs) cannot be functionally profiled, especially for microbial communities from non-human samples cultured in defective media. Therefore, we herein report the development of Hierarchical micrObial functions Prediction by graph aggregated Embedding (HOPE), which utilizes co-occurring patterns and nucleotide sequences to predict microbial functions. HOPE integrates topological structures of microbial co-occurrence networks with k-mer compositions of OTU sequences and embeds them into a lower-dimensional continuous latent space, while maximally preserving topological relationships among OTUs. The high imbalance among KEGG Orthology (KO) functions of microbes is recognized in our framework that usually yields poor performance. A hierarchical multitask learning module is used in HOPE to alleviate the challenge brought by the long-tailed distribution among classes. To test the performance of HOPE, we compare it with HOPE-one, HOPE-seq, and GraphSAGE, respectively, in three microbial metagenomic 16s rRNA sequencing datasets, including abalone gut, human gut, and gut of Penaeus monodon. Experiments demonstrate that HOPE outperforms baselines on almost all indexes in all experiments. Furthermore, HOPE reveals significant generalization ability. HOPE's basic idea is suitable for other related scenarios, such as the prediction of gene function based on gene co-expression networks. The source code of HOPE is freely available at https://github.com/adrift00/HOPE.
| S-EPMC7874084 | biostudies-literature
Project description:Learning accurate drug representations is essential for tasks such as computational drug repositioning and prediction of drug side effects. A drug hierarchy is a valuable source that encodes knowledge of relations among drugs in a tree-like structure where drugs that act on the same organs, treat the same disease, or bind to the same biological target are grouped together. However, its utility in learning drug representations has not yet been explored, and currently described drug representations cannot place novel molecules in a drug hierarchy. Here, we develop a semi-supervised drug embedding that incorporates two sources of information: (1) underlying chemical grammar that is inferred from chemical structures of drugs and drug-like molecules (unsupervised) and (2) hierarchical relations that are encoded in an expert-crafted hierarchy of approved drugs (supervised). We use the Variational Auto-Encoder (VAE) framework to encode the chemical structures of molecules and use the drug-drug similarity information obtained from the hierarchy to induce the clustering of drugs in hyperbolic space. The hyperbolic space is amenable for encoding hierarchical relations. Both quantitative and qualitative results support that the learned drug embedding can accurately reproduce the chemical structure and recapitulate the hierarchical relations among drugs. Furthermore, our approach can infer the pharmacological properties of novel molecules by retrieving similar drugs from the embedding space. We demonstrate that our drug embedding can predict new uses and discover new side effects of existing drugs. We show that it significantly outperforms comparison methods in both tasks.
| S-EPMC7943198 | biostudies-literature