A genomic timescale for the origin of eukaryotes
© Hedges et al; licensee BioMed Central Ltd. 2001
Received: 21 July 2001
Accepted: 12 September 2001
Published: 12 September 2001
Genomic sequence analyses have shown that horizontal gene transfer occurred during the origin of eukaryotes as a consequence of symbiosis. However, details of the timing and number of symbiotic events are unclear. A timescale for the early evolution of eukaryotes would help to better understand the relationship between these biological events and changes in Earth's environment, such as the rise in oxygen. We used refined methods of sequence alignment, site selection, and time estimation to address these questions with protein sequences from complete genomes of prokaryotes and eukaryotes.
Eukaryotes were found to evolve faster than prokaryotes, with those eukaryotes derived from eubacteria evolving faster than those derived from archaebacteria. We found an early time of divergence (~4 billion years ago, Ga) for archaebacteria and the archaebacterial genes in eukaryotes. Our analyses support at least two horizontal gene transfer events in the origin of eukaryotes, at 2.7 Ga and 1.8 Ga. Time estimates for the origin of cyanobacteria (2.6 Ga) and the divergence of an early-branching eukaryote that lacks mitochondria (Giardia) (2.2 Ga) fall between those two events.
We find support for two symbiotic events in the origin of eukaryotes: one premitochondrial and a later mitochondrial event. The appearance of cyanobacteria immediately prior to the earliest undisputed evidence for the presence of oxygen (2.4–2.2 Ga) suggests that the innovation of oxygenic photosynthesis had a relatively rapid impact on the environment as it set the stage for further evolution of the eukaryotic cell.
An emerging pattern found in gene and protein phylogenies that include prokaryotes (archaebacteria and eubacteria) and eukaryotes is the variable position of eukaryotes. In proteins involved in transcription and translation, eukaryotes often cluster with archaebacteria whereas in metabolic proteins they often cluster with eubacteria . Among the latter proteins, eukaryotes sometimes group with α-proteobacteria, presumably reflecting the origin of mitochondria, and plants sometimes cluster with cyanobacteria, reflecting the origin of plastids. These patterns have been interpreted as a general signature of the symbiotic origin of eukaryotes [2, 3] and horizontal gene transfer (HGT) of symbiont genes to the nucleus [4–9]. On the one hand, this complexity resulting from HGT can obscure some aspects of evolutionary history . However, HGT also can provide the means to investigate otherwise difficult questions, such as inferring the number of symbiotic events and estimating the time of those events. This is the approach that we take in this study.
The use of sequence changes to estimate the time of these early events also has its assumptions and limitations [21–23]. Nonetheless, many proteins contain conserved regions of amino acid sequence throughout prokaryotes and eukaryotes that permit alignment and analysis. The most extensive of these analyses have found that all major events related to the origin of eukaryotes occurred about 2.0–2.2 Ga [5, 21]. This includes the divergence of archaebacteria and archaebacterial proteins in eukaryotes, the origin of cyanobacteria, and the divergence of eubacteria and eubacterial proteins in eukaryotes (the latter presumably reflecting symbiosis). However, these times were not adjusted for lineage-specific rate differences that have been discovered subsequently . Here, we estimate the time of these events with protein sequences from complete genomes and consideration of lineage-specific rate variation.
Two other rate comparisons were limited by a small number of proteins: eubacteria versus eukaryotes (KA) and eubacteria versus archaebacteria. Only three proteins were available in the first comparison and all three showed a faster rate in eukaryotes (1.43, 1.12, 1.23; x = 1.26). This result differs from that reported elsewhere , in which the two rates were not significantly different. In the second case, we found that archaebacteria are evolving at a slower rate than eubacteria, as was noted elsewhere . In our case, regression of archaebacterial branch length versus eubacterial branch length, fixed through the origin, resulted in a slope of 0.93 and correlation coefficient of 0.65 (n = 9 proteins). However, in both of these comparisons, rate tests did not yield significant rate differences probably because of the short length of most proteins. Sample size (eight protein sets) also was limited in the Kollman and Doolittle study . Taken together these data suggest the following relative order of rate differences: archaebacteria < eubacteria < eukaryotes (archaebacterial origin) < eukaryotes (eubacterial origin). As additional genomic data become available, more proteins will be useful and greater precision in these rates and rate differences will be possible.
Phylogeny and time estimation
It has been suggested that eukaryotic genes and proteins of archaebacterial origin are more closely related to one lineage of archaebacteria (Crenarchaeota; "eocytes") than the other major lineage (Euryarchaeota) . If true, this would bear on our time estimate for the divergence of archaebacteria and eukaryotes. Thus, we conducted a phylogenetic analysis of 72 proteins containing representatives of the two major groups of archaebacteria, eukaryotes, and eubacteria. At the 95% bootstrap significance level, 19 proteins supported archaebacterial monophyly whereas none supported the eocyte hypothesis (Crenarchaeota + Eukaryota). This indicates that the lineage of archaebacteria leading to the eukaryote nuclear genome diverged prior to the split between the Crenarchaeota and Euryarchaeota. As noted previously , most (in this case, 21 out of 36) eukaryotic proteins with archaebacterial affinity are informational (involved in transcription, translation, and related processes).
Divergence time estimates (billion years ago)
Mean ± SE*
3.50 ± 0.25
3.97 ± 0.32
2.46 ± 0.14
2.73 ± 0.20
1.79 ± 0.29
2.56 ± 0.26
2.50 ± 0.22
2.23 ± 0.12
1.56 ± 0.29
1.84 ± 0.20
The purpose of this study was to examine the temporal relationship between the origin of eukaryotes and events in Earth history. However, some unexpected results required refinement in methodology. These included finding greater among-site rate variation in the calibration group and different rates of sequence change between prokaryotes and eukaryotes, and between eukaryotes derived from different groups of prokaryotes. By taking into account these variables, the resulting time estimates are more robust and have fewer assumptions. For example, the time estimate for the origin of eukaryotes (BK-o) is not based on a general assumption of rate constancy between prokaryotes (or even eubacteria) and eukaryotes because rates are adjusted for each protein and each comparison. Also, the calibration used for BK-o is not a general eukaryotic calibration but one based exclusively on eukaryote sequences derived from eubacteria. A tradeoff in these improved methods was a reduction in the number of proteins that could be used, which increased the variance of the time estimates. Nonetheless, the phylogenies and time estimates obtained in this study have a bearing on current models for the evolution of eukaryotes.
Until about five years ago, it was generally accepted that there was a prior period (before mitochondria) in the history of eukaryotes [2, 26]. The basal position of eukaryotes lacking mitochondria (amitochondriate) in phylogenetic trees  was consistent with this supposition as was evidence from sequence signatures . However, molecular phylogenetic studies of several proteins in recent years have suggested that some or all amitochondriate eukaryotes once possessed mitochondria in the past . Based on this new evidence, most current models for the origin of eukaryotes assume only a single symbiotic or fusion event between an archaebacterium and an α-proteobacterium [8, 28, 29].
Under the single-symbiosis model, eukaryotes should cluster exclusively with an α-proteobacterium (e.g., Rickettsia), among eubacteria. However, our phylogenetic analyses (Fig. 4) instead indicate, significantly, that many eukaryotic proteins originated from one (or more) eubacterial lineages other than α-proteobacteria. The reduced genome of Rickettsia would not explain this result because Rickettsia possesses all of the proteins used in the combined analysis (Fig. 4B). Protein function and location also are consistent with a premitochondrial origin. Only one of the 32 BK-o proteins is restricted to the mitochondrion whereas eight of the nine BK-m proteins are restricted to that organelle. Also, all six of the proteins involved in cellular respiration are in the BK-m group. Based on the serial endosymbiosis theory, the first symbiotic event involved a spirochete . On the other hand, sequence signatures of the heat shock molecular chaperone protein HSP-70 and other evidence have indicated that the first symbiotic event involved a gram-negative eubacterium . Our data are unable to distinguish between these two alternatives but agree with both in implicating an earlier, premitochondrial event. Predation by prokaryotes on early eukaryotes also may have led to HGT.
If two or more symbiotic events were involved, this does not necessarily confirm that any of the living lineages of amitochondriate eukaryotes arose prior to the second (mitochondrial) event. All may have once possessed mitochondria. However, because Giardia arose at an early time (Table 1) and branches near the base of the eukaryote phylogeny, the simplest explanation is that it never possessed mitochondria and is a primary (not secondary) amitochondriate. Although the position of Giardia in some protein phylogenies  has been proposed as evidence that it is a secondary amitochondriate, others have urged caution until additional, more conclusive, data become available .
Extensive glaciations occurred in the Paleoproterozoic (~2.4 Ga), and may have been global in extent . It has been proposed that a major rise in oxygen at this time lowered global temperatures and may have triggered the glaciations . If this is true, and given the time estimates here, the evolutionary innovation of oxygenic photosynthesis may have had a relatively rapid impact on the environment. Moreover, this innovation may have caused a mass extinction of prokaryotes at that time, as a result of the toxic effects of oxygen, as suggested by the virtual absence of lineages prior to ~2.5 Ga and subsequent rapid radiation of lineages (Figs. 4,5).
Our analyses of prokaryotic and eukaryotic genomic sequence data support two symbiotic events in the origin of eukaryotes: one premitochondrial (2.7 billion years ago, Ga) and a later mitochondrial event (1.8 Ga). Our time estimate for the divergence of an early-branching eukaryote (Giardia) that lacks mitochondria, 2.2 Ga, suggests that it is a primary and not secondary amitochondriate organism. Our time estimate for the origin of cyanobacteria (2.6 Ga) is more recent than expected and suggests that earlier fossils claimed to be of cyanobacteria are of other organisms (or artifacts). Moreover, the appearance of cyanobacteria immediately prior to the earliest undisputed evidence for the presence of oxygen (2.4–2.2 Ga) suggests that the innovation of oxygenic photosynthesis had a relatively rapid impact on the environment as it set the stage for further evolution of the eukaryotic cell.
Materials and Methods
Sequence data and alignment
We assembled and aligned protein sequences of all 467 potentially orthologous groups from complete genome databases and these were supplemented with additional eukaryote taxa from the sequence database of the National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov/entrez/. Data from the following species were assembled into presumptive orthology groups (hereafter, proteins) and aligned : Aquifex aeolicus, Bacillus subtilis, Borrelia burgdorferi, Chlamydia trachomatis, Escherichia coli, Haemophilus influenzae, Helicobacter pylori, Mycobacterium tuberculosis, Mycoplasma genitalium, Mycoplasma pneumoniae, Rickettsia prowazekii, Synechocystis sp., and Treponema pallidum (Eubacteria), Archaeoglobus fulgidus, Methanobacterium thermoautotrophicum, Methanococcus jannaschii, Pyrococcus abyssi, and P. horikoshii (Archaebacteria: Euryarchaeota), and Arabidopsis thaliana, Caenorhabditis elegans, Dictyostelium discoideum, Drosophila melanogaster, Gallus gallus, Homo sapiens, Mus musculus, Oryza saticva, Plasmodium falciparum, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Trypanosoma sp., and Xenopus laevis (Eukaryota). The genome sequence of Aeropyrum pernix (Archaebacteria: Crenarchaeota) became available during completion of study and was included in phylogenetic analyses only.
Global alignment algorithms differ from local alignment algorithms in that they sometimes align unrelated (non-homologous) sites together with homologous sites. Using a computational tool, xcons , such unrelated sites were removed from these CLUSTALW  alignments to increase probability of site homology. During construction of protein alignments using the WAT system , short fragmented sequences were manually removed. Of the 204 proteins that could be calibrated for time estimation, the orthology of roughly half (116 proteins) was ambiguous for unknown reasons (e.g., lateral gene transfer, gene loss, or poor phylogenetic resolution) leaving 87 proteins for phylogeny and time estimation. The seven shortest (<75 amino acids) of those were used only in phylogenetic analyses; the remaining proteins averaged 196 amino acids each. Where possible, proteins were rooted by duplicate proteins (duplicate genes); otherwise, they were midpoint-rooted.
Separately, for timing the origin of Giardia, sequences of 17 proteins were obtained from the public databases and aligned  in which the following taxa were available: Giardia and other eukaryotes (including calibration taxa; see below), archaebacteria, and eubacteria. Correspondence and requests for materials should be addressed to S.B.H. (e-mail: firstname.lastname@example.org) or see http://www.evogenomics.org/Publications/data/Eukaryotes/ for alignments and other information.
Methods are described elsewhere  except as follows. Our initial goal was to estimate divergence times for the last common ancestor (LCA), the divergence between archaebacteria and eukaryotes (AK), cyanobacteria and closest eubacterial relatives (origin of cyanobacteria, BC), eubacteria and mitochondrial eukaryotes (origin of mitochondria, BK-m), and Giardia and other eukaryotes (GK) (Fig. 1). The importance of Giardia is its lack of mitochondria and basal location in many phylogenies of eukaryotes [27, 39].
However, our initial phylogenetic analyses revealed that many eukaryotes did not cluster with Rickettsia, the α-proteobacterium, as predicted by current genomic models [8, 25]. Instead, they typically formed a basal lineage among eubacteria in the tree. This result was consistent with the serial endosymbiosis theory  and with other findings  and therefore we designated this divergence as BK-o (origin of eukaryotes). Estimation of the divergence time of the origin of plastids (BK-p) was not a goal of this study, and the LCA was not estimated because of an insufficient number of duplicate proteins needed for reciprocal rooting . Thus, five divergence times were studied: AK, BC, BK-o, BK-m, and GK. Eukaryotes derived from different prokaryotes are referred to herein as KA (from AK), KB-o (from BK-o), and KB-m (from BK-m).
Because of the large amount of sequence conservation in these proteins, it was not possible to calibrate directly by extrapolation from vertebrates , for which an extensive fossil record exists. For example, sequences often were identical among rodents, primates, and birds. Instead, multiple calibrations were used from older divergences among kingdoms (plants, animals, fungi) and animal phyla, derived from analysis of 75 nuclear proteins calibrated with the vertebrate fossil record . This two-step calibration reduced the error involved in extrapolation. Two classes of time estimation methods were used and compared. The multigene (MG) approach uses the mean or mode of many single-gene time estimates [40, 41] whereas the average-distance (AD) approach [42–44] involves the combining of distances and rates among genes or proteins to yield a single time estimate. For the AD approach, we weight each single-gene distance, before combining, by the length of the protein (aligned amino acids).
Protein-specific rates were estimated by regression, fixed through the origin, of these calibration points within eukaryotes: arthropod-chordate (0.993 Ga), chordate-nematode (1.177 Ga), and plant-animal-fungi (1.576 Ga) . During the course of the study, it was discovered that the shape parameter (α) of the gamma distribution used to account for rate variation among sites, estimated by a likelihood method , differed consistently between calibration taxa (average, 1.99) and the overall data set (1.44) for each gene (Fig. 2). Therefore, a dual-gamma approach was taken whereby the eukaryote rate was estimated using the eukaryote gamma parameter and the time estimate (involving prokaryotes) was made using the overall gamma parameter. There is insufficient evidence at present to determine whether or not this difference is biologically based, related to the covarion model , or follows a simple scaling relationship with time or total protein distance. If the relationship is scaled, additional modification in methods may be necessary in the future.
We compared rates of change in archaebacterial versus eubacterial sequences using paralogous sequences (those related by gene duplication) as a root for relative rate tests [47, 48]. To examine rate differences between eukaryotic sequences and their closest prokaryote orthologs (those representing the same gene), we used the more distant prokaryote (archaebacteria or eubacteria) as root. For examining rates in eukaryotic sequences derived from either archaebacteria or eubacteria, we compared pairwise distances of the same taxa present in both locations (e.g., one pair clustering with archaebacteria and the other with eubacteria) in the same protein. The discovery of rate differences among prokaryotes and eukaryotes required rate adjustments for all proteins and comparisons, including those accepted in rate tests. These adjustments were made by estimating time only with the eukaryote lineage, or in the case of BC, using a cyanobacterial rate adjusted by direct comparison of the cyanobacteria branch and eukaryote branch in rate tests. For example, the AK divergence time was estimated only with the KA calibration and the BC, BK-o, and BK-m divergence times were estimated only with the KB-o calibration. These restrictions further reduced the number of proteins available for time estimation to the following: AK (36 total, 21 constant-rate), BK-o (25, 16), BK-m (7, 5), BC (20, 16), and Giardia-eukaryotes (17, 11).
Modes were used in the MG approach, as described previously , except with the BK-m comparison where the median was used because of the small number of proteins. The mode is preferred over the mean or median because it eliminates or reduces the effect of outliers (e.g., unusually high estimates resulting from paralogous comparisons). In this study, a large number of overlapping bins was used initially to better define the distribution of time estimates, followed by use of a smaller number of non-overlapping bins and standard estimation of mode by interpolation. This two-step procedure was found to reduce the influence of bin size on mode estimation. For the AD approach, single-gene distances and rates were weighted by sequence length and then combined distances were divided by combined rates.
Phylogenetic trees  were constructed for each gene from amino acid data to assist in gene selection and interpretation. A gamma distance was used for all trees, with α estimated from the entire data set . An analysis involving combined protein alignments was performed with maximum likelihood , neighbor joining [49, 51], and maximum parsimony , using bootstrapping. Bootstrap consensus trees show branch-lengths estimated by ordinary least-squares method . Bootstrap support ≥ 95% was considered significant.
We thank L. E. Brahmakulam, B. R. Eidell, D. S. Heckman, A. R. Pfaff, J. L. Shoe, and R. L. Stauffer for assistance with analyses, and D. J. Des Marais, J. Kasting, L. Margulis, W. Martin, M. Nei, and M. Rohmer for comments or discussion. S.B.H. was supported by NASA and NSF, and S.K. by NIH and NSF.
- Rivera MC, Jain R, Moore JE, Lake JA: Genomic evidence for two functionally distinct gene classes. Proc Natl Acad Sci USA. 1998, 95: 6239-6244. 10.1073/pnas.95.11.6239.PubMed CentralView ArticlePubMedGoogle Scholar
- Margulis L: Origin of Eukaryotic Cells. New Haven, Connecticut, Yale University Press. 1970Google Scholar
- Margulis L: Archael-eubacterial mergers in the origin of Eukarya: Phylogenetic classification of life. Proc Natl Acad Sci USA. 1996, 93: 1071-1076. 10.1073/pnas.93.3.1071.PubMed CentralView ArticlePubMedGoogle Scholar
- Golding GB, Gupta R: Protein-based phylogenies support a chimeric origin of the eukaryotic genome. Mol Biol Evol. 1995, 12: 1-6.View ArticlePubMedGoogle Scholar
- Feng D-F, Cho G, Doolittle RF: Determining divergence times with a protein clock: Update and reevaluation. Proc Natl Acad Sci USA. 1997, 94: 13028-13033. 10.1073/pnas.94.24.13028.PubMed CentralView ArticlePubMedGoogle Scholar
- Gupta RS: Protein phylogenies and signature sequences: A reappraisal of evolutionary relationships among archaebacteria, eubacteria, and eukaryotes. Micro Mol Biol Rev. 1998, 62: 1435-1491.Google Scholar
- Doolittle WF: A paradigm gets shifty. Nature. 1998, 392: 15-16. 10.1038/32033.View ArticlePubMedGoogle Scholar
- Doolittle WF: Phylogenetic classification and the universal tree. Science. 1999, 284: 2124-2128. 10.1126/science.284.5423.2124.View ArticlePubMedGoogle Scholar
- Lang DF, Gray MW, Burger G: Mitochondrial genome evolution and the origin of eukaryotes. Annu Rev Genet. 1999, 33: 351-397. 10.1146/annurev.genet.33.1.351.View ArticlePubMedGoogle Scholar
- DesMarais DJ: When did photosynthesis emerge on earth?. Science. 2000, 289: 1703-1705.Google Scholar
- Brocks JJ, Logan GA, Buick R, Summons RE: Archean molecular fossils and the early rise of eukaryotes. Science. 1999, 285: 1033-1036. 10.1126/science.285.5430.1033.View ArticlePubMedGoogle Scholar
- Han T-M, Runnegar B: Megascopic eukaryotic algae from the 2.1 billion-year-old Negaunee iron-formation, Michigan. Science. 1992, 257: 232-235.View ArticlePubMedGoogle Scholar
- Schopf JW: Microfossils of the early Archean apex chert: New evidence of the antiquity of life. Science. 1993, 260: 640-646.View ArticlePubMedGoogle Scholar
- Knoll AH: A new molecular window on early life. Science. 1999, 285: 1025-1027. 10.1126/science.285.5430.1025.View ArticlePubMedGoogle Scholar
- Summons RE, Jahnke LL, Hope JM, Logan GA: 2-methylhopanoids as biomarkers for cyanobacterial oxygenic photosynthesis. Nature. 1999, 400: 554-557. 10.1038/23005.View ArticlePubMedGoogle Scholar
- Rohmer M, Bisseret P, Neunlist S: The hopanoids, prokaryotic triterpenoids and precursors of ubiquitous molecular fossils. In: Biological Markers in Sediments and Petroleum : A Tribute to Wolfgang K Seifert (Edited by Moldowan JM, Albrecht P, Philp RP) Englewood Cliffs, New Jersey, Prentice Hall. 1992, 1-17.Google Scholar
- Rohmer M: The biosynthesis of triterpenoids of the hopane series in the Eubacteria: A mine of new enzyme reactions. Pure Appl Chem. 1993, 65: 1293-1298.View ArticleGoogle Scholar
- Ohmoto H: When did the earth's atmosphere become oxic?. The Geochemical News. 1997, 93: 12-26.Google Scholar
- Watanabe Y, Martini JEJ, Ohmoto H: Geochemical evidence for terrestrial ecosystems 2.6 billion years ago. Nature. 2000, 408: 574-578. 10.1038/35046052.View ArticlePubMedGoogle Scholar
- Widdel F, Schnell S, Heising S, Ehrenreich A, Assmus B: Ferrous iron oxidation by anoxygenic phototrophic bacteria. Nature. 1993, 362: 834-838. 10.1038/362834a0.View ArticleGoogle Scholar
- Doolittle RF, Feng D-F, Tsang S, Cho G, Little E: Determining divergence times of the major kingdoms of living organisms with a protein clock. Science. 1996, 271: 470-477.View ArticlePubMedGoogle Scholar
- Hasegawa M, Fitch WM: Dating the cenancestor of organisms. Science. 1996, 274: 1750-10.1126/science.274.5293.1750.View ArticlePubMedGoogle Scholar
- Kollman JM, Doolittle RF: Determining the relative rates of change for prokaryotic and eukaryotic proteins with anciently duplicated paralogs. J Mol Evol. 2000, 51: 173-181.PubMedGoogle Scholar
- Rivera MC, Lake JA: Evidence that eukaryotes and eocyte prokaryotes are immediate relatives. Science. 1992, 257: 74-76.View ArticlePubMedGoogle Scholar
- Andersson SGE, Zomorodipour A, Andersson JO, Sicheritz-Ponten T, Alsmark UCM, Podowski RM, Naslund AK, Eriksson A-S, Winkler HH, Kurland CG: The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature. 1998, 396: 133-140. 10.1038/24094.View ArticlePubMedGoogle Scholar
- Cavalier-Smith T: Kingdom protozoa and its 18 phyla. Microbiol Rev. 1993, 57: 953-994.PubMed CentralPubMedGoogle Scholar
- Sogin ML, Gunderson JH, Elwood HJ, Alonso RA, Peattie DA: Phylogenetic meaning of the kingdom concept: An unusual ribosomal RNA from G iardia lamblia. Science. 1989, 243: 75-77.View ArticlePubMedGoogle Scholar
- Embley TM, Hirt RP: Early branching eukaryotes?. Curr Opin Genet Dev. 1998, 8: 624-629. 10.1016/S0959-437X(98)80029-4.View ArticlePubMedGoogle Scholar
- Martin W, Müller M: The hydrogen hypothesis for the first eukaryote. Nature. 1998, 392: 37-41. 10.1038/32096.View ArticlePubMedGoogle Scholar
- Roger AJ, Svard SG, Tovar J, Clark CG, Smith MW, Gillin FD, Sogin ML: A mitochondrial-like chaperonin 60 gene in G iardia lamblia : Evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria. Proc Natl Acad Sci USA. 1998, 95: 229-234. 10.1073/pnas.95.1.229.PubMed CentralView ArticlePubMedGoogle Scholar
- Mojzsis SJ, Arrhenius G, McKeegan KD, Harrison TM, Nutman AP, Friend CRL: Evidence for life on Earth before 3,800 million years ago. Nature. 1996, 384: 55-59. 10.1038/384055a0.View ArticlePubMedGoogle Scholar
- Gupta RS, Mukhtar T, Singh B: Evolutionary relationships among photosynthetic prokaryotes (Heliobacterium chlorum, Chloroflexus aurantiacus, cyanobacteria, Chlorobium tepidum and proteobacteria): Implications regarding the origin of photosynthesis. Mol Microbiol. 1999, 32: 893-906. 10.1046/j.1365-2958.1999.01417.x.View ArticlePubMedGoogle Scholar
- Xiong J, Fischer WM, Inoue K, Nakahara M, Bauer CE: Molecular evidence for the early evolution of photosynthesis. Science. 2000, 289: 1724-1730. 10.1126/science.289.5485.1724.View ArticlePubMedGoogle Scholar
- Kirschvink JL, Gaidos EJ, Bertani LE, Beukes NJ, Gutzmer J, Maepa LN, Steinberger RE: Paleoproterozoic snowball Earth: Extreme climatic and geochemical global change and its biological consequences. Proc Natl Acad Sci USA. 2000, 97: 1400-1405. 10.1073/pnas.97.4.1400.PubMed CentralView ArticlePubMedGoogle Scholar
- Pavlov AA, Kasting JF, Brown LL: Greenhouse warming by CH4 in the atmosphere of early Earth. J Geophys Res. 2000, 105: 11981-11990. 10.1029/1999JE001134.View ArticlePubMedGoogle Scholar
- Shigenobu S, Watanabe H, Hattori M, Sakaki Y, Ishikawa H: The genome sequence of Buchnera sp. APS, an endocellular bacterial symbiont of aphids. Nature. 2000, 407: 81-86. 10.1016/S0168-9002(97)01373-9.View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: ClustalW: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang DY-C, Kumar S, Hedges SB: Divergence time estimates for the early history of animal phyla and the origin of plants, animals, and fungi. Proc R Soc Lond B. 1999, 266: 163-171. 10.1098/rspb.1999.0617.View ArticleGoogle Scholar
- Sogin ML: History assignment: When was the mitochondrion founded?. Curr Opin Genet Dev. 1997, 7: 792-799. 10.1016/S0959-437X(97)80042-1.View ArticlePubMedGoogle Scholar
- Kumar S, Hedges SB: A molecular timescale for vertebrate evolution. Nature. 1998, 392: 917-920. 10.1038/31927.View ArticlePubMedGoogle Scholar
- Hedges SB, Parker PH, Sibley CG, Kumar S: Continental breakup and the ordinal diversification of birds and mammals. Nature. 1996, 381: 226-229. 10.1038/381226a0.View ArticlePubMedGoogle Scholar
- Fitch WM: Molecular evolutionary clocks. In: Molecular Evolution (Edited by Ayala FJ) Sunderland, MA, Sinauer Associates. 1976, 160-178.Google Scholar
- Gu X: Early metazoan divergence was about 830 million years ago. J Mol Evol. 1998, 47: 369-371.View ArticlePubMedGoogle Scholar
- Lynch M: The age and relationships of the major animal phyla. Evolution. 1999, 53: 319-325.View ArticleGoogle Scholar
- Yang Z: Paml: A program package for phylogenetic analysis by maximum likelihood. CABIOS. 1997, 13: 555-556.PubMedGoogle Scholar
- Fitch WM, Markowitz E: An improved method for determining codon variability in a gene and its application to the rate of fixation of mutations in evolution. Biochem Genet. 1970, 4: 579-593.View ArticlePubMedGoogle Scholar
- Takezaki N, Rzhetsky A, Nei M: Phylogenetic test of the molecular clock and linearized trees. Mol Biol Evol. 1995, 12: 823-833.PubMedGoogle Scholar
- Kumar S: Phyltest: A program for testing phylogenetic hypotheses. University Park, Pennsylvania, Institute of Molecular Evolutionary Genetics, Pennsylvania State University. 1996Google Scholar
- Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Adachi J, Hasegawa M: Molphy version 2.3: Programs for molecular phylogenetics based on maximum likelihood. Tokyo, Institute of Statistical Mathematics. 1996Google Scholar
- Kumar S, Tamura K, Nei M: Mega: Molecular evolutionary genetic analysis. University Park, Pennsylvania State University. 1993Google Scholar
- Holland HD: Early proterozoic atmosphere change. In: Early life on Earth (Edited by Bengston S) New York, Columbia University Press. 1994, 237-244.Google Scholar
- Ohmoto H: Evidence in pre – 2.2 Ga paleosols for the early evolution of atmospheric oxygen and terrestrial biota. Geology. 1996, 24: 1135-1135. 10.1130/0091-7613(1996)024<1135:EIPGPF>2.3.CO;2.View ArticlePubMedGoogle 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.