Comparative study of human mitochondrial proteome reveals extensive protein subcellular relocalization after gene duplications
© Wang et al; licensee BioMed Central Ltd. 2009
Received: 19 May 2009
Accepted: 30 November 2009
Published: 30 November 2009
Gene and genome duplication is the principle creative force in evolution. Recently, protein subcellular relocalization, or neolocalization was proposed as one of the mechanisms responsible for the retention of duplicated genes. This hypothesis received support from the analysis of yeast genomes, but has not been tested thoroughly on animal genomes. In order to evaluate the importance of subcellular relocalizations for retention of duplicated genes in animal genomes, we systematically analyzed nuclear encoded mitochondrial proteins in the human genome by reconstructing phylogenies of mitochondrial multigene families.
The 456 human mitochondrial proteins selected for this study were clustered into 305 gene families including 92 multigene families. Among the multigene families, 59 (64%) consisted of both mitochondrial and cytosolic (non-mitochondrial) proteins (mt-cy families) while the remaining 33 (36%) were composed of mitochondrial proteins (mt-mt families). Phylogenetic analyses of mt-cy families revealed three different scenarios of their neolocalization following gene duplication: 1) relocalization from mitochondria to cytosol, 2) from cytosol to mitochondria and 3) multiple subcellular relocalizations. The neolocalizations were most commonly enabled by the gain or loss of N-terminal mitochondrial targeting signals. The majority of detected subcellular relocalization events occurred early in animal evolution, preceding the evolution of tetrapods. Mt-mt protein families showed a somewhat different pattern, where gene duplication occurred more evenly in time. However, for both types of protein families, most duplication events appear to roughly coincide with two rounds of genome duplications early in vertebrate evolution. Finally, we evaluated the effects of inaccurate and incomplete annotation of mitochondrial proteins and found that our conclusion of the importance of subcellular relocalization after gene duplication on the genomic scale was robust to potential gene misannotation.
Our results suggest that protein subcellular relocalization is an important mechanism for the retention and gain of function of duplicated genes in animal genome evolution.
Gene duplication is an important evolutionary process that plays a key role in generating new genomic information in all the three domains of life: Eubacteria, Archaea and Eukarya [1–5]. Various processes can cause gene duplication on the molecular level, including unequal crossovers, retroposition, or whole chromosome/genome duplication [6–8]. New functional genes resulting from gene duplication are retained in the genome through the processes of subfunctionalization and neofunctionalization [9, 10]. Subfunctionalization refers to a situation when each of the daughter genes adopts only a partial function of the parental gene, while neofunctionalization refers to the gain of new functions by the duplicate, usually related to their ancestor's function .
From the evolutionary perspective, neofunctionalization presents more interest than subfunctionalization because it results in the increase of the total genetic information . Several mechanisms have been invoked to explain the gain of novel gene function for duplicated genes such as dosage compensation, epigenetic complementation, moonlighting, and catalytic promiscuity [12–17]. Interestingly, functional divergence can in some cases precede (and facilitate) gene duplication through allelic divergence . Recently, protein subcellular relocalization or neolocalization has been proposed as a key event for generating new functional genes after duplication [19, 20]. Such neolocalization can be achieved by the gain or loss of N-terminal targeting peptide sequences that can direct the products of duplicated genes from the cytosol to mitochondria, endoplasmic reticulum and chloroplast or vice versa.
The idea of subcellular relocalization underlying the gain of function for duplicated genes has been tested in two yeast genomes by Marques et al. (2008). The authors demonstrated that about one-third of the duplicated genes retained in the yeast genomes had undergone protein subcellular relocalization following whole genome duplication . A few anecdotal observations suggest that neolocalization after duplication also occurs in animal genomes [21, 22], however the magnitude of this process has not been explored. In this study, we performed a systematic survey of subcellular relocalization following gene duplication in the human genome by analyzing nuclear-encoded mitochondrial protein families.
Mitochondria, cell organelles present in nearly all eukaryotes, are instrumental for the production of ATP through oxidative phosphorylation process, and are also involved in heme biosynthesis, cell metabolism, apoptosis, and Fe/S cluster biosynthesis. The complex functions of mitochondria demand a proteome composed of over a thousand of proteins, more than 98% of which are nuclear encoded, which suggests these organelles should play a major role in the process of neolocalization . Hence, the exploration of nuclear encoded mitochondrial gene families is an ideal system to test subcellular relocalization of duplicated genes in the evolution of animal genomes.
Results and Discussion
Subcellular relocalization as a mechanism underlying protein functional divergence
Summary of human mitochondrial multigene families
Total number of proteins
Total number of family
At the root of vertebrates
After the emergence of fish
Mitochondrial- cytosolic families (mt-cy families)
ACAA2, ALDH1B1, CPT2, GPAM, HMGCS2, NDUFA4, SH3BP5, HIBADH, ABCB10, CABC1, ARG2, ATP5A1, BNIP3, CYP11B2, CA5A, CPS1, CDS2, CYB5B, DECR1, DGUOK, DLD, DNAJA3, GFM1, TUFM, FTMT, GK, SHMT2, GPX4, HSPA9, ALAS2, HTRA2, HK1, IDH2, MTIF2, PPA2, AK3, CKMT2, ME2, MFN2, MGST1, MIPEP, NME4, NFS1, OAS2, OPA1, SLC25A15, PCK2, PHB, PPIF, PRDX3, SIRT3, IARS2, SARS2, DNAJC19, TMLHE, TOP1MT, TRAP1, TST, OXR1
Mitochondrial- mitochondrial families (mt-mt families)
AIFM1, MAOA, ATP5G2, BCL2L1, D2HGDH, GLUD2, GRPEL1, LETM1, MCART1, PMPCB, MTCH2, MTERF, ENDOG, OAT, BCKDHA, PDHB, PDK4, PDP2, MTRF1, RHOT1, SCO2, TIMM17A, TOMM40, VDAC3, COX4I2, COX7A1, COX6A1, COX6B1, MCCC2, OXCT2, MRPS18A, MRPS10, MRPS30
For each human protein in the 92 gene families, we retrieved orthologs in mouse (Mus musculus), chicken (Gallus gallus), fish (Danio rerio), fruit fly (Drosophila melanogaster), mosquito (Anopheles gambiae) and nematode (Caenorhabditis elegans) from the Homologene database at NCBI. Phylogenetic analyses were conducted for each family and the time of occurrence of gene duplications in relationship to major divergences in animal evolution was evaluated. Because this study is based on the human mitochondrial proteome, only branches of the phylogenetic tree leading to humans were investigated. For mitochondrial-cytosolic (mt-cy) families that underwent several rounds of duplications, we only considered the duplications that were followed by subcellular relocalizations.
Among the 59 mt-cy families, twenty-seven (45.8%) were inferred to undergo gene duplication prior to the protostome/deuterostome divergence, twenty-nine (49.2%) after the protostome/deuterostome divergence but before that of fish/tetrapods and only three families (5%) within the tetrapod lineage (Table 1). Among the 33 mt-mt families, eight (24%) underwent gene duplication prior to the protostome/deuterostome split, nineteen (57.6%) prior to the fish/tetrapod split and six (18.4%) within the tetrapod lineage. The observation that the majority of investigated families experienced gene duplication between the protostome/deuterostome and fish/tetrapod divergences correlates well with the two rounds of genome duplication at the early stage of vertebrate evolution [24, 25]. However, the scarcity of more recent subcellular relocalization events is surprising, especially considering very high rates of gene birth and death in animal genomes .
Bidirectional relocalization of proteins encoded by duplicated genes in mitochondrial-cytosolic (mt-cy) two-gene families
In order to get insight into the direction of protein subcellular relocalization, we explored mt-cy gene families in which human genes are represented by two copies, one functioning in the cytosol and the other in mitochondria. Among the 24 such families, one third appeared to have its original function in mitochondria with the products of duplicated genes being relocalized to other cellular compartments, another third showed the opposite direction of protein relocalization and for the rest the direction of relocalization could not be determined due to the lack of outgroup information.
Both TOP1 and TOP1MT consist of four domains: a N-terminal localization domain, a core domain, a linker domain and a C-terminal domain . Sequence comparisons show that the N-terminal domain of TOP1MT consists of a mitochondrial targeting signal, while the N-terminus of TOP1 contains a nuclear localization signal. This implies that the change in the N-terminal targeting sequence of TOP1MT helped the protein direct itself to mitochondria and eventually to acquire a new mitochondrial function. It should be noted that the lack of TOP1MT in invertebrates does not mean that type I topoisomerases are not needed for mitochondrial replication and transcription in this group. Recent studies have shown that DNA topoisomerase IIIα from the type IA subfamily has mitochondrial localization in Drosophila melanogaster [29, 30].
Multiple subcellular relocalizations after gene duplications
Expansion of mitochondrial proteome by gene duplications
The presence of 33 mt-mt families among the 92 multigene families supports the notion that gene duplication also contributes to mitochondrial proteome expansion , although the average family size of mt-mt families is smaller than that of mt-cy families (2.4 vs 5.7). Our phylogenetic analyses of these mt-mt families showed that such duplications occurred at different stages in animal evolution, predating the divergence of the protostome/deuterostome lineages, within the vertebrate lineage, and within the mammalian lineage. In general, these families consist of proteins with similar functions that have been retained by subfunctionalization as different subunits or isoforms. Furthermore, the expression of these genes often shows tissue specificity such that one copy in the gene family is expressed ubiquitously, while the other(s) is/are expressed in specific tissues. For example, human SCO2 homolog (SCO cytochrome oxidase deficient homolog 2 (yeast)) is expressed ubiquitously while the SCO1 homolog is predominantly expressed in muscle, heart, and brain, the tissues featured by high rates of oxidative phosphorylation .
Evolutionary modifications of relocalized proteins at sequence level
We further investigated the differences in protein functional domains between mitochondrial and non-mitochondrial proteins based on a suggestion that protein function and/or protein functional efficiency can be modified upon the change in its subcellular location . Here we compared the distribution of mitochondrial Pfam domains that were previously found only in eukaryotic (excluding human) mitochondrial proteins  among members of mt-cy protein families. We found a significant difference in this distribution (figure 5B): 53% of mitochondrial proteins have mitochondrial domains, but only 16% of non-mitochondrial proteins have such domains (chi-square p-value is 3.5e-9). This result indicates that subcellular relocalizations were characterized by the formation of mitochondrial protein domains or their loss in nuclear copies during evolution.
Effects of inaccurate or incomplete cellular component annotation
Limited knowledge of protein subcellular location also prevents us from discovering cases of a different form of subcellular relocalization, called sublocalization, in which duplicated genes become targeted to a subset of their ancestral cellular compartments . For example, the glutamate dehydrogenase family was grouped into mt-mt families since both GLUD1 and GLUD2 are located to mitochondria based on the Genbank and Swissprot annotations. However Rosso et al. recently reported that GLUD1 located to both cytoplasm and mitochondria while GLUD2 became specifically localized to mitochondria owing to a single positively selected amino acid substitution at the N-terminal targeting sequence . In addition, if the ancestral protein had dual localizations but only one of them was annotated, then neolocalization would be inferred instead of sublocalization. The latter problem should be especially pronounced if non-model species are used as outgroups since the localizations of proteins in these species are not thoroughly studied. To investigate the potential effect of this bias, we searched protein localizations in Drosophila melanogaster for the families duplicated in the vertebrate lineage. None of these Drosophila proteins are annotated to have dual localizations. Thus available data suggests that neolocalization rather than sublocalization is the prevalent mode of evolution of duplicated genes studied here.
Protein subcellular relocalization was proposed as an evolutionary mechanism for generating new functional genes after gene duplication . This mechanism was studied in yeast genomes but only received support from individual cases/families in animal genomes [20–22]. Here we systematically investigated human mitochondrial protein families and found that around two thirds of multigene families have protein members that underwent subcellular relocalization after gene duplication. These subcellular relocalizations can occur between mitochondria and another subcellular compartment as well as among several compartments. Comparative sequence analyses showed that the subcellular relocalization processes were primarily enabled via the gain or loss of N-terminal mitochondrial targeting sequences.
After evaluating possible effects of incomplete or incorrect annotations, we conclude that our observation of the subcellular relocalization after gene duplication on the genomic scale was robust to misannotations. Surprisingly our results indicate a scarcity of recent subcellular relocalization events and suggest that protein subcellular relocalization was more important in obtaining new functional genes at the early stages of animal genome evolution. The observation that subcellular relocalization rarely follows recent gene duplications needs to be further investigated because our dataset might be biased in terms of annotation availability. This investigation, together with the previously published finding in yeast , indicates that subcellular relocalization is an important mechanism in the retention and gain of function of duplicated genes over the course of eukaryotic genome evolution.
Retrieving human mitochondrial proteins and identifying the paralogs
We used the MitoP2 database http://www.mitop.de:8080/mitop2/ to retrieve human mitochondrial proteins . MitoP2 is a manually annotated database of mitochondrial proteins that integrates computational predictions, proteome mapping, mutant screening, expression profiling, protein-protein interactions and cellular localization. Among 920 inferred human mitochondrial proteins in MitoP2 at the time of our study, we selected 456 that were also annotated as mitochondrial proteins in the Swissprot database. These 456 protein sequences were used to query the Swissprot database using blastp program. The hit sequences with e-value smaller than 0.001 and bit scores equal to or greater than 100 were kept for further analyses. Single linkage clustering was then carried out to cluster the sequences into 305 families .
Ortholog collection and phylogeny construction
Gene families were removed from analysis if all human proteins in the family were products of alternative splicing of the same gene or the records of the proteins no longer existed in Genbank. For each human protein in the remaining 92 multigene families, we retrieved orthologs in mouse (Mus musculus), chicken (Gallus gallus), fish (Danio rerio), fruit fly (Drosophila melanogaster), mosquito (Anopheles gambiae) and nematode (Caenorhabditis elegans) from the Homologene database at NCBI. The subcellular localizations of proteins in each family were inferred from Swissprot and Genbank databases. The protein sequences in each family were aligned with T-coffee using the default settings . The aligned sequences were manually inspected and then imported into the Mega package to construct neighbor joining and parsimony trees with 1000 bootstrap replicates . For each family, we further performed a maximum likehood search with 100 bootstrap replicates as implemented in the PHYML (v.2.4.4) program using the WAG model with estimated amino acid frequencies and 4 gamma categories. The inferred topologies were congruent among different tree-making methods except for some minor differences (not shown). We selected maximum likelihood trees for illustrations.
Comparative analyses of protein sequences
For human proteins in each family, we retrieved information of N-terminal mitochondrial targeting signal (TARGETP), Pfam domains (MITODOMAIN) and maestro scores developed by Calvo and coauthors for predicting mitochondrial proteins . The distributions of these data were plotted in R package http://www.r-project.org/ for mitochondrial and non-mitochondrial protein members in the analyzed multigene families. The Pearson chi-square tests were applied to test whether the distributional proportions of these genomic criteria were the same for mitochondrial and non-mitochondrial proteins.
We thank Karri Haen, Ehsan Kayal, and two anonymous reviewers for valuable comments on earlier versions of the manuscript, and Pioneer Computational Biology Fellowship for funding.
- Ohno S: Evolution by Gene Duplication. 1970, Springer-VerlagView ArticleGoogle Scholar
- Otto SP, Yong P: The evolution of gene duplicates. Adv Genet. 2002, 46: 451-483. full_text.View ArticlePubMedGoogle Scholar
- Bridges CB: The bar "gene" a duplication. Science. 1936, 83: 210-211. 10.1126/science.83.2148.210.View ArticlePubMedGoogle Scholar
- Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18: 292-298. 10.1016/S0169-5347(03)00033-8.View ArticleGoogle Scholar
- Conant GC, Wolfe KH: Turning a hobby into a job: how duplicated genes find new functions. Nat Rev Genet. 2008, 9: 938-950. 10.1038/nrg2482.View ArticlePubMedGoogle Scholar
- Babushok DV, Ostertag EM, Kazazian HHJ: Current topics in genome evolution: molecular mechanisms of new gene formation. Cell Mol Life Sci. 2007, 64: 542-554. 10.1007/s00018-006-6453-4.View ArticlePubMedGoogle Scholar
- Semon M, Wolfe KH: Consequences of genome duplication. Curr Opin Genet Dev. 2007, 17: 505-512. 10.1016/j.gde.2007.09.007.View ArticlePubMedGoogle Scholar
- Kaessmann H, Vinckenbosch N, Long M: RNA-based gene duplication: mechanistic and evolutionary insights. Nat Rev Genet. 2009, 10: 19-31. 10.1038/nrg2487.PubMed CentralView ArticlePubMedGoogle Scholar
- Katju V, Lynch M: On the formation of novel genes by duplication in the Caenorhabditis elegans genome. Mol Biol Evol. 2006, 23: 1056-1067. 10.1093/molbev/msj114.View ArticlePubMedGoogle Scholar
- Freeling M: The evolutionary position of subfunctionalization, downgraded. Genome Dyn. 2008, 4: 25-40. full_text.View ArticlePubMedGoogle Scholar
- Francino MP: An adaptive radiation model for the origin of new gene functions. Nat Genet. 2005, 37: 573-577. 10.1038/ng1579.View ArticlePubMedGoogle Scholar
- Rodin SN, Parkhomchuk DV, Rodin AS, Holmquist GP, Riggs AD: Repositioning-dependent fate of duplicate genes. DNA Cell Biol. 2005, 24: 529-542. 10.1089/dna.2005.24.529.View ArticlePubMedGoogle Scholar
- Papp B, Pal C, Hurst LD: Dosage sensitivity and the evolution of gene families in yeast. Nature. 2003, 424: 194-197. 10.1038/nature01771.View ArticlePubMedGoogle Scholar
- Qian W, Zhang J: Gene dosage and gene duplicability. Genetics. 2008, 179: 2319-2324. 10.1534/genetics.108.090936.PubMed CentralView ArticlePubMedGoogle Scholar
- McLoughlin SY, Copley SD: A compromise required by gene sharing enables survival: Implications for evolution of new enzyme activities. Proc Natl Acad Sci USA. 2008, 105: 13497-13502. 10.1073/pnas.0804804105.PubMed CentralView ArticlePubMedGoogle Scholar
- Copley SD: Enzymes with extra talents: moonlighting functions and catalytic promiscuity. Curr Opin Chem Biol. 2003, 7: 265-272. 10.1016/S1367-5931(03)00032-2.View ArticlePubMedGoogle Scholar
- Johnston CR, O'Dushlaine C, Fitzpatrick DA, Edwards RJ, Shields DC: Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication. Mol Biol Evol. 2007, 24: 315-323. 10.1093/molbev/msl162.View ArticlePubMedGoogle Scholar
- Proulx SR, Phillips PC: Allelic divergence precedes and promotes gene duplication. Evolution. 2006, 60: 881-892.View ArticlePubMedGoogle Scholar
- Byun-McKay SA, Geeta R: Protein subcellular relocalization: a new perspective on the origin of novel genes. Trends Ecol Evol. 2007, 22: 338-344. 10.1016/j.tree.2007.05.002.View ArticlePubMedGoogle Scholar
- Marques AC, Vinckenbosch N, Brawand D, Kaessmann H: Functional diversification of duplicate genes through subcellular adaptation of encoded proteins. Genome Biol. 2008, 9: R54-10.1186/gb-2008-9-3-r54.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosso L, Marques AC, Reichert AS, Kaessmann H: Mitochondrial targeting adaptation of the hominoid-specific glutamate dehydrogenase driven by positive Darwinian selection. PLoS Genet. 2008, 4: e1000150-10.1371/journal.pgen.1000150.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosso L, Marques AC, Weier M, Lambert N, Lambot MA, Vanderhaeghen P, Kaessmann H: Birth and rapid subcellular adaptation of a hominoid-specific CDC14 protein. PLoS Biol. 2008, 6: e140-10.1371/journal.pbio.0060140.PubMed CentralView ArticlePubMedGoogle Scholar
- Ryan MT, Hoogenraad NJ: Mitochondrial-nuclear communications. Annu Rev Biochem. 2007, 76: 701-722. 10.1146/annurev.biochem.76.052305.091720.View ArticlePubMedGoogle Scholar
- Gu X, Wang Y, Gu J: Age distribution of human gene families shows significant roles of both large- and small-scale duplications in vertebrate evolution. Nat Genet. 2002, 31: 205-209. 10.1038/ng902.View ArticlePubMedGoogle Scholar
- Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005, 3: e314-10.1371/journal.pbio.0030314.PubMed CentralView ArticlePubMedGoogle Scholar
- Demuth JP, Hahn MW: The life and death of gene families. Bioessays. 2009, 31: 29-39. 10.1002/bies.080085.View ArticlePubMedGoogle Scholar
- Wang JC: Cellular roles of DNA topoisomerases: a molecular perspective. Nat Rev Mol Cell Biol. 2002, 3: 430-440. 10.1038/nrm831.View ArticlePubMedGoogle Scholar
- Zhang H, Pommier Y: Mitochondrial topoisomerase I sites in the regulatory D-loop region of mitochondrial DNA. Biochemistry. 2008, 47: 11196-11203. 10.1021/bi800774b.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Lyu YL, Wang JC: Dual localization of human DNA topoisomerase IIIalpha to mitochondria and nucleus. Proc Natl Acad Sci USA. 2002, 99: 12114-12119. 10.1073/pnas.192449499.PubMed CentralView ArticlePubMedGoogle Scholar
- Plank JL, Chu SH, Pohlhaus JR, Wilson-Sali T, Hsieh TS: Drosophila melanogaster topoisomerase IIIalpha preferentially relaxes a positively or negatively supercoiled bubble substrate and is essential during development. J Biol Chem. 2005, 280: 3564-3573. 10.1074/jbc.M411337200.View ArticlePubMedGoogle Scholar
- Frye RA: Phylogenetic classification of prokaryotic and eukaryotic Sir2-like proteins. Biochem Biophys Res Commun. 2000, 273: 793-798. 10.1006/bbrc.2000.3000.View ArticlePubMedGoogle Scholar
- Michishita E, Park JY, Burneskis JM, Barrett JC, Horikawa I: Evolutionarily conserved and nonconserved cellular localizations and functions of human SIRT proteins. Mol Biol Cell. 2005, 16: 4623-4635. 10.1091/mbc.E05-01-0033.PubMed CentralView ArticlePubMedGoogle Scholar
- North BJ, Verdin E: Sirtuins: Sir2-related NAD-dependent protein deacetylases. Genome Biol. 2004, 5: 224-10.1186/gb-2004-5-5-224.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang YH, Chen YH, Zhang CY, Nimmakayalu MA, Ward DC, Weissman S: Cloning and characterization of two mouse genes with homology to the yeast Sir2 gene. Genomics. 2000, 69: 355-369. 10.1006/geno.2000.6360.View ArticlePubMedGoogle Scholar
- Smits P, Smeitink JA, Heuvel van den LP, Huynen MA, Ettema TJ: Reconstructing the evolution of the mitochondrial ribosomal proteome. Nucleic Acids Res. 2007, 35: 4686-4703. 10.1093/nar/gkm441.PubMed CentralView ArticlePubMedGoogle Scholar
- Petruzzella V, Tiranti V, Fernandez P, Ianna P, Carrozzo R, Zeviani M: Identification and characterization of human cDNAs specific to BCS1, PET112, SCO1, COX15, and COX11, five genes involved in the formation and function of the mitochondrial respiratory chain. Genomics. 1998, 54: 494-504. 10.1006/geno.1998.5580.View ArticlePubMedGoogle Scholar
- Omura T: Mitochondria-targeting sequence, a multi-role sorting sequence recognized at all steps of protein import into mitochondria. J Biochem. 1998, 123: 1010-1016.View ArticlePubMedGoogle Scholar
- Emanuelsson O, Nielsen H, Brunak S, von Heijne G: Predicting subcellular localization of proteins based on their N-terminal amino acid sequence. J Mol Biol. 2000, 300: 1005-1016. 10.1006/jmbi.2000.3903.View ArticlePubMedGoogle Scholar
- Dominguez D, Montserrat-Sentis B, Virgos-Soler A, Guaita S, Grueso J, Porta M, Puig I, Baulida J, Franci C, Garcia de Herreros A: Phosphorylation regulates the subcellular location and activity of the snail transcriptional repressor. Mol Cell Biol. 2003, 23: 5078-5089. 10.1128/MCB.23.14.5078-5089.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Calvo S, Jain M, Xie X, Sheth SA, Chang B, Goldberger OA, Spinazzola A, Zeviani M, Carr SA, Mootha VK: Systematic identification of human mitochondrial disease genes through integrative genomics. Nat Genet. 2006, 38: 576-582. 10.1038/ng1776.View ArticlePubMedGoogle Scholar
- Prokisch H, Andreoli C, Ahting U, Heiss K, Ruepp A, Scharfe C, Meitinger T: MitoP2: the mitochondrial proteome database--now including mouse data. Nucleic Acids Res. 2006, 34: D705-11. 10.1093/nar/gkj127.PubMed CentralView ArticlePubMedGoogle Scholar
- Gu Z, Cavalcanti A, Chen FC, Bouman P, Li WH: Extent of gene duplication in the genomes of Drosophila, nematode, and yeast. Mol Biol Evol. 2002, 19: 256-262.View ArticlePubMedGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302: 205-217. 10.1006/jmbi.2000.4042.View ArticlePubMedGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9: 299-306. 10.1093/bib/bbn017.PubMed CentralView ArticlePubMedGoogle Scholar
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.