ClusterMap:用于空间基因表达的多尺度聚类分析 | 空间转录组分析工具推荐

在空间背景下量化RNA是了解复杂组织中基因表达和调控的关键。原位转录组方法可以在完整的组织中产生空间分辨率的RNA图谱。然而,目前还缺乏一个统一的计算工具来综合分析原位转录组数据。2021年10月,Nature Communications发表了一个无监督和无注释的计算工具:ClusterMap,其在二维和三维空间将RNA精确地聚类到亚细胞结构、细胞体和组织区域中,并在不同的组织类型(包括小鼠大脑、胎盘、肠道和人类心脏器官)中表现稳定。ClusterMap广泛适用于各种原位转录组技术,从高维转录组图谱图像中揭示基因表达模式、细胞生态位和组织结构原理。

为了绘制细胞和组织中基因表达的空间异质性,已经开发了大量基于图像的原位转录组学方法(例如STARmap、Fiseq、ISS、MERFISH、seqFISH、osmFISH等),提供了完整组织中亚细胞RNA定位的图谱。然而,从高维空间转录组数据中直接提取生物模式的低维表示仍然具有挑战性。

ClusterMap是什么?

ClusterMap是基于两个关键的生物学现象。首先,细胞内RNA分子的密度高于细胞外;其次,不同基因编码的细胞RNA在不同的亚细胞位置、细胞类型和组织区域富集。因此,开发团队推断,通过对RNA的物理密度和基因身份进行联合聚类,可以直接从原位转录组数据中确定有生物学意义的模式和结构。随后&#x

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
多距离空间聚类分析是一种常见的空间分析方法,它可以帮助我们识别在不同距离阈值下的空间聚类模式。使用arcpy中的Spatial Analyst扩展可以方便地实现这个功能,以下是一个示例代码: ```python import arcpy from arcpy import sa def multi_distance_clustering(input_fc, distances): # 初始化空间参考 spatial_ref = arcpy.Describe(input_fc).spatialReference # 创建空白的分类图像 temp_raster = arcpy.NumPyArrayToRaster(numpy.zeros((spatial_ref.height, spatial_ref.width)), arcpy.Point(spatial_ref.XMin, spatial_ref.YMin), spatial_ref.meanCellWidth, spatial_ref.meanCellHeight) # 循环不同的距离阈值,进行聚类分析 for distance in distances: # 使用FocalStatistics函数计算每个像元在指定距离范围内有多少个邻居 neighborhood = sa.NbrCircle(distance, "MAP") neighbor_count = sa.FocalStatistics(temp_raster, neighborhood, "SUM") # 将邻居数大于等于5的像元分类为1,否则分类为0 cluster_raster = sa.Con(neighbor_count >= 5, 1, 0) # 根据分类结果创建要素类 cluster_fc = arcpy.RasterToPolygon_conversion(cluster_raster, "in_memory\\clusters", "NO_SIMPLIFY") # 将分类结果写入原始要素类的属性表中 arcpy.AddField_management(input_fc, "cluster_" + str(distance), "SHORT") arcpy.JoinField_management(input_fc, "OBJECTID", cluster_fc, "GRIDCODE", "FID") return input_fc ``` 在这个函数中,我们首先使用arcpy.NumPyArrayToRaster函数创建一个空白的分类图像,并循环不同的距离阈值进行聚类分析。对于每个距离阈值,我们使用FocalStatistics函数计算每个像元在指定距离范围内有多少个邻居,并将邻居数大于等于5的像元分类为1,否则分类为0。然后,我们将分类结果写入原始要素类的属性表中,并返回更新后的要素类。请注意,这个示例代码中的距离阈值是在函数参数中作为一个距离列表传递的。如果需要,您可以根据具体情况进行调整。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值