Tag: UMAP

Novel diagnostic biomarkers for keloid based on GEO database

Introduction Keloid is excessive fibrosis of the skin that extends beyond the area of injury and does not regress.1 Keloid can occur in the joints and mouth after several years of severe injury, including burns, chemical injury, wound, and surgical incision.2 Keloids on the joints affect the quality of life,…

Continue Reading Novel diagnostic biomarkers for keloid based on GEO database

Identification of a regulatory pathway inhibiting adipogenesis via RSPO2

Integration of APC scRNA-seq data reveals heterogeneity of adipocyte progenitor cells In a previous study9, we defined Lin−Sca1+CD142+ APCs as adipogenesis regulatory (Areg) cells and demonstrated that these cells are both refractory toward adipogenesis and control adipocyte formation of APCs through paracrine signaling. In contrast, Merrick et. al.4 observed that Lin−CD142+ cells…

Continue Reading Identification of a regulatory pathway inhibiting adipogenesis via RSPO2

Integrating Bulk RNA-seq data with Single cell RNA seq data

Integrating Bulk RNA-seq data with Single cell RNA seq data 0 Hello all, recently, I had been trying to integrate bulk RNAseq data into single-cell data where I treat each sample in my bulk RNAseq data as a single cell and integrate it into the single-cell data based on the…

Continue Reading Integrating Bulk RNA-seq data with Single cell RNA seq data

Adding numbers and characters to legend key in ggplot2 of UMAP clusters

Adding numbers and characters to legend key in ggplot2 of UMAP clusters 0 Hi everyone, I have a UMAP cluster, however there are so many clusters that the descriptions look clunky if i put them on the umap…but then there are too many colors if it just colors. So, i…

Continue Reading Adding numbers and characters to legend key in ggplot2 of UMAP clusters

Non-genetic determinants of malignant clonal fitness at single-cell resolution

1. Turajlic, S., Sottoriva, A., Graham, T. & Swanton, C. Resolving genetic heterogeneity in cancer. Nat. Rev. Genet. 20, 404–416 (2019). CAS  PubMed  Google Scholar  2. Marine, J. C., Dawson, S. J. & Dawson, M. A. Non-genetic mechanisms of therapeutic resistance in cancer. Nat. Rev. Cancer 20, 743–756 (2020). CAS …

Continue Reading Non-genetic determinants of malignant clonal fitness at single-cell resolution

tSNE and UMAP of scATAC-seq data looks like spaghetti

tSNE and UMAP of scATAC-seq data looks like spaghetti 0 I would like to use R to generate cluster my 20k cells from a single cell ATAC-seq experiment. I ran PCA then selected the first 50 components, which were put into tSNE’s normalize_input() then Rtsne(). This is the result I…

Continue Reading tSNE and UMAP of scATAC-seq data looks like spaghetti

Correct usage of FindConservedMarkers() in Seurat

Correct usage of FindConservedMarkers() in Seurat 0 Dear all, I have a Seurat object of a certain cell type with a UMAP of 7 clusters. I also have information about the sample’s origin (primary tumor/metastatic) in my metadata. Looking at the UMAP I can clearly see that clusters 1 and…

Continue Reading Correct usage of FindConservedMarkers() in Seurat

Bioconductor – GEOexplorer (development version)

DOI: 10.18129/B9.bioc.GEOexplorer     This is the development version of GEOexplorer; to use it, please install the devel version of Bioconductor. GEOexplorer: an R/Bioconductor package for gene expression analysis and visualisation Bioconductor version: Development (3.14) GEOexplorer is a Shiny app that enables exploratory data analysis and differential gene expression of…

Continue Reading Bioconductor – GEOexplorer (development version)

WNN in Seurat

Dear all, I am trying to follow the WNN vignette here satijalab.org/seurat/articles/weighted_nearest_neighbor_analysis.html After the steps below, I would like to annotate my clusters, hence I need to know the markers which best represent each cluster. pbmc <- FindMultiModalNeighbors(pbmc, reduction.list = list(“pca”, “lsi”), dims.list = list(1:50, 2:50)) pbmc <- RunUMAP(pbmc, nn.name

Continue Reading WNN in Seurat

The Biostar Herald for Friday, September 03, 2021

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from zx8754, Istvan Albert, and was edited by…

Continue Reading The Biostar Herald for Friday, September 03, 2021

A question about scRNA-seq analysis

A question about scRNA-seq analysis 0 Dear all, please may I ask for a suggestion : I have a scRNA-seq dataset, which has 2 groups (Control and Model). I would like to process all the samples into cluster by tSNE or UMAP , then classification or differentiation the cell types….

Continue Reading A question about scRNA-seq analysis

Best compression for single cell RNA-seq object

Best compression for single cell RNA-seq object 0 I generated a scRNA-seq object (counts, PCA, UMAP embeddings, DEGs etc.) in Scanpy or Seurat. What is the best data structure to store this in to reduce the size of the object? I’m considering H5AD (scanpy/anndata), RDS or H5Seurat (Seurat), or Loom…

Continue Reading Best compression for single cell RNA-seq object

WNN (Seurat v4) vs. totalVI (scvi-tools) for CITE-seq data

WNN (Seurat v4) vs. totalVI (scvi-tools) for CITE-seq data 0 I want to build a UMAP from CITE-seq data with a joint embedding of the scRNA-seq and protein ab data. What’s the ‘best’ method in terms of representing the most accurate embedding? In the totalVI paper, they say totalVI >…

Continue Reading WNN (Seurat v4) vs. totalVI (scvi-tools) for CITE-seq data

How to transform the deg gene list from seurat to a gene list input to clusterProfiler compareCluster ?

Sorry for lateness, I wanted to do something similar. This is what I did for reference: Using a Seurat generated gene list for input into ClusterProfiler to see the GO or KEGG terms per cluster. I’ll keep the meat and potatoes of the Seurat vignette in this tutorial: library(dplyr) library(Seurat)…

Continue Reading How to transform the deg gene list from seurat to a gene list input to clusterProfiler compareCluster ?

VisCello Shiny app code for hosting snRNA-seq data

I am trying to set up VisCello shiny app for hosting some of our single cell data, analyzed in Seurat: github.com/qinzhu/VisCello I am using the following code though running into a problem, the app launches however I dont see the umaps and when I click over to differential expression the…

Continue Reading VisCello Shiny app code for hosting snRNA-seq data

UMAP of TRA/B

Hello, I have output from a single cell sequencing run that has both the VDJ and gene expression data. For the same cells, we also used a hybrid capture approach to sequence the TCR sequences. I have compared the TCR sequences across the two approaches and I have found a…

Continue Reading UMAP of TRA/B

UMAP vs “rigorous” t-SNE

UMAP vs “rigorous” t-SNE 1 I’ve heard a lot of people discussing UMAP recently as though it has essentially superseded t-SNE for visualizing scRNA-seq data. UMAP is certainly impressive, but it seems to me that there are a lot of things one can do to pretty dramatically improve the output…

Continue Reading UMAP vs “rigorous” t-SNE

Experience with Cytobank

Experience with Cytobank 0 Dear all We are medical doctors and immunologists working with 16 color phenotype FCS files created on patient blood PBMC samples. We would like to independently work with a clustering algorithm that is intuitive and easy to manage. We have been exploring Cytobank briefly but we…

Continue Reading Experience with Cytobank

Can someone explain the differences between various 1000 genome project and gnomad call sets? Also any straightforward PCA implementation on them?

I’ve been trying to delve into the data from whole genome sequencing, specifically by looking at the already existing data in the 1000 genome project and gnomad, and I have a lot of questions. Does gnomAD contain the 1000gp samples? I’ve found many vcf including these: ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/supporting/hd_genotype_chip/ ftp.1000genomes.ebi.ac.uk/vol1/ftp/phase3/integrated_sv_map/ ftp.1000genomes.ebi.ac.uk/vol1/ftp/phase1/analysis_results/integrated_call_sets/ gnomad.broadinstitute.org/downloads…

Continue Reading Can someone explain the differences between various 1000 genome project and gnomad call sets? Also any straightforward PCA implementation on them?

Corelate TCR data to clusters/GEX/CITEseq data

Corelate TCR data to clusters/GEX/CITEseq data 1 Hello everyone, I just added my TCR VDJ data as metadata to my Seurat object (as described in the tutorial here). So, I basically ended up with two different collumns of metadata where my barcodes are assigned to the clonotypes and the cdr3…

Continue Reading Corelate TCR data to clusters/GEX/CITEseq data