Multidirectional chromosome painting substantiates the occurrence of extensive genomic reshuffling within Accipitriformes

Background Previous cross-species painting studies with probes from chicken (Gallus gallus) chromosomes 1–10 and a paint pool of nineteen microchromosomes have revealed that the drastic karyotypic reorganization in Accipitridae is due to extensive synteny disruptions and associations. However, the number of synteny association events and identities of microchromosomes involved in such synteny associations remain undefined, due to the lack of paint probes derived from individual chicken microchromosomes. Moreover, no genome-wide homology map between Accipitridae species and other avian species with atypical karyotype organization has been reported till now, and the karyotype evolution within Accipitriformes remains unclear. Results To delineate the synteny-conserved segments in Accipitridae, a set of painting probes for the griffon vulture, Gyps fulvus (2n = 66) was generated from flow-sorted chromosomes. Together with previous generated probes from the stone curlew, Burhinus oedicnemus (2n = 42), a Charadriiformes species with atypical karyotype organization, we conducted multidirectional chromosome painting, including reciprocal chromosome painting between B. oedicnemus and G. fulvus and cross-species chromosome painting between B. oedicnemus and two accipitrid species (the Himalayan griffon, G. himalayensis 2n = 66, and the common buzzard, Buteo buteo, 2n = 68). In doing so, genome-wide homology maps between B. oedicnemus and three Accipitridae species were established. From there, a cladistic analysis using chromosomal characters and mapping of chromosomal changes on a consensus molecular phylogeny were conducted in order to search for cytogenetic signatures for different lineages within Accipitriformes. Conclusion Our study confirmed that the genomes of the diurnal birds of prey, especially the genomes of species in Accipitriformes excluding Cathartidae, have been extensively reshuffled when compared to other bird lineages. The chromosomal rearrangements involved include both fusions and fissions. Our chromosome painting data indicated that the Palearctic common buzzard (BBU) shared several common chromosomal rearrangements with some Old World vultures, and was found to be more closely related to other Accipitridae than to Neotropical buteonine raptors from the karyotypic perspective. Using both a chromosome-based cladistic analysis as well as by mapping of chromosomal differences onto a molecular-based phylogenetic tree, we revealed a number of potential cytogenetic signatures that support the clade of Pandionidae (PHA) + Accipitridae. In addition, our cladistic analysis using chromosomal characters appears to support the placement of osprey (PHA) in Accipitridae. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0484-0) contains supplementary material, which is available to authorized users.


(Continued from previous page)
Conclusion: Our study confirmed that the genomes of the diurnal birds of prey, especially the genomes of species in Accipitriformes excluding Cathartidae, have been extensively reshuffled when compared to other bird lineages. The chromosomal rearrangements involved include both fusions and fissions. Our chromosome painting data indicated that the Palearctic common buzzard (BBU) shared several common chromosomal rearrangements with some Old World vultures, and was found to be more closely related to other Accipitridae than to Neotropical buteonine raptors from the karyotypic perspective. Using both a chromosome-based cladistic analysis as well as by mapping of chromosomal differences onto a molecular-based phylogenetic tree, we revealed a number of potential cytogenetic signatures that support the clade of Pandionidae (PHA) + Accipitridae. In addition, our cladistic analysis using chromosomal characters appears to support the placement of osprey (PHA) in Accipitridae.

Background
Most birds studied so far have highly conserved karyotypes, similar to that of the chicken (Gallus gallus, GGA, 2n = 78), consisting of several pairs of macrochromosomes and a mass of microchromosomes, with diploid chromosome number (2n) varying mainly from 76 to 84 [1,2]. In contrast, the diurnal birds of prey have a unique karyotype organization that differs widely from the usually conserved genomic structure found in other bird lineages [3,4]. In particular, Accipitridae species display the most atypical karyotypes known in Aves [5]. The more than fifty accipitrid species studied so far (reviewed in [6]) share the following karyotypic characteristics: 1) most species have diploid chromosome numbers varying from 66 to 68; 2) they lack large macrochromosomes; and 3) they have many medium-to small-sized bi-armed chromosomes and usually three to six pairs of microchromosomes [3,5]. Only one pair of microchromosomes was found in the black-winged kite (Elanus caeruleus, 2n = 68) [7]. Previous karyotypic comparisons have suggested that such atypical karyotypes of accipitrids could have evolved from a typical bird karyotype by a series of fissions of macrochromosomes and translocations of some macrochromosome segments onto the microchromosomes and small acrocentrics [3]. However, the exact nature of chromosome structural rearrangements that took place in the karyotype evolution of the accipitrid lineage remains largely unclear.
Most recently, two complete sets of painting probes have been generated from two neoavian species: the stone curlew (Burhinus oedicnemus, BOE, 2n = 42, Charadriiformes) [14] and the white hawk (P. albicollis, PAL, 2n = 66, Accipitriformes) [9]. BOE has one of the lowest diploid numbers reported so far in birds. Using the GGA and BOE chromosome painting probes, the first reciprocal chromosome painting between avian species was performed, and a comparative chromosome map between BOE and GGA has been established [14]. The whole set of BOE painting probes was subsequently used to delineate chromosome homology between BOE and representatives of five avian orders [23]. In addition, PAL probes were used to detect chromosomal homologies between GGA 1-10 and that of PAL [9], between PAL and three macrochromosomes of a New World vulture, the turkey vulture (Cathartes aura, 2n = 80) [15], three Buteoninae species (2n = 68) [16] and two Turdidae songbirds (2n = 78, Passeriformes) [17]. Up to now, all cross-species chromosome painting studies in birds have involved either two species with typical avian karyotypes, or one species with a typical and the other with an atypical karyotype. No genome-wide comparative chromosome painting between two bird species with atypical karyotypes has ever been performed.
In the past, the diurnal birds of prey constituted the Order Falconiformes and were classified into five families (i.e. Cathartidae, Falconidae, Accipitridae, and the two monotypic families Pandionidae and Sagittariidae) [24]. However, the status of Falconiformes as a natural (i.e. monophyletic) group was long under debate due to marked heterogeneity in karyological, morphological and mitochondrial (mt) gene order found in different families of Falconiformes (reviewed in [25]). The most recent molecular phylogeny studies retrieved strong support for a clade containing the secretary bird (Sagittariidae), the osprey (Pandionidae) and the traditional accipitrids (e.g., [26][27][28][29][30][31][32][33]). Consequently, Falconiformes now only includes the family Falconidae, while the other diurnal birds of prey constitute a new order, Accipitriformes. Here, we adopted this new taxonomic classification for the diurnal birds of prey. In Accipitriformes, there are some important unresolved issues concerning the rank of the osprey and the status of the Old World vultures as well as the Buteoninae as monophyletic or polyphyletic groups. The osprey was often recognized as a family-level taxa (e.g. [26][27][28][29]32]); however, in some studies, together with the secretary bird, the osprey was assigned the rank of subfamily within Accipitridae (e.g. [30]). Recent molecular studies suggested that the Old World vultures and the Buteoninae were polyphyletic [26,27,29,31,32]. Up to now, no Sagittariidae species (i.e. the secretary bird) has ever been studied using chromosome painting. Besides Accipitridae, three species from other families of Accipitriformes and three species from Falconiformes have been studied by this technique [12,15,18,34]. The availability of chromosome painting data from different diurnal birds of prey families offers us the opportunity to decipher the genomic rearrangements within the diurnal birds of prey and to search for cytogenetic evidence for the phylogenetic relationships within the Accipitriformes.
In this context, we here report the generation of 1) a complete set of painting probes that cover the whole genome of the griffon vulture (G. fulvus, GFU, 2n = 66), an Old World vulture, through bivariate flow sorting and DOP-PCR amplification; and 2) comparative chromosome maps based on reciprocal chromosome painting between BOE and GFU, two species belonging to different orders with atypical karyotypes (namely, Charadriiformes and Accipitriformes). Moreover, the whole set of BOE probes was hybridized to metaphase chromosomes of two species in Accipitridae, Himalayan griffon (G. himalayensis, GHI, 2n = 66) and common buzzard (Buteo buteo, BBU, 2n = 68) to detect karyotypic differences between these species. A cladistic analysis using chromosomal rearrangements as well as the mapping of ZOO-FISH-inferred chromosomal differences on a concensus molecular phylogeny were performed, providing further insight into the origin and evolution of the extensively rearranged Accipitriformes genome.

Results
Characterization of the G. fulvus flow karyotype As reported previously [3], GFU has a 2n = 66 karyotype, including fifteen meta-or sub-metacentric chromosomes, thirteen acrocentric chromosomes, four pairs of dot-like microchromosomes, together with a large submetacentric Fig. 1 Bivariate flow karyotype of G. fulvus (2n = 66) with chromosome assignments Z and a medium-sized submetacentric W. The thirtythree pairs of GFU chromosomes were resolved into twenty-three flow peaks (Fig. 1). Hybridizing paints prepared from individual peaks onto GFU metaphases established the identity of the chromosomes contained in each peak. Reciprocal chromosome painting between GFU and BOE was also performed to resolve any ambiguity in the chromosomal assignment of each flow peak and, in particular, of the peaks that contained multiple chromosomes (see below). The majority of flow peaks each contained only one type of GFU chromosome, with the exception of seven peaks, of which four peaks each contained two GFU chromosomes (GFU 2 + Z, 13 + 20, 14 + 27 and 24 + 25), one peak contained three GFU chromosomes (GFU 11 + 19 + 26) and the remaining two peaks each contained four GFU chromosomes (GFU 5 + 8 + 9 + 10 and 29 + 30 + 31 + 32). Moreover, the two homologues of chromosome 24 were each sorted into a different flow peak. Thus, a set of chromosome painting probes that cover the entire genome of GFU was generated.
Reciprocal chromosome painting between G. fulvus and B. oedicnemus Together with BOE [14] and PAL [9], the set of GFU painting probes represents the third set of paints that have been generated from Neoaves. Reciprocal chromosome painting between GFU and BOE established chromosome homologies between these species and de-  Fig. 2 Reciprocal chromosome painting between BOE and GFU. a BOE 1 probe hybridized to seven pairs of GFU chromosomes. b Probes from GFU 23 (red), 26 (+11 + 19) (green), 7 (red), 22 (green) and 21 (red) painted BOE 1. c BOE 2 probe hybridized to three pairs of GFU chromosomes. d Probes from GFU 1 (red), 24 (+25) (green) and 3 (red) labeled BOE 2. e BOE 3 probe hybridized to four pairs of GFU chromosomes. f Probes from GFU 16 (red), 8 (+5 + 9 + 10) (green) and 19 (green) probes painted BOE 3. g GFU 16 probe labeled GFU16. h GFU 16 probe hybridized to BOE 3 and 15. Note that cross-hybridization signals were also detected on BOE Z and W. Due to the fact that one of the used GFU probes in (b), (d) and (f) contained two or more GFU chromosomes or that one used GFU probe gave signals on two or more BOE chromosomes, besides BOE 1-3, other BOE chromosomes also were painted in (b), (d) and (f) Fig. 3 a DAPI-banded karyotype of GFU with the assignment of homologies to BOE. b DAPI-banded karyotype of BOE with the assignment of homologies to GFU on the right. Homologies to GGA are indicated on the left of each GFU and GHI and BBU chromosome to cross-hybridization signals on the long arm of GFU W (Fig. 2a, c and e). BOE 5 (= GGA 8 + 7)-12 (= two GGA microchromosomes) probes each hybridized to two different GFU chromosomes or chromosome segments (Fig. 4a, b and f ). BOE 4 (= GGA 4q), 13 (= two GGA microchromosomes) and 14 (= two GGA microchromosomes) probes each painted one GFU chromosome or chromosome segment (Fig. 4b). BOE 15 + 16 (= two GGA microchromosomes) probe gave signals on three pairs of GFU chromosomes. A paint derived from four microchromosomes (i.e., BOE 17 + 18 + 19 + 20) hybridized to the four pairs of GFU microchromosomes (29, 30, 31 and 32). In total, BOE autosomal probes revealed forty homologous chromosomal segments in the GFU genome (Fig. 3a).
GFU autosomal probes revealed forty-one segments of conserved synteny in the BOE genome (Fig. 3b). Twenty-four GFU chromosomes (1, 2, 6, 7, 9, 12-14 and 17-32) each corresponded to one BOE chromosome or chromosome segment; eight GFU chromosomes (3-5, 8, 10, 11, 15 and 16) each corresponded to two BOE chromosomes or chromosome segments. In agreement with the finding that BOE 1-3 (= GGA 1-3) probes gave signals on multiple GFU chromosomes, BOE 1-3 (= GGA (See figure on previous page.) Fig. 4 Examples of cross-species chromosome painting with BOE painting probes. a Hybridization of BOE 5 (green) and 11 (red) probes to GFU chromosomes. b Hybridization of BOE 12 (green) and 13 (red) probes to GFU chromosomes. c Hybridization of BOE 1 probe to seven pairs of GHI chromosomes. d Hybridization of BOE 1 probe to eight pairs of BBU chromosomes. e Hybridization of BOE 1 (red) and 9 (green) probes to BBU chromosomes. f Hybridization of BOE 1 (red) and 9 (green) probes to GFU chromosomes. g Hybridization of BOE 10 probe to three pairs of BBU chromosomes. h Hybridization of BOE 10 probe to two pairs of GHI chromosomes. Note that cross-hybridization signal was also detected on the W chromosomes of GFU (a, b, f) and GHI (c) Fig. 5 a DAPI-banded karyotype of GHI with assignment of homologies to BOE. b DAPI-banded karyotype of BBU with assignment of homologies to BOE. Homologies to GGA are indicated on the left of each BOE chromosome 1-3) were each painted by three or more GFU probes (Fig. 2b, d and f ). The GFU 16 probe (Fig. 2g) gave two discrete signals on BOE 3 and one signal on BOE 15 (Fig. 2h). Cross-hybridization signals were also found on BOE Z and W chromosomes (Fig. 2h). Finally, the pooled probe from the four smallest GFU microchromosomes (29)(30)(31)(32) painted the smallest four pairs of BOE microchromosomes (17)(18)(19)(20).
Hybridizing B. oedicnemus probes onto metaphases of G. himalayensis and B. buteo Like other Gyps species, G. himalayensis (GHI) has a 2n = 66 karyotype, which is similar to that of GFU. As expected, the patterns produced by hybridizing each BOE probe to GFU and GHI chromosomes were identical ( Fig. 3a and Fig. 5a). Examples of BOE probes hybridizing to chromosomes of GHI are shown in Fig. 4c and h. Accordingly, paints from BOE autosomes detected forty homologous segments in GHI (Fig. 5a).
The male B. buteo (BBU) shows a 2n = 68 karyotype, comprising fifteen pairs of meta-or sub-metacentric chromosomes, fourteen pairs of acrocentric chromosomes, four pairs of microchromosomes and one pair of large submetacentric Z chromosomes (Fig. 5b), as previously reported [5]. Painting probes from BOE were hybridized to metaphase spreads of BBU. FISH examples are shown in Fig. 4d, e and g. Results are summarized using a DAPI-banded karyotype of BBU (Fig. 5b). Hybridization patterns produced by each BOE probe were very similar in GFU, GHI and BBU. Only two BOE probes gave different numbers of signals in BBU, GFU and GHI. BOE 1 (= GGA 1) probe painted eight pairs of BBU chromosomes ( Fig. 4d and e), but only seven pairs of chromosomes in GFU and GHI (Figs. 2a and 4c, f ); BOE 10 (= GGA 10 + 1 microchromosome) probe hybridized onto three pairs of BBU chromosomes ( Fig. 4g) but two pairs of chromosomes in GFU and GHI (Fig. 4h). In total, probes from BOE autosomes revealed forty-two homologous segments in the BBU genome.
To facilitate comparison of homologous chromosomal segments across bird species, the homologous GGA chromosome segments as inferred from the GGA-BOE comparative chromosomal map [14] were also indicated to the left of GFU, BOE, GHI and BBU chromosomes (Figs. 3 and 5).

Chromosomal phylogeny of the diurnal birds of prey
In total, 59 chromosomal characters were constructed from available inter-species chromosome painting data including 14 species of Accipitriformes and 3 species of Falconiformes (Table 1). The exhaustive search allowed us to retrieve 10 equally most parsimonious trees (L = 67; Additional file 1: Figure S1). These 10 trees were largely congruent in showing three main clades: one represented the New World vultures; one grouped all falconid representatives, while the other gathered all accipitrid species. Furthermore, the grouping of Neotropical Buteoninae (BNI, BME, RMA, PAL) was systematically retrieved while the Palearctic buteonine species BBU was consistently found closer to an assemblage made of accipitrid (GFU and NNI) as well as pandionid (PHA) species. The main topologic instability originated from the uncertainty in the branching of GBA and HHA (see all 10 most parsimonious trees in Additional file 1: Figure S1). The 50 % majority-rule tree is represented in Fig. 6. A posteriori polarization of character changes along this tree using the outgroup comparison criterion showed that more than 50 chromosomal rearrangements occurred during the karyotypic evolution of the diurnal birds of prey studied here (Fig. 6). Among them, 53 could be unambiguously characterized: 38 corresponded to syntenic associations while 15 were syntenic disruptions (Fig. 6). Three characters (characters 18, 33 and 48) could represent convergent or reversal events. Three other characters were hardly interpretable due to the topologic instability of GBA and/or HHA (characters 22, 41 and 42).
As an alternative approach to the above-mentioned parsimony-based cladistic analysis. these chromosomal characters were also mapped onto a consensual molecular tree generated by combining recent phylogenetic data in the literature [29,32,33] (Additional file 2: Figure S2). This treatment allowed the visualization of potential cytogenetic signature rearrangement(s) that could characterize the major lineages observed on the molecular-based consensus phylogeny. In congruence with the cladistic analysis, such an approach shows that different diurnal birds of prey taxa underwent different chromosomal changes, although many chromosomal characters were inferred to appear independently in different lineages. This is an indication of either true homoplasy or potential uncertainty in the branching as suggested by the molecular topology.

Discussion
The generation of probes from GFU, a species of the Old World vultures, the application of both BOE and GFU probes to comparative painting in more accipitrid representatives as well as a cladistics analysis permitting proper a posteriori polarization of evolutionary events have enabled us to define further the types and number of chromosomal rearrangements that have occurred during the karyotypic diversification of Falconiformes and Accipitriforms, with special emphasis on Falconidae, Accipitridae and Pandionidae.
Our painting results with BOE probes have demonstrated that the two 2n = 66 Old World vultures (GFU Table 1 Data matrix used in the PAUP analysis

Rupornis magnirostris
Note: Chromosomal rearrangement characters used in this table were from previous published data and current study (see the section of materials for details). In the column of "Character (GGA)", the individual numbers or the numbers before alphabets represent the numbers of homologous GGA chromosomes or chromosomal segments. "p": the short arm of a given chromosome; "q": the long arm of a given chromosome; "seg": a segment on a given chromosome; "prx": the part which is near to the centromere of a given chromosome or its arm; "med": the part which is in the medial of a given chromosome or its arm; "dis": the part which is distal to the centromere of a given chromosome or its arm; "mic": the homologous GGA microchromosome. The chromosomal rearrangements in the diurnal birds of prey involved many homologous GGA microchromosomes. However, it is impossible to identify each homologous GGA microchromosomes due to the lack of single GGA microchromosome-specific probes. The numbers after "mic" are arbitrary numbers assigned for GGA microchromosomes involved in different species. GGA 1 is homologous to 4-8 pairs of chromosomes in different Accipitriformes species. Whereby the homologues of GGA 1 segments involved in chromosome arrangements could not be identified unambiguously based chromosome banding and painting data, different symbols were used to represent the homologues of GGA 1 segments in different species, such as 1 seg-n1 and 1seg-n2 in HHA, 1seg-nA, 1seg-nB and 1seg-nC in GBA, and 1seg-na to 1seg-ne in PHA and GHI) have identical karyotypes (Figs. 3a and 5a), and that at least twelve syntenic disruptions and nine segmental associations differentiate the karyotype of BOE from those of GFU and GHI. In a previous study, the hybridization patterns produced by GGA 1-10 probes were shown to be identical also between GFU and GRU (G. rueppellii) [20]. These results indicated that extant Gyps vultures have undergone no interchromosomal rearrangements during their divergence. The hybridization patterns of BOE probes in BBU (Fig. 5b) were very similar to those observed in GFU and GHI. Fourteen syntenic disruptions and ten syntenic associations were revealed in the karyotype of BBU by BOE probes. Among these rearrangements, eleven disruptions (characters 1-4, 18-22, 27 and 28) and eight associations (characters 5, 29-33, 41 and 42) detected in GFU and GHI by BOE probes were also shared with BBU. The karyotype differences between BBU (2n = 68) and the 2n = 66 Old World vultures (GFU and GHI) could have resulted from two further syntenic disruptions of homologous segments of BOE 1 and 10 (= GGA 1 and one GGA microchromosome, characters 23 and 50) and one further syntenic association (BOE 1 + 9 = GGA 1 + 6, character 51). Since multiple chromosomes in BBU, GFU and GHI show homologies with BOE 1 (= GGA 1), it is difficult to define which homologous segment of BOE 1 has undergone further fission in BBU based only on painting results from BOE probes. Further hybridization with GFU probes onto chromosomes of GFU, BOE and BBU allowed us to address this uncertainty (Fig. 7). The results indicate that the GFU 23 (Fig. 7a) homologous segment (= BOE 1pter = GGA 1pter, Fig. 7b) has undergone a syntenic disruption in BBU (Fig. 7c). Therefore, further syntenic disruptions of an ancestral macrochromosome segment (= GGA 1pter, character 23) and a microchromosome (= one segment of BOE 10 = one GGA microchromosome, character 50) could have occurred during chromosome evolution in BBU. Moreover, the complete conservation of  in multiple species studied here and in the rock pigeon (Columba livia, 2n = 80, Columbiformes) as revealed by  suggests that these four pairs of dot-shaped microchromosomes have a much more ancient origin, most likely already being present in the ancestral bird karyotype.  Table 1. Published chromosome painting data for Falconidae (Nishida et al., [12]), Pandionidae (Nishida et al., [18]), Accipitridae (de Oliveira et al., [10,16,19]; Nanda et al., [20]; Nishida et al., [13]), Cathartidae (Raudsepp et al., [34]; Tagliarini et al., [15]) and the data in this study were used for this figure. BBU, Buteo buteo; BME, Buteogallus meridionalis; BNI, Buteo nitida Previous comparative cytogenetic analyses indicated that the variation in karyotypes in the diurnal birds of prey was much wider than that in other avian groups, each family having its distinct karyotype features [3]. To date, including the current study, chromosome painting has been performed in fourteen species belonging to three families of the Accipitriformes and three species of Falconiformes. Only the monotypic family Sagittaridae, represented by the Secretary bird (Sagittarius serpentarius), has not been studied yet by chromosome painting. Table 2 lists the correspondence between GGA 1-5 and their counterparts in BOE and seventeen species belonging to Accipitriformes and Falconiformes. These homology data, together with chromosome painting data on other birds, allowed us to deduce the possible process of karyotype evolution in Accipitriformes through a proper cladistics analysis and mapping the chromosomal changes onto a consensus molecular tree.
In our cladistics analysis, two species (HHA and GBA) showed clear unstable locations in the most parsimonious trees that we retrieved here (see Additional file 1: Figure S1). This was probably due to insufficient resolution of some chromosomal changes and further reciprocal painting involving HHA and GBA will be needed to clarify the relationships of these two species with Accipitridae and Pandionidae. Yet, despite HHA and GBA instability, our chromosome-based cladistic analysis retrieved several cytogenetic signatures that support some groupings within Accipitriformes (Fig. 6), with some of them having already been suggested. For instance, seven characters (3, 18, 20, 21, 27 and 28) united the PHA + all accipitrid species (i.e., Pandionidae and Accipitridae, here) in one monophyletic assemblage; three characters (19, 30 and 32) supported the PHA + GFU (hence GRU, GHI) + NNI + BBU clade; one synapomorphy (character 49) linked Neotropical buteoninae species. The three species of Falconiformes shared eight common characters (6-12 and 54), thus leading to karyotypes that differ greatly from those found in Accipitriformes. Two characters (2 and 4, i.e., the syntenic disruptions of 2p/2q and 5a/5b) seemed to link Accipitriformes and Falconiformes. However, the syntenic disruptions of 2p/2q and 5a/5b (characters 2 and 4) were also found in some Galliformes (different from the chicken) [35] and Strigiformes species [8,23,36]. So further analysis is needed to clarify if these characters will be convergent events or not. In addition, species in Cathartidae (the New World vultures) have typical avian karyotypes, and showed a high degree of conservation in chromosomal synteny with GGA, thus differing from other species in Accipitriformes and species in Falconiformes, including those studied here.
Besides our cladistic analysis, we also mapped the chromosomal rearrangements onto a consensus molecular Fig. 7 Examples of cross-species chromosome painting with GFU probes. a Hybridization of GFU 23 probe to GFU chromosomes. b Hybridization of GFU 23 probe to BOE 1pter. c Hybridization of GFU 23 probe to BBU chromosomes tree (Additional file 2: Figure S2). Although the 50 % majority-rule tree of the diurnal birds of prey that we retrieved using chromosomal rearrangements (Fig. 6) had some differences from the molecular phylogeny (Additional file 2: Figure S2), both analyses (i.e. the cladistic analysis and mapping rearrangements onto an independently obtained topology) revealed a number of potential cytogenetic signatures for some clades investigated here. For instance, both analyses suggested that (1) falconids (Falconiformes) had unique chromosomal rearrangements, differing from that of Accipitriformes species; (2) PHA (the osprey, Pandionidae) shared several common chromosomal rearrangements with other Accipitridae species; (3) the Old World vulture of 2n = 60 (GBA) had some distinct chromosomal characters differing from those 2n = 66 Old World vultures, although the number of common characters was different in both analyses. Some common chromosomal rearrangements, such as characters 19, 30, 32, 33, were mapped onto different clades, i.e. PHA, some of Old World vultures (GFU, GHI, GRU), NNI and BBU. However, in our cladistic analysis, characters 19, 30, 32 supported these species in one assemblage. For the placement of Pandionidae (the osprey), most molecular phylogenetic studies indicate that Sagittariidae, Pandionidae and Accipitridae form a clade, with Pandionidae being more closely related to Accipitridae [26,27,29,32], but one molecular phylogenetic study supported that the osprey lies within the Accipitridae [30].
To date, the reason and mechanism explaining why extensive chromosomal fissions and fusions have occurred only in Accipitriformes and Falconiformes species remains unclear. The possible origin of the peculiar characteristics of the Accipitriformes karyotypes was briefly discussed by De Boer [3]. A translocation model of the evolution of karyotype organization in Accipitridae was also proposed based on computer simulations [37]. However, this model had to be verified by analysis of DNA sequences and localization of transposable elements in Accipitridae species [37]. Further high-resolution mapping and sequence analysis of evolutionary synteny breakpoints are required to elucidate the molecular mechanisms that underlie such extensive chromosomal rearrangements. In addition, the species in Sagittariidae and the two accipitrid species with diploid number differing from that of the majority of accipitrids (2n = 66-68), northern goshawk (Accipiter gentilis) with 2n = 78 and crested eagle (Morphnus guianensis) with 2n = 54, remain to be studied by comparative chromosome painting. The chromosome painting data in these species may be the key to understand the whole karyotype evolution in Accipitriformes.

Conclusion
We here established comparative chromosome maps between BOE and the 2n = 66 Old World vultures (GFU and GHI) and the 2n = 68 BBU by cross-species chromosome painting using two sets of probes from avian species with atypical karyotypes (BOE and GFU). Our results indicated that at least eleven syntenic disruptions and eight segmental associations were shared by the 2n = 66 Old World vultures and the 2n = 68 BBU, thus indicating that BBU could be more closely-related to other accipitrids than to the Neotropical buteoninae species (PAL, BNI, BME and RMI). Our investigations within the diurnal birds of prey have revealed some cytogenetic signatures for different lineages and provided support for the groupings of (i) Falconiformes, and (ii) Accipitridae + Pandionidae. Finally, our chromosome-based cladistic analysis reinforced the proposition that the osprey (PHA) might be a member of Accipitridae.

Methods
Cell culture, metaphase preparations and chromosome nomenclature Three accipitrid species belonging to two different genera, namely the griffon vulture (G. fulvus, GFU), Himalayan griffon (G. himalayensis, GHI) and common buzzard (B. buteo, BBU), and one Charadriiformes species, the stone curlew (B. oedicnemus, BOE) were investigated here. The BOE and GFU fibroblast cell lines were obtained from the Paris Natural History Museum and the GHI fibroblast cell line (KCB200021) was provided by Kunming Cell Bank, the Chinese Academy of Sciences. Mitotic chromosomes of BBU were obtained from a BBU chromosomal suspension stored in Kunming Cell Bank, the Chinese Academy of Sciences. We did not undertake any animal research that would require ethics approval to obtain the samples.
Cell culture, metaphase preparation and slide preparation were carried out following conventional methods [14]. The chromosomes of GFU and BOE were numbered according to previously published karyotypes [3,14]. The karyotype of GHI was arranged according to that of GFU. To facilitate inter-specific comparison of the homologous chromosome segments within Accipitridae, the DAPI-banded BBU karyotype was arranged according to GFU and GHI karyotypes, rather than to the karyotype reported previously [5].

FISH, image capture and processing
Cross-species chromosome painting, image capture and processing were carried out following Nie et al. 2009 [14]. Hybridization signals were assigned to specific chromosomes or chromosome regions defined by inverted DAPI (i.e. 4′, 6-diamidino-2-phenylindole)banding patterns.

Chromosomal character matrix
To date, including the current study, chromosome painting has been performed in seventeen species belonging to four families of the diurnal birds of prey using one of four sets of chromosome painting probes developed from the following species: the chicken (G. gallus, GGA), the white hawk (F. albicollis, PAL), the stone curlew (B. oedicnemus, BOE) and the griffon vulture (G. fulvus, GFU). The establishment of chromosomal homologies between GGA and PAL, GGA and BOE, and between BOE and GFU allowed us to integrate all chromosomal homologies that were already available in the diurnal birds of prey using GGA as a reference. Unlike in other avian species, the syntenies of GGA 1, 2 and 3 were each broken into two pieces in Falconidae, and more than three pieces in Accipitridae and Pandionidae. Based on the reciprocal chromosome painting data between BOE and GGA, PAL and GGA, and BOE and GFU as well as chromosome banding comparison among homologues that correspond to subchromosomal regions on GGA chromosomes 1-3, most of the syntenic breakpoints on GGA 1, 2 and 3 seem to be identical in Accipitridae and Pandionidae. Only a couple of further syntenic disruptions in one or two GGA 1 subchromosomal segments resulted in the difference in the numbers of homologous GGA 1 segments among different species. Due to the lack of reverse painting data for most studied Falconiformes and Accipitriformes species and the lack of single GGA microchromosome probes, we could not determine if all chromosomal rearrangements revealed by GGA probes were of the same origin. For such chromosomal rearrangements that we could not define their origins, we treated them as different characters. Table 1 lists the chromosomal rearrangement characters for seventeen studied Falconiformes and Accipitriformes species.

Cladistic analysis and character mapping
A cladistic analysis using 59 chromosomal characters (Table 1) was conducted as recommended by Dobigny et al. (2004) [40]. Briefly, it consisted in scoring syntenic disruptions or associations as unordered and equally weighted characters, and their presence/absence as character states. Gallus gallus (GGA, Galliformes) was used as the only outgroup. The most parsimonious topologies were investigated through an exhaustive search (initial MaxTrees = 100; MulTree option) conducted under PAUP 4.0 b [41]. The character changes were polarized a posteriori using the outgroup comparison criterion [42].
In parallel to the cladistic analysis, the chromosomal differences inferred by ZooFISH were also mapped onto an independently obtained phylogenetic tree. Since no phylogenetic study including all the species studied in our own present work has ever been conducted, we explored the literature and built a consensus topology on the basis of recently published molecular phylogenetic trees [29,30,33]. We took this latter topology as the reference species tree, and mapped the 59 chromosomal changes onto it.