The evolution of the avian genome as revealed by comparative molecular cytogenetics

Interbreed variation in meiotic recombination rate and distribution in the domestic chicken Gallus gallus

The efficiency of natural and artificial selection is critically dependent on the recombination rate. However, interbreed and individual variation in recombination rate in poultry remains unknown. Conventional methods of analysis of recombination such as genetic linkage analysis, sperm genotyping and chiasma count at lampbrush chromosomes are expensive and time-consuming. In this study, we analyzed the number and distribution of recombination nodules in spermatocytes of the roosters of six chicken breeds using immunolocalization of key proteins involved in chromosome pairing and recombination. We revealed significant effects of breed (R 2 = 0.17 ; p < 0.001 ) and individual (R 2 = 0.28 ; p < 0.001 ) on variation in the number of recombination nodules. Both interbreed and individual variations in recombination rate were almost entirely determined by variation in recombination density on macrochromosomes, because almost all microchromosomes in each breed had one recombination nodule. Despite interbreed differences in the density of recombination nodules, the patterns of their distribution along homologous chromosomes were similar. The breeds examined in this study showed a correspondence between the age of the breed and its recombination rate. Those with high recombination rates (Pervomai, Russian White and Brahma) are relatively young breeds created by crossing several local breeds. The breeds displaying low recombination rate are ancient local breeds: Cochin (Indo-China), Brown Leghorn (Tuscany, Italy) and Russian Crested (the European part of Russia).

Evolution of Bird Sex Chromosomes Narrated by Repetitive Sequences: Unusual W Chromosome Enlargement in Gallinula melanops (Aves: Gruiformes: Rallidae)

Among birds, species with the ZZ/ZW sex determination system generally show significant differences in morphology and size between the Z and W chromosomes (with the W usually being smaller than the Z). In the present study, we report for the first time the karyotype of the spot-flanked gallinule (Gallinula melanops) by means of classical and molecular cytogenetics. The spot-flanked gallinule has 2n = 80 (11 pairs of macrochromosomes and 29 pairs of microchromosomes) with an unusual W chromosome that is larger than the Z. Besides being totally heterochromatic, it has a secondary constriction in its long arm corresponding to the nucleolar organizer region, as confirmed by both silver staining and mapping of 18S rDNA probes. This is an unprecedented fact among birds. Additionally, 18S rDNA sites were also observed in 6 microchromosomes, while 5S rDNA was found in just 1 microchromosomal pair. Seven out of the 11 used microsatellite sequences were found to be accumulated in microchromosomes, and 6 microsatellite sequences were found in the W chromosome. In addition to the involvement of heterochromatin and repetitive DNAs in the differentiation of the large W chromosome, the results also show an alternative scenario that highlights the plasticity that shapes the evolutionary history of bird sex chromosomes.

Banding cytogenetics of the vulnerable species Houbara bustard (Otidiformes) and comparative analysis with the Domestic fowl

The Houbara bustard Chlamydotisundulata (Jacquin, 1784) is an emblematic and endangered bird of steppes and desert spaces of North Africa. This species belonging to Otidiformes is recognized as vulnerable by the International Union for Nature Conservation. The critical situation of this species and the revision of its classification on the tree of birds encouraged the authors to start accumulating chromosome data. For that, we propose the GTG- and RBG-banded karyotypes of the Houbara bustard prepared from primary fibroblast cell cultures. The first eight autosomal pairs and sex chromosomes have been described and compared to those of the domestic fowl Gallusdomesticus (Linnaeus, 1758). The diploid number has been estimated as 78 chromosomes with 8 macrochromosomes pairs and 30 microchromosomes pairs, attesting of the stability of chromosome number in avian karyotypes. The description of the karyotype of the Houbara is of crucial importance for the management of the reproduction of this species in captivity. It can be used as a reference in the detection of chromosomal abnormalities, which would be responsible of the early embryonic mortalities.

Determinants of the Efficacy of Natural Selection on Coding and Noncoding Variability in Two Passerine Species

Population genetic theory predicts that selection should be more effective when the effective population size (N_e) is larger, and that the efficacy of selection should correlate positively with recombination rate. Here, we analyzed the genomes of ten great tits and ten zebra finches. Nucleotide diversity at 4-fold degenerate sites indicates that zebra finches have a 2.83-fold larger N_e. We obtained clear evidence that purifying selection is more effective in zebra finches. The proportion of substitutions at 0-fold degenerate sites fixed by positive selection (α) is high in both species (great tit 48%; zebra finch 64%) and is significantly higher in zebra finches. When α was estimated on GC-conservative changes (i.e., between A and T and between G and C), the estimates reduced in both species (great tit 22%; zebra finch 53%). A theoretical model presented herein suggests that failing to control for the effects of GC-biased gene conversion (gBGC) is potentially a contributor to the overestimation of α, and that this effect cannot be alleviated by first fitting a demographic model to neutral variants. We present the first estimates in birds for α in the untranslated regions, and found evidence for substantial adaptive changes. Finally, although purifying selection is stronger in high-recombination regions, we obtained mixed evidence for α increasing with recombination rate, especially after accounting for gBGC. These results highlight that it is important to consider the potential confounding effects of gBGC when quantifying selection and that our understanding of what determines the efficacy of selection is incomplete.

Patterns of microchromosome organization remain highly conserved throughout avian evolution

The structure and organization of a species genome at a karyotypic level, and in interphase nuclei, have broad functional significance. Although regular sized chromosomes are studied extensively in this regard, microchromosomes, which are present in many terrestrial vertebrates, remain poorly explored. Birds have more cytologically indistinguishable microchromosomes (~ 30 pairs) than other vertebrates; however, the degree to which genome organization patterns at a karyotypic and interphase level differ between species is unknown. In species where microchromosomes have fused to other chromosomes, they retain genomic features such as gene density and GC content; however, the extent to which they retain a central nuclear position has not been investigated. In studying 22 avian species from 10 orders, we established that, other than in species where microchromosomal fusion is obvious (Falconiformes and Psittaciformes), there was no evidence of microchromosomal rearrangement, suggesting an evolutionarily stable avian genome (karyotypic) organization. Moreover, in species where microchromosomal fusion has occurred, they retain a central nuclear location, suggesting that the nuclear position of microchromosomes is a function of their genomic features rather than their physical size.

Chromosomal polymorphism and comparative chromosome painting in the rufous-collared sparrow (Zonotrichia capensis)

Zonotrichia capensis is widely distributed in the Neotropics. Previous cytogenetic studies demonstrated the presence of polymorphisms in two chromosome pairs (ZCA2 and ZCA4). Here, we report results based on comparative chromosome painting, using probes derived from Gallus gallus and Leucopternis albicollis, focused on characterizing the chromosome organization of Z. capensis. Our results demonstrate the conservation of ancestral syntenies as observed previously in other species of passerine. Syntenies were rearranged by a series of inversions in the second chromosome as described in other Passeriformes, but in this species, by using probes derived from L. albicollis we observed an extra inversion in the second chromosome that had not previously been reported. We also report a paracentric inversion in pair 3; this chromosome corresponds to the second chromosome in Zonotrichia albicollis and may indicate the presence of ancestral inversions in the genus. The chromosomal inversions we found might be important for understanding the phenotypic variation that exists throughout the distribution of Z. capensis.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes

BACKGROUND

The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions.

RESULTS

We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest.

CONCLUSIONS

The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.

Chromosome Level Genome Assembly and Comparative Genomics between Three Falcon Species Reveals an Unusual Pattern of Genome Organisation

Whole genome assemblies are crucial for understanding a wide range of aspects of falcon biology, including morphology, ecology, and physiology, and are thus essential for their care and conservation. A key aspect of the genome of any species is its karyotype, which can then be linked to the whole genome sequence to generate a so-called chromosome-level assembly. Chromosome-level assemblies are essential for marker assisted selection and genotype-phenotype correlations in breeding regimes, as well as determining patterns of gross genomic evolution. To date, only two falcon species have been sequenced and neither initially were assembled to the chromosome level. Falcons have atypical avian karyotypes with fewer chromosomes than other birds, presumably brought about by wholesale fusion. To date, however, published chromosome preparations are of poor quality, few chromosomes have been distinguished and standard ideograms have not been made. The purposes of this study were to generate analyzable karyotypes and ideograms of peregrine, saker, and gyr falcons, report on our recent generation of chromosome level sequence assemblies of peregrine and saker falcons, and for the first time, sequence the gyr falcon genome. Finally, we aimed to generate comparative genomic data between all three species and the reference chicken genome. Results revealed a diploid number of 2n = 50 for peregrine falcon and 2n = 52 for saker and gyr through high quality banded chromosomes. Standard ideograms that are generated here helped to map predicted chromosomal fragments (PCFs) from the genome sequences directly to chromosomes and thus generate chromosome level sequence assemblies for peregrine and saker falcons. Whole genome sequencing was successful in gyr falcon, but read depth and coverage was not sufficient to generate a chromosome level assembly. Nonetheless, comparative genomics revealed no differences in genome organization between gyr and saker falcons. When compared to peregrine falcon, saker/gyr differed by one interchromosomal and seven intrachromosomal rearrangements (a fusion plus seven inversions), whereas peregrine and saker/gyr differ from the reference chicken genome by 14/13 fusions (11 microchromosomal) and six fissions. The chromosomal differences between the species could potentially provide the basis of a screening test for hybrid animals.

Banding cytogenetics of chimeric hybrids Coturnixcoturnix × Coturnixjaponica and comparative analysis with the domestic fowl

The Common quail Coturnixcoturnix Linnaeus, 1758 is a wild migratory bird which is distributed in Eurasia and North Africa, everywhere with an accelerating decline in population size. This species is protected by the Bonn and Berne conventions (1979) and by annex II/1 of the Birds Directive (2009). In Algeria, its breeding took place at the hunting centre in the west of the country. Breeding errors caused uncontrolled crosses between the Common quail and Japanese quail Coturnixjaponica Temminck & Schlegel, 1849. In order to help to preserve the natural genetic heritage of the Common quail and to lift the ambiguity among the populations of quail raised in Algeria, it seemed essential to begin to describe the chromosomes of this species in the country since no cytogenetic study has been reported to date. Fibroblast cultures from embryo and adult animal were initiated. Double synchronization with excess thymidine allowed us to obtain high resolution chromosomes blocked at prometaphase stage. The karyotype and the idiogram in GTG morphological banding (G-bands obtained with trypsin and Giemsa) corresponding to larger chromosomes 1-12 and ZW pair were thus established. The diploid set of chromosomes was estimated as 2N=78. Cytogenetic analysis of expected hybrid animals revealed the presence of a genetic introgression and cellular chimerism. This technique is effective in distinguishing the two quail taxa. Furthermore, the comparative chromosomal analysis of the two quails and domestic chicken Gallusgallusdomesticus Linnaeus, 1758 has been conducted. Differences in morphology and/or GTG band motifs were observed on 1, 2, 4, 7, 8 and W chromosomes. Neocentromere occurrence was suggested for Common quail chromosome 1 and Chicken chromosomes 4 and W. Double pericentric inversion was observed on the Common quail chromosome 2 while pericentric inversion hypothesis was proposed for Chicken chromosome 8. A deletion on the short arm of the Common quail chromosome 7 was also found. These results suggest that Common quail would be a chromosomally intermediate species between Chicken and Japanese quail. The appearance of only a few intrachromosomal rearrangements that occurred during evolution suggests that the organization of the genome is highly conserved between these three galliform species.

Chromosome Painting in Neotropical Long- and Short-Tailed Parrots (Aves, Psittaciformes): Phylogeny and Proposal for a Putative Ancestral Karyotype for Tribe Arini

Most Neotropical Psittacidae have a diploid number of 2n = 70, and a dichotomy in chromosome patterns. Long-tailed species have biarmed macrochromosomes, while short-tailed species have telo/acrocentric macrochromosomes. However, the use of chromosome painting has demonstrated that karyotype evolution in Psittacidae includes a high number of inter/intrachromosomal rearrangements. To determine the phylogeny of long- and short-tailed species, and to propose a putative ancestral karyotype for this group, we constructed homology maps of Pyrrhura frontalis (PFR) and Amazona aestiva (AAE), belonging to the long- and short-tailed groups, respectively. Chromosomes were analyzed by conventional staining and fluorescent in situ hybridization using whole chromosome paints of Gallusgallus and Leucopternis albicollis. Conventional staining showed a karyotype with 2n = 70 in both species, with biarmed macrochromosomes in PFR and telo/acrocentric chromosomes in AAE. Comparison of the results with the putative avian ancestral karyotype (PAK) showed fusions in PFR of PAK1p/PAK4q (PFR1) and PAK6/PAK7 (PFR6) with a paracentric inversion in PFR6. However, in AAE, there was only the fusion between PAK6/7 (AAE7) with a paracentric inversion. Our results indicate that PFR retained a more basal karyotype than long-tailed species previously studied, and AAE a more basal karyotype for Neotropical Psittacidae analyzed so far.

Reconstruction of avian ancestral karyotypes reveals differences in the evolutionary history of macro- and microchromosomes

Background

Reconstruction of ancestral karyotypes is critical for our understanding of genome evolution, allowing for the identification of the gross changes that shaped extant genomes. The identification of such changes and their time of occurrence can shed light on the biology of each species, clade and their evolutionary history. However, this is impeded by both the fragmented nature of the majority of genome assemblies and the limitations of the available software to work with them. These limitations are particularly apparent in birds, with only 10 chromosome-level assemblies reported thus far. Algorithmic approaches applied to fragmented genome assemblies can nonetheless help define patterns of chromosomal change in defined taxonomic groups.

Results

Here, we make use of the DESCHRAMBLER algorithm to perform the first large-scale study of ancestral chromosome structure and evolution in birds. This algorithm allows us to reconstruct the overall genome structure of 14 key nodes of avian evolution from the Avian ancestor to the ancestor of the Estrildidae, Thraupidae and Fringillidae families.

Conclusions

Analysis of these reconstructions provides important insights into the variability of rearrangement rates during avian evolution and allows the detection of patterns related to the chromosome distribution of evolutionary breakpoint regions. Moreover, the inclusion of microchromosomes in our reconstructions allows us to provide novel insights into the evolution of these avian chromosomes, specifically.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1544-8) contains supplementary material, which is available to authorized users.

Inversions are bigger on the X chromosome

In many insects, X-linked inversions fix at a higher rate and are much less polymorphic than autosomal inversions. Here, we report that in Drosophila, X-linked inversions also capture 67% more genes. We estimated the number of genes captured through an approximate Bayesian computational analysis of gene orders in nine species of Drosophila. X-linked inversions fixed with a significantly larger gene content. Further, X-linked inversions of intermediate size enjoy highest fixation rate, while the fixation rate of autosomal inversions decreases with size. A less detailed analysis in Anopheles suggests a similar pattern holds in mosquitoes. We develop a population genetic model that assumes the fitness effects of inversions scale with the number of genes captured. We show that the same conditions that lead to a higher fixation rate also produce a larger size for inversions on the X.

Chromosome painting in Glyphorynchus spirurus (Vieillot, 1819) detects a new fission in Passeriformes

Glyphorynchus spirurus (GSP), also called the Wedge-billed Woodcreeper (Furnariidae) has an extensive distribution in the Americas, including the Atlantic coast of Brazil. Nevertheless, there is no information about its karyotype or genome organization. To contribute to the knowledge of chromosomal evolution in Passeriformes we analysed the karyotype of Glyphorynchus spirurus by classic and molecular cytogenetics methods. We show that Glyphorynchus spirurus has a 2n = 80 karyotype with a fundamental number (FN) of 84, similar to the avian putative ancestral karyotype (PAK). Glyphorynchus spirurus pair 1 was heteromorphic in the Tapajós population whereby the short arms varied in sizes, possibly due to a pericentric inversion, as described in other Furnariidae birds. FISH with the Histone H5 probe revealed a signal in the pericentromeric region of G. spirurus chromosome 5 and rDNA 18S showed interstitial signal in GSP-1. Chromosome painting with Gallus gallus (GGA) macrochromosomes 1-9 probes showed disruption of chromosome syntenies of GGA-1, 2 and 4 by fission in Glyphorynchus spirurus. Our results confirm that the GGA1 centric fission is a synapomorphic character for the phylogenetic branch composed of Strigiformes, Passeriformes, Columbiformes and Falconiformes. On the other hand, the GGA-2 fission is reported here for the first time in Passeriformes. Chromosome painting with BOE whole chromosome probes confirmed these rearrangements in Glyphorynchus spirurus revealed by Gallus gallus 1-9 probes, in addition to enabling the establishment of genome-wide homology map.

A Genetic Map of Ostrich Z Chromosome and the Role of Inversions in Avian Sex Chromosome Evolution

Recombination arrest is a necessary step for the evolution of distinct sex chromosomes. Structural changes, such as inversions, may represent the mechanistic basis for recombination suppression and comparisons of the structural organization of chromosomes as given by chromosome-level assemblies offer the possibility to infer inversions across species at some detail. In birds, deduction of the process of sex chromosome evolution has been hampered by the lack of a validated chromosome-level assembly from a representative of one of the two basal clades of modern birds, Paleognathae. We therefore developed a high-density genetic linkage map of the ostrich Z chromosome and used this to correct an existing assembly, including correction of a large chimeric superscaffold and the order and orientation of other superscaffolds. We identified the pseudoautosomal region as a 52 Mb segment (≈60% of the Z chromosome) where recombination occurred in both sexes. By comparing the order and location of genes on the ostrich Z chromosome with that of six bird species from the other major clade of birds (Neognathae), and of reptilian outgroup species, 25 Z-linked inversions were inferred in the avian lineages. We defined Z chromosome organization in an early avian ancestor and identified inversions spanning the candidate sex-determining DMRT1 gene in this ancestor, which could potentially have triggered the onset of avian sex chromosome evolution. We conclude that avian sex chromosome evolution has been characterized by a complex process of probably both Z-linked and W-linked inversions (and/or other processes). This study illustrates the need for validated chromosome-level assemblies for inference of genome evolution.

Comparative chromosome painting in Columbidae (Columbiformes) reinforces divergence in Passerea and Columbea

Pigeons and doves (Columbiformes) are one of the oldest and most diverse extant lineages of birds. However, the karyotype evolution within Columbiformes remains unclear. To delineate the synteny-conserved segments and karyotypic differences among four Columbidae species, we used chromosome painting from Gallus gallus (GGA, 2n = 78) and Leucopternis albicollis (LAL, 2n = 68). Besides that, a set of painting probes for the eared dove, Zenaida auriculata (ZAU, 2n = 76), was generated from flow-sorted chromosomes. Chromosome painting with GGA and ZAU probes showed conservation of the first ten ancestral pairs in Z. auriculata, Columba livia, and Columbina picui, while in Leptotila verreauxi, fusion of the ancestral chromosomes 6 and 7 was observed. However, LAL probes revealed a complex reorganization of ancestral chromosome 1, involving paracentric and pericentric inversions. Because of the presence of similar intrachromosomal rearrangements in the chromosomes corresponding to GGA1q in the Columbidae and Passeriformes species but without a common origin, these results are consistent with the recent proposal of divergence within Neoaves (Passerea and Columbea). In addition, inversions in chromosome 2 were identified in C. picui and L. verreauxi. Thus, in four species of distinct genera of the Columbidae family, unique chromosomal rearrangements have occurred during karyotype evolution, confirming that despite conservation of the ancestral syntenic groups, these chromosomes have been modified by the occurrence of intrachromosomal rearrangements.

Banding cytogenetics of the Barbary partridge Alectoris barbara and the Chukar partridge Alectoris chukar (Phasianidae): a large conservation with Domestic fowl Gallus domesticus revealed by high resolution chromosomes

The development of avian cytogenetics is significantly behind that of mammals. In fact, since the advent of cytogenetic techniques, fewer than 1500 karyotypes have been established. The Barbary partridge Alectoris barbara Bonnaterre, 1790 is a bird of economic interest but its genome has not been studied so far. This species is endemic to North Africa and globally declining. The Chukar partridge Alectoris chukar Gray, 1830 is an introduced species which shares the same habitat area as the Barbary partridge and so there could be introgressive hybridisation. A cytogenetic study has been initiated in order to contribute to the Barbary partridge and the Chukar partridge genome analyses. The GTG, RBG and RHG-banded karyotypes of these species have been described. Primary fibroblast cell lines obtained from embryos were harvested after simple and double thymidine synchronisation. The first eight autosomal pairs and Z sex chromosome have been described at high resolution and compared to those of the domestic fowl Gallus domesticus Linnaeus, 1758. The diploid number was established as 2n = 78 for both partridges, as well as for most species belonging to the Galliformes order, underlying the stability of chromosome number in avian karyotypes. Wide homologies were observed for macrochromosomes and gonosome except for chromosome 4, 7, 8 and Z which present differences in morphology and/or banding pattern. Neocentromere occurrence was suggested for both partridges chromosome 4 with an assumed paracentric inversion in the Chukar partridge chromosome 4. Terminal inversion in the long arm of the Barbary partridge chromosome Z was also found. These rearrangements confirm that the avian karyotypes structure is conserved interchromosomally, but not at the intrachromosomal scale.

Karyotype Evolution in Birds: From Conventional Staining to Chromosome Painting

In the last few decades, there have been great efforts to reconstruct the phylogeny of Neoaves based mainly on DNA sequencing. Despite the importance of karyotype data in phylogenetic studies, especially with the advent of fluorescence in situ hybridization (FISH) techniques using different types of probes, the use of chromosomal data to clarify phylogenetic proposals is still minimal. Additionally, comparative chromosome painting in birds is restricted to a few orders, while in mammals, for example, virtually all orders have already been analyzed using this method. Most reports are based on comparisons using Gallus gallus probes, and only a small number of species have been analyzed with more informative sets of probes, such as those from Leucopternis albicollis and Gyps fulvus, which show ancestral macrochromosomes rearranged in alternative patterns. Despite this, it is appropriate to review the available cytogenetic information and possible phylogenetic conclusions. In this report, the authors gather both classical and molecular cytogenetic data and describe some interesting and unique characteristics of karyotype evolution in birds.

Comparative analysis examining patterns of genomic differentiation across multiple episodes of population divergence in birds

Heterogeneous patterns of genomic differentiation are commonly documented between closely related populations and there is considerable interest in identifying factors that contribute to their formation. These factors could include genomic features (e.g., areas of low recombination) that promote processes like linked selection (positive or purifying selection that affects linked neutral sites) at specific genomic regions. Examinations of repeatable patterns of differentiation across population pairs can provide insight into the role of these factors. Birds are well suited for this work, as genome structure is conserved across this group. Accordingly, we reestimated relative (F_ST ) and absolute (d_XY ) differentiation between eight sister pairs of birds that span a broad taxonomic range using a common pipeline. Across pairs, there were modest but significant correlations in window-based estimates of differentiation (up to 3% of variation explained for F_ST and 26% for d_XY ), supporting a role for processes at conserved genomic features in generating heterogeneous patterns of differentiation; processes specific to each episode of population divergence likely explain the remaining variation. The role genomic features play was reinforced by linear models identifying several genomic variables (e.g., gene densities) as significant predictors of F_ST and d_XY repeatability. F_ST repeatability was higher among pairs that were further along the speciation continuum (i.e., more reproductively isolated) providing further insight into how genomic differentiation changes with population divergence; early stages of speciation may be dominated by positive selection that is different between pairs but becomes integrated with processes acting according to shared genomic features as speciation proceeds.

Multidirectional chromosome painting in Synallaxis frontalis (Passeriformes, Furnariidae) reveals high chromosomal reorganization, involving fissions and inversions

In this work we performed comparative chromosome painting using probes from Gallus gallus (GGA) Linnaeus, 1758 and Leucopternis albicollis (LAL) Latham, 1790 in Synallaxis frontalis Pelzeln, 1859 (Passeriformes, Furnariidae), an exclusively Neotropical species, in order to analyze whether the complex pattern of intrachromosomal rearrangements (paracentric and pericentric inversions) proposed for Oscines and Suboscines is shared with more basal species. S. frontalis has 82 chromosomes, similar to most Avian species, with a large number of microchromosomes and a few pairs of macrochromosomes. We found polymorphisms in pairs 1 and 3, where homologues were submetacentric and acrocentric. Hybridization of GGA probes showed syntenies in the majority of ancestral macrochromosomes, except for GGA1 and GGA2, which hybridized to more than one pair of chromosomes each. LAL probes confirmed the occurrence of intrachromosomal rearrangements in the chromosomes corresponding to GGA1q, as previously proposed for species from the order Passeriformes. In addition, LAL probes suggest that pericentric inversions or centromere repositioning were responsible for variations in the morphology of the heteromorphic pairs 1 and 3. Altogether, the analysis of our data on chromosome painting and the data published in other Passeriformes highlights chromosomal changes that have occurred during the evolution of Passeriformes.

Amphibian and Avian Karyotype Evolution: Insights from Lampbrush Chromosome Studies

Amphibian and bird karyotypes typically have a complex organization, which makes them difficult for standard cytogenetic analysis. That is, amphibian chromosomes are generally large, enriched with repetitive elements, and characterized by the absence of informative banding patterns. The majority of avian karyotypes comprise a small number of relatively large macrochromosomes and numerous tiny morphologically undistinguishable microchromosomes. A good progress in investigation of amphibian and avian chromosome evolution became possible with the usage of giant lampbrush chromosomes typical for growing oocytes. Due to the giant size, peculiarities of organization and enrichment with cytological markers, lampbrush chromosomes can serve as an opportune model for comprehensive high-resolution cytogenetic and cytological investigations. Here, we review the main findings on chromosome evolution in amphibians and birds that were obtained using lampbrush chromosomes. In particular, we discuss the data on evolutionary chromosomal rearrangements, accumulation of polymorphisms, evolution of sex chromosomes as well as chromosomal changes during clonal reproduction of interspecies hybrids.

Broad-scale recombination pattern in the primitive bird Rhea americana (Ratites, Palaeognathae)

Birds have genomic and chromosomal features that make them an attractive group to analyze the evolution of recombination rate and the distribution of crossing over. Yet, analyses are biased towards certain species, especially domestic poultry and passerines. Here we analyze for the first time the recombination rate and crossover distribution in the primitive ratite bird, Rhea americana (Rheiformes, Palaeognathae). Using a cytogenetic approach for in situ mapping of crossovers we found that the total genetic map is 3050 cM with a global recombination rate of 2.1 cM/Mb for female rheas. In the five largest macrobivalents there were 3 or more crossovers in most bivalents. Recombination rates for macrobivalents ranges between 1.8-2.1 cM/Mb and the physical length of their synaptonemal complexes is highly predictive of their genetic lengths. The crossover rate at the pseudoautosomal region is 2.1 cM/Mb, similar to those of autosomal pairs 5 and 6 and only slightly higher compared to other macroautosomes. It is suggested that the presence of multiple crossovers on the largest macrobivalents is a feature common to many avian groups, irrespective of their position throughout phylogeny. These data provide new insights to analyze the heterogeneous recombination landscape of birds.

Comparative Cytogenetics between Two Important Songbird, Models: The Zebra Finch and the Canary

Songbird species (order Passeriformes, suborder Oscines) are important models in various experimental fields spanning behavioural genomics to neurobiology. Although the genomes of some songbird species were sequenced recently, the chromosomal organization of these species is mostly unknown. Here we focused on the two most studied songbird species in neuroscience, the zebra finch (Taeniopygia guttata) and the canary (Serinus canaria). In order to clarify these issues and also to integrate chromosome data with their assembled genomes, we used classical and molecular cytogenetics in both zebra finch and canary to define their chromosomal homology, localization of heterochromatic blocks and distribution of rDNA clusters. We confirmed the same diploid number (2n = 80) in both species, as previously reported. FISH experiments confirmed the occurrence of multiple paracentric and pericentric inversions previously found in other species of Passeriformes, providing a cytogenetic signature for this order, and corroborating data from in silico analyses. Additionally, compared to other Passeriformes, we detected differences in the zebra finch karyotype concerning the morphology of some chromosomes, in the distribution of 5S rDNA clusters, and an inversion in chromosome 1.

Upgrading short-read animal genome assemblies to chromosome level using comparative genomics and a universal probe set

Most recent initiatives to sequence and assemble new species’ genomes de novo fail to achieve the ultimate endpoint to produce contigs, each representing one whole chromosome. Even the best-assembled genomes (using contemporary technologies) consist of subchromosomal-sized scaffolds. To circumvent this problem, we developed a novel approach that combines computational algorithms to merge scaffolds into chromosomal fragments, PCR-based scaffold verification, and physical mapping to chromosomes. Multigenome-alignment-guided probe selection led to the development of a set of universal avian BAC clones that permit rapid anchoring of multiple scaffolds to chromosomes on all avian genomes. As proof of principle, we assembled genomes of the pigeon (Columbia livia) and peregrine falcon (Falco peregrinus) to chromosome levels comparable, in continuity, to avian reference genomes. Both species are of interest for breeding, cultural, food, and/or environmental reasons. Pigeon has a typical avian karyotype (2n = 80), while falcon (2n = 50) is highly rearranged compared to the avian ancestor. By using chromosome breakpoint data, we established that avian interchromosomal breakpoints appear in the regions of low density of conserved noncoding elements (CNEs) and that the chromosomal fission sites are further limited to long CNE “deserts.” This corresponds with fission being the rarest type of rearrangement in avian genome evolution. High-throughput multiple hybridization and rapid capture strategies using the current BAC set provide the basis for assembling numerous avian (and possibly other reptilian) species, while the overall strategy for scaffold assembly and mapping provides the basis for an approach that (provided metaphases can be generated) could be applied to any animal genome.

Constitutive heterochromatin in chromosomes of duck hybrids and goose hybrids

Constitutive heterochromatin is a highly condensed fraction of chromatin in chromosomes. It is characterized by a high degree of polymorphism. Heterochromatin is located in the centromeric, telomeric, and interstitial parts of chromosomes. We used the CBG ( C: banding using B: arium hydroxide by G: iemsa) staining technique to identify heterochromatin in chromosomes. Analysis of karyotypes of F1 hybrids resulting from intergeneric hybridization of ducks (A. platyrhynchos × C. moschata) and interspecific crosses of geese (A. anser × A. cygnoides) were used to compare the karyotypes of 2 species of duck and 2 species of geese, as well as to compare the hybrids with the parent species. The localization of C-bands and their size were determined. In the duck hybrid, greater amounts of heterochromatin were noted in the homologous chromosomes from the duck A. platyrhynchos than in the chromosomes from the duck C. moschata. In the goose hybrid more heterochromatin was observed in the homologous chromosomes from the goose A. cygnoides than in the chromosomes from the goose A. anser. Comparison of chromosomes from the duck hybrid with chromosomes of the ducks A. platyrhynchos and C. moschata revealed nearly twice as much constitutive heterochromatin in the chromosomes of the hybrid. When chromosomes from the goose hybrid were compared with those of the geese A. anser and A. cygnoides, differences in the average content of heterochromatin were observed on only a few chromosomes.

Novel Insights into Chromosome Evolution in Birds, Archosaurs, and Reptiles

Homologous synteny blocks (HSBs) and evolutionary breakpoint regions (EBRs) in mammalian chromosomes are enriched for distinct DNA features, contributing to distinct phenotypes. To reveal HSB and EBR roles in avian evolution, we performed a sequence-based comparison of 21 avian and 5 outgroup species using recently sequenced genomes across the avian family tree and a newly-developed algorithm. We identified EBRs and HSBs in ancestral bird, archosaurian (bird, crocodile, and dinosaur), and reptile chromosomes. Genes involved in the regulation of gene expression and biosynthetic processes were preferably located in HSBs, including for example, avian-specific HSBs enriched for genes involved in limb development. Within birds, some lineage-specific EBRs rearranged genes were related to distinct phenotypes, such as forebrain development in parrots. Our findings provide novel evolutionary insights into genome evolution in birds, particularly on how chromosome rearrangements likely contributed to the formation of novel phenotypes.

Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome

Background

Squamates (lizards and snakes) are a speciose lineage of reptiles displaying considerable karyotypic diversity, particularly among lizards. Understanding the evolution of this diversity requires comparison of genome organisation between species. Although the genomes of several squamate species have now been sequenced, only the green anole lizard has any sequence anchored to chromosomes. There is only limited gene mapping data available for five other squamates. This makes it difficult to reconstruct the events that have led to extant squamate karyotypic diversity. The purpose of this study was to anchor the recently sequenced central bearded dragon (Pogona vitticeps) genome to chromosomes to trace the evolution of squamate chromosomes. Assigning sequence to sex chromosomes was of particular interest for identifying candidate sex determining genes.

Results

By using two different approaches to map conserved blocks of genes, we were able to anchor approximately 42 % of the dragon genome sequence to chromosomes. We constructed detailed comparative maps between dragon, anole and chicken genomes, and where possible, made broader comparisons across Squamata using cytogenetic mapping information for five other species. We show that squamate macrochromosomes are relatively well conserved between species, supporting findings from previous molecular cytogenetic studies. Macrochromosome diversity between members of the Toxicofera clade has been generated by intrachromosomal, and a small number of interchromosomal, rearrangements. We reconstructed the ancestral squamate macrochromosomes by drawing upon comparative cytogenetic mapping data from seven squamate species and propose the events leading to the arrangements observed in representative species. In addition, we assigned over 8 Mbp of sequence containing 219 genes to the Z chromosome, providing a list of genes to begin testing as candidate sex determining genes.

Conclusions

Anchoring of the dragon genome has provided substantial insight into the evolution of squamate genomes, enabling us to reconstruct ancestral macrochromosome arrangements at key positions in the squamate phylogeny, demonstrating that fusions between macrochromosomes or fusions of macrochromosomes and microchromosomes, have played an important role during the evolution of squamate genomes. Assigning sequence to the sex chromosomes has identified NR5A1 as a promising candidate sex determining gene in the dragon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2774-3) contains supplementary material, which is available to authorized users.

Marsupials as models for understanding the role of chromosome rearrangements in evolution and disease

Chromosome rearrangements have been implicated in diseases, such as cancer, and speciation, but it remains unclear whether rearrangements are causal or merely a consequence of these processes. Two marsupial families with very different rates of karyotype evolution provide excellent models in which to study the role of chromosome rearrangements in a disease and evolutionary context. The speciose family Dasyuridae displays remarkable karyotypic conservation, with all species examined to date possessing nearly identical karyotypes. Despite the seemingly high degree of chromosome stability within this family, they appear prone to developing tumours, including transmissible devil facial tumours. In contrast, chromosome rearrangements have been frequent in the evolution of the species-rich family Macropodidae, which displays a high level of karyotypic diversity. In particular, the genus Petrogale (rock-wallabies) displays an extraordinary level of chromosome rearrangement among species. For six parapatric Petrogale species, it appears that speciation has essentially been caught in the act, providing an opportunity to determine whether chromosomal rearrangements are a cause or consequence of speciation in this system. This review highlights the reasons that these two marsupial families are excellent models for testing hypotheses for hotspots of chromosome rearrangement and deciphering the role of chromosome rearrangements in disease and speciation.

Gene Regulatory Evolution During Speciation in a Songbird

Over the last decade, tremendous progress has been made toward a comparative understanding of gene regulatory evolution. However, we know little about how gene regulation evolves in birds, and how divergent genomes interact in their hybrids. Because of the unique features of birds – female heterogamety, a highly conserved karyotype, and the slow evolution of reproductive incompatibilities – an understanding of regulatory evolution in birds is critical to a comprehensive understanding of regulatory evolution and its implications for speciation. Using a novel complement of analyses of replicated RNA-seq libraries, we demonstrate abundant divergence in brain gene expression between zebra finch (Taeniopygia guttata) subspecies. By comparing parental populations and their F1 hybrids, we also show that gene misexpression is relatively rare among brain-expressed transcripts in male birds. If this pattern is consistent across tissues and sexes, it may partially explain the slow buildup of postzygotic reproductive isolation observed in birds relative to other taxa. Although we expected that the action of genetic drift on the island-dwelling zebra finch subspecies would be manifested in a higher rate of trans regulatory divergence, we found that most divergence was in cis regulation, following a pattern commonly observed in other taxa. Thus, our study highlights both unique and shared features of avian regulatory evolution.

Cytotaxonomy of Eurypyga helias (Gruiformes, Eurypygidae): First Karyotypic Description and Phylogenetic Proximity with Rynochetidae

The sunbittern (Eurypyga helias) is a South American Gruiformes, the only member of Family Eurypigidae. In most phylogenetic proposals, it is placed in a more distant position than other families of the so-called "core Gruiformes". Different studies based on molecular, morphological and biogeographical data suggest that the Eurypigidae is closely related to the kagu (Rhynochetos jubatus), the only species in Rynochetidae, another family not included in the core Gruiformes. Here, the karyotype of the sunbittern is described for the first time, by classical and molecular cytogenetics, using whole chromosome probes derived from Gallus gallus and Leucopternis albicollis. We found a diploid number of 80, with only one pair of biarmed autosomal macrochromosomes, similar to that observed in the kagu. Chromosome painting revealed that most syntenies found in the avian putative ancestral karyotype (PAK) were conserved in the sunbittern. However, PAK1, PAK2, and PAK5 corresponded to two chromosome pairs each. Probes derived from L. albicollis confirm that fissions in PAK1 and PAK2 were centric, whereas in PAK5 the fission is interstitial. In addition, there is fusion of segments homologous to PAK2q and PAK5. From a phylogenetic point of view, comparisons of our results with two other Gruiformes belonging to family Rallidae suggest that the PAK5q fission might be a synapomorphy for Gruiformes. Fissions in PAK1 and PAK2 are found only in Eurypigidae, and might also occur in Rynochetidae, in view of the similar chromosomal morphology between the sunbittern and the kagu. This suggests a close phylogenetic relationship between Eurypigidae and Rynochetidae, whose common ancestor was separated by the Gondwana vicariancy in South America and New Caledonia, respectively.

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.

High-density sex-specific linkage maps of a European tree frog (Hyla arborea) identify the sex chromosome without information on offspring sex

Identifying homology between sex chromosomes of different species is essential to understanding the evolution of sex determination. Here, we show that the identity of a homomorphic sex chromosome pair can be established using a linkage map, without information on offspring sex. By comparing sex-specific maps of the European tree frog Hyla arborea, we find that the sex chromosome (linkage group 1) shows a threefold difference in marker number between the male and female maps. In contrast, the number of markers on each autosome is similar between the two maps. We also find strongly conserved synteny between H. arborea and Xenopus tropicalis across 200 million years of evolution, suggesting that the rate of chromosomal rearrangement in anurans is low. Finally, we show that recombination in males is greatly reduced at the centers of large chromosomes, consistent with previous cytogenetic findings. Our research shows the importance of high-density linkage maps for studies of recombination, chromosomal rearrangement and the genetic architecture of ecologically or economically important traits.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.

Chromosomal Diversity and Karyotype Evolution in South American Macaws (Psittaciformes, Psittacidae)

Most species of macaws, which represent the largest species of Neotropical Psittacidae, characterized by their long tails and exuberant colours, are endangered, mainly because of hunting, illegal trade and habitat destruction. Long tailed species seem to represent a monophyletic group within Psittacidae, supported by cytogenetic data. Hence, these species show karyotypes with predominance of biarmed macrochromosomes, in contrast to short tailed species, with a predominance of acro/telocentric macrochromosomes. Because of their similar karyotypes, it has been proposed that inversions and translocations may be the main types of rearrangements occurring during the evolution of this group. However, only one species of macaw, Ara macao, that has had its genome sequenced was analyzed by means of molecular cytogenetics. Hence, in order to verify the rearrangements, we analyzed the karyotype of two species of macaws, Ara chloropterus and Anodorhynchus hyacinthinus, using cross-species chromosome painting with two different sets of probes from chicken and white hawk. Both intra- and interchromosomal rearrangements were observed. Chicken probes revealed the occurrence of fusions, fissions and inversions in both species, while the probes from white hawk determined the correct breakpoints or chromosome segments involved in the rearrangements. Some of these rearrangements were common for both species of macaws (fission of GGA1 and fusions of GGA1p/GGA4q, GGA6/GGA7 and GGA8/GGA9), while the fissions of GGA 2 and 4p were found only in A. chloropterus. These results confirm that despite apparent chromosomal similarity, macaws have very diverse karyotypes, which differ from each other not only by inversions and translocations as postulated before, but also by fissions and fusions.

Chromosome Painting in Vanellus chilensis: Detection of a Fusion Common to Clade Charadrii (Charadriiformes)

The Southern lapwing (Vanellus chilensis) is endemic to America and is well-known because of the vast expansion of its geographical distribution and its involvement in air accidents. Despite its popularity, there is no information concerning the genomic organization and karyotype of this species. Hence, because other species of the genus Vanellus have variable diploid numbers from 2n = 58 to 76, the aim of this report was to analyze the karyotype of V. chilensis by means of classical and molecular cytogenetics. We found that 2n = 78 and chromosome painting using probes of Gallus gallus (GGA) and Leucopternis albicollis revealed an organization similar to the avian putative ancestral karyotype, except for the fusion of GGA7 and GGA8, also found in Burhinus oedicnemus, the only Charadriiforme species analyzed by FISH so far. This rearrangement may represent a cytogenetic signature for this group and, in addition, must be responsible for the difference between the diploid number found in the avian putative ancestral karyotype (2n = 80) and V. chilensis (2n = 78).

Highly variable recombinational landscape modulates efficacy of natural selection in birds

Determining the rate of protein evolution and identifying the causes of its variation across the genome are powerful ways to understand forces that are important for genome evolution. By using a multitissue transcriptome data set from great tit (Parus major), we analyzed patterns of molecular evolution between two passerine birds, great tit and zebra finch (Taeniopygia guttata), using the chicken genome (Gallus gallus) as an outgroup. We investigated whether a special feature of avian genomes, the highly variable recombinational landscape, modulates the efficacy of natural selection through the effects of Hill-Robertson interference, which predicts that selection should be more effective in removing deleterious mutations and incorporating beneficial mutations in high-recombination regions than in low-recombination regions. In agreement with these predictions, genes located in low-recombination regions tend to have a high proportion of neutrally evolving sites and relaxed selective constraint on sites subject to purifying selection, whereas genes that show strong support for past episodes of positive selection appear disproportionally in high-recombination regions. There is also evidence that genes located in high-recombination regions tend to have higher gene expression specificity than those located in low-recombination regions. Furthermore, more compact genes (i.e., those with fewer/shorter introns or shorter proteins) evolve faster than less compact ones. In sum, our results demonstrate that transcriptome sequencing is a powerful method to answer fundamental questions about genome evolution in nonmodel organisms.

Gene turnover in the avian globin gene families and evolutionary changes in hemoglobin isoform expression

The apparent stasis in the evolution of avian chromosomes suggests that birds may have experienced relatively low rates of gene gain and loss in multigene families. To investigate this possibility and to explore the phenotypic consequences of variation in gene copy number, we examined evolutionary changes in the families of genes that encode the α- and β-type subunits of hemoglobin (Hb), the tetrameric α2β2 protein responsible for blood-O2 transport. A comparative genomic analysis of 52 bird species revealed that the size and membership composition of the α- and β-globin gene families have remained remarkably constant during approximately 100 My of avian evolution. Most interspecific variation in gene content is attributable to multiple independent inactivations of the α(D)-globin gene, which encodes the α-chain subunit of a functionally distinct Hb isoform (HbD) that is expressed in both embryonic and definitive erythrocytes. Due to consistent differences in O2-binding properties between HbD and the major adult-expressed Hb isoform, HbA (which incorporates products of the α(A)-globin gene), recurrent losses of α(D)-globin contribute to among-species variation in blood-O2 affinity. Analysis of HbA/HbD expression levels in the red blood cells of 122 bird species revealed high variability among lineages and strong phylogenetic signal. In comparison with the homologous gene clusters in mammals, the low retention rate for lineage-specific gene duplicates in the avian globin gene clusters suggests that the developmental regulation of Hb synthesis in birds may be more highly conserved, with orthologous genes having similar stage-specific expression profiles and similar functional properties in disparate taxa.

Conservation and losses of non-coding RNAs in avian genomes

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous “losses” of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

Low rate of interchromosomal rearrangements during old radiation of gekkotan lizards (Squamata: Gekkota)

Gekkotan lizards are a highly specious (∼1600 described species) clade of squamate lizards with nearly cosmopolitan distribution in warmer areas. The clade is primarily nocturnal and forms an ecologically dominant part of the world nocturnal herpetofauna. However, molecular cytogenetic methods to study the evolution of karyotypes have not been widely applied in geckos. Our aim here was to uncover the extent of chromosomal rearrangements across the whole group Gekkota and to search for putative synapomorphies supporting the newly proposed phylogenetic relationships within this clade. We applied cross-species chromosome painting with the recently derived whole-chromosomal probes from the gekkonid species Gekko japonicus to members of the major gekkotan lineages. We included members of the families Diplodactylidae, Carphodactylidae, Pygopodidae, Eublepharidae, Phyllodactylidae and Gekkonidae. Our study demonstrates relatively high chromosome conservatism across the ancient group of gekkotan lizards. We documented that many changes in chromosomal shape across geckos can be attributed to intrachromosomal rearrangements. The documented rearrangements are not totally in agreement with the recently newly erected family Phyllodactylidae. The results also pointed to homoplasy, particularly in the reuse of chromosome breakpoints, in the evolution of gecko karyotypes.

Genomic resources for the endangered Hawaiian honeycreepers

Background

The Hawaiian honeycreepers are an avian adaptive radiation containing many endangered and extinct species. They display a dramatic range of phenotypic variation and are a model system for studies of evolution, conservation, disease dynamics and population genetics. Development of a genome-scale resources for this group would augment the quality of research focusing on Hawaiian honeycreepers and facilitate comparative avian genomic research.

Results

We assembled the genome sequence of a Hawaii amakihi (Hemignathus virens),and identified ~3.9 million single nucleotide polymorphisms (SNPs) in the genome. Using the amakihi genome as a reference, we also identified ~156,000 SNPs in RAD tag (restriction site associated DNA) sequencing of five honeycreeper species (palila [Loxioides bailleui], Nihoa finch [Telespiza ultima], iiwi [Vestiaria coccinea], apapane [Himatione sanguinea], and amakihi). SNPs are distributed throughout the amakihi genome, and the individual sequenced shows several large regions of low heterozygosity on chromosomes 1, 5, 6, 8 and 11. SNPs from RAD tag sequencing were also found throughout the genome but were found to be more densely located on microchromosomes, apparently a result of differential distribution of the particular site recognized by restriction enzyme BseXI.

Conclusions

The amakihi genome sequence will be useful for comparative avian genomics research and provides a significant resource for studies in such areas as disease ecology, evolution, and conservation genetics. The genome sequences will enable mapping of transcriptome data for honeycreepers and comparison of gene sequences between avian taxa. Researchers will be able to use the large number of SNP markers to genotype honeycreepers in regions of interest or across the whole genome. There are enough markers to enable use of methods such as genome-wide association studies (GWAS) that will allow researchers to make connections between phenotypic diversity of honeycreepers and specific genetic variants. Genome-wide markers will also help resolve phylogenetic and population genetic questions in honeycreepers.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1098) contains supplementary material, which is available to authorized users.

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Reconstruction of gross avian genome structure, organization and evolution suggests that the chicken lineage most closely resembles the dinosaur avian ancestor

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Global patterns of apparent copy number variation in birds revealed by cross-species comparative genomic hybridization

There is a growing interest in copy number variation (CNV) and the recognition of its importance in phenotype, disease, adaptation and speciation. CNV data is usually ascertained by array-CGH within-species, but similar inter-species comparisons have also been made in primates, mice and domestic mammals. Here, we conducted a broad appraisal of putative cross-species CNVs in birds, 16 species in all, using the standard array-CGH approach. Using a chicken oligonucleotide microarray, we detected 790 apparent CNVs within 135 unique regions and developed a bioinformatic tool 'CNV Analyser' for analysing and visualising cross-species data sets. We successfully addressed four hypotheses as follows: (a) Cross-species CNVs (compared to chicken) are, as suggested from preliminary evidence, smaller and fewer in number than in mammals; this 'dogma' was rejected in the light of the new evidence. (b) CNVs in birds are likely to have a functional effect through an association with genes; a large proportion of detected regions (70 %) were indeed associated with genes (suggesting functional significance), however, not necessarily more so than in mammals. (c) There are more CNVs in birds with more rearranged karyotypes; this hypothesis was rejected. Indeed, Falco species contained fewer than most with relatively standard (chicken-like) karyotypes. (d) There are more CNVs per megabase on micro-chromosomes than macrochromosomes; this hypothesis was accepted. Indeed, in species with rearranged karyotypes characterised by chromosomal fusions, the fused former microchromosomes still 'behaved' as though they were their microchromosomal ancestors. Gene ontology analysis of CNVRs revealed enrichment in immune response and antigen presentation genes and five CNVRs were perfectly correlated with the unique loss of sexual dichromatism in one Galliformes species.

Comparative cytogenomics of poultry: mapping of single gene and repeat loci in the Japanese quail (Coturnix japonica)

Well-characterized molecular and cytogenetic maps are yet to be established in Japanese quail (Coturnix japonica). The aim of the current study was to cytogenetically map and determine linkage of specific genes and gene complexes in Japanese quail through the use of chicken (Gallus gallus) and turkey (Meleagris gallopavo) genomic DNA probes and conduct a comparative study among the three genomes. Chicken and turkey clones were used as probes on mitotic metaphase and meiotic pachytene stage chromosomes of the three species for the purpose of high-resolution fluorescence in situ hybridization (FISH). The genes and complexes studied included telomerase RNA (TR), telomerase reverse transcriptase (TERT), 5S rDNA, 18S-5.8S-28S rDNA (i.e., nucleolus organizer region (NOR)), and the major histocompatibility complex (MHC). The telomeric profile of Japanese quail was investigated through the use of FISH with a TTAGGG-PNA probe. A range of telomeric array sizes were confirmed as found for the other poultry species. Three NOR loci were identified in Japanese quail, and single loci each for TR, TERT, 5S rDNA and the MHC-B. The MHC-B and one NOR locus were linked on a microchromosome in Japanese quail. We confirmed physical linkage of 5S rDNA and the TR gene on an intermediate-sized chromosome in quail, similar to both chicken and turkey. TERT localized to CJA 2 in quail and the orthologous chromosome region in chicken (GGA 2) and in turkey (MGA 3). The cytogenetic profile of Japanese quail was further developed by this study and synteny was identified among the three poultry species.

Novel tools for characterising inter and intra chromosomal rearrangements in avian microchromosomes

Avian genome organisation is characterised, in part, by a set of microchromosomes that are unusually small in size and unusually large in number. Although containing about a quarter of the genome, they contain around half the genes and three quarters of the total chromosome number. Nonetheless, they continue to belie analysis by cytogenetic means. Chromosomal rearrangements play a key role in genome evolution, fertility and genetic disease and thus tools for analysis of the microchromosomes are essential to analyse such phenomena in birds. Here, we report the development of chicken microchromosomal paint pools, generation of pairs of specific microchromosome BAC clones in chicken, and computational tools for in silico comparison of the genomes of microchromosomes. We demonstrate the use of these molecular and computational tools across species, suggesting their use to generate a clear picture of microchromosomal rearrangements between avian species. With increasing numbers of avian genome sequences that are emerging, tools such as these will find great utility in assembling genomes de novo and for asking fundamental questions about genome evolution from a chromosomal perspective.