Volume 10 Supplement 3

Eighth International Conference on Bioinformatics (InCoB2009): Computational Biology

Open Access

MitoVariome: a variome database of human mitochondrial DNA

  • Yong Seok Lee1,
  • Woo-Yeon Kim1,
  • Mihyun Ji1,
  • Ji Han Kim1 and
  • Jong Bhak1Email author
Contributed equally
BMC Genomics200910(Suppl 3):S12

DOI: 10.1186/1471-2164-10-S3-S12

Published: 3 December 2009



Mitochondrial sequence variation provides critical information for studying human evolution and variation. Mitochondrial DNA provides information on the origin of humans, and plays a substantial role in forensics, degenerative diseases, cancers, and aging process. Typically, human mitochondrial DNA has various features such as HVSI, HVSII, single-nucleotide polymorphism (SNP), restriction enzyme sites, and short tandem repeat (STR).


We present a variome database (MitoVariome) of human mitochondrial DNA sequences. Queries against MitoVariome can be made using accession numbers or haplogroup/continent. Query results are presented not only in text but also in HTML tables to report extensive mitochondrial sequence variation information. The variation information includes repeat pattern, restriction enzyme site polymorphism, short tandem repeat, disease information as well as single nucleotide polymorphism. It also provides a graphical interface as Gbrowse displaying all variations at a glance. The web interface also provides the tool for assigning haplogroup based on the haplogroup-diagnostic system with complete human mitochondrial SNP position list and for retrieving sequences that users query against by using accession numbers.


MitoVariome is a freely accessible web application and database that enables human mitochondrial genome researchers to study genetic variation in mitochondrial genome with textual and graphical views accompanied by assignment function of haplogrouping if users submit their own data. Hence, the MitoVariome containing many kinds of variation features in the human mitochondrial genome will be useful for understanding mitochondrial variations of each individual, haplogroup, or geographical location to elucidate the history of human evolution.


Human mitochondrial genome is 16,569 base pairs (bp) in length, a closed, circular molecule located within the cytoplasmic mitochondria [1]. It has 37 genes: 22 transfer RNA genes, 13 protein coding genes, two ribosomal RNA genes, and one non-coding region containing the origin of replication. Twenty-eight of these genes are encoded by the heavy strand, and nine by the light strand.

Mitochondrial DNA is known for maternal clonal inheritance, rapid evolutionary rate, lack of introns, absence of recombination events, and haploidy [24]. Mitochondrial variations are linked to the origin of humans, and play a substantial role in forensics, degenerative diseases, cancers, and the aging process [5]. These mitochondrial DNA roles caused the development of various mitochondrial databases, i.e. MITOMAP containing information related to human evolution, disease [6], mtDB having the human mitochondrial DNAs which have not been deposited in a publicly available database such as GenBank [7], Mitome containing data on whole mitochondrial genomes for metazoan species [8], and MitoRes containing data on nuclear-encoded mitochondrial-related products for metazoan species. However, these databases haven't provided the integrated variation information. The environment of these web-based databases typically lacks a user-friendly and interactive platform for presenting all kinds of variation information.

Here we present an extensive web-based database, MitoVariome. The database will help researchers to identify mitochondrial DNA sequence variation, which is useful for understanding human evolution. It has integrated information such as gene and disease annotation. It provides a user-friendly interface to accept user-defined rules and display graphical views of variations and gene structures.

Methods and results

The MitoVariome database was described below (Figure 1).
Figure 1

Flowchart of MitoVariome database construction.

Data source

The database consists of over 5,000 complete human sequences of mitochondrial DNA. They were obtained from the NCBI GenBank database, MITOMAP at http://www.mitomap.org, and mtDB at http://www.genpat.uu.se/mtDB/. The nucleotide position and gene annotation corresponding to those of rCRS (revised Cambridge Reference Sequence) [9] were used. Geographic origin and disease information of human mitochondrial DNA were parsed from the above mentioned databases.

Assignment of haplogroup

Each human mitochondrial DNA sequence was aligned with the rCRS by ClustalW software and haplogroups were assigned defined by sequence differences from the rCRS. The haplogroup-diagnostic nucleotide sequence variants were assembled from the information from the references [10, 11]. The distribution of each haplogroup is shown in statistics page http://variome.kobic.re.kr/MitoVariome/statistics.jsp. The mitochondrial haplogroups in this study are divided into twenty groups broadly. We defined mitochondrial DNA sequences sharing over 95 percent of SNPs of each haplogroup as each haplogroup.

Finding of short tandem repeats, restriction enzyme sites

We used mreps [12] tool to detect short tandem repeats in mitochondrial DNA. We found all the mitochondrial tandem repeats which have at least two bases in a repeated unit and two copies as the minimum threshold. From the tandem repeat analysis, we found that each human mitochondrial DNA had on average 929 tandem repeat regions. Also, around 31.2% of each mitochondrial genome was detected as tandem repeat regions.

In order to predict restriction enzyme cleavage sites in human mitochondrial DNA, we used the 'restrict' program from the EMBOSS package. These analyses used 19 restriction enzymes which are frequently used in several studies on mitochondrial genomes. Each human mitochondrial DNA had on average 506 restriction enzyme sites http://variome.kobic.re.kr/MitoVariome/view_enzymelist.jsp. On average there were 99 MseI restriction enzyme sites (TTAA) at the most and one XhoI restriction enzyme site (CTCGAG) at the least.

Mapping disease information

Mapping disease information to mitochondrial DNA was performed by using mitochondrial disease information (Leber Hereditary Optic Neuropathy, Alzeimer's Disease, Parkinsons's Disease, and so on) on coding and non-coding regions. This information was obtained from Pubmed and MITOMAP. There are forty-two reported mitochondrial DNA base substitution diseases in our database (disease listing in http://variome.kobic.re.kr/MitoVariome/view_diseaselist.jsp, and 235 sites were related to base substitution diseases.

Visualization of results

Tables and two-dimensional view are provided to the users as visualize the final resultants. The tables provide detailed information of genome/gene and diseases associated with variation data. This information linked to GenBank and Pubmed databases with ID number and users can download sequences in FASTA files. The two-dimensional view result is shown by a method developed using Gbrowse [13]. This browser provides a graphical view of the mitochondrial genome/gene structure and the results obtained from above mentioned analyses.

Web interface

MitoVariome provides a web interface that allows researchers to browse the integrated information of mitochondrial DNA variation and disease. As shown in Figure 2, users can search using two entries; accession number (from NCBI), and haplogroup/continent. By submitting an accession number or haplogroup/continent information with haplogroup matching percent optionally (default 95%), the user can obtain genome or gene information related to the query and information for variations within the queried genome. The summary shows geographic information, haplogroup, sequence, reference and the number of SNP, disease related SNP, short tandem repeat, SNP within short tandem repeat, restriction enzyme site, and SNP within restriction enzyme site. MitoVariome also provides detailed disease information with associated variations in a table. For an advanced study, this database provides nucleotide sequences in a FASTA format for all the entries.
Figure 2

Search interface of MitoVariome. There are three functions; Database search, Haplogrouping and Get sequences. Users can search the MitoVariome database to get integrated information for variations using accession number and haplogroup/continent in Database search. By entering a SNPs list user provided, user data is matched with SNP positions of each haplogroup in Haplogrouping. By typing accession numbers, users can get the human mitochondrial genome sequences in Get sequences.

There are two additional applications in the database for further study, Haplogrouping and Get sequences. In order to assign haplogroup with user provided data, MitoVariome has the Haplogrouping application. By entering a user provided SNPs list, MitoVariome is matching with SNP positions in each haplogroup. It provides a table with over 50% matched haplogroup(s); haplogroup name, matched SNPs, unmatched SNPs, the count of matched SNPs and SNPs in each group, and percentage of matched over SNPs in each group. The second application, Get sequence, is provided in MitoVariome. Users can download the human mitochondrial genome sequences by typing accession numbers.

Here, we present an example of MitoVariome's usage for studying genetic variations in Korean mitochondrial genome. A Korean genome (SJK) derived mitochondrial genome was used as input for MitoVariome. SJK is the first full length Korean individual genome sequence [14]. We found 43 novel SNPs by comparing SNPs with rCRS. Two insertions and one deletion were detected with six non-synonymous SNPs. STR has 919 tandem repeat regions and 496 restriction enzyme sites showing similar numbers with the average of all human mitochondrial genome. 15 of 44 mtDNA variations are referred as the marker of "D4" haplogroup which is prevalent in Korea. Figure 3 provided variation information and a corresponding two-dimensional image.
Figure 3

Example of MitoVariome ('SJK'). Variation information and a corresponding two-dimensional image are shown. The top box is variation information. The top line with small bars is a scale bar for mitochondrial genome. The light brown horizontal bar is the Cambridge Reference Sequence. The yellow bars and blocks are CDS. The red vertical bars below the yellow bars are SNPs. Small black vertical bars are STRs. The last small black vertical bars are restriction enzyme sites.


MitoVariome was constructed as a database and analysis server for human mitochondrial DNA variation. It provides a platform to retrieve all kinds of variation information found in human mitochondrial DNA. MitoVariome provides assignment for the haplogroup with user provided SNP list by using rCRS information. MitoVariome database can be accessed at http://variome.kobic.re.kr/MitoVariome.


Other papers from the meeting have been published as part of BMC Bioinformatics Volume 10 Supplement 15, 2009: Eighth International Conference on Bioinformatics (InCoB2009): Bioinformatics, available online at http://www.biomedcentral.com/1471-2105/10?issue=S15.




The authors would like to thank Byun, Ha-na who designed the logo of MitoVariome. This research was supported by a grant from the Korean Research Institute of Bioscience and Biotechnology (KRIBB) Research Initiative Program and the Korean Ministry of Science and Technology (MOST) under grant number (M10757020001-07N5702-00110). We thank Maryana Bhak for editing the manuscript.

This article has been published as part of BMC Genomics Volume 10 Supplement 3, 2009: Eighth International Conference on Bioinformatics (InCoB2009): Computational Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2164/10?issue=S3.

Authors’ Affiliations

Korean Bioinformation Center (KOBIC), KRIBB


  1. Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, et al: Sequence and organization of the human mitochondrial genome. Nature. 1981, 290 (5806): 457-465. 10.1038/290457a0.View ArticlePubMedGoogle Scholar
  2. Brown WM, George M, Wilson AC: Rapid evolution of animal mitochondrial DNA. Proceedings of the National Academy of Sciences of the United States of America. 1979, 76 (4): 1967-1971. 10.1073/pnas.76.4.1967.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Giles RE, Blanc H, Cann HM, Wallace DC: Maternal inheritance of human mitochondrial DNA. Proceedings of the National Academy of Sciences of the United States of America. 1980, 77 (11): 6715-6719. 10.1073/pnas.77.11.6715.PubMed CentralView ArticlePubMedGoogle Scholar
  4. Olivo PD, Walle Van de MJ, Laipis PJ, Hauswirth WW: Nucleotide sequence evidence for rapid genotypic shifts in the bovine mitochondrial DNA D-loop. Nature. 1983, 306 (5941): 400-402. 10.1038/306400a0.View ArticlePubMedGoogle Scholar
  5. Taylor RW, Turnbull DM: Mitochondrial DNA mutations in human disease. Nature reviews. 2005, 6 (5): 389-402.PubMed CentralView ArticlePubMedGoogle Scholar
  6. Kogelnik AM, Lott MT, Brown MD, Navathe SB, Wallace DC: MITOMAP: a human mitochondrial genome database. Nucleic acids research. 1996, 24 (1): 177-179. 10.1093/nar/24.1.177.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Ingman M, Gyllensten U: mtDB: Human Mitochondrial Genome Database, a resource for population genetics and medical sciences. Nucleic acids research. 2006, D749-751. 10.1093/nar/gkj010. 34 Database
  8. Lee YS, Oh J, Kim YU, Kim N, Yang S, Hwang UW: Mitome: dynamic and interactive database for comparative mitochondrial genomics in metazoan animals. Nucleic acids research. 2008, D938-942. 36 Database
  9. Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N: Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nature genetics. 1999, 23 (2): 147-10.1038/13779.View ArticlePubMedGoogle Scholar
  10. Kong QP, Bandelt HJ, Sun C, Yao YG, Salas A, Achilli A, Wang CY, Zhong L, Zhu CL, Wu SF, et al: Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations. Human molecular genetics. 2006, 15 (13): 2076-2086. 10.1093/hmg/ddl130.View ArticlePubMedGoogle Scholar
  11. Brandon MC, Lott MT, Nguyen KC, Spolim S, Navathe SB, Baldi P, Wallace DC: MITOMAP: a human mitochondrial genome database--2004 update. Nucleic acids research. 2005, D611-613. 33 Database
  12. Kolpakov R, Bana G, Kucherov G: mreps: Efficient and flexible detection of tandem repeats in DNA. Nucleic acids research. 2003, 31 (13): 3672-3678. 10.1093/nar/gkg617.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Stein LD, Mungall C, Shu S, Caudy M, Mangone M, Day A, Nickerson E, Stajich JE, Harris TW, Arva A, et al: The generic genome browser: a building block for a model organism system database. Genome research. 2002, 12 (10): 1599-1610. 10.1101/gr.403602.PubMed CentralView ArticlePubMedGoogle Scholar
  14. Ahn SM, Kim TH, Lee S, Kim D, Ghang H, Kim DS, Kim BC, Kim SY, Kim WY, Kim C, et al: The first Korean genome sequence and analysis: Full genome sequencing for a socio-ethnic group. Genome research. 2009Google Scholar


© Lee et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.