The genomic landscape of species divergence in Ficedula flycatchers

Advancing avian behavioral neuroendocrinology through genomics

Genome technologies are transforming all areas of biology, including the study of hormones, brain and behavior. Annotated reference genome assemblies are rapidly being produced for many avian species. Here we briefly review the basic concepts and tools used in genomics. We then consider how these are informing the study of avian behavioral neuroendocrinology, focusing in particular on lessons from the study of songbirds. We discuss the impact of having a complete "parts list" for an organism; the transformational potential of studying large sets of genes at once instead one gene at a time; the growing recognition that environmental and behavioral signals trigger massive shifts in gene expression in the brain; and the prospects for using comparative genomics to uncover the genetic roots of behavioral variation. Throughout, we identify promising new directions for bolstering the application of genomic information to further advance the study of avian brain and behavior.

Experimental evidence for ecological selection on genome variation in the wild

Understanding natural selection's effect on genetic variation is a major goal in biology, but the genome-scale consequences of contemporary selection are not well known. In a release and recapture field experiment we transplanted stick insects to native and novel host plants and directly measured allele frequency changes within a generation at 186 576 genetic loci. We observed substantial, genome-wide allele frequency changes during the experiment, most of which could be attributed to random mortality (genetic drift). However, we also documented that selection affected multiple genetic loci distributed across the genome, particularly in transplants to the novel host. Host-associated selection affecting the genome acted on both a known colour-pattern trait as well as other (unmeasured) phenotypes. We also found evidence that selection associated with elevation affected genome variation, although our experiment was not designed to test this. Our results illustrate how genomic data can identify previously underappreciated ecological sources and phenotypic targets of selection.

Transcriptome sequencing reveals the character of incomplete dosage compensation across multiple tissues in flycatchers

Sex chromosome divergence, which follows the cessation of recombination and degeneration of the sex-limited chromosome, can cause a reduction in expression level for sex-linked genes in the heterozygous sex, unless some mechanisms of dosage compensation develops to counter the reduction in gene dose. Because large-scale perturbations in expression levels arising from changes in gene dose might have strong deleterious effects, the evolutionary response should be strong. However, in birds and in at least some other female heterogametic organisms, wholesale sex chromosome dosage compensation does not seem to occur. Using RNA-seq of multiple tissues and individuals, we investigated male and female expression levels of Z-linked and autosomal genes in the collared flycatcher, a bird for which a draft genome sequence recently has been reported. We found that male expression of Z-linked genes was on average 50% higher than female expression, although there was considerable variation in the male-to-female ratio among genes. The ratio for individual genes was well correlated among tissues and there was also a correlation in the extent of compensation between flycatcher and chicken orthologs. The relative excess of male expression was positively correlated with expression breadth, expression level, and number of interacting proteins (protein connectivity), and negatively correlated with variance in expression. These observations lead to a model of compensation occurring on a gene-by-gene basis, supported by an absence of clustering of genes on the Z chromosome with respect to the extent of compensation. Equal mean expression level of autosomal and Z-linked genes in males, and 50% higher expression of autosomal than Z-linked genes in females, is compatible with that partial compensation is achieved by hypertranscription from females' single Z chromosome. A comparison with male-to-female expression ratios in orthologous Z-linked genes of ostriches, where Z-W recombination still occurs, suggests that male-biased expression of Z-linked genes is a derived trait after avian sex chromosome divergence.

Spatiotemporal dynamics of gene flow and hybrid fitness between the M and S forms of the malaria mosquito, Anopheles gambiae

The M and S forms of Anopheles gambiae have been the focus of intense study by malaria researchers and evolutionary biologists interested in ecological speciation. Divergence occurs at three discrete islands in genomes that are otherwise nearly identical. An "islands of speciation" model proposes that diverged regions contain genes that are maintained by selection in the face of gene flow. An alternative "incidental island" model maintains that gene flow between M and S is effectively zero and that divergence islands are unrelated to speciation. A "divergence island SNP" assay was used to explore the spatial and temporal distributions of hybrid genotypes. Results revealed that hybrid individuals occur at frequencies ranging between 5% and 97% in every population examined. A temporal analysis revealed that assortative mating is unstable and periodically breaks down, resulting in extensive hybridization. Results suggest that hybrids suffer a fitness disadvantage, but at least some hybrid genotypes are viable. Stable introgression of the 2L speciation island occurred at one site following a hybridization event.

Demographic divergence history of pied flycatcher and collared flycatcher inferred from whole-genome re-sequencing data

Profound knowledge of demographic history is a prerequisite for the understanding and inference of processes involved in the evolution of population differentiation and speciation. Together with new coalescent-based methods, the recent availability of genome-wide data enables investigation of differentiation and divergence processes at unprecedented depth. We combined two powerful approaches, full Approximate Bayesian Computation analysis (ABC) and pairwise sequentially Markovian coalescent modeling (PSMC), to reconstruct the demographic history of the split between two avian speciation model species, the pied flycatcher and collared flycatcher. Using whole-genome re-sequencing data from 20 individuals, we investigated 15 demographic models including different levels and patterns of gene flow, and changes in effective population size over time. ABC provided high support for recent (mode 0.3 my, range <0.7 my) species divergence, declines in effective population size of both species since their initial divergence, and unidirectional recent gene flow from pied flycatcher into collared flycatcher. The estimated divergence time and population size changes, supported by PSMC results, suggest that the ancestral species persisted through one of the glacial periods of middle Pleistocene and then split into two large populations that first increased in size before going through severe bottlenecks and expanding into their current ranges. Secondary contact appears to have been established after the last glacial maximum. The severity of the bottlenecks at the last glacial maximum is indicated by the discrepancy between current effective population sizes (20,000-80,000) and census sizes (5-50 million birds) of the two species. The recent divergence time challenges the supposition that avian speciation is a relatively slow process with extended times for intrinsic postzygotic reproductive barriers to evolve. Our study emphasizes the importance of using genome-wide data to unravel tangled demographic histories. Moreover, it constitutes one of the first examples of the inference of divergence history from genome-wide data in non-model species.

Validation of SNP allele frequencies determined by pooled next-generation sequencing in natural populations of a non-model plant species

Sequencing of pooled samples (Pool-Seq) using next-generation sequencing technologies has become increasingly popular, because it represents a rapid and cost-effective method to determine allele frequencies for single nucleotide polymorphisms (SNPs) in population pools. Validation of allele frequencies determined by Pool-Seq has been attempted using an individual genotyping approach, but these studies tend to use samples from existing model organism databases or DNA stores, and do not validate a realistic setup for sampling natural populations. Here we used pyrosequencing to validate allele frequencies determined by Pool-Seq in three natural populations of Arabidopsis halleri (Brassicaceae). The allele frequency estimates of the pooled population samples (consisting of 20 individual plant DNA samples) were determined after mapping Illumina reads to (i) the publicly available, high-quality reference genome of a closely related species (Arabidopsis thaliana) and (ii) our own de novo draft genome assembly of A. halleri. We then pyrosequenced nine selected SNPs using the same individuals from each population, resulting in a total of 540 samples. Our results show a highly significant and accurate relationship between pooled and individually determined allele frequencies, irrespective of the reference genome used. Allele frequencies differed on average by less than 4%. There was no tendency that either the Pool-Seq or the individual-based approach resulted in higher or lower estimates of allele frequencies. Moreover, the rather high coverage in the mapping to the two reference genomes, ranging from 55 to 284x, had no significant effect on the accuracy of the Pool-Seq. A resampling analysis showed that only very low coverage values (below 10-20x) would substantially reduce the precision of the method. We therefore conclude that a pooled re-sequencing approach is well suited for analyses of genetic variation in natural populations.

Differentiation in neutral genes and a candidate gene in the pied flycatcher: using biological archives to track global climate change

Global climate change is one of the major driving forces for adaptive shifts in migration and breeding phenology and possibly impacts demographic changes if a species fails to adapt sufficiently. In Western Europe, pied flycatchers (Ficedula hypoleuca) have insufficiently adapted their breeding phenology to the ongoing advance of food peaks within their breeding area and consequently suffered local population declines. We address the question whether this population decline led to a loss of genetic variation, using two neutral marker sets (mitochondrial control region and microsatellites), and one potentially selectively non-neutral marker (avian Clock gene). We report temporal changes in genetic diversity in extant populations and biological archives over more than a century, using samples from sites differing in the extent of climate change. Comparing genetic differentiation over this period revealed that only the recent Dutch population, which underwent population declines, showed slightly lower genetic variation than the historic Dutch population. As that loss of variation was only moderate and not observed in all markers, current gene flow across Western and Central European populations might have compensated local loss of variation over the last decades. A comparison of genetic differentiation in neutral loci versus the Clock gene locus provided evidence for stabilizing selection. Furthermore, in all genetic markers, we found a greater genetic differentiation in space than in time. This pattern suggests that local adaptation or historic processes might have a stronger effect on the population structure and genetic variation in the pied flycatcher than recent global climate changes.

Hybridization reveals the evolving genomic architecture of speciation

The rate at which genomes diverge during speciation is unknown, as are the physical dynamics of the process. Here, we compare full genome sequences of 32 butterflies, representing five species from a hybridizing Heliconius butterfly community, to examine genome-wide patterns of introgression and infer how divergence evolves during the speciation process. Our analyses reveal that initial divergence is restricted to a small fraction of the genome, largely clustered around known wing-patterning genes. Over time, divergence evolves rapidly, due primarily to the origin of new divergent regions. Furthermore, divergent genomic regions display signatures of both selection and adaptive introgression, demonstrating the link between microevolutionary processes acting within species and the origin of species across macroevolutionary timescales. Our results provide a uniquely comprehensive portrait of the evolving species boundary due to the role that hybridization plays in reducing the background accumulation of divergence at neutral sites.

The evolution of novelty in conserved genes; evidence of positive selection in the Drosophila fruitless gene is localised to alternatively spliced exons

There has been much debate concerning whether cis-regulatory or coding changes are more likely to produce evolutionary innovation or adaptation in gene function, but an additional complication is that some genes can dramatically diverge through alternative splicing, increasing the diversity of gene function within a locus. The fruitless gene is a major transcription factor with a wide range of pleiotropic functions, including a fundamental conserved role in sexual differentiation, species-specific morphology and an important influence on male sexual behaviour. Here, we examine the structure of fruitless in multiple species of Drosophila, and determine the patterns of selective constraint acting across the coding region. We found that the pattern of selection, estimated from the ratio of non-synonymous to synonymous substitutions, varied considerably across the gene, with most regions of the gene evolutionarily conserved but with several regions showing evidence of divergence as a result of positive selection. The regions that showed evidence of positive selection were found to be localised to relatively consistent regions across multiple speciation events, and are associated with alternative splicing. Alternative splicing may thus provide a route to gene diversification in key regulatory loci.

Why time matters: codon evolution and the temporal dynamics of dN/dS

The ratio of divergence at nonsynonymous and synonymous sites, dN/dS, is a widely used measure in evolutionary genetic studies to investigate the extent to which selection modulates gene sequence evolution. Originally tailored to codon sequences of distantly related lineages, dN/dS represents the ratio of fixed nonsynonymous to synonymous differences. The impact of ancestral and lineage-specific polymorphisms on dN/dS, which we here show to be substantial for closely related lineages, is generally neglected in estimation techniques of dN/dS. To address this issue, we formulate a codon model that is firmly anchored in population genetic theory, derive analytical expressions for the dN/dS measure by Poisson random field approximation in a Markovian framework and validate the derivations by simulations. In good agreement, simulations and analytical derivations demonstrate that dN/dS is biased by polymorphisms at short time scales and that it can take substantial time for the expected value to settle at its time limit where only fixed differences are considered. We further show that in any attempt to estimate the dN/dS ratio from empirical data the effect of the intrinsic fluctuations of a ratio of stochastic variables, can even under neutrality yield extreme values of dN/dS at short time scales or in regions of low mutation rate. Taken together, our results have significant implications for the interpretation of dN/dS estimates, the McDonald-Kreitman test and other related statistics, in particular for closely related lineages.

Genome-wide evidence for speciation with gene flow in Heliconius butterflies

Most speciation events probably occur gradually, without complete and immediate reproductive isolation, but the full extent of gene flow between diverging species has rarely been characterized on a genome-wide scale. Documenting the extent and timing of admixture between diverging species can clarify the role of geographic isolation in speciation. Here we use new methodology to quantify admixture at different stages of divergence in Heliconius butterflies, based on whole-genome sequences of 31 individuals. Comparisons between sympatric and allopatric populations of H. melpomene, H. cydno, and H. timareta revealed a genome-wide trend of increased shared variation in sympatry, indicative of pervasive interspecific gene flow. Up to 40% of 100-kb genomic windows clustered by geography rather than by species, demonstrating that a very substantial fraction of the genome has been shared between sympatric species. Analyses of genetic variation shared over different time intervals suggested that admixture between these species has continued since early in speciation. Alleles shared between species during recent time intervals displayed higher levels of linkage disequilibrium than those shared over longer time intervals, suggesting that this admixture took place at multiple points during divergence and is probably ongoing. The signal of admixture was significantly reduced around loci controlling divergent wing patterns, as well as throughout the Z chromosome, consistent with strong selection for Müllerian mimicry and with known Z-linked hybrid incompatibility. Overall these results show that species divergence can occur in the face of persistent and genome-wide admixture over long periods of time.

POPBAM: Tools for Evolutionary Analysis of Short Read Sequence Alignments

BACKGROUND

While many bioinformatics tools currently exist for assembling and discovering variants from next-generation sequence data, there are very few tools available for performing evolutionary analyses from these data. Evolutionary and population genomics studies hold great promise for providing valuable insights into natural selection, the effect of mutations on phenotypes, and the origin of species. Thus, there is a need for an extensible and flexible computational tool that can function into a growing number of evolutionary bioinformatics pipelines.

RESULTS

This paper describes the POPBAM software, which is a comprehensive set of computational tools for evolutionary analysis of whole-genome alignments consisting of multiple individuals, from multiple populations or species. POPBAM works directly from BAM-formatted assembly files, calls variant sites, and calculates a variety of commonly used evolutionary sequence statistics. POPBAM is designed primarily to perform analyses in sliding windows across chromosomes or scaffolds. POPBAM accurately measures nucleotide diversity, population divergence, linkage disequilibrium, and the frequency spectrum of mutations from two or more populations. POPBAM can also produce phylogenetic trees of all samples in a BAM file. Finally, I demonstrate that the implementation of POPBAM is both fast and memory-efficient, and also can feasibly scale to the analysis of large BAM files with many individuals and populations. Software: The POPBAM program is written in C/C++ and is available from http://dgarriga.github.io/POPBAM. The program has few dependencies and can be built on a variety of Linux platforms. The program is open-source and users are encouraged to participate in the development of this resource.

The faster-X effect: integrating theory and data

Population genetics theory predicts that X (or Z) chromosomes could play disproportionate roles in speciation and evolutionary divergence, and recent genome-wide analyses have identified situations in which X or Z-linked divergence exceeds that on the autosomes (the so-called 'faster-X effect'). Here, we summarize the current state of both the theory and data surrounding the study of faster-X evolution. Our survey indicates that the faster-X effect is pervasive across a taxonomically diverse array of evolutionary lineages. These patterns could be informative of the dominance or recessivity of beneficial mutations and the nature of genetic variation acted upon by natural selection. We also identify several aspects of disagreement between these empirical results and the population genetic models used to interpret them. However, there are clearly delineated aspects of the problem for which additional modeling and collection of genomic data will address these discrepancies and provide novel insights into the population genetics of adaptation.

Segmenting the human genome based on states of neutral genetic divergence

Many studies have demonstrated that divergence levels generated by different mutation types vary and covary across the human genome. To improve our still-incomplete understanding of the mechanistic basis of this phenomenon, we analyze several mutation types simultaneously, anchoring their variation to specific regions of the genome. Using hidden Markov models on insertion, deletion, nucleotide substitution, and microsatellite divergence estimates inferred from human-orangutan alignments of neutrally evolving genomic sequences, we segment the human genome into regions corresponding to different divergence states--each uniquely characterized by specific combinations of divergence levels. We then parsed the mutagenic contributions of various biochemical processes associating divergence states with a broad range of genomic landscape features. We find that high divergence states inhabit guanine- and cytosine (GC)-rich, highly recombining subtelomeric regions; low divergence states cover inner parts of autosomes; chromosome X forms its own state with lowest divergence; and a state of elevated microsatellite mutability is interspersed across the genome. These general trends are mirrored in human diversity data from the 1000 Genomes Project, and departures from them highlight the evolutionary history of primate chromosomes. We also find that genes and noncoding functional marks [annotations from the Encyclopedia of DNA Elements (ENCODE)] are concentrated in high divergence states. Our results provide a powerful tool for biomedical data analysis: segmentations can be used to screen personal genome variants--including those associated with cancer and other diseases--and to improve computational predictions of noncoding functional elements.

Inferring the geographic mode of speciation by contrasting autosomal and sex-linked genetic diversity

When geographic isolation drives speciation, concurrent termination of gene flow among genomic regions will occur immediately after the formation of the barrier between diverging populations. Alternatively, if speciation is driven by ecologically divergent selection, gene flow of selectively neutral genomic regions may go on between diverging populations until the completion of reproductive isolation. It may also lead to an unsynchronized termination of gene flow between genomic regions with different roles in the speciation process. Here, we developed a novel Approximate Bayesian Computation pipeline to infer the geographic mode of speciation by testing for a lack of postdivergence gene flow and a concurrent termination of gene flow in autosomal and sex-linked markers jointly. We applied this approach to infer the geographic mode of speciation for two allopatric highland rosefinches, the vinaceous rosefinch Carpodacus vinaceus and the Taiwan rosefinch C. formosanus from DNA polymorphisms of both autosomal and Z-linked loci. Our results suggest that the two rosefinch species diverged allopatrically approximately 0.5 Ma. Our approach allowed us further to infer that female effective population sizes are about five times larger than those of males, an estimate potentially useful when comparing the intensity of sexual selection across species.

Lessons learned from implementing a national infrastructure in Sweden for storage and analysis of next-generation sequencing data

: Analyzing and storing data and results from next-generation sequencing (NGS) experiments is a challenging task, hampered by ever-increasing data volumes and frequent updates of analysis methods and tools. Storage and computation have grown beyond the capacity of personal computers and there is a need for suitable e-infrastructures for processing. Here we describe UPPNEX, an implementation of such an infrastructure, tailored to the needs of data storage and analysis of NGS data in Sweden serving various labs and multiple instruments from the major sequencing technology platforms. UPPNEX comprises resources for high-performance computing, large-scale and high-availability storage, an extensive bioinformatics software suite, up-to-date reference genomes and annotations, a support function with system and application experts as well as a web portal and support ticket system. UPPNEX applications are numerous and diverse, and include whole genome-, de novo- and exome sequencing, targeted resequencing, SNP discovery, RNASeq, and methylation analysis. There are over 300 projects that utilize UPPNEX and include large undertakings such as the sequencing of the flycatcher and Norwegian spruce. We describe the strategic decisions made when investing in hardware, setting up maintenance and support, allocating resources, and illustrate major challenges such as managing data growth. We conclude with summarizing our experiences and observations with UPPNEX to date, providing insights into the successful and less successful decisions made.

Speciation: clash of the genomes

Complete genomes of hybridizing bird species demonstrate the importance of the sex chromosomes, telomeres and centromeres to the initial stages of speciation.

Specific-gene studies of evolutionary mechanisms in an age of genome-wide surveying

The molecular tools of genomics have great power to reveal patterns of genetic difference within or among species, but must be complemented by the mechanistic study of the genetic variants found if these variants' evolutionary meaning is to be well understood. Central to this purpose is knowledge of the organisms' genotype-phenotype-environment interactions, which embody biological adaptation and constraint and thus drive natural selection. The history of this approach is briefly reviewed. Strategies embracing the complementarity of genomics and specific-gene studies in evolution are considered. Implementation of these strategies, and examples showing their feasibility and power, are discussed. Initial generalizations emphasize: (1) reproducibility of adaptive mechanisms; (2) evolutionary co-importance of variation in protein sequences and expression; (3) refinement of rudimentary molecular functions as an origin of evolutionary innovations; (4) identification of specific-gene mechanisms as underpinnings of genomic or quantitative genetic variation; and (5) multiple forms of adaptive or constraining epistasis among genes. Progress along these lines will advance understanding of evolution and support its use in addressing urgent medical and environmental applications.

Characterisation of a transcriptome to find sequence differences between two differentially migrating subspecies of the willow warbler Phylloscopus trochilus

BACKGROUND

Animal migration requires adaptations in morphological, physiological and behavioural traits. Several of these traits have been shown to possess a strong heritable component in birds, but little is known about their genetic architecture. Here we used 454 sequencing of brain-derived transcriptomes from two differentially migrating subspecies of the willow warbler Phylloscopus trochilus to detect genes potentially underlying traits associated with migration.

RESULTS

The transcriptome sequencing resulted in 1.8 million reads following filtering steps. Most of the reads (84%) were successfully mapped to the genome of the zebra finch Taeniopygia gutatta. The mapped reads were situated within at least 12,101 predicted zebra finch genes, with the greatest sequencing depth in exons. Reads that were mapped to intergenic regions were generally located close to predicted genes and possibly located in uncharacterized untranslated regions (UTRs). Out of 85,000 single nucleotide polymorphisms (SNPs) with a minimum sequencing depth of eight reads from each of two subspecies-specific pools, only 55 showed high differentiation, confirming previous studies showing that most of the genetic variation is shared between the subspecies. Validation of a subset of the most highly differentiated SNPs using Sanger sequencing demonstrated that several of them also were differentiated between an independent set of individuals of each subspecies. These SNPs were clustered in two chromosome regions that are likely to be influenced by divergent selection between the subspecies and that could potentially be associated with adaptations to their different migratory strategies.

CONCLUSIONS

Our study represents the first large-scale sequencing analysis aiming at detecting genes underlying migratory phenotypes in birds and provides new candidates for genes potentially involved in migration.

Understanding the origin of species with genome-scale data: modelling gene flow

As it becomes easier to sequence multiple genomes from closely related species, evolutionary biologists working on speciation are struggling to get the most out of very large population genomic data sets. Such data hold the potential to resolve long-standing questions in evolutionary biology about the role of gene exchange in species formation. In principle, the new population genomic data can be used to disentangle the conflicting roles of natural selection and gene flow during the divergence process. However, there are great challenges in taking full advantage of such data, especially with regard to including recombination in genetic models of the divergence process. Current data, models, methods and the potential pitfalls in using them will be considered here.

A multi-platform draft de novo genome assembly and comparative analysis for the Scarlet Macaw (Ara macao)

Data deposition to NCBI Genomes: This Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession AMXX00000000 (SMACv1.0, unscaffolded genome assembly). The version described in this paper is the first version (AMXX01000000). The scaffolded assembly (SMACv1.1) has been deposited at DDBJ/EMBL/GenBank under the accession AOUJ00000000, and is also the first version (AOUJ01000000). Strong biological interest in traits such as the acquisition and utilization of speech, cognitive abilities, and longevity catalyzed the utilization of two next-generation sequencing platforms to provide the first-draft de novo genome assembly for the large, new world parrot Ara macao (Scarlet Macaw). Despite the challenges associated with genome assembly for an outbred avian species, including 951,507 high-quality putative single nucleotide polymorphisms, the final genome assembly (>1.035 Gb) includes more than 997 Mb of unambiguous sequence data (excluding N's). Cytogenetic analyses including ZooFISH revealed complex rearrangements associated with two scarlet macaw macrochromosomes (AMA6, AMA7), which supports the hypothesis that translocations, fusions, and intragenomic rearrangements are key factors associated with karyotype evolution among parrots. In silico annotation of the scarlet macaw genome provided robust evidence for 14,405 nuclear gene annotation models, their predicted transcripts and proteins, and a complete mitochondrial genome. Comparative analyses involving the scarlet macaw, chicken, and zebra finch genomes revealed high levels of nucleotide-based conservation as well as evidence for overall genome stability among the three highly divergent species. Application of a new whole-genome analysis of divergence involving all three species yielded prioritized candidate genes and noncoding regions for parrot traits of interest (i.e., speech, intelligence, longevity) which were independently supported by the results of previous human GWAS studies. We also observed evidence for genes and noncoding loci that displayed extreme conservation across the three avian lineages, thereby reflecting their likely biological and developmental importance among birds.

Low fertility of wild hybrid male flycatchers despite recent divergence

Postzygotic isolation may be important for maintaining species boundaries, particularly when premating barriers are incomplete. Little is known about the course of events leading from minor environmental mismatches affecting hybrid fitness to severe genetic incompatibilities causing sterility or inviability. We investigated whether reduced reproductive success of hybrid males was caused by suboptimal sperm traits or by more severe genetic incompatibilities in a hybrid zone of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) on the island of Öland, Sweden. About 4 per cent hybridization is observed in this population and all female hybrids are sterile. We found no sperm in the ejaculates of most sampled hybrid males, and sperm with abnormal morphology in two hybrids. Furthermore, none of the hybrids sired any offspring because of high levels of hatching failure and extra-pair paternity in their nests. These results from a natural hybrid zone suggest that the spermatogenesis of hybrid males may become disrupted despite little genetic divergence between the parental species.

Insights into the evolution of Darwin's finches from comparative analysis of the Geospiza magnirostris genome sequence

Background

A classical example of repeated speciation coupled with ecological diversification is the evolution of 14 closely related species of Darwin’s (Galápagos) finches (Thraupidae, Passeriformes). Their adaptive radiation in the Galápagos archipelago took place in the last 2–3 million years and some of the molecular mechanisms that led to their diversification are now being elucidated. Here we report evolutionary analyses of genome of the large ground finch, Geospiza magnirostris.

Results

13,291 protein-coding genes were predicted from a 991.0 Mb G. magnirostris genome assembly. We then defined gene orthology relationships and constructed whole genome alignments between the G. magnirostris and other vertebrate genomes. We estimate that 15% of genomic sequence is functionally constrained between G. magnirostris and zebra finch. Genic evolutionary rate comparisons indicate that similar selective pressures acted along the G. magnirostris and zebra finch lineages suggesting that historical effective population size values have been similar in both lineages. 21 otherwise highly conserved genes were identified that each show evidence for positive selection on amino acid changes in the Darwin's finch lineage. Two of these genes (Igf2r and Pou1f1) have been implicated in beak morphology changes in Darwin’s finches. Five of 47 genes showing evidence of positive selection in early passerine evolution have cilia related functions, and may be examples of adaptively evolving reproductive proteins.

Conclusions

These results provide insights into past evolutionary processes that have shaped G. magnirostris genes and its genome, and provide the necessary foundation upon which to build population genomics resources that will shed light on more contemporaneous adaptive and non-adaptive processes that have contributed to the evolution of the Darwin’s finches.

Relative performance of hybrid nestlings in Ficedula flycatchers: a translocation experiment

Ecological speciation predicts that hybrids should experience relatively low fitness in the local environments of their parental species. In this study, we performed a translocation experiment of nestling hybrids between collared and pied flycatchers into the nests of conspecific pairs of their parental species. Our aim was to compare the performance of hybrids with purebred nestlings. Nestling collared flycatchers are known to beg and grow faster than nestling pied flycatchers under favorable conditions, but to experience higher mortality than nestling pied flycatchers under food limitation. The experiment was performed relatively late in the breeding season when food is limited. If hybrid nestlings have an intermediate growth potential and begging intensity, we expected them to beg and grow faster, but also to experience lower survival than pied flycatchers. In comparison with nestling collared flycatchers, we expected them to beg and grow slower, but to survive better. We found that nestling collared flycatchers indeed begged significantly faster and experienced higher mortality than nestling hybrids. Moreover, nestling hybrids had higher weight and tended to beg faster than nestling pied flycatchers, but we did not detect a difference in survival between the latter two groups of nestlings. We conclude that hybrid Ficedula nestlings appear to have a better intrinsic adaptation to food limitation late in the breeding season compared with nestling collared flycatchers. We discuss possible implications for gene flow between the two species.

Genomic diversity and evolution of the head crest in the rock pigeon

The geographic origins of breeds and the genetic basis of variation within the widely distributed and phenotypically diverse domestic rock pigeon (Columba livia) remain largely unknown. We generated a rock pigeon reference genome and additional genome sequences representing domestic and feral populations. We found evidence for the origins of major breed groups in the Middle East and contributions from a racing breed to North American feral populations. We identified the gene EphB2 as a strong candidate for the derived head crest phenotype shared by numerous breeds, an important trait in mate selection in many avian species. We also found evidence that this trait evolved just once and spread throughout the species, and that the crest originates early in development by the localized molecular reversal of feather bud polarity.

Lack of dosage compensation accompanies the arrested stage of sex chromosome evolution in ostriches

Sex chromosome evolution is usually seen as a process that, once initiated, will inevitably progress toward an advanced stage of degeneration of the nonrecombining chromosome. However, despite evidence that avian sex chromosome evolution was initiated >100 Ma, ratite birds have been trapped in an arrested stage of sex chromosome divergence. We performed RNA sequencing of several tissues from male and female ostriches and assembled the transcriptome de novo. A total of 315 Z-linked genes fell into two categories: those that have equal expression level in the two sexes (for which Z–W recombination still occurs) and those that have a 2-fold excess of male expression (for which Z–W recombination has ceased). We suggest that failure to evolve dosage compensation has constrained sex chromosome divergence in this basal avian lineage. Our results indicate that dosage compensation is a prerequisite for, not only a consequence of, sex chromosome evolution.