【生物信息学学习】第一天：生物数据库使用

本文内容均来自山东大学生物信息学课程

生物数据库

这一篇文章主要介绍生物信息学需要用到的各个数据库

由于数据库的数量众多，本文将对数据库进行一个分类，并选取其中几个重要数据库进行分析。

核酸数据库	蛋白质数据库	专用数据库
NCBI-GenBank		PubMed
EEMBL-ENA
DDBJ

一、文献搜索（PubMed）

你的导师给你扔来了一大段基因序列让你去研究，该如何进行呢？
你头大，你去问你师兄。你的师兄告诉你，这一段序列十分的像dUTPase。
但什么是dUTPase呢？你的师兄不告诉你

这时候我们就需要借助自己的力量，使用PubMed来解决！

1. 什么是PubMed

PubMed是拥有超过两百四十万的生物医学文献。它们来源于MEDLINE（生物医学文献数据库）、生命科学领域数据杂志以及在线的专业书籍。
PubMed链接: http://www.ncbi.nlm.nih.gov/pubmed

在PubMed中进行搜索：

关于PubMed搜索的相关内容，通过WoS的知识进行理解。（后续会更新WoS搜索相关文章）

2. PubMed存在的问题

（1）搜索1995年前文献中排名是为以后的作者
（2）搜索1976年以前的文献是没有摘要的
（3）1965年前的文献较难搜索

通过PubMed我们了解到dUTPase到底是什么：
dUTP焦磷酸酶（dUTPase）是DNA合成中的一种关键酶，广泛分布于真核、原核细胞以及病毒等生命有机体内。该酶能够水解细胞质中的dUTP，从而最大限度地减少尿嘧啶在DNA合成中的错误插入，降低细胞中dUTP/dTTP的比例，维持基因组复制的保真度和顺利进行。

二、一级核酸数据库

当你对导师的任务有了一个初步的了解，之后就该去干活惹。那我们在哪里获取核酸的数据呢？

1. GenBank

GenBank是由美国国家生物技术信息中心（NCBI）开发并负责维护，隶属于美国国立卫生研究院（NIH）。
GenBank数据库: http://www.ncbi.nlm.nih.gov/

注意：在NCBI网站中需选择Nucleotide选项（GenBank）

解读GenBank

恭喜！到了这一阶段，说明你已经开始渐渐习惯打工人的生活了。我们现在开始学习你该学会的第一个技能：解读一级核酸数据库

在这一章，我们用真核生物和原核生物作为例子，来了解两者的基因在数据库中不同的存储以及注释方式

首先，我们得明白，为什么他们之间为何存在不同？

	Prokaryotes	Eukaryotes
基因组大小	0.5-91 million bp	10-670,000 million bp
基因密度	one gene / 1,000 bp	one gene / 100,000 bp
编码区含量	70%	5%
基因是否线性分布	是	否
mRNA有无内含子	无	有

原核生物核酸数据库解读

当你了解到了真核生物和原核生物的区别之后。你的导师又跑来找你，让你搜索一下大肠杆菌（原核生物）dUTPase的DNA序列：X01714

那这时候我们就可以通过上面的几个数据库，轻松地得到对应的信息，X01714的DNA序列如下:

LOCUS       X01714                  1609 bp    DNA     linear   BCT 23-OCT-2008
DEFINITION  E. coli dut gene for dUTPase (EC 3.6.1.23) (deoxyuridine
            5'-triphosphate nucleotidohydrolase).
ACCESSION   X01714
VERSION     X01714.1
KEYWORDS    dUTPase; unidentified reading frame.
SOURCE      Escherichia coli
  ORGANISM  Escherichia coli
            Bacteria; Pseudomonadota; Gammaproteobacteria; Enterobacterales;
            Enterobacteriaceae; Escherichia.
REFERENCE   1  (bases 1 to 1609)
  AUTHORS   Lundberg,L.G., Thoresson,H.O., Karlstrom,O.H. and Nyman,P.O.
  TITLE     Nucleotide sequence of the structural gene for dUTPase of
            Escherichia coli K-12
  JOURNAL   EMBO J. 2 (6), 967-971 (1983)
   PUBMED   6139280
COMMENT     Data kindly reviewed (25-NOV-1985) by L. Lundberg.
FEATURES             Location/Qualifiers
     source          1..1609
                     /organism="Escherichia coli"
                     /mol_type="genomic DNA"
                     /db_xref="taxon:562"
     regulatory      286..291
                     /regulatory_class="promoter"
                     /note="-35 region"
     regulatory      310..316
                     /regulatory_class="promoter"
                     /note="-10 region"
     misc_feature    322..324
                     /note="put. transcription start region"
     regulatory      330..333
                     /regulatory_class="ribosome_binding_site"
                     /note="put. rRNA binding site"
     CDS             343..798
                     /note="unnamed protein product; dUTP-ase (aa 1-151)"
                     /codon_start=1
                     /transl_table=11
                     /protein_id="CAA25859.1"
                     /db_xref="GOA:P06968"
                     /db_xref="InterPro:IPR008180"
                     /db_xref="InterPro:IPR008181"
                     /db_xref="PDB:1DUD"
                     /db_xref="PDB:1DUP"
                     /db_xref="PDB:1EU5"
                     /db_xref="PDB:1EUW"
                     /db_xref="PDB:1RN8"
                     /db_xref="PDB:1RNJ"
                     /db_xref="PDB:1SEH"
                     /db_xref="PDB:1SYL"
                     /db_xref="PDB:2HR6"
                     /db_xref="PDB:2HRM"
                     /db_xref="UniProtKB/Swiss-Prot:P06968"
                     /translation="MKKIDVKILDPRVGKEFPLPTYATSGSAGLDLRACLNDAVELAP
                     GDTTLVPTGLAIHIADPSLAAMMLPRSGLGHKHGIVLGNLVGLIDSDYQGQLMISVWN
                     RGQDSFTIQPGERIAQMIFVPVVQAEFNLVEDFDATDRGEGGFGHSGRQ"
     misc_feature    831..851
                     /note="put.stem-loop structure"
     repeat_region   831..838
                     /note="inverted repeat A"
     repeat_region   844..851
                     /note="inverted repeat A'"
     misc_feature    866..893
                     /note="put. stem-loop structure"
     repeat_region   866..872
                     /note="imp. inverted repeat B"
     repeat_region   888..893
                     /note="imp. inverted repeat B'"
     regulatory      889..895
                     /regulatory_class="ribosome_binding_site"
                     /note="pot. rRNA binding site"
     CDS             905..1540
                     /note="unnamed protein product; unidentified reading
                     frame"
                     /codon_start=1
                     /transl_table=11
                     /protein_id="CAA25860.1"
                     /db_xref="GOA:P0C093"
                     /db_xref="InterPro:IPR001647"
                     /db_xref="InterPro:IPR009057"
                     /db_xref="InterPro:IPR011075"
                     /db_xref="InterPro:IPR015893"
                     /db_xref="UniProtKB/Swiss-Prot:P0C093"
                     /translation="MAEKQTAKRNRREEILQSLALMLESSDGSQRITTAKLAASVGVS
                     EAALYRHFPSKTRMFDSLIEFIEDSLITRINLILKDEKDTTARLRLIVLLLLGFGERN
                     PGLTRILTGHALMFEQDRLQGRINQLFERIEAQLRQVLREKRMREGEGYTTDETLLAS
                     QILAFCEGMLSRFVRSEFKYRPTDDFDARWPLIAASCSNMTPDDFSSGEFL"
ORIGIN      
        1 cagagaaaat caaaaagcag gccacgcagg gtgatgaatt aacaataaaa atggttaaaa
       61 accccgatat cgtcgcaggc gttgccgcac taaaagacca tcgaccctac gtcgttggat
      121 ttgccgccga aacaaataat gtggaagaat acgcccggca aaaacgtatc cgtaaaaacc
      181 ttgatctgat ctgcgcgaac gatgtttccc agccaactca aggatttaac agcgacaaca
      241 acgcattaca ccttttctgg caggacggag ataaagtctt accgcttgag cgcaaagagc
      301 tccttggcca attattactc gacgagatcg tgacccgtta tgatgaaaaa aatcgacgtt
      361 aagattctgg acccgcgcgt tgggaaggaa tttccgctcc cgacttatgc cacctctggc
      421 tctgccggac ttgacctgcg tgcctgtctc aacgacgccg tagaactggc tccgggtgac
      481 actacgctgg ttccgaccgg gctggcgatt catattgccg atccttcact ggcggcaatg
      541 atgctgccgc gctccggatt gggacataag cacggtatcg tgcttggtaa cctggtagga
      601 ttgatcgatt ctgactatca gggccagttg atgatttccg tgtggaaccg tggtcaggac
      661 agcttcacca ttcaacctgg cgaacgcatc gcccagatga tttttgttcc ggtagtacag
      721 gctgaattta atctggtgga agatttcgac gccaccgacc gcggtgaagg cggctttggt
      781 cactctggtc gtcagtaaca catacgcatc cgaataacgt cataacatag ccgcaaacat
      841 ttcgtttgcg gtcatagcgt gggtgccgcc tggcaagtgc ttattttcag gggtattttg
      901 taacatggca gaaaaacaaa ctgcgaaaag gaaccgtcgc gaggaaatac ttcagtctct
      961 ggcgctgatg ctggaatcca gcgatggaag ccaacgtatc acgacggcaa aactggccgc
     1021 ctctgtcggc gtttccgaag cggcactgta tcgccacttc cccagtaaga cccgcatgtt
     1081 cgatagcctg attgagttta tcgaagatag cctgattact cgcatcaacc tgattctgaa
     1141 agatgagaaa gacaccacag cgcgcctgcg tctgattgtg ttgctgcttc tcggttttgg
     1201 tgagcgtaat cctggcctga cccgcatcct cactggtcat gcgctaatgt ttgaacagga
     1261 tcgcctgcaa gggcgcatca accagctgtt cgagcgtatt gaagcgcagc tgcgccaggt
     1321 attgcgtgaa aagagaatgc gtgagggtga aggttacacc accgatgaaa ccctgctggc
     1381 aagccagatc ctggccttct gtgaaggtat gctgtcacgt tttgtccgca gcgaatttaa
     1441 ataccgcccg acggatgatt ttgacgcccg ctggccgcta attgcggcca gttgcagtaa
     1501 tatgacgccg gatgactttt catccggcga gtttctttaa acgccaaact cttcgcgata
     1561 ggccttaacc gccgccagat gttccgccat ttccggcttc tcttccagg
//

观察数据的第一行，是不是一时有些懵？这一段代表的是什么意思呢？用中文标记一下，或许更有利于理解：

LOCUS	X01714	1609bp	DNA	linear	BCT 23-0CT-2000
	基因座名	核酸序列长度	分子的类别	拓扑类型	更新时间

以下是剩余的几行相关术语的注解：

ACCESSION：检索号在数据库中是唯一旦不变的，即使数据提交者改变数据内容。ACCESSION不一定和IOCUS相同。（原因：LOCUS是真实姓名，ACCESSION是编号。同一个基因只有一个名字，但可以在不同的数据库中有不同的编号。）

VERSION：版本号的格式是“检索号.版本编号”。版本号于1999年2月由三大数据库采纳使用。主要用于识别数据库中一条单一的特定核苷酸序列。在数据库中，如果某条序列数据发生了变化，即使是单碱基的改变，它的版本号都将增加，而它的检索号保持不变。
（例如：由U12345.1变为U12345.2。）

GI (Genlnfo Identifier）号：与前面的版本号系统是平行运行的。当一条序列改变后，它将被赋予一个新的GI号，同时它的版本号將增加。

KEYWORDS:能够大致描述该条目的几个关键词。

SOURCE：基因序列所属物种的俗名。

ORGANISM：对所属物种的更详细定义，包括他的科学分类。

REFERENCE：基因序列来源的科学文献（一条基因序列的不同片段可能来源于不同的文献）。文献具体分为作者、题目和刊物。刊物还包括PubMedID作为其子条目。

COMMENT：自由撰写内容，比如致谢或者无法归入前几类的内容。

FEATURES（重要信息）： 描迷核酸序列中各个已确定的片段区域，包含很多子条目，比如来源(source)，启动子 （ptomotet）等。
source：说明了核酸序列的来源，据此可以容易地分辦出该序列是来源于克隆載体还是基因组。当前序列（全长）来源于大肠杆茵的基因组DNA。
promotet：列出了启动子的位置。复习：细菌有两个启
动子区，一个-35区（5’-TTGACA-3〞）位置在第286
个碱基到第291个碱基，一个-10区（5’-TATAAT-3〞）
位置在第310个碱基到第316个碱基。
misc_ featute：混合内容。比如，这条说明了从第322个
碱基到第324个碱基是一个推测的(putative，推定，但无实验证实） 转录起始位置。
RBS (Ribosome Binding Site)：核糖体结合位点。
CDs(CodingSegment)：记录了一个ORF(open reading frame），从第343个碱基开始的ATG（起始密码子）到第798个碱基结束的TAA（结束密码子）。除了第一行的位置信息，还包括翻译产物（蛋白质）的诸多信息。

  CDS             343..798
                     /note="unnamed protein product; dUTP-ase (aa 1-151)"
                     /*翻译产物蛋白的名字:该ORF编码其1-151个氨基酸*/
                     /codon_start=1
                     /transl_table=11
                     /*翻译起始位置和使用的密码本*/
                     /protein_id="CAA25859.1"
                     /db_xref="GOA:P06968"
                     /db_xref="InterPro:IPR008180"
                     /db_xref="InterPro:IPR008181"
                     /db_xref="PDB:1DUD"
                     /db_xref="PDB:1DUP"
                     /db_xref="PDB:1EU5"
                     /db_xref="PDB:1EUW"
                     /db_xref="PDB:1RN8"
                     /db_xref="PDB:1RNJ"
                     /db_xref="PDB:1SEH"
                     /db_xref="PDB:1SYL"
                     /db_xref="PDB:2HR6"
                     /db_xref="PDB:2HRM"
                     /db_xref="UniProtKB/Swiss-Prot:P06968"
                     /*中间一大段是该蛋白质序列在各蛋白质数据库中的检索号*/
                     /translation="MKKIDVKILDPRVGKEFPLPTYATSGSAGLDLRACLNDAVELAP
                     GDTTLVPTGLAIHIADPSLAAMMLPRSGLGHKHGIVLGNLVGLIDSDYQGQLMISVWN
                     RGQDSFTIQPGERIAQMIFVPVVQAEFNLVEDFDATDRGEGGFGHSGRQ"               		
                     /*计算机使用翻译密码本根据核酸序列翻译出的蛋白质序列（不是实际获得）*/

/*该核酸序列还存在着潜在基因unidentified reading frame，是由计算机预测出的基因，编码蛋白在系统中无明确记录*/
  CDS             905..1540
                     /note="unnamed protein product; unidentified reading
                     frame"
                     /codon_start=1
                     /transl_table=11
                     /protein_id="CAA25860.1"
                     /db_xref="GOA:P0C093"
                     /db_xref="InterPro:IPR001647"
                     /db_xref="InterPro:IPR009057"
                     /db_xref="InterPro:IPR011075"
                     /db_xref="InterPro:IPR015893"
                     /db_xref="UniProtKB/Swiss-Prot:P0C093"
                     /translation="MAEKQTAKRNRREEILQSLALMLESSDGSQRITTAKLAASVGVS
                     EAALYRHFPSKTRMFDSLIEFIEDSLITRINLILKDEKDTTARLRLIVLLLLGFGERN
                     PGLTRILTGHALMFEQDRLQGRINQLFERIEAQLRQVLREKRMREGEGYTTDETLLAS
                     QILAFCEGMLSRFVRSEFKYRPTDDFDARWPLIAASCSNMTPDDFSSGEFL"

/*核酸序列*/
ORIGIN      
        1 cagagaaaat caaaaagcag gccacgcagg gtgatgaatt aacaataaaa atggttaaaa
       61 accccgatat cgtcgcaggc gttgccgcac taaaagacca tcgaccctac gtcgttggat
      121 ttgccgccga aacaaataat gtggaagaat acgcccggca aaaacgtatc cgtaaaaacc
      181 ttgatctgat ctgcgcgaac gatgtttccc agccaactca aggatttaac agcgacaaca
      241 acgcattaca ccttttctgg caggacggag ataaagtctt accgcttgag cgcaaagagc
      301 tccttggcca attattactc gacgagatcg tgacccgtta tgatgaaaaa aatcgacgtt
      361 aagattctgg acccgcgcgt tgggaaggaa tttccgctcc cgacttatgc cacctctggc
      421 tctgccggac ttgacctgcg tgcctgtctc aacgacgccg tagaactggc tccgggtgac
      481 actacgctgg ttccgaccgg gctggcgatt catattgccg atccttcact ggcggcaatg
      541 atgctgccgc gctccggatt gggacataag cacggtatcg tgcttggtaa cctggtagga
      601 ttgatcgatt ctgactatca gggccagttg atgatttccg tgtggaaccg tggtcaggac
      661 agcttcacca ttcaacctgg cgaacgcatc gcccagatga tttttgttcc ggtagtacag
      721 gctgaattta atctggtgga agatttcgac gccaccgacc gcggtgaagg cggctttggt
      781 cactctggtc gtcagtaaca catacgcatc cgaataacgt cataacatag ccgcaaacat
      841 ttcgtttgcg gtcatagcgt gggtgccgcc tggcaagtgc ttattttcag gggtattttg
      901 taacatggca gaaaaacaaa ctgcgaaaag gaaccgtcgc gaggaaatac ttcagtctct
      961 ggcgctgatg ctggaatcca gcgatggaag ccaacgtatc acgacggcaa aactggccgc
     1021 ctctgtcggc gtttccgaag cggcactgta tcgccacttc cccagtaaga cccgcatgtt
     1081 cgatagcctg attgagttta tcgaagatag cctgattact cgcatcaacc tgattctgaa
     1141 agatgagaaa gacaccacag cgcgcctgcg tctgattgtg ttgctgcttc tcggttttgg
     1201 tgagcgtaat cctggcctga cccgcatcct cactggtcat gcgctaatgt ttgaacagga
     1261 tcgcctgcaa gggcgcatca accagctgtt cgagcgtatt gaagcgcagc tgcgccaggt
     1321 attgcgtgaa aagagaatgc gtgagggtga aggttacacc accgatgaaa ccctgctggc
     1381 aagccagatc ctggccttct gtgaaggtat gctgtcacgt tttgtccgca gcgaatttaa
     1441 ataccgcccg acggatgatt ttgacgcccg ctggccgcta attgcggcca gttgcagtaa
     1501 tatgacgccg gatgactttt catccggcga gtttctttaa acgccaaact cttcgcgata
     1561 ggccttaacc gccgccagat gttccgccat ttccggcttc tcttccagg

终于，我们把X01714的信息给解读完了，但此时还不能休息

我们还需要学会如何下载其信息，我们直接采用ctrlcv大法是不行的。我们把界面翻到最上面，可以看到两个选项：FASTA和Graphics。

点击FASTA我们可以获取FASTA格式的核酸序列信息，点击Graphics可以获取序列的图形概览。

真核生物核酸数据库解读

你的导师又派给了你一个任务，去研究一下人（真核生物）dUTPase的成熟mRNA（U90223）序列信息。

根据上文所述的技巧，我们可以得到真核生物的序列信息：

LOCUS       HSU90223                 960 bp    mRNA    linear   PRI 03-JAN-1998
DEFINITION  Human deoxyuridine triphosphate nucleotidohydrolase precursor mRNA,
            nuclear gene encoding mitochondrial protein, complete cds.
ACCESSION   U90223
VERSION     U90223.1
KEYWORDS    .
SOURCE      Homo sapiens (human)
  ORGANISM  Homo sapiens
            Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi;
            Mammalia; Eutheria; Euarchontoglires; Primates; Haplorrhini;
            Catarrhini; Hominidae; Homo.
REFERENCE   1  (bases 1 to 960)
  AUTHORS   Ladner,R.D. and Caradonna,S.J.
  TITLE     The Human dUTPase Gene Encodes Both Nuclear and Mitochondrial
            Isoforms: Differential Expression of the Isoforms and
            Characterization of a cDNA Encoding the Mitochondrial Species
  JOURNAL   Unpublished
REFERENCE   2  (bases 1 to 960)
  AUTHORS   Ladner,R.D. and Caradonna,S.J.
  TITLE     Direct Submission
  JOURNAL   Submitted (19-FEB-1997) Dept. of Molecular Biology, Univ. of Med.
            and Dent. of NJ-School of Osteopathic Medicine, 2 Medical Center
            Drive, Stratford, NJ 08084, USA
FEATURES             Location/Qualifiers
     source          1..960
                     /organism="Homo sapiens"
                     /mol_type="mRNA"
                     /db_xref="taxon:9606"
     CDS             63..821
                     /note="mitochondrial dUTPase isoform; DUT-M"
                     /codon_start=1
                     /product="deoxyuridine triphosphate nucleotidohydrolase
                     precursor"
                     /protein_id="AAB94642.1"
                     /translation="MTPLCPRPALCYHFLTSLLRSAMQNARGTAEGRSRGTLRARPAP
                     RPPAAQHGIPRPLSSAGRLSQGCRGASTVGAAGWKGELPKAGGSPAPGPETPAISPSK
                     RARPAEVGGMQLRFARLSEHATAPTRGSARAAGYDLYSAYDYTIPPMEKAVVKTDIQI
                     ALPSGCYGRVAPRSGLAAKHFIDVGAGVIDEDYRGNVGVVLFNFGKEKFEVKKGDRIA
                     QLICERIFYPEIEEVQALDDTERGSGGFGSTGKN"
     sig_peptide     63..269
                     /note="mitochondrial targeting presequence"
     mat_peptide     270..818
                     /product="deoxyuridine triphosphate nucleotidohydrolase"
ORIGIN      
        1 ggtggaagcc tggcgcacgt ccggaggtgc cgaggaccca accagcccaa actctggggg
       61 aaatgactcc cctctgccct cgccccgcgc tctgctacca tttccttacg tctctgcttc
      121 gctcagcgat gcaaaacgcg cgaggcacgg cagagggccg aagccgcggt actctccggg
      181 ccaggcccgc ccctcggccg ccggcggcgc agcacgggat tccccggccg ctgtccagcg
      241 ctggccgcct gagccaaggc tgccgcggag ccagtacagt cggggccgct ggctggaagg
      301 gcgagcttcc taaggcgggg ggaagcccgg cgccggggcc ggagacaccc gccatttcac
      361 ccagtaagcg ggcccggcct gcggaggtgg gcggcatgca gctccgcttt gcccggctct
      421 ccgagcacgc cacggccccc acccggggct ccgcgcgcgc cgcgggctac gacctgtaca
      481 gtgcctatga ttacacaata ccacctatgg agaaagctgt tgtgaaaacg gacattcaga
      541 tagcgctccc ttctgggtgt tatggaagag tggctccacg gtcaggcttg gctgcaaaac
      601 actttattga tgtaggagct ggtgtcatag atgaagatta tagaggaaat gttggtgttg
      661 tactgtttaa ttttggcaaa gaaaagtttg aagtcaaaaa aggtgatcga attgcacagc
      721 tcatttgcga acggattttt tatccagaaa tagaagaagt tcaagccttg gatgacaccg
      781 aaaggggttc aggaggtttt ggttccactg gaaagaatta aaatttatgc caagaacaga
      841 aaacaagaag tcataccttt ttcttaaaaa aaaaaaaagt ttttgcttca agtgttttgg
      901 tgttttgcac ttctgtaaac ttactagctt taccttctaa aagtactgca ttttttactt
//

观察这一大串符号我们可以发现，真核生物和原核生物在数据库中的数据排列方式很类似，但是却出现了两个不一样的条目：

  sig_peptide     63..269
                     /note="mitochondrial targeting presequence"
  /*该段指出了编码用于亚细胞定位的信号肽的碱基位置*/
  mat_peptide     270..818
                     /product="deoxyuridine triphosphate nucleotidohydrolase"
  /*该段指出了编码成熟蛋白的碱基位置 */

mat_peptide 270…818的编码位置与 CDS 343…798相差了三个碱基，为什么？
回答：在这一段序列最后的三个氨基酸为编码区的终止密码子（不翻译）

新任务！：研究一下人（真核生物）dUTPase的基因组DNA序列信息（AH005568）

LOCUS       AH005568                5507 bp    DNA     linear   PRI 10-JUN-2016
DEFINITION  Homo sapiens dUTPase (DUT) gene, complete cds, alternatively
            spliced.
ACCESSION   AH005568 AF018429 AF018430 AF018431 AF018432
VERSION     AH005568.2
KEYWORDS    .
SOURCE      Homo sapiens (human)
  ORGANISM  Homo sapiens
            Eukaryota; Metazoa; Chordata; Craniata; Vertebrata; Euteleostomi;
            Mammalia; Eutheria; Euarchontoglires; Primates; Haplorrhini;
            Catarrhini; Hominidae; Homo.
REFERENCE   1  (bases 1 to 5507)
  AUTHORS   Pearlman,R.E.
  TITLE     Human genomic nuclear and mitochondria dUTPase gene
  JOURNAL   Unpublished
REFERENCE   2  (bases 1 to 5507)
  AUTHORS   Pearlman,R.E.
  TITLE     Direct Submission
  JOURNAL   Submitted (11-AUG-1997) Biology, York University, 4700 Keele St.,
            North York, ONT M3J 1P3, Canada
COMMENT     On or before Jun 10, 2016 this sequence version replaced
            AF018429.1, AF018430.1, AF018431.1, AF018432.1, AH005568.1.
FEATURES             Location/Qualifiers
     source          1..5507
                     /organism="Homo sapiens"
                     /mol_type="genomic DNA"
                     /db_xref="taxon:9606"
                     /map="15q15-q21.1"
     gene            <1..>5236
                     /gene="DUT"
     mRNA            join(<282..561,1034..1172,2395..2486,3113..3157,
                     4447..4521,4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes mitochondrial form
                     of the protein"
     CDS             join(282..561,1034..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..5236)
                     /gene="DUT"
                     /note="DUT-M; alternatively spliced; mitochondrial form of
                     the protein; similar to H. sapiens dUTPase encoded by
                     GenBank Accession Number U90224"
                     /codon_start=1
                     /product="dUTPase"
                     /protein_id="AAB71393.1"
                     /translation="MTPLCPRPALCYHFLTSLLRSAMQNARGTAEGRSRGTLRARPAP
                     RPPAAQHGIPRPLSSAGRLSQGCRGASTVGAAGWKGELPKAGGSPAPGPETPAISPSK
                     RARPAEVGGMQLRFARLSEHATAPTRGSARAAGYDLYSAYDYTIPPMEKAVVKTDIQI
                     ALPSGCYGRVAPRSGLAAKHFIDVGAGVIDEDYRGNVGVVLFNFGKEKFEVKKGDRIA
                     QLICERIFYPEIEEVQALDDTERGSGGFGSTGKN"
     exon            <282..561
                     /gene="DUT"
                     /note="exon used only in the mitochondrial form of the
                     protein"
                     /number=1
     protein_bind    898..903
                     /gene="DUT"
                     /bound_moiety="SP-1"
     protein_bind    935..943
                     /gene="DUT"
                     /bound_moiety="E2F"
     mRNA            join(<1018..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes nuclear form of the
                     protein"
     CDS             join(1018..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..5236)
                     /gene="DUT"
                     /note="DUT-N; alternatively spliced; nuclear form of the
                     protein; similar to H. sapiens dUTPase encoded by GenBank
                     Accession Number U90224"
                     /codon_start=1
                     /product="dUTPase"
                     /protein_id="AAB71394.1"
                     /translation="MPCSEETPAISPSKRARPAEVGGMQLRFARLSEHATAPTRGSAR
                     AAGYDLYSAYDYTIPPMEKAVVKTDIQIALPSGCYGRVAPRSGLAAKHFIDVGAGVID
                     EDYRGNVGVVLFNFGKEKFEVKKGDRIAQLICERIFYPEIEEVQALDDTERGSGGFGS
                     TGKN"
     exon            <1018..1172
                     /gene="DUT"
                     /note="alternative exon used in nuclear form of the
                     protein"
                     /number=2
     exon            1034..1172
                     /gene="DUT"
                     /note="alternative exon used in mitochondrial form of the
                     protein"
                     /number=2
     gap             1736..1835
                     /estimated_length=unknown
     exon            2395..2486
                     /gene="DUT"
                     /number=3
     gap             3013..3112
                     /estimated_length=unknown
     exon            3113..3157
                     /gene="DUT"
                     /number=4
     gap             3690..3789
                     /estimated_length=unknown
     exon            4447..4521
                     /gene="DUT"
                     /number=5
     exon            4673..4743
                     /gene="DUT"
                     /number=6
     exon            5180..5236
                     /gene="DUT"
                     /number=7
ORIGIN      
        1 tcggaaaaat gggggccaga gcaaacaaga agagcgaaag caagagggct aggcagccag
       61 aggcggcagc aagactcaag acgccaacgg cgctgtttcc tggggcccca gggcttgcgc
      121 catccctggg ctgccggggc accgcctctc cacgcccctc gtccggcggc ggctgcgact
      181 gcttccgagg tcatgttccc aggacgggcg cgtcttcagg gtggaagcct ggcgcacgtc
      241 cggaggtgcc gaggacccaa ccagcccaaa ctctggggga aatgactccc ctctgccctc
      301 gccccgcgct ctgctaccat ttccttacgt ctctgcttcg ctcagcgatg caaaacgcgc
      361 gaggcacggc agagggccga agccgcggta ctctccgggc caggcccgcc cctcggccgc
      421 cggccgcgca gcacgggatt ccccggccgc tgtccagcgc tggccgcctg agccaaggct
      481 gccgcggagc cagtacagtc ggggccgctg gctggaaggg cgagcttcct aaggcggggg
      541 gaagcccggc gccggggccg ggtaggaaag gcgggggagg ggctccggcc gtctggaagg
      601 aatccaacgc ggcttgaggc tgtggggaag tagggtggcg agcggtcctt ctgcgcgcgg
      661 ggggccgggg gggtggggtg gtccattagg gtcccctggc gagggggcgg ctttctagtg
      721 tgtgaaggcg acgccctaga agctcccctt caaagttggc cccacgcgct gaatgtggaa
      781 agttgactgg gacccagtag tttcccatcc caaacctgct ttccgagaag ggcttcaaac
      841 ccaaaatgtg aatcccgcct cccctctcag ccagaactgt ggactcgtcc cggggagggg
      901 cggtgggtgg ggcggggctg gcgggaaatt tcggttttgg cgcgctccct gcggcgacgc
      961 tcatcgtgcg ctctcctctt cccccggtgg tctcctcgct cgccttctgg ctctgccatg
     1021 ccctgctctg aagagacacc cgccatttca cccagtaagc gggcccggcc tgcggaggtg
     1081 ggcggcatgc agctccgctt tgcccggctc tccgagcacg ccacggcccc cacccggggc
     1141 tccgcgcgcg ccgcgggcta cgacctgtac aggtgagcgg ggacctgccg gcgaggaggc
     1201 tgggaagggc cggccgtccg ctgccacagc tagaaacagt caccggagag atcacaggaa
     1261 cacactagct ataaatagga tttctgcctt tttcgtgttt aaaattttag ctttcatctt
     1321 tggcataaat taaatagaga tttgggcaaa gactgcagaa taagtaaaat agctatacgg
     1381 tgtctagcaa ggcgttactt tgcaacgttt attgtgccct tcctaaatag aagatagaga
     1441 ggaaggccca tggtggcttt cgaagtggcc cgagggtgat gctgtgctca atagaaaaac
     1501 caaggtgaga gcctagatgt gagcgtgaaa atacctaaga aggatgaacg aagatgcatc
     1561 tgccttaaaa agttatttct atacattcat ccggcccagg gcggaatttg agaagcatct
     1621 gaaaacgaag gcagactgcc tgtatctacc acactttcat ctctacagca cgttttactg
     1681 tactaaaact tcccgtatgc tgttgtatag tcctccacac atccctaact agata     
          [gap 100 bp]    Expand Ns
     1836                                       tccct aaatcaacac agatcatgtg
     1861 gaggaataaa atggggttaa tatatgtaaa accaattagg aaactgtttc tggggcaaca
     1921 cagtaaaggg cttattcaat ggataggcta gtattattag ttagtaattg ggcccttttt
     1981 ttctttgttt cttttcttca tttttttcct tttcaaacta tgggttgtaa agcatccacc
     2041 ttttgaaagt ttgcctttct gccctttcac gctgataagt acctcagttt ccaataaact
     2101 tttgttcagg ggcaaacatt tacaatgttg acatctcttc acaccaccaa aaatattcat
     2161 ggagaattat tttatctaaa gctgtctttt taataataaa atagccacct ctaccttctt
     2221 cataaacttt taagatgaat tggtaattca tcatagcaag gttgatttta gaaactaaag
     2281 ttgcattaat tcattaaata cactgaaagt aattttgtat gcttggtcac aaagaaaata
     2341 taaaaacaat tttataaata gatttgcagt tattttcttt caatattttc ttagtgccta
     2401 tgattacaca ataccaccta tggagaaagc tgttgtgaaa acggacattc agatagcgct
     2461 cccttctggg tgttatggaa gagtgggtaa gtcatttaag aaacaggtaa ctatttgtca
     2521 agttctcctt tgtgatagat tcttcatgtt tcatttgggg taataagcag gcaatattgc
     2581 ttgggctgtg tcctaaaaga agcaccattt gtgatagcaa atgcactctt tgaaaggctt
     2641 tatttacatc tctgctttgc ctctttttga cccttttatt tttctccttc ctcactggag
     2701 cttttaggct cacactggcc tagaaggctg ttctcagaac atggcatttt atattatgag
     2761 agtaaaactt ctgacctgtt ggtcccagaa tgtgtaagcc tacttaacct tttcttgttt
     2821 ggccatgggg tttagggtaa gggatactct tcagtgtttg tagaggcact gggaggaagc
     2881 taggacaaaa tggagttaca cgtcaacagg tttgattttt cctggaagcg aattcagtgt
     2941 ttaccagaca gttcctttgc agagcgttag ttcctttttg actacttcca agttaactta
     3001 aggaggcatg ga                                                    
          [gap 100 bp]    Expand Ns
     3113                                                          ctccacgg
     3121 tcaggcttgg ctgcaaaaca ctttattgat gtaggaggta atatatttcc ttttttattc
     3181 tgtaaatgtt tgcaagtatt tactttgtct ttaaaaggta atattcaaat gacagatttt
     3241 atttttaaga aaagaaaatg attagaggaa agcttgttat aataggagga aaagctttgt
     3301 ggttatttta agtaattata ttttgtttaa ctactaactt aatttttaaa ataatgatta
     3361 cctaattaat aatgacagat ttaatgatat aataattaag taattattta agcaacaaaa
     3421 gttctagaaa acacaggtct atcttgtaag aattttcttc ataagagagc tatctgtctt
     3481 tagcaaaatt taagaagaaa atctttccct gaagtcagga gagaaaacac ttttctctaa
     3541 atttgttaag acagtccctc tgtgttacct gtttcccctt tataataaag cctttctcat
     3601 gtttggtcct caaagttaac ccctacctcc tgttggtacc ctttcatgaa gaaccctcag
     3661 tgtccccccg gccaactaac ccctgatac                                  
          [gap 100 bp]    Expand Ns
     3790          t cccccgggtc atattctaac agtatataca caaaaagtgt actcaggaag
     3841 cctatgtgct acagtctctg agaaccatgt ctaaggaagg atggtggttc agattaaaag
     3901 gagaggaaga ttatgaactg aaaaatacag ctcattttta gtgtaaaaat ggaatattct
     3961 aatatttgca agactgctgc ttaacaatgt cagcgttgat ttactccact gcagagcttt
     4021 aggctttttg acatttttat ttatgactca attgtatttc aatctgggta ttagcgggaa
     4081 gctgatacca ggaaaaacct ttaaaaccca aattgtaatt ttctcctagt ctttttggca
     4141 agtgtgccgt gggtagacag taggttggag atggccattg ctactgccag cttaaaatag
     4201 gaaatggtcg ccatttgata aagagtgaca ggttggaaaa ttaggacctg ttttctaaat
     4261 ttttatttgt attttagtaa tgtcatttct tggttttgag gctttacata gtaccagggc
     4321 aagattcact tattcaattg taagaataca aatctcagag ctacatgata atgattttga
     4381 gataatctta cacctctagt attccaatag caaaaattga gataatatta cttttctttt
     4441 ctctagctgg tgtcatagat gaagattata gaggaaatgt tggtgttgta ctgtttaatt
     4501 ttggcaaaga aaagtttgaa ggtatgttaa atatatacat tcacataatt ttagtgaatt
     4561 ttcagagtca tgtatgtgta aattaatatt gactccttta attctcattg aataagacag
     4621 gatatggcga atgtgtcagt aacgtcagta ataaactatt ctttctttga agtcaaaaaa
     4681 ggtgatcgaa ttgcacagct catttgcgaa cggatttttt atccagaaat agaagaagtt
     4741 caagtaagta ttacaaagga agatacagaa taagtaatat aacatcttaa gtgaagaaat
     4801 atatataatc ttgagaattt aatatgctgt ttgtaactaa atagtatata tgactaaact
     4861 tattttaagc aaatttaaaa tactagtttt agaatttctt taaatgtttt tcatgtagct
     4921 attatgtagt attactttgg ataataagtt atttaaacat actgtgaact tctaatattt
     4981 tattagaatt ttttaaagtt atccagtatt ctaatttatg gagcttttta gaatttaatt
     5041 ttctttctgt aatctccctt ttgaaaagat gatatagcaa gagtagaatt ctggctgtat
     5101 ttttcttagg agctggagag gaaaactgaa agatgctctt aaaaaaaact gtgaagctta
     5161 ctacctttct atctttcagg ccttggatga caccgaaagg ggttcaggag gttttggttc
     5221 cactggaaag aattaaaatt tatgccaaga acagaaaaca agaagtcata cctttttctt
     5281 aaaaaaaaaa aaaaagtttt tgcttcaagt gttttggtgt tttgcacttc tgtaaactta
     5341 ctagctttac cttctaaaag tactgcattt tttacttttt tttatgatca aggaaaagat
     5401 cattaaaaaa aaacacaaag aagtttttcc tttgtgttgg gatcaaaaga aactttgttt
     5461 tcccgcaatt gaagggttgt atgtaaatcg ccttgttgtg aaccaca
//

我们观察FEATURES段：

FEATURES             Location/Qualifiers
     source          1..5507
                     /organism="Homo sapiens"
                     /mol_type="genomic DNA"
                     /db_xref="taxon:9606"
                     /map="15q15-q21.1"
     gene            <1..>5236
                     /gene="DUT"
     mRNA            join(<282..561,1034..1172,2395..2486,3113..3157,
                     4447..4521,4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes mitochondrial form
                     of the protein"
     CDS             join(282..561,1034..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..5236)
                     /gene="DUT"
                     /note="DUT-M; alternatively spliced; mitochondrial form of
                     the protein; similar to H. sapiens dUTPase encoded by
                     GenBank Accession Number U90224"
                     /codon_start=1
                     /product="dUTPase"
                     /protein_id="AAB71393.1"
                     /translation="MTPLCPRPALCYHFLTSLLRSAMQNARGTAEGRSRGTLRARPAP
                     RPPAAQHGIPRPLSSAGRLSQGCRGASTVGAAGWKGELPKAGGSPAPGPETPAISPSK
                     RARPAEVGGMQLRFARLSEHATAPTRGSARAAGYDLYSAYDYTIPPMEKAVVKTDIQI
                     ALPSGCYGRVAPRSGLAAKHFIDVGAGVIDEDYRGNVGVVLFNFGKEKFEVKKGDRIA
                     QLICERIFYPEIEEVQALDDTERGSGGFGSTGKN"
     exon            <282..561
                     /gene="DUT"
                     /note="exon used only in the mitochondrial form of the
                     protein"
                     /number=1
     protein_bind    898..903
                     /gene="DUT"
                     /bound_moiety="SP-1"
     protein_bind    935..943
                     /gene="DUT"
                     /bound_moiety="E2F"
     mRNA            join(<1018..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes nuclear form of the
                     protein"
     CDS             join(1018..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..5236)
                     /gene="DUT"
                     /note="DUT-N; alternatively spliced; nuclear form of the
                     protein; similar to H. sapiens dUTPase encoded by GenBank
                     Accession Number U90224"
                     /codon_start=1
                     /product="dUTPase"
                     /protein_id="AAB71394.1"
                     /translation="MPCSEETPAISPSKRARPAEVGGMQLRFARLSEHATAPTRGSAR
                     AAGYDLYSAYDYTIPPMEKAVVKTDIQIALPSGCYGRVAPRSGLAAKHFIDVGAGVID
                     EDYRGNVGVVLFNFGKEKFEVKKGDRIAQLICERIFYPEIEEVQALDDTERGSGGFGS
                     TGKN"
     exon            <1018..1172
                     /gene="DUT"
                     /note="alternative exon used in nuclear form of the
                     protein"
                     /number=2
     exon            1034..1172
                     /gene="DUT"
                     /note="alternative exon used in mitochondrial form of the
                     protein"
                     /number=2
     gap             1736..1835
                     /estimated_length=unknown
     exon            2395..2486
                     /gene="DUT"
                     /number=3
     gap             3013..3112
                     /estimated_length=unknown
     exon            3113..3157
                     /gene="DUT"
                     /number=4
     gap             3690..3789
                     /estimated_length=unknown
     exon            4447..4521
                     /gene="DUT"
                     /number=5
     exon            4673..4743
                     /gene="DUT"
                     /number=6
     exon            5180..5236
                     /gene="DUT"
                     /number=7

第五行的 /map=“15q15-q21.1”，指出了该序列属于第15号染色体，更精确地说是该染色体的长臂q上的q21.1条带内。

  gene            <1..>5236
                     /gene="DUT"

gen指出了完整基因的具体位置。

mRNA            join(<282..561,1034..1172,2395..2486,3113..3157,
                     4447..4521,4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes mitochondrial form
                     of the protein"

mRNA指出了所有外显子的具体位置，mRNA中的基因串起来相当于完成了剪切的过程，就是成熟的mRNA。

值得注意的是，剪切后的mRNA有两种：

mRNA            join(<282..561,1034..1172,2395..2486,3113..3157,
                     4447..4521,4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes mitochondrial form
                     of the protein"

 mRNA            join(<1018..1172,2395..2486,3113..3157,4447..4521,
                     4673..4743,5180..>5236)
                     /gene="DUT"
                     /product="dUTPase"
                     /note="alternatively spliced; encodes nuclear form of the
                     protein"

上方的mRNA在前端多一个外显子，这一块将被翻译为定位线粒体的信号肽，从而翻译出线粒体型的蛋白质。下一种没有信号肽的，将形成细胞核型蛋白质。

     exon            <282..561
                     /gene="DUT"
                     /note="exon used only in the mitochondrial form of the
                     protein"
                     /number=1

exon段表示的含义是：当前序列所包含的282-561号碱基是“DUT”基因的第一个外显子。

2. ENA

欧洲核苷酸序列数据集（ENA）是由欧洲分子生物学研究室（EMBL）开发并维护。
ENA数据库: http://www.ebi.ac.uk/ena/

3. DDBJ

日本DNA数据库（DDBJ），由位于日本国立遗传学研究所（NIG）开发并负责维护。
DDBJ数据库： http://www.ddbj.nig.ac.jp/

4. INSDC

INSDC是三家数据库共同组成的国际核酸序列数据库合作联盟，通过每日的数据共通，使得他们几乎在任何时候都享有相同的数据。
INSDC数据库： http://www.insdc.org/

5. 基因组数据库Ensemble

这部分我们来学习使用Ensemble基因组数据库

Ensemble链接: http://www.ensembl.org

Ensembl是由欧洲生物信息学研究所（EBI）和英国桑格研究院合作开发。Esemble收入了各种动物的基因组，特别是那些离我们人类相近的动物。这些基因组的注释都是通过配套开发的软件自动添加的。

通过点击上方的连接我们可以前往Ensemble的官网。

从中我们可以获取到染色体一览图等信息：

点击DNA条段，可以进入查看条段的详细信息。从上文的学习可以得知，DUT基因位于15号染色体的长臂条段21.1附近，我们进入这条条段查看一下。

点击条段进入

进入之后可以看到该条段的详细信息

我们在基因图谱中找到DUT基因，并以其为中心放大

点击DUT基因在Ensemble数据库中的检索号，便可以得到该基因在数据库中的详细记录。

6. 微生物宏基因组数据库JCVI

JCVI链接: http://www.jcvi.org

美国基因组研究所（TIGR）致力于为生物基因组的研究，也有部分植物基因组项目。它是克莱格·凡特研究所（JCVI）的一部分，自1995年成立之初的两个基因组，至今已有超过700个基因组。TIGR是NCBI基因组资源的有力补充，因为它不仅拥有已经完成测序的基因组，还有哪些测序中的基因组信息。

三、二级核酸数据库

二级核酸数据库包含的内容很多，经常会用到的几个数据库有：NCBI下属的RefSeq数据库，dbEST数据库以及Gene数据库。

RefSeq数据库：参考序列数据库，是通过自动及人工精选出的非冗余数据库，包括基因组序列、转录序列和蛋白质序列。

dbEST数据库：表达序列标签数据库，包含来源于不同物种的表达序列标签（EST）

Gene数据库：为用户提供基因序列注释和检索服务，收录了来自5300多个物种的430万条基因记录

ncRNAdb：非编码RNA数据库，提供非编码RNA的序列和功能信息。包含来源于99种细菌，古细菌和真核生物的3万多条序列。

ncRNA链接: http://biobases.ibch.poznan.pl/ncRNA/

miRBase：主要存放已发表的microRNA序列和注释。可以分析microRNA在基因组中的定位和挖掘miRNA序列间的关系。

miRBase链接: http://www.mirbase.org/

至此，今天对基因库的学习基本上就告一段落。之后需要的就是多搜索几个DNA序列，来巩固自己的所学。