首页 > 其他分享 >空间转录组cell niche（流程试用版）

空间转录组cell niche（流程试用版）

时间：2024-06-10 16:33:33浏览次数：14  
标签：niche title text 试用版 element cell var fill

作者，Evil Genius

放假了，简单画画图。

最近粉丝问的问题简直就是送命题，空间培训班要不要上，我细数了一下，如果从空间基础分析开始将讲，R，python版本，当然是更新后的code，包括所有的2024年更新的新方法，新的绘图方法，再加上HD、Xenium，文献解读，我初步排了一下课，没40节课根本上不完。

即使去年的课程不包含在内，新更新的方法 + python版本讲起来也够30节课了。

但也不得不说，现在的分析越来越高端了。

而且似乎分子niche和细胞niche是必要的分析

spatial niche analyses

这里备注一下cell niche的脚本示例，以备后续使用

library(Seurat)
library(SeuratObject)
library(tidyverse)
library(gplots)
library(tibble)


## CREATE COLOR LISTS ----
color_list <- list(`aCap` = "#aa4e26",
                   `Arteriole` = "#d59445",
                   `gCap` = "#cdbf5e",
                   `DCs/Lymphatic` = "#e46300",
                   `Venous` = "#975f47",
                   `Proliferating Endothelial` = "#dacf3c",  
                   
                   `Activated FBs (CTHRC1+/FAP+)` = "#765fa5",        
                   `Activated FBs` = "#765fa5",  
                   `Adventitial FBs` = "#01b0ea",                     
                   `Alveolar FBs` = "#6d4bce",        
                   `Fibroblasts` = "#7cb3f5",                        
                   `HAS1+ Fibroblasts` = "#0065f5",  
                   `Lipofibroblasts` = "#c6adf2",    
                   `Mesothelial` = "#374ca3",    
                   `Peribronchial FBs` = "#9da7ff",
                   `PLIN2+ Fibroblasts` = "#4c4e8c",
                   `SMCs` = "#0192fe",    
                   `Proliferating Mesenchymal` = "#7895ff",
                   
                   `B cells` = "#d84265",
                   `NK cells` = "#b16f6f",
                   `T-cells` = "#e62f4c",  
                   `Plasma` = "#df4738",
                   `pDCs` = "#7b39c7",
                   `Proliferating Lymphoid` = "#993432",
                   
                   `cDCs` =  "#7a4ca1",
                   `DCs` = "#ac67dc",
                   `FABP4+ Macrophages` = "#cc85d9",
                   `Macrophages` = "#dc48bd",
                   `Mast` = "#bf4ce3",
                   `Monocyte-derived Macrophages` = "#963095",
                   `Monocyte-derived Macrophages (FCN1+)`  = "#d497ba",                          
                   `SPP1+ Macrophages` = "#9e3c6d",
                   `Proliferating Myeloid` = "#e25292",  
                   
                   `AT1` = "#6cd66d",                                 
                   `AT2` = "#34491d",                                                        
                   `Basal` = "#acde47",                                                     
                   `Ciliated` = "#b9b700",                          
                   `Differentiating Ciliated` = "#01eaaa",           
                   `KRT5-/KRT17+` = "#a0cd9e",                                        
                   `MUC5B+ Secretory` = "#31a200",      
                   `SCGB3A2+` = "#ccd87a",
                   `SCGB3A2+/SCGB1A1+` = "#006e3e",  
                   `Transitional AT2` = "#778d39",   
                   `Proliferating Epithelial` = "#586551")

nuclei_niche_color_list <- list(`1` = "#ffa26d",
                                `2` = "#1036b5",
                                `3` = "#cabd00",
                                `4` = "#a6513c",
                                `5` = "#84cd5d",
                                `6` = "#72c7ff",
                                `7` = "#d10040",
                                `8` = "black",
                                `9` = "#71c8a5",
                                `10` = "#6f774b",
                                `11` = "#8E76DB",
                                `12` = "#924373")



## SET ENVIRONMENT AND LOAD DATA ----
library(Seurat)
library(ggplot2)
library(tidyverse)
library(Cairo)

# Better UMAP function
pretty_umap <- theme(axis.text = element_blank(), axis.ticks = element_blank(),
                     axis.title = element_text(hjust = 1))

# Load data
xenium <- readRDS("/scratch/avannan/late_IPF_spatial/xenium/rds_files/xenium_full_CT_051923.rds")

# Set directories for saving plots
scratch_dir <- "/Users/avannan/Documents/nuclei_niches/"


## CREATE COLOR LISTS ----
color_list <- list(`aCap` = "#aa4e26",
                   `Arteriole` = "#d59445",
                   `gCap` = "#cdbf5e",
                   `DCs/Lymphatic` = "#e46300",
                   `Venous` = "#975f47",
                   `Proliferating Endothelial` = "#dacf3c",  
                   
                   `Activated FBs (CTHRC1+/FAP+)` = "#765fa5",        
                   `Activated FBs` = "#765fa5",  
                   `Adventitial FBs` = "#01b0ea",                     
                   `Alveolar FBs` = "#6d4bce",        
                   `Fibroblasts` = "#7cb3f5",                        
                   `HAS1+ Fibroblasts` = "#0065f5",  
                   `Lipofibroblasts` = "#c6adf2",    
                   `Mesothelial` = "#374ca3",    
                   `Peribronchial FBs` = "#9da7ff",
                   `PLIN2+ Fibroblasts` = "#4c4e8c",
                   `SMCs` = "#0192fe",    
                   `Proliferating Mesenchymal` = "#7895ff",
                   
                   `B cells` = "#d84265",
                   `NK cells` = "#b16f6f",
                   `T-cells` = "#e62f4c",  
                   `Plasma` = "#df4738",
                   `pDCs` = "#7b39c7",
                   `Proliferating Lymphoid` = "#993432",
                   
                   `cDCs` =  "#7a4ca1",
                   `DCs` = "#ac67dc",
                   `FABP4+ Macrophages` = "#cc85d9",
                   `Macrophages` = "#dc48bd",
                   `Mast` = "#bf4ce3",
                   `Monocyte-derived Macrophages` = "#963095",
                   `Monocyte-derived Macrophages (FCN1+)`  = "#d497ba",                          
                   `SPP1+ Macrophages` = "#9e3c6d",
                   `Proliferating Myeloid` = "#e25292",  
                   
                   `AT1` = "#6cd66d",                                 
                   `AT2` = "#34491d",                                                        
                   `Basal` = "#acde47",                                                     
                   `Ciliated` = "#b9b700",                          
                   `Differentiating Ciliated` = "#01eaaa",           
                   `KRT5-/KRT17+` = "#a0cd9e",                                        
                   `MUC5B+ Secretory` = "#31a200",      
                   `SCGB3A2+` = "#ccd87a",
                   `SCGB3A2+/SCGB1A1+` = "#006e3e",  
                   `Transitional AT2` = "#778d39",   
                   `Proliferating Epithelial` = "#586551")

nuclei_niche_color_list <- list(`1` = "#ffa26d",
                                `2` = "#1036b5",
                                `3` = "#cabd00",
                                `4` = "#a6513c",
                                `5` = "#84cd5d",
                                `6` = "#72c7ff",
                                `7` = "#d10040",
                                `8` = "black",
                                `9` = "#71c8a5",
                                `10` = "#6f774b",
                                `11` = "#8E76DB",
                                `12` = "#924373")


## FIX COORDINATES ----
# Set spatial rows and columns
xenium$spatial_row <- ""
xenium$spatial_row[xenium$sample %in% c("VUHD116A", "VUILD102LF", "VUHD095", "VUILD105MF", "VUHD069")] <- "R1"
xenium$spatial_row[xenium$sample %in% c("VUILD96MF", "VUILD102MF", "VUHD116B", "VUILD104LF", "VUHD113", "VUILD105LF")] <- "R2"
xenium$spatial_row[xenium$sample %in% c("VUILD96LF", "VUILD107MF", "VUILD48MF", "VUILD104MF", "VUILD48LF")] <- "R3"
xenium$spatial_row[xenium$sample %in% c("VUILD115", "VUILD110", "THD0011", "TILD117MF")] <- "R4"
xenium$spatial_row[xenium$sample %in% c("VUILD78MF", "TILD117LF", "VUILD91MF")] <- "R5"
xenium$spatial_row[xenium$sample %in% c("VUILD106", "THD0008", "VUILD91LF", "VUILD78LF", "TILD175")] <- "R6"
xenium$spatial_column <- ""
xenium$spatial_column[xenium$sample %in% c("VUILD96MF", "VUILD115", "VUILD106")] <- "C1"
xenium$spatial_column[xenium$sample %in% c("VUHD116A", "VUILD102MF", "VUILD96LF")] <- "C2"
xenium$spatial_column[xenium$sample %in% c("VUILD102LF", "VUHD116B", "VUILD107MF", "VUILD110", "THD0008")] <- "C3"
xenium$spatial_column[xenium$sample %in% c("VUHD095", "VUILD104LF", "VUILD48MF", "VUILD78MF", "VUILD91LF")] <- "C4"
xenium$spatial_column[xenium$sample %in% c("VUILD105MF", "VUHD113", "VUILD104MF", "THD0011", "TILD117LF", "VUILD78LF")] <- "C5"
xenium$spatial_column[xenium$sample %in% c("VUHD069", "VUILD105LF", "VUILD48LF", "TILD117MF", "VUILD91MF", "TILD175")] <- "C6"

# Adjust coordinates based on row and column positions
xenium$super_adj_y_centroid <- xenium$adj_y_centroid
xenium$super_adj_y_centroid[xenium$spatial_row == "R1"] <- xenium$adj_y_centroid[xenium$spatial_row == "R1"] + 4400
xenium$super_adj_y_centroid[xenium$spatial_row == "R2"] <- xenium$adj_y_centroid[xenium$spatial_row == "R2"] + 2200
xenium$super_adj_y_centroid[xenium$spatial_row == "R4"] <- xenium$adj_y_centroid[xenium$spatial_row == "R4"] -1000
xenium$super_adj_y_centroid[xenium$spatial_row == "R5"] <- xenium$adj_y_centroid[xenium$spatial_row == "R5"] - 2200
xenium$super_adj_y_centroid[xenium$spatial_row == "R6"] <- xenium$adj_y_centroid[xenium$spatial_row == "R6"] - 4400
xenium$super_adj_x_centroid <- xenium$adj_x_centroid
xenium$super_adj_x_centroid[xenium$spatial_column == "C1"] <- xenium$adj_x_centroid[xenium$spatial_column == "C1"] - 3000
xenium$super_adj_x_centroid[xenium$spatial_column == "C2"] <- xenium$adj_x_centroid[xenium$spatial_column == "C2"] - 2000
xenium$super_adj_x_centroid[xenium$spatial_column == "C3"] <- xenium$adj_x_centroid[xenium$spatial_column == "C3"] - 1000
xenium$super_adj_x_centroid[xenium$spatial_column == "C5"] <- xenium$adj_x_centroid[xenium$spatial_column == "C5"] + 2200
xenium$super_adj_x_centroid[xenium$spatial_column == "C6"] <- xenium$adj_x_centroid[xenium$spatial_column == "C6"] + 4400

# Add "super" spatial dimension reduction object
position_xy <- cbind(xenium$super_adj_x_centroid,
                     xenium$super_adj_y_centroid)
row.names(position_xy) <- row.names([email protected])
colnames(position_xy) <- c("Super_SP_1", "Super_SP_2")
xenium[["super_sp"]] <- CreateDimReducObject(
  embeddings = position_xy, key = "Super_SP_", assay = DefaultAssay(xenium))



# SET VARIABLE NAMES ----
# Broad, 25, 8-12 niches
ct_var = xenium$broad_CT5
ct_name_var = "broad_CT5"
num_neighbors_list = 25
num_niches_list = c(8, 9, 10, 11, 12)


## RUN NICHE ANALYSIS ----
for (j in 1:length(num_neighbors_list)) {
  num_neighbors = num_neighbors_list[j]
  print(paste0("Neighbors K: ", num_neighbors))
  
  for (i in num_niches_list) {
    num_niches = i
    print(paste0("Niches K: ", num_niches))
    
    # Set directory names for saving files
    niche_name_var <-  paste0(ct_name_var, "_nichek", num_niches, "_neighborsk",
                              num_neighbors)
    plot_save_dir <- paste0(scratch_dir, niche_name_var, "/")
    
    print(paste0("Saving plots to ", plot_save_dir))

    ## CREATE FOV AND FIND CELL NICHES ----
    # Create FOV
    coords <- [email protected][, c("super_adj_x_centroid", "super_adj_y_centroid")]
    xenium[["fov"]]  <- CreateFOV(coords, assay = "RNA", type = "centroids")

    Sys.time()
    xenium <- BuildNicheAssay(object = xenium, fov = "fov", group.by = ct_name_var,
                              niches.k = num_niches, neighbors.k = num_neighbors)
    Sys.time()

    # Set more variables
    niche_var <- xenium$niches
    [email protected][[niche_name_var]] <- xenium$niches

    ## GENERAL SUMMARY FIGURES ----
    # Niche by Cell Type
    CairoPNG(paste0(plot_save_dir, niche_name_var, "_by_ct.png"), width = 10, height = 6, units = "in", dpi = 300)
    a <- table(ct_var, niche_var) %>%
      as.data.frame() %>%
      ggplot(aes(x = niche_var, y = Freq, fill = ct_var)) +
      geom_col(position = position_fill()) +
      scale_y_continuous(expand = c(0, 0)) +
      scale_fill_manual(values = color_list) +
      labs(x = "Nuclei Niche", y = "Cell Type Proportion", fill = "", title = niche_name_var) +
      theme_bw() +
      theme(axis.text = element_text(color = "black"),
            plot.title = element_text(hjust = 0.5))
    print(a)
    dev.off()
    
    # Cell Type by Niche
    CairoPNG(paste0(plot_save_dir, "ct_by_", niche_name_var, ".png"), width = 10, height = 6, units = "in", dpi = 300)
    a <- table(ct_var, niche_var) %>%
      as.data.frame() %>%
      ggplot(aes(x = ct_var, y = Freq, fill = niche_var)) +
      geom_col(position = position_fill(), width = 0.8) +
      scale_y_continuous(expand = c(0, 0)) +
      scale_fill_manual(values = nuclei_niche_color_list) +
      labs(y = "Nuclei Niche Proportion", fill = "", title = niche_name_var) +
      theme_bw(base_size = 12) +
      theme(axis.text = element_text(color = "black"),
            axis.text.x = element_text(angle = 45, hjust = 1),
            axis.title.x = element_blank(),
            plot.title = element_text(hjust = 0.5))
    print(a)
    dev.off()
    
    # Niche by Sample Type
    CairoPNG(paste0(plot_save_dir, niche_name_var, "_by_sample_type.png"), width = 5, height = 5, units = "in", dpi = 300)
    a <- table(xenium$sample_type, niche_var) %>% 
      as.data.frame() %>%
      mutate(Var1 = ordered(Var1, levels = c("Unaffected", "LF", "MF", "ILD"))) %>%
      ggplot(aes(x = Var1, y = Freq, fill = niche_var)) +
      geom_col(position = position_fill()) +
      scale_y_continuous(expand = c(0, 0)) +
      scale_fill_manual(values = nuclei_niche_color_list) +
      theme_classic(base_size = 14) +
      labs(x = "", y = "Niche Proportions", title = niche_name_var) +
      theme(legend.title = element_blank(), 
            axis.text = element_text(color = "black"),
            plot.title = element_text(hjust = 0.5))
    print(a)
    dev.off()
    
    # Niche by Sample
    CairoPNG(paste0(plot_save_dir, niche_name_var, "_by_sample.png"), width = 5, height = 5, units = "in", dpi = 300)
    a <- table(xenium$sample, niche_var) %>%
      as.data.frame() %>%
      mutate(Var1 = ordered(Var1, levels = c("VUHD116A", "VUHD116B", "VUHD095", "VUHD113", 
                                             "VUHD069", "THD0011", "THD0008", "TILD117LF",
                                             "VUILD96LF", "VUILD91LF", "VUILD78LF", "VUILD105LF",
                                             "VUILD48LF", "VUILD102LF", "VUILD104LF", "VUILD78MF",
                                             "VUILD96MF", "VUILD48MF", "VUILD104MF", "TILD117MF", 
                                             "VUILD102MF", "VUILD107MF", "VUILD91MF",  "VUILD110", 
                                             "VUILD105MF", "TILD175", "VUILD115", "VUILD106"))) %>%
      ggplot(aes(x = Var1, y = Freq, fill = niche_var)) +
      geom_col(position = position_fill()) +
      scale_y_continuous(expand = c(0, 0)) +
      scale_fill_manual(values = nuclei_niche_color_list) +
      theme_classic(base_size = 14) +
      labs(x = "", y = "Niche Proportions", title = niche_name_var) +
      theme(legend.title = element_blank(),
            axis.text = element_text(color = "black"),
            axis.text.x = element_text(angle = 45, hjust = 1),
            plot.title = element_text(hjust = 0.5))
    print(a)
    dev.off()
    
    ## SPLIT SAMPLE PLOTS ----
    # Split samples
    sample_ids <- unique(xenium$sample)
    sample_list <- lapply(sample_ids, function(XX) { subset(xenium, subset = sample == XX) })
    names(sample_list) <- sample_ids
    
    # Create DimPlots for each sample
    dimplot_list <- lapply(sample_ids, function(XX) {
      DimPlot(sample_list[[XX]], reduction = "sp", 
              group.by = "niches", cols = nuclei_niche_color_list, 
              raster = FALSE) + pretty_umap + 
        ggtitle(paste0(XX, "\n", niche_name_var))
    })
    names(dimplot_list) <- sample_ids
    
    # Create PNGs for each figure
    for (sm in sample_ids) {
      CairoPNG(paste0(plot_save_dir, sm, "_", niche_name_var, ".png"), width = 7.5, height = 7.5, units = "in", dpi = 300)
      print(dimplot_list[[sm]])
      dev.off()
    }
    # ave object
    saveRDS(xenium, paste0(scratch_dir, "xenium_", niche_name_var, ".rds"))
  }
}

saveRDS(xenium, paste0(scratch_dir, "xenium_", paste0(ct_name_var, "_neighborsk", num_neighbors_list), ".rds"))

spatial niche analyses

生活很好，有你更好

标签：niche,title,text,试用版,element,cell,var,fill
From： https://blog.csdn.net/weixin_53637133/article/details/139552012

相关文章

  • 数据结构和矩阵细节用法：double、cell和complex #matlab
    建立矩阵用[]；矩阵的同一行内的元素用逗号或者空格隔开矩阵的不同行的元素用分号隔开eg.       矩阵     A=123456789在matlab中矩阵A表示为：clc;clear;A=[1,2,3;4,5,6;7,8,9]; %或者A=[123;456;789];小矩阵构建大矩阵假设：A=[......
  • 覆盖路径规划经典算法 The Boustrophedon Cellular Decomposition 详解
    2000年一篇论文CoverageofKnownSpaces:TheBoustrophedonCellularDecomposition横空出世，解决了很多计算机和机器人领域的覆盖路径问题，今天我来详细解读这个算法。TheBoustrophedonCellularDecomposition算法详解这篇论文标题为"CoveragePathPlanning:TheB......
  • Xenium In Situ数据的质量评估及优化方案（试用版）
    作者，EvilGenius作为公司的分析人员，多方法比较和优化也是必修课。今天我们来看看Xenium的一些分析点及可能得优化方案Xenium本身的细胞分割更加适合于规整的、圆形或者其他形状不太奇怪的细胞，这一点得到了上海10X技术支持的确认，不过对于大多数样本，这个方法就够用了。这一篇......
  • ISS空间转录组的细胞分割算法汇总(stardist、cellpose、QuPath、SCS)
    作者，EvilGenius我们来了解一下关于HE图片细胞分割的一些算法，以备后续的使用我们在做Xenium或者其他ISS技术的数据时候，通常都要获得一个polygons文件，直译过来就是多边形文件，其实就是我们进行的图像细胞分割文件。其实ISS技术由来已久，不过现在由于10Xxenium的发布，对于空间原......
  • 易基因： ChIP-seq等揭示酿酒酵母Tho2介导Nrd1调控衰老相关基因表达｜Aging Cell
    大家好，这里是专注表观组学十余年，领跑多组学科研服务的易基因。衰老的特点是分子、细胞和器官损伤的逐渐积累，导致生物功能下降，对疾病和死亡的敏感性增加。酿酒酵母（Saccharomycescerevisiae）作为首个实现完整基因组测序的模式真核细胞，具有繁殖迅速、培养条件简单、基因组小等优点，并......
  • Flutter 中的 TableCell 小部件：全面指南
    Flutter中的TableCell小部件：全面指南Flutter是一个功能强大的UI框架，由Google开发，允许开发者使用Dart语言构建跨平台的移动、Web和桌面应用。在Flutter的丰富组件库中，TableCell是一个用于创建表格单元格的组件，它提供了灵活的方式来自定义表格的内容和布局。本......
  • 关于前端上传excell时间的问题
    当前端导入excell里的数据时，只能获取到下面类似的这种数据Excel存储的日期是从1900年1月1日开始按天数来计算的，也就是说1900年1月1日在Excel中是1。转化的思路和对Excel中的错误的处理因为时间戳是从1970年1月1日算起的（时间戳为0的时候是1970年1月1日）也就是说newDate(0).toLo......
  • 浅谈C#中取消令牌CancellationTokenSource
    基础操作CancellationTokenSourcecancellationTokenSource=newCancellationTokenSource();cancellationTokenSource.Token.Register(()=>{Console.WriteLine("取消了1111");});......
  • cellchat安装问题
    运行linux安装踩坑失败，留个痕迹install-cbiocondabioconductor-biobase=2.54.0-y报错ledwithinitialfrozensolve.Retryingwithflexiblesolve.Solvingenvironment:failedwithrepodatafromcurrent_repodata.json,willretrywithnextrepodatasource.Ree......
  • windows安装cellchat
    检查包位置.libPaths('D:/R/R-4.3.1/library')确认R版本和包路径一致用R=4.3.1下载https://github.com/jinworks/CellChat安装cellchat前根据说明安装依赖https://github.com/jinworks/CellChat切换到CellChat-master.zip所在位置,删掉zip/src里的.o.sosetwd('D:/R/R-4.3......

赞助商

阅读排行

网站主页关于我们联系我们网站地图

本网站内容转载自其他媒体，侵权联系[admin##ips99.com]。

Copyright © 2020-2023 IPS99 版权所有 IPS99