目前基因组测序和组装成本几乎已经到任何一个课题组都可以单独负担的价码,大量物种的基因组序列被测定和释放。与此同时,对应的基因结构注释信息文件,如GTF或GFF3文件等,也可公开下载。
对于绝大多数要使用这些公共资源的研究人员而言,有了这两个文件就足够了。但想象总是美好,现实却常常骨感。物种基因组很多,基因组序列质量不高的同时,基因结构注释信息文件更加不规范,往往直接影响下游数据分析。我最近在看一些数据,也就遇到类似的问题。如果规范化一个GFF3文件,补充对应的信息,使得该文件能够尽可能适应各类下游数据分析,更或者方便进一步分析使用。
做了大概检索,发现几乎没有同类功能的工具,除了一个perl脚本集合(github搜索 aget )或许可以做类似工作外,确实没发现相关功能的软件实现,尤其是...又还有界面。(Emmm,我还是会乐于承认功能已有实现的软件存在,当然可能是就算啥也没有,我还有 JIGplot 绘图引擎。不会像某两三个团队,似乎有意回避 TBtools,无视他人科研贡献。 :D - PS: 好彩 TBtools 文稿,尤其是预印本、github和一些早期视频讲演早早放到网上去了,不然就要被“李鬼打李逵”了)
回到主题,花了点时间,写了 GXF Fix,然后放到 TBtools 主程序,主要用于修复基因结构注释文件中的缺失部分(基于文件中已有信息)。
基因结构注释信息文件,
所有这些,都可能影响下游分析。如果有一个方便快捷的功能就好了。
我已经很久没有写出这么符合 TBtools 开发理念的功能了。我们的口号是 -IOS~:
于是,功能界面如下,通过菜单跳转。
Emmm,拟南芥的基因结构注释信息文件,可以说是非常全面....我们可以跑跑看。
内容对比如下
整体行数,完全没变化...
因为其实就没啥好修复的,这个注释信息很好。
香蕉的基因结构注释质量有待提升,不过具体内容上都还可以。
内容对比如下
整体行数变化较多,毕竟 UTR 还是有不少。
有一些注释信息文件,还是比较粗糙的。比如缺少 gene,甚至 mRNA feature都没有。
# stringtie --merge -G F:\RNA-seq\algnment\sme.gff3 -o merged.stringtie.gtf WT.cd1_1.fq_clean.fq.pos.sorted.assembly.gtf WT.cd2_1.fq_clean.fq.pos.sorted.assembly.gtf WT.unc1_1.fq_clean.fq.pos.sorted.assembly.gtf WT.unc2_1.fq_clean.fq.pos.sorted.assembly.gtf
# StringTie version 2.1.4
SMEL3Ch00.00004 StringTie transcript 42 1079 1000 - . gene_id "MSTRG.1"; transcript_id "MSTRG.1.1";
SMEL3Ch00.00004 StringTie exon 42 717 1000 - . gene_id "MSTRG.1"; transcript_id "MSTRG.1.1"; exon_number "1";
SMEL3Ch00.00004 StringTie exon 836 1079 1000 - . gene_id "MSTRG.1"; transcript_id "MSTRG.1.1"; exon_number "2";
SMEL3Ch00.00004 StringTie exon 97 210 1000 - . gene_id "MSTRG.1"; transcript_id "SMEL_000g000120.1.01"; exon_number "1";
SMEL3Ch00.00004 StringTie exon 484 717 1000 - . gene_id "MSTRG.1"; transcript_id "SMEL_000g000120.1.01"; exon_number "2";
SMEL3Ch00.00004 StringTie exon 836 1079 1000 - . gene_id "MSTRG.1"; transcript_id "SMEL_000g000120.1.01"; exon_number "3";
直接进行修复
Emmm,功能还是很不错的。我觉得这个应该会成为一个常用功能吧。尤其是,越来越多不够规范的 GTF 或者 GFF3 文件出现之后。
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。
原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。
如有侵权,请联系 cloudcommunity@tencent.com 删除。