The genomic mosaicism of hybrid speciation

Consistent scaling of inbreeding depression in space and time in a house sparrow metapopulation

Inbreeding may increase the extinction risk of small populations. Yet, studies using modern genomic tools to investigate inbreeding depression in nature have been limited to single populations, and little is known about the dynamics of inbreeding depression in subdivided populations over time. Natural populations often experience different environmental conditions and differ in demographic history and genetic composition, characteristics that can affect the severity of inbreeding depression. We utilized extensive long-term data on more than 3,100 individuals from eight islands in an insular house sparrow metapopulation to examine the generality of inbreeding effects. Using genomic estimates of realized inbreeding, we discovered that inbred individuals had lower survival probabilities and produced fewer recruiting offspring than noninbred individuals. Inbreeding depression, measured as the decline in fitness-related traits per unit inbreeding, did not vary appreciably among populations or with time. As a consequence, populations with more resident inbreeding (due to their demographic history) paid a higher total fitness cost, evidenced by a larger variance in fitness explained by inbreeding within these populations. Our results are in contrast to the idea that effects of inbreeding generally depend on ecological factors and genetic differences among populations, and expand the understanding of inbreeding depression in natural subdivided populations.

Recent hybrids recapitulate ancient hybrid outcomes

Genomic outcomes of hybridization depend on selection and recombination in hybrids. Whether these processes have similar effects on hybrid genome composition in contemporary hybrid zones versus ancient hybrid lineages is unknown. Here we show that patterns of introgression in a contemporary hybrid zone in Lycaeides butterflies predict patterns of ancestry in geographically adjacent, older hybrid populations. We find a particularly striking lack of ancestry from one of the hybridizing taxa, Lycaeides melissa, on the Z chromosome in both the old and contemporary hybrids. The same pattern of reduced L. melissa ancestry on the Z chromosome is seen in two other ancient hybrid lineages. More generally, we find that patterns of ancestry in old or ancient hybrids are remarkably predictable from contemporary hybrids, which suggests selection and recombination affect hybrid genomes in a similar way across disparate time scales and during distinct stages of speciation and species breakdown.

The house sparrow in the service of basic and applied biology

From the northernmost tip of Scandinavia to the southernmost corner of Patagonia, and across six continents, house sparrows (Passer domesticus) inhabit most human-modified habitats of the globe. With over 7,000 articles published, the species has become a workhorse for not only the study of self-urbanized wildlife, but also for understanding life history and body size evolution, sexual selection and many other biological phenomena. Traditionally, house sparrows were studied for their adaptations to local biotic and climatic conditions, but more recently, the species has come to serve as a focus for studies seeking to reveal the genomic, epigenetic and physiological underpinnings of success among invasive vertebrate species. Here, we review the natural history of house sparrows, highlight what the study of these birds has meant to bioscience generally, and describe the many resources available for future work on this species.

Genetic Variation Related to High Elevation Adaptation Revealed by Common Garden Experiments in Pinus yunnanensis

Local adaptation, adaptation to specialized niches and environmental clines have been extensively reported for forest trees. Investigation of the adaptive genetic variation is crucial for forest resource management and breeding, especially in the context of global climate change. Here, we utilized a Pinus yunnanensis common garden experiments established at high and low elevation sites to assess the differences in growth and survival among populations and between the two common garden sites. The studied traits showed significant variation between the two test sites and among populations, suggesting adaptive divergence. To detect genetic variation related to environment, we captured 103,608 high quality SNPs based on RNA sequencing, and used them to assess the genetic diversity and population structure. We identified 321 outlier SNPs from 131 genes showing significant divergence in allelic frequency between survival populations of two sites. Functional categories associated with adaptation to high elevation were found to be related to flavonoid biosynthesis, response to UV, DNA repair, response to reactive oxygen species, and membrane lipid metabolic process. Further investigation of the outlier genes showed overrepresentation of the flavonoid biosynthesis pathway, suggesting that this pathway may play a key role in P. yunnanensis adaptation to high elevation environments. The outlier genes identified, and their variants, provide a basic reference for advanced investigations.

The Genome of the Great Gerbil Reveals Species-Specific Duplication of an MHCII Gene

The great gerbil (Rhombomys opimus) is a social rodent living in permanent, complex burrow systems distributed throughout Central Asia, where it serves as the main host of several important vector-borne infectious pathogens including the well-known plague bacterium (Yersinia pestis). Here, we present a continuous annotated genome assembly of the great gerbil, covering over 96% of the estimated 2.47-Gb genome. Taking advantage of the recent genome assemblies of the sand rat (Psammomys obesus) and the Mongolian gerbil (Meriones unguiculatus), comparative immunogenomic analyses reveal shared gene losses within TLR gene families (i.e., TLR8, TLR10, and the entire TLR11-subfamily) for Gerbillinae, accompanied with signs of diversifying selection of TLR7 and TLR9. Most notably, we find a great gerbil-specific duplication of the MHCII DRB locus. In silico analyses suggest that the duplicated gene provides high peptide binding affinity for Yersiniae epitopes as well as Leishmania and Leptospira epitopes, putatively leading to increased capability to withstand infections by these pathogens. Our study demonstrates the power of whole-genome sequencing combined with comparative genomic analyses to gain deeper insight into the immunogenomic landscape of the great gerbil and its close relatives.

Eukaryote hybrid genomes

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve.

Dissimilar molecular and morphological patterns in an introgressed peripheral population of a sand dune species (Armeria pungens, Plumbaginaceae)

Introgression is a poorly understood evolutionary outcome of hybridisation because it may remain largely undetected whenever it involves the transfer of small parts of the genome from one species to another. Aiming to understand the early stages of this process, a putative case from the southernmost border of the Armeria pungens range from its congener A. macrophylla is revisited following the discovery of a subpopulation that does not show phenotypic signs of introgression and resembles typical A. pungens. We analysed morphometrics, nuclear ribosomal DNA ITS and plastid DNA (trnL-trnF) sequences, genome size, 45S and 5S rDNA loci-FISH data and nrDNA IGS sequences. Within the study site, most individuals match morphologies of either of the two hybridising species, particularly the new subpopulation, with intermediate phenotypes being scarce. This pattern does not fully fit molecular evidence revealing two ITS ribotypes co-occurring intragenomically in most plants from the study site and one single plastid haplotype. Genome size and structural features of the IGS sequences both indicate that A. pungens from the study site is genetically more similar to its sympatric congener than to the remainder of its conspecifics. Introgression of A. macrophylla into A. pungens and plastid capture explain all the evidence analysed. However, important features to understand the origin and fate of the introgressed population, such as the degree and direction of introgression, which are important for understanding early stages of hybridisation in plants with low reproductive barriers, should be addressed with new data.

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves, Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early-breeding house sparrows and hybrids and rural late-breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z-chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land-use changes in a mosaic landscape.

Hybridization increases population variation during adaptive radiation

Adaptive radiations are prominent components of the world's biodiversity. They comprise many species derived from one or a small number of ancestral species in a geologically short time that have diversified into a variety of ecological niches. Several authors have proposed that introgressive hybridization has been important in the generation of new morphologies and even new species, but how that happens throughout evolutionary history is not known. Interspecific gene exchange is expected to have greatest impact on variation if it occurs after species have diverged genetically and phenotypically but before genetic incompatibilities arise. We use a dated phylogeny to infer that populations of Darwin's finches in the Galápagos became more variable in morphological traits through time, consistent with the hybridization hypothesis, and then declined in variation after reaching a peak. Some species vary substantially more than others. Phylogenetic inferences of hybridization are supported by field observations of contemporary hybridization. Morphological effects of hybridization have been investigated on the small island of Daphne Major by documenting changes in hybridizing populations of Geospiza fortis and Geospiza scandens over a 30-y period. G. scandens showed more evidence of admixture than G. fortis Beaks of G. scandens became progressively blunter, and while variation in length increased, variation in depth decreased. These changes imply independent effects of introgression on 2, genetically correlated, beak dimensions. Our study shows how introgressive hybridization can alter ecologically important traits, increase morphological variation as a radiation proceeds, and enhance the potential for future evolution in changing environments.

Mitochondrial heteroplasmy in an avian hybrid form (Passer italiae: Aves, Passeriformes)

Mitochondrial heteroplasmy is the result from biparental transmission of mitochondrial DNA (mtDNA) to the offspring. In such rare cases, maternal and paternal mtDNA is present in the same individual. Though recent studies suggested that mtDNA heteroplasmy might be more common than previously anticipated, that phenomenon is still poorly documented and was mostly detected in case studies on hybrid populations. The Italian sparrow, Passer italiae is a homoploid hybrid form that occurs all across the Italian Peninsula mostly under strict absence of either of its parent species, the house sparrow (P. domesticus) and the Spanish sparrow (P. hispaniolensis). In this study, we document a new case of mitochondrial heteroplasmy from two island populations of P. italiae (Ustica and Lipari). Our analysis was based on the mitochondrial NADH dehydrogenase subunit 2 (ND2) that allows for a clear distinction between mitochondrial lineages of the two parental species. We amplified and sequenced the mitochondrial ND2 gene with specifically designed primer combinations for each of the two parental species. In two of our study populations, a single individual carried two different ND2 haplotypes from each of the two parental lineages. These findings contribute to current knowledge on the still poorly documented phenomenon of paternal leakage in vertebrates.

A Multireference-Based Whole Genome Assembly for the Obligate Ant-Following Antbird, Rhegmatorhina melanosticta (Thamnophilidae)

Current generation high-throughput sequencing technology has facilitated the generation of more genomic-scale data than ever before, thus greatly improving our understanding of avian biology across a range of disciplines. Recent developments in linked-read sequencing (Chromium 10×) and reference-based whole-genome assembly offer an exciting prospect of more accessible chromosome-level genome sequencing in the near future. We sequenced and assembled a genome of the Hairy-crested Antbird (Rhegmatorhina melanosticta), which represents the first publicly available genome for any antbird (Thamnophilidae). Our objectives were to (1) assemble scaffolds to chromosome level based on multiple reference genomes, and report on differences relative to other genomes, (2) assess genome completeness and compare content to other related genomes, and (3) assess the suitability of linked-read sequencing technology for future studies in comparative phylogenomics and population genomics studies. Our R. melanosticta assembly was both highly contiguous (de novo scaffold N50 = 3.3 Mb, reference based N50 = 53.3 Mb) and relatively complete (contained close to 90% of evolutionarily conserved single-copy avian genes and known tetrapod ultraconserved elements). The high contiguity and completeness of this assembly enabled the genome to be successfully mapped to the chromosome level, which uncovered a consistent structural difference between R. melanosticta and other avian genomes. Our results are consistent with the observation that avian genomes are structurally conserved. Additionally, our results demonstrate the utility of linked-read sequencing for non-model genomics. Finally, we demonstrate the value of our R. melanosticta genome for future researchers by mapping reduced representation sequencing data, and by accurately reconstructing the phylogenetic relationships among a sample of thamnophilid species.

Ancient introgression drives adaptation to cooler and drier mountain habitats in a cypress species complex

Introgression may act as an important source of new genetic variation to facilitate the adaptation of organisms to new environments, yet how introgression might enable tree species to adapt to higher latitudes and elevations remains unclear. Applying whole-transcriptome sequencing and population genetic analyses, we present an example of ancient introgression from a cypress species (Cupressus gigantea) that occurs at higher latitude and elevation on the Qinghai-Tibet Plateau into a related species (C. duclouxiana), which has likely aided the latter species to extend its range by colonizing cooler and drier mountain habitats during postglacial periods. We show that 16 introgressed candidate adaptive loci could have played pivotal roles in response to diverse stresses experienced in a high-elevation environment. Our findings provide new insights into the evolutionary history of Qinghai-Tibet Plateau plants and the importance of introgression in the adaptation of species to climate change.

Phylogenomics Reveals an Ancient Hybrid Origin of the Persian Walnut

Persian walnut (Juglans regia) is cultivated worldwide for its high-quality wood and nuts, but its origin has remained mysterious because in phylogenies it occupies an unresolved position between American black walnuts and Asian butternuts. Equally unclear is the origin of the only American butternut, J. cinerea. We resequenced the whole genome of 80 individuals from 19 of the 22 species of Juglans and assembled the genome of its relatives Pterocarya stenoptera and Platycarya strobilacea. Using phylogenetic-network analysis of single-copy nuclear genes, genome-wide site pattern probabilities, and Approximate Bayesian Computation, we discovered that J. regia (and its landrace J. sigillata) arose as a hybrid between the American and the Asian lineages and that J. cinerea resulted from massive introgression from an immigrating Asian butternut into the genome of an American black walnut. Approximate Bayesian Computation modeling placed the hybrid origin in the late Pliocene, ∼3.45 My, with both parental lineages since having gone extinct in Europe.

Reconciling the Mitonuclear Compatibility Species Concept with Rampant Mitochondrial Introgression

The mitonuclear compatibility species concept defines a species as a population that is genetically isolated from other populations by uniquely coadapted mitochondrial (mt) and nuclear genes. A key prediction of this hypothesis is that the mt genotype of each species will be functionally distinct and that introgression of mt genomes will be prevented by mitonuclear incompatibilities that arise when heterospecific mt and nuclear genes attempt to cofunction to enable aerobic respiration. It has been proposed, therefore, that the observation of rampant introgression of mt genotypes from one species to another constitutes a strong refutation of the mitonuclear speciation. The displacement of a mt genotype from a nuclear background with which it co-evolved to a foreign nuclear background will necessarily lead to fitness loss due to mitonuclear incompatibilities. Here I consider two potential benefits of mt introgression between species that may, in some cases, overcome fitness losses arising from mitonuclear incompatibilities. First, the introgressed mt genotype may be better adapted to the local environment than the native mt genotype such that higher fitness is achieved through improved adaptation via introgression. Second, if the mitochondria of the recipient taxa carry a high mutational load, then introgression of a foreign, less corrupt mt genome may enable the recipient taxa to escape its mutational load and gain a fitness advantage. Under both scenarios, fitness gains from novel mt genotypes could theoretically compensate for the fitness that is lost via mitonuclear incompatibility. I also consider the role of endosymbionts in non-adaptive rampant introgression of mt genomes. I conclude that rampant introgression is not necessarily evidence against the idea of tight mitonuclear coadaptation or the mitonuclear compatibility species concept. Rampant mt introgression will typically lead to erasure of species but in some cases could lead to hybrid speciation.

Mature sperm small-RNA profile in the sparrow: implications for transgenerational effects of age on fitness

Mammalian sperm RNA has recently received a lot of interest due to its involvement in epigenetic germline inheritance. Studies of epigenetic germline inheritance have shown that environmental exposures can induce effects in the offspring without altering the DNA sequence of germ cells. Most mechanistic studies were conducted in laboratory rodents and C.elegans while observational studies confirm the phenotypic phenomenon in wild populations of humans and other species including birds. Prominently, paternal age in house sparrows affects offspring fitness, yet the mechanism is unknown. This study provides a first reference of house sparrow sperm small RNA as an attempt to uncover their role in the transmission of the effects of paternal age on the offspring. In this small-scale pilot, we found no statistically significant differences between miRNA and tRNA fragments in aged and prime sparrow sperm. These results indicate a role of other epigenetic information carriers, such as distinct RNA classes, RNA modifications, DNA methylation and retained histones, and a clear necessity of future studies in wild populations.

Single-Cell Omics Analyses Enabled by Microchip Technologies

Single-cell omics studies provide unique information regarding cellular heterogeneity at various levels of the molecular biology central dogma. This knowledge facilitates a deeper understanding of how underlying molecular and architectural changes alter cell behavior, development, and disease processes. The emerging microchip-based tools for single-cell omics analysis are enabling the evaluation of cellular omics with high throughput, improved sensitivity, and reduced cost. We review state-of-the-art microchip platforms for profiling genomics, epigenomics, transcriptomics, proteomics, metabolomics, and multi-omics at single-cell resolution. We also discuss the background of and challenges in the analysis of each molecular layer and integration of multiple levels of omics data, as well as how microchip-based methodologies benefit these fields. Additionally, we examine the advantages and limitations of these approaches. Looking forward, we describe additional challenges and future opportunities that will facilitate the improvement and broad adoption of single-cell omics in life science and medicine.

Signs of adaptation to trace metal contamination in a common urban bird

Metals and metalloids at elevated concentrations can be toxic to both humans and wildlife. In particular, lead exposure can act as a stressor to wildlife and cause negative effects on fitness. Any ability to adapt to stress caused by the negative effects of trace metal exposure would be beneficial for species living in contaminated environments. However, mechanisms for responding adaptively to metal contamination are not fully understood in free-living organisms. The Australian populations of the house sparrow (Passer domesticus) provides an excellent opportunity to study potential adaptation to environmental lead contamination because they have a commensal relationship with humans and are distributed broadly across Australian settlements including many long-term mining and smelting communities. To examine the potential for an evolutionary response to long-term lead exposure, we collected genomic SNP data using the house sparrow 200 K SNP array, from 11 localities across the Australian distribution including two mining sites (Broken Hill and Mount Isa, which are two genetically independent populations) that have well-established elevated levels of lead contamination as well as trace metals and metalloids. We contrast these known contaminated locations to other lesser-contaminated environments. Using an ecological association genome scan method to identify genomic differentiation associated with estimates of lead contamination we identified 60 outlier loci across three tests. A total of 39 genes were found to be physically linked (within 20 kbps) of all outliers in the house sparrow reference genome. The linked candidate genes included 12 genes relevant to lead exposure, such as two metal transporters that can transport metals including lead and zinc across cell membranes. These candidate genes provide targets for follow up experiments comparing resilience to lead exposure between populations exposed to varied levels of lead contamination.

Predictable genome-wide sorting of standing genetic variation during parallel adaptation to basic versus acidic environments in stickleback fish

Genomic studies of parallel (or convergent) evolution often compare multiple populations diverged into two ecologically different habitats to search for loci repeatedly involved in adaptation. Because the shared ancestor of these populations is generally unavailable, the source of the alleles at adaptation loci, and the direction in which their frequencies were shifted during evolution, remain elusive. To shed light on these issues, we here use multiple populations of threespine stickleback fish adapted to two different types of derived freshwater habitats-basic and acidic lakes on the island of North Uist, Outer Hebrides, Scotland-and the present-day proxy of their marine ancestor. In a first step, we combine genome-wide pooled sequencing and targeted individual-level sequencing to demonstrate that ecological and phenotypic parallelism in basic-acidic divergence is reflected by genomic parallelism in dozens of genome regions. Exploiting data from the ancestor, we next show that the acidic populations, residing in ecologically more extreme derived habitats, have adapted by accumulating alleles rare in the ancestor, whereas the basic populations have retained alleles common in the ancestor. Genomic responses to selection are thus predictable from the ecological difference of each derived habitat type from the ancestral one. This asymmetric sorting of standing genetic variation at loci important to basic-acidic divergence has further resulted in more numerous selective sweeps in the acidic populations. Finally, our data suggest that the maintenance in marine fish of standing variation important to adaptive basic-acidic differentiation does not require extensive hybridization between the marine and freshwater populations. Overall, our study reveals striking genome-wide determinism in both the loci involved in parallel divergence, and in the direction in which alleles at these loci have been selected.

The Timing and Direction of Introgression Under the Multispecies Network Coalescent

Introgression is a pervasive biological process, and many statistical methods have been developed to infer its presence from genomic data. However, many of the consequences and genomic signatures of introgression remain unexplored from a methodological standpoint. Here, we develop a model for the timing and direction of introgression based on the multispecies network coalescent, and from it suggest new approaches for testing introgression hypotheses. We suggest two new statistics, D ₁ and D ₂, which can be used in conjunction with other information to test hypotheses relating to the timing and direction of introgression, respectively. D ₁ may find use in evaluating cases of homoploid hybrid speciation (HHS), while D ₂ provides a four-taxon test for polarizing introgression. Although analytical expectations for our statistics require a number of assumptions to be met, we show how simulations can be used to test hypotheses about introgression when these assumptions are violated. We apply the D ₁ statistic to genomic data from the wild yeast Saccharomyces paradoxus-a proposed example of HHS-demonstrating its use as a test of this model. These methods provide new and powerful ways to address questions relating to the timing and direction of introgression.

Exploring the hybrid speciation continuum in birds

Hybridization is increasingly recognized as a creative evolutionary force contributing to adaptation and speciation. Homoploid hybrid speciation-the process in which hybridization results in a stable, fertile, and reproductively isolated hybrid lineage where there is no change in ploidy-has been documented in several taxa. Hybridization can directly contribute to reproductive isolation or reinforce it at a later stage. Alternatively, hybridization might not be related to the evolution of reproductive isolation. To account for these different scenarios, I propose to discriminate between two types of hybrid speciation: type I where reproductive isolation is a direct consequence of hybridization and type II where it is the by-product of other processes. I illustrate the applicability of this classification scheme with avian examples. To my knowledge, seven hybrid bird species have been proposed: Italian sparrow, Audubon's warbler, Genovesa mockingbird, Hawaiian duck, red-breasted goose, golden-crowned manakin, and a recent lineage of Darwin's finches on the island of Daphne Major ("Big Bird"). All studies provide convincing evidence for hybridization, but do not always confidently discriminate between scenarios of hybrid speciation and recurrent introgressive hybridization. The build-up of reproductive isolation between the hybrid species and their parental taxa is mainly driven by premating isolation mechanisms and comparable to classical speciation events. One hybrid species can be classified as type I ("Big Bird") while three species constitute type II hybrid species (Italian sparrow, Audubon's warbler, and golden-crowned manakin). The diversity in hybrid bird species across a range of divergence times also provides an excellent opportunity to study the evolution of hybrid genomes in terms of genome stabilization and adaptation.

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.

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.

Signatures of human-commensalism in the house sparrow genome

House sparrows (Passer domesticus) are a hugely successful anthrodependent species; occurring on nearly every continent. Yet, despite their ubiquity and familiarity to humans, surprisingly little is known about their origins. We sought to investigate the evolutionary history of the house sparrow and identify the processes involved in its transition to a human-commensal niche. We used a whole genome resequencing dataset of 120 individuals from three Eurasian species, including three populations of Bactrianus sparrows, a non-commensal, divergent house sparrow lineage occurring in the Near East. Coalescent modelling supports a split between house and Bactrianus sparrow 11 Kya and an expansion in the house sparrow at 6 Kya, consistent with the spread of agriculture following the Neolithic revolution. Commensal house sparrows therefore likely moved into Europe with the spread of agriculture following this period. Using the Bactrianus sparrow as a proxy for a pre-commensal, ancestral house population, we performed a comparative genome scan to identify genes potentially involved with adaptation to an anthropogenic niche. We identified potential signatures of recent, positive selection in the genome of the commensal house sparrow that are absent in Bactrianus populations. The strongest selected region encompasses two major candidate genes; COL11A-which regulates craniofacial and skull development and AMY2A, part of the amylase gene family which has previously been linked to adaptation to high-starch diets in humans and dogs. Our work examines human-commensalism in an evolutionary framework, identifies genomic regions likely involved in rapid adaptation to this new niche and ties the evolution of this species to the development of modern human civilization.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

New insights into the phylogenetics and population structure of the prairie falcon (Falco mexicanus)

Background

Management requires a robust understanding of between- and within-species genetic variability, however such data are still lacking in many species. For example, although multiple population genetics studies of the peregrine falcon (Falco peregrinus) have been conducted, no similar studies have been done of the closely-related prairie falcon (F. mexicanus) and it is unclear how much genetic variation and population structure exists across the species’ range. Furthermore, the phylogenetic relationship of F. mexicanus relative to other falcon species is contested. We utilized a genomics approach (i.e., genome sequencing and assembly followed by single nucleotide polymorphism genotyping) to rapidly address these gaps in knowledge.

Results

We sequenced the genome of a single female prairie falcon and generated a 1.17 Gb (gigabases) draft genome assembly. We generated maximum likelihood phylogenetic trees using complete mitochondrial genomes as well as nuclear protein-coding genes. This process provided evidence that F. mexicanus is an outgroup to the clade that includes the peregrine falcon and members of the subgenus Hierofalco. We annotated > 16,000 genes and almost 600,000 high-quality single nucleotide polymorphisms (SNPs) in the nuclear genome, providing the raw material for a SNP assay design featuring > 140 gene-associated markers and a molecular-sexing marker. We subsequently genotyped ~ 100 individuals from California (including the San Francisco East Bay Area, Pinnacles National Park and the Mojave Desert) and Idaho (Snake River Birds of Prey National Conservation Area). We tested for population structure and found evidence that individuals sampled in California and Idaho represent a single panmictic population.

Conclusions

Our study illustrates how genomic resources can rapidly shed light on genetic variability in understudied species and resolve phylogenetic relationships. Furthermore, we found evidence of a single, randomly mating population of prairie falcons across our sampling locations. Prairie falcons are highly mobile and relatively rare long-distance dispersal events may promote gene flow throughout the range. As such, California’s prairie falcons might be managed as a single population, indicating that management actions undertaken to benefit the species at the local level have the potential to influence the species as a whole.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4615-z) contains supplementary material, which is available to authorized users.

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.

Genomics of adaptive divergence with chromosome-scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome-scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low-crossover centres relative to the high-crossover peripheries of chromosomes ("Chromosome Centre-Biased Differentiation", CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome-scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome-scale heterogeneity in crossover rate in evolutionary genomics.

Interpreting the genomic landscape of introgression

Introgression, the transfer of genetic material between species through hybridisation, occurs in many taxa and has important consequences. Genomic studies allow us to characterise the landscape of introgression across the genome, shedding light on both its adaptive benefits and the incompatibilities that help to maintain species barriers. Studies taking a genome-wide view suggest that adaptive introgression may be common, but that introgressed variation between many species is selected against throughout much of the genome. Confounding factors can complicate interpretations from these data, and computational simulations have proved vital to illustrate expected patterns under different scenarios. Future developments will move beyond correlative evidence to explicit models that account for how selection and genetic drift influence introgressed variation.

A history of hybrids? Genomic patterns of introgression in the True Geese

Background

The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression.

Results

We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow.

Conclusions

Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1048-2) contains supplementary material, which is available to authorized users.