<<

Received 24 September 2001 Accepted 3 December 2001 Published online 9 April 2002

A -level supertree of the Dinosauria Davide Pisani1,2*, Adam M. Yates1, Max C. Langer1 and Michael J. Benton1 1Department of Earth Sciences, University of Bristol, Bristol BS8 1RJ, UK 2Department of Zoology, The Natural History Museum, Cromwell Road, SW7 5BD, UK One of the ultimate aims of systematics is the reconstruction of the tree of life. This is a huge undertaking that is inhibited by the existence of a computational limit to the inclusiveness of phylogenetic analyses. Supertree methods have been developed to overcome, or at least to go around this problem by combining smaller, partially overlapping . Here, we present a very inclusive generic-level supertree of Dinosauria (covering a total of 277 genera), which is remarkably well resolved and provides some clarity in many contentious areas of systematics. Keywords: supertrees; phylogeny; ; ;

1. INTRODUCTION tations of the information conveyed by the trees in the input set. In the last few decades, advances in theoretical system- The original MRP method of Baum (1992) and Ragan atics, as well as in computer sciences, have provided bio- (1992), referred to as component coding-MRP (CC- logists and palaeontologists with tools to investigate the MRP) in this paper, was used to combine 126 `partial phy- phylogenetic relationships among organisms, yet phylo- logenies’ of dinosaurs (covering a total of 277 genera) with genetic inference is computationally expensive and analy- the aim of obtaining a single, genus-level supertree. CC- ses of large datasets hardly feasible (Graham & Foulds MRP makes use of the fact that a tree can be represented 1982). This problem is exacerbated in palaeontology, in a matrix by means of a set of binary `characters’. Fol- where missing data can lead to an exponential increase in lowing Baum & Ragan (1993), we refer to these `charac- the number of most-parsimonious solutions permitted by ters’ as matrix elements. Matrix representations of a given dataset. However, `partial phylogenies’ (i.e. phy- different source trees can be readily combined with miss- logenies including only a subset of the taxa known to ing entries in the matrix elements corresponding to taxa belong to a certain group) are accumulating at an increas- that are not present in the matching source tree. Parsi- ing rate, and the Dinosauria represent a typical, if not mony analysis of the combined matrix representations extreme, example of this trend. yields one or more most parsimonious trees (MPTs), the Supertree methods, formally introduced by Gordon CC-MRP supertree(s), which comprise all the included (1986), can be implemented to combine such `partial phy- taxa (see Baum 1992; Ragan 1992 for more details). logenies’ and obtain more inclusive estimates without the need to pool the original datasets. Accordingly, whenever a collection of `partial phylogenies’ exists, supertrees can 2. MATERIAL AND METHODS represent an alternative to a direct analysis of the primary Potential source trees were identi® ed from online searches. data (i.e. the pooled dataset). Two electronic literature databases were searched for publications Several supertree methods have been developed (e.g. with relevant source trees, BFV (Bibliography of Fossil Verte- Sanderson et al. 1998; Wilkinson et al. 2001 for review), brates) online (http://eteweb.lscf.ucsb.edu/bfv/bfv_form.html) but only those using matrix representation with parsimony and Web of Science (http://wos.mimas.ac.uk/). BFV online does (MRP methods) have software implementation (Swofford not allow the enquirer to con® ne his search to a limited number 1998; Thorley & Page 2000) and these are currently the of and it was thus screened over the whole time-span that only viable choice in supertree reconstruction. Neverthe- it covers (1509± 1993), while Web of Science was searched for less, their characterization is still incomplete. In particular, the years 1993± 2001. All publications that were likely to include when the source trees are incompatible (i.e. they disagree a were examined. Because the aim of this analysis about the relationships of the taxa that overlap), their was to obtain a consensus representation of cladistic views on behaviour is not completely understood. In these cases, dinosaur systematics, the earliest of which date back to the MRP methods can be inaccurate. It is debatable whether beginning of the 1980s (Benton 1990), publications presented MRP estimates are accurate enough to be useful tree- before the 1980 have been excluded a priori. Publications reconstruction methods (Purvis 1995a; Wilkinson et al. after this date were retained when they presented cladogram(s) 2001). Yet, there is evidence (Bininda-Emonds & San- resulting from an original study, or from the modi® cation of a derson 2001) that, if the set of source trees is large pre-existing dataset. The list of all retained publications is avail- enough, the MRP supertree(s) are accurate represen- able (see electronic appendix). For each selected cladogram, a Nexus (Swofford 1998) tree ® le was produced using MacClade (Maddison & Maddison 1997). Some terminal genera were omitted a priori (see electronic appendix, available on The Royal * Author for correspondence ([email protected]). Society’s Publications Web site, for the list of the omitted taxa).

Proc. R. Soc. Lond. B (2002) 269, 915± 921 915 Ó 2002 The Royal Society DOI 10.1098/rspb.2001.1942 916 D. Pisani and others A supertree of Dinosauria

These were either nomina dubia, unconstrained taxa, or both. of support, following Bininda-Emonds et al. (1999) and Liu et The former are taxa based on material inadequate for diagnosis al. (2001), decay indices (Bremer 1988) are provided. Decay and for proper phylogenetic analysis. The latter may be valid values have been calculated for Matrices II± IV, but because of taxa but they appear only in one source tree and as members the size of the resulting supertree, they have not been calculated of a polytomy. Accordingly, their inclusion adds nothing to the for Matrix I. However, it is our opinion that this measure cannot analysis but an increase in the number of possible MPTs. Most be considered really representative of the support for the original studies included a variable proportion of suprageneric recovered by the supertree. Accordingly, decay values should be terminal taxa. For analytical proposes, their was interpreted with care. Consider, for example, that novel clades assumed and a standard substitution (Wilkinson et al. (sensu Bininda-Emonds & Bryant 1998), despite being unsup- 2001) was performed. Each suprageneric taxon was substituted ported (Pisani & Wilkinson 2002), have a positive decay index. with a star tree (i.e. an unresolved polytomy) including either: Thus, new support measures for MRP supertrees are needed (i) the taxa embodied in each of them when explicitly given in and are under development (D. Pisani, unpublished data), but the original study; or (ii) only its non-controversial members in they will not be considered any further in this paper. the case that no such de® nition was given (see electronic appen- dix for the minimal non-controversial de® nition of each supra- 3. RESULTS generic taxon). When multiple MPTs were presented, strict reduced consensus (RC; Wilkinson 1994) was implemented A total of 34 900 MPTs was found from the analysis of using RadCon (Thorley & Page 2000) and the tree in the RC Matrix I before the search had to be aborted because of pro® le with the highest cladistic information content (Thorley memory limitations. We suggest that this sample of et al. 1998) was retained for the CC-MRP analysis. MPTs, although less than the set of all possible supertrees, Five CC-MRP matrices (see electronic appendix), one scoring is representative of the whole. This is because the topology all of the 277 included taxa (Matrix I), one for Sauropodomor- of their strict consensus agrees perfectly with those of the pha (Matrix II), one for (Matrix III) and two for strict consensus trees summarizing the results obtained (Matrices IV and V) were produced using RadCon. from the compartmentalized analyses (Matrices II, III and These have been analysed separately, resulting in one compre- V), which have been run to completion. It is thus possible hensive (Matrix I) and four compartmentalized (Matrices II± V) to argue that increasing the number of retained MPTs supertrees. It is important to note that Matrices II± V are matr- from the analysis of Matrix I will result in a more compre- ices in their own right, not subsamples of Matrix I. They have hensive collection of alternative solutions for the polytom- not been obtained by pruning the excluded taxa, e.g. all thero- ies seen in ® gure 2 and not in the collapse of its resolved pods and all ornithischians (in the case of Matrix II) from Matrix areas. The strict consensus of the 34 900 MPTs from the I, but by pruning the source trees and then recoding them in a analysis of Matrix I (® gure 2) therefore represents the new CC-MRP matrix. This approach, implemented here for the most inclusive estimate that has ever, to our knowledge, ® rst time, to our knowledge, is important because each matrix been presented of the phylogeny of the Dinosauria. element is representative of a node in a tree. Thus, pruning a Contentious issues in dinosaur systematics are resolved taxon from an MRP matrix will create a matrix that is not rep- in the present supertree. is composed resentative of the real topology of the pruned tree (® gure 1). of Saturnalia and monophyletic Prosauropoda and Sauro- Because objective functions for implementing a correct, poda (Sereno 1999; Benton et al. 2000; contra Gauthier unequal character-weighting scheme in MRP analyses are still 1986). Within Prosauropoda, a monophyletic Melanoro- unknown (Bininda-Emonds et al. 1999; Bininda-Emonds & sauridae (, and ) Sanderson 2001), equally weighted parsimony was used to ana- emerges as the of (Masso- lyse the ® ve CC-MRP matrices. spondylus, and ). Within

Paup¤ (Swofford 1998) was used to analyse the data matrices , is the sister group to all other with parsimony and, in to avoid the problem of being (Upchurch 1998; contra Wilson & Sereno `bogged down’ in sub-optimal tree islands, the following tree- 1998) and is more closely related to Neosauro- search protocol was applied: (1) 100 heuristic searches were poda than it is to (Wilson & Sereno 1998; carried out using tree bisection± reconnection (TBR) branch contra Upchurch 1998). The early dinosaurs and swapping with the `multrees’ option not in effect; random step- the are theropods (Sereno et al. wise addition was used to obtain the starting trees for these repli- 1993; contra Langer et al. 1999). is a mono- cates; (2) the MPTs found were retained and swapped using phyletic group containing and Neocerato- TBR branch swapping with the `multrees’ option in effect. (Sereno 1999; Holtz 2000; contra Forster 1999), A potentially serious problem of MRP supertree methods is while torvosaurids are more closely related to derived their ability to resolve con¯ ict among the source trees by the tetanurans () than to spinosaurids (Holtz generation of novel clades (sensu Bininda-Emonds & Bryant 2000; contra Sereno et al. 1994). Therizinosauroidea is the 1998). Bininda-Emonds et al. (1999) and Bininda-Emonds & sister group of (Russell & Dong 1993; Sanderson (2001) have likened the appearance of such clades in Holtz 2000; contra Sereno 1999), but is a MRP to their appearance in Total Evidence analyses. However, basal member of : it is excluded from both Ovira- as shown by Pisani & Wilkinson (2002), this cannot be the case ptorosauria (contra Sereno 1999) and (contra Ji and such novel clades should always be collapsed. No novel et al. 1998), while Metornithes is monophyletic and con- clades (sensu Bininda-Emonds & Bryant 1998) are present in the tains , (Chatterjee 1991; contra dinosaur supertree. Holtz 1994) and (Chiappe et al. 1996; A poorly understood problem in MRP supertree building is contra Sereno 1999). how to evaluate the support for the resulting clades. For those Within Ornithischia, Minmi and Gargoyleosaurus are who feel that an MRP supertree is naked without some measure basal Ankylosauridae (Sereno 1999; contra Kirkland

Proc. R. Soc. Lond. B (2002) A supertree of Dinosauria D. Pisani and others 917

(tree A) (tree B) A B C D E F A B C E F

node 1 node 1 node 2 node 3 node 3 pruning taxon D node 4 node 4

CC-MRP representation of tree A CC-MRP representation of tree B

nodes: N1 N2 N3 N4

A 0 0 0 0 B 0 0 0 1 C 0 0 1 1 The deletion of taxon D implies the D 0 1 1 1 E 1 1 1 1 suppression of node 2. F 1 1 1 1 Accordingly, the corresponding matrix element (N2) is not present in the correct CC-MRP representation Matrix obtained when pruning of tree B. taxon D from the CC-MRP representation of tree A (i.e. using the PAUP option ‘Delete-

Restore taxa’). Note that the

new matrix is not a correct CC- MRP representation of tree B.

nodes: N1 N2 N3 N4 nodes: N1 N3 N4

A 0 0 0 0 A 0 0 0 B 0 0 0 1 B 0 0 1 C 0 0 1 1 C 0 1 1 E 1 1 1 1 E 1 1 1 F 1 1 1 1 F 1 1 1

Figure 1. A tree on six (tree A) and its corresponding pruned version (tree B), from which taxon D has been deleted. Pruning taxon D from the CC-MRP representation of tree A results in a matrix that is not the correct CC-MRP matrix representing tree B.

1998). A monophyletic (Norman 1998; three long-necked Chinese taxa, Omeisaurus, Mamenchi- contra Sereno 1999) containing , Altirhinus and saurus and , and to the poorly constrained pos- is recovered, while the Hadrosauridae ition of . The ® rst problem re¯ ects differences appears more closely related to than to in the results obtained by different authors on the relation- Iguanodontidae. Leptoceratops and Udanoceratops are ships of those taxa. Among the source trees, two classes excluded from Coronosauria (Protoceratopsidae plus of cladogram exist. Upchurch (1995, 1998) includes all Ceratopsoidea) and a monophyletic Protoceratopsidae three Chinese taxa in a nested outside Neosauro- containing Bagaceratops, Breviceratops, Graciliceratops and poda, while Wilson & Sereno (1998; also Sereno 1999) is recovered (Sereno 2000; contra Chinnery & do not include and consider Omeisaurus Weishampel 1998). and Euhelopus to be distantly related. In the supertree, The present supertree contains three important poly- Omeisaurus and Euhelopus are nested according to the tomies. One lies within Sauropoda and the other two second of these two hypotheses, but Mamenchisaurus is within Theropoda: one at the base of and free to cluster with either of the two. Implementing RC, the other within Coelurosaria, at the base of the Eumanir- two `rogue’ taxa (Cetiosaurus and Mamenchisaurus) are aptora. identi® ed. With their exclusion (® gure 3a), the following The polytomy within Sauropoda involves the taxa more taxonomic statements are true in all the supertrees: (i) derived than Shunosaurus and is due to both the existence falls within (de® ned by of con¯ icting hypotheses concerning the relationships of the node connecting and ) and is

Proc. R. Soc. Lond. B (2002) 918 D. Pisani and others A supertree of Dinosauria

Saturnalia 1 P

Melanorosaurus r 1 1 o 1 Riojasaurus s 7 a

1 Camelotia u

1 Yunnanosaurus r o

Massospondylus p o 1 d

Lufengosaurus a 1 Riojasaurus 1 Sellosaurus 1 1 1 Barapasaurus Shunosaurus 1 Omeisaurus Mamenchisaurus 1 Cetiosaurus 1 Haplocantosaurus Gongbusaurus 2 1 Euhelopus D 2 i a p i 1 l r

Dacentrurus 1 o

u

3 d

a Toujiangosaurus

1 o s 3

Lexovisaurus c o 1 2 o g Stegosaurus 2 i e 1 d t Wuerhosaurus

Barosaurus e S 1

Kentrosaurus a 1 Diplodocus Apatosaurus Pawpawsaurus 2 2 1 3 Silvisaurus 1

2 T a 4

i 2 1 i r t

Gargoyleosaurus a u n a Minmi 1 Opisthocoelocaudia

s 1 o 1 Peiropolis titanosaur o Shamosaurus 1 s l Tsagantegia a

y 1 1 1 u k r i n 1 Mymoorapelta a A 2 1 1 Saltasaurus Opistocelicaudia 1 Pinacosaurus 1 Nequensaurus Euoplocephalus Eoraptor 2 1 Tarchia 9 1 2 Nodocephalosaurus a 3 1 i r 2 2 u a 1 C s 2

o 1 e

l 1 7 4 Syntarsus r a Ornatotholus 3 a h 3 t o p s e 1 a c 2 1 i Dwarf pachicephalosaur 2 3 u

h r i

c Prenocephale 1 1 a a

P Noasaurus Stygimoloch 1 1 15 Abelisaurus Chaoyangsaurus S e 4 Majungatholus p a Archaeoceratops 1 i

d Carnotaurus n

i Leptoceratops 2 2 o

s

Udanoceratops s

p 2 1 a

o Bagaceratops

1 u t r a Breviceratops 1 4 i r 2 1 d e Protoceratops 1 a

c Graciliceratops 1 9 e

Protoceratops o t 2 o

r 1 3 P Turanoceratops 1 Baryonyx 1 Anchiceratops 1 2 1 2 3 1

e 1 a Arrhinoceratops 1 d C i 1 2 s a

p Avaceratops 1

Sinraptor r o n t Szechuanosaurus o a 1 1 r Monoclonius

Yangchuanosaurus s e

2 a Centrosaurus C u

Styracosaurus 3 2 1 r

1 i Einiosaurus a Achelousaurus 1 2 2 1 1

a

i t

n o 1

d Othnielia

o Yandusaurus 2 h 1 3 Bagarataan

p Zephyrosaurus o 1 l 1 i

s

p 4 T Parksosaurus 1 y Thescelosaurus y

Tenontosaurus r H a n e 2 1 n a o

d 1 1 1 s i Horseshoe Canyon tyrannosaurid a t

1 u n Iguanodon 1 r o 1 i

d Altirhinus d

2 2 1 O o Ouranosaurus 1 1 Tyrannosaurus a e n r

Harpymimus n a Probactrosaurus 3 i u

Protohadros Pelecanimimus t

Ouranosaurus h g

I 1

Telmatosaurus 3 o

Gilmoreosaurus 2 1 m

Eolambia 1 i m 1 1 1 o Pararabdodon 1 s 1 a 1 1 u

Charonosaurus 1 1 r

1 i a e 1 Dromiceomimus a Quarry 9 coelurosaur d 2 i 2 Galimimus r 4 u O a

s v 2 1 1 o Aralosaurus i r r 2

Parasaurolophus a d 2 1 p a 1

Shantungosaurus t

1 El Brete oviraptosaur o

H 2

Anatotitan s 1 a u

Lophorhothon 2 r 1 1 i 1 Ingenia a 1 1 Kritosaurus 1 1 Caudipteryx Therizinosaurus 1 1 2 1 1 1 E u

1 Dornogov troodontid m

Sinornithoides a n

1 i r

Troodon a 1 p 1 t o

1 r

Ornithothoraces a Avimimus Rahonavis 1 1 1 Hell Creek alvaresaurid 1 Parvicursor 1

Proc. R. Soc. Lond. B (2002) A supertree of Dinosauria D. Pisani and others 919

(a) Shunosaurus Omeisaurus Patagosaurus Camarasaurus Haplocanthosaurus Euhelopus Phuwiangosaurus Nemegtosaurus Quaesitosaurus Chubutisaurus

(b) Coelurus Proceratosaurus Ornitholestes Scipionyx Bagarataan Sinosauroptreryx Compsognathus Oviraptosauria Eumaniraptoria

Figure 3. (a) Reduced consensus of the 14 323 MPTs of length 356 obtained from the analyses of Matrix II. (b) Reduced consensus of the 85 900 MPTs of length 910 obtained from the analysis of Matrix V. more closely related to Camarasaurus, Brachiosaurus and RC identi® es these taxa and their omission allows extra Titanosauria than to ; (ii) Euhelopus, Chubut- resolution of the supertree (® gure 3b). Conversely, the isaurus, Phuwiangosaurus, Nemegtosaurus, Quaesitosaurus variable placement of eumaniraptorians in the source trees and Titanosauria are more closely related to each other seems to be due to real character incongruences (Holtz than they are to any other retained sauropod; (iii) Nemeg- 2000) and leads to an unresolved polytomy in the super- tosaurus and Quaesitosaurus are more closely related to the tree topology. derived Titanosauria than they are to Diplodocidae. Within Theropoda, the polytomy at the base of Coelu- rosauria is simply caused by a few poorly constrained gen- 4. DISCUSSION era (i.e. genera that are present in few trees and are not tightly bracketed by taxa of a relatively stable position) Supertree techniques are now becoming commonplace in biology (e.g. Purvis 1995b; Bininda-Emonds et al. 1999; Deltadromeus, Dryptosaurus and Nqwebasaurus. As above, Liu et al. 2001), but this is the ® rst time, to our know- ledge, that a supertree method has been implemented to Figure 2. (Opposite.) The strict consensus of the 34 900 disentangle the phylogeny of a fossil group. MPTs of length 1988 obtained from the analysis of Matrix Supertree methods resemble meta-analysis in several I. Numbers along the branches correspond to the decay respects (see also Sanderson et al. 1998). In meta-analysis, values of the supertrees obtained in the compartmentalized formal statistical techniques are implemented to sum up a analyses (which are not reported). body of separate (but similar) experiments (Mann 1990).

Proc. R. Soc. Lond. B (2002) 920 D. Pisani and others A supertree of Dinosauria

Similarly, supertree methods exploit either the existence The authors thank O.R.P. Bininda-Emonds and M.J. San- of a one-to-one relation between a tree and its matrix rep- derson for providing a preprint of their simulation work and resentation (Ragan 1992), or subtrees± tree reconstruction three anonymous referees for their useful comments on the manuscript. A.M.Y. and M.J.B. were funded by the Lever- algorithms (e.g. Aho et al. 1981; Semple & Steel 2000) hulme Trust, D.P. was partially funded by the University of to combine a body of different (but partially overlapping) Bristol through the Bob Savage Memorial Fund, and M.C.L. phylogenies. Meta-analysis is like an ordinary scienti® c was supported by an ORS grant. review of research, except that ordinary reviews provide a qualitative, and often subjective, assessment of few stud- REFERENCES ies, while meta-analysis promises a quantitative synthesis of all available data (Mann 1990). In the same way, super- Aho, A. V., Savig, Y., Szymansky, T. G. & Ullman, J. D. 1981 tree methods introduce objectivity to phylogenetic reviews Inferring a tree from lowest common ancestors with an by quantitatively synthesizing results of previous phylo- application to the optimization of relational expressions. genetic studies. SIAM J. Comput. 10, 405± 421. Baum, B. R. 1992 Combining trees as a way of combining data Phylogenetic trees should be used as frameworks to test sets for phylogenetic inference, and the desirability of com- evolutionary hypotheses (Felsenstein 1985; MacLeod bining gene trees. Taxon 41, 3± 10. 2001). Such frameworks should be (i) obtained indepen- Baum, B. R. & Ragan, M. 1993 Reply to A.G. Rodrigo’s `A dently from the characters where evolution is under inves- comment on Baum’s method for combining phylogenetic tigation; (ii) objective (i.e. they should not represent the trees’. Taxon 42, 637± 640. personal view of the author); (iii) inclusive of all relevant Benton, M. J. 1990 Origin and interrelationships of dinosaurs. taxa. Biologists often solve point (i) by testing the evol- In The Dinosauria (ed. D. B. Weishampel, P. Dodson & H. ution of morphological characters on molecular phy- Osmo lska), pp. 11± 30. Berkeley, CA: University of Califor- logenies, an option unavailable to palaeontologists. nia Press. Benton, M. J., Juul, L., Storrs, G. W. & Galton, P. M. 2000 However, since phylogenetic supertrees, including MRP Anatomy and systematics of the prosauropod dinosaur Thec- supertrees, represent consensus solutions (Pisani & Wilk- odontosaurus antiquus from the Upper of southern inson 2002), they are only indirectly based on any parti- . J. Vertebr. Palaeontol. 20, 77± 108. cular morphological characters. Moreover, they are Bininda-Emonds, O. R. P. & Bryant, H. N. 1998 Properties of objective (point (ii)) and can be very inclusive (point (iii)). matrix representation with parsimony analyses. Syst. Biol. Therefore, they should be used as frameworks (see also 47, 497± 508. Lanyon 1993) for evolutionary palaeobiological studies, Bininda-Emonds, O. R. P. & Sanderson, M. J. 2001 Assess- such as implementation of the comparative method of ment of the accuracy of matrix representation with parsi- Felsenstein (1985), or the phylogenetic bracket of Witmer mony analysis supertree construction. Syst. Biol. 50, 565± 579. (1995). Recent evolutionary studies within Dinosauria, Bininda-Emonds, O. R. P., Gittleman, J. L. & Purvis, A. 1999 (e.g. Carrano 2000; Hutchinson 2001a,b) have generally Building large trees by combining phylogenetic information: used those authors’ own subjective estimation of the con- a complete phylogeny of the extant Carnivora (Mammalia). sensus. Biol. Rev. 74, 143± 175. With regard to the supertree under discussion, it is Bremer, K. 1988 The limits of amino acid sequence data in important to note that some of its structure can be angiosperm phylogenetic reconstruction. Evolution 42, expected to change in the near future as works in the 795± 803. press, or in progress, are published. In particular, the Carrano, M. T. 2000 Homoplasy and the evolution of dino- saur locomotion. Palaeobiology 26, 489± 512. monophyly of Prosauropoda, Ceratosauria and Metor- Chatterjee, S. 1991 Cranial anatomy and relationships of a new nithes may not be sustained, while Eoraptor and Herrera- Triassic from . Phil. Trans. R. Soc. Lond. B 332, sauridae may be found to be more basal within 277± 346. Dinosauria. Chiappe, L. M., Norrell, M. A. & Clark, J. M. 1996 Phylogen- As previously stated, Bremer support values must be etic position of Mononykus (Aves: Alvarezsauridae) from the interpreted with care when used as a measure of support Late of the . Mem. Queensland Mus. for a supertree. However, the low values obtained for a 39, 557± 582. high proportion of nodes shows that there is a certain Chinnery, B. J. & Weishampel, D. B. 1998 Montanoceratops amount of discordance among the source trees, i.e. there cerorhynchus (Dinosauria: ) and relationships among basal neoceratopsians. J. Vertebr. Palaeontol. 18, is still disagreement among dinosaur systematicians rela- 569± 585. tive to most areas of the trees. This is more noticeable in Felsenstein, J. 1985 Phylogenies and the comparative method. the that have, on average, lower decay values Am. Nat. 125, 1± 15. than the Ornithischia. This might indicate that the phylo- Forster, C. T. 1999 Gondwanan dinosaur evolution and biog- genetic signal of ornithischians is less masked by homo- eographic analysis. J. Afr. Earth Sci. 28, 169± 185. plasy or may simply re¯ ect greater attention focused upon Gauthier, J. 1986 Saurischian monophyly and the origin of saurischian phylogeny. Never the less, we are con® dent . Mem. Calif. Acad. Sci. 8, 1± 55. with the resolution of con¯ icts of our supertree. Gordon, A. D. 1986 Consensus supertrees: the synthesis of In conclusion, we hope the dinosaur supertree rooted trees containing overlapping sets of labeled leaves. J. Classi® cation 3, 335± 348. presented will represent a solid framework for future study Graham, R. L. & Foulds, L. R. 1982 Unlikelihood that mini- of dinosaurian evolution, and will stimulate and direct mal phylogenies for a realistic biological study can be con- further systematic studies towards the less well understood structed in a reasonable computational time. Math. Biosci. areas of dinosaur systematics that are highlighted as poly- 60, 133± 142. tomies in the supertree. Holtz Jr, T. R. 1994 The phylogenetic position of Tyrannosau-

Proc. R. Soc. Lond. B (2002) A supertree of Dinosauria D. Pisani and others 921

ridae: implications for theropod systematics. J. Palaeontol. Sanderson, M. J., Purvis, A. & Henze, C. 1998 Phylogenetic 68, 1100± 1117. supertrees: assembling the tree of life. Trends Ecol. Evol. 13, Holtz Jr, T. R. 2000 A new phylogeny of carnivorous dino- 105± 109. saurs. Gaia 15, 5± 61. Sereno, P. C. 1999 The evolution of dinosaurs. Science 284, Hutchinson, J. R. 2001a The evolution of the pelvic osteology 2137± 2147. and soft tissues on the line to extant birds (Neornithes). Zool. Sereno, P. C. 2000 The fossil record, systematics and evolu- J. Linn. Soc. 131, 123± 168. tion of pachycephalosaurs and ceratopsians from . In Hutchinson, J. R. 2001b The evolution of the femoral oste- The age of dinosaurs in Russia and (ed. M. J. Benton, ology and soft tissues on the line to extant birds M. A. Shishkin, D. M. Unwin & E. N. Kurochkin), pp. 480± (Neornithes). Zool. J. Linn. Soc. 131, 169± 197. 516. Cambridge University Press. Ji, Q., Currie, P. J., Norrell, M. A. & Ji, S. A. 1998 Two feath- Sereno, P. C., Forster, C. A., Rogers, R. R. & Monetta, A. M. ered dinosaurs from northeastern . 393, 753± 1993 Primitive dinosaur from and the 761. early evolution of the Dinosauria. Nature 361, 64± 66. Kirkland, J. I. 1998 A polacanthine ankylosaur (Ornithischia: Sereno, P. C., Wilson, J. A., Larsson, H. C. E., Duthiel, Dinosauria) from the () of east- D. B. & Sues, H.-D. 1994 Early Cretaceous dinosaurs from ern . Mus. Nat. Hist. Sci. Bull. 14, 271± the Sahara. Science 265, 267± 271. 282. Swofford, D. L. 1998 Paup¤ . Phylogenetic analysis using parsi- Langer, M. C., Abdala, F., Richter, M. & Benton, M. J. 1999 mony (¤ and other methods). Version 4. Sunderland, MA: A sauropodomorph dinosaur from the Upper Triassic Sinauer. () of southern . Crit. Rev. Acad. Sci. II Thorley, J. L. & Page, R. D. M. 2000 RadCon: phylogenetic 329, 511± 517. tree comparison and consensus. Bioinformatics 16, 486± 487. Lanyon, S. M. 1993 Phylogenetic frameworks: towards a Thorley, J. L., Wilkinson, M. & Charleston, M. 1998 The ® rmer foundation for the comparative approach. Biol. J. information content of consensus trees. In Advances in data Linn. Soc. 49, 45± 61. science and classi® cation (ed. A. Rizzi, M. Vichi & H. H. Liu, F. G. R., Miyamoto, M. M., Freire, N. P., Ong, P. Q., Bock), pp. 91± 98. Berlin: Springer. Tennant, M. R., Young, T. S. & Gugel, K. F. 2001 Molecu- Upchurch, P. 1995 Evolutionary history of sauropod dino- lar and morphological supertrees for eutherian (placental) saurs. Phil. Trans. R. Soc. Lond. B 349, 365± 390. . Science 291, 1786± 1789. Upchurch, P. 1998 The phylogenetic relationships of sauropod MacLeod, N. 2001 The role of phylogeny in quantitative dinosaurs. Zool. J. Linn. Soc. 124, 43± 103. paleobiological data analysis. 27, 226± 240. Wilkinson, M. 1994 Common cladistic information and its Maddison, W. P. & Maddison, D. R. 1997 MacClade, v. 3.07. consensus representation: reduced Adams and reduced Sunderland, MA: Sinauer. cladistic consensus trees and pro® les. Syst. Biol. 43, 343± Mann, C. 1990 Meta-analysis in the breech. Science 249, 368. 476± 480. Wilkinson, M., Thorley, J. L., Littlewood, D. T. & Bray, R. A. Norman, D. B. 1998 On Asian ornithopods (Dinosauria: 2001 Towards a phylogenetic supertree for the Platyhel- Ornithischia). A new of iguanodontid dinosaur. Zool. minthes? In Interrelationships of the Platyhelminthes (ed. D. T. J. Linn. Soc. 122, 291± 348. Littlewood & R. A. Bray), pp. 292± 301. London: Chap- Pisani, D. & Wilkinson, M. 2002 Matrix representation with man & Hall. parsimony, taxonomic congruence and total evidence. Syst. Wilson, J. A. & Sereno, P. C. 1998 Early evolution and higher- Biol. 51, 151± 155. level phylogeny of sauropod dinosaurs. Soc. Vertebr. Purvis, A. 1995a A modi® cation to Baum and Ragan’s method Palaeontol. Mem. 5, 1± 68. for combining phylogenetic trees. Syst. Biol. 44, 251± 255. Witmer, L. M. 1995 The extant phylogenetic bracket and the Purvis, A. 1995b A composite estimate of primate phylogeny. importance of reconstructing soft tissues in . In Func- Phil. Trans. R. Soc. Lond. B 348, 405± 421. tional morphology in vertebrate (ed. J. J. Ragan, M. 1992 Phylogenetic inference based on matrix rep- Thomason), pp. 19± 33. Cambridge University Press. resentation of trees. Mol. Phylogenet. Evol. 1, 113± 126. Russell, D. A. & Dong, Z. M. 1993 The af® nities of a new theropod dinosaur from the Alxa Desert As this paper exceeds the maximum length normally permitted, the authors have agreed to contribute to production costs. People’s Republic of China. Can. J. Earth Sci. 30, 2107± 2127. Semple, C. & Steel, M. 2000 A supertree method for rooted Visit http://www.pubs.royalsoc.ac.uk to see an electronic appendix to trees. Discrete Appl. Math. 105, 147± 158. this paper.

Proc. R. Soc. Lond. B (2002)