The zebra finch is an important model organism in several fields1,2 with unique relevance to human neuroscience3,4. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken5—the only bird with a sequenced genome until now6. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes7. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation 1 . The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models 2-7 that show pre-as well as postzygotic isolation 8,9 . We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map 10 and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.As lineages diverge, a combination of pre-as well as postzygotic reproductive isolation barriers will eventually arise 1 . Divergence is likely to start from specific loci that may precede and cause the evolution of reproductive incompatibility. Hybridization between diverging lineages may therefore create a genomic mosaic of regions where interspecific gene flow occurs at different rates (the genic view of speciation 11 ), with introgression expected to be weak in genomic regions involved in speciation. Revealing the genomic regions with elevated levels of divergence will eventually deepen our knowledge of the speciation process. However, more than 150 years after the publication of On the Origin of Species 12 , the genetic basis of speciation is still largely unresolved 13,14 . We know little about the identity, number and effect size of loci involved in population divergence, their genomic distribution and the type of mutations involved. Advances in sequencing technology now open a promising avenue for the study of genomic divergence, even for non-model vertebrate species with gigabase (Gb)-sized genomes.The collared flycatcher Ficedula albicollis and the pied flycatcher Ficedula hypoleuca (Fig. 1) are important study organisms for key aspects of evolutionary ecology and biology 2-7 . Diverged less than 2 million years ago, their history has been shaped by repeated cycles of glaciation in Eurasia where periods of allopatric divergence in refugia probably alternated with periods of secondary contact...

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cm and a mean genetic distance between markers of 0.12 cm. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cm/Mb and correlated closely with chromosome size, from 2 cm/Mb for chromosomes >100 Mb to >10 cm/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cm was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Understanding the causes and consequences of variation in the rate of recombination is essential since this parameter is considered to affect levels of genetic variability, the efficacy of selection, and the design of association and linkage mapping studies. However, there is limited knowledge about the factors governing recombination rate variation. We genotyped 1920 single nucleotide polymorphisms in a multigeneration pedigree of more than 1000 zebra finches (Taeniopygia guttata) to develop a genetic linkage map, and then we used these map data together with the recently available draft genome sequence of the zebra finch to estimate recombination rates in 1 Mb intervals across the genome. The average zebra finch recombination rate (1.5 cM/Mb) is higher than in humans, but significantly lower than in chicken. The local rates of recombination in chicken and zebra finch were only weakly correlated, demonstrating evolutionary turnover of the recombination landscape in birds. The distribution of recombination events was heavily biased toward ends of chromosomes, with a stronger telomere effect than so far seen in any organism. In fact, the recombination rate was as low as 0.1 cM/Mb in intervals up to 100 Mb long in the middle of the larger chromosomes. We found a positive correlation between recombination rate and GC content, as well as GC-rich sequence motifs. Levels of linkage disequilibrium (LD) were significantly higher in regions of low recombination, showing that heterogeneity in recombination rates have left a footprint on the genomic landscape of LD in zebra finch populations.

Although there is growing interest to take genomics into the complex realms of natural populations, there is a general shortage of genomic resources and tools available for wild species. This applies not at least to birds, for which genomic approaches should be helpful to questions such as adaptation, speciation and population genetics. In this study, we describe a genome-wide reference set of conserved avian gene markers, broadly applicable across birds. By aligning protein-coding sequences from the recently assembled chicken genome with orthologous sequences in zebra finch, we identified particularly conserved exonic regions flanking introns of suitable size for subsequent amplification and sequencing. Primers were designed for 242 gene markers evenly distributed across the chicken genome, with a mean inter-marker interval of 4.2 Mb. Between 78% and 93% of the markers amplified a specific product in five species tested (chicken, peregrine falcon, collared flycatcher, great reed warbler and blue tit). Two hundred markers were sequenced in collared flycatcher, yielding a total of 122.41 kb of genomic DNA sequence (12096 bp coding sequence and 110 314 bp noncoding). Intron size of collared flycatcher and chicken was highly correlated, as was GC content. A polymorphism screening using these markers in a panel of 10 unrelated collared flycatchers identified 871 single nucleotide polymorphisms (pi = 0.0029) and 33 indels (mainly very short). Avian genome characteristics such as uniform genome size and low rate of syntenic rearrangements suggest that this marker set will find broad utility as a genome-wide reference resource for molecular ecological and population genomic analysis of birds. We envision that it will be particularly useful for obtaining large-scale orthologous targets in different species--important in, for instance, phylogenetics--and for large-scale identification of evenly distributed single nucleotide polymorphisms needed in linkage mapping or in studies of gene flow and hybridization.

The view of species as entities subjected to natural selection and amenable to change put forth by Charles Darwin and Alfred Wallace laid the conceptual foundation for understanding speciation. Initially marred by a rudimental understanding of hereditary principles, evolutionists gained appreciation of the mechanistic underpinnings of speciation following the merger of Mendelian genetic principles with Darwinian evolution. Only recently have we entered an era where deciphering the molecular basis of speciation is within reach. Much focus has been devoted to the genetic basis of intrinsic postzygotic isolation in model organisms and several hybrid incompatibility genes have been successfully identified. However, concomitant with the recent technological advancements in genome analysis and a newfound interest in the role of ecology in the differentiation process, speciation genetic research is becoming increasingly open to non-model organisms. This development will expand speciation research beyond the traditional boundaries and unveil the genetic basis of speciation from manifold perspectives and at various stages of the splitting process. This review aims at providing an extensive overview of speciation genetics. Starting from key historical developments and core concepts of speciation genetics, we focus much of our attention on evolving approaches and introduce promising methodological approaches for future research venues.

Chromosomal rearrangements and copy number variants (CNVs) play key roles in genome evolution and genetic disease; however, the molecular mechanisms underlying these types of structural genomic variation are not fully understood. The availability of complete genome sequences for two bird species, the chicken and the zebra finch, provides, for the first time, an ideal opportunity to analyze the relationship between structural genomic variation (chromosomal and CNV) and recombination on a genome-wide level. The aims of this study were therefore threefold: (1) to combine bioinformatics, physical mapping to produce comprehensive comparative maps of the genomes of chicken and zebra finch. In so doing, this allowed the identification of evolutionary chromosomal rearrangements distinguishing them. The previously reported interchromosomal conservation of synteny was confirmed, but a larger than expected number of intrachromosomal rearrangements were reported; (2) to hybridize zebra finch genomic DNA to a chicken tiling path microarray and identify CNVs in the zebra finch genome relative to chicken; 32 interspecific CNVs were identified; and (3) to test the hypothesis that there is an association between CNV, chromosomal rearrangements, and recombination by correlating data from (1) and (2) with recombination rate data from a high-resolution genetic linkage map of the zebra finch. We found a highly significant association of both chromosomal rearrangements and CNVs with elevated recombination rates. The results thus provide support for the notion of recombination-based processes playing a major role in avian genome evolution.

No abstract