- Open Access
The rhomboids: a nearly ubiquitous family of intramembrane serine proteases that probably evolved by multiple ancient horizontal gene transfers
© Koonin et al.; licensee BioMed Central Ltd. 2003
Received: 30 September 2002
Accepted: 3 February 2003
Published: 28 February 2003
The rhomboid family of polytopic membrane proteins shows a level of evolutionary conservation unique among membrane proteins. They are present in nearly all the sequenced genomes of archaea, bacteria and eukaryotes, with the exception of several species with small genomes. On the basis of experimental studies with the developmental regulator rhomboid from Drosophila and the AarA protein from the bacterium Providencia stuartii, the rhomboids are thought to be intramembrane serine proteases whose signaling function is conserved in eukaryotes and prokaryotes.
Phylogenetic tree analysis carried out using several independent methods for tree constructions and the corresponding statistical tests suggests that, despite its broad distribution in all three superkingdoms, the rhomboid family was not present in the last universal common ancestor of extant life forms. Instead, we propose that rhomboids evolved in bacteria and have been acquired by archaea and eukaryotes through several independent horizontal gene transfers. In eukaryotes, two distinct, ancient acquisitions apparently gave rise to the two major subfamilies, typified by rhomboid and PARL (presenilins-associated rhomboid-like protein), respectively. Subsequent evolution of the rhomboid family in eukaryotes proceeded by multiple duplications and functional diversification through the addition of extra transmembrane helices and other domains in different orientations relative to the conserved core that harbors the protease activity.
Although the near-universal presence of the rhomboid family in bacteria, archaea and eukaryotes appears to suggest that this protein is part of the heritage of the last universal common ancestor, phylogenetic tree analysis indicates a likely bacterial origin with subsequent dissemination by horizontal gene transfer. This emphasizes the importance of explicit phylogenetic analysis for the reconstruction of ancestral life forms. A hypothetical scenario for the origin of intracellular membrane proteases from membrane transporters is proposed.
Polytopic transmembrane proteins are, in general, not particularly strongly conserved during evolution. Inspection of the database of Clusters of Orthologous Groups of proteins (COGs)  revealed only one family of such proteins that is represented in most of the sequenced bacterial, archaeal and eukaryotic genomes. The prototype of this family is the rhomboid (RHO) protein from Drosophila melanogaster, a developmental regulator involved in epidermal growth factor (EGF)-dependent signaling pathways [2–4]. Not only were homologs of rhomboid detected in prokaryotes and eukaryotes, but the pattern of sequence conservation in this family appeared uncharacteristic of nonenzymatic membrane proteins, such as transporters [5, 6]. Specifically, several polar amino-acid residues are conserved in nearly all members of the rhomboid family, suggesting the possibility of an enzymatic activity. As three of these conserved residues were histidines, it has been hypothesized that rhomboid-family proteins could function as metal-dependent membrane proteases [5, 6]. Recently, however, it has been shown that RHO cleaves a transmembrane helix (TMH) in the membrane-bound precursor of the TGFα-like growth factor Spitz, enabling the released Spitz to activate the EGF receptor, and that a conserved serine and a conserved histidine in RHO are essential for this cleavage [7, 8]. Thus, it appears that rhomboid-family proteins are a distinct group of intramembrane serine proteases. Altogether, the genome of Drosophila encodes seven RHO paralogs (now designated RHO1-7, with the original rhomboid becoming RHO-1), at least three of which are involved in distinct EGF-dependent pathways, apparently through proteolytic activation of diverse ligands of the EGF receptor [9, 10].
The newly discovered intramembrane proteolytic activity of RHO places the rhomboid family within the framework of regulated intramembrane proteolysis (RIP), a new paradigm of signal transduction, which appears to be prominent in all forms of life [11, 12]. Under RIP, signaling proteins undergo site-specific proteolysis within TMH, resulting in the release of active fragments, which are the actual effectors in signal tranduction cascades. Until recently, the only characterized cases of RIP in eukaryotes involved presenilin-1, an aspartyl protease, which cleaves a transmembrane helix in type-1 membrane proteins such as amyloid β-precursor protein (AβPP), Notch and Ire1 , and the metalloprotease S2P, which cleaves a TMH in a type-2 transmembrane protein, the sterol-dependent transcription factor SREBP . Notably, S2P has highly conserved bacterial homologs, and the protease domain of presenilins also might be homologous to bacterial and archaeal type IV prepilin peptidases, although, in this case, the sequence similarity is low [14, 15].
In the case of the rhomboid family, the existence of homologs of RHO in most prokaryotes is particularly remarkable because animal RHO proteins are involved in signaling pathways that are not found outside metazoa, which seems to make functional conservation in prokaryotes a remote possibility. The only prokaryotic protein of the rhomboid family that has been characterized experimentally in considerable detail is AarA from the bacterium Providencia stuartii [16, 17]. This protein is involved in the export of a quorum-sensing peptide, a function that, in physiological terms, resembles that of RHO, although the signaling molecules, other than RHO and AarA, are obviously unrelated . In a striking recent development, two independent research groups have shown that several bacterial rhomboid-family proteins, including AarA, can cleave the EGF receptor ligands (Spitz, Keren and Gurken) that are normally cleaved by RHO paralogs [19, 20]. The cleavage depended on the conserved serine and histidine residues  and, moreover, transgenic flies that expressed AarA developed a phenotype indistinguishable from that induced by overexpression of RHO, whereas RHO could substitute for AarA in Providencia stuartii . These unexpected findings demonstrated the conservation of a RIP mechanism producing extracellular signals in eukaryotes and prokaryotes. Eukaryotic rhomboid family proteins seem to show considerable functional variability; in particular, cross-talk might exist between different RIP pathways. A distinct representative of the rhomboid family has been shown to physically interact with presinilins 1 and 2, and was accordingly named presenilins-associated rhomboid-like protein (PARL) . The yeast ortholog of PARL has been suggested to participate in the processing of cytochrome c peroxidase precursor during its import into the mitochondrion .
The near ubiquity of the rhomboid family among bacteria, archaea and eukaryotes, along with the remarkable functional conservation, suggests that a signaling mechanism mediated by rhomboids might have functioned already in the last common ancestor of all extant life forms, with subsequent loss in several lineages. To address this possibility, we performed a detailed phylogenetic analysis of the rhomboid family.
Results and discussion
Sequence and structural features and phyletic distribution of the rhomboid family
When examining the multiple alignment of the rhomboid superfamily proteins, we noticed that several eukaryotic members appear to be inactivated proteases, as indicated by the loss of the predicted catalytic serine or histidine (Figure 1, and data not shown); these inactivated forms could be regulators of active rhomboid proteases. Several other proteins lack one or more of the conserved residues in TMH3; it remains unclear whether or not these are active proteases.
Bacterial and archaeal members of the rhomboid superfamily contain six TMH, whereas the eukaryotic members typically have an additional seventh TMH, which may be attached to the core either from the amino terminus or from the carboxyl terminus as discussed below.
The phyletic distribution pattern of the rhomboid family shows that this intramembrane protease is extremely common in all three kingdoms of life, but is not necessarily essential for cell function. Rhomboids are missing in the microsporidian Encephalitozoon cuniculi, a eukaryotic intracellular parasite with a highly degraded genome, the archaea Methanothermobacter thermoautotrophicus and Thermoplasma volcanium, and several bacterial species, primarily parasites with small genomes but also species with moderately sized genomes, such as Xylella fastidiosum (see COG0705 at ). In two instances, a representative of the rhomboid family is present in only one of a pair of relatively close genomes (present in T. acidophilum but missing in T. volcanium; present in the spirochete Treponema pallidum but missing in the related bacterium Borrelia burgdorferi), which suggests relatively recent, repeated losses of this gene. Most of the prokaryotic species have a single gene coding for a rhomboid-family protein, although some have two or three paralogs (see COG0705 ); in contrast, eukaryotes show expansion of the rhomboid family, with seven members in Drosophila, and as many as 13 in Arabidopsis.
Phylogeny and evolutionary history of the rhomboid family
The second eukaryotic subfamily, which we designated the PARL subfamily, after PARL, the human ortholog of Drosophila RHO7 , resides within a large, heterogeneous prokaryotic cluster (Figure 2). Within this subfamily, PARL and its orthologs from other animals and from fungi have distinct domain architecture, with an extra TMH added to the amino terminus of the core, whereas the rest have only the core (a carboxy-terminal TMH and a ubiquitin-associated domain are appended in one Arabidopsis protein; Figure 2). Thus, the existence of two distinct subfamilies of eukaryotic rhomboids is supported by features of domain architectures that appear to comprise shared derived characters. Within these two major eukaryotic subfamilies, evolution apparently proceeded by both ancient and more recent duplications. Several lineage-specific expansions of paralogs  are noticeable, in insects, mammals and plants (Figure 2).
Archaeal rhomboids are scattered over the phylogenetic tree, with two major clusters and, in addition, three isolated proteins joining different bacterial branches (Figure 2). There is no indication of an affinity between any of the archaeal and eukaryotic rhomboids. Although many of the bacterial rhomboids form phylogenetically coherent clusters corresponding to the established bacterial lineages, there are also several clusters that have an odd composition, such as the grouping of proteobacterial and Gram-positive species; some of these clusters are well supported by bootstrap (see clusters 1-4 in Figure 2).
Log-likelihood analysis of possible placements of selected branches of maximum likelihood trees for the proteins analyzed
A → B
B → A
A → C
A → D
In addition, a tree of the rhomboid family was constructed using the Bayesian inference method, which has recently become a practical alternative to the more traditional methods of phylogenetic analysis [24, 25]. The tree produced using the MRBAYES package  showed the same major clades as the tree in Figure 2 (data not shown); moreover, clustering of the RHO and PARL subfamilies of eukaryotic rhomboids with the respective prokaryotic clades was supported by high posterior probabilities (Figure 2).
We also attempted to construct a phylogenetic tree of the rhomboid family by using the maximum parsimony method . The resulting tree contained the same major clades as the trees constructed using ML and MRBAYES; however, the number of parsimony-informative sites was insufficient to obtain high bootstrap support with this approach (data not shown).
Statistical comparisons of the best neighbor-joining tree with the hypothesis 1 and hypothesis 2 trees
Templeton (Wilcoxon signed-ranks) test
Winning-sites (sign) test
The concordance of the results obtained with several independent methods for phylogenetic tree construction and statistical analysis specifically aimed at testing the alternative hypothesis of monophyletic origin of eukaryotic rhomboids shows strong support for the major aspects of the tree topology in Figure 2 and, in particular, for the polyphyly of eukaryotic rhomboids.
The phylogenetic tree of the rhomboid family shown in Figure 2 and supported by the additional tests described above follows neither the 'standard model' scenario [28, 29], with the major split between the archaeo-eukaryotic and bacterial lineages nor the 'mitochondrial' scenario, which postulates acquisition of a gene by eukaryotes from the pro-mitochondrial endosymbiont. Neither can this tree be explained by postulating a small number of lineage-specific gene losses. The parsimonious interpretation of the rhomboid family tree seems to be that the evolutionary history of this family had been replete with horizontal gene transfer (HGT) and lineage-specific gene loss events. In particular, in spite of the presence of rhomboids in the majority of modern life forms from all three primary superkingdoms, phylogenetic analysis suggests that this family has not been inherited from the last universal common ancestor (LUCA). Instead, the tree topology seems to indicate that this family emerged in some bacterial lineage and afterwards had been widely disseminated by HGT, and then lost in some lineages. Both archaea and eukaryotes seem to have acquired rhomboids on several independent occasions. In particular, at least two HGT events seem to have contributed to the origin of eukaryotic rhomboids, one of them yielding the RHO subfamily and the other one the PARL subfamily, with a possible additional HGT in plants (Figures 2,3).
Given the broad phyletic representation of both subfamilies of eukaryotic rhomboids, both the RHO subfamily and the PARL subfamily must have been acquired through HGT at an early stage of eukaryotic evolution, definitely before the divergence of the major crown-group lineages. This early epoch in eukaryotic evolution is thought to have been dominated by HGT from multiple bacterial symbionts [30, 31].
An alternative to this multiple-HGT scenario is that LUCA already had multiple, paralogous rhomboids, which evolved by a series of ancient gene duplications, and the odd topology of the phylogenetic tree is due primarily to differential loss of these ancient paralogs. Although this cannot be ruled out formally, this hypothesis implies the existence of an elaborate signaling system in LUCA and, accordingly, suggests that LUCA was a complex organism, which might have had as many genes as modern bacteria. Theoretical analysis of evolutionary scenarios constructed on the basis of the phyletic patterns of COGs by applying the parsimony principle shows that the complexity of the inferred gene set of LUCA critically depends on the relative rates of gene loss and HGT at the early stages of evolution . A complex LUCA with around 2,000 genes is predicted only when one assumes that the rate of gene loss is an order of magnitude greater than the rate of HGT. However, explicit reconstruction of the gene set of LUCA under the assumption of equal rates of gene loss and HGT leads to a hypothetical genome that consists of only around 600 genes but appears to be 'compatible with life', that is, it includes genes responsible for most, if not all, essential cellular functions . We currently believe that this is the most realistic, albeit inevitably imprecise, reconstruction of LUCA's gene set. With respect to the rhomboid family and other families whose phylogenetic trees show similar patterns, this makes the multiple-HGT interpretation the scenario of choice. Further theoretical, comparative-genomic and experimental analyses aimed at determining relative rates of gene loss and HGT will help in a more objective assessment of the validity of this argument.
The rhomboid family might be the most widespread and conserved group of integral membrane proteins. In and by itself, this would suggest that this family is part of the gene repertoire of LUCA. However, phylogenetic analysis suggests a different scenario, one of emergence in a bacterial lineage with subsequent multiple, independent HGT events and gene losses. Although caution is due in the evolutionary interpretation of phylogenetic trees for large families, particularly when membrane proteins with a relatively small number of conserved positions, such as the rhomboids, are involved, the multiple-HGT scenario seemed to be supported by several methods of tree analysis and statistical tests.
Eukaryotes probably acquired their two major rhomboid subfamilies, RHO and PARL, as the result of two independent, early HGT events. These events, which might have introduced RIP as a means of intercellular communication, could have been pivotal in the evolution of eukaryotic multicellularity along the lines discussed previously with regard to the apparent bacterial origin of key components of eukaryotic programmed cell death machinery . Subsequent evolution of rhomboids in eukaryotes proceeded by lineage-specific expansion of paralogs  followed by diversification through the addition of an extra TMH in different positions relative to the catalytic core, some limited domain accretion (see Figure 2) and sequence divergence.
Phylogenetic analysis of the rhomboid family described here carries a general message for studies aimed at the reconstruction of ancestral life forms, particularly LUCA. Although most of the (nearly) ubiquitous protein families probably do derive from LUCA, explicit phylogenetic analysis is required to ascertain this in each case.
Materials and methods
The nonredundant (NR) protein sequence database at the National Center for Biotechnology Information (NIH, Bethesda) was searched iteratively using the PSI-BLAST program with multiple starting queries . PSI-BLAST was normally run with expectation (E) value of 0.01 as the cut-off for inclusion of sequences into the position-specific scoring matrix. Multiple alignments of protein sequences were constructed using the ClustalW program  and manually adjusted on the basis of the examination of PSI-BLAST search outputs and the superposition of the predicted TMHs, which were identified using the programs TMpred  and TMAP .
Phylogenetic trees were built using the least-squares method  implemented in the FITCH program of the PHYLIP package , with subsequent local rearrangement using the PROTML program of the MOLPHY package to obtain the maximum likelihood tree . The reliability of the tree topology was assessed using the RELL (resampling of estimated log-likelihoods) bootstrap method of MOLPHY, with 10,000 replications . Alternative placements of selected clades in maximum-likelihood trees were compared by using the rearrangement optimization method (Kishino-Hasegawa test) as implemented in the ProtML program [43–45]. Maximum parsimony trees were constructed using the heuristic search option of PAUP* . In addition, trees were constructed by Bayesian inference using the Markov chain Monte Carlo method as implemented in the MRBAYES package [24, 26]. The complete alignment information, including columns with gaps, was used for the MRBAYES analysis.
Constraint trees for phylogenetic hypothesis testing were generated using the TreeView program . Constraint trees were imported into PAUP*  and subjected to neighbor-joining search to generate the phylogenies corresponding to alternative hypotheses. These phylogenies were compared using the KH , Templeton (Wilcoxon signed-ranks)  and Winning-sites (sign)  tests implemented in PAUP*.
L.P. is supported by a grant from the Natural Sciences and Engineering Research Council of Canada.
- Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS, Kiryutin B, Galperin MY, Fedorova ND, Koonin EV: The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res. 2001, 29: 22-28. 10.1093/nar/29.1.22.PubMedPubMed CentralView ArticleGoogle Scholar
- Sturtevant MA, Roark M, Bier E: The Drosophila rhomboid gene mediates the localized formation of wing veins and interacts genetically with components of the EGF-R signaling pathway. Genes Dev. 1993, 7: 961-973.PubMedView ArticleGoogle Scholar
- Sturtevant MA, Roark M, O'Neill JW, Biehs B, Colley N, Bier E: The Drosophila rhomboid protein is concentrated in patches at the apical cell surface. Dev Biol. 1996, 174: 298-309. 10.1006/dbio.1996.0075.PubMedView ArticleGoogle Scholar
- Guichard A, Biehs B, Sturtevant MA, Wickline L, Chacko J, Howard K, Bier E: rhomboid and Star interact synergistically to promote EGFR/MAPK signaling during Drosophila wing vein development. Development. 1999, 126: 2663-2676.PubMedGoogle Scholar
- Mushegian AR, Koonin EV: Sequence analysis of eukaryotic developmental proteins: ancient and novel domains. Genetics. 1996, 144: 817-828.PubMedPubMed CentralGoogle Scholar
- Pellegrini L, Passer BJ, Canelles M, Lefterov I, Ganjei JK, Fowlkes BJ, Koonin EV, D'Adamio L: PAMP and PARL, two novel putative metalloproteases interacting with the COOH-terminus of Presenilin-1 and -2. J Alzheimers Dis. 2001, 3: 181-190.PubMedGoogle Scholar
- Urban S, Lee JR, Freeman M: Drosophila rhomboid-1 defines a family of putative intramembrane serine proteases. Cell. 2001, 107: 173-182.PubMedView ArticleGoogle Scholar
- Klambt C: EGF receptor signalling: roles of star and rhomboid revealed. Curr Biol. 2002, 12: R21-R23. 10.1016/S0960-9822(01)00642-X.PubMedView ArticleGoogle Scholar
- Guichard A, Roark M, Ronshaugen M, Bier E: brother of rhomboid, a rhomboid-related gene expressed during early Drosophila oogenesis, promotes EGF-R/MAPK signaling. Dev Biol. 2000, 226: 255-266. 10.1006/dbio.2000.9851.PubMedView ArticleGoogle Scholar
- Wasserman JD, Urban S, Freeman M: A family of rhomboid-like genes: Drosophila rhomboid-1 and roughoid/rhomboid-3 cooperate to activate EGF receptor signaling. Genes Dev. 2000, 14: 1651-1663.PubMedPubMed CentralGoogle Scholar
- Brown MS, Ye J, Rawson RB, Goldstein JL: Regulated intramembrane proteolysis: a control mechanism conserved from bacteria to humans. Cell. 2000, 100: 391-398.PubMedView ArticleGoogle Scholar
- Urban S, Freeman M: Intramembrane proteolysis controls diverse signalling pathways throughout evolution. Curr Opin Genet Dev. 2002, 12: 512-518. 10.1016/S0959-437X(02)00334-9.PubMedView ArticleGoogle Scholar
- Wolfe MS, Xia W, Ostaszewski BL, Diehl TS, Kimberly WT, Selkoe DJ: Two transmembrane aspartates in presenilin-1 required for presenilin endoproteolysis and gamma-secretase activity. Nature. 1999, 398: 513-517. 10.1038/19077.PubMedView ArticleGoogle Scholar
- Steiner H, Kostka M, Romig H, Basset G, Pesold B, Hardy J, Capell A, Meyn L, Grim ML, Baumeister R, et al: Glycine 384 is required for presenilin-1 function and is conserved in bacterial polytopic aspartyl proteases. Nat Cell Biol. 2000, 2: 848-851. 10.1038/35041097.PubMedView ArticleGoogle Scholar
- Sreekumar KR, Aravind L, Koonin EV: Computational analysis of human disease-associated genes and their protein products. Curr Opin Genet Dev. 2001, 11: 247-257. 10.1016/S0959-437X(00)00186-6.PubMedView ArticleGoogle Scholar
- Rather PN, Orosz E: Characterization of aarA, a pleiotrophic negative regulator of the 2'-N-acetyltransferase in Providencia stuartii. J Bacteriol. 1994, 176: 5140-5144.PubMedPubMed CentralGoogle Scholar
- Rather PN, Ding X, Baca-DeLancey RR, Siddiqui S: Providencia stuartii genes activated by cell-to-cell signaling and identification of a gene required for production or activity of an extracellular factor. J Bacteriol. 1999, 181: 7185-7191.PubMedPubMed CentralGoogle Scholar
- Gallio M, Kylsten P: Providencia may help find a function for a novel, widespread protein family. Curr Biol. 2000, 10: R693-R694. 10.1016/S0960-9822(00)00722-3.PubMedView ArticleGoogle Scholar
- Urban S, Schlieper D, Freeman M: Conservation of intramembrane proteolytic activity and substrate specificity in prokaryotic and eukaryotic rhomboids. Curr Biol. 2002, 12: 1507-1512. 10.1016/S0960-9822(02)01092-8.PubMedView ArticleGoogle Scholar
- Gallio M, Sturgill G, Rather P, Kylsten P: A conserved mechanism for extracellular signaling in eukaryotes and prokaryotes. Proc Natl Acad Sci USA. 2002, 99: 12208-12213. 10.1073/pnas.192138799.PubMedPubMed CentralView ArticleGoogle Scholar
- Esser K, Tursun B, Ingenhoven M, Michaelis G, Pratje E: A novel two-step mechanism for removal of a mitochondrial signal sequence involves the mAAA complex and the putative rhomboid protease Pcp1. J Mol Biol. 2002, 323: 835-843. 10.1016/S0022-2836(02)01000-8.PubMedView ArticleGoogle Scholar
- COGS: phylogenetic classification of proteins encoded in complete genomes. [http://www.ncbi.nlm.nih.gov/COG]
- Lespinet O, Wolf YI, Koonin EV, Aravind L: The role of lineage-specific gene family expansion in the evolution of eukaryotes. Genome Res. 2002, 12: 1048-1059. 10.1101/gr.174302.PubMedPubMed CentralView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001, 294: 2310-2314. 10.1126/science.1065889.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Larget B, Miller RE, Ronquist F: Potential applications and pitfalls of bayesian inference of phylogeny. Syst Biol. 2002, 51: 673-688. 10.1080/10635150290102366.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and other methods). 1998, Sunderland, MA: SinauerGoogle Scholar
- Brown JR, Doolittle WF: Archaea and the prokaryote-to-eukaryote transition. Microbiol Mol Biol Rev. 1997, 61: 456-502.PubMedPubMed CentralGoogle Scholar
- Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA. 1990, 87: 4576-4579.PubMedPubMed CentralView ArticleGoogle Scholar
- Doolittle WF: You are what you eat: a gene transfer ratchet could account for bacterial genes in eukaryotic nuclear genomes. Trends Genet. 1998, 14: 307-311. 10.1016/S0168-9525(98)01494-2.PubMedView ArticleGoogle Scholar
- Koonin EV, Galperin MY: Sequence - Evolution - Function. Computational Approaches in Comparative Genomics. 2002, Boston: KluwerGoogle Scholar
- Mirkin BG, Fenner TI, Galperin MY, Koonin EV: Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes. BMC Evol Biol. 2003, 3: 2-10.1186/1471-2148-3-2.PubMedPubMed CentralView ArticleGoogle Scholar
- Snel B, Bork P, Huynen MA: Genomes in flux: the evolution of archaeal and proteobacterial gene content. Genome Res. 2002, 12: 17-25. 10.1101/gr.176501.PubMedView ArticleGoogle Scholar
- Gogarten JP, Doolittle WF, Lawrence JG: Prokaryotic evolution in light of gene transfer. Mol Biol Evol. 2002, 19: 2226-2238.PubMedView ArticleGoogle Scholar
- Weihofen A, Binns K, Lemberg MK, Ashman K, Martoglio B: Identification of signal peptide peptidase, a presenilin-type aspartic protease. Science. 2002, 296: 2215-2218. 10.1126/science.1070925.PubMedView ArticleGoogle Scholar
- Koonin EV, Aravind L: Origin and evolution of eukaryotic apoptosis: the bacterial connection. Cell Death Differ. 2002, 9: 394-404. 10.1038/sj.cdd.4400991.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.PubMedPubMed CentralView ArticleGoogle Scholar
- Hofmann K, Stoffel W: TMbase - A database of membrane spanning protein segments. Biol Chem Hoppe-Seyler. 1993, 374: 166-Google Scholar
- Persson B, Argos P: Prediction of membrane protein topology utilizing multiple sequence alignments. J Protein Chem. 1997, 16: 453-457. 10.1023/A:1026353225758.PubMedView ArticleGoogle Scholar
- Fitch WM, Margoliash E: Construction of phylogenetic trees. Science. 1967, 155: 279-284.PubMedView ArticleGoogle Scholar
- Felsenstein J: Inferring phylogenies from protein sequences by parsimony, distance, and likelihood methods. Methods Enzymol. 1996, 266: 418-427.PubMedView ArticleGoogle Scholar
- Adachi J, Hasegawa M: MOLPHY: Programs for Molecular Phylogenetics. 1992, Tokyo: Institute of Statistical MathematicsGoogle Scholar
- Kishino H, Miyata T, Hasegawa M: Maximum likelihood inference of protein phylogeny and the origin of chloroplasts. J Mol Evol. 1990, 31: 151-160.View ArticleGoogle Scholar
- Kishino H, Hasegawa M: Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea. J Mol Evol. 1989, 29: 170-179.PubMedView ArticleGoogle Scholar
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
- Templeton AR: Phylogenetic inference from restriction endonuclease cleavage site maps with particular reference to the humans and apes. Evolution. 1983, 37: 221-244.View ArticleGoogle Scholar
- Prager EM, Wilson AC: Ancient origin of lactalbumin from lysozyme: analysis of DNA and amino acid sequences. J Mol Evol. 1988, 27: 326-335.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.