Inter- and intrafamilial relationships within the New World Piciformes were examined through an electrophoretic analysis of 20 protein-coding loci (19 of which varied between taxa). One individual from each of 26 species representing 25 genera and 5 families was analyzed; Momotus momota (Coraciiformes, Momotidae) was used as an outgroup. Although levels of genetic differentiation were high (the mean Nei's unbiased distance was 1.07), the data proved useful for phylogenetic inference. The jackknife technique was used to estimate the robustness of phylogenetic hypotheses. At the interfamilial level, the results suggest the following groupings: [[[(Bucconidae) (Galbulidae)]{(Picidae)[(Capitonidae) (Ramphastidae)]}]]. These results were consistent with hypotheses of familial relationships proposed by two recent cladistic analyses of morphological character complexes (Simpson and Cracraft 1981, Swierczewski and Raikow 1981). Our data challenge the currently accepted monophyly of the Piciformes, however, in much the same way as do DNA-DNA hybridization data. Agreement among independently derived hypotheses of interfamilial relationships suggests confidence in our knowledge of evolutionary patterns among piciform taxa. Hypotheses of intrafamilial relationships, some of which agreed with morphological patterns obtained in other studies, were presented. This study shows that starch-gel electrophoresis may be useful at higher taxonomic levels. Received 10 November 1986, accepted 4 June 1987.

Museum of Zoology, Louisiana State University, Baton Rouge, Louisiana 70803 USA RECENT years have witnessed a plethora of biochemical systematic studies (Avise and Aquadro 1982). These studies have demonstrat- ed that molecular characters have both advan- tages and limitations for phylogeny reconstruc- tion (Lanyon 1985b). There are several types of limitations. For example, molecular methods are scale dependent. That is, one would generally not use starch-gel electrophoresis of proteins to investigate relationships of higher taxa (Buth 1984) because often taxa share no alleles, which contributes no phylogenetic information. In other instances molecular methods may be un- able to resolve evolutionary patterns because of the relationship between the nature and rate of character evolution (e.g. allelic substitutions oc- curring more or less uniformly over time) and the nature of the evolutionary history itself (Fi- ala and Sokal 1985). For example, if cladogenetic events are close in time and lineages are short- lived before fragmentation (speciation), then there is a low probability of the origin (and retention through evolutionary time, given a  Present address: Division of Birds, Field Museum of Natural History, Roosevelt Road at Lake Shore Drive, Chicago, Illinois 60605 USA. fast rate of change) of a synapomorphic char- acter state, such as a particular allele. Nonetheless, the potential of biochemical methods for the inference of evolutionary his- tory is widely appreciated (Buth 1984, Wilson et al. 1985, Sibley and Ahlquist 1986) because they allow direct access to genetic information. A framework of quantitative models provides objective methods for estimating phylogenies (Felsenstein 1982). Probably few phylogenetic estimates, however, irrespective of the data or methods used to construct them, are immune to bias. Recently, statistical or quantitative tests have been developed to assess the accuracy of a particular phylogenetic pattern (e.g. Temple- ton 1983, Felsenstein 1985, Lanyon 1985a). Fur- thermore, systematic studies that compare pat- terns of morphological, biochemical, and behavioral variation may lead to more robust estimates of evolutionary history (Wagner 1961, Miyamoto 1981, Lanyon and Lanyon 1986). Such multidisciplinary approaches can exploit posi- tive aspects of each suite of characters, and per- mit evaluation of whether consistent patterns emerge. A disadvantage of this approach is that it requires investigators to become familiar with many different techniques, which can increase greatly the time required to complete studies of even small assemblages. It is possible to sim- ulate such a multidisciplinary study by inte- grating independently or collaboratively gath- ered data sets of several investigators. We examined phylogenetic patterns in elec- trophoretic characters (allozymes) in the Pici- formes. Recent hypotheses (Fig. 1) of piciform relationships have been based on hindlimb musculature (Swierczewski and Raikow 1981), osteological characters (Simpson and Cracraft 1981), and the anatomy of the feeding apparatus (Burton 1984). We compare phylogenetic pat- terns derived from analysis of our allozyme data set with those derived from analyses of muscle and skeletal variation. We illustrate the impor- tance of examining the robustness of phylo- genetic estimates through the use of the jack- knife technique (Lanyon 1985a). Furthermore, this is one of the few systematic applications of protein electrophoresis at higher taxonomic levels in birds (Lanyon 1985b); we suggest that the limits of the technique for avian systematics are unknown. METHODS Taxon selection.--The 27 tissue samples (represent- ing 25 genera and 5 families and 1 outgroup; see Table 1) analyzed in this study were collected over a period of 3 yr by personnel of the Louisiana State University Museum of Zoology. From 2 to 7 individuals, each in a different genus, were selected to represent each of the five traditionally recognized New World piciform families [Bucconidae (puffbirds), Galbulidae (jaca- mars), Capitonidae (barbers), Ramphastidae (toucans), and Picidae (woodpeckers)]. We maximized the num- ber of species-level taxa used, instead of individuals, to estimate levels and patterns of genetic variation within and among families. The use of a single in- dividual/taxon to estimate genetic distances is suffi- cient when investigating relationships at higher taxo- nomic levels (Gorman and Renzi 1979, Lanyon 1985b). The outgroup (Momotus momota, Momotidae) was se- lected from the presumed sister group, the Coraci- iformes (Simpson and Cracraft 1981, Swierczewski and Raikow 1981; but see Olson 1983). Electrophoresis.--Homogenates were prepared from pooled samples of liver and pectoral muscle: 0.5-1.0 g of minced tissue and an equal volume of Cleland's Reagent were centrifuged at 12,000 rpm for 60 min, the supernatant removed and preserved at - 70øC, and the pellet discarded. Homogenates were applied to 12% horizontal starch gels using filter-paper wicks. Two discontinuous and two continuous buffer sys- tems were used to isolate proteins using horizontal starch-gel electrophoresis [buffers 8-10 of Aquadro and Avise (1982); and Tris-maleate of Selander et al. Fig. 1. Cladogram of higher-level relationships supported by Simpson and C racraft (1981) and Swier- czewski and Raikow (1981). A = Ramphastoidea, B = Picoidea, C = Pici, D = Galbulae, and E = Piciformes. Derived characters support each node. BUCCONIDAE GALBULIDAE CAPITONIDAE RAMPHASTIDAE PICIDAE INDICATORIDAE (1971), adjusted to pH 6.5]. Twenty proteins (see Table 1), presumed homologous across taxa, were identified using protein-specific assays outlined by Harris and Hopkinson (1976). During the initial survey of loci, the mobility of each character (=electromorph) across the 27 taxa was recorded relative to a standard. All alleles of similar mobility were then compared in side-by-side tests on additional gels to ensure the ac- curacy of character-state designations. Data analysis.--We used the computer program BIOSYS-1 (Swofford and Selander 1981) to compute Rogers' (1972) and Nei's (1978) genetic distances and to estimate phylogenetic patterns using the UPGMA and distance Wagner procedures. The computer pro- gram PHYLIP (Felsenstein 1985) was used to produce Fitch-Margoliash (F-M) trees. These distance analyses yield phylogenetic information because of the pre- dominantly stochastic manner in which the charac- ters evolve (Kimura 1983). We note, however, that rates of nucleotide substitutions may differ between lineages (Avise and Aquadro 1982, Britten 1986), which results in reduced accuracy of phylogenetic hypoth- eses derived using distance analyses (Felsenstein 1978, 1982). To estimate the degree of error in phylogenies con- structed by our analyses of genetic distances, we em- ployed a jackknife analysis (Lanyon 1985a), in which pseudoreplicate distance matrices were created (each consisting of all possible combinations of n - 1 taxa, where n = the total number of taxa under consider- ation), and F-M trees generated. A single strict con- sensus tree was then produced to identify the nodes that were consistent with all analyses. Values reported are means ñ SD. RESULTS Allelic frequencies.--Inspection of the pattern of electromorphic variation at the 20 loci (Table 1) reveals considerable variation among the taxa. An average of 5.3 + 2.1 alleles/locus was ob- T^BLE 1. Distribution of electromorphs (denoted by lowercase letters) in 26 piciform taxa and 1 outgroup (Mornotus rnornota, Momotidae). Taxa 2-9 represent the family Bucconidae, 10-13 the Galbulidae, 14-20 the Picidae, 21-22 the Capitonidae, and 23-27 the Ramphastidae. Genus LOCUS a,b 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 1. Momotus momota 2. Bucco capensis 3. Nystalus radiatus 4. Malacoptila fusca 5. Micromonacha lanceolata 6. Nonnula rubecula 7. Hapaloptila castanea 8. Monasa nigrifrons 9. Chelidoptera tenebrosa 10. Galbula aIbirostris 11. Jacamerops aurea 12. Brachygalba salmoni 13. Galbalcyrhynchus leucotis 14. Carnpephilus haematogaster 15. Picoides scalaris 16. Sphyrapicus varius 17. Picoides villosus 18. Melanerpes aurifrons 19. Colaptes auratus 20. Picumnus borbae 21. Capito niger 22. Eubucco bourcierii 23. Selenidera spectabilis 24. Aulacorhynchus prasinus 25. Andigena hypoglauca 26. Pteroglossus castanotis 27. Rarnphastos sulfuratus g b b e c f b e b g e f h cd c c b ab b a e a a d c f c cf f e e d gh a e c ab c b a e a a d c f c e g e e d i a e c d bc b a e a d d a f c e d e e d i c e a d c b a e a c d c f c e f e f d i c e a d c b a e a d d c a c e f e e d h c e a d c b a e a d d c f a e f e e d i c e a d bc c a e a a d c f c e f e c d i c d b d c b a e a a d c f d e f e c d i b d c d c b a f a d d c g a e e e f f d c e a d b b a f b e d c b a e a d g g h d e a d b b a f a d e c e a e e g g d j c e a d b b a f a a d c b a g g e e f f ac e a d b b a b c f f b g c b e a b d g c f c c d b a b c f c b g c e h f a k c c d c c c a a b c f c c g c b h fad f a c d c c c b a a c f c b g c bd h f a k c c d c c c b a ab c f c d g d b h b ac d bc c d c c c b a b c f c d g c g h f a g cd c e c c c b a b c d c b g c g h h b f d d f c bc c b a h c f d b g c g f f c k i c b a b c b a h c f d b g c g i a c k h c b b b c b a h c f b a g c g c a e d d c b a b c b a h c f b a g c g f c c f d d b a b c b a dh c f b a g e g c c ce f f c b b b c b a h c f b a g c g c c c g d c b a b c a a h c d a a g c g c e c f d a b b b c b a  Loci in order are: Pgm-1 (E.C. 2.7.5.1), Mdh-l,2 (1.1.1.37), Got-I (2.6.1.1), Ck-2 (2.7.3.2), Eap (3.1.3.2), Pgi (5.3.1.9), Est "D" (3.1.1.1), Sod-I (1.15.1.1), Gda (3.5.1.2), 6-Pgd (1.1.1.44), Peptidase C, B (3.4.11), Mpi (5.3.1.8), Ldh-l,2 (L1.1.27), ldh-l,2 (1.1.1.42), Fum (4.2.1.2), and Lap (3.4.11.1). b Two letters signify a heterozygous genotype. served; only one locus (Lap) was monomorphic across all taxa. Ten other loci were examined but were too variable to warrant continued analysis from a logistical standpoint. Genetic distances.--A summary of Nei's (1978) genetic distances (Table 2) shows average ge- netic distances within and between major groups (the complete matrix is available from the au- thors). Because the distance values in Table 2 were based on conservativeloci, theyare under- estimates of the actual level of genetic differ- entiation in piciform birds; this should be con- sidered when comparing our results with values reported in the literature. The distances are, nevertheless, greater than those reported for comparisons at comparable taxonomic levels in other birds (Barrowclough 1980). The average Nei (1978) genetic distance among all piciform taxa was 1.07. Distance analysis and interfamilial phylogenetic patterns.--The UPGMA and distance Wagner analyses (Fig. 2) and Fitch-Margoliash analysis (not shown) produced the same branching se- quences for the five families. Furthermore, the branching sequence was completely stable to the jackknife analysis. Three hypotheses of in- terfamilial relationships proposed indepen- dently by Swierczewski and Raikow (1981) and by Simpson and Cracraft (1981)were supported. The Bucconidae and Galbulidae were identified as sister taxa, supporting a hypothesis previ- ously defined by four osteological characters (Simpson and Cracraft 1981), six myological characters (Swierczewski and Raikow 1981), and similarities in conalbumins and ovalbumin (Sibley and Ahlquist 1972). Although these taxa have consistently been identified as close tel- atives, the evidence presented here reveals marked genetic dissimilarity (D = 0.80 + 0.22), and a considerable amount of time probably has elapsed since they shared a common ancestor. This is in contrast to the Capitonidae and Ram- TABLE 2. Nei's (1978) genetic distance (+ SD) between major groupings of piciform taxa. Entries in the diagonal are comparisons between taxa within each grouping. Sample sizes are the number of pairwise comparisons. Motmot Puffbirds Jacamars Barbets Toucans Woodpeckers Motmot -- Puffbirds 1.15 + 0.10 0.32 + 0.09 (Bucconidae) n = 7 n = 28 Jacamars 1.08 + 0.08 0.80 + 0.22 0.52 + 0.11 (Galbulidae) n = 4 n = 32 n = 6 Barbets 1.57 1.12 + 0.24 1.50 + 0.20 (Capitonidae) n = 2 n = 16 n = 8 Toucans 1.55 + 0.25 1.38 + 0.29 1.59 + 0.37 (Ramphastidae) n = 5 n = 40 n = 20 Woodpeckers 1.55 + 0.19 1.38 + 0.22 1.75 + 0.36 (Picidae) n = 7 n = 56 n = 28 0.22 n=l 0.48 + 0.07 0.32 + 0.07 n = 10 n = 10 1.26 + 0.49 0.98 _+ 0.19 n = 14 n = 40 0.44 + 0.18 n = 21 phastidae, also recognized as sister taxa, be- tween which the average intergeneric distance was considerably smaller (D = 0.48 + 0.07). The monophyly of the Capitonidae and Ram- phastidae, supported by two myological char- acters (Swierczewski and Raikow 1981) and one osteological synapomorphy (Simpson and Cra- craft 1981), was also supported by the electro- phoretic analysis. We suggest that capitonids and ramphastids are as closely related as are the jacamar genera (D = 0.52 + 0.11). The electrophoretic data supported a mono- phyletic assemblage consisting of the Picidae, Capitonidae, and Ramphastidae. This hypoth- esis has been supported previously by four osteological characters, six myological syn- apomorphies, and the mobility of an s-Mdh zymogram (Avise and Aquadro 1987). A close relationship between capitonids and picids was suggested by analysis of egg-white proteins (Sibley and Ahlquist 1972), but ramphastids were not analyzed. In contrast, Peters (1964) placed the capitonids and ramphastids with the bucconids and galbulids rather than with the picids. Sibley and Ahlquist (1985) suggested that the New World capitonids were more closely related to ramphastids than to Old World cap- itonids. We had tissue samples only from New World capitonids, and we cannot comment on the monophyly of the Capitonidae. To evaluate the familial relationships sup- ported by this study, one specimen from each of the five families (plus Momotus) was selected at random and a Fitch-Margoliash analysis con- ducted. Ten iterations were performed, and in all cases the results (not shown) supported the branching sequence of families shown in Fig. 2. IntrafamiIiaI phylogenetic relationships.--To in- vestigate relationships within piciform fami- lies, we assumed that each family represented a monophyletic assemblage. A jackknifed F-M strict consensus tree was produced for each fam- ily and was rooted using the remaining taxa. At this lower level of taxonomic investigation, we detected tree instability. Jackknifing allowed us to expose inconsistencies within the distance matrix, identified as unstable nodes, and to iso- late reliable nodes within families (Figs. 3-5). Four nodes were retained for the eight buc- conid genera (Fig. 3a). The hypothesis that Buc- co and NystaIus represent sister taxa was sup- ported [Peters (1964) and Cottrell (1968) considered them to be congeners]. Swierczew- ski and Raikow (1981) identified NystaIus and CheIidoptera as being more closely related to each other than either is to MaIacoptiIa or NonnuIa (Fig. 3b). This was neither supported nor re- futed by the electrophoretic investigation, which failed to resolve the branching sequence for these taxa. The electrophoretic characters also identified MaIacoptiIa, NonnuIa, HapaIoptiIa, and Micromonacha as a monophyletic assem- blage within which Micromonacha is the sister taxon to the remaining three genera. Finally, we note that the distance between Bucco and NystaIus, which are currently considered con- generic (A.O.U. 1983), is comparable to the dis- tance between CheIidoptera and Monasa (D = 0.418 and 0.387, respectively). The latter two taxa have not previously been suggested as sis- ter taxa (or congeners) but have been placed .O8 .21 .66 .13 Momotus .10 .20 Bucco .0  Nystalus .15 Monasa ß 2  .15 Ch(c)lldoptera .20 .03  Malacoptlla  .20 Mlcromonacha Nonnula .23 Hapaloptlla ß 10 .30 Galbula  4 Brachygalba .10 ß Galbalcyrbynchus .45 Jacamerops .47 Campephllus .12 P. scalarIs .01 P. vlllosus .05 Sphyraplcus Colapt(c)s Melanerpes .48 Plcumnus .20 Caplto  Eubucco Selenldera Pteroglossus .08 Aulacorbynchus .30 Andlgena .3 Ramphastos b .07 .O8 Momotus .05 Bucco ß 01 Ny$ talus .05 Chelldoptera .09 Mlcromonacha .09 Nonnula .07 .26 dacamerops Galbalcyrhynchus .25 Camp(c)phllus .07 Sphyraplcus Melanerpes .23 Plcumnus .06 Caplto  .12 Euhucco ß 14 Aulacorhynchus  .17 Ramphastos  .14 Pteroglossus Fig. 2. UPGMA (a) and optimized distance Wagner (b) trees derived using Rogers' genetic distance mea- sures calculated from the matrix of electromorphs (Table 1). The cophenetic correlation coefficient for the UPGMA phenogram is 0.93, indicating a reasonably good fit of the branching diagram to the distance matrix. Farris' "f" for the distance Wagner tree is 14.2. The values on the diagrams represent branch lengths. The distance Wagner tree was rooted using Momotus as an outgroup. See legend to Table 1 for family membership of taxa. together at the end of the family in recent clas- sifications (Peters 1964, A.O.U. 1983). Although genetic distance cannot be an absolute measure of taxonomic status, our data suggest scrutiny of the postulated close relationship between Bucco and Nystalus. None of the galbulid phylogenetic hypoth- eses generated by the distance analyses (Fig. 2) was retained after jackknifing. Despite this po- tential lack of confidence in the branching pat- tern depicted (it might, after all, be correct), the allelic frequency data do provide insight into galbulid evolution. The relatively great genetic distance between galbulid genera, 0.52 +_ 0.11, suggests that the four genera analyzed herein are from distinct, relatively old lineages. These lineages might have originated close in time. Two nodes were consistently supported for -  Mlcromonacha I . ...... Malacoptlla I - ..... Hapaloptlla %  -- Bucco  Nonnula / / Malacoptlla   Chelldop tera Nystalu$ Fig. 3. (a) Jackknifed strict consensus Fitch-Mar- goliash tree for the Bucconidae. All other taxa were used as a composite outgroup to root the tree. Dotted lines indicate portions of the topology that were un- stable to the jackknife manipulation. (b) Phylogenetic relationships within the Bucconidae supported by Swierczewski and Raikow's (1981) analysis of hind- limb musculature. the five ramphastid taxa examined (Fig. 4a). Ramphastos was identified as the outgroup to the remaining genera, and Selenidera and Pteroglos- sus were identified as sister taxa (note, however, that these results differ from the distance Wag- ner tree shown in Fig. 2b). Our findings conflict with those of Swierczewski and Raikow (1981), who concluded that Pteroglossus and Ramphastos were sister taxa (Fig. 4b), but are consistent with Haffer's (1974) phylogeny for part of the family (Fig. 4c). Haffer (1974) suggested that the low- land genera Selenidera and Pteroglossus were sis- ter taxa on the basis of vocal similarities. Three apparently robust phylogenetic hy- potheses were generated for the seven wood- pecker taxa analyzed (Fig. 5a). The identifica- tion of Picumnus as the sister group to the remaining forms is consistent with Short's (1982) classification and the findings of Swierczewski and Raikow (1981) (Fig. 5b). Campephilus was shown to be the next lineage to arise (relative to those surveyed), as suggested by Swierczew- ski and Raikow (1981). In the classification of Short (1982) the relationships among the Me- lanerpini, Campetherini, and Campephilini / Andlgena / Aulecorhynchus  %_ .... Selenldere Pteroglossus Ramhastos b  Aulacorhynchus Pteroglossus  Ramphastos 1/2 Aulacorhynchus --  Selenldera Pteroglossus Fig. 4. (a) Jackknifed strict consensus Fitch-Mar- goliash tree for the Ramphastidae. All other taxa were used as a composite outgroup to root the tree. Dotted lines indicate portions of the topology that were un- stable to the jackknife manipulation. (b) Phylogenetic relationships within the Ramphastidae supported by Swierczewski and Raikow's (1981) analysis of hind- limb musculature. (c) Phylogenetic relationships within the Ramphastidae supported by Haffer's (1974) analysis of plumage, vocalizations, and biogeogra- phy. were unresolved. Lastly, Picoides villosus and P. scalaris were identified as sister taxa, as one would expect from their congeneric status. Sev- eral hypotheses suggested by Short's classifi- cation and the cladogram developed by Swier- czewski and Raikow (1981) cannot be addressed here. Specifically, we lack a consistently sup- ported branching sequence for Melanerpes, Sphyrapicus, Colaptes, and Picoides. Additional study of patterns of protein variation promises insights into picid relationships. DISCUSSION Polyphyly of the Piciformes.--Sibley and Ahl- quist (1972, 1985), Olson (1983), and Burton (1984) have suggested that the Piciformes might be polyphyletic. These authors suggest that the / P. calarls ......... Melanerpes Campephllus Plcumnus ___ P. scalarIs f $phyraplcus / i ' Melanerpes Colapres Campephllus Plcumnus Fig. 5. (a) Jackknifed strict consensus Fitch-Mar- goliash tree for the Picidae. All other taxa were used as a composite outgroup to root the tree. (b) Phylo- genetic relationships within the Picidae supported by Swierczewski and Raikow's (1981) analysis of hind- limb musculature. Galbulae (Bucconidae and Galbulidae) might be more closely related to some Coraciiformes than either is to the Picae (Picidae, Ramphastidae, and Capitonidae). We did not specifically ad- dress this question. We can consider the ques- tion of monophyly, however, under the as- sumption of a strong positive correlation between genetic distance and time since diver- gence. The UPGMA phenogram, which as- sumes a strong time correlation, places Momotus as the sister taxon to the Galbulae. The reason for this topology can be seen from the summary of genetic distances (Table 2). The average dis- tances between Momotus and the galbulid and bucconid genera (1.08 + 0.08 and 1.15 + 0.10, respectively) are considerably smaller than the average distance between the Picae genera and Momotus (D = 1.55 + 0.19). Therefore, our find- ings are consistent with the possibility that the order Piciformes (sensu Peters 1964) is polyphy- letic. Testing the stability of phylogenetic trees.--Sev- eral sources of bias cause phylogenetic trees to be unstable or unreliable. A particular phylo- genetic estimate is dependent on the samples of individuals, characters, and taxa used in its construction; different samples of each might result in different phylogenetic estimates. In addition, the particular algorithm used has in- herent assumptions, such as an average uniform rate of change (UPGMA) or parsimony (dis- tance Wagner). It usually is not possible to know the nature of character evolution and therefore to select the appropriate algorithm. For some algorithms (e.g. distance Wagner), several equivalent or nearly equivalent trees can result. Thus, methods of testing robustness of trees are being developed (Templeton 1983, Felsenstein 1985, Lanyon 1985a). It is sometimes not appreciated that methods of producing branching diagrams do so without any implied confidence in the branching pat- tern of either the entire tree or parts thereof. Often, a tree's structure depends on few char- acters, and sampling additional characters or individuals shifts branching patterns. System- atists should evaluate routinely the robustness of particular nodes in phylogenetic estimates. Several factors merit comment, however. The patterns obtained from the first pass of an al- gorithm (e.g. UPGMA) might be correct in their entirety; such a conclusion would be enhanced if a similar branching pattern was obtained from an independent data set. Because some nodes of the phylogeny might be based on a single or very few characters, such nodes might disap- pear wrongly with some resampling method, such as jackknifing. In other words, an ancestral branch might be very short lived, with time for few character-state transitions. If one happened to analyze the (derived) characters that define such a branch, the fact that only a few characters (could) support the branch might cause it to be collapsed in tests of robustness. In addition, other methods of testing the stability or confi- dence of branching diagrams, such as boot- strapping (Felsenstein 1985), might indicate in- stability in a consensus tree obtained by jack- knifing. Therefore, we note that our branching patterns for all taxa might be correct, but we limit interpretation to those patterns for which some confidence is indicated by jackknifing. We suggest that at least one test of the robustness of branching patterns be employed, along with conservative interpretation of patterns. Few avian systematic studies have used explicit, ob- jective methods to examine the robustness of phylogenetic hypotheses (Lanyon 1985b). Lastly, we note the ongoing debate concern- ing the inference of phylogenies from distance matrices (e.g. Farris 1986, Felsenstein 1986). Resolution of this issue, we suspect, will not alter our primary conclusions. Suggestions for further work on piciform relation- ships.--As a result of this study and indepen- dent analyses of two very different suites of characters (Simpson and Cracraft 1981, Swier- czewski and Raikow 1981), a set of reliable hy- potheses of systematic relationships has been identified for the Piciformes at the family level and above. With the exception of the Picidae, however, for which Short (1982) constructed a relatively detailed classification, relationships within the piciform families remain obscure. Using the reliable hypotheses of familial rela- tionships to identify proper outgroups and the preliminary intrafamilial hypotheses presented here and by Swierczewski and Raikow (1981) as a foundation, studies of phylogenetic rela- tionships within each piciform family may be designed. For example, it will be useful to de- termine whether New World capitonids are more closely related to ramphastids or to Old World capitonids (see Sibley and Ahlquist 1985). Several taxa representing the range of variation within the Picidae should be used as a com- posite outgroup. A critical higher-level systematic question that remains to be answered definitively is whether the Piciformes are monophyletic. Allozymic data suggest that the traditional Piciformes may be polyphyletic (Fig. 2a), although further assess- ment of this question must await a study in which a composite outgroup is constructed from closely related orders. We note also the useful- ness of electrophoresis of proteins (allozymes) for estimating phylogenetic relationships at higher taxonomic levels, levels at which the technique supposedly fails (Buth 1984). Given that this approach is less expensive than others (e.g. DNA-DNA hybridization), it deserves ap- plication at higher taxonomic levels in birds, especially when it can complement other mo- lecular and nonmolecular data sets. ACKNOWLEDGMENTS We thank the Louisiana State University Museum of Zoology for financial support and the personnel of the Museum for collecting the tissue samples used in this study. We also thank J. C. Avise and C. F. Aquadro for providing a prepublication copy of their manuscript. R. L. Cann, A. P. Capparella, J. Cracraft, J. A. Gerwin, M. S. Hafner, S. J. Hackett, and J. V. Remsen, Jr., commented constructively on the manu- script. LITERATURE CITED AMERICAN ORNITHOLOGISTS' UNIONß 1983. Check-list of North American birds, 6th ed. Washington, D.C., Amer. Ornithol. Union. AQUADRO, C. F., &: J. C. AVISE. 1982. Evolutionary genetics of birds. VI. A reexamination of protein divergence using varied electrophoretic condi- tions. Evolution 36: 1003-1019. AVISE, J. C., &: C. F. AQUADRO. 1982. A comparative summary of genetic distances in the vertebrates. Pp. 151-185 in Evolutionary biology, vol. 15 (M. Hecht, B. Wallace, and R. Prans, Eds.). New York, Plenum Press. , &--. 1987. Malate dehydrogenase iso- zymes provide a phylogenetic marker for the Pic- iformes (woodpeckers and allies). Auk 104: 324- 328. BARROWCLOUGH, G.F. 1980. Genetic and phenotypic differentiation in a wood warbler (genus Den- droica) hybrid zone. Auk 97: 789-798. BRITTEN, R. J. 1986. Rates of DNA sequence evolu- tion differ between taxonomic groups. Science 231: 1393-1398. BURTON, P. J. K. 1984. Anatomy and evolution of the feeding apparatus in the avian orders Cora- ciiformes and Piciformes. Bull. Brit. Mus. Nat. Hist. (Zool.) 47: 331-443. BUTH, D.G. 1984. The application of electrophoretic data in systematic studies. Ann. Rev. Ecol. Syst. 15: 501-522. COTTRELL, G.W. 1968. The genera of puffbirds. Bre- viora, No. 285. FARRIS, J. S. 1986. Distances and statistics. Cladistics 2: 144-157. FELSENSTEIN, J. 1978. Cases in which parsimony or compatibility methods will be positively mis- leading. Syst. Zool. 27: 401-410. ß 1982. Numerical methods for inferring evo- lutionary trees. Q. Rev. Biol. 57: 379-404. ß 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783- 791. 1986. Distance methods: a reply to Farris. Cladistics 2: 130-143. FIALA, K. L., & R. R. SOKAL. 1985. Factors determin- ing the accuracy of cladogram estimation: eval- uation using computer simulation. Evolution 39: 609-622. GORMAN, G. C., & J. R. RENZ1. 1979. Genetic distance and heterozygosity estimates in electrophoretic studies: effects of sample size. Copeia 1979: 242- 249. HAFFER, J. 1974. Avian speciation in tropical South America. Publ. Nuttall Ornithol. Club No. 14. HARRIS, H., & D. A. HOPKINSON. 1976. Handbook of electrophoresis in human genetics. Amsterdam, North-Holland Publ. Co. KIMURA, M. 1983. The neutral theory of molecular evolution. Cambridge, England, Cambridge Univ. Press. LANYON, S. g. 1985a. Detecting internal inconsis- tencies in distance data. Syst. Zool. 34: 397-403. --. 1985b. Molecular perspective on higher-level relationships in the Tyrannoidea (Aves). Syst. Zool. 34: 404-418. LANYON, W. E., & S. g. LANYON. 1986. Generic status of Euler's Flycatcher: a morphological and bio- chemical study. Auk 103: 341-350. MIYAMOTO, g. 1981. Congruence among character sets in phylogenetic studies of the frog genus Leptodactylus. Syst. Zool. 30: 281-290. NEI, M. 1978. Estimation of average heterozygosity and genetic distance from a small number of in- dividuals. Genetics 89: 583-590. OLSON, S.L. 1983. Evidence for a polyphyletic origin of the Piciformes. Auk 100: 126-133. PETERS, J.L. 1964. Check-list of birds of the world, vol. VI. Cambridge, Massachusetts, Harvard Univ. Press. ROGERS, J. S. 1972. Measures of genetic similarity and genetic distance. Univ. Texas Publ. 7213: 145- 153. SELANDER, R. K., M. H. SMITH, S. Y. YANG, W. E. JOHNSON, & J. B. GENTRY. 1971. Biochemical polymorphism and systematics in the genus Pero- myscus. I. Variation in the old-field mouse (Pero- myscus polionotus). Univ. Texas Publ. 7103: 49-90. SHORT, L.L. 1982. Woodpeckers of the world. Green- ville, Delaware Mus. Nat. Hist. SIBLEY, C. G., & J. E. AHLQUIST. 1972. A comparative study of the egg white proteins of non-passerine birds. Bull. Peabody Mus. Nat. Hist., Yale Univ., No. 39. --, & ---. 1986. Reconstructing bird phy- logeny by comparing DNA's. Sci. Amer. 254: 82- 92. SIMPSON, S. F.,  J. CRACRAFT. 1981. The phyloge- netic relationships of the Piciformes (class Aves). Auk 98: 481-494. SWIERCZEWSKI, E. V., & R. J. RAIKOW. 1981. Hind limb morphology, phylogeny, and classification of the Piciformes. Auk 98: 466-480. SWOFFORD, D. L., & R. K. SELANDER. 1981. A computer program for the analysis of allelic variation in genetics. J. Hered. 72: 281-283. TEMPLETON, A.g. 1983. Phylogenetic inference from restriction endonuclease cleavage site maps with particular reference to the evolution of humans and apes. Evolution 37: 221-244. WAGNER, W. H., JR. 1961. Problems in the classifi- cation of ferns. Recent Adv. Botany 1: 841-844. WILSON, A. C., R. L. CANN, S. M. CARR, M. GEORGE, U. B. GYLLENSTEN, K. M. HELM-BYCHOWSKI, R. G. HIGUCHI, S. R. PALUMBI, E. M. PRAGER, R. D. SAGE, & G. M. STONEKING. 1985. Mitochondrial DNA and two perspectives on evolutionary genetics. Biol. J. Linn. Soc. 26: 375-400.