Host-derived population genomics data provides insights into bacterial and diatom composition of the killer whale skin

Recent exploration into the interactions and relationship between hosts and their microbiota has revealed a connection between many aspects of the host's biology, health and associated micro-organisms. Whereas amplicon sequencing has traditionally been used to characterize the microbiome, the increasing number of published population genomics data sets offers an underexploited opportunity to study microbial profiles from the host shotgun sequencing data. Here, we use sequence data originally generated from killer whale Orcinus orca skin biopsies for population genomics, to characterize the skin microbiome and investigate how host social and geographical factors influence the microbial community composition. Having identified 845 microbial taxa from 2.4 million reads that did not map to the killer whale reference genome, we found that both ecotypic and geographical factors influence community composition of killer whale skin microbiomes. Furthermore, we uncovered key taxa that drive the microbiome community composition and showed that they are embedded in unique networks, one of which is tentatively linked to diatom presence and poor skin condition. Community composition differed between Antarctic killer whales with and without diatom coverage, suggesting that the previously reported episodic migrations of Antarctic killer whales to warmer waters associated with skin turnover may control the effects of potentially pathogenic bacteria such as Tenacibaculum dicentrarchi. Our work demonstrates the feasibility of microbiome studies from host shotgun sequencing data and highlights the importance of metagenomics in understanding the relationship between host and microbial ecology.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci.

Behavioral Isolation and Incipient Speciation in Birds

Behavioral changes, such as those involved in mating, foraging, and migration, can generate reproductive barriers between populations. Birds, in particular, are known for their great diversity in these behaviors, and so behavioral isolation is often proposed to be the major driver of speciation. Here, we review empirical evidence to evaluate the importance of behavioral isolation in the early stages of avian speciation. Experimentally measured mating preferences indicate that changes in mating behavior can result in premating barriers, with their strength depending on the extent of divergence in mating signals. Differences in migratory and foraging behavior also can play important roles in generating reproductive barriers in the early stages of speciation. However, because premating behavioral isolation is imperfect, extrinsic postzygotic barriers, in the form of selection against hybrids having intermediate phenotypes, also play an important role in avian diversification, especially in completing the speciation process.

Illumina Library Preparation for Sequencing the GC-Rich Fraction of Heterogeneous Genomic DNA

Standard Illumina libraries are biased toward sequences of intermediate GC-content. This results in an underrepresentation of GC-rich regions in sequencing projects of genomes with heterogeneous base composition, such as mammals and birds. We developed a simple, cost-effective protocol to enrich sheared genomic DNA in its GC-rich fraction by subtracting AT-rich DNA. This was achieved by heating DNA up to 90 °C before applying Illumina library preparation. We tested the new approach on chicken DNA and found that heated DNA increased average coverage in the GC-richest chromosomes by a factor up to six. Using a Taq polymerase supposedly appropriate for PCR amplification of GC-rich sequences had a much weaker effect. Our protocol should greatly facilitate sequencing and resequencing of the GC-richest regions of heterogeneous genomes, in combination with standard short-read and long-read technologies.

Selection and gene flow shape genomic islands that control floral guides

Populations often show “islands of divergence” in the genome. Analysis of divergence between subspecies of Antirrhinum that differ in flower color patterns shows that sharp peaks in relative divergence occur at two causal loci. The island is shaped by a combination of gene flow and multiple selective sweeps, showing how divergence and barriers between populations can arise and be maintained.

Genomes of closely-related species or populations often display localized regions of enhanced relative sequence divergence, termed genomic islands. It has been proposed that these islands arise through selective sweeps and/or barriers to gene flow. Here, we genetically dissect a genomic island that controls flower color pattern differences between two subspecies of Antirrhinum majus, A.m.striatum and A.m.pseudomajus, and relate it to clinal variation across a natural hybrid zone. We show that selective sweeps likely raised relative divergence at two tightly-linked MYB-like transcription factors, leading to distinct flower patterns in the two subspecies. The two patterns provide alternate floral guides and create a strong barrier to gene flow where populations come into contact. This barrier affects the selected flower color genes and tightly-linked loci, but does not extend outside of this domain, allowing gene flow to lower relative divergence for the rest of the chromosome. Thus, both selective sweeps and barriers to gene flow play a role in shaping genomic islands: sweeps cause elevation in relative divergence, while heterogeneous gene flow flattens the surrounding “sea,” making the island of divergence stand out. By showing how selective sweeps establish alternative adaptive phenotypes that lead to barriers to gene flow, our study sheds light on possible mechanisms leading to reproductive isolation and speciation.

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.

A non-coding region near Follistatin controls head colour polymorphism in the Gouldian finch

Discrete colour morphs coexisting within a single population are common in nature. In a broad range of organisms, sympatric colour morphs often display major differences in other traits, including morphology, physiology or behaviour. Despite the repeated occurrence of this phenomenon, our understanding of the genetics that underlie multi-trait differences and the factors that promote the long-term maintenance of phenotypic variability within a freely interbreeding population are incomplete. Here, we investigated the genetic basis of red and black head colour in the Gouldian finch (Erythrura gouldiae), a classic polymorphic system in which naturally occurring colour morphs also display differences in aggressivity and reproductive success. We show that the candidate locus is a small (approx. 70 kb) non-coding region mapping to the Z chromosome near the Follistatin (FST) gene. Unlike recent findings in other systems where phenotypic morphs are explained by large inversions containing hundreds of genes (so-called supergenes), we did not identify any structural rearrangements between the two haplotypes using linked-read sequencing technology. Nucleotide divergence between the red and black alleles was high when compared to the remainder of the Z chromosome, consistent with their maintenance as balanced polymorphisms over several million years. Our results illustrate how pleiotropic phenotypes can arise from simple genetic variation, probably regulatory in nature.

Don't throw out the sympatric speciation with the crater lake water: fine-scale investigation of introgression provides equivocal support for causal role of secondary gene flow in one of the clearest examples of sympatric speciation

Genomic data has revealed complex histories of colonization and repeated gene flow previously unrecognized in some of the most celebrated examples of sympatric speciation and radiation. However, much of the evidence for secondary gene flow into these radiations comes from summary statistics calculated from sparse genomic sampling without knowledge of which specific genomic regions introgressed. This tells us little about how gene flow potentially influenced sympatric diversification. Here, we investigated whole genomes of Barombi Mbo crater lake cichlids for fine-scale patterns of introgression with neighboring riverine cichlid populations. We found evidence of secondary gene flow into the radiation scattered across <0.24% of the genome; however, from our analyses, it is not clear if the functional diversity in these regions contributed to the ecological, sexual, and morphological diversity found in the lake. Unlike similar studies, we found no obvious candidate genes for adaptive introgression and we cannot rule out that secondary gene flow was predominantly neutral with respect to the diversification process. We also found evidence for differential assortment of ancestral polymorphisms found in riverine populations between sympatric sister species, suggesting the presence of an ancestral hybrid swarm. Although the history of gene flow and colonization is more complicated than previously assumed, the lack of compelling evidence for secondary gene flow's role in species diversification suggests that we should not yet rule out one of the most celebrated examples of sympatric speciation in nature without a more thorough investigation of the timing and functional role of each introgressed region.

Comparative genomics and transcriptomics of Chrysolophus provide insights into the evolution of complex plumage coloration

Background

As one of the most recognizable characteristics in birds, plumage color has a high impact on understanding the evolution and mechanisms of coloration. Feather and skin are ideal tissues to explore the genomics and complexity of color patterns in vertebrates. Two species of the genus Chrysolophus, golden pheasant (Chrysolophus pictus) and Lady Amherst's pheasant (Chrysolophus amherstiae), exhibit brilliant colors in their plumage, but with extreme phenotypic differences, making these two species great models to investigate plumage coloration mechanisms in birds.

Results

We sequenced and assembled a genome of golden pheasant with high coverage and annotated 15,552 protein-coding genes. The genome of Lady Amherst's pheasant is sequenced with low coverage. Based on the feather pigment identification, a series of genomic and transcriptomic comparisons were conducted to investigate the complex features of plumage coloration. By identifying the lineage-specific sequence variations in Chrysolophus and golden pheasant against different backgrounds, we found that four melanogenesis biosynthesis genes and some lipid-related genes might be candidate genomic factors for the evolution of melanin and carotenoid pigmentation, respectively. In addition, a study among 47 birds showed some candidate genes related to carotenoid coloration in a broad range of birds. The transcriptome data further reveal important regulators of the two colorations, particularly one splicing transcript of the microphthalmia-associated transcription factor gene for pheomelanin synthesis.

Conclusions

Analysis of the golden pheasant and its sister pheasant genomes, as well as comparison with other avian genomes, are helpful to reveal the underlying regulation of their plumage coloration. The present study provides important genomic information and insights for further studies of avian plumage evolution and diversity.

Genetic and Molecular Basis of Feather Diversity in Birds

Feather diversity is striking in many aspects. Although the development of feather has been studied for decades, genetic and genomic studies of feather diversity have begun only recently. Many questions remain to be answered by multidisciplinary approaches. In this review, we discuss three levels of feather diversity: Feather morphotypes, intraspecific variations, and interspecific variations. We summarize recent studies of feather evolution in terms of genetics, genomics, and developmental biology and provide perspectives for future research. Specifically, this review includes the following topics: 1) Diversity of feather morphotype; 2) feather diversity among different breeds of domesticated birds, including variations in pigmentation pattern, in feather length or regional identity, in feather orientation, in feather distribution, and in feather structure; and 3) diversity of feathers among avian species, including plumage color and morph differences between species and the regulatory differences in downy feather development between altricial and precocial birds. Finally, we discussed future research directions.

Bidirectional adaptive introgression between two ecologically divergent sparrow species

Natural hybrid zones can be used to dissect the mechanisms driving key evolutionary processes by allowing us to identify genomic regions important for establishing reproductive isolation and that allow for transfer of adaptive variation. We leverage whole-genome data in a system where two bird species, the saltmarsh (Ammospiza caudacuta) and Nelson's (A. nelsoni) sparrow, hybridize despite their relatively high background genetic differentiation and past ecological divergence. Adaptive introgression is plausible in this system because Nelson's sparrows are recent colonists of saltwater marshes, compared to the specialized saltmarsh sparrow that has a longer history of saltmarsh adaptation. Comparisons among whole-genome sequences of 34 individuals from allopatric and sympatric populations show that ongoing gene flow is shaping the genomic landscape, with allopatric populations exhibiting genome-wide F_ST estimates close to double of that observed in sympatry. We characterized patterns of introgression across the genome and identify regions that exhibit biased introgression into hybrids from one parental species. These regions offer compelling candidates for genes related to tidal marsh adaptations suggesting that adaptive introgression may be an important consequence of hybridization. These findings highlight the value of considering the landscapes of both genome-wide introgression and divergence when characterizing the evolutionary forces that drive speciation.

Resource stability and geographic isolation are associated with genome divergence in western Palearctic crossbills

While many conifers produce annually variable seed crops, serotinous species (which hold seeds in cones for multiple years) represent unusually stable food resources for seed predators. Such stability is conducive to residency and potentially population divergence of consumers as exemplified by the Cassia crossbill (Loxia sinesciuris) in North America. We used genotyping by sequencing (GBS) to test whether three Mediterranean subspecies of common crossbills (L. curvirostra) associated with the serotinous Aleppo pine (Pinus halepensis) were more genetically distinct than European crossbills associated with nonserotinous conifers. We assembled a Cassia crossbill draft genome as a reference for mapping GBS reads and as a first step towards a more contiguous genome assembly. We found clear patterns of genetic divergence for each of the P. halepensis-associated subspecies. Geographic isolation, as promoted by resource stability and residency, is associated with genetic divergence of two of these subspecies. However, geographic isolation cannot account for divergence of L. c. hispana. Instead, resource stability likely contributed to divergence by reducing dispersal and increasing resource competition that may limit breeding by immigrants. In contrast, we found no differentiation among common crossbills associated with less stable resources, and only slight differentiation between common crossbills and parrot crossbills (L. pytyopsittacus). The substantial morphological divergence between common and parrot crossbills has likely originated or been maintained by selection despite gene flow generated by spatiotemporal resource fluctuation. Our results indicate that phenological as well as morphological characteristics of conifers have influenced crossbill diversification, and suggest a possible link between resource stability and population divergence.

Field studies reveal a close relative of C. elegans thrives in the fresh figs of Ficus septica and disperses on its Ceratosolen pollinating wasps

BACKGROUND

Biotic interactions are ubiquitous and require information from ecology, evolutionary biology, and functional genetics in order to be understood. However, study systems that are amenable to investigations across such disparate fields are rare. Figs and fig wasps are a classic system for ecology and evolutionary biology with poor functional genetics; Caenorhabditis elegans is a classic system for functional genetics with poor ecology. In order to help bridge these disciplines, here we describe the natural history of a close relative of C. elegans, Caenorhabditis inopinata, that is associated with the fig Ficus septica and its pollinating Ceratosolen wasps.

RESULTS

To understand the natural context of fig-associated Caenorhabditis, fresh F. septica figs from four Okinawan islands were sampled, dissected, and observed under microscopy. C. inopinata was found in all islands where F. septica figs were found. C.i nopinata was routinely found in the fig interior and almost never observed on the outside surface. C. inopinata was only found in pollinated figs, and C. inopinata was more likely to be observed in figs with more foundress pollinating wasps. Actively reproducing C. inopinata dominated early phase figs, whereas late phase figs with emerging wasp progeny harbored C. inopinata dauer larvae. Additionally, C. inopinata was observed dismounting from Ceratosolen pollinating wasps that were placed on agar plates. C. inopinata was not found on non-pollinating, parasitic Philotrypesis wasps. Finally, C. inopinata was only observed in F. septica figs among five Okinawan Ficus species sampled.

CONCLUSION

These are the first detailed field observations of C. inopinata, and they suggest a natural history where this species proliferates in early phase F. septica figs and disperses from late phase figs on Ceratosolen pollinating fig wasps. While consistent with other examples of nematode diversification in the fig microcosm, the fig and wasp host specificity of C. inopinata is highly divergent from the life histories of its close relatives and frames hypotheses for future investigations. This natural co-occurrence of the fig/fig wasp and C. inopinata study systems sets the stage for an integrated research program that can help to explain the evolution of interspecific interactions.

A High-Quality, Long-Read De Novo Genome Assembly to Aid Conservation of Hawaii's Last Remaining Crow Species

Genome-level data can provide researchers with unprecedented precision to examine the causes and genetic consequences of population declines, which can inform conservation management. Here, we present a high-quality, long-read, de novo genome assembly for one of the world’s most endangered bird species, the ʻAlalā (Corvus hawaiiensis; Hawaiian crow). As the only remaining native crow species in Hawaiʻi, the ʻAlalā survived solely in a captive-breeding program from 2002 until 2016, at which point a long-term reintroduction program was initiated. The high-quality genome assembly was generated to lay the foundation for both comparative genomics studies and the development of population-level genomic tools that will aid conservation and recovery efforts. We illustrate how the quality of this assembly places it amongst the very best avian genomes assembled to date, comparable to intensively studied model systems. We describe the genome architecture in terms of repetitive elements and runs of homozygosity, and we show that compared with more outbred species, the ʻAlalā genome is substantially more homozygous. We also provide annotations for a subset of immunity genes that are likely to be important in conservation management, and we discuss how this genome is currently being used as a roadmap for downstream conservation applications.

Speciation with gene flow via cycles of isolation and migration: insights from multiple mangrove taxa

Allopatric speciation requiring an unbroken period of geographical isolation has been the standard model of neo-Darwinism. While doubts have been repeatedly raised, strict allopatry without any gene flow remains a plausible mechanism in most cases. To rigorously reject strict allopatry, genomic sequences superimposed on the geological records of a well-delineated geographical barrier are necessary. The Strait of Malacca, narrowly connecting the Pacific and Indian Ocean coasts, serves at different times either as a geographical barrier or a conduit of gene flow for coastal/marine species. We surveyed 1700 plants from 29 populations of 5 common mangrove species by large-scale DNA sequencing and added several whole-genome assemblies. Speciation between the two oceans is driven by cycles of isolation and gene flow due to the fluctuations in sea level leading to the opening/closing of the Strait to ocean currents. Because the time required for speciation in mangroves is longer than the isolation phases, speciation in these mangroves has proceeded through many cycles of mixing-isolation-mixing, or MIM, cycles. The MIM mechanism, by relaxing the condition of no gene flow, can promote speciation in many more geographical features than strict allopatry can. Finally, the MIM mechanism of speciation is also efficient, potentially yielding mⁿ (m > 1) species after n cycles.

Repeated Selection of Alternatively Adapted Haplotypes Creates Sweeping Genomic Remodeling in Stickleback

Heterogeneous genetic divergence can accumulate across the genome when populations adapt to different habitats while still exchanging alleles. How long does diversification take and how much of the genome is affected? When divergence occurs in parallel from standing genetic variation, how often are the same haplotypes involved? We explore these questions using restriction site-associated DNA sequencing genotyping data and show that broad-scale genomic repatterning, fueled by copious standing variation, can emerge in just dozens of generations in replicate natural populations of threespine stickleback fish (Gasterosteus aculeatus). After the catastrophic 1964 Alaskan earthquake, marine stickleback colonized newly created ponds on seismically uplifted islands. We find that freshwater fish in these young ponds differ from their marine ancestors across the same genomic segments previously shown to have diverged in much older lake populations. Outside of these core divergent regions the genome shows no population structure across the ocean-freshwater divide, consistent with strong local selection acting in alternative environments on stickleback populations still connected by significant gene flow. Reinforcing this inference, a majority of divergent haplotypes that are at high frequency in ponds are detectable in the sea, even across great geographic distances. Building upon previous population genomics work in this model species, our data suggest that a long history of divergent selection and gene flow among stickleback populations in oceanic and freshwater habitats has maintained polymorphisms of alternatively adapted DNA sequences that facilitate parallel evolution.

Whole genome and transcriptome maps of the entirely black native Korean chicken breed Yeonsan Ogye

Background

Yeonsan Ogye (YO), an indigenous Korean chicken breed (Gallus gallus domesticus), has entirely black external features and internal organs. In this study, the draft genome of YO was assembled using a hybrid de novo assembly method that takes advantage of high-depth Illumina short reads (376.6X) and low-depth Pacific Biosciences (PacBio) long reads (9.7X).

Findings

The contig and scaffold NG50s of the hybrid de novo assembly were 362.3 Kbp and 16.8 Mbp, respectively. The completeness (97.6%) of the draft genome (Ogye_1.1) was evaluated with single-copy orthologous genes using Benchmarking Universal Single-Copy Orthologs and found to be comparable to the current chicken reference genome (galGal5; 97.4%; contigs were assembled with high-depth PacBio long reads (50X) and scaffolded with short reads) and superior to other avian genomes (92%-93%; assembled with short read-only or hybrid methods). Compared to galGal4 and galGal5, the draft genome included 551 structural variations including the fibromelanosis (FM) locus duplication, related to hyperpigmentation. To comprehensively reconstruct transcriptome maps, RNA sequencing and reduced representation bisulfite sequencing data were analyzed from 20 tissues, including 4 black tissues (skin, shank, comb, and fascia). The maps included 15,766 protein-coding and 6,900 long noncoding RNA genes, many of which were tissue-specifically expressed and displayed tissue-specific DNA methylation patterns in the promoter regions.

Conclusions

We expect that the resulting genome sequence and transcriptome maps will be valuable resources for studying domestic chicken breeds, including black-skinned chickens, as well as for understanding genomic differences between breeds and the evolution of hyperpigmented chickens and functional elements related to hyperpigmentation.

Extensive phenotypic diversification coexists with little genetic divergence and a lack of population structure in the White Wagtail subspecies complex (Motacilla alba)

Geographically clustered phenotypes often demonstrate consistent patterns in molecular markers, particularly mitochondrial DNA (mtDNA) traditionally used in phylogeographic studies. However, distinct evolutionary trajectories among traits and markers can lead to their discordance. First, geographic structure in phenotypic traits and nuclear molecular markers can be co-aligned but inconsistent with mtDNA (mito-nuclear discordance). Alternatively, phenotypic variation can have little to do with patterns in neither mtDNA nor nuclear markers. Disentangling between these distinct patterns can provide insight into the role of selection, demography and gene flow in population divergence. Here, we examined a previously reported case of strong inconsistency between geographic structure in mtDNA and plumage traits in a widespread polytypic bird species, the White Wagtail (Motacilla alba). We tested whether this pattern is due to mito-nuclear discordance or discrepancy between morphological evolution and both nuclear and mtDNA markers. We analysed population differentiation and structure across six out of nine commonly recognized subspecies using 17 microsatellite loci and a combination of microsatellites and plumage indices in a comprehensively sampled region of a contact between two subspecies. We did not find support for the mito-nuclear discordance hypothesis: nuclear markers indicated a subtle signal of genetic clustering only partially consistent with plumage groups, similar to previous findings that relied on mtDNA. We discuss evolutionary factors that could have shaped the intricate patterns of phenotypic diversification in the White wagtail and the role that repeated selection on plumage 'hotspots' and hybridization may have played.

Adaptive introgression underlies polymorphic seasonal camouflage in snowshoe hares

Snowshoe hares (Lepus americanus) maintain seasonal camouflage by molting to a white winter coat, but some hares remain brown during the winter in regions with low snow cover. We show that cis-regulatory variation controlling seasonal expression of the Agouti gene underlies this adaptive winter camouflage polymorphism. Genetic variation at Agouti clustered by winter coat color across multiple hare and jackrabbit species, revealing a history of recurrent interspecific gene flow. Brown winter coats in snowshoe hares likely originated from an introgressed black-tailed jackrabbit allele that has swept to high frequency in mild winter environments. These discoveries show that introgression of genetic variants that underlie key ecological traits can seed past and ongoing adaptation to rapidly changing environments.

Transgressive phenotypes and evidence of weak postzygotic isolation in F1 hybrids between closely related capuchino seedeaters

Postzygotic reproductive isolation may become strong only once the process of speciation is in its advanced stages. For taxa in the early stages of speciation, prezygotic reproductive isolation barriers may play a predominant role in maintaining species boundaries. Here, we study the recent capuchino seedeater biological radiation, a group of highly sympatric species from the genus Sporophila that have diversified during the Pleistocene in Neotropical grasslands. Capuchinos can be diagnosed by adult male coloration patterns and song, two sets of characters known to contribute to pre-mating reproductive isolation. However, it remains unknown whether potzygotic incompatibilities contribute to maintaining species limits in this group. Here we use existing breeding records from captive individuals to test for patterns consistent with F1 inviability. We compare hatching success, fledging success, and the sex ratio at adulthood between conspecific and hybrid capuchino pairs. We observed a trend towards lower numbers of the heterogametic sex among adult hybrids, consistent Haldane's rule, but this was supported by only one of our statistical tests. Our study is the first to document hybrid male capuchino phenotypes based on known crosses. We observed phenotypes that were similar or intermediate to those of the parental species, as well as novel plumage patterns that have not been described in the wild. One cross produced a plumage pattern that has been observed at low frequencies in natural populations. We discuss the implications of our results for understanding the relative importance of the mechanisms of reproductive isolation in capuchino seedeaters.

Varied Genomic Responses to Maladaptive Gene Flow and Their Evidence

Adaptation to a local environment often occurs in the face of maladaptive gene flow. In this perspective, I discuss several ideas on how a genome may respond to maladaptive gene flow during adaptation. On the one hand, selection can build clusters of locally adaptive alleles at fortuitously co-localized loci within a genome, thereby facilitating local adaptation with gene flow ('allele-only clustering'). On the other hand, the selective pressure to link adaptive alleles may drive co-localization of the actual loci relevant for local adaptation within a genome through structural genome changes or an evolving intra-genomic crossover rate ('locus clustering'). While the expected outcome is, in both cases, a higher frequency of locally adaptive alleles in some genome regions than others, the molecular units evolving in response to gene flow differ (i.e., alleles versus loci). I argue that, although making this distinction is important, we commonly lack the critical empirical evidence to do so. This is mainly because many current approaches are biased towards detecting local adaptation in genome regions with low crossover rates. The importance of low-crossover genome regions for adaptation with gene flow, such as in co-localizing relevant loci within a genome, thus remains unclear. Future empirical investigations should address these questions by making use of comparative genomics, where multiple de novo genome assemblies from species evolved under different degrees of genetic exchange are compared. This research promises to advance our understanding of how a genome adapts to maladaptive gene flow, thereby promoting adaptive divergence and reproductive isolation.

Admixture mapping in a hybrid zone reveals loci associated with avian feather coloration

Identifying the genetic bases for colour patterns has provided important insights into the control and expression of pigmentation and how these characteristics influence fitness. However, much more is known about the genetic bases for traits based on melanin pigments than for traits based on another major class of pigments, carotenoids. Here, we use natural admixture in a hybrid zone between Audubon's and myrtle warblers (Setophaga coronata auduboni/S. c. coronata) to identify genomic regions associated with both types of pigmentation. Warblers are known for rapid speciation and dramatic differences in plumage. For each of five plumage coloration traits, we found highly significant associations with multiple single-nucleotide polymorphisms (SNPs) across the genome and these were clustered in discrete regions. Regions near significantly associated SNPs were enriched for genes associated with keratin filaments, fibrils that make up feathers. A carotenoid-based trait that differs between the taxa-throat colour-had more than a dozen genomic regions of association. One cluster of SNPs for this trait overlaps the Scavenger Receptor Class F Member 2 (SCARF2) gene. Other scavenger receptors are presumed to be expressed at target tissues and involved in the selective movement of carotenoids into the target cells, making SCARF2 a plausible new candidate for carotenoid processing. In addition, two melanin-based plumage traits-colours of the eye line and eye spot-show very strong associations with a single genomic region mapping to chromosome 20 in the zebra finch. These findings indicate that only a subset of the genomic regions differentiated between these two warblers are associated with the plumage differences between them and demonstrate the utility of reduced-representation genomic scans in hybrid zones.

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Genomewide signatures of selection in Epichloë reveal candidate genes for host specialization

Host specialization is a key process in ecological divergence and speciation of plant-associated fungi. The underlying determinants of host specialization are generally poorly understood, especially in endophytes, which constitute one of the most abundant components of the plant microbiome. We addressed the genetic basis of host specialization in two sympatric subspecies of grass-endophytic fungi from the Epichloë typhina complex: subsp. typhina and clarkii. The life cycle of these fungi entails unrestricted dispersal of gametes and sexual reproduction before infection of a new host, implying that the host imposes a selective barrier on viability of the progeny. We aimed to detect genes under divergent selection between subspecies, experiencing restricted gene flow due to adaptation to different hosts. Using pooled whole-genome sequencing data, we combined F_ST and D_XY population statistics in genome scans and detected 57 outlier genes showing strong differentiation between the two subspecies. Genomewide analyses of nucleotide diversity (π), Tajima's D and dN/dS ratios indicated that these genes have evolved under positive selection. Genes encoding secreted proteins were enriched among the genes showing evidence of positive selection, suggesting that molecular plant-fungus interactions are strong drivers of endophyte divergence. We focused on five genes encoding secreted proteins, which were further sequenced in 28 additional isolates collected across Europe to assess genetic variation in a larger sample size. Signature of positive selection in these isolates and putative identification of pathogenic function supports our findings that these genes represent strong candidates for host specialization determinants in Epichloë endophytes. Our results highlight the role of secreted proteins as key determinants of host specialization.

Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography

Sex chromosomes are disproportionately involved in reproductive isolation and adaptation. In support of such a “large‐X” effect, genome scans between recently diverged populations and species pairs often identify distinct patterns of divergence on the sex chromosome compared to autosomes. When measures of divergence between populations are higher on the sex chromosome compared to autosomes, such patterns could be interpreted as evidence for faster divergence on the sex chromosome, that is “faster‐X”, barriers to gene flow on the sex chromosome. However, demographic changes can strongly skew divergence estimates and are not always taken into consideration. We used 224 whole‐genome sequences representing 36 populations from two Heliconius butterfly clades (H. erato and H. melpomene) to explore patterns of Z chromosome divergence. We show that increased divergence compared to equilibrium expectations can in many cases be explained by demographic change. Among Heliconius erato populations, for instance, population size increase in the ancestral population can explain increased absolute divergence measures on the Z chromosome compared to the autosomes, as a result of increased ancestral Z chromosome genetic diversity. Nonetheless, we do identify increased divergence on the Z chromosome relative to the autosomes in parapatric or sympatric species comparisons that imply postzygotic reproductive barriers. Using simulations, we show that this is consistent with reduced gene flow on the Z chromosome, perhaps due to greater accumulation of incompatibilities. Our work demonstrates the importance of taking demography into account to interpret patterns of divergence on the Z chromosome, but nonetheless provides evidence to support the Z chromosome as a strong barrier to gene flow in incipient Heliconius butterfly species.

Genomic data reveal cryptic lineage diversification and introgression in Californian golden cup oaks (section Protobalanus)

Here we study hybridization, introgression and lineage diversification in the widely distributed canyon live oak (Quercus chrysolepis) and the relict island oak (Q. tomentella), two Californian golden cup oaks with an intriguing biogeographical history. We employed restriction-site-associated DNA sequencing and integrated phylogenomic and population genomic analyses to study hybridization and reconstruct the evolutionary past of these taxa. Our analyses revealed the presence of two cryptic lineages within Q. chrysolepis. One of these lineages shares its most recent common ancestor with Q. tomentella, supporting the paraphyly of Q. chrysolepis. The split of these lineages was estimated to take place during the late Pliocene or the early Pleistocene, a time corresponding well with the common presence of Q. tomentella in the fossil records of continental California. Analyses also revealed historical hybridization among lineages, high introgression from Q. tomentella into Q. chrysolepis in their current area of sympatry, and widespread admixture between the two lineages of Q. chrysolepis in contact zones. Our results support that the two lineages of Q. chrysolepis behave as a single functional species phenotypically and ecologically well differentiated from Q. tomentella, a situation that can be only accommodated considering hybridization and speciation as a continuum with diffuse limits.

Parallel evolution of gene expression between trophic specialists despite divergent genotypes and morphologies

Parallel evolution of gene expression commonly underlies convergent niche specialization, but parallel changes in expression could also underlie divergent specialization. We investigated divergence in gene expression and whole-genome genetic variation across three sympatric Cyprinodon pupfishes endemic to San Salvador Island, Bahamas. This recent radiation consists of a generalist and two derived specialists adapted to novel niches: a scale-eating and a snail-eating pupfish. We sampled total mRNA from all three species at two early developmental stages and compared gene expression with whole-genome genetic differentiation among all three species in 42 resequenced genomes. Eighty percent of genes that were differentially expressed between snail-eaters and generalists were up or down regulated in the same direction between scale-eaters and generalists; however, there were no fixed variants shared between species underlying these parallel changes in expression. Genes showing parallel evolution of expression were enriched for effects on metabolic processes, whereas genes showing divergent expression were enriched for effects on cranial skeleton development and pigment biosynthesis, reflecting the most divergent phenotypes observed between specialist species. Our findings reveal that even divergent niche specialists may exhibit convergent adaptation to higher trophic levels through shared genetic pathways. This counterintuitive result suggests that parallel evolution in gene expression can accompany divergent ecological speciation during adaptive radiation.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

A genomic map of clinal variation across the European rabbit hybrid zone

Speciation is a process proceeding from weak to complete reproductive isolation. In this continuum, naturally hybridizing taxa provide a promising avenue for revealing the genetic changes associated with the incipient stages of speciation. To identify such changes between two subspecies of rabbits that display partial reproductive isolation, we studied patterns of allele frequency change across their hybrid zone using whole-genome sequencing. To connect levels and patterns of genetic differentiation with phenotypic manifestations of subfertility in hybrid rabbits, we further investigated patterns of gene expression in testis. Geographic cline analysis revealed 253 regions characterized by steep changes in allele frequency across their natural region of contact. This catalog of regions is likely to be enriched for loci implicated in reproductive barriers and yielded several insights into the evolution of hybrid dysfunction in rabbits: (i) incomplete reproductive isolation is likely governed by the effects of many loci, (ii) protein-protein interaction analysis suggest that genes within these loci interact more than expected by chance, (iii) regulatory variation is likely the primary driver of incompatibilities, and (iv) large chromosomal rearrangements appear not to be a major mechanism underlying incompatibilities or promoting isolation in the face of gene flow. We detected extensive misregulation of gene expression in testis of hybrid males, but not a statistical overrepresentation of differentially expressed genes in candidate regions. Our results also did not support an X chromosome-wide disruption of expression as observed in mice and cats, suggesting variation in the mechanistic basis of hybrid male reduced fertility among mammals.

Hidden endemism, deep polyphyly, and repeated dispersal across the Isthmus of Tehuantepec: Diversification of the White-collared Seedeater complex (Thraupidae: Sporophila torqueola)

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.

Biodiversity and the Species Concept-Lineages are not Enough

The nature and definition of species continue to be matters of debate. Current views of species often focus on their nature as lineages-maximal reproductive communities through time. Whereas many authors point to the Evolutionary Species Concept as optimal, in its original form it stressed the ecological role of species as well as their history as lineages, but most recent authors have ignored the role aspect of the concept, making it difficult to apply unambiguously in a time-extended way. This trend has been exacerbated by the application of methods and concepts emphasizing the notion of monophyly, originally applied only at higher levels, to the level of individuals, as well as by the current emphasis on molecular data. Hence, some current authors recognize units that are no more than probable exclusive lineages as species. We argue that biodiversity is inherently a phenotypic concept and that role, as manifested in the organismal extended phenotype, is a necessary component of the species concept. Viewing species as historically connected populations with unique role brings together the temporal and phenotypic natures of species, providing a clear way to view species both in a time-limited and time-extended way. Doing so alleviates perceived issues with "paraphyletic species" and returns the focus of species to units that are most relevant for biodiversity.

Genomic islands of divergence or opportunities for introgression?

In animals, introgression between species is often perceived as the breakdown of reproductive isolating mechanisms, but gene flow between incipient species can also represent a source for potentially beneficial alleles. Recently, genome-wide datasets have revealed clusters of differentiated loci ('genomic islands of divergence') that are thought to play a role in reproductive isolation and therefore have reduced gene flow. We use simulations to further examine the evolutionary forces that shape and maintain genomic islands of divergence between two subspecies of the migratory songbird, Swainson's thrush (Catharus ustulatus), which have come into secondary contact since the last glacial maximum. We find that, contrary to expectation, gene flow is high within islands and is highly asymmetric. In addition, patterns of nucleotide diversity at highly differentiated loci suggest selection was more frequent in a single ecotype. We propose a mechanism whereby beneficial alleles spread via selective sweeps following a post-glacial demographic expansion in one subspecies and move preferentially across the hybrid zone. We find no evidence that genomic islands are the result of divergent selection or reproductive isolation, rather our results suggest that differentiated loci both within and outside islands could provide opportunities for adaptive introgression across porous species boundaries.

Androgen and estrogen sensitivity of bird song: a comparative view on gene regulatory levels

Singing of songbirds is sensitive to testosterone and its androgenic and estrogenic metabolites in a species-specific way. The hormonal effects on song pattern are likely mediated by androgen receptors (AR) and estrogen receptor alpha (ERα), ligand activated transcription factors that are expressed in neurons of various areas of the songbirds' vocal control circuit. The distribution of AR in this circuit is rather similar between species while that of ERα is species variant and concerns a key vocal control area, the HVC (proper name). We discuss the regulation of the expression of the cognate AR and ERα and putative splice variants. In particular, we suggest that transcription factor binding sites in the promoter of these receptors differ between bird species. Further, we suggest that AR- and ERα-dependent gene regulation in vocal areas differs between species due to species-specific DNA binding sites of putative target genes that are required for the transcriptional activity of the receptors. We suggest that species differences in the distribution of AR and ERα in vocal areas and in the genomic sensitivity to these receptors contribute to species-specific hormonal regulation of the song.

Urban landscape genomics identifies fine-scale gene flow patterns in an avian invasive

Invasive species exert a serious impact on native fauna and flora and have been the target of many eradication and management efforts worldwide. However, a lack of data on population structure and history, exacerbated by the recency of many species introductions, limits the efficiency with which such species can be kept at bay. In this study we generated a novel genome of high assembly quality and genotyped 4735 genome-wide single nucleotide polymorphic (SNP) markers from 78 individuals of an invasive population of the Javan Myna Acridotheres javanicus across the island of Singapore. We inferred limited population subdivision at a micro-geographic level, a genetic patch size (~13-14 km) indicative of a pronounced dispersal ability, and barely an increase in effective population size since introduction despite an increase of four to five orders of magnitude in actual population size, suggesting that low population-genetic diversity following a bottleneck has not impeded establishment success. Landscape genomic analyses identified urban features, such as low-rise neighborhoods, that constitute pronounced barriers to gene flow. Based on our data, we consider an approach targeting the complete eradication of Javan Mynas across Singapore to be unfeasible. Instead, a mixed approach of localized mitigation measures taking into account urban geographic features and planning policy may be the most promising avenue to reducing the adverse impacts of this urban pest. Our study demonstrates how genomic methods can directly inform the management and control of invasive species, even in geographically limited datasets with high gene flow rates.

Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of divergence within a poplar species complex

One of the outstanding questions in understanding how new species form is how reproductive isolation arises. In particular, the relative roles of gene flow and natural selection in creating two separate species remains open for debate. Here we show within the four continuously speciating lineages of a poplar that local genomic differentiation of populations is not associated with either rate of recent gene flow or time of species divergence. By contrast, we found that these genomic islands of divergence most likely came about by selective processes—sorting of ancient genetic polymorphisms and the incidental hitchhiking of linked variations. These findings substantially enhance our understanding of genomic changes in speciation.

How genome divergence eventually leads to speciation is a topic of prime evolutionary interest. Genomic islands of elevated divergence are frequently reported between diverging lineages, and their size is expected to increase with time and gene flow under the speciation-with-gene-flow model. However, such islands can also result from divergent sorting of ancient polymorphisms, recent ecological selection regardless of gene flow, and/or recurrent background selection and selective sweeps in low-recombination regions. It is challenging to disentangle these nonexclusive alternatives, but here we attempt to do this in an analysis of what drove genomic divergence between four lineages comprising a species complex of desert poplar trees. Within this complex we found that two morphologically delimited species, Populus euphratica and Populus pruinosa, were paraphyletic while the four lineages exhibited contrasting levels of gene flow and divergence times, providing a good system for testing hypotheses on the origin of divergence islands. We show that the size and number of genomic islands that distinguish lineages are not associated with either rate of recent gene flow or time of divergence. Instead, they are most likely derived from divergent sorting of ancient polymorphisms and divergence hitchhiking. We found that highly diverged genes under lineage-specific selection and putatively involved in ecological and morphological divergence occur both within and outside these islands. Our results highlight the need to incorporate demography, absolute divergence measurement, and gene flow rate to explain the formation of genomic islands and to identify potential genomic regions involved in speciation.