使用ANNOVAR匹配 variant rsid

发表于 2023-9-24 16:37:47

最近一直在尝试使用ANNOVAR给我从网上下载的GWAS summary data 进行注释，主要就是别人给了CHR:POS, 但没有rsid, 我试着用avsnp150/avsnp147/snp138 注释，但发现成功匹配的都是ins/del，没有位点突变那种（e.g., C/T），我想知道是我文件格式不对还是这个软件的问题，还是哪个步骤错了，谢谢！

发表于 2023-9-24 16:37:48

ANNOVAR是一个很有用的工具，可以用于注释基因组中的突变位点。根据你描述的情况，可能有几个原因导致无法匹配到 rsid。

首先，确保你的 GWAS summary data 文件与 ANNOVAR 能够兼容。GWAS summary data 一般以 tab 分隔的文本文件形式呈现，其中包含列名和对应的信息。确保你的文件格式正确，并且包含 CHR、POS 和其他必要的列。

其次，确认你使用的数据库版本是否适合你的数据。不同版本的数据库可能会存在差异，特别是对于一些较新的变异位点。你提到使用了 avsnp150/avsnp147/snp138 数据库，但仍然无法匹配到 rsid。可以尝试使用更新的数据库版本，例如 avsnp151 或最新的数据库版本。

另外，如果你在注释中只看到了插入/删除类型的变异位点而没有突变位点，则可能是由于数据中的变异位点没有被正确识别或注释。这可能是由于数据源的问题，例如数据文件中的变异位点命名方式与数据库中的命名方式不匹配。你可以检查数据文件中的变异位点命名方式是否与数据库中的一致，或者尝试使用其他的注释软件来进行比对，例如 VEP（Variant Effect Predictor）。

最后，确保你的输入数据中所列的 CHR 和 POS 是正确的，以免导致注释不准确。如果确定输入数据正确无误，你可以尝试使用其他注释工具或数据库来进行匹配，以确定是否存在问题。

总而言之，通过检查文件格式、使用适当的数据库版本、确认变异位点命名方式的一致性以及尝试其他注释工具，你应该能够解决无法匹配 rsid 的问题。如果仍然遇到困难，可以提供更多详细的信息，以便我能够提供更具体的帮助。

发表于 2023-9-24 17:27:43

根据你的描述，你使用ANNOVAR对从网上下载的GWAS摘要数据进行注释时遇到一些问题。你提到已经尝试使用了avsnp150、avsnp147和snp138数据库进行注释，但发现成功匹配的只有插入/删除（ins/del）变异，而没有位点突变（例如C/T）。你想知道是文件格式有问题还是软件的问题，或者可能是其他步骤出错了。

首先，请确保你使用的GWAS摘要数据文件格式正确，并且含有准确的染色体位置信息（CHR:POS）。通常来说，GWAS摘要数据文件中会包含每个位点（SNP）的染色体编号、位置信息以及其他相关的GWAS统计信息。

然后，针对你遇到的问题，值得注意的是，AVSAnno和ANNOVAR所使用的数据库中的标识符通常是基于rsID。因此，如果你的GWAS摘要数据中没有rsID，并且通过注释时只能成功匹配到插入/删除变异，可能是因为无法直接将CHR:POS转化为rsID。

解决这个问题的一种方法是使用NCBI的dbSNP数据库进行注释，该数据库包含了大量的SNP标识符（包括rsID）。你可以尝试下载最新版本的dbSNP数据库，并将其与ANNOVAR一起使用，以便更好地匹配和注释你的GWAS摘要数据。

此外，还有一些其他策略可以尝试。例如，你可以使用NCBI的dbSNP网站或其他在线工具，根据CHR:POS信息进行查询，以获取相应的rsID。然后，将这些rsID与ANNOVAR一起使用，进行后续的注释分析。

最后，如果以上方法仍然无法解决问题，建议参考ANNOVAR的官方文档、用户手册或在线帮助文档，查找关于该软件版本和数据库使用的详细信息。你也可以参考其他研究论文或在相关的生物信息学社区寻求帮助，以获得更多专业的建议和指导。

[已完结] 使用ANNOVAR匹配 variant rsid

已采纳

发表回复

悬赏问答

最新回复

站点公告

今日热心研友

（热心研友QQ群）

[已完结] 使用ANNOVAR匹配 variant rsid 一键复制

已采纳