GO 分析
基因本体论（Gene Ontology, GO）是一个用于描述基因和基因产品属性的标准术语体系。它提供了一个有组织的方式来表示基因在生物体内的各种角色。基因本体论通常从三个层面对基因进行描述：细胞成分（ Cellular Component ， CC ）、生物学过程（ Biological Process ， BP ）和分子功能（ Molecular Function ， MF ）。
示例：
细胞成分（ Cellular Component ， CC ）：这个层面描述了基因产物（如蛋白质）在细胞内的定位。例如，它们可能位于细胞核、细胞质、线粒体膜或其他细胞器上。这有助于了解基因产物在细胞内的作用和功能。
生物学过程（ Biological Process ， BP ）：这个层面描述了基因参与的生物学过程。这些过程可能包括细胞生长、信号传导、基因表达调控、代谢途径等。通过了解基因参与的生物学过程，我们可以更好地理解生物体的生理功能和疾病发生机制。
分子功能（ Molecular Function ， MF ）：这个层面描述了基因产物在分子层面的功能，通常涉及到与其他分子的相互作用或催化生化反应。例如，某个基因产物可能是一个酶，它能够催化某个特定的生化反应；或者它可能是一个结构蛋白，参与细胞骨架的组装和维持。
本期将介绍使用R clusterProfiler包对基因进行GO功能富集分析、OrgDb包制作以及结果可视化。使用R clusterProfiler包对基因进行富集，需要导入目的基因（前景基因）相对应物种的参考基因组（背景基因），现阶段"Bioconductor"已有十几种常见动物，如人类、小鼠等物种的OrgDb。但仍然有许多物种不在Bioconductor的OrgDb列表里，但存在参考基因组，如山羊，绵羊等，因此大家可以选择性使用。
上篇文章我们已经用Perl脚本将蛋白名批量转化为了基因名，下面首先用R 语言将基因名找到基因 ID，以作后续分析。
1.在上篇文章转换的结果targetSymbol文件中提取基因名，如图所示，建立"基因.txt"文件
2.导入文件，运行代码，得到id.txt文件
library("org.Hs.eg.db") ##载入物种包 ##id转换 rt=read.table("基因.txt",sep="\t",check.names=F,header=F) ##读取数据 genes=as.vector(rt[,1]) entrezIDs <- mget(genes, org.Hs.egSYMBOL2EG, ifnotfound=NA) entrezIDs <- as.character(entrezIDs) out=cbind(rt,1,entrezID=entrezIDs) colnames(out)=c("symbol","logFC","entrezID") write.table(out,file="id.txt",sep="\t",quote=F,row.names=F)
3.GO分析，将上述分析得到的id.txt文件导入，运行以下代码以做GO分析
4.得到GO文件，绘制气泡图和柱状图
下篇笔记将会展示如何应用R语言作KEGG分析，敬请关注与期待。整理不易用你们发正刊（CNS）的金手指为本篇内容点赞支持，本文代码与文件在以下链接。
链接：https://pan.baidu.com/s/1CW6Trynp9DNPhyCncFBLfA
提取码：0jsj
文章转载自公众号：皮蛋笔记，欢迎关注，随时获取第一手文章内容。