使用Biopython.Entrez返回与基因符号列表关联的pubmed记录

时间:2014-10-27 15:56:03

标签: python search biopython pubmed

我想在pubmed数据库中搜索一个基因符号列表(下面命名为t),以便(最终)检索相关基因的DNA序列。我想限制我的搜索只限于人类,但我现在的代码给了我除了人类以外的有机体。

from Bio import Entrez
Entrez.email = '...'           #my email: always tell Entrez who you are

t = ['FOXO3']
for i in range(len(t)):
    search = 'human[orgn]'+t[i]
    handle = Entrez.esearch(db='gene',term=search)
    record = Entrez.read(handle)
    t = record[u'IdList']
    handle = Entrez.efetch('nucleotide',id=t[0],rettype='gb',retmode='text')
    print handle.read()

有人可以看到我出错的地方吗?

1 个答案:

答案 0 :(得分:0)

你弄乱了数据库。在您使用db=gene的搜索中,但在efetch中,您将其更改为db=nucleotide。它们是不同的东西: