Research article | Open | Published:
Polyphyly of the hawk genera Leucopternis and Buteogallus (Aves, Accipitridae): multiple habitat shifts during the Neotropical buteonine diversification
BMC Evolutionary Biologyvolume 6, Article number: 10 (2006)
The family Accipitridae (hawks, eagles and Old World vultures) represents a large radiation of predatory birds with an almost global distribution, although most species of this family occur in the Neotropics. Despite great morphological and ecological diversity, the evolutionary relationships in the family have been poorly explored at all taxonomic levels. Using sequences from four mitochondrial genes (12S, ATP8, ATP6, and ND6), we reconstructed the phylogeny of the Neotropical forest hawk genus Leucopternis and most of the allied genera of Neotropical buteonines. Our goals were to infer the evolutionary relationships among species of Leucopternis, estimate their relationships to other buteonine genera, evaluate the phylogenetic significance of the white and black plumage patterns common to most Leucopternis species, and assess general patterns of diversification of the group with respect to species' affiliations with Neotropical regions and habitats.
Our molecular phylogeny for the genus Leucopternis and its allies disagrees sharply with traditional taxonomic arrangements for the group, and we present new hypotheses of relationships for a number of species. The mtDNA phylogenetic trees derived from analysis of the combined data posit a polyphyletic relationship among species of Leucopternis, Buteogallus and Buteo. Three highly supported clades containing Leucopternis species were recovered in our phylogenetic reconstructions. The first clade consisted of the sister pairs L. lacernulatus and Buteogallus meridionalis, and Buteogallus urubitinga and Harpyhaliaetus coronatus, in addition to L. schistaceus and L. plumbeus. The second clade included the sister pair Leucopternis albicollis and L. occidentalis as well as L. polionotus. The third lineage comprised the sister pair L. melanops and L. kuhli, in addition to L. semiplumbeus and Buteo buteo. According to our results, the white and black plumage patterns have evolved at least twice in the group. Furthermore, species found to the east and west of the Andes (cis- Andean and trans- Andean, respectively) are not reciprocally monophyletic, nor are forest and non-forest species.
The polyphyly of Leucopternis, Buteogallus and Buteo establishes a lack of concordance of current Accipitridae taxonomy with the mtDNA phylogeny for the group, and points to the need for further phylogenetic analysis at all taxonomic levels in the family as also suggested by other recent analyses. Habitat shifts, as well as cis- and trans-Andean disjunctions, took place more than once during buteonine diversification in the Neotropical region. Overemphasis of the black and white plumage patterns has led to questionable conclusions regarding the relationships of Leucopternis species, and suggests more generally that plumage characters should be used with considerable caution in the taxonomic evaluation of the Accipitridae.
The family Accipitridae comprises approximately 237 species of predatory birds distributed worldwide except Antarctica , with diversity concentrated in the Neotropics [1, 2]. Despite numerous taxonomic revisions (e.g., [3, 4]), the evolutionary history of the family has not been sufficiently explored using methods of phylogenetic inference, and current classifications are mainly based on plumage and ecological resemblance between taxa . Current taxonomy is still highly provisional at all taxonomic levels [1, 6], and does not appear to reflect phylogenetic relationships in several cases [7, 8], thus retarding biogeographic analysis, morphological trait mapping and the general understanding of the evolutionary history of the Accipitridae.
The Accipitridae morphological diversity has been traditionally represented in sub-groups of similar or supposedly closely related species, such as "kites", "harriers", "booted eagles" and "buteonines" . The buteonine hawks are represented by the large cosmopolitan genus Buteo and several related genera, called "sub-buteonines" by Amadon , which includes the predominantly Neotropical genera Buteogallus, Parabuteo, Asturina, Leucopternis, Busarellus, Geranoaetus, Geranospiza and Harpyhaliaetus. Two old world genera, Kaupifalco and Butastur, were formerly included as part of the "sub-buteonines" group, but were subsequently removed from this division . Some authors consider the buteonine as a sub-family (Buteonineae, e.g. Friedman , Grossman and Hamlet ), but formal sub-familial division of Accipitridae has been a contentious issue due to a lack of knowledge of the evolutionary history of the family (see [9, 12]).
Evolutionary biologists have long sought to understand the processes responsible for the generation of the high species richness found in the Neotropics, and several models of biotic diversification have been invoked to explain such patterns, for example forest refuges resulting from climatic fluctuations [13–15], rivers as barriers to gene flow , river dynamics , sea level oscillations [18–20], geotectonic vicariance  and ecological factors  (see Moritz et al.  for a revision). However, these models have only rarely been tested with organisms capable of long-distance dispersal (e.g., [24, 25]), such as hawks and eagles capable of soaring and gliding flight. Numerous flocks of migrant hawk species as Buteo platypterus and Buteo swainsoni, for example, cross the Andes as part of their yearly migrations , and call into question the degree to which the geographical barriers to gene flow identified in many models of Neotropical diversification have been important in Accipitridae speciation.
The genus Leucopternis is a morphologically heterogeneous group of 10 buteonine species distributed in forested habitats from southern Mexico to Paraguay and Uruguay , and offers an opportunity to explore the diversification of an Accipitridae group distributed throughout the Neotropical region. Species in the genus vary from the small L. semiplumbeus (250 g) to the large L. princeps (1 kg) , and are hawks with broad wings and medium to short tails. Two species, L. schistaceus and L. plumbeus, are entirely dark slate; however, most Leucopternis have primarily white plumage and vary in the amount of black, grey or slate black on the back, wings and/or head. Those which we here refer to as "black-and-white" Leucopternis species are L. albicollis, L. polionotus, L. occidentalis, L. lacernulatus, L. melanops, L. kuhli, L. semiplumbeus, and L. princeps [1, 6]. While Leucopternis are found exclusively in forest habitats, other Neotropical buteonine species occur in a variety of habitats, such as mangroves (Buteogallus aequinoctialis), savannahs (Harpyhaliaetus coronatus, Buteogallus meridionalis) and wetlands (Busarellus nigricollis) [1, 6], which makes this group suitable to analysis of the evolutionary relationships of forest and non-forest species.
The buteonine phylogeny has been partially explored recently using morphological and molecular data [7, 27, 28], but Neotropical species have not been well represented. Incomplete taxon sampling notwithstanding, these analyses have called into question the monophyly of Leucopternis, Buteo and Buteogallus [7, 28]. The present work constitutes an effort to clarify the relationships among all Leucopternis species, and their position relative to other Neotropical buteonine genera. We address the following questions: (1) Is Leucopternis as currently recognized monophyletic? (2) What are the relationships among species of Leucopternis to other genera of buteonine hawks? (3) Is the black and white plumage pattern a synapomorphic trait uniting the majority of species in the genus Leucopternis? (4) Are phylogenetic relationships among Neotropical buteonines predicted by biogeography or habitat?
Datasets, molecular variation
Our final alignment of the total dataset (12S, ATP8 and 6 and ND6) without gaps totalled 2179 base pairs, with 651 variable and 505 parsimony informative sites. Uncorrected distances ranged from 0 to 7.1% for 12S (without gaps), 0 to 21.4% for ATP8, 0 to 11.9% for ATP6, and 0 to 13.9% for ND6. Deviations from linearity were found in third position plots of ATP8, ATP6 and ND6. We did not detect significant departures from homogeneity of base frequencies across taxa in any dataset (P > 0.05, data not shown). We are confident of mitochondrial origin of our sequences because: (1) most of our samples were represented by mitochondrial-rich tissues (feathers, muscle or liver); (2) most samples (comprising 12 of 20 species) had all regions sequenced using independent fragments amplified with different primer sets (with sequence overlap ranging from 67 to more than 400bp, in highly variable regions), and sequences were identical; (3) sequences were easily aligned to published sequences of other Accipitridae species; (4) electropherograms were carefully checked for double peaks; (5) coding regions did not show unexpected stop codons; and 6) gene specific phylogenetic analyses revealed similar relationships to those inferred from the combined data, indicating that a mitochondrial translocation to the nucleus would have to have been more than eight kilobases in length.
A single nucleotide site in the 12S sequence of the muscle sample LGEMA F39 (L. lacernulatus) presented a strong "C" peak with a lower "A" peak at the position 593, and this same pattern persisted in sequences obtained from amplifications using three different primer combinations, with sizes ranging from approximately 800 to 2700 bp. The sequence can be easily aligned and has a base composition similar to other published sequences for the Accipitridae. We could not find any evidence of pseudogene amplification, and since it has been suggested that PCR amplifications larger than 1.5 kilobases are likely to represent true mitochondrial amplifications , this site was coded as "M" (IUPAC code representing C and A) in all analyses, and it may represent an example of mitochondrial heteroplasmy.
Phylogenetic analysis, single and combined datasets
None of the phylogenetic analyses supported the monophyly of the hawk genera Leucopternis, Buteogallus or Buteo. Phylogenies inferred from the subsets of mtDNA genes and combined data were largely congruent under all optimization criteria, and topologies differed mainly in resolution and nodal support. Although phylogenetic analysis of single gene subsets resulted in poorly resolved trees (data not shown), the nodes identified with high bootstrap support (> 75) or posterior probabilities (> 0.95) were entirely congruent with those identified in the combined analyses.
The maximum likelihood (ML) inference of the total dataset resulted in one completely resolved tree with likelihood -ln 9813.2993 (figure 1), which was identical to the majority rule consensus topology of the Bayesian analysis (BA). Most nodes were highly supported by both bootstrap proportions and posterior probabilities. Maximum parsimony (MP) analyses of the total dataset resulted in two equally parsimonious trees (1531 steps, CI = 0.4899, RI= 0.6938) (trees not shown), which were also highly congruent with the ML and BA inferences. The main differences between the MP trees and the other combined analyses were the position of L. princeps basal to Clades 2 and 3 plus the remaining Buteo species and Parabuteo unicinctus (bootstrap support 68, data not shown), and the position of L. schistaceus basal to a poorly supported clade containing the sister pairs Buteogallus meridionalis and L. lacernulatus, and H. coronatus and Buteogallus urubitinga (bootstrap support < 50, data not shown).
Three major mtDNA clades containing Leucopternis species (figure 1) were recovered in all combined analysis, with high levels of statistical support in most cases. Clade 1 included Buteogallus meridionalis, L. lacernulatus, Buteogallus urubitinga, H. coronatus, L. schistaceus and L. plumbeus; Leucopternis princeps was basal to this Clade 1 in the ML and BA analyses, but with low support. The ML and BA analyses strongly supported a sister relationship between L. lacernulatus and Buteogallus meridionalis, and between H. coronatus and B. urubitinga; the latter sister pair was also highly supported by parsimony bootstrap (bootstrap support 98, data not shown). There was weaker support for a sister relationship between L. schistaceus and Buteogallus meridionalis plus L. lacernulatus, obtained in the ML and BA analyses. All phylogenetic trees inferred from the combined dataset establish an early split of L. plumbeus from the rest of the Clade 1.
Clade 2 comprised the two sampled L. albicollis subspecies (which carried identical mtDNA sequences), L. occidentalis and L. polionotus. Clade 3 was sister to Clade 2 and included L. semiplumbeus, L. melanops and L. kuhli sister to Buteo buteo. Buteo albicaudatus, Buteo magnirostris, Parabuteo unicinctus and Buteo leucorrhous were outside the sister relationship of Clade 2 and Clade 3 in all analyses. There was strong support in the ML, MP and BA trees obtained from the combined dataset for a sister relationship between Buteo leucorrhous and P. unicinctus.
Shimodaira-Hasegawa tests were conducted on topologies constrained by the monophyly of all species in the genus Leucopternis, monophyly of species with the black and white plumage pattern, monophyly of all forest species, and reciprocal monophyly of the cis- and trans-Andean species. In all tests the constraint trees had a significantly poorer fit to the data than the unconstrained ML tree (P < 0.001). Ancestral state reconstructions of habitat and plumage traits onto the ML tree using unordered parsimony (Figure 2) showed that the black-and-white plumage characteristic to most Leucopternis species evolved at least twice in Neotropical buteonines, and shifts between forest and open habitats occurred at least four times.
Phylogenetic relationships, novel groupings, and evolution of common plumage patterns
The mtDNA-based phylogenies for Neotropical buteonines reject the monophyly of the genera Leucopternis, Buteogallus and Buteo. Our analysis provides another example of the lack of correspondence between classic taxonomic groupings within the Accipitridae and molecular phylogenies [7, 8, 28]. The genus Leucopternis is a composite of three independent lineages, and several Leucopternis species are more closely related to other buteonine taxa than to each other.
One of the novel phylogenetic arrangements presented here is the sister relationship between L. lacernulatu s and Buteogallus meridionalis, rather than between L. lacernulatus and the sympatric Atlantic forest endemic L. polionotus, despite overall similarities in plumage pattern between the latter pair (see ). A proposed close relationship [1, 4, 12, 30] between L. lacernulatus and the L. melanops/L. kuhli complex was not supported in any of our analyses.
Although close association of L. lacernulatus and B. meridionalis was unanticipated, the distinctiveness of the latter from other Buteogallus species has been recognized by its placement by some authors in the monotypic genus Heterospizias (e.g.,[10, 11, 30–34]). A close relationship between Buteogallus urubitinga, and Harpyhaliaetus eagles has been previously suggested [35, 36], and the retention of the former in Buteogallus has been justified only by the bigger size of the Harpyhaliaetus species . A recent higher-level molecular analysis of hawks and eagles  posited the paraphyly of the genus Buteogallus, and established a closer relationship of Buteogallus urubitinga and Harpyhaliaetus species than between the two Buteogallus species sampled (B. urubitinga and B. anthracinus). Given the polyphyletic relationships of Buteogallus as presented here it is clear that a complete re-assessment of the genus is required. It is worth noting that Amadon , Grossman and Hamlet , Brown and Amadon  and Ridgway  predicted the relatively close relationship of L. schistaceus and L. plumbeus to Buteogallus established in the mtDNA phylogeny presented here. However, we did not recover a sister relationship between the latter pair, sometimes considered subspecies  or members of a super-species complex (e.g., [9, 38, 39]).
The proximity of L. polionotus, L. albicollis and L. occidentalis is fully supported by our data, and has been strongly suggested in taxonomic revisions (e.g., [1, 3, 4, 6, 9, 12, 30, 38, 39]). Species limits in the L. albicollis complex have been a contentious taxonomic issue, and revision of this complex is needed. Similarly, a sister relationship between L. melanops and L. kuhli has been long suggested (e.g.,[1, 4, 6, 30, 38–40]), as well as the close relationship of both to L. semiplumbeus , based on similarities of plumage and external morphology. However, their close relationship to Buteo buteo is novel (see also ). In our study, Buteo buteo represents a genus containing several species of North American and Old World hawks (see ), and the nested position of this species within the Neotropical buteonines corroborates Amadon's  hypothesis of Neotropical origins for the entire buteonine lineage. The polyphyly of Buteo species in our study corroborates the results of Riesing et al. .
Our study includes all recognized Leucopternis species and represents the largest Neotropical buteonine taxonomic sample investigated to date. The polyphyly of Leucopternis, Buteogallus and Buteo indicates that current taxonomy does not reflect the phylogenetic history of the group. Only a complete buteonine phylogeny is likely to provide sufficient guidance regarding the nomenclatural issues concerning Buteo, Leucopternis, Buteogallus and Harpyhaliaetus, as well as fine scale biogeographic inferences. The predominantly black and white plumage pattern shared by several Leucopternis species (L. polionotus, L. occidentalis, L. lacernulatus, L. melanops, L. kuhli, L. semiplumbeus, L. albicollis, L. princeps) has evolved at least twice (figure 2), and the widespread occurrence of this pattern may result from plumage convergence in forested habitats. Similarly, L. plumbeus and L. schistaceus posses a slate plumage pattern very similar to that of Rostrhamus species, and all those species are associated with riparian habitats both within (L. schistaceus, L. plumbeus, R. hamatus) and outside forests (R. sociabilis). In order to test the hypothesis that the slate plumage results from selection in riparian habitats, additional analysis including Rostrhamus species is required. It has been shown that plumage characters under strong selection may evolve rapidly , and in some cases may represent evolutionary convergence instead of reflecting shared phylogenetic or phylogeographic history . An overemphasis on the black and white plumage pattern influenced the grouping of Leucopternis species, and more generally our results indicate that plumage patterns alone may not be reliable taxonomic markers among Accipitridae species.
Biogeography and habitat shifts
We can confidently reject the reciprocal monophyly of cis- and trans- Andean distributed buteonines (figure 2). At least three cis-trans (east-west) disjunctions were identified in all phylogenetic trees inferred from the combined dataset: (1) L. semiplumbeus (trans-) versus its Amazonian (cis-) sister clade, L. melanops plus L. kuhli, (2) L. polionotus (cis-) versus the two sampled subspecies of L. albicollis plus L. occidentalis (trans-), and (3) L. plumbeus (trans-) and the rest of Clade 1 (mostly cis - Andean, but with Buteogallus species occurring on both sides). The uncertain position of L. princeps may hide a possible fourth disjunction. Three major hypothesis have been suggested to explain the differentiation of ancestral populations into cis- and trans- Andean lineages: the Andean Uplift Hypothesis , advocating separation of populations on either side of the mountains as a consequence of Andean orogeny; the Across Andes Dispersal Hypothesis [43–45] proposing long distance dispersal across the Andes as the cause of diversification, and the Forest Refugia Hypothesis, with historical expansions and contractions of forest corridors that linked forested lowlands on either side of the Andes  (see Brumfield and Capparella , Ribas et al. ). Lack of fossil calibrations for raptors, as well as rejection of clock-like evolution for the ATP8 and ATP6 dataset using a likelihood ratio test (data not shown) precluded calculation of divergence times for the hawk species analyzed here, and thus without a temporal framework we are unable to reject any of the Andean biogeography hypotheses based solely on mtDNA phylogenetic inference.
Few geographic barriers besides the Andes seem to define ranges of Neotropical buteonine species. The Amazon River has been suggested to delineate the distributions of terra-firme forest species [48, 49], and may have played a role in the separation of L. melanops and L. kuhli north and south of the river, respectively (see ). However, there are recent and historical records of sympatry between L. melanops and L. kuhli, with observations of both species south of the Amazon River (Barlow et al. , Amaral et al. pers. obs.). Because of this fact it would appear that the Amazon River does not currently impose a strong barrier blocking L. melanops from occupying southern Amazon forests.
The putative contraction and expansion of lowland tropical forests associated with climate change has been promoted as one of the factors that explain avian areas of endemism (see  for limits and details on those areas), but most Neotropical buteonine species have geographic ranges that cross many of the proposed areas. Nonetheless, one might posit that the two Atlantic forest endemics, L. lacernulatus and L. polionotus represent species formed by the contraction and isolation of forest fragments during glacial episodes, but the mtDNA phylogenetic tree establishes that they are not sister species (Figure 1). The distant phylogenetic separation of these two Leucopternis species supports the proposal that the Atlantic Forest biota has complex origins .
Our mtDNA-based phylogenetic analysis of Neotropical buteonines (figure 2) also permits strong inference that forest and open-vegetation species are not reciprocally monophyletic. Furthermore, ecological shifts between forest and non-forest habitats occurred early and late in the Neotropical buteonine diversification. The phylogenetic tree presented in Figure 2 indicates that separation of the open-vegetation Buteo albicaudatus from Leucopternis forest species occurred early in the radiation compared to the more recent divergence of Buteogallus meridionalis from L. lacernulatus and Parabuteo unicinctus from Buteo leucorrhous. The geographic and altitudinal distributions of these more recently derived sister pairs overlap greatly, but the sister species occupy different habitats, either open-country or forest. This same pattern is reflected in the sister relationship between Harpyhaliaetus coronatus (savannah) and Harpyhaliaetus solitarius (forest, unavailable for this study) , although these species do not overlap much in their ranges.
Habitat shifts between sister species or groups of closely related species have been poorly explored in studies of avian historical biogeography in the Neotropics, mainly due to the predominance of forest lowland bird species in such studies (but see Garcia-Moreno and Cardoso da Silva  and Ribas et al. ). The pattern of habitat shifts between sister pairs can be consistent with both allopatric (Theory of Vanishing Refuges ) and parapatric models of speciation (Gradient Hypothesis ). Occurrence of parapatric speciation and divergence with gene flow has been a controversial issue (see Brown , Cracraft and Prum ), but explicit tests with vertebrates [58–61], including birds  have pointed to patterns consistent with between-habitat divergence in the presence of gene flow in tropical habitats. High mobility due to soaring and gliding flight, occurrence of several species in ecotones, and the pattern of sister relationships between forest and non-forest species suggest that buteonine species offer future opportunities to test alternative models of diversification in the Neotropics using phylogeographic data.
The determination of geographically structured areas of endemism in the Neotropics [51, 63] has promoted allopatry as the principal mode of speciation in the Neotropics, with the Andean orogeny, rivers, and changes in forest cover serving as the principal vicariant events separating populations. Accumulation of phylogenetic and phylogeographic data for Neotropical birds suggests that the process of diversification is more complex [e. g. [42, 64]], and that parapatric and sympatric models of speciation must be properly tested [19, 23, 58–62]. The family Accipitridae is extensively represented in most Neotropical habitats, and offers opportunities to explore the radiation of an ecologically diverse group with high dispersal capabilities. Further phylogenetic and phylogeographic studies of diurnal raptors, as well as other groups representing varying degrees of vagility and occurrence in habitats other than forests, will permit more explicit tests of the role of alternative modes of speciation acting on Neotropical birds, and refinement of general explanations for the origin and maintenance of Neotropical biodiversity.
Our mtDNA-based inference of Neotropical buteonine phylogeny establishes a polyphyletic relationship among the hawk genera Leucopternis, Buteogallus and Buteo. Thus the phylogeny indicates that the current taxonomy of the Accipitridae is not a good guide to the evolutionary relationships of species in the group, and identifies a need for further systematic analysis of the family at all taxonomic levels. We do not propose nomenclatural modifications, since only a complete buteonine analysis would permit such taxonomic changes. Nonetheless, our results coupled to earlier work predict some of the nomenclatural changes that will undoubtedly be forthcoming, and also establish that plumage has been overemphasized in defining the taxonomy of the Accipitridae. Finally, we conclude that shifts between forest and non-forest habitats, as well as movement across the Andes, have occurred more than once during the Neotropical buteonine diversification.
Taxon sampling, DNA extraction, amplification and sequencing
We sampled a total of 31 specimens, comprising all 10 recognized species of the genus Leucopternis (L. lacernulatus, L. polionotus, L. semiplumbeus, L. plumbeus, L. occidentalis, L. schistaceus, L. princeps, L. melanops, L. kuhli, L. albicollis), and including two sub-species of L. albicollis (L. a. costaricensis and L. a. ghiesbreghti), as well as almost all Neotropical buteonine genera (representing 15 of the 21 species of the Neotropical "sub-buteonines" sensu Amadon ), plus four Buteo species. We chose Geranospiza caerulescens as an outgroup based on a recent higher-level analysis of Accipitridae . When possible, we included two individuals per taxon. Nomenclature follows the South American Classification Committee of the American Ornithologists' Union . Sequences of Buteo buteo were obtained from Genbank [NC_003128]. Tissue, feather and blood samples were obtained from specimens collected in the field, museum tissue collections, and captive birds (see Additional file 1: Table 1). Known localities of origin, feathers and photographs are available for most captive specimens.
DNA extraction, amplification and sequencing were performed at the Universidade de São Paulo (Brazil), the Smithsonian Tropical Institute (Panama) and the Royal Ontario Museum (Canada) based on earlier protocols . DNA extraction followed Bruford et al. , or via the DNeasy kit (Qiagen); for feather samples, we added 30 ug of dithiothreitol to the digestion buffer. We sequenced four mitochondrial genes: a portion of the 12S ribosomal RNA gene (12S, longest sequence of 833 bp), the complete ATP synthase F0 subunit 8 (ATP8, 168 bp) and subunit 6 (ATP6, 684 bp) genes, as well as the complete NADH dehydrogenase subunit 6 (ND6, 519 bp) using several primer pair combinations via polymerase chain reaction (PCR) (table 2). In few cases, weak amplification products were re-amplified using internal primers. Both strands of the amplified products were sequenced.
Alignment and phylogenetic analysis
Multiple strands obtained for each specimen were assembled in CodonCode Aligner v. 1.3.4 (CodonCode Corporation) or Sequencher v. 4 (Gene Codes Corporation). Contigs were exported and alignment performed in Clustal X 1.83  with default parameters. The 12S alignment had 21 indels, which consisted mostly of autapomorphies and sites of ambiguous alignment; these were removed from all analyses. All single marker and combined datasets were tested for significant departures from average base frequencies with PAUP* 4b10 , using only variable sites. Uncorrected codon-based (ATP8, ATP6, ND6) and total (12S) transition and transversion distances were plotted against Kimura-2-parameters distances using the software Dambe v4.2  to evaluate the effect of multiple substitutions in each dataset. We implemented a partition homogeneity test in PAUP*, using only variable sites with 1000 replicates with random additions, to evaluate the congruence of the phylogenetic signal between the different genes. Because the latter test did not detect significantly-different phylogenetic signal among the partitions (P = 0.29), all four genes were combined into a single combined dataset. Separated analyses of subsets of single gene region were also performed to evaluate the concordance among those datasets. The ATP8 and ATP6 genes overlap by 10 bp; however in all phylogenetic analyses of the combined dataset and the single gene region of ATP (subunit 8 plus subunit 6) this region was considered only once. We performed phylogenetic reconstructions using maximum likelihood (ML) and maximum parsimony (MP) implemented in PAUP*, and Bayesian analysis (BA) with MrBayes v3.1.1 , to evaluate the concordance of topologies obtained under different optimization criteria. ML and MP heuristic searches were performed using 1000 and 10 random additions of sequences, respectively. Nonparametric bootstrapping was performed to assess branch support (100 replicates with single random additions for the ML analysis and 1,000 replicates with 10 random additions for the MP analysis). Modeltest v3.7  was used to choose among evolutionary models of DNA substitution for ML and BA analyses using a hierarchical likelihood ratio test. Modeltest determined that the TrN+I+G model was the best fit for the total dataset with base frequencies of A = 0.3217, C = 0.3521, G = 0.1235, T = 0.2027, a gamma shape parameter of 1.0617 and proportion of invariable sites of 0.5631. The BA analyses of the combined dataset was run with individual likelihood for each of the three gene regions (12S, ATP 6 and ATP 8, and ND6) as selected by Modeltest (TrN+I+G, HKY+I+G, TrN+G), which were the same models used for subset ML and BA analyses. MrBayes was run with four chains for 4,000,000 generations with trees sampled every 100 generations, replicated four times. All runs reached stationarity around 400 sampled generations, so we discarded the first 40000 generations as a "burnin"; a consensus topology was created with all the remaining sampled generations.
To determine whether our data support monophyly of Leucopternis, monophyly of forest species, monophyly of species presenting black and white plumage patterns, or reciprocal monophyly of trans- and cis- Andean species, we compared alternative constraint topologies to the ML tree using the nonparametric Shimodaira-Hasegawa test  implemented in PAUP*. Ancestral states of habitat and general plumage pattern (according to Thiollay , Fergusson-Lees and Christie  and Sibley and Monroe ) were mapped onto the ML tree inferred from the combined dataset using unordered parsimony in Mesquite v1.05 .
Thiollay JM: Family Accipitridae (Hawks and Eagles). Handbook of the birds of the world. New World Vultures to Guineafowl. Edited by: Hoyo J del, Elliott A, Sargatal J. 1994, Barcelona: Lynx Edicions, 2: 52-205.
Sick H: Ornitologia Brasileira. 1997, Rio de Janeiro: Nova Fronteira
Amadon D: Taxonomic notes on birds of prey. Amer Mus Novit. 1964, 2166: 1-24.
Amadon D: A revision of the sub-buteonine hawks (Accipitridae, Aves). Amer Mus Novit. 1982, 2741: 1-20.
Olson SL: The fossil record of birds. Avian biology. Edited by: Farner D, King J R, Parkes K. 1985, New York:Academic Press, VIII: 79-239.
Ferguson-Lees J, Christie DA: Raptors of the world. 2001, London: Christopher Helm
Riesing MJ, Kruckenhauser L, Gamauf A, Haring E: Molecular phylogeny of the genus Buteo (Aves: Accipitridae) based on mitochondrial marker sequences. Mol Phylogenet Evol. 2003, 27: 328-342. 10.1016/S1055-7903(02)00450-5.
Helbig AJ, Kocum A, Seibold I, Braun MJ: A multi-gene phylogeny of aquiline eagles (Aves: Accipitriformes) reveals extensive paraphyly at the genus level. Mol Phyl Evol. 2005, 35: 147-164. 10.1016/j.ympev.2004.10.003.
Amadon D, Bull J: Hawks and owls of the world. A distributional and taxonomic list. Proc West Found Vert Zool. 1988, 3: 295-357.
Friedman H: The birds of North and Middle America. Part XI. Cathartidae to Falconidae. 1950, United States National Museum Bulletin 50. Washington D.C.: Smithsonian Institution, 50:
Grossman ML, Hamlet J: Birds of prey of the world. 1964, New York: Bonanza Books
Brown LH, Amadon D: Eagles, hawks and falcons of the world. 1968, Feltham: Country Life Books
Haffer J: Speciation in Amazonian forest birds. Science. 1969, 165: 131-137.
Haffer J: Time's cycle and time's arrow in the history of Amazonia. Biogeographica. 1993, 69: 15-45.
Vanzolini PE, Williams EE: South American anoles: the geographical variation and and evolution of the Anolis chrysolepis species group (Sauria, Iguanidae). Arq Zool São Paulo. 1970, 19: 1-298.
Wallace AR: On the monkeys of the Amazon. Proc Zool Soc London. 1852, 20: 107-110.
Salo J, Kalliola R, Häkkinen I, Mäkinen Y, Niemelä P, Pubakka M, Coley PD: River dynamics and the diversity of Amazon lowland forest. Nature. 1986, 322: 254-258. 10.1038/322254a0.
Räsänen ME, Linna AM, Santos JCR, Negri FR: Late Miocene tidal deposits in the Amazonian foreland basin. Science. 1995, 269: 386-390.
Marroig G, Cerqueira R: Plio-Pleistocene South American history and the Amazon lagoon hypothesis: a piece in the puzzle of Amazonian diversification. J Comp Biol. 1997, 2: 103-119.
Nores M: An alternative hypothesis for the origin of Amazonian bird diversity. J Biogeography. 1999, 26: 475-485. 10.1046/j.1365-2699.1999.t01-1-00311.x.
Platnick NI, Nelson GJ: A method of analysis for historical biogeography. Syst Zool. 1978, 27: 1-16.
Endler JA: Pleistocene forest refuges: fact or fancy. Biological Diversification in the Tropics. Edited by: Prance GT. 1982, New York: Columbia University Press, 641-657.
Moritz C, Patton JL, Schneider CJ, Smith TB: Diversification of rainforest faunas: an integrated molecular approach. Annu Rev Ecol Syst. 2000, 31: 533-563. 10.1146/annurev.ecolsys.31.1.533.
Dick CW, Abdul-Salim K, Bermingham E: Molecular systematic analysis reveals cryptic Tertiary diversification of a widespread tropical rain forest tree. Am Nat. 2003, 162: 691-703. 10.1086/379795.
Dick CW, Roubik DW, Gruber KF, Bermingham E: Long-distance gene flow and cross-Andean dispersal of lowland rainforest bees (Apidae: Euglossini) revealed by comparative mitochondrial DNA phylogeography. Mol Ecol. 2004, 13: 3775-3885. 10.1111/j.1365-294X.2004.02374.x.
Bildstein KL: Raptor migration in the Neotropics: patterns, processes and consequences. Ornitol Neotrop. 2004, 83-99. Suppl 15
Holdaway R: An exploratory phylogenetic analysis of the genera of the Accipitridae, with notes on the biogeography of the family. Raptor conservation today: proceedings of the IV World Conference on Birds of Prey and Owls. 10–17 May 1992, Berlin. Edited by: Meyburg BU, Chancellor RD. 1994, Berlin, London, Paris: Pica Press and World Working Group on Birds of Prey and Owls, 601-637.
Lerner HRL, Mindell DP: Phylogeny of eagles, Old World vultures, and other Accipitridae based on nuclear and mitochondrial DNA. Mol Phylogenet Evol. 2005, 37: 327-346. 10.1016/j.ympev.2005.04.010.
Pereira SL, Baker AJ: Low number of mitochondrial pseudogenes in the chicken (Gallus gallus) nuclear genome: implications for molecular inference of population history and phylogenetics. BMC Evol Biol. 2004, 4: 17-10.1186/1471-2148-4-17.
Hellmayr CE, Conover B: Catalogue of Birds of the Americas. Field Mus Nat Hist Zool Ser. 1949, 13-
Pinto OMO: Catálogo das aves do Brasil e lista dos exemplares que as representam no Museu Paulista, 1.a Parte. Aves não Passeriformes e Passeriformes não Oscines excluida a Fam. Tyrannidae e seguintes. Rev Mus Paulista. 1938, 22: 1-566.
Plótnick R: Posicíon sistemática del género Heterospizias. Hornero. 1956, 10: 136-139.
Phelps WH, Phelps WH: Lista de las aves de Venezuela con su distribucion. Tomo II. parte I. No Passeriformes. Bol Soc Venez Cienc Nat. 1958, 19: 1-317.
Meyer de Schauensee R: A guide to the birds of South America. 1970, Wynnewood, Pennsylvania: Livingston Publishing
Amadon D: Notes on Harpyhaliaetus. Auk. 1949, 66: 53-56.
Amadon D, Eckelberry DR: Observations of Mexican Birds. Condor. 1955, 57: 65-80.
Ridgway R: Studies of the American Falconidae. Bull U S Geol Geog Surv Terr. 1876, 2: 91-182.
Stressemann E, Amadon D: Order Falconiformes. Check-List of Birds of the World. Edited by: Mayr E, Cottrell GW. 1979, Cambridge, Massachusetts: Museum of Comparative Zoology, 1: 271-425. 2
Sibley CG, Monroe BL: Distribution and taxonomy of birds of the world. 1990, New Haven and London: Yale Univ Press
Haffer J: Biogeography of Neotropical birds. Biogeography and quaternary history in tropical America. Edited by: Whitmore TC, Prance GT. 1987, Oxford:Clarendon Press
Omland KE, Lanyon SM: Reconstructing plumage evolution in orioles (Icterus): repeated convergence and reversal in patterns. Evolution. 2000, 54: 2119-2133.
Cheviron ZA, Hackett SJ, Capparella AP: Complex evolutionary history of a Neotropical lowland forest bird (Lepidothrix coronata) and its implications for historical hypotheses of the origin of Neotropical avian diversity. Mol Phylogenet Evol. 2005, 36: 338-357. 10.1016/j.ympev.2005.01.015.
Chapman FM: The distribution of bird life in Colombia. Bull Am Mus Nat Hist. 1917, 36: 1-729.
Chapman FM: The distribution of bird life in Colombia. Bull Am Mus Nat Hist. 1926, 55: 1-784.
Haffer J: Speciation in Colombian forest birds west of the Andes. American Museum Novitates. 1967, 294: 1-57.
Brumfiled RT, Capparella AP: Historical diversification of birds in northwestern South America: a molecular perspective on the role of vicariant events. Evolution. 1996, 1607-1624.
Ribas CC, Gaban-Lima R, Miyaki CY, Cracraft J: Historical biogeography and diversification within the Neotropical parrot genus Pionopsitta (Aves:Psittacidae). J Biogeography. 2005, 32: 1409-1427.
Capparella AP: Genetic variation in Neotropical birds: implications for the speciation process. Acta XIX Congr Int Ornithol. 1988, 19: 1658-1673.
Aleixo A: Historical diversification of a Terra-firme forest bird superspecies: a phylogeographic perspective on the role of different hypotheses of Amazonian diversification. Evolution. 2004, 58: 1303-1317.
Barlow J, Haugaasen T, Peres C: Sympatry of the Black-faced Hawk Leucopternis melanops and the White-browed Hawk Leucopternis kuhli in the Lower Rio Tapajós, Pará, Brazil. Cotinga. 2002, 18: 77-79.
Cracraft J: Historical biogeography and patterns of differentiation within the South American avifauna: areas of endemism. Ornith Monogr. 1985, 36: 49-84.
Costa L: The historical bridge between Amazon and the Atlantic Forest of Brazil: a study of molecular phylogeography with small mammals. J Biogeography. 2003, 30: 71-86. 10.1046/j.1365-2699.2003.00792.x.
Garcia-Moreno J, Cardoso da Silva JM: An interplay between forest and non-forest South American avifaunas suggested by a phylogeny of Lepidocolaptes woodcreepers (Dendrocolaptidae). Stud Neotrop Fauna & Environm. 1997, 32: 164-173.
Ribas CC, Miyaki CY: Molecular systematics in Aratinga parakeets: species limits and historical biogeography in the 'solstitialis' group, and the systematic position of Nandayus nenday. Mol Phylogenet Evol. 2004, 30: 663-675. 10.1016/S1055-7903(03)00223-9.
Vanzolini PE, Williams EE: The vanishing refuge: a mechanism for ecogeographic speciation. Pap Avulsos Zoo. 1981, 34: 251-255.
Brown KS: Conclusions, synthesis, and alternative hypothesis. Biogeography and quaternary history in tropical America. Edited by: Whitmore TC, Prance GT. 1987, New York: Clarendon Press, 175-196.
Cracraft J, Prum RO: Patterns and processes of diversification: speciation and historical congruence in some Neotropical birds. Evolution. 1988, 42: 603-620.
Schneider CJ, Smith TB, Larison B, Moritz C: A test of alternative models of diversification in tropical rainforests: ecological gradients vs rainforest refugia. Proc Nat Acad Sci USA. 1999, 24: 13869-13873. 10.1073/pnas.96.24.13869.
Smith TB, Wayne RK, Girman DJ, Bruford MW: A role for ecotones in generating rainforest biodiversity. Science. 1997, 276: 1855-1857. 10.1126/science.276.5320.1855.
Smith TB, Schneider CJ, Holder K: Refugial isolation versus ecological gradients. Genetica. 2001, 112–113: 383-398. 10.1023/A:1013312510860.
Ogden R, Thorpe RS: Molecular evidence for ecological speciation in tropical habitats. Proc Natl Acad Sci USA. 2002, 99: 13612-13615. 10.1073/pnas.212248499.
Smith TB, Calsbeek R, Wayne RK, Holder KH, Pires D, Bardeleben C: Testing alternative mechanisms of evolutionary divergence in an African rain forest passerine bird. J Evol Biol. 2005, 18: 257-268. 10.1111/j.1420-9101.2004.00825.x.
Bates JM, Hackett SJ, Cracraft J: Area-relationships in the Neotropical lowlands: a hypothesis based on raw distributions of Passerine birds. J Biogeography. 1998, 25: 783-793. 10.1046/j.1365-2699.1998.2540783.x.
Marks BD, Hackett SJ, Capparella AP: Historical relationships among Neotropical lowland forest areas of endemism as determined by mitochondrial DNA sequence variation within the Wedge-billed Woodcreeper (Aves: Dendrocolaptidae: Glyphorynchus spirurus). Mol Phylogenet Evol. 2002, 24: 153-167. 10.1016/S1055-7903(02)00233-6.
Remsen JV, Jaramillo A, Nores M, Pacheco JF, Robbins MB, Schulenberg TS, Stiles FG, da Silva JMC, Stotz DF, Zimmer KJ: A classification of the bird species of South America. American Ornithologists' Union. [Version 30 november 2005], [http://www.museum.lsu.edu/~Remsen/SACCBaseline.html]
Tavares ES, Yamashita C, Miyaki CY: Phylogenetic relationships among some Neotropical parrot genera (Psittacidae) based on mitochondrial sequences. Auk. 2004, 121: 230-242.
Bruford MW, Hanotte O, Brookfield JFY, Burke T: Multi and single-locus DNA fingerprinting. Molecular Genetic Analysis of Populations – a Practical Approach. Edited by: Hoelzel AR. 1992, New York: IRL Press, 225-269.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.
Swofford DL: PAUP* Phylogenetic analysis using parsimony (*and other methods). 2002, Sunderland, Massachusetts: Sinauer Associates
Xia X, Xie Z: DAMBE: Data analysis in molecular biology and evolution. Journal of Heredity. 2001, 92: 371-373. 10.1093/jhered/92.4.371.
Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Posada D, Crandall KA: MODELTEST: testing the model of nucleotide substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116.
Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 1.05. 1999, [http://mesquiteproject.org]
Miyaki CY, Matioli S, Burke T, Wajntal A: Parrot evolution and paleogeographic events: mitochondrial DNA evidence. Mol Biol Evol. 1998, 15: 544-551.
Eberhard JR, Wright TF, Bermingham E: Duplication and concerted evolution of the mitochondrial control region in the parrot genus Amazona. Mol Biol Evol. 2001, 18: 1330-1342.
Sorenson MD, Ast JC, Dimcheff DE, Yuri T, Mindell DP: Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates. Mol Phylogenet Evol. 1999, 12: 105-114. 10.1006/mpev.1998.0602.
Haring E, Kruckenhauser L, Gamauf A, Riesing MJ, Pinsker W: The complete sequence of the mitochondrial genome of Buteo buteo (Aves, Accipitridae) indicates an early split in the phylogeny of raptors. Mol Biol Evol. 2001, 18: 1892-1904.
Eberhard JR, Bermingham E: Phylogeny and biogeography of the Amazona ocrocephala (Aves:Psittacidae) complex. Auk. 2004, 121: 318-332.
Fleischer RC, Olson SL, James HF, Cooper AC: Identification of the extinct Hawaiian eagle (Haliaeetus) by mtDNA sequence analysis. Auk. 2000, 117: 1051-1056.
We thank the following individuals and institutions for kindly providing samples, without which this work could not have been performed: Donna Dittman, Robb Brumfield and Fred Sheldon (LSUMZ); Leo Joseph (ANSP); Alexandre Aleixo and Marcos Pérsio D. Santos (MPEG); Marcelo Soares (UFPA); Carlos E. A. Carvalho; Eduardo P. M. de Carvalho Filho; Giancarlo Zorzin, Gustavo D. M. de Carvalho and Marcus Canuto (SOS Falconiformes), Pedro Scherer Neto (Museu de História Natural Capão do Imbuía), Zoológico de Paulínia, Ricardo Pereira (Zoológico de São Paulo), Leo Fukui and Jorge Lisboa (ABPFAR); David Whitacre, Rick Watson and Martin Gilbert (Peregrine Fund); Tammo Hoeksema (ZOOMAT); Sérgio Aguilar (Veracruz Pronatura); Alexandre Miranda, José H. Fontenele, Greicelene R. Pedro (Orquidário de Santos); Renato Gaban-Lima and Guilherme Renzo (USP); Adriana Joppert (DEPAVE/SP); Robson Silva e Silva and Fábio Olmos. We are grateful to Erika S. Tavares for kindly sequencing the ANSP samples at the Royal Ontario Museum and helping with the analyses, and Allan Baker for allowing use of the ROM laboratory for sequencing. We also thank Sérgio Pereira and Martin Riesing for useful suggestions and support, Erwin T. Grau for invaluable help with the laboratory methods, Robb Brumfield for help with the analyses, and Camila Ribas, Cibele Biondo, Gustavo S. Cabanne, Rodrigo Pessoa and three anonymous referees for comments and several contributions to the manuscript. Financial support was provided by Fundação de Amparo a Pesquisa do Estado de São Paulo, CNPq and CAPES.
FSRA conceived the study, carried out most of the data collection and phylogenetic analysis, and drafted the manuscript. MJM carried out the data collection from the LSUMZ samples, and made substantial contributions to the manuscript. LFS obtained part of the samples, and made substantial contributions to the manuscript. EB made substantial contributions to the manuscript. AW helped to conceive the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.