如何在R语言中通过染色体坐标（Chromosome:BP）获取SNP名称？

完全可行，下面是两种常用的实现方案：

方法一：使用biomaRt包（在线查询，数据更新及时） #

biomaRt可对接Ensembl数据库，在线查询最新的SNP信息，适合需要实时数据的场景：

• 安装并加载包：

install.packages("biomaRt")
library(biomaRt)

• 连接对应物种的Ensembl数据集（以人类hg38版本为例，其他物种或版本可替换数据集名称）：

ensembl <- useMart("ensembl", dataset = "hsapiens_gene_ensembl")

• 准备目标坐标信息（示例为chr1:1000000，注意染色体号无需带chr前缀，与Ensembl格式统一）：

target_chr <- "1"
target_pos <- 1000000

• 执行查询，获取SNP ID及对应坐标：

snp_result <- getBM(
  attributes = c("rsid", "chromosome_name", "start_position"),
  filters = c("chromosome_name", "start"),
  values = list(chromosome_name = target_chr, start = target_pos),
  mart = ensembl
)

查询结果snp_result中的rsid列就是对应的SNP名称。如果需要查询坐标区间，只需增加end过滤器并传入区间终点值即可。

方法二：使用SNPlocs系列离线包（无网络环境适用） #

如果没有网络，可使用Bioconductor提供的离线SNP位置包，以人类dbSNP155版本为例：

• 安装并加载包（需先安装BiocManager）：

if (!require("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("SNPlocs.Hsapiens.dbSNP155")
library(SNPlocs.Hsapiens.dbSNP155)

• 构建基因组范围对象，指定目标坐标：

library(GenomicRanges)
target_gr <- GRanges(seqnames = "chr1", ranges = IRanges(start = 1000000, end = 1000000))

• 查询对应位置的SNP：

matched_snps <- snpsByOverlaps(SNPlocs.Hsapiens.dbSNP155, target_gr)
snp_rsids <- rsids(matched_snps)

snp_rsids即为该位置对应的所有SNP名称列表。

注意事项 #

• 确保使用的基因组版本（如hg19/hg38）与你的坐标数据一致，不同版本的坐标对应关系可能存在差异；
• 染色体格式要统一：部分数据库要求带chr前缀，部分不需要，查询前需确认包或数据库的格式要求；
• 单个位置可能对应多个SNP，查询结果会返回所有匹配条目。

内容的提问来源于stack exchange，提问作者user169605