29+ Evidences for Macroevolution

Part 1:
The Unique Universal Phylogenetic Tree

Copyright © 1999-2013 by Douglas Theobald, Ph.D.

"As buds give rise by growth to fresh buds, and these, if vigorous, branch out and overtop on all sides many a feebler branch, so by generation I believe it has been with the great Tree of Life, which fills with its dead and broken branches the crust of the earth, and covers the surface with its ever-branching and beautiful ramifications."

Charles Darwin
The Origin of Species, p. 171

Part 1 Outline

  1. The fundamental unity of life
  2. A nested hierarchy of species
  3. Independent determination of the historical phylogeny
  4. Intermediate and transitional forms
  5. Chronology of common ancestors

Prediction 1.1: The fundamental unity of life

"OH JEHOVA, quam ampla sunt opera Tua."
– Carolus Linnaeus
at the beginning of Systema Naturae, 1757

According to the theory of common descent, modern living organisms, with all their incredible differences, are the progeny of one single species in the distant past. In spite of the extensive variation of form and function among organisms, several fundamental criteria characterize all life. Some of the macroscopic properties that characterize all of life are (1) replication, (2) heritability (characteristics of descendents are correlated with those of ancestors), (3) catalysis, and (4) energy utilization (metabolism). At a very minimum, these four functions are required to generate a physical historical process that can be described by a phylogenetic tree.

If every living species descended from an original species that had these four obligate functions, then all living species today should necessarily have these functions (a somewhat trivial conclusion). Most importantly, however, all modern species should have inherited the structures that perform these functions. Thus, a basic prediction of the genealogical relatedness of all life, combined with the constraint of gradualism, is that organisms should be very similar in the particular mechanisms and structures that execute these four basic life processes.


The common polymers of life

The structures that all known organisms use to perform these four basic processes are all quite similar, in spite of the odds. All known living things use polymers to perform these four basic functions. Organic chemists have synthesized hundreds of different polymers, yet the only ones used by life, irrespective of species, are polynucleotides, polypeptides, and polysaccharides. Regardless of the species, the DNA, RNA and proteins used in known living systems all have the same chirality, even though there are at least two chemically equivalent choices of chirality for each of these molecules. For example, RNA has four chiral centers in its ribose ring, which means that it has 16 possible stereoisomers—but only one of these stereoisomers is found in the RNA of known living organisms.

Nucleic acids are the genetic material of life

Ten years after the publication of The Origin of Species, nucleic acids were first isolated by Friedrich Miescher in 1869. It took another 75 years after this discovery before DNA was identified as the genetic material of life (Avery et al. 1944). It is quite conceivable that we could have found a different genetic material for each species. In fact, it is still possible that newly identified species might have unknown genetic materials. However, all known life uses the same polymer, polynucleotide (DNA or RNA), for storing species specific information. All known organisms base replication on the duplication of this molecule. The DNA used by living organisms is synthesized using only four nucleosides (deoxyadenosine, deoxythymidine, deoxycytidine, and deoxyguanosine) out of the dozens known (at least 102 occur naturally and many more have been artificially synthesized) (Rozenski et al. 1999; Voet and Voet 1995, p. 969).

Protein catalysis

In order to perform the functions necessary for life, organisms must catalyze chemical reactions. In all known organisms, enzymatic catalysis is based on the abilities provided by protein molecules (and in relatively rare, yet important, cases by RNA molecules). There are over 390 naturally occurring amino acids known (Voet and Voet 1995, p. 69; Garavelli et al. 2001); however, the protein molecules used by all known living organisms are constructed with the same subset of 22 amino acids.

The univeral genetic code

There must be a mechanism for transmitting information from the genetic material to the catalytic material. All known organisms, with extremely rare exceptions, use the same genetic code for this. The few known exceptions are, nevertheless, simple and minor variations from the "universal" genetic code (see Figure 1.1.1) (Lehman 2001; Voet and Voet 1995, p. 967), exactly as predicted by evolutionary biologists based on the theory of common descent, years before the genetic code was finally solved (Brenner 1957; Crick et al. 1961; Hinegardner and Engelberg 1963; Judson 1996, p. 280-281).

The scientists who cracked the genetic code in the 1950's and 1960's worked under the assumption that the code was universal or nearly so (Judson 1996, p. 280-281). These scientists (which included Francis Crick, Sydney Brenner, George Gamow, and several others) all made this assumption and justified it based upon evolutionary reasoning, even in the complete absence of any direct experimental evidence for a universal code.


Figure 1.1.1. The standard genetic code and known variant nuclear codes. (1) Candida, a unicellular yeast. (2) Micrococcus. (3) ciliated protozoans and green algae. (4) Mycoplasma. (5) suppressor codon in bacteria. (6) Euplotes. (7) the selenocysteine codon (8) Spiroplasma. (9) Micrococcus. (10) resume codon in ssrA RNA (Lehman 2001).

"Crick urged on his companions two other simplifying assumptions of great audacity. ... they assumed, with some apprehension, that the genetic code would be the same for all living things. There was no evidence whatever for this; .... Yet universality seemed inevitable for an obvious reason: since a mutation that changed even one word or letter of the code would alter most of a creature's proteins, it looked sure to be lethal." (Judson 1996, p. 280-281)

In fact, the assumption of a universal genetic code was instrumental in their success in solving the code. For instance, in 1957, nearly ten years before the genetic code was finally solved, Sydney Brenner published an influential paper in which he concluded that all overlapping triplet codes were impossible if the code was universal (Brenner 1957). This paper was widely considered a landmark work, since many researchers were leaning towards an overlapping code. Of course, it turned out that Brenner was correct about the nature of the true code.

In 1961, five years before the code was deciphered, Crick referenced Brenner's work in his landmark report in the journal Nature, "General nature of the genetic code for proteins" (Crick et al. 1961). Although the organism used in the paper was the bacterium E. coli, Crick titled the paper "the genetic code for proteins", not "a genetic code" or "the genetic code of E. coli". In this paper, Crick and others concluded that the code was (1) a triplet code, (2) non-overlapping, and (3) that the code is read from a fixed starting point (i.e. the "start" codon) (Crick et al. 1961). These conclusions were explicitly based on the assumption that the code was essentially the same in tobacco, humans, and bacteria, though there was no direct empirical support for this assumption. These conclusions, when applied to organisms from bacteria to humans, turned out to be correct. Thus, experimental work also assumed a universal code due to common descent.

In fact, in 1963—three years before the code was finally solved—Hinegardner and Engelberg published a paper in Science formally explaining the evolutionary rationale for why the code must be universal (or nearly so) if universal common descent were true, since most mutations in the code would likely be lethal to all living things. Note that, although these early researchers predicted a universal genetic code based on common descent, they also predicted that minor variations could likely be found. Hinegardner and Engelberg allowed for some variation in the genetic code, and predicted how such variation should be distributed if found:

"... if different codes do exist they should be associated with major taxonomic groups such as phyla or kingdoms that have their roots far in the past." (Hinegardner and Engelberg 1963)

Similarly, before alternate codes were found, Francis Crick and Leslie Orgel expressed surprise that minor variants of the code had not been observed yet:

"It is a little surprising that organisms with somewhat different codes do not coexist." (Crick and Orgel 1973, p. 344)

Crick and Orgel were correct in their surprise, and today we know of about a dozen minor variants of the standard, universal genetic code (see the grey, red, and green codons in Figure 1.1.1). As Hinegardner and Engelberg predicted, the minor variations in the standard genetic code are indeed associated with major taxonomic groups (vertebrates vs. plants vs. single-celled ciliates, etc.).

Common metabolism

All known organisms use extremely similar, if not the same, metabolic pathways and metabolic enzymes in processing energy-containing molecules. For example, the fundamental metabolic systems in living organisms are glycolysis, the citric acid cycle, and oxidative phosphorylation. In all eukaryotes and in the majority of prokaryotes, glycolysis is performed in the same ten steps, in the same order, using the same ten enzymes (Voet and Voet 1995, p. 445). In addition, the most basic unit of energy storage, the adenosine triphosphate molecule (ATP), is the same in all species that have been studied.

Potential Falsification:

Thousands of new species are discovered yearly, and new DNA and protein sequences are determined daily from previously unexamined species (Wilson 1992, Ch. 8). At the current rate, which is increasing exponentially, nearly 30,000 new sequences are deposited at GenBank every day, amounting to over 38 million new bases sequenced every day. Each and every one is a test of the theory of common descent. When I first wrote these words in 1999, the rate was less than one tenth what it is today (in 2006), and we now have 20 times the amount of DNA sequenced.

Based solely on the theory of common descent and the genetics of known organisms, we strongly predict that we will never find any modern species from known phyla on this Earth with a foreign, non-nucleic acid genetic material. We also make the strong prediction that all newly discovered species that belong to the known phyla will use the "standard genetic code" or a close derivative thereof. For example, according to the theory, none of the thousands of new and previously unknown insects that are constantly being discovered in the Brazilian rainforest will have non-nucleic acid genomes. Nor will these yet undiscovered species of insects have genetic codes which are not close derivatives of the standard genetic code. In the absence of the theory of common descent, it is quite possible that every species could have a very different genetic code, specific to it only, since there are 1.4 x 1070 informationally equivalent genetic codes, all of which use the same codons and amino acids as the standard genetic code (Yockey 1992). This possibility could be extremely useful for organisms, as it would preclude interspecific viral infections. However, this has not been observed, and the theory of common descent effectively prohibits such an observation.

As another example, nine new lemur and two marmoset species (all primates) were discovered in the forests of Madagascar and Brazil in 2000 (Groves 2000; Rasoloarison et al. 2000; Thalmann and Geissmann 2000). Ten new monkey species have been discovered in Brazil alone since 1990 (Van Roosmalen et al. 2000). Nothing in biology prevents these various species from having a hitherto unknown genetic material or a previously unused genetic code—nothing, that is, except for the theory of common descent. However, we now know definitively that the new lemurs use DNA with the standard genetic code (Yoder et al. 2000); the marmosets have yet to be tested.

Furthermore, each species could use a different polymer for catalysis. The polymers that are used could still be chemically identical yet have different chiralities in different species. There are thousands of thermodynamically equivalent glycolysis pathways (even using the same ten reaction steps but in different orders), so it is possible that every species could have its own specific glycolysis pathway, tailored to its own unique needs. The same reasoning applies to other core metabolic pathways, such as the citric acid cycle and oxidative phosphorylation.

Finally, many molecules besides ATP could serve equally well as the common currency for energy in various species (CTP, TTP, UTP, ITP, or any ATP-like molecule with one of the 293 known amino acids or one of the dozens of other bases replacing the adenosine moiety immediately come to mind). Discovering any new animals or plants that contained any of the anomalous examples proffered above would be potential falsifications of common ancestry, but they have not been found.

Prediction 1.2: A nested hierarchy of species

As seen from the phylogeny in Figure 1, the predicted pattern of organisms at any given point in time can be described as "groups within groups", otherwise known as a nested hierarchy. The only known processes that specifically generate unique, nested, hierarchical patterns are branching evolutionary processes. Common descent is a genetic process in which the state of the present generation/individual is dependent only upon genetic changes that have occurred since the most recent ancestral population/individual. Therefore, gradual evolution from common ancestors must conform to the mathematics of Markov processes and Markov chains. Using Markovian mathematics, it can be rigorously proven that branching Markovian replicating systems produce nested hierarchies (Givnish and Sytsma 1997; Harris 1989; Norris 1997). For these reasons, biologists routinely use branching Markov chains to effectively model evolutionary processes, including complex genetic processes, the temporal distributions of surnames in populations (Galton and Watson 1874), and the behavior of pathogens in epidemics.

The nested hierarchical organization of species contrasts sharply with other possible biological patterns, such as the continuum of "the great chain of being" and the continuums predicted by Lamarck's theory of organic progression (Darwin 1872, pp. 552-553; Futuyma 1998, pp. 88-92). Mere similarity between organisms is not enough to support macroevolution; the nested classification pattern produced by a branching evolutionary process, such as common descent, is much more specific than simple similarity. Real world examples that cannot be objectively classified in nested hierarchies are the elementary particles (which are described by quantum chromodynamics), the elements (whose organization is described by quantum mechanics and illustrated by the periodic table), the planets in our Solar System, books in a library, or specially designed objects like buildings, furniture, cars, etc.

Although it is trivial to classify anything subjectively in a hierarchical manner, only certain things can be classified objectively in a consistent, unique nested hierarchy. The difference drawn here between "subjective" and "objective" is crucial and requires some elaboration, and it is best illustrated by example. Different models of cars certainly could be classified hierarchically—perhaps one could classify cars first by color, then within each color by number of wheels, then within each wheel number by manufacturer, etc. However, another individual may classify the same cars first by manufacturer, then by size, then by year, then by color, etc. The particular classification scheme chosen for the cars is subjective. In contrast, human languages, which have common ancestors and are derived by descent with modification, generally can be classified in objective nested hierarchies (Pei 1949; Ringe 1999). Nobody would reasonably argue that Spanish should be categorized with German instead of with Portugese.

The difference between classifying cars and classifying languages lies in the fact that, with cars, certain characters (for example, color or manufacturer) must be considered more important than other characters in order for the classification to work. Which types of car characters are more important depends upon the personal preference of the individual who is performing the classification. In other words, certain types of characters must be weighted subjectively in order to classify cars in nested hierarchies; cars do not fall into natural, unique, objective nested hierarchies.

Because of these facts, a cladistic analysis of cars will not produce a unique, consistent, well-supported tree that displays nested hierarchies. A cladistic analysis of cars (or, alternatively, a cladistic analysis of imaginary organisms with randomly assigned characters) will of course result in a phylogeny, but there will be a very large number of other phylogenies, many of them with very different topologies, that are as well-supported by the same data. In contrast, a cladistic analysis of organisms or languages will generally result in a well-supported nested hierarchy, without arbitrarily weighting certain characters (Ringe 1999). Cladistic analysis of a true genealogical process produces one or relatively few phylogenetic trees that are much more well-supported by the data than the other possible trees.

Interestingly, Linnaeus, who originally discovered the objective hierarchical classification of living organisms, also tried to classify rocks and minerals hierarchically. However, his classification for non-living objects eventually failed, as it was found to be very subjective. Hierarchical classifications for inanimate objects don't work for the very reason that unlike organisms, rocks and minerals do not evolve by descent with modification from common ancestors.

The degree to which a given phylogeny displays a unique, well-supported, objective nested hierarchy can be rigorously quantified. Several different statistical tests have been developed for determining whether a phylogeny has a subjective or objective nested hierarchy, or whether a given nested hierarchy could have been generated by a chance process instead of a genealogical process (Swofford 1996, p. 504). These tests measure the degree of "cladistic hierarchical structure" (also known as the "phylogenetic signal") in a phylogeny, and phylogenies based upon true genealogical processes give high values of hierarchical structure, whereas subjective phylogenies that have only apparent hierarchical structure (like a phylogeny of cars, for example) give low values (Archie 1989; Faith and Cranston 1991; Farris 1989; Felsenstein 1985; Hillis 1991; Hillis and Huelsenbeck 1992; Huelsenbeck et al. 2001; Klassen et al. 1991).

There is one caveat to consider with this prediction: if rates of evolution are fast, then cladistic information can be lost over time since it would be essentially randomized. The faster the rate, the less time needed to obliterate information about the historical branching pattern of evolution. Slowly evolving characters let us see farther back into time; faster evolving characters restrict that view to more recent events. If the rate of evolution for a certain character is extremely slow, a nested hierarchy will be observed for that character only for very distantly related taxa. However, "rate of evolution" vs. "time since divergence" is relative; if common descent is true, then in some time frame we will always be able to observe a nested hierarchy for any given character. Furthermore, we know empirically that different characters evolve at different rates (e.g. some genes have higher background mutation rates than others). Thus, if common descent is true, we should observe nested hierarchies over a broad range of time at various biological levels.

Therefore, since common descent is a genealogical process, common descent should produce organisms that can be organized into objective nested hierarchies. Equivalently, we predict that, in general, cladistic analyses of organisms should produce phylogenies that have large, statistically significant values of hierarchical structure (in standard scientific practice, a result with "high statistical significance" is a result that has a 1% probability or less of occurring by chance [P < 0.01]). As a representation of universal common descent, the universal tree of life should have very high, very significant hierarchical structure and phylogenetic signal.


Most existing species can be organized rather easily in a nested hierarchical classification. This is evident in the use of the Linnaean classification scheme. Based on shared derived characters, closely related organisms can be placed in one group (such as a genus), several genera can be grouped together into one family, several families can be grouped together into an order, etc.

As a specific example (see Figure 1), plants can be classified as vascular and nonvascular (i.e. they have or lack xylem and phloem). Nested within the vascular group, there are two divisions, seed and non-seed plants. Further nested within the seed plants are two more groups, the angiosperms (which have enclosed, protected seeds) and the gymnosperms (having non-enclosed seeds). Within the angiosperm group are the monocotyledons and the dicotyledons.

Most importantly, the standard phylogenetic tree and nearly all less inclusive evolutionary phylogenies have statistically significant, high values of hierarchical structure (Baldauf et al. 2000; Brown et al. 2001; Hillis 1991; Hillis and Huelsenbeck 1992; Klassen et al. 1991).


Figure 1.2.1. A plot of the CI values of cladograms versus the number of taxa in the cladograms. CI values are on the y-axis; taxa number are on the x-axis. The 95% confidence limits are shown in light turquoise. All points above and to the right of the turquoise region are statistically significant high CI values. Similarly, all points below and to the left of the turquoise region are statistically significant low values of CI. (reproduced from Klassen et al. 1991, Figure 6).

Potential Falsification:

It would be very problematic if many species were found that combined characteristics of different nested groupings. Proceeding with the previous example, some nonvascular plants could have seeds or flowers, like vascular plants, but they do not. Gymnosperms (e.g. conifers or pines) occasionally could be found with flowers, but they never are. Non-seed plants, like ferns, could be found with woody stems; however, only some angiosperms have woody stems. Conceivably, some birds could have mammary glands or hair; some mammals could have feathers (they are an excellent means of insulation). Certain fish or amphibians could have differentiated or cusped teeth, but these are only characteristics of mammals. A mix and match of characters like this would make it extremely difficult to objectively organize species into nested hierarchies. Unlike organisms, cars do have a mix and match of characters, and this is precisely why a nested hierarchy does not flow naturally from classification of cars.

If it were impossible, or very problematic, to place species in an objective nested classification scheme (as it is for the car, chair, book, atomic element, and elementary particle examples mentioned above), macroevolution would be effectively disproven. More precisely, if the phylogenetic tree of all life gave statistically significant low values of phylogenetic signal (hierarchical structure), common descent would be resolutely falsified.

In fact, it is possible to have a "reciprocal" pattern from nested hierarchies. Mathematically, a nested hierarchy is the result of specific correlations between certain characters of organisms. When evolutionary rates are fast, characters become randomly distributed with respect to one another, and the correlations are weakened. However, the characters can also be anti-correlated—it is possible for them to be correlated in the opposite direction from what produces nested hierarchies (Archie 1989; Faith and Cranston 1991; Hillis 1991; Hillis and Huelsenbeck 1992; Klassen et al. 1991). The observation of such an anti-correlated pattern would be a strong falsification of common descent, regardless of evolutionary rates.

One widely used measure of cladistic hierarchical structure is the consistency index (CI). The statistical properties of the CI measure were investigated in a frequently cited paper by Klassen et al. (Klassen et al. 1991; see Figure 1.2.1). The exact CI value is dependent upon the number of taxa in the phylogenetic tree under consideration. In this paper, the authors calculated what values of CI were statistically significant for various numbers of taxa. Higher values of CI indicate a greater degree of hierarchical structure.

As an example, a CI of 0.2 is expected from random data for 20 taxa. A value of 0.3 is, however, highly statistically significant. Most interesting for the present point is the fact that a CI of 0.1 for 20 taxa is also highly statistically significant, but it is too low—it is indicative of anti-cladistic structure. Klassen et al. took 75 CI values from published cladograms in 1989 (combined from three papers) and noted how they fared in terms of statistical significance. The cladograms used from 5 to 49 different taxa (i.e. different species). Three of the 75 cladograms fell within the 95% confidence limits for random data, which means that they were indistinguishable from random data. All the rest exhibited highly statistically significant values of CI. None exhibited significant low values; none displayed an anti-correlated, anti-hierarchical pattern.

Note, this study was performed before there were measures of statistical significance which would allow researchers to "weed out" the bad cladograms. Predictably, the three "bad" data sets considered under ten taxa—it is of course more difficult to determine statistical significance with very little data. Seventy-five independent studies from different researchers, on different organisms and genes, with high values of CI (P < 0.01) is an incredible confirmation with an astronomical degree of combined statistical significance (P << 10-300, Bailey and Gribskov 1998; Fisher 1990). If the reverse were true—if studies such as this gave statistically significant values of CI (i.e. cladistic hierarchical structure) which were lower than that expected from random data—common descent would have been firmly falsified.

Keep in mind that about 1.5 million species are known currently, and that the majority of these species has been discovered since Darwin first stated his hypothesis of common ancestry. Even so, they all have fit the correct hierarchical pattern within the error of our methods. Furthermore, it is estimated that only 1 to 10% of all living species has even been catalogued, let alone studied in detail. New species discoveries pour in daily, and each one is a test of the theory of common descent (Wilson 1992, Ch. 8).

Prediction 1.3: Consilience of independent phylogenies

"It will be determined to what extent the phylogenetic tree, as derived from molecular data in complete independence from the results of organismal biology, coincides with the phylogenetic tree constructed on the basis of organismal biology. If the two phylogenetic trees are mostly in agreement with respect to the topology of branching, the best available single proof of the reality of macro-evolution would be furnished. Indeed, only the theory of evolution, combined with the realization that events at any supramolecular level are consistent with molecular events, could reasonably account for such a congruence between lines of evidence obtained independently, namely amino acid sequences of homologous polypeptide chains on the one hand, and the finds of organismal taxonomy and paleontology on the other hand. Besides offering an intellectual satisfaction to some, the advertising of such evidence would of course amount to beating a dead horse. Some beating of dead horses may be ethical, when here and there they display unexpected twitches that look like life."

Emile Zuckerkandl and Linus Pauling, discussing the possibility of the twin nested hierarchy before the first molecular phylogenies had been made.
(1965) "Evolutionary Divergence and Convergence in Proteins." in Evolving Genes and Proteins, p. 101.

Here we commence to beat Pauling's poor 40-year dead horse. If there is one historical phylogenetic tree which unites all species in an objective genealogy, all separate lines of evidence should converge on the same tree (Penny et al. 1982; Penny et al. 1991; Zuckerkandl and Pauling 1965). Independently derived phylogenetic trees of all organisms should match each other with a high degree of statistical significance.


Well-determined phylogenetic trees inferred from the independent evidence of morphology and molecular sequences match with an extremely high degree of statistical significance. Many genes with very basic cellular functions are ubiquitous—they occur in the genomes of most or all organisms. An oft-cited example is the cytochrome c gene. Since all eukaryotes contain the gene for this essential protein, neither its presence nor its function correlates with organismal morphology. Additionally, because of the fact of DNA coding redundancy, parts of certain DNA sequences have absolutely no correlation with phenotype (e.g. certain introns or the four-fold degenerate third-base position of most DNA codons). Due to these two aspects of certain DNA sequences, ubiquity and redundancy, DNA sequences can be carefully chosen that constitute completely independent data from morphology. (See point 17 and 18 for more background about the molecular sequence evidence and for more detail about how it is independent of morphology.) The degree of phylogenetic congruence between these independent data sets is nothing short of incredible.

In science, independent measurements of theoretical values are never exact. When inferring any value (such as a physical constant like the charge of the electron, the mass of the proton, or the speed of light) some error always exists in the measurement, and all independent measurements are incongruent to some extent. Of course, the true value of something is never known for certain in science—all we have are measurements that we hope approximate the true value. Scientifically, then, the important relevant questions are "When comparing two measurements, how much of a discrepancy does it take to be a problem?" and "How close must the measurements be in order to give a strong confirmation?" Scientists answer these questions quantitatively with probability and statistics (Box 1978; Fisher 1990; Wadsworth 1997). To be scientifically rigorous we require statistical significance. Some measurements of a given value match with statistical significance (good), and some do not (bad), even though no measurements match exactly (reality).

Table 1.3.1: Number of possible phylogenetic trees for a given number of taxa
Number of taxa Number of possible trees
NR = (2n-3)!! = (2n-3)!/(2n-2(n-2)!)
rooted unrooted
2 3 1
3 4 3
4 5 15
5 6 105
6 7 945
7 8 10,395
8 9 135,135
9 10 2,027,025
10 11 34,459,425
11 12 654,729,075
12 13 13,749,310,575
13 14 316,234,143,225
14 15 7,905,853,580,625
15 16 213,458,046,676,875
16 17 6,190,283,353,629,375
17 18 191,898,783,962,510,625
18 19 6,332,659,870,762,850,625
19 20 221,643,095,476,699,771,875
20 21 8,200,794,532,637,891,559,375

Phylogenetic Trees Calculator

rooted   unrooted
Number of Taxa:

Number of different trees =
1.2e+50 means 1.2 × 1050

So, how well do phylogenetic trees from morphological studies match the trees made from independent molecular studies? There are over 1038 different possible ways to arrange the 30 major taxa represented in Figure 1 into a phylogenetic tree (see Table 1.3.1; Felsenstein 1982; Li 1997, p. 102). In spite of these odds, the relationships given in Figure 1, as determined from morphological characters, are completely congruent with the relationships determined independently from cytochrome c molecular studies (for consensus phylogenies from pre-molecular studies see Carter 1954, Figure 1, p. 13; Dodson 1960, Figures 43, p. 125, and Figure 50, p. 150; Osborn 1918, Figure 42, p. 161; Haeckel 1898, p. 55; Gregory 1951, Fig. opposite title page; for phylogenies from the early cytochrome c studies see McLaughlin and Dayhoff 1973; Dickerson and Timkovich 1975, pp. 438-439). Speaking quantitatively, independent morphological and molecular measurements such as these have determined the standard phylogenetic tree, as shown in Figure 1, to better than 38 decimal places. This phenomenal corroboration of universal common descent is referred to as the "twin nested hierarchy". This term is something of a misnomer, however, since there are in reality multiple nested hierarchies, independently determined from many sources of data.

When two independently determined trees mismatch by some branches, they are called "incongruent". In general, phylogenetic trees may be very incongruent and still match with an extremely high degree of statistical significance (Hendy et al. 1984; Penny et al. 1982; Penny and Hendy 1986; Steel and Penny 1993). Even for a phylogeny with a small number of organisms, the total number of possible trees is extremely large. For example, there are about a thousand different possible phylogenies for only six organisms; for nine organisms, there are millions of possible phylogenies; for 12 organisms, there are nearly 14 trillion different possible phylogenies (Table 1.3.1; Felsenstein 1982; Li 1997, p. 102). Thus, the probability of finding two similar trees by chance via two independent methods is extremely small in most cases. In fact, two different trees of 16 organisms that mismatch by as many as 10 branches still match with high statistical significance (Hendy et al. 1984, Table 4; Steel and Penny 1993). For more information on the statistical significance of trees that do not match exactly, see "Statistics of Incongruent Phylogenetic Trees".

The stunning degree of match between even the most incongruent phylogenetic trees found in the biological literature is widely unappreciated, mainly because most people (including many biologists) are unaware of the mathematics involved (Bryant et al. 2002; Penny et al. 1982; Penny and Hendy 1986). Penny and Hendy have performed a series of detailed statistical analyses of the significance of incongruent phylogenetic trees, and here is their conclusion:

"Biologists seem to seek the 'The One Tree' and appear not to be satisfied by a range of options. However, there is no logical difficulty in having a range of trees. There are 34,459,425 possible [unrooted] trees for 11 taxa (Penny et al. 1982), and to reduce this to the order of 10-50 trees is analogous to an accuracy of measurement of approximately one part in 106." (Penny and Hendy 1986, p. 414)

For a more realistic universal phylogenetic tree with dozens of taxa including all known phyla, the accuracy is better by many orders of magnitude. To put the significance of this incredible confirmation in perspective, consider the modern theory of gravity. Both Newton's Theory of Universal Gravitation and Einstein's General Theory of Relativity rely upon a fundamental physical constant, G, the gravitational constant. If these theories of gravity are correct, independent methods should determine similar values for G. However, to date, very precise independent measurements of the gravitational constant G disagree by nearly 1% (Kestenbaum 1998; Quinn 2000). Here is how David Kestenbaum describes the current scientific status of the theory of gravity, as reported in the prestigious journal Science:

"While the charge of the electron is known to seven decimal places, physicists lose track of G after only the third. For some, that's an embarrassment. 'It grates on me like a burr in the saddle,' says Alvin Sanders, a physicist at the University of Virginia in Charlottesville. Over the past few decades, he and a handful of other physicists have dedicated themselves to measuring G more accurately. To their dismay, they've come up with wildly different values. 'You might say we've had negative progress,' says Barry Taylor, a physicist at the National Institute of Standards and Technology (NIST) in Gaithersburg, Maryland. ...

'Nobody understands it [the far-out results of the PTB, the German standards lab in Braunschweig],' says Meyer. 'They must have made an unbelievable mistake, but we cannot find it.' ... says Terry Quinn, 'we may just have to throw the PTB result out.'" (Kestenbaum 1998)

Over two years later, the same Terry Quinn (of the International Bureau of Weights and Measures [BIPM] in Sèvres, France) summarized the situation in a review for the journal Nature:

"The current interest in measuring G was stimulated by the publication in 1996 of a value for G that differed by 0.6% from the accepted value given in the previous 1986 CODATA report. To take account of this, the 1998 CODATA report recommends a value for G ... with an uncertainty of 0.15%, some ten times worse than in 1986. Whereas the other fundamental constants were more accurately known in 1998 than in 1986, the uncertainty in G increased dramatically. The G community appeared to be going backwards rather than forwards." (Quinn 2000)

Nevertheless, a precision of just under 1% is still pretty good; it is not enough, at this point, to cause us to cast much doubt upon the validity and usefulness of modern theories of gravity. However, if tests of the theory of common descent performed that poorly, different phylogenetic trees, as shown in Figure 1, would have to differ by 18 of the 30 branches! In their quest for scientific perfection, some biologists are rightly rankled at the obvious discrepancies between some phylogenetic trees (Gura 2000; Patterson et al. 1993; Maley and Marshall 1998). However, as illustrated in Figure 1, the standard phylogenetic tree is known to 38 decimal places, which is a much greater precision than that of even the most well-determined physical constants. For comparison, the charge of the electron is known to only seven decimal places, the Planck constant is known to only eight decimal places, the mass of the neutron, proton, and electron are all known to only nine decimal places, and the universal gravitational constant has been determined to only three decimal places.

Furthermore, if common descent is true, we expect that including more data in phylogenetic analyses will increase the correspondence between phylogenetic trees. As explained in the phylogenetic caveats sidebar, gene trees are not equivalent to species trees (Avise and Wollenberg 1997; Fitch 1970; Hudson 1992; Nichols 2001; Wu 1991). Genetics and heredity are stochastic (i.e. probabilistic) processes, and consequently we expect that phylogenies constructed with single genes will be partially incongruent. However, including multiple independent genes in a phylogenetic analysis should circumvent this difficulty; in general more than five independent genes are needed to accurately reconstruct a species phylogeny (Wu 1991). Phylogenetic trees constructed with multiple genes should thus be more accurate than those constructed with single genes, and indeed combined gene trees are more congruent (Baldauf et al. 2000; Hedges 1994; Hedges and Poling 1999; Penny et al. 1982).

Potential Falsification:

When it became possible to sequence biological molecules, the realization of a markedly different tree based on the independent molecular evidence would have been a fatal blow to the theory of evolution, even though that is by far the most likely result. More precisely, the common descent hypothesis would have been falsified if the universal phylogenetic trees determined from the independent molecular and morphological evidence did not match with statistical significance. Furthermore, we are now in a position to begin construction of phylogenetic trees based on other independent lines of data, such as chromosomal organization. In a very general sense, chromosome number and length and the chromosomal position of genes are all causally independent of both morphology and of sequence identity. Phylogenies constructed from these data should recapitulate the standard phylogenetic tree as well (Hillis et al. 1996; Li 1997).


One common objection is the assertion that anatomy is not independent of biochemistry, and thus anatomically similar organisms are likely to be similar biochemically (e.g. in their molecular sequences) simply for functional reasons. According to this argument, then, we should expect phylogenies based on molecular sequences to be similar to phylogenies based on morphology even if organisms are not related by common descent. This argument is very wrong. There is no known biological reason, besides common descent, to suppose that similar morphologies must have similar biochemistry. Though this logic may seem quite reasonable initially, all of molecular biology refutes this "common sense" correlation. In general, similar DNA and biochemistry give similar morphology and function, but the converse is not true—similar morphology and function is not necessarily the result of similar DNA or biochemistry. The reason is easily understood once explained; many very different DNA sequences or biochemical structures can result in the same functions and the same morphologies (see predictions 4.1 and 4.2 for a detailed explanation).

As a close analogy, consider computer programs. Netscape works essentially the same on a Macintosh, an IBM, or a Unix machine, but the binary code for each program is quite different. Computer programs that perform the same functions can be written in most any computer language—Basic, Fortran, C, C++, Java, Pascal, etc. and identical programs can be compiled into binary code many different ways. Furthermore, even using the same computer language, there are many different ways to write any specific computer program, even using the same algorithms and subroutines. In the end, there is no reason to suspect that similar computer programs are written with similar code, based solely on the function of the program. This is the reason why software companies keep their source code secret, but they don't care that competitors can use their programs—it is essentially impossible to deduce the program code from the function and operation of the software. The same conclusion applies to biological organisms, for very similar reasons.

To reiterate, although similar genotypes (e.g. molecular sequences) often give similar phenotypes (e.g. morphological characters), similar phenotypes are not necessarily the result of similar genotypes. Thus, it is entirely possible that phylogenetic trees constructed from genotypic data could be radically different from phylogenetic trees constructed from phenotypic data. In fact, in the absence of common descent or any other reason to suppose that these two types of trees should be similar, the most likely result by far is that they will be radically different. This is precisely why it is possible to falsify the macroevolutionary prediction that independently derived phylogenies should be similar.

Prediction 1.4: Intermediate and transitional forms: the possible morphologies of predicted common ancestors

All fossilized animals found should conform to the standard phylogenetic tree. If all organisms are united by descent from a common ancestor, then there is one single true historical phylogeny for all organisms. Similarly, there is one single true historical genealogy for any individual human. It directly follows that if there is a unique universal phylogeny, then all organisms, both past and present, fit in that phylogeny uniquely. Since the standard phylogenetic tree is the best approximation of the true historical phylogeny, we expect that all fossilized animals should conform to the standard phylogenetic tree within the error of our scientific methods.

Every node shared between two branches in a phylogeny or cladogram represents a predicted common ancestor; thus there are ~29 common ancestors predicted from the tree shown in Figure 1. Our standard tree shows that the bird grouping is most closely related to the reptilian grouping, with a node linking the two (A in Figure 1); thus we predict the possibility of finding fossil intermediates between birds and reptiles. The same reasoning applies to mammals and reptiles (B in Figure 1). However, we predict that we should never find fossil intermediates between birds and mammals.

It should be pointed out that there is no requirement for intermediate organisms to go extinct. In fact, all living organisms can be thought of as intermediate between adjacent taxa in a phylogenetic tree. For instance, modern reptiles are intermediate between amphibians and mammals, and reptiles are also intermediate between amphibians and birds. As far as macroevolutionary predictions of morphology are concerned, this point is trivial, as it is essentially just a restatement of the concept of a nested hierarchy.

However, a phylogenetic tree does make significant predictions about the morphology of intermediates which no longer exist or which have yet to be discovered. Each predicted common ancestor has a set of explicitly specified morphological characteristics, based on each of the most common derived characters of its descendants and based upon the transitions that must have occurred to transform one taxa into another (Cunningham et al. 1998; Futuyma 1998, pp. 107-108). From the knowledge of avian and reptilian morphology, it is possible to predict some of the characteristics that a reptile-bird intermediate should have, if found. Therefore, we expect the possibility of finding reptile-like fossils with feathers, bird-like fossils with teeth, or bird-like fossils with long reptilian tails. However, we do not expect transitional fossils between birds and mammals, like mammalian fossils with feathers or bird-like fossils with mammalian-style middle ear bones.


Example 1: bird-reptiles

In the case just mentioned, we have found a quite complete set of dinosaur-to-bird transitional fossils with no morphological "gaps" (Sereno 1999), represented by Eoraptor, Herrerasaurus, Ceratosaurus, Allosaurus, Compsognathus, Sinosauropteryx, Protarchaeopteryx, Caudipteryx, Velociraptor, Sinovenator, Beipiaosaurus, Sinornithosaurus, Microraptor, Archaeopteryx, Rahonavis, Confuciusornis, Sinornis, Patagopteryx, Hesperornis, Apsaravis, Ichthyornis, and Columba, among many others (Carroll 1997, pp. 306-323; Norell and Clarke 2001; Sereno 1999; Xu et al. 1999; Xu et al. 2000; Xu et al. 2002). All have the expected possible morphologies (see Figure 3.1.1 from Prediction 3.1 for a few examples), including organisms such as Protarchaeopteryx, Caudipteryx, and the famous "BPM 1 3-13" (a dromaeosaur from China now named Cryptovolans pauli; Czerkas et al. 2002 ) which are flightless bipedal dinosaurs with modern-style feathers (Chen et al. 1998 ; Qiang et al. 1998; Norell et al. 2002). Additionally, several similar flightless dinosaurs have been found covered with nascent evolutionary precursors to modern feathers (branched feather-like integument indistinguishable from the contour feathers of true birds), including Sinornithosaurus ("Bambiraptor"), Sinosauropteryx, Beipiaosaurus, Microraptor, and an unnamed dromaeosaur specimen, NGMC 91, informally called "Dave" (Ji et al. 2001). The All About Archaeopteryx FAQ gives a detailed listing of the various characters of Archaeopteryx which are intermediate between reptiles and modern birds.

Example 2: reptile-mammals

[Figure1.4.1 (cartoon of vertebrate jaws)]

Figure 1.4.1. The jaws of three vertebrates—mammal, therapsid, and pelycosaur. A side view of three idealized skulls of mammals, therapsids (mammal-like reptiles), and pelycosaurs (early reptiles). The figure shows the differences between mammal and reptilian jaws and ear-bone structures. The jaw joint is shown as a large black dot, the quadrate (mammalian anvil or incus) is in turquoise, the articular (mammalian hammer or malleus) is in yellow, and the angular (mammalian tympanic annulus) is in pink. Note how, in the reptile, the jaw joint is formed between the blue quadrate and the yellow articular (with the pink angular close by), and how, in the mammal, the jaw joint is formed between the squamosal above and the dentary below. In the reptile, the squamosal is just above and contacting the quadrate. Advanced therapsids have two jaw joints: a reptile-like joint and a mammal-like joint (Figure based on Kardong 2002, pp. 275, reproduced with permission from the publisher, Copyright © 2002 McGraw-Hill)

We also have an exquisitely complete series of fossils for the reptile-mammal intermediates, ranging from the pelycosauria, therapsida, cynodonta, up to primitive mammalia (Carroll 1988, pp. 392-396; Futuyma 1998, pp. 146-151; Gould 1990; Kardong 2002, pp. 255-275). As mentioned above, the standard phylogenetic tree indicates that mammals gradually evolved from a reptile-like ancestor, and that transitional species must have existed which were morphologically intermediate between reptiles and mammals—even though none are found living today. However, there are significant morphological differences between modern reptiles and modern mammals. Bones, of course, are what fossilize most readily, and that is where we look for transitional species from the past. Osteologically, two major striking differences exist between reptiles and mammals: (1) reptiles have at least four bones in the lower jaw (e.g. the dentary, articular, angular, surangular, and coronoid), while mammals have only one (the dentary), and (2) reptiles have only one middle ear bone (the stapes), while mammals have three (the hammer, anvil, and stapes) (see Figure 1.4.1).

Early in the 20th century, developmental biologists discovered something that further complicates the picture. In the reptilian fetus, two developing bones from the head eventually form two bones in the reptilian lower jaw, the quadrate and the articular (see the Pelycosaur in Figure 1.4.1). Surprisingly, the corresponding developing bones in the mammalian fetus eventually form the anvil and hammer of the unique mammalian middle ear (also known more formally as the incus and malleus, respectively; see Figure 1.4.2) (Gilbert 1997, pp. 894-896). These facts strongly indicated that the hammer and anvil had evolved from these reptilian jawbones—that is, if common descent was in fact true. This result was so striking, and the required intermediates so outlandish, that many anatomists had extreme trouble imagining how transitional forms bridging these morphologies could have existed while retaining function. Young-earth creationist Duane Gish stated the problem this way:

"All mammals, living or fossil, have a single bone, the dentary, on each side of the lower jaw, and all mammals, living or fossil, have three auditory ossicles or ear bones, the malleus, incus and stapes. ... Every reptile, living or fossil, however, has at least four bones in the lower jaw and only one auditory ossicle, the stapes. ... There are no transitional fossil forms showing, for instance, three or two jawbones, or two ear bones. No one has explained yet, for that matter, how the transitional form would have managed to chew while his jaw was being unhinged and rearticulated, or how he would hear while dragging two of his jaw bones up into his ear." (Gish 1978, p. 80)

[Figure1.4.2a (cartoon of vertebrate ears)] [Figure1.4.2b (cartoon of vertebrate ears)]

Figure 1.4.2. A comparison of the ears of reptiles and mammals. The reptile ear is shown on the left, the mammal ear on the right. As in Figure 1.4.1, the quadrate (mammalian anvil or incus) is in turquoise and the articular (mammalian hammer or malleus) is in yellow. The stapes is shown in brown. Note how the relative arrangement of these bones is similar in both taxa, in the order of inner ear-stapes-quadrate-articular.

Gish was incorrect in stating that there were no transitional fossil forms, and he has been corrected on this gaffe numerous times since he wrote these words. However, Gish's statements nicely delineate the morphological conundrum at hand. Let's review the required evolutionary conclusion. During their evolution, two mammalian middle ear bones (the hammer and anvil, aka malleus and incus) were derived from two reptilian jawbones. Thus there was a major evolutionary transition in which several reptilian jawbones (the quadrate, articular, and angular) were extensively reduced and modified gradually to form the modern mammalian middle ear. At the same time, the dentary bone, a part of the reptilian jaw, was expanded to form the major mammalian lower jawbone. During the course of this change, the bones that form the hinge joint of the jaw changed identity. Importantly, the reptilian jaw joint is formed at the intersection of the quadrate and articular whereas the mammalian jaw joint is formed at the intersection of the squamosal and dentary (see Figure 1.4.1).

How could hearing and jaw articulation be preserved during this transition? As clearly shown from the many transitional fossils that have been found (see Figure 1.4.3), the bones that transfer sound in the reptilian and mammalian ear were in contact with each other throughout the evolution of this transition. In reptiles, the stapes contacts the quadrate, which in turn contacts the articular. In mammals, the stapes contacts the incus, which in turn contacts the malleus (see Figure 1.4.2). Since the quadrate evolved into the incus, and the articular evolved into the malleus, these three bones were in constant contact during this impressive evolutionary change. Furthermore, a functional jaw joint was maintained by redundancy—several of the intermediate fossils have both a reptilian jaw joint (from the quadrate and articular) and a mammalian jaw joint (from the dentary and squamosal). Several late cynodonts and Morganucodon clearly have a double-jointed jaw. In this way, the reptilian-style jaw joint was freed to evolve a new specialized function in the middle ear. It is worthy of note that some modern species of snakes have a double-jointed jaw involving different bones, so such a mechanical arrangement is certainly possible and functional.

Since Figure 1.4.3 was made, several important intermediate fossils have been discovered that fit between Morganucodon and the earliest mammals. These new discoveries include a complete skull of Hadrocodium wui (Luo et al. 2001) and cranial and jaw material from Repenomamus and Gobiconodon (Wang et al. 2001). These new fossil finds clarify exactly when and how the malleus, incus, and angular completely detached from the lower jaw and became solely auditory ear ossicles.

Recall that Gish stated: "There are no transitional fossil forms showing, for instance, three or two jawbones, or two ear bones" (Gish 1978, p. 80). Gish simply does not understand how gradual transitions happen (something he should understand, obviously, if he intends to criticize evolutionary theory). These fossil intermediates illustrate why Gish's statement is a gross mischaracterization of how a transitional form should look. In several of the known intermediates, the bones have overlapping functions, and one bone can be called both an ear bone and a jaw bone; these bones serve two functions. Thus, there is no reason to expect transitional forms with intermediate numbers of jaw bones or ear bones. For example, in Morganucodon, the quadrate (anvil) and the articular (hammer) serve as mammalian-style ear bones and reptilian jaw bones simultaneously. In fact, even in modern reptiles the quadrate and articular serve to transmit sound to the stapes and the inner ear (see Figure 1.4.2). The relevant transition, then, is a process where the ear bones, initially located in the lower jaw, become specialized in function by eventually detaching from the lower jaw and moving closer to the inner ear.

[Figure1.4.3 (cartoon of vertebrate jaws)]

Figure 1.4.3. A comparison of the jawbones and ear-bones of several transitional forms in the evolution of mammals. Approximate stratigraphic ranges of the various taxa are indicated at the far left (more recent on top). The left column of jawbones shows the view of the left jawbone from the inside of the mouth. The right column is the view of the right jawbone from the right side (outside of the skull). As in Figure 1.4.1, the quadrate (mammalian anvil or incus) is in turquoise, the articular (mammalian hammer or malleus) is in yellow, and the angular (mammalian tympanic annulus) is in pink. For clarity, the teeth are not shown, and the squamosal upper jawbone is omitted (it replaces the quadrate in the mammalian jaw joint, and forms part of the jaw joint in advanced cynodonts and Morganucodon). Q = quadrate, Ar = articular, An = angular, I = incus (anvil), Ma = malleus (hammer), Ty = tympanic annulus, D = dentary. (Reproduced from Kardong 2002, pp. 274, with permission from the publisher, Copyright © 2002 McGraw-Hill)

Example 3: human-apes

"All advances by degrees in Nature, and nothing by leaps, and this law as applied to each, is part of my doctrine of Continuity. Although there may exist in some other world species intermediate between Man and the Apes, Nature has thought it best to remove them from us, in order to establish our superiority beyond question. I speak of intermediate species, and by no means limit myself to those leading to Man."

– Gottfried Wilhelm Leibniz
Protogaea 1749

One of the most celebrated examples of transitional fossils is our collection of fossil hominids (see Figure 1.4.4 below). Based upon the consensus of numerous phylogenetic analyses, Pan troglodytes (the chimpanzee) is the closest living relative of humans. Thus, we expect that organisms lived in the past which were intermediate in morphology between humans and chimpanzees. Over the past century, many spectacular paleontological finds have identified such transitional hominid fossils.

[Figure 1.4.4: Hominid skulls]

Figure 1.4.4. Fossil hominid skulls. (Images © 2000 Smithsonian Institution.) (larger 76K JPG version)

Example 4: legged fossil whales

"Finally, and most glaringly obvious, if random evolution is true there must have been a large number of transitional forms between the mesonychid and the ancient whale: Where are they? It seems like quite a coincidence that of all the intermediate species that must have existed between the mesonychid and whale, only species that are very similar to the end species have been found. (Behe 1994)"

– Michael J. Behe

Anti-Darwinian, Intelligent Design conjecturist,
writing against the validity of evolution less than a year before three transitional species between whales and land-dwelling Eocene Mesonychids were found.

Another impressive example of incontrovertible transitional forms predicted to exist by evolutionary biologists is the collection of land mammal-to-whale fossil intermediates. Whales, of course, are sea animals with flippers, lacking external hindlimbs. Since they are also mammals, the consensus phylogeny indicates that whales and dolphins evolved from land mammals with legs. In recent years, we have found several transitional forms of whales with legs, both capable and incapable of terrestrial locomotion (Gingerich et al. 1983; Gingerich et al. 1990; Gingerich et al. 1994; Gingerich et al. 2001; Thewissen et al. 2001).

Example 5: legged seacows

[Figure 1.4.5: Seacow with legs]

Figure 1.4.5. An intermediate sirenian fossil: a seacow with legs. Reconstructed skeleton of Pezosiren portelli. Length is approximately 7 feet. Gray coloring represents extant fossils; white elements are partly conjectural. (Domning 2001; Image © 2001 Macmillan Magazines Ltd.)

Seacows (manatees and dugongs) are fully aquatic mammals with flippers for forelimbs and no hindlimbs. Evolutionary theory predicts that seacows evolved from terrestrial ancestors with legs, and that thus we could find seacow intermediates with legs. Recently, a new transitional fossil has been found in Jamaica, a seacow with four legs (Domning 2001).

There are many other examples such as these—most can be found in the excellent Transitional Vertebrate Fossils FAQ.

Potential Falsification:

Any finding of a striking half-mammal, half-bird intermediate would be highly inconsistent with common descent. Many other examples of prohibited intermediates can be thought of, based on the standard tree (Kemp 1982; Stanley 1993; Carroll 1997; Chaterjee 1997).

A subtle, yet important point is that a strict cladistic evolutionary interpretation precludes the possibility of identifying true ancestors; only intermediates or transitionals can be positively identified. (For the purposes of this article, transitionals and intermediates are considered synonymous.) The only incontrovertible evidence for an ancestor-descendant relationship is the observation of a birth; obviously this is normally rather improbable in the fossil record. Intermediates are not necessarily the same as the exact predicted ancestors; in fact, it is rather unlikely that they would be the same. Simply due to probability considerations, the intermediates that we find will most likely not be the true ancestors of any modern species, but will be closely related to a predicted common ancestor. Therefore, the intermediates we do find will likely have additional derived characters besides the characters that identified them as intermediates. Because of these considerations, when a new and important intermediate fossil species is discovered, careful paleontologists will often note that the transitional species under study is probably not an ancestor, but rather is "representative of a common ancestor" or is an evolutionary "side-branch". The fewer extra dervied characters that an intermediate fossil has, the higher the probability that an intermediate fossil is an actual ancestor. For further clarification see prediction 5.4.

Prediction 1.5: Chronological order of intermediates

Fossilized intermediates should appear in the correct general chronological order based on the standard tree. Any phylogenetic tree predicts a relative chronological order of the evolution of hypothetical common ancestors and intermediates between these ancestors. For instance, in our current example, the reptile-mammal common ancestor (B) and intermediates should be older than the reptile-bird common ancestor (A) and intermediates.

Note, however, that there is some "play" within the temporal constraints demanded by any phylogeny, for two primary reasons: (1) the statistical confidence (or conversely, the error) associated with a phylogeny and its specific internal branches, and (2) the inherent resolution of the fossil record (ultimately stemming from the vagaries of the fossilization process). As mentioned earlier, most phylogenetic trees have some branches with high confidence, because they are well-supported by the data, and other branches in which we have less confidence, because they are statistically less significant and poorly supported by the data. See also the caveats associated with phylogenetic analysis.

When evaluating the geological order of fossils, remember that once a transitional species appears there is no reason why it must become extinct and be replaced. For instance, some organisms have undergone little change in as much as 100 to 200 million years in rare cases. Some familiar examples are the "living fossils", such as the coelacanth, which has persisted for approximately 80 million years; the bat, which has not changed much in the past 50 million years; and even the modern tree squirrel, which has not changed in 35 million years. In fact, paleontological studies indicate the average longevity of 21 living families of vertebrates is approximately 70 million years (Carroll 1997, p. 167).

Furthermore, the fossil record is demonstrably incomplete; species appear in the fossil record, then disappear, then reappear later. An exceptional instance is the coelacanth, which last appeared in the fossil record 80 million years ago, yet it is alive today. During the Cretaceous (a critical time in bird evolution), there is a 50 million-year gap in the diplodocoidean record, greater than a 40 million-year gap in the pachycephalosaurian record, greater than a 20 million-year gap in the trodontidiae, and about a 15 million-year gap in the oviraptosaurian fossil record (both of these last two orders of dinosaurs are maniraptoran coelurosaurian theropods, which figure significantly in the evolution of birds). During the Jurassic, there is a 40 million-year gap in the fossil record of the heterodontosauridae (Sereno 1999). Most organisms do not fossilize, and there is no reason why a representative of some species must be found in the fossil record. As every graduate student in scientific research knows (or eventually learns, perhaps the hard way), arguments based upon negative evidence are very weak scientific arguments, especially in the absence of proper positive controls. Thus, based on the fossil remains of modern species and the known gaps in the current paleontological records of extinct species, the observation of transitional species "out of order" by 40 million years should be fairly common. This degree of "play" in the fossil record is actually rather minor, considering that the fossil record of life spans between 2 to 3.8 billion years and that of multicellular organisms encompasses a total of ~660 million years. An uncertainty of 40 million years is equivalent to about a 1% or 6% relative error, respectively—rather small overall.


The reptile-bird intermediates mentioned above date from the Upper Jurassic and Lower Cretaceous (about 150 million years ago), whereas pelycosauria and therapsida (reptile-mammal intermediates) are older and date from the Carboniferous and the Permian (about 250 to 350 million years ago, see the Geological Time Scale). This is precisely what should be observed if the fossil record matches the standard phylogenetic tree.

The most scientifically rigorous method of confirming this prediction is to demonstrate a positive corellation between phylogeny and stratigraphy, i.e. a positive corellation between the order of taxa in a phylogenetic tree and the geological order in which those taxa first appear and last appear (whether for living or extinct intermediates). For instance, within the error inherent in the fossil record, prokaryotes should appear first, followed by simple multicellular animals like sponges and starfish, then lampreys, fish, amphibians, reptiles, mammals, etc., as shown in Figure 1. Contrary to the erroneous (and unreferenced) opinions of some anti-evolutionists (e.g. Wise 1994, p. 225-226), studies from the past ten years addressing this very issue have confirmed that there is indeed a positive corellation between phylogeny and stratigraphy, with statistical significance (Benton 1998; Benton and Hitchin 1996; Benton and Hitchin 1997; Benton et al. 1999; Benton et al. 2000; Benton and Storrs 1994; Clyde and Fisher 1997; Hitchin and Benton 1997; Huelsenbeck 1994; Norell and Novacek 1992a; Norell and Novacek 1992b; Wills 1999). Using three different measures of phylogeny-stratigraphy correlation [the RCI, GER, and SCI (Ghosts 2.4 software, Wills 1999)], a high positive correlation was found between the standard phylogenetic tree portrayed in Figure 1 and the stratigraphic range of the same taxa, with very high statistical significance (P < 0.0001) (this work, Ghosts input file available upon request).

As another specific example, an early analysis published in Science by Mark Norell and Michael Novacek (Norell and Novacek 1992b) examined 24 different taxa of vertebrates (teleosts, amniotes, reptiles, synapsids, diapsids, lepidosaurs, squamates, two orders of dinosaurs, two orders of hadrosaurs, pachycephalosaurs, higher mammals, primates, rodents, ungulates, artiodactyls, ruminants, elephantiformes, brontotheres, tapiroids, chalicotheres, Chalicotheriinae, and equids). For each taxa, the phylogenetic position of known fossils was compared with the stratigraphic position of the same fossils. A positive correlation was found for all of the 24 taxa, 18 of which were statistically significant.

As a third example, Michael Benton and Rebecca Hitchin published a more recent, greatly expanded, and detailed stratigraphic analysis of 384 published cladograms of various multicellular organisms (Benton and Hitchin 1997). Using three measures of congruence between the fossil record and phylogeny (the RCI, SRC, and SCI), these researchers observed values "skewed so far from a normal distribution [i.e. randomness] that they provide evidence for strong congruence of the two datasets [fossils and cladograms]." Furthermore, Benton and Hitchin's analysis was extremely conservative, since they made no effort to exclude cladograms with poor resolution, to exclude cladograms with very small numbers of taxa, or to use only fossils with reliable dates. Including these types of data will add confounding random elements to the analysis and will decrease the apparent concordance between stratigraphy and cladograms. Even so, the results were overall extremely statistically significant (P < 0.0005). As the authors comment in their discussion:

"... the RCI and SCI metrics showed impressive left-skewing; the majority of cladograms tested show good congruence between cladistic and stratigraphic information. Cladists and stratigraphers may breathe easy: the cladistic method appears, on the whole, to be finding phylogenies that may be close to the true phylogeny of life, and the sequence of fossils in the rocks is not misleading. ... it would be hard to explain why the independent evidence of the stratigraphic occurrence of fossils and the patterns of cladograms should show such striking levels of congruence if the fossil record and the cladistic method were hopelessly misleading." (Benton and Hitchin 1997, p. 889)

Additionally, if the correlation between phylogeny and stratigraphy is due to common descent, we would expect the correlation to improve over longer geological time frames (since the relative error associated with the fossil record decreases). This is in fact observed (Benton et al. 1999). We also would expect the correlation to improve, not to get worse, as more fossils are discovered, and this has also been observed (Benton and Storrs 1994).

Potential Falsification:

It would be highly inconsistent if the chronological order were reversed in the reptile-bird and reptile-mammal example. Even the finding that there was no overall correlation between stratigraphy and the consensus phylogeny of the major taxa would be very problematic for the theory of common descent. In addition, the observed correlation could decrease over longer time frames or as we acquire more paleontological data—but neither is the case (Benton et al. 1999; Benton and Storrs 1994).

Based on the high confidence in certain branches of phylogenetic trees, some temporal constraints are extremely rigid. For example, we should never find mammalian or avian fossils in or before Devonian deposits, before reptiles had diverged from the amphibian tetrapod line. This excludes Precambrian, Cambrian, Ordovician, and Silurian deposits, encompassing 92% of the earth's geological history and 65% of the biological history of multicellular organisms. Even one incontrovertible find of any pre-Devonian mammal, bird, or flower would shatter the theory of common descent (Kemp 1982; Carroll 1988; Stanley 1993; Chaterjee 1997).



Albert, J., Wahlberg, J., Leitner, T., Escanilla, D. and Uhlen, M. (1994) "Analysis of a rape case by direct sequencing of the human immunodeficiency virus type 1 pol and gag genes." J Virol 68: 5918-24. [PubMed]

Archie, J. W. (1989) "A randomization test for phylogenetic information in systematic data." Systematic Zoology 38: 219-252.

Arnold, C., Balfe, P. and Clewley, J. P. (1995) "Sequence distances between env genes of HIV-1 from individuals infected from the same source: implications for the investigation of possible transmission events." Virology 211: 198-203. [PubMed]

Atchely, W. R., and Fitch, W. M. (1991) "Gene trees and the origins of inbred strains of mice." Science 254: 554-558. [PubMed]

Avery, O. T., MacLeod, C. M. and McCarty, M. (1944) "Studies on the chemical nature of the substance inducing transformation of pneumococcal types." J. Exp. Med. 79:137-158.

Avise, J. C., and Wollenberg, K. (1997) "Phylogenetics and the origin of species." PNAS 94: 7748-7755. http://www.pnas.org/cgi/ content/full/94/15/7748

Bailey, T. L., and Gribskov, M. (1998) "Combining evidence using p-values: application to sequence homology searches. Bioinformatics. 14: 48-54. [PubMed]

Bajpai, S. and Gingerich, P. D. (1998) "A new Eocene archaeocete (Mammalia, Cetacea) from India and the time of origin of whales." Proc Natl Acad Sci U S A. 95: 15464-15468. [PubMed]

Baldauf, S. L., Roger, A. J., Wenk-Siefert, I., and Doolittle, W. F. (2000) "A kingdom-level phylogeny of eukaryotes based on combined protein data." Science 290: 972-7. [PubMed]

Behe, M.J. (1994) "Experimental support for regarding functional classes of proteins to be highly isolated from each other." In Darwinism, Science or Philisophy? Eds. Jon Buell and Virginia Hearn, Foundation for Thought and Ethics: Houston, Texas. [ARN Michael Behe Files]

Benton, M. J. (1998) "Molecular and morphological phylogenies of mammals: Congruence with stratigraphic data." Molecular Phylogenetics and Evolution 9: 398-407. [PubMed]

Benton, M. J., and Hitchin, R. (1996) "Testing the quality of the fossil record by groups and by major habitats." Historical Biol. 12: 111-157.

Benton, M. J., and Hitchin, R. (1997) "Congruence between phylogenetic and stratigraphic data on the history of life." Proc. R. Soc. Lond. B. 264: 885-890.

Benton, M. J., Hitchin, R., and Wills, M. A. (1999) "Assessing congruence between cladistic and stratigraphic data." Syst. Biol. 48: 581-596.

Benton, M. J., and Storrs, G. W. (1994) "Testing the quality of the fossil record: paleontological knowledge is improving." Geology 22: 111-114.

Benton, M. J., Wills, M. A., and Hitchin, R. (2000) "Quality of the fossil record through time." Nature 403: 534-537. [PubMed]

Birch, C. J., McCaw, R. F., Bulach, D. M., Revill, P. A., Carter, J. T., Tomnay, J., Hatch, B., Middleton, T. V., Chibo, D., Catton, M. G., Pankhurst, J. L., Breschkin, A. M., Locarnini, S. A. and Bowden, D. S. (2000) "Molecular analysis of human immunodeficiency virus strains associated with a case of criminal transmission of the virus." J Infect Dis 182: 941-4. http://www.journals.uchicago.edu/JID/journal/issues/v182n3/000154/000154.html

Blanchard, A., Ferris, S., Chamaret, S., Guetard, D. and Montagnier, L. (1998) "Molecular evidence for nosocomial transmission of human immunodeficiency virus from a surgeon to one of his patients." J Virol 72: 4537-40. http://jvi.asm.org/cgi/content/full/72/5/4537?view=full&pmid=9557756

Box, G. E. P. (1978) Statistics for experimenters : an introduction to design, data analysis, and model building. New York: John Wiley and Sons.

Brenner, S. (1957) "On the impossibility of all overlapping triplet codes in information transfer from nucleic acid to proteins." PNAS 43:687-694.

Brown, J. R., Douady, C. J., Italia, M. J., Marshall, W. E., and Stanhope, M. J. (2001) "Universal trees based on large combined protein sequence data sets." Nature Genetics 28: 281-285. [PubMed]

Bryant, D., MacKenzie, A. and Steel, M. (2002). "The size of a maximum agreement subtree for random binary trees." In: Bioconsensus II. DIMACS Series in Discrete Mathematics and Theoretical Computer Science (American Mathematical Society). ed., M.F. Janowitz.

Carroll, R. L. (1988) Vertebrate Paleontology and Evolution. New York: W. H. Freeman and Co.

Carroll, R. L. (1997) Patterns and Processes of Vertebrate Evolution. Cambridge: Cambridge University Press.

Carter, G. S. (1954) Animal Evolution. London: Sidgwick and Jackson.

Czerkas, S.A., Zhang, D., Li, J. and Li, Y. (1954) "Flying dromaeosaurs." Dinosaur Museum Journal 1:97-126.

Chaterjee, S. (1997) The Rise of Birds: 225 million years of evolution. Baltimore, MD: Johns Hopkins University Press.

Chen, P.-j., Z.-m. Dong, et al. (1998) "An exceptionally well-preserved theropod dinosaur from the Yixian formation of China." Nature 391: 147-152.

Clyde, W. C., and Fisher, D. C. (1997) "Compaing the fit of stratigraphic and morphologic data in phylogenetic analysis." Paleobiology 23: 1-19.

Crick, F. H. C., Barnett, L., Brenner, S., and Watts-Tobin, R. J. (1961) "General nature of the genetic code for proteins." Nature 192:1227-1232.

Crick, F. H. C. and Orgel, L. E. (1973) "Directed panspermia." Icarus 19:341-346.

Cunningham, C. W., Omland, K. E., and Oakley, T. H. (1998) "Reconstructing ancestral character states: a critical reappraisal." Trends in Ecology and Evolution 13: 361-366.

Darwin, C. (1872) The Origin of Species. Sixth Edition. The Modern Library, New York. [Project Gutenberg text file]

Dickerson, R. E., and Timkovich, R. (1975) cytochrome c. The Enzymes. P. D. Boyer. New York: Academic Press. 11: 397-547.

Dodson, E. O. (1960) Evolution: Process and Product. New York: Reinhold Publishers.

Domning, D. P. (2001) "The earliest known fully quadrupedal sirenian." Nature. 413: 625-627. [PubMed]

Doolittle, W. F. (1999) "Phylogenetic Classification and the Universal Tree." Science 284: 2124. [PubMed]

Doolittle, W. F. (2000) "The nature of the universal ancestor and the evolution of the proteome." Current Opinion in Structural Biology 10: 355-358. [PubMed]

Faith, D. P., and Cranston, P. S. (1991) "Could a cladogram this short have arisen by chance alone?: on permutation tests for cladistic structure." Cladistics 7: 1-28.

Farris, J. S. (1989) "The retention index and the rescaled consistency index." Cladistics 5:417-419.

Felsenstein, J. (1982) "Numerical methods for inferring evolutionary trees." Quart. Rev. Biol. 57: 379-404.

Felsenstein, J. (1985) "Confidence limits on phylogenies: an approach using the bootstrap." Evolution 39: 783-791.

Fisher, R. A. (1990) Statistical methods, experimental design, and scientific inference. Bennet, J. H. ed., Oxford: Oxford University Press.

Fitch, W. M. (1970) "Distinguishing homologous from analogous proteins." Syst. Zool. 28: 132-163.

Futuyma, D. (1998) Evolutionary Biology. Third edition. Sunderland, MA: Sinauer Associates.

Galton, F. and Watson, H.W. (1874) "On the probability of the extinction of families." Journal of the Royal Anthropological Institute 4: 138-144.

Garavelli, J. S., Hou, Z., et al. (2001) "The RESID Database of protein structure modifications and the NRL-3D Sequence-Structure Database." Nucleic Acids Research 29: 199-201. http://www.ncifcrf.gov/RESID/

Gish, D. T. (1978) Evolution? The Fossils Say No! Public School Edition, San Diego: Creation-Life Publishers.

Gilbert, S. F. (1997) Developmental Biology. Fifth edition. Sinauer Associates.

Gingerich, P., Wells, N. A., Russell, D. E., and Shah, S. M. I. (1983) "Origin of whales in epicontinental remnant seas: new evidence from the early Eocene of Pakistan." Science 220: 403-406.

Gingerich, P., Smith, B. H., and Simons, E. L. (1990) "Hind limbs of Eocene Basilosaurus: evidence of feet in whales." Science 249: 154-157.

Gingerich, P., Raza, S. M., Arif, M., Anwar, M., and Zhou, X. (1994) "New whale from the Eocene of Pakistan and the origin of cetacean swimming." Nature 368: 844-847.

Gingerich, P. D., Haq Mu, Zalmout, I. S., Khan, I. H., and Malkani, M. S. (2001). "Origin of whales from early artiodactyls: hands and feet of Eocene Protocetidae from Pakistan." Science. 293: 2239-2242. [PubMed]

Givnish, T.J., and Sytsma, K.J. (1997) "Consistency, characters, and the likelihood of correct phylogenetic inference." Mol Phylogenet Evol. 7: 320-330. [PubMed]

Goujon, C. P., Schneider, V. M., Grofti, J., Montigny, J., Jeantils, V., Astagneau, P., Rozenbaum, W., Lot, F., Frocrain-Herchkovitch, C., Delphin, N., Le Gal, F., Nicolas, J. C., Milinkovitch, M. C. and Deny, P. (2000) "Phylogenetic analyses indicate an atypical nurse-to-patient transmission of human immunodeficiency virus type 1." J Virol 74: 2525-32. http://jvi.asm.org/cgi/content/full/74/6/2525?view=full&pmid=10684266

Gould, S. J. (1990) "An earful of jaw." Natural History 3: 12-23.

Gregory, W. K. (1951) Evolution Emerging: A survey of changing patterns from primeval life to man. Volume I. The Macmillan Co., New York.

Groves, C. P. (2000) "The genus Cheirogaleus: Unrecognized biodiversity in dwarf lemurs." International Journal of Primatology. 21(6): 943-962.

Gura, T. (2000) "Bones, molecules ... or both?" Nature 406: 230-233. [PubMed]

Haeckel, E. (1898) The Last Link. London: Adam and Charles Black.

Harris, T. E. (1989) The Theory of Branching Processes. New York: Dover.

Hedges, S. B. (1994) "Molecular evidence for the origin of birds." PNAS 91: 2621-2624. http://www.pnas.org/cgi/content/abstract/91/7/2621

Hedges, S. B., and Poling, L. L. (1999) "A molecular phylogeny of reptiles." Science 283: 998-1001. [PubMed]

Hendy, M. D., Little, C. H. C., and Penny, D. (1984) "Comparing trees with pendant vertices labelled." SIAM J. Appl. Math. 44: 1054-1065.

Hillis, D. M. (1991) "Discriminating between phylogenetic signal and random noise in DNA sequences." In Phylogenetic analysis of DNA sequences. pp. 278-294 M. M. Miyamoto and J. Cracraft, eds. New York: Oxford University Press.

Hillis, D. M., and Bull, J. J. (1993) "An empirical test of bootstrapping as a method for assessing confidence on phylogenetic analysis." Syst. Biol. 42: 182-192.

Hillis, D. M., and Huelsenbeck, J. P. (1992) "Signal, noise, and reliability in molecular phylogenetic analyses." Journal of Heredity 83: 189-195. [PubMed]

Hillis, D. M., Moritz, C. and Mable, B. K. Eds. (1996) Molecular systematics. Sunderland, MA: Sinauer Associates.

Hinegardner, T. T., and Engelberg, J. (1963) "Rationale for a Universal Genetic Code." Science 142: 1083-1085.

Hitchin, R., and Benton, M. J. (1997) "Congruence between parsimony and stratigraphy: comparisons of three indices." Paleobiology 23: 20-32.

Hofreiter, M., Serre, D., Poinar, H. N., Kuch, M., and Paabo, S. (2001) "Ancient DNA." Nat Rev Genet 2: 353-359.

Holmes, E. C., Zhang, L. Q., Simmonds, P., Rogers, A. S. and Brown, A. J. (1993) "Molecular investigation of human immunodeficiency virus (HIV) infection in a patient of an HIV-infected surgeon." J Infect Dis 167: 1411-4. [PubMed]

Hudson, R. R. (1992) "Gene trees, species trees and the segregation of ancestral alleles." Genetics 131: 509-513. [PubMed]

Huelsenbeck, J. P. (1994) "Comparing the stratigraphic record to estimates of phylogeny." Palaeobiology 20: 470-483.

Huelsenbeck, J. P., Ronquist, F., Nielsen, R., and Bollback, J. P. (2001) "Bayesian inference of phylogeny and its impact on evolutionary biology." Science 294: 2310-2314. [PubMed]

Ji, Q., Norell, M. A., Gao, K-Q., Ji, S-A., and Ren, D. (2001) "The distribution of integumentary structures in a feathered dinosaur." Nature 410: 1084-1088. [PubMed]

Judson, H. F. (1996) The Eighth Day of Creation: Makers of the Revolution in Biology. Expanded edition. Cold Spring Harbor Press.

Kardong, K. V. (2002) Vertebrates: Comparative Anatomy, Function, Evolution. Third ed. New York: McGraw Hill.

Klassen, G. J., Mooi, R. D., and Locke, A. (1991) "Consistency indices and random data." Syst. Zool. 40:446-457.

Kemp, J. S. (1982) Mammal-like reptiles and the origin of mammals. New York: Academic Press.

Kestenbaum, D. (1998) "Gravity Measurements Close in on Big G." Science 282: 2180-2181. http://www.sciencemag.org/cgi/content/full/282/5397/2180

Krings, M., Capelli, C., Tschentscher, F., Geisert, H., Meyer, S., von Haeseler, A., Grossschmidt, K., Possnert, G., Paunovic, M., and Paabo, S. (2000) "A view of Neandertal genetic diversity." Nat Genet 26: 144-146.

Lehman, N. (2001) "Please release me, genetic code." Current Biology 11: R63-R66. [PubMed]

Li, W.-H. (1997) Molecular Evolution. Sunderland, MA: Sinauer Associates.

Luo, Z.-X., Crompton, A. W. and Sun, A.-L. (2001) "A new mammaliaform from the early Jurassic and evolution of mammalian characteristics." Science 292: 1535-1539. [PubMed]

Lyons-Weiler, J., Hoeker, G. A., and Tausch, R. J. (1996) "Relative Apparent Synapomorphy Analysis (RASA) I: The Statistical Measurement of Phylogenetic Signal." Mol. Biol. Evol. 13: 749-757. http://www. molbiolevol.org/cgi/content/abstract/13/6/749

Machuca, R., Jorgensen, L. B., Theilade, P. and Nielsen, C. (2001) "Molecular investigation of transmission of human immunodeficiency virus type 1 in a criminal case." Clin Diagn Lab Immunol 8: 884-90. [PubMed]

Maddison, W. P., and Maddison, D. R. (1992) MacClade. Sunderland, MA: Sinauer Associates.

Maley, L. E., and Marshall, C. R. (1998) "The coming age of molecular systematics." Science 279: 505-506. [PubMed]

McLaughlin, P. J., and Dayhoff, M. O. (1973) "Eukaryote evolution: a view based on cytochrome c sequence data." Journal of Molecular Evolution 2: 99-116.

Nichols, R. (2001) "Gene trees and species trees are not the same." Trends Ecol Evol. 16: 358-364. [PubMed]

Norell, M. A., and Clarke, J. A. (2001) "Fossil that fills a critical gap in avian evolution." Nature 409: 181-184. [PubMed]

Norell, M. A., Ji, Q., Gao, K., Yuan, C., Zhao, Y., and Wang, L. (2002) "'Modern' feathers on a non-avian dinosaur." Nature 416: 36-37. [PubMed]

Norell, M. A. and Novacek, M. J. (1992a) "Congruence between superpositional and phylogenetic patterns: Comparing cladistic patterns with fossil records." Cladistics 8: 319-337.

Norell, M. A. and Novacek, M. J. (1992b) "The fossil record and evolution: Comparing cladistic and paleontologic evidence for vertebrate history." Science 255: 1690-93.

Norris, J. R. (1997) Markov Chains. Cambridge: Cambridge University Press.

Osborn, H. F. (1918) The Origin and Evolution of Life. New York: Charles Scribner's Sons.

Ou, C. Y., Ciesielski, C. A., Myers, G., Bandea, C. I., Luo, C. C., Korber, B. T., Mullins, J. I., Schochetman, G., Berkelman, R. L., Economou, A. N. and et al. (1992) "Molecular epidemiology of HIV transmission in a dental practice." Science 256: 1165-71. [PubMed]

Patterson, C., Williams, D., and Humphries, C. (1993) "Congruence between molecular and morphological phylogenies." Annual Review of Ecology and Systematics 24: 153-188.

Pei, M. (1949) The Story of Language. Philadelphia: Lippincott.

Penny, D., and Hendy, M. D. (1986) "Estimating the reliability of phylogenetic trees." Mol. Biol. Evol. 3: 403-417. http://www. molbiolevol.org/cgi/content/abstract/3/5/403

Penny, D., Foulds, L. R., and Hendy, M. D. (1982) "Testing the theory of evolution by comparing phylogenetic trees constructed from five different protein sequences." Nature 297: 197-200.

Penny, D., Hendy, M. D., and Steel, M. A. (1991) "Testing the theory of descent." In Phylogenetic Analysis of DNA Sequences, eds. Miyamoto, M. and Cracraft, J., New York: Oxford University Press. pp. 155-183.

Qiang, J., Currie, P. J. et al. (1998) "Two feathered dinosaurs from northeastern China." Nature 393: 753-761.

Quinn, T. (2000) "Measuring big G." Nature 408: 919-921.

Rasoloarison, R. M., Goodman, S. M., and Ganzhorn, J. U. (2000) "Taxonomic revision of mouse lemurs (Microcebus) in the western portions of Madagascar." International Journal of Primatology. 21(6): 963-1019.

Ringe, D. (1999) "Language classification: scientific and unscientific methods." in The Human Inheritance, ed. B. Sykes. Oxford: Oxford University Press, pp. 45-74.

Rozenski, J., Crain, P. F. et al. (1999) "The RNA Modification Database: 1999 update." Nucleic Acids Research 27: 196-197. http://medlib.med.utah.edu/RNAmods/

Sereno, P. C. (1999) "The Evolution of Dinosaurs." Science 284: 2137-2147. [PubMed]

Stanley, S. (1993) Earth and Life Through Time. New York: W. H. Freeman.

Steel, M. A., and Penny,D. (1993) "Distributions of tree comparison metrics—some new results." Systematic Biology 42: 126-141.

Strait, D. S., Grine, F. E., and Moniz, M. A. (1997) "A reappraisal of early hominid phylogeny." J Hum Evol 32: 17-82.

Swofford, D. L., Olsen, G. J., Waddell, P. J., and Hillis, D. M. (1996) "Phylogenetic inference." In Molecular Systematics, pp 407-514. Hillis, D. M., Moritiz, C. and Mable, B. K. eds., Sunderland, Massachusetts: Sinauer.

Thalmann, . U., and Geissmann, T. (2000) "Distribution and geographic variation in the western woolly lemur (Avahi occidentalis) with description of a new species (A. unicolor)." International Journal of Primatology. 21(6): 915-941.

Thewissen, J. G, Williams, E. M., Roe, L. J., and Hussain, S. T. (2001) "Skeletons of terrestrial cetaceans and the relationship of whales to artiodactyls." Nature. 413: 277-281. [PubMed]

van Roosmalen, M. G. M., van Roosmalen, T., Mittermeier, R. A., and Rylands, A. B. (2000) "Two new species of marmoset, genus Callithrix Erxleben, 1777 (Callitrichidae, Primates), from the Tapajos/Madeira interfluvium, south central Amazonia." Neotropical Primates 8(1): 2-18.

Veenstra, J., Schuurman, R., Cornelissen, M., van't Wout, A. B., Boucher, C. A., Schuitemaker, H., Goudsmit, J. and Coutinho, R. A. (1995) "Transmission of zidovudine-resistant human immunodeficiency virus type 1 variants following deliberate injection of blood from a patient with AIDS: characteristics and natural history of the virus." Clin Infect Dis 21: 556-60. [PubMed]

Voet, D., and Voet, J. (1995) Biochemistry. New York: John Wiley and Sons.

Vogel, G. (1997) "Phylogenetic analysis: getting its day in court." Science 275: 1559-60. [PubMed]

Wadsworth, H. M. (1997) Handbook of statistical methods for engineers and scientists. New York: McGraw-Hill Pub. Co.

Wang, Y. Hu, Y., Meng, J. and Li,C. (2001) "An ossified meckel's cartilage in two Cretaceous mammals and origin of the mammalian middle ear." Science 294: 357-361. [PubMed]

Wills, M. A. (1999) "Congruence between phylogeny and stratigraphy: Randomization tests and the gap excess ratio." Syst. Biol. 48: 559-580.

Wilson, E. O. (1992) The Diversity of Life. Cambridge, MA, Harvard University Press.

Wise, K. P. (1994) "The Origin of Life's Major Groups." In The Creation Hypothesis, pp. 211-234. Moreland, J. P. ed., Downers Grove, IL: InterVarsity Press.

Woese, C. (1998) "The universal ancestor." PNAS 95: 6854-6859. http://www.pnas.org/cgi/ content/full/95/12/6854

Wu, C. I. (1991) "Inferences of species phylogeny in relation to segregation of ancient polymorphisms." Genetics 127: 429-435. [PubMed]

Xu, X., Tang, Z-T., Wang, X-l. (1999) "A therizinosauroid dinosaur with integumentary structures from China." Nature 399: 350-354.

Xu, X., Zhou, Z., and Wang ,X. (2000) "The smallest known non-avian theropod dinosaur." Nature 408: 705-708 [PubMed]

Xu, X., Norell, M. A., Wang, X-l., Makovicky, P. J., and Wu, X.C. (2002) "A basal troodontid from the Early Cretaceous of China." Nature 415: 780-784. [PubMed]

Yirrell, D. L., Robertson, P., Goldberg, D. J., McMenamin, J., Cameron, S. and Leigh Brown, A. J. (1997) "Molecular investigation into outbreak of HIV in a Scottish prison." Bmj 314: 1446-50. http://bmj.com/cgi/content/full/314/7092/1446?view=full&pmid=9167560

Yockey, H. P. (1992) Information Theory and Molecular Biology. New York: Cambridge University Press.

Yoder, A. D., Rasoloarison, R. M., Goodman, S. M., Irwin, J. A., Atsalis, S., Ravosa, M. J., and Ganzhorn, J. U. (2000) "Remarkable species diversity in Malagasy mouse lemurs (primates, Microcebus)." Proc Natl Acad Sci 97(21): 11325-30. http://www.pnas.org/cgi/ content/full/97/21/11325

Zhu, T., B. Korber, et al. (1998) "An African HIV-1 sequence from 1959 and implications for the origin of the epidemic." Nature 391: 594-597. [PubMed]

Zuckerkandl, E. and Pauling, L. (1965) "Evolutionary Divergence and Convergence in Proteins." in Evolving Genes and Proteins: a symposium held at the Institute of Microbiology of Rutgers, with support from the National Science Foundation. Eds Vernon Bryson and Henry J. Vogel. New York: Academic Press.

Full references