python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子相关的知识,希望对你有一定的参考价值。

import xmlrpclib
import time

url = "http://deepblue.mpi-inf.mpg.de/xmlrpc"
user_key = "anonymous_key"

server = xmlrpclib.Server(url, allow_none=True)

(status, exps) = server.select_regions("", "GRCh38", "H3k27ac", "", "", 
                                       "BLUEPRINT Epigenome", "chr1", None,
                                        None, user_key)
(status, exps_peaks) = server.query_experiment_type(exps, "peaks", user_key)


(status, promoters) = server.select_annotations("promoters", "GRCh38","chr1", 
                                                None, None, user_key)

(status, exps_promoters) = server.intersection(exps_peaks, promoters, user_key)

(status, tf) = server.select_regions("","hg19", ["SP1"], "", 
                                     "", "ENCODE", "chr1", None, None, user_key)
(status, ts_signals) = server.query_experiment_type(tf, "signal", user_key)

(status, final) = server.intersection(tf, exps_promoters, user_key)

(status, request) = server.get_regions(final, 
                      "CHROMOSOME,START,END,@NAME,@EPIGENETIC_MARK,@BIOSOURCE",
                      user_key)

# Wait the processing
(status, info) = server.info(request, user_key)
while info[0]["state"] != "done" and info[0]["state"] != "error":
    time.sleep(5)
    print info
    (status, info) = server.info(request, user_key)

(status, regions) = server.get_request_data(request, user_key)
print regions

以上是关于python 鉴定H3k27ac峰,其与所有BLUEPRINT数据集中的启动子重叠并随后鉴定转录因子的主要内容,如果未能解决你的问题,请参考以下文章

R语言学习 - 线图绘制

python 通过H-ESC H3K4me3峰聚集肝脏数据实验的甲基化数据信号

python 总结来自人肝细胞的H3K4me3峰区域的肝组织中的DNA甲基化水平

硬件开源电路STC8H3K64S4开发板

R语言学习 - 线图一步法

Journal of Proteomics Research | 构建用于鉴定蓖麻毒素的串联质谱库