Natural hybrid zones provide powerful opportunities for identifying the mechanisms that facilitate and inhibit speciation. Documenting the extent of genomic admixture allows us to discern the architecture of reproductive isolation through the identification of isolating barriers. This approach is particularly powerful for characterizing the accumulation of isolating barriers in systems exhibiting varying levels of genomic divergence. Here, we use a hybrid zone between two species-the Baltimore (Icterus galbula) and Bullock's (I. bullockii) orioles-to investigate this architecture of reproductive isolation. We combine whole genome re-sequencing with data from an additional 313 individuals amplityped at ancestry-informative markers to characterize fine-scale patterns of admixture, and to quantify links between genes and the plumage traits. On a genome-wide scale, we document several putative barriers to reproduction, including elevated peaks of divergence above a generally high genomic baseline, a large putative inversion on the Z chromosome, and complex interactions between melanogenesis-pathway candidate genes. Concordant and coincident clines for these different genomic regions further suggest the coupling of pre- and post-mating barriers. Our findings of complex and coupled interactions between pre- and post-mating barriers suggest a relatively rapid accumulation of barriers between these species, and they demonstrate the complexities of the speciation process.

Colour polymorphisms play a key role in sexual selection and speciation, yet the mechanisms that generate and maintain them are not fully understood. Here, we use genomic and transcriptomic tools to identify the precise genetic architecture and evolutionary history of a sex-linked colour polymorphism in the Gouldian finch Erythrura gouldiae that is also accompanied by remarkable differences in behaviour and physiology. We find that differences in colour are associated with an ~72-kbp region of the Z chromosome in a putative regulatory region for follistatin, an antagonist of the TGF-β superfamily genes. The region is highly differentiated between morphs, unlike the rest of the genome, yet we find no evidence that an inversion is involved in maintaining the distinct haplotypes. Coalescent simulations confirm that there is elevated nucleotide diversity and an excess of intermediate frequency alleles at this locus. We conclude that this pleiotropic colour polymorphism is most probably maintained by balancing selection.

Instances of recent and rapid speciation are suitable for associating phenotypes with their causal genotypes, especially if gene flow homogenizes areas of the genome that are not under divergent selection. We study a rapid radiation of nine sympatric bird species known as capuchino seedeaters, which are differentiated in sexually selected characters of male plumage and song. We sequenced the genomes of a phenotypically diverse set of species to search for differentiated genomic regions. Capuchinos show differences in a small proportion of their genomes, yet selection has acted independently on the same targets in different members of this radiation. Many divergent regions contain genes involved in the melanogenesis pathway, with the strongest signal originating from putative regulatory regions. Selection has acted on these same genomic regions in different lineages, likely shaping the evolution of cis-regulatory elements, which control how more conserved genes are expressed and thereby generate diversity in classically sexually selected traits.

Speciation is one of the most important processes in biology, yet the study of the genomic changes underlying this process is in its infancy. North American warbler species Setophaga townsendi and Setophaga occidentalis hybridize in a stable hybrid zone, following a period of geographic separation. Genomic differentiation accumulated during geographic isolation can be homogenized by introgression at secondary contact, whereas genetic regions that cause low hybrid fitness can be shielded from such introgression. Here, we examined the genomic underpinning of speciation by investigating (1) the genetic basis of divergent pigmentation traits between species, (2) variation in differentiation across the genome, and (3) the evidence for selection maintaining differentiation in the pigmentation genes. Using tens of thousands of single nucleotide polymorphisms (SNPs) genotyped in hundreds of individuals within and near the hybrid zone, genome-wide association mapping revealed a single SNP associated with cheek, crown, breast coloration, and flank streaking, reflecting pleiotropy (one gene affecting multiple traits) or close physical linkage of different genes affecting different traits. This SNP is within an intron of the RALY gene, hence we refer to it as the RALY SNP. We then examined between-species genomic differentiation, using both genotyping-by-sequencing and whole genome sequencing. We found that the RALY SNP is within one of the highest peaks of differentiation, which contains three genes known to influence pigmentation: ASIP, EIF2S2, and RALY (the ASIP-RALY gene block). Heterozygotes at this gene block are likely of reduced fitness, as the geographic cline of the RALY SNP has been narrow over two decades. Together, these results reflect at least one barrier to gene flow within this narrow (∼200 kb) genomic region that modulates plumage difference between species. Despite extensive gene flow between species across the genome, this study provides evidence that selection on a phenotype-associated genomic region maintains a stable species boundary.

Many vertebrates breed in cooperative groups in which more than two members provide care for young. Studies of cooperative breeding behavior within species have long highlighted the importance of environmental factors in mediating the paradox of why some such individuals delay independent breeding to help raise the offspring of others. In contrast, studies involving comparisons among species have not shown a similarly clear evolutionary-scale relationship between the interspecific incidence of cooperative breeding and any environmental factors. Here, we use a phylogenetically controlled comparative analysis of a complete, socially diverse group of birds-45 species of African starlings-to show that cooperative breeding is positively associated with living in semiarid savanna habitats and with temporal variability in rainfall. Savanna habitats are not only highly seasonal, but also temporally variable and unpredictable, and this temporal variability directly influences individual reproductive decisions in starlings and helps explain interspecific patterns of sociality. Cooperative breeding is likely to be adaptive in temporally variable environments because it allows for both reproduction in harsh years and sustained breeding during benign years. This "temporal variability" hypothesis might help explain the phylogenetic and geographic concentrations of cooperatively breeding vertebrates in savanna-like habitats and other temporally variable environments worldwide.

Animal migration demands an interconnected suite of adaptations for individuals to navigate over long distances. This trait complex is crucial for small birds whose migratory behaviors-such as directionality-are more likely innate, rather than being learned as in many longer-lived birds. Identifying causal genes has been a central goal of migration ecology, and this endeavor has been furthered by genome-scale comparisons. However, even the most successful studies of migration genetics have achieved low-resolution associations, identifying large chromosomal regions that encompass hundreds of genes, one or more of which might be causal. Here we leverage the genomic similarity among golden-winged (Vermivora chrysoptera) and blue-winged (V. cyanoptera) warblers to identify a single gene-vacuolar protein sorting 13A (VPS13A)-that is associated with distinct differences in migration to Central American (CA) or South American (SA) wintering areas. We find reduced sequence variation in this gene region for SA wintering birds, and show this is the likely result of natural selection on this locus. In humans, variants of VPS13A are linked to the neurodegenerative disorder chorea-acanthocytosis. This association provides one of the strongest gene-level associations with avian migration differences.

The Horned Lark (Eremophila alpestris) is a small songbird that exhibits remarkable geographic variation in appearance and habitat across an expansive distribution. While E. alpestris has been the focus of many ecological and evolutionary studies, we still lack a highly contiguous genome assembly for the Horned Lark and related taxa (Alaudidae). Here, we present CLO_EAlp_1.0, a highly contiguous assembly for E. alpestris generated from a blood sample of a wild, male bird captured in the Altiplano Cundiboyacense of Colombia. By combining short-insert and mate-pair libraries with the ALLPATHS-LG genome assembly pipeline, we generated a 1.04 Gb assembly comprised of 2713 scaffolds, with a largest scaffold size of 31.81 Mb, a scaffold N50 of 9.42 Mb, and a scaffold L50 of 30. These scaffolds were assembled from 23685 contigs, with a largest contig size of 1.69 Mb, a contig N50 of 193.81 kb, and a contig L50 of 1429. Our assembly pipeline also produced a single mitochondrial DNA contig of 14.00 kb. After polishing the genome, we identified 94.5% of single-copy gene orthologs from an Aves data set and 97.7% of single-copy gene orthologs from a vertebrata data set, which further demonstrates the high quality of our assembly. We anticipate that this genomic resource will be useful to the broader ornithological community and those interested in studying the evolutionary history and ecological interactions of larks, which comprise a widespread, yet understudied lineage of songbirds.

Phenotypic and genetic variation are present in all species, but lineages differ in how variation is partitioned among populations. Examining phenotypic clustering and genetic structure within a phylogeographic framework can clarify which biological processes have contributed to extant biodiversity in a given lineage. Here, we investigate genetic and phenotypic variation among populations and subspecies within a Neotropical songbird complex, the White-collared Seedeater (Sporophila torqueola) of Central America and Mexico. We combine measurements of morphology and plumage patterning with thousands of nuclear loci derived from ultraconserved elements (UCEs) and mitochondrial DNA to evaluate population differentiation. We find deep levels of molecular divergence between two S. torqueola lineages that are phenotypically diagnosable: One corresponds to S. t. torqueola along the Pacific coast of Mexico, and the other includes S. t. morelleti and S. t. sharpei from the Gulf Coast of Mexico and Central America. Surprisingly, these two lineages are strongly differentiated in both nuclear and mitochondrial markers, and each is more closely related to other Sporophila species than to one another. We infer low levels of gene flow between these two groups based on demographic models, suggesting multiple independent evolutionary lineages within S. torqueola have been obscured by coarse-scale similarity in plumage patterning. These findings improve our understanding of the biogeographic history of this lineage, which includes multiple dispersal events out of South America and across the Isthmus of Tehuantepec into Mesoamerica. Finally, the phenotypic and genetic distinctiveness of the range-restricted S. t. torqueola highlights the Pacific Coast of Mexico as an important region of endemism and conservation priority.

Many animals use photoperiod cues to synchronize reproduction with environmental conditions and thereby improve their reproductive success. The circadian clock, which creates endogenous behavioral and physiological rhythms typically entrained to photoperiod, is well characterized at the molecular level. Recent work provided evidence for an association between Clock poly-Q length polymorphism and latitude and, within a population, an association with the date of laying and the length of the incubation period. Despite relatively high overall breeding synchrony, the timing of clutch initiation has a large impact on the fitness of swallows in the genus Tachycineta. We compared length polymorphism in the Clock poly-Q region among five populations from five different Tachycineta species that breed across a hemisphere-wide latitudinal gradient (Fig. 1). Clock poly-Q variation was not associated with latitude; however, there was an association between Clock poly-Q allele diversity and the degree of clutch size decline within breeding seasons. We did not find evidence for an association between Clock poly-Q variation and date of clutch initiation in for any of the five Tachycineta species, nor did we found a relationship between incubation duration and Clock genotype. Thus, there is no general association between latitude, breeding phenology, and Clock polymorphism in this clade of closely related birds.Figure 1Photos of Tachycineta swallows that were used in this study: A) T. bicolor from Ithaca, New York, B) T. leucorrhoa from Chascomús, Argentina, C) T. albilinea from Hill Bank, Belize, D) T. meyeni from Puerto Varas, Chile, and E) T. thalassina from Mono Lake, California, Photographers: B: Valentina Ferretti; A, C-E: David Winkler.

Recent studies of several species have reported a latitudinal cline in the circadian clock gene, Clock, which influences rhythms in both physiology and behavior. Latitudinal variation in this gene may hence reflect local adaptation to seasonal variation. In some bird populations, there is also an among-individual association between Clock poly-Q genotype and clutch initiation date and incubation period. We examined Clock poly-Q allele variation in the Barn Swallow (Hirundo rustica), a species with a cosmopolitan geographic distribution and considerable variation in life-history traits that may be influenced by the circadian clock. We genotyped Barn Swallows from five populations (from three subspecies) and compared variation at the Clock locus to that at microsatellite loci and mitochondrial DNA (mtDNA). We found very low variation in the Clock poly-Q region, as >96% of individuals were homozygous, and the two other alleles at this locus were globally rare. Genetic differentiation based on the Clock poly-Q locus was not correlated with genetic differentiation based on either microsatellite loci or mtDNA sequences. Our results show that high diversity in Clock poly-Q is not general across avian species. The low Clock variation in the background of heterogeneity in microsatellite and mtDNA loci in Barn Swallows may be an outcome of stabilizing selection on the Clock locus.

Theory suggests that different taxa having colonized a similar, challenging environment will show parallel or lineage-specific adaptations to shared selection pressures, but empirical examples of parallel evolution in independent taxa are exceedingly rare. We employed comparative genomics to identify parallel and lineage-specific responses to selection within and among four species of North American sparrows that represent four independent, post-Pleistocene colonization events by an ancestral, upland subspecies and a derived salt marsh specialist. We identified multiple cases of parallel adaptation in these independent comparisons following salt marsh colonization, including selection of 12 candidate genes linked to osmoregulation. In addition to detecting shared genetic targets of selection across multiple comparisons, we found many novel, species-specific signatures of selection, including evidence of selection of loci associated with both physiological and behavioral mechanisms of osmoregulation. Demographic reconstructions of all four species highlighted their recent divergence and small effective population sizes, as expected given their rapid radiation into saline environments. Our results highlight the interplay of both shared and lineage-specific selection pressures in the colonization of a biotically and abiotically challenging habitat and confirm theoretical expectations that steep environmental clines can drive repeated and rapid evolutionary diversification in birds.

Species radiations have long served as model systems in evolutionary biology.1,2 However, it has only recently become possible to study the genetic bases of the traits responsible for diversification and only in a small number of model systems.3 Here, we use genomes of 36 species of North, Central, and South American warblers to highlight the role of pigmentation genes-involved in melanin and carotenoid processing-in the diversification of this group. We show that agouti signaling protein (ASIP) and beta-carotene oxygenase 2 (BCO2) are predictably divergent between species that differ in the distribution of melanin and carotenoid in their plumages, respectively. Among species, sequence variation at ASIP broadly mirrors the species' phylogenetic history, consistent with repeated, independent mutations generating melanin-based variation. In contrast, BCO2 variation is highly discordant from the species tree, with evidence of cross-lineage introgression among species like the yellow warbler (Setophaga petechia) and magnolia warbler (S. magnolia) with extensive carotenoid-based coloration. We also detect introgression of a small part of the BCO2 coding region (<3 kb) in S. discolor and S. vitellina, including an amino acid substitution that is unique to warblers but otherwise highly conserved across birds. Lateral transfer of carotenoid-processing genes has been documented in arthropods, but introgression of BCO2 as demonstrated here-presumably adaptive-represents the first example of carotenoid gene transfer among vertebrates. These contrasting genomic patterns show that both independent evolution in a common set of genes and past hybridization have fueled plumage diversification in this colorful avian radiation.

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity1-4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference5, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.