- Research article
- Open Access
A supertree approach to shorebird phylogeny
BMC Evolutionary Biology volume 4, Article number: 28 (2004)
Order Charadriiformes (shorebirds) is an ideal model group in which to study a wide range of behavioural, ecological and macroevolutionary processes across species. However, comparative studies depend on phylogeny to control for the effects of shared evolutionary history. Although numerous hypotheses have been presented for subsets of the Charadriiformes none to date include all recognised species. Here we use the matrix representation with parsimony method to produce the first fully inclusive supertree of Charadriiformes. We also provide preliminary estimates of ages for all nodes in the tree.
Three main lineages are revealed: i) the plovers and allies; ii) the gulls and allies; and iii) the sandpipers and allies. The relative position of these clades is unresolved in the strict consensus tree but a 50% majority-rule consensus tree indicates that the sandpiper clade is sister group to the gulls and allies whilst the plover group is placed at the base of the tree. The overall topology is highly consistent with recent molecular hypotheses of shorebird phylogeny.
The supertree hypothesis presented herein is (to our knowledge) the only complete phylogenetic hypothesis of all extant shorebirds. Despite concerns over the robustness of supertrees (see Discussion), we believe that it provides a valuable framework for testing numerous evolutionary hypotheses relating to the diversity of behaviour, ecology and life-history of the Charadriiformes.
The shorebirds and allies (Aves: Charadriiformes; ) present an exceptional group for studying numerous evolutionary hypotheses. Their remarkable diversity of social mating system, parental care, sexual dimorphism, ecology and life-history make them an ideal group for unravelling the mechanisms of, for example, sexual selection and sexual conflict. Previous comparative studies have made significant contributions to our understanding of the evolution of mating systems , parental care [3, 4], sexual size dimorphism [5–7], locomotion and morphology , migratory behaviour , egg size , and plumage colouration . The importance of phylogeny in cross-species comparative studies is well documented [12–14]. Large and well-resolved phylogenies that incorporate divergence times provide powerful tests of a wide range of hypotheses whilst accounting for the effects of shared evolutionary history [13, 15]. However, the shorebird studies listed above were limited by the lack of a complete phylogeny for the group. Most of these studies are based on derivations of the seminal work of Sibley and Ahlquist , yet this study included less than a quarter of extant and recently extinct shorebird species. Recently extinct taxa (according to Monroe and Sibley ) are: the Tahitian sandpiper Prosobonia leucoptera, the Canary Islands oystercatcher Haematopus maedewaldoi, and the Great auk Pinguinus impennis.
Recent molecular studies covering a wide range of shorebird families have drawn attention to conflict in the reconstruction of the deep basal nodes of shorebird phylogeny (figure 1; reviewed by van Tuinen et al. ). For example, morphological data [18, 19] places Alcinae (auks, puffins, murres) at the base of the shorebird tree whilst sequence [20–22] and DNA-DNA hybridisation  data suggests that they are a highly derived sister group to Stercorariini (skuas and jaegers), Larini (gulls), Sternini (terns), and Rynchopini (skimmers). It is important to note that taxon coverage differs between these studies and this may be an important factor in determining the tree topology. Specific phylogenies have been derived, for example, for sandpipers , the genus Charadrius , and jacanas  using DNA sequence data. In contrast, morphological evidence provided the basis for Chu's  study of gull phylogeny. Strauch  presented the most complete data set of 227 Charadriiformes species. However, despite the plethora of cladograms for particular shorebird groups (see reviews by Sibley and Ahlquist ; Thomas et al. ), those that address relationships across the whole clade use either sparse taxon sampling [16, 27], or are based on reassessments of Strauch's  data [19, 28–30]. Note that Dove  included a feather microstructural analysis in addition to her reanalysis of Strauch's  data.
Combining phylogenetic data
Numerous methods and types of data can be used to infer phylogeny. Frequently, as in Charadriiformes, a single analysis incorporating all taxa of interest is absent. Under the principle of total evidence , all sources of phylogenetic information should be combined to maximize their explanatory power. Eernisse and Kluge  define total evidence as a method for seeking the best fitting phylogenetic hypothesis for an unpartitioned set of synapomorphies (shared derived characters) using character congruence (characters combined in a supermatrix). Hence, this method combines the primary data (molecular, morphological and behavioural characters) into a single analysis. The approach is powerful because weak signals in the partitioned data sets may be enhanced when combined, and previously obscured relationships may be revealed .
The total evidence approach has both practical and theoretical problems. First, only certain types of data can be combined. For example, nucleotide sequences and morphological traits can be readily assessed together as characters, but it is not generally possible to include nucleotide sequences and genetic distance data in a single analysis . We acknowledge that Lapointe et al.  suggest a distance based approach to combine otherwise incompatible data in a total evidence analysis, although this method has not been tested beyond a single application. The consequence is that it is rarely possible to combine all sources of data in practice and the lack of overlap in combinable data sets may result in a reduction of the number of taxa included. Second, Miyamoto and Fitch  contend that combining data sets is rarely justified because partitions of phylogenetic data are real and unequivocal. They argue that several partitions producing similar topologies provide multiple lines of independent evidence supporting that topology.
Theoretical arguments over the benefits of total evidence will undoubtedly continue, but perhaps the major barriers to its use are the often very high computational demands of large matrices, and the a priori exclusion of certain data types. This is particularly true of Charadriiformes phylogeny, where one of the most significant contributions to the field – DNA-DNA hybridisation – cannot be included. An alternative set of techniques, collectively termed supertrees (e.g., Matrix Representation with Parsimony, MRP; [37, 38]), enables combination of trees (rather than raw data) from otherwise incompatible sources. MRP methods code source phylogenies based on the presence and absence of taxa at each node of the tree [37–39] and are thus one step removed from the primary data. It is important to recognise that supertrees should not be regarded as a replacement for exhaustive phylogenetic studies of the primary data and there are drawbacks to the methods (see Discussion). However, they do enable very large phylogenies to be constructed rapidly . Supertrees have been constructed successfully for a wide variety of taxa including carnivores , primates , seabirds , dinosaurs , and grasses .
Shorebirds are particularly well suited for supertree treatment, since there are numerous incomplete phylogenies available and a broader phylogeny is desirable to facilitate powerful analyses of numerous evolutionary hypotheses (see above). Here, we present the first complete composite phylogeny of extant and recently extinct  shorebirds using the MRP approach. We are therefore combining data on tree topologies, and not conducting a simultaneous analysis on the original data. We also use fossil and molecular data to estimate divergence times (see Methods). The combination of complete taxonomic coverage and the inclusion of branch lengths provide the basis for future comparative analyses of Charadriiformes evolution. In addition, conflicting and unresolved areas of Charadriiformes phylogeny are revealed.
Results and Discussion
Supertree resolution and topology
We found 1469 equally short trees of length 1847 steps using the parsimony ratchet approach (see Methods). This compares favourably to a standard heuristic search that yielded shortest trees of 1853 steps. All subsequent results and discussion refer to the parsimony ratchet analyses. Figure 2 shows the family and subfamily level relationships of shorebirds based on the strict and 50 % majority-rule consensus tree (see additional file 1 for branch length estimates). Figures 3,4,5,6,7,8,9 show the species level phylogeny. The full 50% majority rule consensus and the strict consensus trees are available as additional file 2 and 3 respectively. The 50% majority-rule consensus tree is well resolved (73.1%; 255 nodes out of a possible 349 in a fully bifurcating tree), although the strict consensus tree is only 49.6% resolved (173 from 349 possible nodes). The majority rule tree includes nine novel clades (numbers 20, 29, 57, 85, 89, 108, 122, 139, 140) that do not appear in any of the source trees; all of these occur towards the tips of the tree. This is a general problem in supertree construction and such clades should be collapsed as they have no support . To demonstrate where the MRP method has performed badly we have included the novel clades in all figures and list details in the figure legends. In addition, 58 nodes are supported by only one character (see additional file 1). Each of these nodes is left over from a single source tree. Assessing the support for such nodes is problematic because this may simply reflect a lack of research directed at the taxa in question. A major challenge for supertree construction is to develop measures of support that reflect the robustness of nodes in the source trees. We list the number of characters supporting each node (additional file 1) but stress that these are not measures of tree robustness and may not be directly comparable even within the same tree. This is because the taxon coverage across source trees is highly variable so some nodes have more potential support than others. Furthermore, because measures of support used in the source trees differ between studies (some source trees include no measures of support), it is impractical and of dubious value to use these measures to assess the robustness of the supertree.
The majority of unresolved nodes in the shorebird supertree are located towards the tips of the phylogeny. For example, the genus Gallinago forms a monophyletic clade but only two pairs of species are resolved from 14 species (G. megala and G. negripennis; G. macrodactyla and G. media) in the majority-rule tree. Only the latter relationship remains in the strict consensus tree. In addition, clades including the genera Charadrius and Vanellus, Calidris and Tringa, Sterna, and Scolopax are poorly resolved. This may reflect a bias in phylogenetic studies of shorebirds. For instance, we found six source trees for Alcinae [43–48] but none devoted to Scolopax or Gallinago. Thomas et al.  indicate that this may be a problem for shorebird studies in general and reported a strong skew favouring research on northern hemisphere species.
In contrast to the within genera relationships, the generic and family levels are generally well resolved. The supertree indicates three monophyletic Charadriiformes lineages (figure 2). Family and subfamily resolution within each lineage is high, however the relative position of each group is unresolved in the strict consensus tree. This is an important point because the deepest relationships of shorebird phylogeny are contentious . The 50% majority-rule consensus tree indicates that the gulls and allies (Larini, Sternini, Rynchopini, Stercorariini, Dromas, Alcinae, and Glareolidae) are sister to the sandpipers and allies (Scolopacidae, Jacanidae, Rostratulidae, Thinocoridae, Pedionomidae). The most basal lineage includes the plovers and allies (Charadriinae, Pluvianellidae, Chionidae, Burhinidae, Haematopodini and Recurvirostrini). The gulls and allies clade is most consistent with DNA-DNA hybridisation , indicating that Larini are sister to Sternini and that Rynchopini are sister to this group. This conflicts with morphology-based topologies where Stercorariini are sister to Larini and Sternini with Rynchopini basal to both. Indeed, the position of Stercorariini remains controversial and most recently they were placed as sister to Alcinae [20–22]. In contrast, morphological evidence [18, 19] places Alcinae at the base of the whole Charadriiformes tree with Stercorariini sister to Larini. Thus, the position of Alcinae is uncertain and appears to be dependent on the type of data, with fundamental differences between molecular based analyses and morphological analyses. The taxon sampling of previous morphological and molecular studies varies considerably and it may be this, rather than genuine differences in the phylogenetic signal of different data types, that is the cause of conflict in resolving the phylogenetic position of Alcinae. However, it is encouraging that van Tuinen et al.  suggested that new unpublished osteological data are consistent with the more derived position indicated by molecular data. The supertree resolves Glareolidae outside the Larini, Sternini, Rynchopini, Stercorariini, Dromas, Alcinae clade. This is also the case with recent molecular and previous DNA-DNA hybridisation studies. Morphological studies have failed to resolve the position of Glareolidae, placing the family in a large polytomy with all other major groups except Alcinae and the sandpipers and allies (fig. 1). A novel development in shorebird phylogeny is the placement of the black-rumped buttonquail Turnix hottentotta as a sister to the gulls and allies (Larini, Sternini, Rynchopini, Stercorariini, Dromas, Alcinae, and Glareolidae) based on the nuclear RAG-1 gene . We did not include this species in the supertree because to date Paton et al.  remains the only study to reveal an apparently robust relationship. More diverse sampling of the buttonquails (Turnicidae) is essential to corroborate the general affinities of this family.
The relationships within the plover clade appear to be reasonably stable. Morphological, molecular, and DNA-DNA hybridisation all place Charadriinae as sister to Haematopodini and Recurvirostrini; our supertree is consistent with these relationships. However, it is not clear whether Burhinidae and Chionidae are sister to each other [20–22] or whether Chionidae are sister to a Charadriinae, Haematopodini, Recurvirostrini, and Burhinidae clade . Our supertree also included Pluvianellidae, a family consisting of only one species (magellanic plover Pluvianellus socialis) and places this as sister to Chionidae. If Pluvianellidae are excluded, the supertree is consistent with the sister group relationship of Burhinidae and Chionidae.
The sister group relationship of Jacanidae to Rostratulidae is well established [16, 18–22] and is found in our supertree. The supertree resolves the Thinocoridae and Pedionomidae as sister taxa and this group is sister to the Jacanidae and Rostratulidae. The large Scolopacidae clade is at the base of the sandpiper clade consistent with recent molecular studies [20–22] and the DNA-DNA hybridisation tapestry .
Taken together, it is evident that the supertree is generally more consistent with molecular data (both recent sequence studies and DNA-DNA hybridisation) than with analyses based on morphology. However, it is of course possible that this reflects the greater number of molecular source trees available rather than indicating that molecular data is actually better at resolving shorebird phylogeny. We included several large morphological phylogenies [e.g [18, 19, 26, 30, 43]] but the majority of source trees (29 out of 51) were based on molecular evidence (see additional file 5).
The higher resolution of the majority-rule tree means it is more likely to be of use in comparative studies. We therefore estimated node ages for this topology only (see additional file 1 and 2). We stress that our estimates of node dates are a first attempt at dating the whole tree and have several limitations. First, the fossils used to calibrate seven nodes in the tree are unlikely to be the earliest members of their respective families thus these dates will be underestimates. Second, we assumed that the fossils are grouped with the extant members of the family but this requires formal testing in a phylogenetic framework. Third, the pure birth model assumes that no extinction occurs but this may be unrealistic and it is likely that extinction processes have reduced the representation of older lineages . Furthermore, this model is derived from the topological structure of the tree so errors in tree reconstruction will likely lead to errors in branch length estimation. However, this approach has been employed previously in supertrees of primates  and carnivores  explicitly to facilitate comparative analyses. Despite these caveats, simulation studies have demonstrated that comparative methods such as independent contrasts are robust to errors in branch length  and no viable alternative for dating supertrees has been proposed. Nonetheless, we urge that alternative branch length assumptions are explored if the shorebird supertree is used in future comparative studies. At present, the calibrated RAG-1 tree of Paton et al.  remains arguably the most thorough and reliable measure of divergence times for Charadriiformes.
A fuller understanding of the phylogenetic affinities of fossil shorebirds will probably improve estimates of node ages for the group. For example, the extinct form Graculavidae, is represented by fossils from the Maastrichtian of New Jersey  and Cretaceous of Wyoming  but its position within the shorebird clade is unclear. Feduccia  suggests that it may be basal and a formal corroboration of this would support proposals for a late Cretaceous origin of shorebirds. The difficulties in dating the shorebird tree are further illustrated by fossil representatives of Recurvirostrini and Burhinidae which are much older than current estimates suggests. The earliest record of the Recurvirostrini is estimated to be over 50 million years old  whilst recent discoveries of a possible member of the Burhinidae are dated to around 70 mya [55, 56]. There is clearly a need for an integrated phylogenetic study including both extinct and extant shorebirds.
Supertrees are still at an early stage of development and many aspects of MRP, and supertree methods in general, are not yet clearly understood. Steps can be taken to ensure that the supertree includes the most appropriate sets of sources trees, such as only using trees from explicitly phylogenetic studies. This is not always straightforward and could result in the exclusion of important information. For instance, in our shorebird supertree, we included Sibley and Ahlquist's DNA-DNA hybridisation tapestry  although this is based on distance measures rather than more rigorous phylogenetic methods. Even if very strict tree selection criteria are applied, there are still likely to be biases in the data set. For example, not all source trees are equally well supported, yet in most supertree analyses each tree is treated equally . This is a problem for supertree construction because whilst it is theoretically possible, and indeed beneficial, to weight source trees based on support values  it is rarely possible in practice. Many source trees do not have support values and those that do may use different methods, (e.g, bootstrapping or decay indices) which cannot be directly compared with each other. An additional problem that has not been fully resolved relates to correlations between source trees . Several source trees based on the same data set may unduly increase the influence of that data set on the supertree analysis. However, there is no formal way of determining how much overlap to allow and the choice of source trees that go into supertree construction inevitably involves some degree of subjective reasoning. For the shorebird supertree we used strict Reduced Cladistic Consensus trees to summarise potential source trees that were from the same data set but based on different methods. For example, Thomas et al.  based their phylogeny on cytochrome-b but used a range of methods including parsimony and Bayesian analyses. We therefore combined these trees to minimise bias. In contrast, Ericson et al.  used two types of data: sequences from the nuclear RAG 1 gene and sequences from the myoglobin intron II. They carried out three analyses: each gene separately and then the two combined in a single analysis. In this case, we used three source trees. It could be argued that the combined analysis of Ericson et al.  should be excluded because of the possible overlap with the individual analyses. However, under the principle of total evidence, the combined data set may result in novel relationships being revealed [31, 33] and therefore could contribute important information to the supertree. Simulation and empirical studies are required to fully understand these and other possible biases in supertree construction (e.g., the influence of source tree size and shape) and formal protocols for the selection of source trees are desirable. For transparency, we include a summary of the source trees used, data type, and the main taxa included in the study (additional file 5).
Our shorebird supertree is highly consistent with recent advances in the molecular phylogenetics Charadriiformes. However, we urge caution when using the tree in comparative analyses and encourage the additional use of alternative phylogenies and branch length assumptions. It is particularly important to note that the position of some groups such as the Alcinae remains controversial and that although the majority rule tree is consistent with recent molecular studies, the strict consensus tree fails to resolve the deepest nodes.
The supertree presented here is, to our knowledge, the first attempt to reconstruct the phylogeny of the entire order Charadriiformes. Overall, the supertree is highly consistent with recent molecular hypotheses of shorebird phylogeny. However, it is apparent that fresh attempts to resolve both the phylogeny and estimates of age will be dependent on further gene sequencing and new fossil discoveries. The affinities of the Alcinae and the relationships between the three major shorebird clades require further corroboration, and studies of several genera such as Gallinago and Vanellus are desirable. Furthermore, additional work is required to establish the true affinities of the Turnicidae. Nonetheless, it appears that shorebird phylogeny is gradually approaching a consensus view. The broad taxonomic scope and consistency of the supertree mean that is of potentially great value to future comparative studies (accepting the caveats discussed above) of the behaviour, life-history, ecology and conservation of this diverse group.
Possible source trees were identified from online searches of Web of Science http://wos.mimas.ac.uk covering the years 1981 to 2004. We used the single key strings phylogen*, cladistic*, clado*, classif*, systematic*, and taxonom* (where the asterisks allow variations such as "phylogeny" or "phylogenetics") in the topic field, in conjunction with a major Charadriiformes taxon name (scientific or common). As supertree methods have been criticized for being biased towards historical trends, we preferred those studies that explicitly set out to derive a phylogenetic hypothesis and so exclude purely (and typically older) descriptive taxonomic works. The Sibley and Ahlquist  DNA-DNA hybridisation tapestry may be viewed as non-cladistic, but it was clearly the authors' intention to reconstruct the phylogeny of birds. Furthermore, it provided a vital catalyst for subsequent studies of avian (including shorebird) phylogeny. We therefore included the DNA-DNA hybridisation hypothesis as a source tree in our analyses. Simulation studies have demonstrated that the performance of supertree methods is improved by including at least one taxonomically complete (or near complete) source tree . We therefore make an exception to our self-imposed rule, and in addition use the taxonomic hierarchy of Monroe and Sibley  as a source tree as this includes all extant Charadriiformes species. We acknowledge that this taxonomy is based largely on Sibley and Ahlquist's  DNA-DNA hybridisation tapestry.
The initial search identified 78 source trees from 44 publications. Each source tree was typed as a text file in Nexus format . We coded trees to the species level with species names taken from Monroe and Sibley , but note that contra Monroe and Sibley , we use Charadriiformes not Charadrii to refer to the whole group. Several studies included the gull Larus thayeri [26, 60–63] either as a subspecies of Larus glaucoides (Larus glaucoides thayeri in Monroe and Sibley ) or a species in its own right. In recognition of this, we included Larus glaucoides thayeri as the only subspecies in our data set thus increasing the total taxa to 366. Monroe and Sibley  include 16 species of the family Pteroclidae within the Charadriiformes. However, the relationship of this family to the Charadriiformes is uncertain and they have recently been placed in their own order . We include the Pteroclidae in our analyses only as a means of rooting the tree. Where there were multiple most parsimonious trees (MPTs), or where source trees had been derived from predominantly overlapping data (e.g., from the same data but using alternative methods), we used RadCon  to produce strict Reduced Cladistic Consensus trees (RCC [66, 67]). The output is in the form of a reduced consensus profile and from this we selected the tree with the highest Cladistic Information Content (CIC) [65, 68]. This resulted in a total of 51 source trees from which our supertree is derived and these are summarised in additional file 5.
We produced an MRP matrix of the 51 Nexus  source trees in RadCon  (see additional file 6 for the MRP file). We used the original MRP coding method of Baum  and Ragan . Weighting source trees based on node support such as bootstrapping improves the accuracy of MRP supertrees . However, this is only possible if all source trees can be weighted on the same criteria . The absence of branch support measures in many of the shorebird source trees precludes this approach from the present study; hence, subsequent analyses were conducted using equally weighted parsimony.
The tendency of large data sets to produce many sub-optimal trees that are close in length and topology to the shortest tree is a serious problem in phylogenetics. Standard heuristic searches frequently are trapped searching within globally sub-optimal "islands" and the tree search is often aborted before completion. Nixon  proposed a new method to avoid this problem. The "Parsimony Ratchet" reweights a random set of characters from the data set. This may result in the tree island no longer representing a local optimum and the heuristic search continues until a new optimum is reached. The algorithm then reverts to the original weighting and the search continues. Nixon  demonstrated the efficacy of the method on a 500-taxon data set, where the ratchet-based search found a tree two steps shorter than standard heuristic searches. We used PAUPRat  to implement a parsimony ratchet in PAUP* . The default settings of 200 iterations and 15% perturbation of characters for reweighting were used and we carried out 20 replicates. Equally parsimonious trees were summarized using both strict and 50% majority-rule consensus methods.
We did not calculate any measures of branch support for two reasons. First, their validity and meaning is questionable in MRP supertrees . Second, the number of taxa included in our data set is too large to allow practical calculation of any branch support indices (e.g., decay indices ) on a desktop computer.
Dating the supertree
Following Purvis  and Bininda-Emonds et al.  we dated the supertree using both absolute and relative dates. We used data from the Fossil Record 2  as the source of fossil-based absolute dates. This yielded estimates for Jacanidae (Nupharanassa tolutaria, Rupellian), Phalaropus (Phalaropus elenorae, Middle Pliocene), Burhinidae (Burhinus lucorum, Lower Miocene), Glareolidae (Paractiornis perpusillus, Lower Miocene), Alcinae (Petralca austrica, Rupellian), Stercoariini (Stercorarius sp., Middle Miocene), and Larini (undetermined, Rupellian). We took the midpoint of the range from the Fossil Record 2  as our date estimate. More recent publications of fossil Charadriiformes were not included because they either represent specimens that are younger or have not been assigned to families that are represented amongst the extant Charadriiformes (such as Turnipacidae ). We assumed that fossil dates represent the earliest occurrence for each group which inevitably resulted in underestimates of clade age. The fossil record of Charadriiformes is amongst the best of the modern bird groups  in terms of the numbers of taxa, but many specimens are fragmentary and reliable estimates of divergence dates are dependent on a limited number of exceptional specimens . The phylogenetic affinities of the fossil shorebirds in relation to their extant relatives have not yet been fully established, hence have implicitly assumed that fossil representatives of extant groups would be resolved amongst their living relatives.
Source trees may include estimates of relative branch lengths (e.g., genetic distances). This allows further dating of the supertree but is problematic because different relative estimates are not comparable and cannot be applied directly to the supertree . However, where a source trees shares a node that has an absolute date in the supertree (a node dated from fossil evidence), the relative branch lengths can easily be converted to estimates of age. All taxa in our supertree are either extant, or very recently extinct; hence, the tips of the calibrated supertree should be equidistant from the root of the tree. In source trees where the relative branch lengths are not equidistant from the root, we followed the protocol of Purvis [; p.407–8]. We estimated relative dates using the local molecular clock logic  as implemented by Purvis  and Bininda-Emonds et al. . For example, consider three taxa A, B, and C where A and B are sister taxa and C is sister to A and B. The root is dated to 10 million years (myr) from fossil evidence, and independent molecular data provides estimates of divergence based on the number of substitutions per site. The molecular estimates of branch lengths are as follows: A, 6 substitutions; B, 8 substitutions; C, 20 substitutions; A and B are 11 substitutions from the root. A and B are therefore separated from their common node by a mean of 7 substitutions. The total length from A and B to the root is thus 18 substitutions compared to 20 for C (a mean of 19). This can be converted to date estimates such that 19 substitutions are equivalent to 10 myr. The dates of the tree are then: ((A: 3.68, B: 3.68), C: 10)). There were no cases where multiple source trees with molecular divergence dates were able to provide estimates for the same node. We estimated relative dates from multiple nodes rather than a single dated node to minimise correlative errors in estimates.
To provide date estimates for all nodes in the tree we employed a pure birth model to date nodes for which absolute and relative dates could not be attained . Pure birth models infer that a clade's age is proportional to the logarithm of the number of species within the clade:
date of daughter = date of ancestor *(log daughter clade size/log parent clade size)
For example, the age of a daughter node that subtends 12 taxa, estimated from its immediate ancestor dated to 20 myr and which subtends 19 taxa is:
20*(log(12)/log(19)) = 16.879
We applied this approach to estimate the ages of daughter nodes based on dates (absolute or calibrated) of ancestral nodes. We had no ancestral node on which to base estimates of the most basal clade. In this case, we rearranged the pure birth formula and calculated the age of the ancestral node from its two daughter nodes, taking the mean as our "best estimate". Finally, to estimate the ages of nodes between daughter and ancestor nodes of known age we spaced the nodes evenly along the branches length .
Monroe BL, Sibley CG: A World Checklist of Birds. 1993, New Haven; Yale University Press
Whitfield DP, Tomkovich PS: Mating system and timing of breeding in Holarctic waders. Biol J Linn Soc. 1996, 57: 277-290. 10.1006/bijl.1996.0015.
Reynolds JD, Székely T: The evolution of parental care in shorebirds: life histories, ecology, and sexual selection. Behav Ecol. 1997, 8: 126-134.
Borowik OA, Mclennan DA: Phylogenetic patterns of parental care in Calidrine sandpipers. Auk. 1999, 116: 1107-1117.
Figuerola J: A comparative study on the evolution of reversed size dimorphism in monogamous waders. Biol J Linn Soc. 1999, 67: 1-18. 10.1006/bijl.1998.0254.
Székely T, Reynolds JD, Figuerola J: Sexual size dimorphism in shorebirds, gulls, and alcids: the influence of sexual and natural selection. Evolution. 2000, 54: 1404-1413.
Lindenfors P, Székely T, Reynolds JD: Directional changes in sexual size dimorphism in shorebirds, gulls and alcids. J Evol Biol. 2003, 16: 1-9. 10.1046/j.1420-9101.2003.00506.x.
Barbosa A, Moreno E: Hindlimb morphology and locomotor performance in waders: an evolutionary approach. Biol J Linn Soc. 1999, 67: 313-330. 10.1006/bijl.1998.0282.
Myers JP: Cross-seasonal interactions in the evolution of sandpiper social systems. Behav Ecol Sociobiol. 1981, 8: 195-202.
Liker A, Reynolds JD, Székely T: The evolution of egg size in socially polyandrous shorebirds. Oikos. 2001, 95: 3-14.
Bókony V, Liker A, Székely T, Kis J: Melanin-based plumage colouration and flight displays in plovers and allies. Proc Roy Soc Lond B. 2003, 270: 2491-2497. 10.1098/rspb.2003.2506.
Felsenstein J: Phylogenies and the comparative method. Am Nat. 1985, 125: 1-15. 10.1086/284325.
Harvey PH, Pagel MD: The Comparative Method in Evolutionary Biology. 1991, Oxford; Oxford University Press
Martins EP: Adaptation and the comparative method. TREE. 2000, 15: 295-299.
Bininda-Emonds ORP, Gittleman JL, Purvis A: Building large trees by combining phylogenetic information: a complete phylogeny of the extant Carnivora (Mammalia). Biol Rev. 1999, 74: 143-175. 10.1017/S0006323199005307.
Sibley CG, Ahlquist JE: Phylogeny and Classification of Birds. 1990, New Haven; Yale University Press
van Tuinen M, Waterhouse D, Dyke GJ: Avian molecular systematics on the rebound: a fresh look at modern shorebird phylogenetic relationships. J Avian Biol. 2004, 35: 191-194. 10.1111/j.0908-8857.2004.03362.x.
Strauch JG: The phylogeny of the Charadriiformes (Aves): a new estimate using the method of character compatibility analysis. Trans Zool Soc Lond. 1978, 34: 263-345.
Chu PC: Phylogenetic reanalysis of Strauch's osteological data set for the Charadriiformes. Condor. 1995, 97: 174-196.
Paton TA, Baker AJ, Groth JG, Barrowclough GF: RAG-1 sequences resolve phylogenetic relationships within Charadriiform birds. Mol Phyl Evol. 2003, 29: 268-278. 10.1016/S1055-7903(03)00098-8.
Ericson PGP, Envall I, Irestedt M, Norman JA: Inter-familial relationships of the shorebirds (Aves: Charadriiformes) based on nuclear DNA sequence data. BMC Evol Biol. 2003, 3: 16-10.1186/1471-2148-3-16.
Thomas GH, Wills MA, Székely T: Phylogeny of shorebirds, gulls, and alcids (Aves: Charadrii) from the cytochrome-b gene: parsimony, Bayesian inference, minimum evolution and quartet puzzling. Mol Phyl Evol. 2004, 30: 516-526. 10.1016/S1055-7903(03)00222-7.
Christidis L, Davies K, Westerman M, Christian PD, Schodde R: Molecular assessment of the taxonomic status of Cox's sandpiper. Condor. 1996, 98: 459-463.
Joseph L, Lessa EP, Christidis L: Phylogeny and biogeography in the evolution of migration: shorebirds of the Charadrius complex. J Biogeogr. 1999, 26: 329-342. 10.1046/j.1365-2699.1999.00269.x.
Whittingham LA, Sheldon FH, Emlen ST: Molecular phylogeny of jacanas and its implications for morphologic and biogeographic evolution. Auk. 2000, 117: 22-32.
Chu PC: A phylogeny of the gulls (Aves: Larinae) inferred from osteological and integumentary characters. Cladistics. 1998, 14: 1-43. 10.1006/clad.1997.0051.
Christian PD, Christidis L, Schodde R: Biochemical systematics of the Charadriiformes (shorebirds): relationships between the Charadrii, Scolopaci, and Lari. Aus J Zool. 1992, 40: 291-302.
Mickevich MF, Parenti LR: [Review of] The phylogeny of the Charadriiformes (Aves): a new estimate using the method of character compatibility analysis. Syst Biol. 1980, 29: 108-113.
Björklund M: Phylogenetic relationships among Charadriiformes: reanalysis of previous data. Auk. 1994, 111: 825-832.
Dove CJ: A descriptive and phylogenetic analysis of plumalaceous feather characters in Charadriiformes. Ornith Monogr. 2000, 51: 1-163.
Kluge AG: A concern for evidence and a phylogenetic hypothesis of relationships among Epicrates (Boidae, Serpentes). Syst Zool. 1989, 38: 7-25.
Eernisse DJ, Kluge AG: Taxonomic congruence versus total evidence, and amniote phylogeny inferred from fossils, molecules, and morphology. Mol Biol Evol. 1993, 10: 1170-1195.
Barrett M, Donoghue MJ, Sober E: Against consensus. Syst Zool. 1991, 40: 486-493.
de Queiroz A, Donoghue MJ, Kim J: Separate versus combined analysis of phylogenetic evidence. Ann Rev Ecol Syst. 1995, 26: 657-681. 10.1146/annurev.es.26.110195.003301.
Lapointe F-J, Kirsch JAW, Hutcheon JM: Total evidence, consensus, and bat phylogeny: a distance-based approach. Mol Phyl Evol. 1999, 11: 55-66. 10.1006/mpev.1998.0561.
Miyamoto MM, Fitch WM: Testing species phylogenies and phylogenetic methods with congruence. Syst Biol. 1995, 44: 64-76.
Baum BR: Combining trees as a way of combining data sets for phylogenetic inference, and the desirability of combining gene trees. Taxon. 1992, 42: 637-640.
Ragan MA: Phylogenetic inference based on matrix representation of trees. Mol Phyl Evol. 1992, 8: 53-58.
Purvis A: A composite estimate of primate phylogeny. Phil Trans Roy Soc Lond B. 1995, 348: 405-421.
Kennedy M, Page RDM: Seabird supertrees: combing partial estimates of procellariform phylogeny. Auk. 2002, 119: 88-108.
Pisani D, Yates AM, Langer MC, Benton MJ: A genus-level supertree of the Dinosauria. Proc Roy Soc Lond B. 2002, 269: 915-921. 10.1098/rspb.2001.1942.
Salamin N, Hodkinson TR, Savolainen V: Building supertrees: an empirical assessment using the grass family (Poaceae). Syst Biol. 2002, 51: 136-150. 10.1080/106351502753475916.
Strauch JG: The phylogeny of the Alcidae. Auk. 1985, 102: 520-539.
Friesen VL, Baker AJ, Piatt JF: Phylogenetic relationships within the Alcidae (Charadriiformes: Aves) inferred from total molecular evidence. Mol Biol Evol. 1996, 13: 359-367.
Friesen VL, Montevecchi WA, Davidson WS: Cytochrome-b nucleotide-sequence variation among the atlantic Alcidae. Hereditas. 1993, 119: 245-252. 10.1111/j.1601-5223.1993.00245.x.
Moum T, Arnason U, Arnason E: Mitochondrial DNA sequence evolution and phylogeny of the Atlantic Alcidae, including the extinct great auk (Pinguinus impennis). Mol Biol Evol. 2002, 19: 1434-1439.
Moum T, Johansen S, Erikstad KE, Piatt JF: Phylogeny and evolution of the auks (subfamily Alcinae) based on mitochondrial DNA sequences. PNAS. 1994, 91: 7912-7916.
Watada M, Kakizawa R, Kuroda N, Utida S: Genetic differentiation and phylogenetic relationships of an avian family, Alcidae (auks). J Yamashina Inst Ornith. 1987, 19: 79-88.
Thomas GH, Székely T, Sutherland WJ: Publication bias in waders. Wader Study Group Bull. 2003, 100: 216-223.
Diaz-Uriarte R, Garland T: Effects of branch length errors on the performance of phylogenetically independent contrasts. Syst Biol. 1998, 47: 654-672. 10.1080/106351598260653.
Olsen SL, Parris DC: The Cretaceous birds of New Jersey. Smithsonian Cont Paleobiol. 1987, 63: 1-22.
Hope S: A new species of Graculavus from the Cretaceous of Wyoming (Aves: Neornithes). Smithsonian Cont Paleobiol. 1999, 89: 261-266.
Feduccia A: The Origin and Evolution of Birds. 1996, New Haven: Yale Univeristy Press
Unwin DM: Aves. In The Fossil Record 2. Edited by: Benton MJ. 1990, London: Chapman and Hall, London, 717-737.
Case JA: Latest Cretaceous record of modern birds from Antarctica: center of origin or fortuitous occurrence?. PaleoBios. 2001, 21: 40A-
Cordes AH: A new Charadriiform avian specimen from the early Maastrichtian of Cape Lamb, Vega Island, Antarctic Peninsula. J Vertebr Paleontol. 2002, 22: 99A-
Bininda-Emonds ORP, Sanderson MJ: Assessment of the accuracy of matrix representation with parsimony analysis supertree construction. Syst Biol. 2001, 50: 565-579. 10.1080/106351501750435112.
Springer MS, de Jong WW: Phylogenetics – Which mammalian supertree to bark up?. 2001, 291: 1709-10.1126/science.1059434.
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony and Other Methods (software). 2000, Sunderland: Sinauer Associates
Crochet P-A, Desmaris E: Slow rate of evolution in the mitochondrial control region of gulls (Aves: Laridae). Mol Biol Evol. 2000, 17: 1797-1806.
Crochet P-A, Lebreton J-D, Bonhomme F: Systematics of large white-headed gulls: patterns of mitochondrial DNA variation in Western European taxa. Auk. 2002, 119: 603-620.
Hoffman W: Phylogeny, feeding behaviour, and wing structure in gulls, terns, and allies (Laroidea). PhD thesis. 1984, University of South Florida
Schnell GD: A phenetic study of the suborder Lari (Aves) II. Phenograms, discussion, and conclusions. Syst Zool. 1970, 19: 264-302.
Perrins C: The New Encyclopedia of Birds. 2003, Oxford: Oxford University Press
Thorley JL, Page RDM: RADCON: phylogenetic tree comparison and consensus. Bioinformatics. 2000, 16: 486-487. 10.1093/bioinformatics/16.5.486.
Wilkinson M: Common cladistic information and its consensus representation: reduced Adams and reduced cladistic consensus trees and profiles. Syst Biol. 1994, 43: 343-368.
Wilkinson M: Missing entries and multiple trees: Instability, relationships, and support in parsimony analysis. J Vertebr Paleont. 2003, 23: 311-323.
Thorley JL, Wilkinson M, Charleston M: The information content of consensus trees. In Advances in Data Science and Classification. Edited by: Rizzi A, Vichi M, Bock HH. 1998, Berlin: Springer, 91-98.
Nixon KC: The parsimony ratchet, a new method for rapid parsimony analysis. Cladistics. 1999, 15: 407-414. 10.1006/clad.1999.0121.
Sikes DS, Lewis PO: beta software, version 1. PAUPRat: PAUP* implementation of the parsimony ratchet. Distributed by the authors. Department of Ecology and Evolutionary Biology, University of Connecticut, Storrs, USA;. 2001
Bremer K: The limits of amino acid sequence data in angiosperm phylogenetic reconstruction. Evolution. 1988, 9: 427-436.
Mayr G: Charadriiform birds from the early Oligocene of Céreste (France) and the Middle Eocene of Messel (Hessn, Germany). Geobios. 2000, 33: 625-636.
Olsen SL: The fossil record of birds. In Avian Biology. Edited by: Farner DS, King JR, Parkes KC. 1985, Orlando: Academic Press, 8: 80-238.
Bailey WJ, Fitch DHA, Tagle DA, Czelusniak J, Slightom JA, Goodman M: Molecular evolution of the _-globin gene locus: gibbon phylogeny and the hominid slowdown. Mol Biol Evol. 1991, 8: 155-184.
Losos JB: Ecomorphology, performance capability and scaling of West Indian Anolis lizards: an evolutionary analysis. Ecol Monogr. 1990, 60: 369-388.
Baker AJ, Strauch JG: Genetic variation and differentiation in shorebirds. In Proceedings of the 19th International Ornithological Congress. 1988, 1639-1645.
Baker AJ: Molecular genetics of Calidris, with special reference to Knots. Wader Study Group Bull. 1992, 64 (Supplement): 29-35.
Byrkjedal I, Thompson D: Tundra Plovers: The Eurasian, Pacific and American Golden Plovers and Grey Plover. 1998, Princeton: Princeton University Press
Christian PD, Christidis L, Schodde R: Biochemical systematics of the Australian dotterels and plovers (Charadriiformes: Charadriidae). Aus J Zool. 1992, 40: 225-233.
Cohen BL, Baker AJ, Blechschmidt K, Dittmann DL, Furness RW, Gerwin JA, Helbig AJ, De Korte J, Marshall HD, Palma RL, Peter H-U, Ramli R, Siebold I, Willcox SM, Wilson RH, Zink RM: Enigmatic phylogeny of skuas (Aves: Stercorariidae). Proc Roy Soc Lond B. 1997, 264: 181-190. 10.1098/rspb.1997.0026.
Crochet P-A, Bonhomme F, Lebreton J-D: Molecular phylogeny and plumage evolution in gulls (Larini). J Evol Biol. 2000, 13: 47-57. 10.1046/j.1420-9101.2000.00135.x.
de Knijff P, Denkers F, van Swelm ND, Kuiper M: Genetic affinities within the herring gull Larus argentatus assemblage revealed by AFLP genotyping. J Mol Evol. 2001, 52: 85-93.
Dittmann DL, Zink RM: Mitochondrial DNA variation among phalaropes and allies. Auk. 1991, 108: 771-779.
Dittmann DL, Zink RM, Gerwin JA: Evolutionary genetics among phalaropes. Auk. 1989, 106: 327-331.
Ericson PGP: Systematic relationships of the palaeogene family Presbyornithidae (Aves: Anseriformes). Zool J Linn Soc. 1997, 121: 429-483. 10.1006/zjls.1997.0098.
Hacket SJ: Effects of varied electrophoretic conditions on detection of evolutionary patterns in the Laridae. Condor. 1989, 91: 73-90.
Haffer J: Systematik und Taxonomie der Larus argentatus Argtengruppe. In Handbuch der Vogel Mitteleuropas. Edited by: Glutz von Blotzheim UM, Bauer KM. 1982, Wiesbaden: Akademische Verlagsgesellschaft, 8: 502-515.
Heidrich P, Ristow D, Wink M: Molecular differentiation of Cory's and Manx shearwaters (Calonedris diomedea, Puffinus puffinus, P. yelkouan) and the herring gull complex (Larus argentatus, L. fuscus, L. cachinnans). J Ornith. 1996, 137: 281-294.
Helbig AJ: Mitochondrial genetic variation in gulls and shearwaters: a comment. J Ornith. 1994, 135: 609-615.
Jehl JR: Relationships in the Charadrii (shorebirds): a taxonomic study based on color patterns of the downy young. San Diego Soc Nat Hist. 1968, Memoir 3: 1-54.
Mindell DP, Sites JW: Tissue expression patterns of avian isozymes: a preliminary study of phylogenetic applications. Syst Zool. 1987, 2: 137-152.
Randi E, Spina F: An electrophoretic approach to the systematics of Italian gulls and terns (Aves, Laridae and Sternidae). Ital J Zool. 1987, 21: 317-344.
Ward D: The behavioural and morphological affinities of some vanelline plovers (Vanellinae: Charadriiformes: Aves). J Zool. 1992, 228: 625-640.
Wink M, Kahl U, Heidrich P: Genetic distinction of Larus argentatus, L. fuscus, and L. cachinnans. J Ornith. 1994, 135: 73-80.
We thank for Davide Pisani for stimulating discussion of supertree methods during the early planning of this work, Gareth Dyke for information on recent shorebird fossil discoveries, and three anonymous reviewers for helpful comments on an earlier version of the manuscript. Gavin Thomas was supported by a University of Bath Research Studentship.
GHT assisted in the design of the study, carried out the phylogenetic analyses and node dating, and drafted the manuscript in partial fulfillment of a doctoral degree at the University of Bath. MAW assisted in the design of the study and with editing and revision of the manuscript. TS assisted in the design of the study, collection of source trees, and editing and revision of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Estimates of node ages and node support (branch lengths.xls) Node numbers correspond to figures 2-9. Five types of estimate were used: a) absolute dates from the fossil record; b) absolute dates from molecular point estimates; c) relative dates based on branch length estimates from molecular studies; d) estimates based on a pure birth model (see text for details); and e) even spacing of nodes along branches with daughters and ancestors of known age. The numbers of characters supporting each node are provided (column D), this is equivalent to the number of source trees that share the equivalent node (see text for details). (XLS 32 KB)
Additional File 2: Shorebird supertree (50% majority-rule consensus; majrulesupertree.tiff) Shorebird supertree based on 50% majority-rule consensus of 1496 shortest trees with calibrated branch lengths. Scale bar indicates time from the present in millions of years. (TIFF 313 KB)
Additional File 3: Shorebird supertree (strict consensus; strictsupertree.tif) Shorebird supertree based on 50% majority-rule consensus of 1496 shortest trees. (TIFF 305 KB)
Additional File 4: Calibrated supertree (shorebirdsupertree.txt) The supertree in nexus format including branch length estimates. (TXT 24 KB)
Additional File 5: Source trees (source trees.xls) A summary of each tree used is given including the data type and main taxa studied. This is a brief summary and the original papers should be consulted for full details. (XLS 14 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Thomas, G.H., Wills, M.A. & Székely, T. A supertree approach to shorebird phylogeny. BMC Evol Biol 4, 28 (2004). https://doi.org/10.1186/1471-2148-4-28
- Branch Length
- Source Tree
- Total Evidence
- Strict Consensus Tree
- Recent Molecular Study