美文网首页
制作小麦的GO注释文件

制作小麦的GO注释文件

作者: 无言_俗人 | 来源:发表于2022-12-04 12:09 被阅读0次

参考资料:

由于在GO数据库中并没有小麦的GO注释信息所以需要自己整理获得小麦中gene与GO对应关系。当然也可以利用TGT数据库网站的GOEnrichment项目进行GO富集分析。
自制GO注释文件过程具体如下:

1.下载小麦1.0版本基因注释信息表

#此处选的是高可信基因版本
wget https://urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Annotations/v1.0/iwgsc_refseqv1.0_FunctionalAnnotation_v1.zip
#解压
unzip iwgsc_refseqv1.0_FunctionalAnnotation_v1.zip

2. 编写perl脚本实现GO注释信息提取

Taes_GO_file_from_1.0HCTAB.pl内容如下:

#!/usr/bin/perl -w
use strict;
 "usage: perl $0 TAB.file  |sort |uniq > GO_file.txt\n" unless @ARGV==1;
open TAB,"$ARGV[0]" or die $!;
print "gene\tGO\tprocess\tdesp\n";
while(<TAB>){
chomp;
my @line=split/\t/;
#将那些没有注释到任何GO通路的基因过滤掉
next if $line[7] !~/GO:/;
#仅使用基因id,而不用转录本id
my $gene=(split/\./,$line[0])[0];
#将那些有多个GO注释过程的基因分割成多行,每行仅含一个过程
$line[7]=~s/;/\n$gene\t/g;
$line[7]=~s/\sMF:\s/\tMF\t/g;
$line[7]=~s/\sCC:\s/\tCC\t/g;
$line[7]=~s/\sBP:\s/\tBP\t/g;
print $gene,"\t",$line[7],"\n"; 
}
close TAB;

运行脚本获得GO_file.txt

#此处需要去除重复的行
perl Taes_GO_file_from_1.0HCTAB.pl iwgsc_refseqv1.0_FunctionalAnnotation_v1__HCgenes_v1.0.TAB |sort |uniq > GO_file.txt

GO_file.txt结果展示如下:

gene    GO      process desp
TraesCS1A01G001800      GO:0003735      MF      structural constituent of ribosome
TraesCS1A01G001800      GO:0005622      CC      intracellular
TraesCS1A01G001800      GO:0005840      CC      ribosome
TraesCS1A01G001800      GO:0006412      BP      translation
TraesCS1A01G002000      GO:0005515      MF      protein binding
TraesCS1A01G002000      GO:0008270      MF      zinc ion binding
TraesCS1A01G002400      GO:0005506      MF      iron ion binding

相关文章

  • python从基因组注释文件中提取GO信息

    问题 在转录组或基因组注释文件中有基因的GO注释,为了做GO富集分析,需要从中提取每个有GO注释的基因对应的GO号...

  • 代码规范

    文件注释 file - File and code templetes -Go File 添加快捷键 函数注释(...

  • 初识Go语言

    Go 语言结构 初识Go语言 Go 语言结构 代码说明: ① 定义包名,必须在源文件中非注释的第一行指明这个文件...

  • 10 注释

    GO 注释 KEGG注释 Go注释批量导出结果

  • 3.6. 生成代码文档

    go doc工具会从 Go 程序和包文件中提取顶级声明的首行注释以及每个对象的相关注释,并生成相关文档。 用法 g...

  • Golang学习笔记-基础知识

    包和文件名 Go源文件中非注释的第一行需要指明这个文件属于哪个包,每个 Go 应用程序都包含一个名为 main 的...

  • 02-Go语言Hello World

    Go与C语法相似 Go注释和C语言一样// 表示单行注释/**/ 表示多行注释 Go也是由函数构成 Go语言的入口...

  • awk 提取

    今日分析遇到的问题GWAS结果合并很多SNP,想在注释文件中提取某些SNP的注释信息,本来想用R分析,奈何小麦的注...

  • golang Grpc入门Demo

    1、制作中间文件 安装protoc并配置好环境变量,用来编译中间文件 安装protoc-gen-go来生成go相关...

  • Android Studio添加文件模板

    Android开发中常用的注释有文件注释,类注释,方法注释,类的成员变量注释,常量注释,xml文件注释等,而文件模...

网友评论

      本文标题:制作小麦的GO注释文件

      本文链接:https://www.haomeiwen.com/subject/nhlvfdtx.html