跳到目录

利用R进行基因组数据检索

动机:

这个软件包诞生于我自己的挫败感,我想让基因组数据检索过程自动化,为大规模基因组研究创建可计算的可复制脚本。由于我找不到易于使用和完全可复制的软件库,我坐下来尝试实现一个框架,使任何人都能自动化和标准化基因组数据检索过程。我希望这个软件包对其他人也有用,并有助于促进基因组学研究中的可复制研究。

我很高兴欢迎任何希望为这个项目做出贡献的人:)给我发封电子邮件就行了。

请查找详细信息此处显示文档.

引用

请引用生物标记如果它对你的研究有帮助。这将允许我在未来继续维护这个项目。

Drost HG、Paszkowski J。Biomartr:利用R检索基因组数据.生物信息学(2017) 33(8): 1216-1217.doi:10.1093/bioinformatics/btw821.

短包装说明:

测序基因组数量的急剧增长使我们能够进行一种新型的生物研究。使用比较方法,这些基因组为我们提供了关于生物信息如何在分子水平上编码以及这些信息如何随进化时间变化的新见解。

然而,任何基于基因组的研究的第一步都是从数据库中检索基因组及其注释。为了在元基因组尺度上自动检索此信息生物标记该软件包为基因组序列检索和功能注释检索提供了接口功能。的主要目标生物标记是为了促进(元)基因组分析中基因组数据的计算再现性和大规模处理。此外,生物标记旨在解决基因组版本危机。使用生物标记用户现在可以控制并了解他们自动检索的基因组版本。许多大规模基因组学研究缺乏这一信息,因此,当基因组版本信息的文档被忽略时,再现性和数据解释几乎不可能实现。

具体来说,生物标记自动化基因组、蛋白质组、CDS、RNA、重复序列、GFF/GTF(注释)、基因组组装质量和从主要生物数据库(如

此外Ensemb生物集市数据库允许用户使用一种新颖的以生物体为中心的搜索策略检索基因组位点的功能注释。此外,用户可以下载整个数据库例如

  • NCBI参考序列
  • NCBI编号
  • NCBI nt公司
  • NCBI基因库
  • ENSEMBL公司

只有一个命令。

类似工作

The main difference between theBiomaRt公司包和生物标记包裹是这样的生物标记扩展功能注释检索程序BiomaRt公司此外为基因组、蛋白质组、编码序列、gff文件、RNA序列、Repeat Masker注释文件提供有用的检索功能,以及检索整个数据库的功能,例如NCBI编号等。

请咨询教程部分了解更多详细信息。

在功能注释检索的上下文中这个生物标记包允许用户仅使用感兴趣生物体的学名筛选可用的mart,而不是首先搜索支持特定感兴趣生物体(使用BiomaRt公司包装)。此外,生物标记允许您在搜索属性和筛选器时搜索特定主题。我知道,软件包的类似命名很不幸,但这是由于历史原因引起的(请在此处找到详细解释:https://github.com/ropensci/biomartr/blob/master/FASKS.md还有这里#11).

我也致力于要比较的整个小插曲这个BiomaRt公司生物标记上下文中的包功能功能注释(它们的功能重叠,仅占生物rtr包整体功能的20%左右)。

反馈

我真的很重视你的意见和改进建议。因此,如果您能参加这项1分钟3个问题的调查,我将不胜感激(https://goo.gl/forms/Qaoxjb1EnNSLpM02)这样我就可以学习如何改进生物标记以最好的方式。提前致谢。

安装

这个生物标记包依赖于一些生物导体工具,因此需要安装以下软件包:

#安装核心生物导体包
如果 (requireNamespace(必需命名空间)(“生物经理”))
    安装.包(“生物经理”)
生物技术经理::安装()
#安装包依赖项
生物技术经理::安装(“生物串”)
生物技术经理::安装(“生物反应”)

现在用户可以安装生物标记来自CRAN:

#从CRAN安装biomartr 1.0.7
安装.包(“生物标记”,依赖项= 真的)

#安装包含最新功能的开发人员版本
生物管理器::安装(“ropensci/生物标志物”)

使用Bioconda安装

使用激活的Bioconda通道(参见2。设置通道),安装时使用:

康达安装r-biomertr

并更新为:

康达更新r-biomertr

或者使用docker容器:

码头工人拉码头。io/生物容器/r-biomartr:<tag>

(检查r-生物标记对于的有效值)

例子

集合检索

自动检索集合(=基因组、蛋白质组、CDS、RNA、GFF、Repeat Masker、AssemblyStats文件)将确保生物体的基因组文件与CDS、蛋白质组,RNA、GFF等文件匹配,并使用相同的基因组组装版本生成。基因组学研究在计算和生物再现性方面失败的原因之一是,尚不清楚拟议分析中使用的CDS、蛋白质组、RNA、GFF等文件是否是使用表示相同基因组组装版本的相同基因组组装文件生成的。为了避免这个看似微不足道的错误,我们鼓励用户使用生物标记功能获取集合()并将相应的输出作为补充数据附加到相应的基因组学研究中,以确保计算和生物再现性。

#酿酒酵母下载集
生物标记::获取收藏(数据库= “参考序列”,有机体= “酿酒酵母”)

在内部获取集合()函数现在将生成一个名为refseq/收集/酵母菌_cerevisiae并将存储所有基因组和注释文件酿酒酵母在同一文件夹中。此外,准确的基因和注释版本将记录在文件文件夹。

在内部,一个名为doc_Saccharomyces_cerevisiae_db_refseq.txt文件生成。此日志文件中存储的信息结构如下:

文件名:Saccharomyces_cerevisiae_assembly_stats_refseq.txt生物体名称:酵母菌_cerevisiae数据库:NCBI参考网址:ftp://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/000/146/045/GCF_000146045.2_R64/GCF_000146045.2_R64_assembly_stats.txt下载日期:2018年6月27日星期三15:21:51refseq_category:参考基因组装配_加入:GCF_000146045.2生物项目:PRJNA128生物样本:不适用出租车:559292次特异性名称:菌株=S288Cversion_status:最新release_type:主要genome_rep:完整序号_日期:2014-12-17提交人:酵母基因组数据库

在理想情况下,该参考文件可以作为补充信息包含在任何依赖基因组信息的生命科学出版物中,从而实现实验和分析的再现性。

仅使用一个命令检索数百个基因组

下载所有哺乳动物脊椎动物基因组NCBI参考序列通过:

#下载所有脊椎动物基因组
元检索(王国= “脊椎动物_哺乳动物”,分贝= “参考序列”,类型= “基因组”)

所有基因都存储在根据王国命名的文件夹中。在这种情况下脊椎动物哺乳动物。或者,用户可以指定出文件夹参数定义自定义输出文件夹路径。

常见问题(FAQ)

请查找此处有所有常见问题解答.

讨论和错误报告

我很乐意了解更多关于此包中提供的概念和功能的潜在改进。

此外,如果您发现一些错误或需要此软件包部分的其他(更灵活的)功能,请告诉我:

https://github.com/HajkD/bioratr/issues(https://github.com/HajkD/biorartr/issues)

教程

入门生物标记:

用户还可以阅读中的教程(Posit(前RStudio)) :

#获取biomertr包
图书馆(生物标记)

#查找biomartr软件包中提供的所有教程(小插曲)
#这将打开您的web浏览器
浏览幻影(“生物标记”)

新闻

包的当前状态以及每个版本的功能的详细历史记录生物标记可以在中找到新闻部分。

安装开发人员版本

一些错误修复或新功能在CRAN上尚不可用,但在GitHub上的开发人员版本中可用。下载并安装最新版本的生物标记运行:

#在系统上安装当前版本的biometr
如果 (requireNamespace(必需命名空间)(“生物经理”,悄悄地= 真的))
    安装.包(“生物经理”)

生物技术经理::安装(“ropensci/biomartr”)

基因组数据检索

元基因组检索

基因组检索

导入下载的文件

数据库检索

BioMart查询

执行基因本体查询

基因本体论

在Windows系统上下载开发人员版本

#在Windows上,这不起作用-看到了吗?构建github_devtools
安装github(“HajkD/biomartr”,构建的小插曲= 真的,依赖项= 真的)

#使用Windows时,首先需要安装
#R包:rtools->install.packages(“rtools”)

#然后,您可以安装devtools->install.packages(“devtools”)
#然后您可以运行:

开发工具::安装github(“HajkD/biomartr”,内部版本(_V)= 真的,依赖项= 真的)

#然后从图书馆打电话
图书馆(“生物标记”,库.loc= “C:/Program Files/R/R-3.1.1/library”)

行为准则

请注意,此项目发布时带有贡献者行为准则。参与此项目即表示您同意遵守其条款。