如何优雅展示GO富集结果

作者: TOP生物信息 | 来源:发表于2021-03-17 22:10 被阅读0次

    图片好像加载不出来,可移步:https://blog.csdn.net/qq_38774801/article/details/114952565


    不知道分析过转录组的小伙伴有没有碰到过这样的问题:转录组后续做GO富集的时候,发现最前面的那些term其实对应的基因都差不多,如果选前几个画图,可能说明的只是一件事。我一般是手动选,在足够显著的前提下,选择几个不一样的term来画图(当然还得符合预期)。今天就发现了一个网页工具,可以对(语义)相似的GO term进行聚类,然后选出每个cluster中的代表term,同时不得不提该工具的图形展示,个人觉得相当好看。

    先看图:

    该工具支持大部分的模式物种:Homo sapiens (human), Gallus gallus (chicken), Bos taurus (cow), Canis familiaris (dog), Mus musculus (mouse), Rattus norvegicus (rat), Caenorhabditis elegans (nematode), Arabidopsis thaliana (thale cress), Drosophila melanogaster (fruit fly), Saccharomyces cerevisiae (budding yeast),and Danio rerio (zebrafish)

    使用很简单,只需额外提供差异分析得到的基因集,以HGNC symbols命名。

    网址在:http://funset.uno/
    用法如下:

    #我用的时候不能上传差异基因文件,所以就转换了一下,复制粘贴上去的
    $ less -S 0317.txt | awk '{ORS=",";$1=$1; print $0}' | cat
    GZMK,CCL4L2,CCL4,GZMH,CMC1,NKG7,CST7,TUBA4A,CD8B,ITM2C,CD8A
    ......
    

    整个富集、聚类、绘图过程,300个差异基因不到5分钟,在可接受范围之内
    最终结果展示是这样的

    左边图片显示区可以放大、缩小、拖拽,右上角可以保存svg图片,以及富集结果的表格

    后续可以在svg文件的基础上做一些细节的修改,如果觉得网页找的代表性term不准确,也可以自己根据表格修改。表格长这样


    整体来看,这个工具的绘图很赞,值得一试!

    因水平有限,有错误的地方,欢迎批评指正!

    相关文章

      网友评论

        本文标题:如何优雅展示GO富集结果

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