Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor
- Michael N Romanov†1,
- Marta Farré†2,
- Pamela E Lithgow1,
- Katie E Fowler1, 10,
- Benjamin M Skinner3,
- Rebecca O’Connor1,
- Gothami Fonseka1,
- Niclas Backström4,
- Yoichi Matsuda5,
- Chizuko Nishida6,
- Peter Houde7,
- Erich D Jarvis8,
- Hans Ellegren4,
- David W Burt9,
- Denis M Larkin†2Email author and
- Darren K Griffin†1Email author
© Romanov et al.; licensee BioMed Central. 2014
Received: 11 November 2014
Accepted: 27 November 2014
Published: 11 December 2014
The availability of multiple avian genome sequence assemblies greatly improves our ability to define overall genome organization and reconstruct evolutionary changes. In birds, this has previously been impeded by a near intractable karyotype and relied almost exclusively on comparative molecular cytogenetics of only the largest chromosomes. Here, novel whole genome sequence information from 21 avian genome sequences (most newly assembled) made available on an interactive browser (Evolution Highway) was analyzed.
Focusing on the six best-assembled genomes allowed us to assemble a putative karyotype of the dinosaur ancestor for each chromosome. Reconstructing evolutionary events that led to each species’ genome organization, we determined that the fastest rate of change occurred in the zebra finch and budgerigar, consistent with rapid speciation events in the Passeriformes and Psittaciformes. Intra- and interchromosomal changes were explained most parsimoniously by a series of inversions and translocations respectively, with breakpoint reuse being commonplace. Analyzing chicken and zebra finch, we found little evidence to support the hypothesis of an association of evolutionary breakpoint regions with recombination hotspots but some evidence to support the hypothesis that microchromosomes largely represent conserved blocks of synteny in the majority of the 21 species analyzed. All but one species showed the expected number of microchromosomal rearrangements predicted by the haploid chromosome count. Ostrich, however, appeared to retain an overall karyotype structure of 2n = 80 despite undergoing a large number (26) of hitherto un-described interchromosomal changes.
Results suggest that mechanisms exist to preserve a static overall avian karyotype/genomic structure, including the microchromosomes, with widespread interchromosomal change occurring rarely (e.g., in ostrich and budgerigar lineages). Of the species analyzed, the chicken lineage appeared to have undergone the fewest changes compared to the dinosaur ancestor.
KeywordsAncestral karyotype Avian genome Chromosome evolution Dinosaur
The mechanisms of genome evolution are most often considered from the perspective of individual genes or gene families; there is nonetheless increasing evidence supporting the functional role and significance of events at a chromosomal (cytogenetic) level . To date, bird genomes remain relatively understudied from an overall genome organization perspective; however, the recent availability of multiple avian genome sequence assemblies  allows us to consider the role of chromosomal change in the evolution of Aves from their dinosaur ancestors. Chromosome rearrangements between species can cause or reinforce reproductive isolation through reduced fitness of hybrid offspring due to a compromised ability to synapse and segregate chromosomes at meiosis [3, 4]. Moreover, reduced interspecific recombination in rearranged regions is thought to promote the accumulation of incompatibility loci in such regions [5–7]. The purpose of this study was to gain further insight into the mechanism of bird evolution through the multiple comparative analyses of chromosomal segments and breakpoints.
Unraveling the mechanisms and relevance of bird karyotype evolution has hitherto been impeded by a karyotype that is difficult to define because of indistinct banding on the macrochromosomes and a preponderance of cytogenetically indistinguishable microchromosomes. Indeed, to date, only a single avian karyotype (chicken) has been fully defined using a combination of BAC/cosmid clones and chromosome paints generated by flow cytometry and microdissection . Moreover, karyotypes are broadly similar in overall pattern from species to species. For instance, at a cytogenetic level, two thirds of bird species have a chromosome number of around 2n = 80 with similar numbers of macro- and microchromosomes suggesting little interchromosomal changes between species . Molecular insights into interchromosomal differences between species (and the evolutionary events that have led to them) have focused mostly on the largest macrochromosomes. These studies applied chicken chromosome paints  to the chromosomes of numerous other species (reviewed in ) in zoo-FISH experiments. Such investigations have provided much insight into inter-macrochromosomal rearrangements between birds with the underlying message that the ancestral pattern has remained largely unaltered in the majority of species. Rare exceptions include significant chromosome rearrangement in Psittaciformes (parrots etc.), Falconiformes (falcons) and Sphenisciformes (penguins) . There are also individual changes associated with representative orders, e.g., fission of chromosome 1 in Passeriformes (songbirds) and of chromosome 2 in certain Galliformes (land fowl) (reviewed in ). Studies of interchromosomal changes involving the microchromosomes are much more limited as the flow cytometry methods used to generate the chromosome paints  do not have the resolution to isolate individual microchromosomes.
Using chicken BAC clones, studies provide a low-resolution appraisal of intrachromosomal rearrangements between chicken and other species [12–14] (turkey, duck, zebra finch, respectively). This approach, however, is limited in its ability to identify the molecular coordinates of evolutionary breakpoints. The availability of whole assembled genomes [15–17] allows comparative genomics at a much more detailed level of resolution than can be achieved by cross-species FISH. Burt et al. were the first to use bioinformatics to define inter-species analysis of whole avian chromosomes at a genomic level (chicken-human). The publication of the chicken genome sequence  provided more detailed information, establishing conserved synteny between chicken and human whole genome assemblies. In the ten years since, only conserved synteny comparisons have been made between the chromosomes of two [14, 19], or at most three [20, 21] avian species.
The use of whole genome assemblies to study cytogenetic phenomena has raised interest in the study of comparative cytogenetics from the perspective of evolutionary breakpoint regions (EBRs) and homologous synteny blocks (HSBs). To date, the majority of such studies have focused on mammals , however, analysis of other groups, such as birds, is essential in order to establish whether mammalian systems are representative of, or an exception to, general patterns observed in other animal groups. Larkin et al. found that, in mammals, EBRs can lie in gene-dense regions. In the human genome EBRs also lie in regions with more zinc finger protein genes, more genes whose function is associated with environmental stimulus response, as well as more segmental duplications, CNVs, SNPs and retrotransposed genes. Such “EBR genes” appear to be related to lineage-specific biology and adaptive features [22–24]. EBRs are also frequently reused, i.e. there are regions of the genome that are prone to chromosomal breakage leading to translocations, inversions and fissions [25, 26]. Comparison of sequence assemblies in chicken, zebra finch and turkey suggests that breakpoint reuse is higher in birds than in mammals [20, 21]. The data in birds also suggests a key role for recombination-based mechanisms in the generation of chromosome rearrangements in that EBR location is consistent with elevated levels of genetic recombination at these loci . This is consistent with the notion that, if recombination drives chromosomal rearrangements and assuming an evolutionarily conserved recombination landscape [27–29], EBRs might be enriched in genomic regions with elevated recombination rates. Not all species show an association of chromosomal breakage and elevated recombination however, e.g., insects [30, 31] and mammals. Indeed, in mammals Larkin et al. suggested that the highest levels of recombination are located between the EBRs rather than in association with them.
HSBs have been defined in all animal species thus far examined for conserved chromosomal synteny . Larkin et al. argue that the continued presence of HSBs in all species may indicate a selective advantage to the retention of gene combinations in close proximity. Supporting evidence is found in the fact that multispecies HSBs (msHSBs) involving nine mammals plus chicken, unlike EBRs, are enriched in gene ontology (GO) terms for organismal development, central nervous system, and brain function in the human genome. Others argue that the idea of close proximity and any resulting correlation in expression patterns (if present) are not necessarily adaptive or required (e.g., [33, 34]). Given that around three quarters of avian chromosomes are small, cytogenetically indistinguishable microchromosomes, and that overall karyotype structure appears broadly similar between at least two thirds of bird species, a high degree of conserved chromosomal synteny is inferred . This raises the hypothesis that avian karyotypes are evolutionarily static; however, for this to be tested, we would first need to establish that inter-microchromosomal rearrangements are rare or absent in most birds. If true, we would subsequently hypothesize that, like HSBs in mammals, individual whole microchromosomes are enriched for functional GO terms (regardless of any intrachromosomal rearrangements between them).
A detailed account of the chromosomal differences and changes that have occurred during the evolution of avian chromosomes is an essential prerequisite for any further insights into functional and/or mechanistic relevance. The combination of comparative analysis by bioinformatics and chromosome painting has the potential to do this, provided the appropriate tools are developed and used. The purpose of this study was thus to examine multiple avian genomes recently sequenced [2, 35], reconstruct the common ancestral karyotype and thence the evolutionary events that led to extant karyotypes. Furthermore, we tested the hypothesis that EBRs occurring in two lineages (chicken and zebra finch) are associated with elevated levels of genetic recombination and assessed the degree to which EBRs are reused in avian evolution. Finally, we tested the hypothesis that whole microchromosomes essentially constitute interchromosomal HSBs (i.e. that rearrangements between them are rare or absent) and that each microchromosome consists of functionally enriched GO terms.
Genomic data and visualization of HSBs and EBRs
Reconstructions of scaffold-based assemblies also relied, in part, on previously published zoo-FISH (BAC and chromosome painting) data for the macro- and microchromosomes of chicken, turkey, duck and zebra finch [12–14] as well as newly generated data in this study as follows: we used seven new chicken microchromosomal paints A–G , verifying their assignments with chicken BACs (see Additional file 1) by dual color FISH and painting them onto ostrich and budgerigar metaphases.
Comparative mapping of chicken chromosome paints A–G, and their ostrich and budgerigar orthologs
Chromosome paint ID
Ostrich orthologs (all microchromosomes)
Fusion as part of chromosome 5
10 and 12
2 pairs of microchromosomes (no apparent fissions/fusions at this resolution)
1 pair of microchromosomes (no apparent fissions/fusions at this resolution)
13 and 14
1 microchromosome pair +1 arm of chromosome 8 = fission and fusion at this resolution
10 and 12
1 pair = fusion
16, 17 and 18
2 pairs = fusion
~5 pairs smaller than 18
3 pairs = 2 fusions (although some signals are weak so may be failure of hybridization)
Reconstruction of ancestral karyotypes and chromosomal changes
Total numbers of inter- and intrachromosomal rearrangements since divergence from avian ancestor 100 MYA
No. of interchromosomal changes (as determined by FISH) from avian ancestor
No. of interchromosomal changes (determined using bioinformatics) from avian ancestor
No. of intrachromosomal changes from avian ancestor in chromosomes 1–5 (excluding 4p)
No. of intrachromosomal changes from Neognathae ancestor in chromosomes 6–28 + 4p + Z
Closer analysis of the breakpoints to address the question of breakpoint reuse (see Background) identified, in chicken chromosomes 1–5 (and their turkey, duck, zebra finch, budgerigar and ostrich orthologs), 620 segment ends, of which 421 were involved in rearrangements. The most parsimonious predicted pathways from the common avian ancestor suggested that 100 breakpoint regions (23.8%) recurred in different lineages, whereas 214 breakpoint regions (50.8%) recurred in either the same or different lineages. In chicken chromosomes 4p, 6–28 and Z, and their turkey, duck, zebra finch and budgerigar orthologs, 560 segment ends were identified, of which 428 were involved in rearrangements. The most parsimonious predicted pathways from the common avian ancestor suggested that 109 breakpoint regions (25.5%) recurred in different lineages, whereas 210 breakpoint regions (49.1%) recurred in either the same or different lineages.
EBRs and recombination in chicken and zebra finch
As also mentioned in the Background section, we tested the hypothesis that the presence of EBRs was related to the regional recombination rate. Given the quality of the genetic maps and the data available in this study, this could be achieved for the chicken and zebra finch only.
In zebra finch, 1 Mb non-overlapping windows with EBRs (n = 31) had a slightly higher recombination rate than windows without (n = 952; 1.60 vs. 1.29 cM/Mb), although this was not statistically significant (Wilcoxon’s test, P = 0.1; randomization test, empirical difference in mean between classes = 0.31, P = 0.1; Figure 5).
Interchromosomal changes in multiple species and GO of microchromosomes
Total number of interchromosomal rearrangements involving microchromosomes in 21 avian species compared to chicken
Total number of interchromosomal changes involving macro- and microchromosomes
Interchromosomal changes between microchromosomes only
Haploid chromosome number (difference fromn = 40)
Adélie penguin (Pygoscelis adeliae)
American crow (Corvus brachyrhynchos)
Common cuckoo (Cuculus canorus)
Pekin duck (Anas platyrhynchos)
Little egret (Egretta garzetta)
Emperor penguin (Aptenodytes forsteri)
Peregrine falcon (Falco peregrinus)
Zebra finch (Taeniopygia guttata)
Hoatzin (Ophisthocomus hoazin)
Anna’s hummingbird (Calypte anna)
Crested ibis (Nipponia nippon)
Killdeer (Charadrius vociferous)
Golden collared manakin (Manacus vitellinus)
Medium ground finch (Geospiza fortis)
Ostrich (Struthio camelus)
Budgerigar (Melopsittacus undulates)
Rock dove (Columba livia)
Chimney swift (Chaetura pelagica)
Turkey (Meleagris gallopavo)
Downy woodpecker (Picoides pubescens)
Chicken (Gallus gallus)
Once we had established (above) that rearrangements were rare in the microchromosomes, then this led to the hypothesis that each microchromosome contained functionally enriched GO categories (see Background). We found evidence to support this hypothesis only for chromosome 16 (enriched for immune function) when P < 0.05 and a false discovery rate (FDR) threshold of 0.05 were applied. Nonetheless several chromosomes had a significant P value but did not pass the FDR threshold: for chromosome 11 enrichment categories were apparent for drug/caffeine metabolism as well as hemophilic cell adhesion; for chromosome 12 genes for nucleotide binding were clustered together; for chromosome 13 there were enrichment categories for GTPase regulator activity; phosphatase activity in chromosome 15; chromosome 17 for glycosylation and glycoprotein related processes; chromosome 18 for cytoskeletal and motor protein related genes; and chromosome 20 for genes involved in apoptosis and cell death.
We thus find evidence to support our hypothesis that microchromosomes represent highly conserved blocks of interchromosomal synteny but find limited evidence to support the hypothesis that one possible explanation for this is a clustering of genes of associated function on the same chromosome.
The results presented here signify the most comprehensive appraisal of avian comparative cytogenetics to date. They provide a more detailed reconstruction of avian genome evolution than could be achieved by zoo-FISH analysis alone and demonstrate proof of principle from which further studies of genome evolution and comparative genomics can ensue.
We used a highly interactive avian genome dataset from the Evolution Highway comparative chromosome browser [37, 38] that, as has already been demonstrated in mammals, can be applied to compare the chromosome organization of individual or multiple species. The ultimate aim for this browser is that, in chromosomes for all avian species uploaded, HSBs will be displayed with reference to the chromosome number, as is currently the case for turkey, zebra finch and duck, or to specific scaffolds for other birds. In future, this will be achieved by a number of strategies: (a) by improved scaffold sizes, e.g., using optical mapping such as has been achieved to some degree in ostrich and budgerigar in this study; (b) by linkage to radiation hybrid (RH) maps such as was achieved for duck in this study (see also ); (c) by association with known linkage and other physical maps (e.g., [39, 40]); d) by use of novel algorithms to order and orient scaffolds into longer chromosomal fragments or whole chromosomes using comparative genome information and pair-end reads (reference-assisted chromosome assembly; ); (e) by systematic FISH mapping to chromosomes of orthologous clones derived from the individual scaffolds. We are currently concentrating our efforts on the development of FISH probes that will identify not only on which chromosomes the scaffolds lie in the species of interest, but also the order in which they appear on the chromosome. With current technology, however, even the best-assembled genomes (e.g., assisted with optical mapping) require a degree of intervention by molecular cytogenetics in order to generate a complete picture of overall genome organization. Given the efforts that have been made to sequence the genomes of the birds recently by current technologies , it is questionable how many of them will be re-sequenced using newer technologies that generate large scaffolds. A note of caution is relevant here: no genome assembly is “perfect” - the results reported here and elsewhere represent the state of the art in terms of what can be reasonably gleaned with the current technology available. Our future studies will focus on the systematic molecular characterization by zoo-FISH of as many scaffolds and EBRs as time and resources allow.
Earlier cytogenetic data suggested that, for the majority of bird species, karyotypic patterns are broadly similar to one another [9, 11, 14, 20]. This purportedly extends to ratite birds [42–44]; however, further analysis presented in this study challenges this notion. That is, we identified 26 interchromosomal rearrangements in ostrich compared to the ancestor. Moreover, the question of whether the conserved interchromosomal synteny seen in the macrochromosomes applies to the microchromosomes has hitherto been beyond the resolution of contemporary methodology. This study is the first to classify inter-microchromosomal rearrangements in any species; we provide evidence that interchromosomal rearrangements are nonetheless rare, except in cases (around 1/3 of species) where we already knew that karyotypes were highly rearranged . Ostrich is the exception and it will be interesting to note whether this applies to other ratite birds.
Microchromosomes are not a uniquely avian feature. They are also found in some primitive amphibians (Cryptobranchidae and Hynobiidae have 14–19 pairs [45, 46]), most (but not all) reptiles (snakes have around 20 pairs ), but paradoxically not Crocodylia  – the closest phylogenetic lineage to birds. Indeed microchromosomes are typical of most amniotes (mammals and crocodilians being exceptions); however, the greatest number and smallest size of microchromosomes are typically found among birds. Burt  in a “fission-fusion” hypothesis suggested that most microchromosomes were already present in the common dinosaur ancestor that gave rise to birds (which probably had already evolved a small genome size and karyotype of around 2n = 60 including 20 pairs of microchromosomes) but that chromosome fission created the remainder, presumably including the smallest ones. In the current study, the similar number of chromosomes amongst most species but relatively large number of rearrangements between ostrich and all the other birds studied suggest that a basic pattern of 2n = 80 (~30 pairs of microchromosomes) became fixed before the Palaeognathae-Neognathae divergence 100 MYA but that interchromosomal rearrangement was still relatively common in birds at the time. Another alternative is that ratite birds underwent further adaptive changes that may be associated with the very different phenotypes present in this clade alone. The paucity of inter-microchromosomal rearrangements between most Neognathae (if the evidence presented here is representative, this would presumably include the 2/3 of Neognathae species where 2n = ~80) supports our hypothesis that the microchromosomes represent blocks of conserved synteny at an interchromosomal level. An absence of interchromosomal rearrangement could either suggest an evolutionary advantage to retaining this particular configuration or a lack of opportunity for chromosome rearrangement. The latter might be explained by few recombination hotspots, transposable elements or endogenous retroviruses, all of which have been associated with chromosomal change. Both inter- and intrachromosomal change can arise via these mechanisms, and thus the rapid amount of intrachromosomal but not interchromosomal change in our representative passeriform species, the zebra finch, suggest that there may be an evolutionary advantage to keeping microchromosomes numerous, gene dense, compact and evolutionarily static. Stasis in evolution can, however, arise via alternative interpretations; it may be that the mutational mechanisms underlying chromosomal changes are different in birds or that lack of adaptive value, rather than purifying selection, slows down the rate of chromosomal changes. At the time of writing no sequences have yet been associated with the very smallest of the avian microchromosomes (29–38) and this is an issue that will require rectifying in future avian genome projects using more sophisticated technologies.
The rate of chromosomal change in any eukaryotic organism, and the speciation that ultimately arises from it, is dependent on two factors: the rate of mutation and the rate of fixation . The mutation rate of chromosomes is, in turn, related to the frequency of homologous sites . Repeat structures in general, and transposable elements in particular, provide substrates for chromosomal rearrangement. In a genome that is constrained by size (perhaps, as has been suggested, because of the energy requirements associated with flight [50, 51]), the opportunity for mutation is reduced and only fission (or intrachromosomal rearrangement such as inversion) can occur. This would explain first why the avian genome is the most fragmented of any vertebrate genome (i.e. birds have the most chromosomes) and second why there have been few interchromosomal rearrangements in most species. There are also possible advantages of multiple chromosomes in a karyotype in terms of generating variation, the driver of natural selection. That is, more chromosomes lead to more combinations of gametes as well as an increase in recombination rate as there has to be at least one obligatory chiasma per chromosome. The absence of positive selection for much change in chromosome number is a possible explanation of why there was little fixation of any interchromosomal changes among birds although inbreeding and genetic drift may play a role [18, 49, 52, 53]. Burt  suggested that a higher recombination rate is another constraint that has resulted in the properties we most associate with microchromosomes (e.g., high GC-content, low repeats, high gene-density) and led to the maintenance of the typical avian karyotype with both macro- and microchromosomes and few rearrangements between them.
A constraint of overall karyotype structure does not preclude intrachromosomal rearrangements. Indeed there is a correlation between the rates of speciation and intrachromosomal rearrangement . In the current study, the rapid rate of intrachromosomal rearrangement in the zebra finch would argue for a relationship between intrachromosomal rearrangement and speciation in birds given the Passeriformes represent over half of all species. Such mechanisms could be mediated through an increase in localized repeat content. Hotspots of recombination have previously been reported to also play a role  and in this study we tested the hypothesis further utilizing “zebra finch only” and “chicken only” breakpoints comparing them to previously reported genetic maps of each species [37, 54, 55]. In chicken, recombination rates were near identical in regions with breakpoints compared to those without. In zebra finch, the difference in rates between regions containing EBRs and regions without EBRs, although similar in magnitude to that previously reported , failed to reach statistical significance (at P < 0.05). This therefore casts doubt on our original findings, thereby either suggesting that our hypothesis should be rejected or that the numbers in the study were not sufficiently large to reach statistical significance. A further alternative explanation is that the available recombination maps have too low marker density (typically Mb scale) to pick up local recombination rate variation at a sufficiently detailed scale (Kb scale) to detect associations with EBRs. Study of a greater number of species in this manner using high-density linkage maps or population based recombination rate estimates may resolve the paradoxical difference between  and the current study.
Some avian species undergo a radical departure from the typical (2n = ~80) avian genome organization. The presence of an unusually high chromosome number in the Adélie penguin (2n = 96) and a lower than average number in the emperor penguin (2n = 72) (but both associated with high degrees of inter-microchromosomal rearrangement) suggest that similar mechanisms can act to either reduce or increase chromosome number rapidly. Evidence from the penguins and the rearranged karyotypes of the Falconiformes and the Psittaciformes suggest that these changes can happen in a relatively short time. Mammals, reptiles and amphibians with larger, repeat-rich genomes have the potential to undergo rapid intra- and interchromosomal rearrangements and the results presented here suggest that birds too can undergo similar changes in certain groups. We are not, however, aware of any evidence to suggest that highly rearranged avian genomes are especially large, or significantly more repeat-rich than other avian genomes. Comparisons of the zebra finch and the budgerigar suggest that mutation rates of chromosomes may well be similarly high in both groups but that they are features associated with exploiting evolutionary niches in certain groups that serve to fix interchromosomal rearrangements, while in others such fixation is prevented and the overall avian karyotype maintained. Such processes are, to date, undiscovered but possible clues might lie in the study of GO terms present in EBRs. In an associated study, a correlation between EBRs and specific avian adaptive features in individual species has been demonstrated. This included forebrain development in budgerigar, one of the six species focused upon in this study and consistent with this species being not only vocal-learner but having distinctive neuronal connections compared to other vocal-learners . As more genomes become available with better assemblies, these analyses may well point to adaptive phenotypic features of individual orders and families.
Finally, we observed that it appears to be the chicken that seems to have undergone the fewest chromosomal changes compared to the ancestor. There are interesting parallels between this study and another study  examining sex chromosome evolution. While our data demonstrates that autosomes have been reorganized least in chicken chromosomes 1–5 in comparison to the common avian ancestor, Zhou et al. conclude that the ancestral sex chromosome organization is observed closer to that of the Palaeognathae (ostrich and emu). Zhou et al. show less degradation of the sex chromosomes and a closer synteny to the lizard. As, in this study, we only examined the Z chromosome in the Neognathae (for the reasons given), further studies will be required to establish whether sex chromosomes and autosomes preserve their ancestry differently in the different lineages. The question also arises of whether chicken and related species, having undergone the fewest chromosomal changes, have undergone the fewest adaptive changes compared to the avian ancestor. Most authors agree that the dinosaur ancestors of birds were bipedal and terrestrial, relatively small (small size being an immediate pre-adaptation to flight) and had limited flying ability, not unlike Galliformes . On the other hand, the earliest known Ornithurae along the presumed direct line to modern birds were either fully aquatic or amphibious (e.g., Gansus) and details of their anatomy, including webbed feet, have been likened to ducks [59, 60]. The oldest relatively certain fossil representative of Neornithes (modern birds) is aquatic, and identified as a Galloanseres (e.g., Vegavis). However, the fossil record may be difficult to interpret due to geographic and depositional sampling biases, limited understanding of functional anatomy, and the uncertainty that avian ancestors were ecologically and behaviorally typical of the larger groups to which they belonged. As an independent record of the actual substance of inheritance of living birds, genomic characteristics such as chromosomal arrangement complement a fossil record that may imperfectly represent actual neornithine forebears. Thus, chromosomal rearrangements may provide information on the ecological adaptations of avian ancestors that the fossil record may never be able to establish unambiguously .
In summary, this study represents the most comprehensive appraisal of changes in overall avian genome structure hitherto reported. We provide further insight on previously reported roles of genetic recombination in chromosome rearrangement and on the functional significance of karyotype stability in the avian genome. Here, we establish that the chicken lineage contains the fewest number of chromosomal changes compared to the dinosaur ancestor relative to the other five species studied. At this stage it would be unwise automatically to infer that this means that the chicken has the fewest number of adaptive changes also. This will nonetheless be the topic of future study.
Presentation of multiple avian genome assemblies
In order to present and visualize comparative cytogenetics and identify HSBs and EBRs in multiple avian species, an interactive, comparative chromosome browser Evolution Highway was used . All blocks of synteny were identified and displayed relative to chromosomes of the reference chicken genome (ICGSC Gallus_gallus-4.0/galGal4). Evolution Highway was used to display the sequence coordinates of all syntenic fragments (SF) and HSBs in each genome ). We made use of the set of HSBs and SFs that contained rearrangements that are ≥ 300 Kb in the reference genome. This set, together with two other separate sets that visualize HSBs and SFs that are larger than 100 Kb and 500 Kb in the reference genome, is publicly available from the Evolution Highway website  (Figure 1) and are further described in .
For the purposes of this study, 21 avian genomes plus one outgroup species were utilized to address the questions set out in the Background section and made up of the following: of these 21, 17 were recently sequenced and presented  including common cuckoo, peregrine falcon, American crow, little egret, crested ibis, domestic pigeon, hoatzin, golden-collared manakin, medium ground finch, downy woodpecker, Adélie penguin, emperor penguin, Anna’s hummingbird, chimney swift, killdeer, budgerigar and ostrich. Conserved blocks of synteny are presented as scaffolds (scaffold 1 being the largest and the rest numbered accordingly to size) in relation to chicken chromosomes. Chromosome-level assembly and analysis of conserved synteny had been previously reported for the largest (macro-) chromosomes of chicken, turkey and zebra finch [14, 20, 21]. Thus, the turkey (TGC Turkey_2.01/melGal1) and zebra finch (WUGSC 3.2.4/taeGut1) genomes were presented in Evolution Highway with reference to published chromosome number (e.g., chromosome 11 in chicken corresponds to chromosome 12 in duck and 13 in turkey; see Figure 1). Chromosome-level assembly of the Pekin duck genome was constructed from available genome scaffolds  using an original RH mapping approach through hybrid sequencing (Faraut et al., personal communication). Pekin duck was added and presented with reference to published chromosome number. The Carolina anole was the only reptile outgroup genome available with reference to whole chromosomes and therefore this was chosen for this study as the outgroup for reconstruction of the ancestral chromosomes (see the sub-section Establishment of ancestral avian karyotypes).
Of the 17 newly sequenced species, two (ostrich and budgerigar) were selected for studies involving reconstruction of the ancestral chromosomes. These species, thanks to optical mapping, had the largest N50 (>10 Mb) and were also the species on which we performed zoo-FISH studies due to the availability of material for chromosome preparation. These and the remaining 15 species were used for defining EBRs to compare with recombination rate and for establishing interchromosomal conserved synteny among the microchromosomes .
Karyotype and zoo-FISH analysis
For chromosome analysis, rapidly dividing embryonic fibroblasts or white blood cells were arrested in metaphase using colchicine (Sigma), swollen using 75 mM KCl and fixed to glass slides using 3:1 methanol : acetic acid mix. Metaphases were stained with a combination of DAPI and propidium iodide in VECTASHIELD® antifade medium (Vector Laboratories). Image capture involved an Olympus BX61 epifluorescence microscope with cooled CCD camera; SmartCapture system and SmartType software (Digital Scientific UK) were used for capturing and karyotyping purposes, respectively. Microchromosome paints described elsewhere  were generated by flow cytometry, then amplified and directly labeled with FITC using DOP-PCR. BAC clone DNAs were used to verify chromosome paint alignment and were extracted by miniprep (QIAprep Spin Miniprep Kit, QIAGEN), then directly labeled by nick translation with FITC or Cy3.5.
For FISH, metaphases were probed with chicken chromosome paints and BACs generated above. Briefly, probes were dissolved in a formamide buffer and applied, under a coverslip, and then sealed using rubber cement. Simultaneous denaturation of probe and genomic DNA on a 75°C hotplate preceded hybridization at 37°C (overnight for same species FISH, three days for zoo-FISH). Post-hybridization washes (2 minutes in 0.4 × SSC at 73°C; 30 seconds in 2 × SSC/0.5% Tween 20 at room temperature) were followed by chromosome counterstaining using VECTASHIELD® anti-fade medium with DAPI and viewed as above using epifluorescence and SmartCapture (Digital Scientific UK).
Establishment of ancestral avian karyotypes
In total six avian species (chicken, turkey, duck, zebra finch, ostrich and budgerigar) plus one lizard outgroup species (Carolina anole) were chosen for reconstruction of the ancestral karyotypes (for the reasons given in the sub-section Presentation of multiple avian genome assemblies). A combination of bioinformatics, zoo-FISH and karyotyping allowed us to make reconstructions of the order and orientation of scaffolds and thence the ancestral chromosomes. To reconstruct a putative avian ancestor as inferred from orthology maps the Multiple Genomes Rearrangements and Ancestors (MGRA) tool on the Algorithmic Biology Lab web server at St. Petersburg Academic University of the Russian Academy of Sciences [64, 65] was used as follows: using Evolution Highway, pairwise alignments for turkey, duck, zebra finch, budgerigar and ostrich were visualized relative to the chicken whole genome sequence as a reference at the 300 Kb resolution. The orthology map of the Carolina anole, also visualized by Evolution Highway, was used as an input for the MGRA program and included in the analysis as an outgroup. Orthologous regions observed in all the species compared were defined as msHSBs and served as MGRA inputs for individual genomes. The hypothetical ancestral genome was determined using the phylogenetic tree information for this set of six species .For chromosomes 1–5, 80% of the avian genomes were also represented by orthologous sequences in the Carolina anole outgroup. In this case we could therefore reconstruct the ancestral chromosomes for all birds. For chromosomes 6–28 and Z, we used ostrich as the outgroup (thus only drawing conclusions about the Neognathae), as only ~9% of the genome had orthologous sequences represented in the lizard outgroup. Where the ostrich and Neognathae ancestor had the same arrangement of HSBs, we could infer the avian ancestor (as with chromosome 11, Figure 3).
In order to test the robustness of our analysis in a series of additional MGRA simulations and iterations, we established if exclusion of one species at a time from the set of six species would affect the overall pattern of the reconstructed avian ancestor genome organization.
Reconstruction of evolutionary events guided by MGRA
The positions of CARs and HSBs or SFs within each species genome were noted, allowing correlation with our previously published FISH based physical mapping data in chicken turkey, duck and zebra finch [12–14] and that derived by cross-species chromosome painting in former publications [66, 67] and in the current study. These data were previously acquired by cross-species FISH of chicken BACs and chromosome paints onto turkey, duck, ostrich and budgerigar chromosomes, and same-species FISH of orthologous zebra finch BACs onto zebra finch chromosomes.
The available karyotypic, FISH and bioinformatic data were combined to generate the “best-fit” model for chromosomal evolution in the six avian species of interest, i.e. the one with the minimum number of rearrangements. The MGRA tool was used on the whole genome datasets to reconstruct the evolutionary events that, most parsimoniously, led to the arrangement seen in the extant species. For the most part, the changes suggested by MGRA were accepted as the most parsimonious involving the minimum inversions for intrachromosomal rearrangements and fissions/fusions for interchromosomal rearrangements (the process of defining the inversions is illustrated in Figure 3; see also ). In cases where apparent interchromosomal rearrangements (such as translocations) had occurred, the MGRA solution was cross-referenced with the reconstructions on a chromosome-by-chromosome basis using the Multiple Genome Rearrangements (MGR) tool [68, 69] and with zoo-FISH data. In cases of disagreement on the pattern of rearrangements, three independent observers with extensive cytogenetic expertise manually checked and decided the rearrangement pattern. When a whole, otherwise independent, block (scaffold or chromosome) was classed as inverted, this was counted in the analysis as a true inversion if a different orientation was recovered for two or more species (example shown in Figure 3b for chromosome 11 in zebra finch).
Identification of EBRs and breakpoint reuse
We used the EBRs defined in  that involved a single reference chromosome (intrachromosomal EBRs) and more than one reference chromosome (interchromosomal EBRs) in target species’ chromosomes or scaffolds . Interchromosomal EBRs delineated interchromosomal rearrangements, which were then compared with published chromosome number , or more specifically deviation from n = 40; correlation coefficient R 2 was calculated using Microsoft Excel. In order to determine breakpoint reuse, the series of possible rearrangements from the common avian ancestor (with lizard as the outgroup, chromosomes 1–5) or Neognathae ancestor (with ostrich as the outgroup, chromosomes 4p, Z and 6–28) to each species was considered, and for each rearrangement, the segment ends flanking the breakpoints were noted. Within each lineage, the number of times a segment end was involved in a rearrangement was counted and reuse classified if it occurred more than once in any lineage or between lineages.
Recombination rate analyses
We used the chicken- and finch-specific EBRs defined in  to compare with chicken-specific recombination rates and zebra finch-specific EBRs with zebra-finch recombination rates. This differed from our previous approach  in which we examined all EBRs between three species compared to the zebra finch genetic map. Zebra finch-specific EBRs coordinates initially identified in chicken chromosomes were translated into zebra finch chromosome coordinates (WUGSC 3.2.4/taeGut1) using the correspondence between coordinates of finch HSB boundaries in the chicken and finch chromosome assemblies . In this way all chicken-specific and zebra finch-specific EBRs identified at 300 Kb resolution were compared directly with genetic maps in chicken and zebra finch, respectively.
We obtained sex-averaged recombination rate estimates for 1 Mb non-overlapping windows by comparing genetic and physical positions of SNPs distributed along the chicken and zebra finch genomes (data from [54, 55]). To assess if the recombination rate differed between regions with and without chromosomal breakpoints, we partitioned the recombination data into two classes, one with windows containing at least one breakpoint and one with windows without breakpoints, using the zebra finch and chicken breakpoint data . We applied a non-parametric test (Wilcoxon’s rank sum test with continuity correction as implemented in R ) to assess the level of significance for the difference in recombination rates between classes. Since the sample size differed considerably between classes (i.e. windows not containing EBRs vastly exceeded those that contained EBRs) we also applied a randomization test in R . We randomly sampled the same number of windows as those containing EBRs in each respective taxon (n = 31 for zebra finch, n = 35 for chicken) from the entire sample 10,000 times. Lastly, we calculated the average recombination rate in the random sample of windows for each iteration to obtain an expected distribution.
GO analysis of microchromosomes
In order to ask whether individual microchromosomes were enriched for specific GO categories, whole gene sets for each microchromosome were collated and loaded both into DAVID [72, 73] and GOEAST [74, 75]. Specifically, Ensembl gene ID data and gene name for each microchromosome were extracted from the BioMart Ensembl Genes 75 Database [76, 77], using galGal4 as the dataset. In order to eliminate any “significant” results arising through the presence of multiple copies of genes in the same family being present on the same chromosome, gene families were reduced to a single representative member. Downloaded gene IDs and gene names were then copied into a spreadsheet for further analysis using DAVID and GOEAST. Gene IDs for each microchromosome were uploaded into DAVID Bioinformatics Resources 6.7, using Ensembl Gene ID as the list identifier and subsequently analyzed using the Functional Annotation Clustering tool. Cluster data from each microchromosome gene list output was downloaded into Microsoft Excel and filtered using an enrichment score of 1.3 and above and a P value less than 0.05 to edit the list for clusters considered to be significant. BioMart (Ensembl) derived gene names for each microchromosome were also uploaded into GOEAST using Gallus gallus as the reference. Batch-gene analysis was performed by GOEAST, and enriched GO term outputs with a P value less than 0.05 were considered to be significant. The GO results obtained from GOEAST were downloaded into Microsoft Excel and presented with graphic files created directly from GOEAST for each microchromosome where results were available. Finally, in order to correct for multiple sampling error, an FDR threshold of 0.05 was used.
Michael N Romanov and Marta Farré, joint first authors.
Denis M Larkin and Darren K Griffin, joint last and corresponding authors.
Bacterial artificial chromosome
Contiguous ancestral region
Copy number variation
Evolutionary breakpoint region
False discovery rate
Fluorescent in situ hybridization
Homologous synteny block
Multispecies homologous synteny block
Multiple Genome Rearrangements
Multiple Genomes Rearrangements and Ancestors
Million years ago
Single nucleotide polymorphism.
The authors would like to thank Alain Vignal and Thomas Faraut of INRA Toulouse (France) for access to the duck chromosome assembly data. This research was funded in part by PL-Grid Infrastructure (DML), Biotechnology and Biological Sciences Research Council BB/K008161 (DML, DKG), BB/K008226/1 (DML), BB/J010170/1 (DML, MF) and a knowledge transfer partnership award (DKG and Cytocell Ltd). The authors are grateful to Malcolm Ferguson-Smith’s lab (Cambridge, UK) for producing the flow-sorted chicken microchromosome paints. We also thank Cytocell Ltd (Cambridge, UK) for technical support in FISH technologies.
- Lewin HA, Larkin DM, Pontius J, O'Brien SJ: Every genome sequence needs a good map. Genome Res. 2009, 19: 1925-1928. 10.1101/gr.094557.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhang G, Li C, Li Q, Li B, Larkin DM, Lee C, Storz JF, Antunes A, Greenwold MJ, Meredith RW, Ödeen A, Cui J, Zhou Q, Xu L, Pan H, Wang Z, Jin L, Zhang P, Hu H, Yang W, Hu J, Xiao J, Yang Z, Liu Y, Xie Q, Yu H, Lian J, Wen P, Zhang F, Li H, et al: Comparative genomics reveals insights into avian genome evolution and adaptation. Science. 2014, 346: 1311-1320. 10.1126/science.1251385.PubMed CentralPubMedView ArticleGoogle Scholar
- White BJ, Crandall C, Raveche ES, Hjio JH: Laboratory mice carrying three pairs of Robertsonian translocations: establishment of a strain and analysis of meiotic segregation. Cytogenet Cell Genet. 1978, 21: 113-138. 10.1159/000130886.PubMedView ArticleGoogle Scholar
- King M: Species evolution. The Role of Chromosome Change. 1993, Cambridge: Cambridge University PressGoogle Scholar
- Delneri D, Colson I, Grammenoudi S, Roberts IN, Louis EJ, Oliver SG: Engineering evolution to study speciation in yeasts. Nature. 2003, 422: 68-72. 10.1038/nature01418.PubMedView ArticleGoogle Scholar
- Noor MA, Grams KL, Bertucci LA, Reiland J: Chromosomal inversions and the reproductive isolation of species. Proc Natl Acad Sci U S A. 2001, 98: 12084-12088. 10.1073/pnas.221274498.PubMed CentralPubMedView ArticleGoogle Scholar
- Rieseberg LH: Chromosomal rearrangements and speciation. Trends Ecol Evol. 2001, 16: 351-358. 10.1016/S0169-5347(01)02187-5.PubMedView ArticleGoogle Scholar
- Masabanda JS, Burt DW, O'Brien PC, Vignal A, Fillon V, Walsh PS, Cox H, Tempest HG, Smith J, Habermann F, Schmid M, Matsuda Y, Ferguson-Smith MA, Crooijmans RP, Groenen MA, Griffin DK: Molecular cytogenetic definition of the chicken genome: the first complete avian karyotype. Genetics. 2004, 166: 1367-1373. 10.1534/genetics.166.3.1367.PubMed CentralPubMedView ArticleGoogle Scholar
- Christidis L: Aves. Animal Cytogenetics. Volume 4: Chordata 3 B. Edited by: John B, Kayano H, Levan A. 1990, Berlin: Gebrüder BorntraegerGoogle Scholar
- Griffin DK, Haberman F, Masabanda J, O'Brien P, Bagga M, Sazanov A, Smith J, Burt DW, Ferguson-Smith M, Wienberg J: Micro- and macrochromosome paints generated by flow cytometry and microdissection: tools for mapping the chicken genome. Cytogenet Cell Genet. 1999, 87: 278-281. 10.1159/000015449.PubMedView ArticleGoogle Scholar
- Griffin DK, Robertson LB, Tempest HG, Skinner BM: The evolution of the avian genome as revealed by comparative molecular cytogenetics. Cytogenet Genome Res. 2007, 117: 64-77. 10.1159/000103166.PubMedView ArticleGoogle Scholar
- Griffin DK, Robertson LB, Tempest HG, Vignal A, Fillon V, Crooijmans RP, Groenen MA, Deryusheva S, Gaginskaya E, Carré W, Waddington D, Talbot R, Völker M, Masabanda JS, Burt DW: Whole genome comparative studies between chicken and turkey and their implications for avian genome evolution. BMC Genomics. 2008, 9: 168-10.1186/1471-2164-9-168.PubMed CentralPubMedView ArticleGoogle Scholar
- Skinner BM, Völker M, Ellis M, Griffin DK: An appraisal of nuclear organisation in interphase embryonic fibroblasts of chicken, turkey and duck. Cytogenet Genome Res. 2009, 126: 156-164. 10.1159/000245915.PubMedView ArticleGoogle Scholar
- Völker M, Backström N, Skinner BM, Langley EJ, Bunzey SK, Ellegren H, Griffin DK: Copy number variation, chromosome rearrangement, and their association with recombination during avian evolution. Genome Res. 2010, 20: 503-511. 10.1101/gr.103663.109.PubMed CentralPubMedView ArticleGoogle Scholar
- International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.View ArticleGoogle Scholar
- Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Künstner A, Searle S, White S, Vilella AJ, Fairley S, Heger A, Kong L, Ponting CP, Jarvis ED, Mello CV, Minx P, Lovell P, Velho TA, Ferris M, Balakrishnan CN, Sinha S, Blatti C, London SE, Li Y, Lin YC, George J, Sweedler J, Southey B, Gunaratne P, Watson M, et al: The genome of a songbird. Nature. 2010, 464: 757-762. 10.1038/nature08819.PubMed CentralPubMedView ArticleGoogle Scholar
- Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Blomberg LA, Bouffard P, Burt DW, Crasta O, Crooijmans RP, Cooper K, Coulombe RA, De S, Delany ME, Dodgson JB, Dong JJ, Evans C, Frederickson KM, Flicek P, Florea L, Folkerts O, Groenen MA, Harkins TT, Herrero J, Hoffmann S, Megens HJ, Jiang A, de Jong P, Kaiser P, Kim H, et al: Multi-platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): genome assembly and analysis. PLoS Biol. 2010, 8: e1000475-10.1371/journal.pbio.1000475.PubMed CentralPubMedView ArticleGoogle Scholar
- Burt DW, Bruley C, Dunn IC, Jones CT, Ramage A, Law AS, Morrice DR, Paton IR, Smith J, Windsor D, Sazanov A, Fries R, Waddington D: The dynamics of chromosome evolution in birds and mammals. Nature. 1999, 402: 411-413. 10.1038/46555.PubMedView ArticleGoogle Scholar
- Rao M, Morisson M, Faraut T, Bardes S, Fève K, Labarthe E, Fillon V, Huang Y, Li N, Vignal A: A duck RH panel and its potential for assisting NGS genome assembly. BMC Genomics. 2012, 13: 513-10.1186/1471-2164-13-513.PubMed CentralPubMedView ArticleGoogle Scholar
- Skinner BM, Griffin DK: Intrachromosomal rearrangements in avian genome evolution: evidence for regions prone to breakpoints. Heredity (Edinb). 2012, 108: 37-41. 10.1038/hdy.2011.99.View ArticleGoogle Scholar
- Lithgow PE, O’Connor R, Smith D, Fonseka G, Al Mutery A, Rathje C, Frodsham R, O’Brien P, Ferguson-Smith MA, Skinner BM, Griffin DK: Novel tools for characterising inter and intra chromosomal rearrangements in avian microchromosomes. Chromosome Res. 2014, 22: 85-97.PubMedGoogle Scholar
- Larkin DM, Pape G, Donthu R, Auvil L, Welge M, Lewin HA: Breakpoint regions and homologous synteny blocks in chromosomes have different evolutionary histories. Genome Res. 2009, 19: 770-777. 10.1101/gr.086546.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigó R, Hamernik DL, Kappes SM, Lewin HA, Lynn DJ, Nicholas FW, Reymond A, Rijnkels M, Skow LC, Zdobnov EM, Schook L, Womack J, Alioto T, Antonarakis SE, Astashyn A, Chapple CE, Chen HC, Chrast J, Câmara F, Ermolaeva O, Bovine Genome Sequencing and Analysis Consortium, et al: The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009, 324: 522-528.PubMed CentralPubMedView ArticleGoogle Scholar
- Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, Rogel-Gaillard C, Park C, Milan D, Megens HJ, Li S, Larkin DM, Kim H, Frantz LA, Caccamo M, Ahn H, Aken BL, Anselmo A, Anthon C, Auvil L, Badaoui B, Beattie CW, Bendixen C, Berman D, Blecha F, Blomberg J, Bolund L, Bosse M, Botti S, Bujie Z, et al: Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012, 491: 393-398. 10.1038/nature11622.PubMed CentralPubMedView ArticleGoogle Scholar
- Sankoff D: Genome rearrangement with gene families. Bioinformatics. 1999, 15: 909-917. 10.1093/bioinformatics/15.11.909.PubMedView ArticleGoogle Scholar
- Stankiewicz P, Lupski JR: Molecular-evolutionary mechanisms for genomic disorders. Curr Opin Genet Dev. 2002, 12: 312-319. 10.1016/S0959-437X(02)00304-0.PubMedView ArticleGoogle Scholar
- Dumont BL, Payseur BA: Evolution of the genomic recombination rate in murid rodents. Genetics. 2011, 187: 643-657. 10.1534/genetics.110.123851.PubMed CentralPubMedView ArticleGoogle Scholar
- Garcia-Cruz R, Pacheco S, Brieno MA, Steinberg ER, Mudry MD, Ruiz-Herrera A, Garcia-Caldes M: A comparative study of the recombination pattern in three species of Platyrrhini monkeys (primates). Chromosoma. 2011, 120: 521-530. 10.1007/s00412-011-0329-6.PubMedView ArticleGoogle Scholar
- Segura J, Ferretti L, Ramos-Onsins S, Capilla L, Farré M, Reis F, Oliver-Bonet M, Fernández-Bellón H, Garcia F, Garcia-Caldés M, Robinson TJ, Ruiz-Herrera A: Evolution of recombination in eutherian mammals: insights into mechanisms that affect recombination rates and crossover interference. Proc Biol Sci. 2013, 280: 20131945-10.1098/rspb.2013.1945.PubMed CentralPubMedView ArticleGoogle Scholar
- Ranz JM, Maurin D, Chan YS, von Grotthuss M, Hillier LW, Roote J, Ashburner M, Bergman CM: Principles of genome evolution in the Drosophila melanogaster species group. PLoS Biol. 2007, 5: e152-10.1371/journal.pbio.0050152.PubMed CentralPubMedView ArticleGoogle Scholar
- Stevison LS, Hoehn KB, Noor MA: Effects of inversions on within- and between-species recombination and divergence. Genome Biol Evol. 2011, 3: 830-841. 10.1093/gbe/evr081.PubMed CentralPubMedView ArticleGoogle Scholar
- Eichler EE, Sankoff D: Structural dynamics of eukaryotic chromosome evolution. Science. 2003, 301: 793-797. 10.1126/science.1086132.PubMedView ArticleGoogle Scholar
- Singer GA, Lloyd AT, Huminiecki LB, Wolfe KH: Clusters of co-expressed genes in mammalian genomes are conserved by natural selection. Mol Biol Evol. 2005, 22: 767-775.PubMedView ArticleGoogle Scholar
- Sémon M, Duret L: Evolutionary origin and maintenance of coexpressed gene clusters in mammals. Mol Biol Evol. 2006, 23: 1715-1723. 10.1093/molbev/msl034.PubMedView ArticleGoogle Scholar
- Jarvis ED, Mirarab S, Aberer AJ, Li B, Houde P, Li C, Ho SYW, Faircloth BC, Nabholz B, Howard JT, Suh A, Weber CC, da Fonseca RR, Li J, Zhang F, Li H, Zhou L, Narula N, Liu L, Ganapathy G, Boussau B, Bayzid MS, Zavidovych V, Subramanian S, Gabaldón T, Capella-Gutiérrez S, Huerta-Cepas J, Rekepalli B, Munch K, Schierup M, et al: Whole-genome analyses resolve the early branches in the tree of life of modern birds. Science. 2014, 346: 1320-1331. 10.1126/science.1253451.PubMed CentralPubMedView ArticleGoogle Scholar
- Evolution Highway. [http://evolutionhighway.ncsa.uiuc.edu]
- Farré M, Narayan J, Slavov G, Auvil L, Li C, Jarvis ED, Burt DW, Griffin DK, Larkin DM: Chromosome dynamics is associated with ancestral and lineage-specific phenotypes in birds, archosaurians, and other reptiles. Proc Natl Acad Sci U S A.
- Murphy WJ, Larkin DM, Everts van der Wind A, Bourque G, Tesler G, Auvil L, Beever JE, Chowdhary BP, Galibert F, Gatzke L, Hitte C, Meyers SN, Milan D, Ostrander EA, Pape G, Parker HG, Raudsepp T, Rogatcheva MB, Schook LB, Skow LC, Welge M, Womack JE, O'Brien SJ, Pevzner PA, Lewin HA: Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps. Science. 2005, 309: 613-617. 10.1126/science.1111387.PubMedView ArticleGoogle Scholar
- Wallis JW, Aerts J, Groenen MA, Crooijmans RP, Layman D, Graves TA, Scheer DE, Kremitzki C, Fedele MJ, Mudd NK, Cardenas M, Higginbotham J, Carter J, McGrane R, Gaige T, Mead K, Walker J, Albracht D, Davito J, Yang SP, Leong S, Chinwalla A, Sekhon M, Wylie K, Dodgson J, Romanov MN, Cheng H, de Jong PJ, Osoegawa K, Nefedov M, et al: A physical map of the chicken genome. Nature. 2004, 432: 761-764. 10.1038/nature03030.PubMedView ArticleGoogle Scholar
- Romanov MN, Dodgson JB: Cross-species overgo hybridization and comparative physical mapping within avian genomes. Animal Genet. 2006, 37: 397-399. 10.1111/j.1365-2052.2006.01463.x.PubMedView ArticleGoogle Scholar
- Kim J, Larkin DM, Cai Q, Cai Q, Zhang Y, Ge R-L, Auvil L, Capitanu B, Zhang G, Lewin HA, Ma J: Reference-assisted chromosome assembly. Proc Natl Acad Sci U S A. 2013, 110: 1785-1790. 10.1073/pnas.1220349110.PubMed CentralPubMedView ArticleGoogle Scholar
- Shetty S, Griffin DK, Graves JA: Comparative painting reveals strong chromosome homology over 80 million years of bird evolution. Chromosome Res. 1999, 7: 289-295. 10.1023/A:1009278914829.PubMedView ArticleGoogle Scholar
- Guttenbach M, Nanda I, Feichtinger W, Masabanda JS, Griffin DK, Schmid M: Comparative chromosome painting of chicken autosomal paints 1–9 in nine different bird species. Cytogenet Genome Res. 2003, 103: 173-184. 10.1159/000076309.PubMedView ArticleGoogle Scholar
- Nishida-Umehara C, Tsuda Y, Ishijima J, Ando J, Fujiwara A, Matsuda Y, Griffin DK: The molecular basis of chromosome orthologies and sex chromosomal differentiation in palaeognathous birds. Chromosome Res. 2007, 15: 721-734. 10.1007/s10577-007-1157-7.PubMedView ArticleGoogle Scholar
- Morescalchi A: Phylogenetic aspects of karyological evidence. Major Patterns in Vertebrate Evolution. Edited by: Hecht MK, Goody PC, Hecht BM. 1977, New York, London: Plenum Press, 149-167.View ArticleGoogle Scholar
- Morescalchi A, Odierna G, Olmo E: Karyology of the primitive salamanders, family Hynobiidae. Experientia. 1979, 35: 1434-1436. 10.1007/BF01962768.PubMedView ArticleGoogle Scholar
- Mengden GA, Stock AD: Chromosomal evolution in Serpentes: a comparison of G and C chromosome banding patterns of some colubrid and boid genera. Chromosoma. 1980, 79: 53-64. 10.1007/BF00328472.View ArticleGoogle Scholar
- King M, Honeycutt R, Contreras N: Chromosomal repatterning in crocodiles: C, G and N-banding and the in situ hybridization of 18S and 26S rRNA cistrons. Genetica. 1986, 70: 191-201. 10.1007/BF00122186.View ArticleGoogle Scholar
- Burt DW: Origin and evolution of avian microchromosomes. Cytogenet Genome Res. 2002, 96: 97-112. 10.1159/000063018.PubMedView ArticleGoogle Scholar
- Gregory TR: A bird’s-eye view of the C-value enigma: genome size, cell size, and metabolic rate in the class Aves. Evolution. 2002, 56: 121-130. 10.1111/j.0014-3820.2002.tb00854.x.PubMedView ArticleGoogle Scholar
- Andrews CB, Mackenzie SA, Gregory TR: Genome size and wing parameters in passerine birds. Proc Biol Sci. 2009, 276: 55-61. 10.1098/rspb.2008.1012.PubMed CentralPubMedView ArticleGoogle Scholar
- Bush GL, Case SM, Wilson AC, Patton JL: Rapid speciation and chromosomal evolution in mammals. Proc Natl Acad Sci U S A. 1977, 74: 3942-3946. 10.1073/pnas.74.9.3942.PubMed CentralPubMedView ArticleGoogle Scholar
- Fontdevila A, Ruiz A, Ocaña J, Alonso G: Evolutionary history of Drosophila buzzatii. II. How much has chromosomal polymorphism changed in colonization?. Evolution. 1982, 36: 843-851. 10.2307/2407897.View ArticleGoogle Scholar
- Groenen MA, Wahlberg P, Foglio M, Cheng HH, Megens HJ, Crooijmans RP, Besnier F, Lathrop M, Muir WM, Wong GK, Gut I, Andersson L: A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate. Genome Res. 2009, 19: 510-519.PubMed CentralPubMedView ArticleGoogle Scholar
- Backström N, Forstmeier W, Schielzeth H, Mellenius H, Nam K, Bolund E, Webster MT, Öst T, Schneider M, Kempenaers B, Ellegren H: The recombination landscape of the zebra finch Taeniopygia guttata genome. Genome Res. 2010, 20: 485-495. 10.1101/gr.101410.109.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou Q, Zhang J, Bachtrog D, An N, Huang Q, Jarvis ED, Gilbert MT, Zhang G: Complex evolutionary trajectories of sex chromosomes across bird taxa. Science. 2014, 346: 1246338-10.1126/science.1246338.PubMedView ArticleGoogle Scholar
- Witmer LM: The debate on avian ancestry: phylogeny, function, and fossils. Mesozoic Birds: Above the Heads of Dinosaurs. Edited by: Chiappe LM, Witmer LM. 2002, Berkeley: University of California Press, 3-30.Google Scholar
- Chiappe LM, Dyke GJ: The early evolutionary history of birds. J Paleont Soc Korea. 2006, 22: 133-151.Google Scholar
- You HL, Lamanna MC, Harris JD, Chiappe LM, O'Connor J, Ji SA, Lü JC, Yuan CX, Li DQ, Zhang X, Lacovara KJ, Dodson P, Ji Q: A nearly modern amphibious bird from the Early Cretaceous of northwestern China. Science. 2006, 312: 1640-1643. 10.1126/science.1126377.PubMedView ArticleGoogle Scholar
- Nudds RL, Atterholt J, Wang X, You HL, Dyke GJ: Locomotory abilities and habitat of the Cretaceous bird Gansus yumenensis inferred from limb length proportions. J Evol Biol. 2013, 26: 150-154. 10.1111/jeb.12036.PubMedView ArticleGoogle Scholar
- Clarke JA, Tambussi CP, Noriega JI, Erickson GM, Ketcham RA: Definitive fossil evidence for the extant avian radiation in the Cretaceous. Nature. 2005, 433: 305-308. 10.1038/nature03150.PubMedView ArticleGoogle Scholar
- Gatesy SM: Locomotor evolution on the line to modern birds. Mesozoic Birds: Above the Heads of Dinosaurs. Edited by: Chiappe LM, Witmer LM. 2002, Berkeley: University of California Press, 432-447.Google Scholar
- Huang Y, Li Y, Burt DW, Chen H, Zhang Y, Qian W, Kim H, Gan S, Zhao Y, Li J, Yi K, Feng H, Zhu P, Li B, Liu Q, Fairley S, Magor KE, Du Z, Hu X, Goodman L, Tafer H, Vignal A, Lee T, Kim KW, Sheng Z, An Y, Searle S, Herrero J, Groenen MA, Crooijmans RP, et al: The duck genome and transcriptome provide insight into an avian influenza virus reservoir species. Nat Genet. 2013, 45: 776-783. 10.1038/ng.2657.PubMed CentralPubMedView ArticleGoogle Scholar
- MGRA (Multiple Genome Rearrangements and Ancestors) web server, beta version. [http://mgra.bioinf.spbau.ru/]
- Alekseyev MA, Pevzner PA: Breakpoint graphs and ancestral genome reconstructions. Genome Res. 2009, 19: 943-957. 10.1101/gr.082784.108.PubMed CentralPubMedView ArticleGoogle Scholar
- Nanda I, Karl E, Volobouev V, Griffin DK, Schartl M, Schmid M: Extensive gross genomic rearrangements between chicken and Old World vultures (Falconiformes: Accipitridae). Cytogenet Genome Res. 2006, 112: 286-295. 10.1159/000089883.PubMedView ArticleGoogle Scholar
- Nishida C, Ishijima J, Kosaka A, Tanabe H, Habermann FA, Griffin DK, Matsuda Y: Characterization of chromosome structures of Falconinae (Falconidae, Falconiformes, Aves) by chromosome painting and delineation of chromosome rearrangements during their differentiation. Chromosome Res. 2008, 16: 171-181. 10.1007/s10577-007-1210-6.PubMedView ArticleGoogle Scholar
- MGR: Multiple Genome Rearrangements. [http://grimm.ucsd.edu/MGR/]
- Bourque G, Pevzner PA: Genome-Scale Evolution: Reconstructing Gene Orders in the Ancestral Species. Genome Res. 2002, 12: 26-36.PubMed CentralPubMedGoogle Scholar
- Ma J, Zhang L, Suh BB, Raney BJ, Burhans RC, Kent JW, Blanchette M, Haussler D, Miller W: Reconstructing contiguous regions of an ancestral genome. Genome Res. 2006, 16: 1557-1565. 10.1101/gr.5383506.PubMed CentralPubMedView ArticleGoogle Scholar
- The R Project for Statistical Computing. [http://www.r-project.org/]
- DAVID Functional Annotation Bioinformatics Microarray Analysis. [http://david.abcc.ncifcrf.gov/]
- Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nature Protoc. 2009, 4: 44-57.View ArticleGoogle Scholar
- Gene Ontology Enrichment Analysis Software Toolkit (GOEAST). [http://omicslab.genetics.ac.cn/GOEAST/]
- Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36 (Web Server issue): W358-W363.PubMed CentralPubMedView ArticleGoogle Scholar
- BioMart – MartView. [http://www.biomart.org/biomart/martview]
- Kasprzyk A: BioMart: driving a paradigm change in biological data management. Database (Oxford). 2011, 2011: bar049-View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.