Evidence for mito-nuclear and sex-linked reproductive barriers between the hybrid Italian sparrow and its parent species

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Conservation Mitonuclear Replacement: Facilitated mitochondrial adaptation for a changing world

Most species will not be able to migrate fast enough to cope with climate change, nor evolve quickly enough with current levels of genetic variation. Exacerbating the problem are anthropogenic influences on adaptive potential, including the prevention of gene flow through habitat fragmentation and the erosion of genetic diversity in small, bottlenecked populations. Facilitated adaptation, or assisted evolution, offers a way to augment adaptive genetic variation via artificial selection, induced hybridization, or genetic engineering. One key source of genetic variation, particularly for climatic adaptation, are the core metabolic genes encoded by the mitochondrial genome. These genes influence environmental tolerance to heat, drought, and hypoxia, but must interact intimately and co-evolve with a suite of important nuclear genes. These coadapted mitonuclear genes form some of the important reproductive barriers between species. Mitochondrial genomes can and do introgress between species in an adaptive manner, and they may co-introgress with nuclear genes important for maintaining mitonuclear compatibility. Managers should consider the relevance of mitonuclear genetic variability in conservation decision-making, including as a tool for facilitating adaptation. I propose a novel technique dubbed Conservation Mitonuclear Replacement (CmNR), which entails replacing the core metabolic machinery of a threatened species-the mitochondrial genome and key nuclear loci-with those from a closely related species or a divergent population, which may be better-adapted to climatic changes or carry a lower genetic load. The most feasible route to CmNR is to combine CRISPR-based nuclear genetic editing with mitochondrial replacement and assisted reproductive technologies. This method preserves much of an organism's phenotype and could allow populations to persist in the wild when no other suitable conservation options exist. The technique could be particularly important on mountaintops, where rising temperatures threaten an alarming number of species with almost certain extinction in the next century.

Hybridization and gene expression: Beyond differentially expressed genes

Gene expression has a key role in reproductive isolation, and studies of hybrid gene expression have identified mechanisms causing hybrid sterility. Here, we review the evidence for altered gene expression following hybridization and outline the mechanisms shown to contribute to altered gene expression in hybrids. Transgressive gene expression, transcending that of both parental species, is pervasive in early generation sterile hybrids, but also frequently observed in viable, fertile hybrids. We highlight studies showing that hybridization can result in transgressive gene expression, also in established hybrid lineages or species. Such extreme patterns of gene expression in stabilized hybrid taxa suggest that altered hybrid gene expression may result in hybridization-derived evolutionary novelty. We also conclude that while patterns of misexpression in hybrids are well documented, the understanding of the mechanisms causing misexpression is lagging. We argue that jointly assessing differences in cell composition and cell-specific changes in gene expression in hybrids, in addition to assessing changes in chromatin and methylation, will significantly advance our understanding of the basis of altered gene expression. Moreover, uncovering to what extent evolution of gene expression results in altered expression for individual genes, or entire networks of genes, will advance our understanding of how selection moulds gene expression. Finally, we argue that jointly studying the dual roles of altered hybrid gene expression, serving both as a mechanism for reproductive isolation and as a substrate for hybrid ecological adaptation, will lead to significant advances in our understanding of the evolution of gene expression.

Bayesian hybrid index and genomic cline estimation with the R package gghybrid

Admixture, the interbreeding of individuals from differentiated source populations, is now known to be a widespread phenomenon. Genomic studies of natural hybridisation can help to answer many questions on the impacts of admixture on adaptive evolution, reproductive isolation, and speciation. When a large variety of admixture proportions between two source populations exist, both geographic and genomic cline analysis are suitable methods for inferring biased, restricted or excessive gene flow at individual loci into the foreign genomic background, providing evidence for reproductive isolation, selection across an environmental transition, balancing selection, and adaptive introgression. Genomic cline analysis replaces geographic location with genome-wide hybrid index and is therefore useable in circumstances that violate geographic cline assumptions. Here, I introduce gghybrid, an R package for simple and flexible Bayesian estimation of Buerkle's hybrid index and Fitzpatrick's logit-logistic genomic clines using bi-allelic data, suitable for both small and large datasets. gghybrid allows any ploidy and uses Structure input file format. It has separate functions for hybrid index and cline estimation, treating each individual and locus respectively as an independent analysis, making it highly parallelisable. Admixture proportions from other software can alternatively be used in cline analysis, alongside parental allele frequencies. Parameters can be fixed and samples pooled for statistical model comparison with AIC or waic. Here, I describe the functions, pipeline, and statistical properties of gghybrid. Simulations reveal that model comparison with waic is preferred, and use of Bayesian posterior distributions and p values to select candidate non-null loci is problematic and should be avoided.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds

BACKGROUND

The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome.

RESULTS

We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle.

CONCLUSIONS

To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Hybridization and polyploidization generate evolutionary diversification of the genus Buddleja in the Sino-Himalayan region. A commentary on 'Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region'

This article comments on:

Fengmao Yang, Jia Ge, Yongjie Guo, Richard Olmstead and Weibang Sun. Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region, Annals of Botany, Volume 132, Issue 1, 01 July 2023, Pages 15–28, https://doi.org/10.1093/aob/mcad022

Diet and mitonuclear haplotype interactions affect growth rate in a slime mould

Trait expression in metazoans is strongly influenced by the balance of macronutrients (i.e. protein, carbohydrate and fat) in the diet. At the same time, an individual's genetic background seems to regulate the magnitude of phenotypic response to a particular diet. It needs to be better understood whether interactions between diet, genetic background and trait expression are found in unicellular eukaryotes. A protist-the slime mould, Physarum polycephalum can choose diets based on protein-to-carbohydrate (P:C) content to support optimal growth rate. Yet, the role of genetic background (variation in the mitochondrial and nuclear DNAs) in mediating growth rate response to dietary P:C ratios in the slime mould is unknown. Here, we studied the effects of interactions between mitochondrial and nuclear DNA haplotypes and diet (i.e. G × G × E interactions) on the growth rate of P. polycephalum. A genetic panel of six distinct strains of P. polycephalum that differ in their mitochondrial and nuclear DNA haplotypes was used to measure growth rate across five diets that varied in their P:C ratio and total calories. We first determined the strains' growth rate (total biomass and surface area) when grown on a set menu with access to a particular diet. We then assessed whether the growth rate of strains increased on a buffet menu with access to all diets. Our findings show that the growth rate of P. polycephalum is generally higher on diets containing more carbohydrates than protein and that total calories negatively affect the growth rate. Three-way interactions between mitochondrial, nuclear haplotypes and dietary P:C ratios affected the strains' surface area of growth but not biomass. Intriguingly, strains did not increase their surface area and biomass when they had access to all diets on the buffet menu. Our findings have broad implications for our understanding of the effect of mitonuclear interactions on trait expression across diverse eukaryotic lineages.

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-023-09612-5.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Extensive transgressive gene expression in testis but not ovary in the homoploid hybrid Italian sparrow

Hybridization can result in novel allelic combinations which can impact the hybrid phenotype through changes in gene expression. While misexpression in F₁ hybrids is well documented, how gene expression evolves in stabilized hybrid taxa remains an open question. As gene expression evolves in a stabilizing manner, break‐up of co‐evolved cis‐ and trans‐regulatory elements could lead to transgressive patterns of gene expression in hybrids. Here, we address to what extent gonad gene expression has evolved in an established and stable homoploid hybrid, the Italian sparrow (Passer italiae). Through comparison of gene expression in gonads from individuals of the two parental species (i.e., house and Spanish sparrow) to that of Italian sparrows, we find evidence for strongly transgressive expression in male Italian sparrows—2530 genes (22% of testis genes tested for inheritance) exhibit expression patterns outside the range of both parent species. In contrast, Italian sparrow ovary expression was similar to that of one of the parent species, the house sparrow (Passer domesticus). Moreover, the Italian sparrow testis transcriptome is 26 times as diverged from those of the parent species as the parental transcriptomes are from each other, despite being genetically intermediate. This highlights the potential for regulation of gene expression to produce novel variation following hybridization. Genes involved in mitochondrial respiratory chain complexes and protein synthesis are enriched in the subset that is over‐dominantly expressed in Italian sparrow testis, suggesting that selection on key functions has moulded the hybrid Italian sparrow transcriptome.

Genome-wide local ancestry and evidence for mitonuclear coadaptation in African hybrid cattle populations (Bos taurus/indicus)

The phenotypic diversity of African cattle reflects adaptation to a wide range of agroecological conditions, human-mediated selection preferences, and complex patterns of admixture between the humpless Bos taurus (taurine) and humped Bos indicus (zebu) subspecies, which diverged 150-500 thousand years ago. Despite extensive admixture, all African cattle possess taurine mitochondrial haplotypes, even populations with significant zebu biparental and male uniparental nuclear ancestry. This has been interpreted as the result of human-mediated dispersal ultimately stemming from zebu bulls imported from South Asia during the last three millennia. Here, we assess whether ancestry at mitochondrially targeted nuclear genes in African admixed cattle is impacted by mitonuclear functional interactions. Using high-density SNP data, we find evidence for mitonuclear coevolution across hybrid African cattle populations with a significant increase of taurine ancestry at mitochondrially targeted nuclear genes. Our results, therefore, support the hypothesis of incompatibility between the taurine mitochondrial genome and the zebu nuclear genome.

Parent-of-Origin Effects on Seed Size Modify Heterosis Responses in Arabidopsis thaliana

Parent-of-origin effects arise when a phenotype depends on whether it is inherited maternally or paternally. Parent-of-origin effects can exert a strong influence on F1 seed size in flowering plants, an important agronomic and life-history trait that can contribute to biomass heterosis. Here we investigate the natural variation in the relative contributions of the maternal and paternal genomes to F1 seed size across 71 reciprocal pairs of F1 hybrid diploids and the parental effect on F1 seed size heterosis. We demonstrate that the paternally derived genome influences F1 seed size more significantly than previously appreciated. We further demonstrate (by disruption of parental genome dosage balance in F1 triploid seeds) that hybridity acts as an enhancer of genome dosage effects on F1 seed size, beyond that observed from hybridity or genome dosage effects on their own. Our findings indicate that interactions between genetic hybridity and parental genome dosage can enhance heterosis effects in plants, opening new avenues for boosting heterosis breeding in crop plants.

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.

Genomes of hybrid lineages are mosaics of those of their parent species and harbour variation that has the potential to facilitate adaptation when hybrids encounter diverse environments. However, genetic incompatibilities between parental species can also act to limit possible combinations of parental alleles, constraining hybrid genome formation. What is the relative importance of selection and constraints in form of admixture proportions and genomic architecture in this process? We investigated this in the Italian sparrow, a hybrid species resulting from past hybridization between the house and Spanish sparrow. Using three independent hybrid lineages, we addressed how their genomes, harbouring different parental combinations, have evolved. We examined the roles of selection due to divergent local adaptation, recombination and purging of genetic incompatibilities in predicting differentiation. We found that selection against incompatibilities may constrain hybrid genome composition. In addition, signals of local selection as well as estimates of differentiation were correlated across populations, and outliers were shared among the hybrid lineages more often than expected by chance. Overall, our results suggest that in the Italian sparrow selection interacts with constraints linked to genetic incompatibilities affecting which sections of the genome can readily diverge among hybrid lineages.

The role of mitonuclear incompatibilities in allopatric speciation

Aerobic metabolism in eukaryotic cells requires extensive interactions between products of the nuclear and mitochondrial genomes. Rapid evolution of the mitochondrial genome, including fixation of both adaptive and deleterious mutations, creates intrinsic selection pressures favoring nuclear gene mutations that maintain mitochondrial function. As this process occurs independently in allopatry, the resulting divergence between conspecific populations can subsequently be manifest in mitonuclear incompatibilities in inter-population hybrids. Such incompatibilities, mitonuclear versions of Bateson-Dobzhansky-Muller incompatibilities that form the standard model for allopatric speciation, can potentially restrict gene flow between populations, ultimately resulting in varying degrees of reproductive isolation. The potential role of mitonuclear incompatibilities in speciation is further enhanced where mtDNA substitution rates are elevated compared to the nuclear genome and where population structure maintains allopatry for adequate time to evolve multiple mitonuclear incompatibilities. However, the fact that mitochondrial introgression occurs across species boundaries has raised questions regarding the efficacy of mitonuclear incompatibilities in reducing gene flow. Several scenarios now appear to satisfactorily explain this phenomenon, including cases where differences in mtDNA genetic load may drive introgression or where co-introgression of coadapted nuclear genes may support the function of introgressed mtDNA. Although asymmetries in reproductive isolation between taxa are consistent with mitonuclear incompatibilities, interactions between autosomes and sex chromosomes yield similar predictions that are difficult to disentangle. With regard to establishing reproductive isolation while in allopatry, existing studies clearly suggest that mitonuclear incompatibilities can contribute to the evolution of barriers to gene flow. However, there is to date relatively little definitive evidence supporting a primary role for mitonuclear incompatibilities in the speciation process.

Predictability and parallelism in the contemporary evolution of hybrid genomes

Hybridization between species is widespread across the tree of life. As a result, many species, including our own, harbor regions of their genome derived from hybridization. Despite the recognition that this process is widespread, we understand little about how the genome stabilizes following hybridization, and whether the mechanisms driving this stabilization tend to be shared across species. Here, we dissect the drivers of variation in local ancestry across the genome in replicated hybridization events between two species pairs of swordtail fish: Xiphophorus birchmanni × X. cortezi and X. birchmanni × X. malinche. We find unexpectedly high levels of repeatability in local ancestry across the two types of hybrid populations. This repeatability is attributable in part to the fact that the recombination landscape and locations of functionally important elements play a major role in driving variation in local ancestry in both types of hybrid populations. Beyond these broad scale patterns, we identify dozens of regions of the genome where minor parent ancestry is unusually low or high across species pairs. Analysis of these regions points to shared sites under selection across species pairs, and in some cases, shared mechanisms of selection. We show that one such region is a previously unknown hybrid incompatibility that is shared across X. birchmanni × X. cortezi and X. birchmanni × X. malinche hybrid populations.

Mitochondrial Transfer in Cardiovascular Disease: From Mechanisms to Therapeutic Implications

Mitochondrial dysfunction has been proven to play a critical role in the pathogenesis of cardiovascular diseases. The phenomenon of intercellular mitochondrial transfer has been discovered in the cardiovascular system. Studies have shown that cell-to-cell mitochondrial transfer plays an essential role in regulating cardiovascular system development and maintaining normal tissue homeostasis under physiological conditions. In pathological conditions, damaged cells transfer dysfunctional mitochondria toward recipient cells to ask for help and take up exogenous functional mitochondria to alleviate injury. In this review, we summarized the mechanism of mitochondrial transfer in the cardiovascular system and outlined the fate and functional role of donor mitochondria. We also discussed the advantage and challenges of mitochondrial transfer strategies, including cell-based mitochondrial transplantation, extracellular vesicle-based mitochondrial transplantation, and naked mitochondrial transplantation, for the treatment of cardiovascular disorders. We hope this review will provide perspectives on mitochondrial-targeted therapeutics in cardiovascular diseases.

The avian W chromosome is a refugium for endogenous retroviruses with likely effects on female-biased mutational load and genetic incompatibilities

It is a broadly observed pattern that the non-recombining regions of sex-limited chromosomes (Y and W) accumulate more repeats than the rest of the genome, even in species like birds with a low genome-wide repeat content. Here, we show that in birds with highly heteromorphic sex chromosomes, the W chromosome has a transposable element (TE) density of greater than 55% compared to the genome-wide density of less than 10%, and contains over half of all full-length (thus potentially active) endogenous retroviruses (ERVs) of the entire genome. Using RNA-seq and protein mass spectrometry data, we were able to detect signatures of female-specific ERV expression. We hypothesize that the avian W chromosome acts as a refugium for active ERVs, probably leading to female-biased mutational load that may influence female physiology similar to the 'toxic-Y' effect in Drosophila males. Furthermore, Haldane's rule predicts that the heterogametic sex has reduced fertility in hybrids. We propose that the excess of W-linked active ERVs over the rest of the genome may be an additional explanatory variable for Haldane's rule, with consequences for genetic incompatibilities between species through TE/repressor mismatches in hybrids. Together, our results suggest that the sequence content of female-specific W chromosomes can have effects far beyond sex determination and gene dosage. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

The rates of introgression and barriers to genetic exchange between hybridizing species: sex chromosomes vs autosomes

Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.

Mitonuclear conflict in a macaque species exhibiting phylogenomic discordance

Speciation and hybridization are intertwined processes in the study of evolution. Hybridization between sufficiently diverged populations can result in genomic conflict within offspring, causing reduced viability and fertility, thus increasing divergence between populations. Conflicts between mitochondrial and nuclear genes are increasingly found to play a role in this process in various systems. We examine the possibility of this conflict in the bear macaque, Macaca arctoides (Primates: Cercopithecidae), a primate species exhibiting mitonuclear discordance due to extensive hybridization with species in the sinica and fascicularis groups. Here, divergence, introgression and natural selection of mitonuclear genes (N = 160) relative to nuclear control genes (N = 144) were analysed to determine whether there are evolutionary processes involved in resolving the potential conflict caused by mitonuclear discordance. Nucleotide divergence of mitonuclear genes is increased relative to control nuclear genes between M. arctoides and the species sharing its nuclear ancestry (p = 0.007), consistent with genetic conflict. However, measures of introgression and selection do not identify large-scale co-introgression or co-evolution as means to resolve mitonuclear conflict. Nonetheless, mitochondrial tRNA synthetases stand out in analyses using dN/dS and extended branch lengths as potential targets of selection. The methodology implemented provides a framework that can be used to examine the effects of mitonuclear co-introgression and co-evolution on a genomic scale in a variety of systems.

The genomic consequences of hybridization

In the past decade, advances in genome sequencing have allowed researchers to uncover the history of hybridization in diverse groups of species, including our own. Although the field has made impressive progress in documenting the extent of natural hybridization, both historical and recent, there are still many unanswered questions about its genetic and evolutionary consequences. Recent work has suggested that the outcomes of hybridization in the genome may be in part predictable, but many open questions about the nature of selection on hybrids and the biological variables that shape such selection have hampered progress in this area. We synthesize what is known about the mechanisms that drive changes in ancestry in the genome after hybridization, highlight major unresolved questions, and discuss their implications for the predictability of genome evolution after hybridization.

Common patterns in the molecular phylogeography of western palearctic birds: a comprehensive review

A plethora of studies have offered crucial insights in the phylogeographic status of Western Palearctic bird species. However, an overview integrating all this information and analyzing the combined results is still missing. In this study, we compiled all published peer-reviewed and grey literature available on the phylogeography of Western Palearctic bird species. Our literature review indicates a total number of 198 studies, with the overwhelming majority published as journal articles (n = 186). In total, these literature items offer information on 145 bird species. 85 of these species are characterized by low genetic differentiation, 46 species indicate genetic variation but no geographic structuring i.e. panmixia, while 14 species show geographically distinct lineages and haplotypes. Majority of bird species inhabiting the Western Palearctic display genetic admixture. The glaciation cycles in the past few million years were pivotal factors in shaping this situation: during warm periods many species expanded their distribution range to the north over wide areas of Eurasia; whereas, during ice ages most areas were no longer suitable and species retreated to refugia, where lineages mixed.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s10336-021-01893-x.

Locus-specific introgression in young hybrid swarms: Drift may dominate selection

Closely related species that have previously inhabited geographically separated ranges are hybridizing at an increasing rate due to human disruptions. These human-mediated hybrid zones can be used to study reproductive isolation between species at secondary contact, including examining locus-specific rates of introgression. Introgression is expected to be heterogenous across the genome, reflecting variation in selection. Those loci that introgress especially slowly are good candidates for being involved in reproductive isolation, while those loci that introgress quickly may be involved in adaptive introgression. In the context of conservation, policy makers are especially concerned about introduced alleles moving quickly into the background of a native or endemic species, as these alleles could replace the native alleles in the population, leading to extinction via hybridization. We applied genomic cline analyses to 44,997 SNPs to identify loci introgressing more or less when compared to the genome wide expectation in a human-mediated hybridizing population of red deer and sika in Kintyre Scotland. We found 11.4% of SNPs had cline centres that were significantly different from the genome wide expectation, and 17.6% of all SNPs had excess rates of introgression. Based on simulations, we believe that many of these markers have diverged from the genome-wide average due to drift, rather than because of selection, and we suggest that these simulations can be useful as a null distribution for future studies of genomic clines. Future work on red deer and sika could determine the policy implications of allelic-replacement due to drift rather than selection, and could use replicate, geographically distinct hybrid zones to narrow down those loci that are responding to selection.

Intraspecific genomic variation and local adaptation in a young hybrid species

Hybridization increases genetic variation, hence hybrid species may have greater evolutionary potential once their admixed genomes have stabilized and incompatibilities have been purged. Yet, little is known about how such hybrid lineages evolve at the genomic level following their formation, in particular their adaptive potential. Here we investigate how the Italian sparrow (Passer italiae), a homoploid hybrid species, has evolved and locally adapted to its variable environment. Using restriction site-associated DNA sequencing (RAD-seq) on several populations across the Italian peninsula, we evaluate how genomic constraints and novel genetic variation have influenced population divergence and adaptation. We show that population divergence within this hybrid species has evolved in response to climatic variation, suggesting ongoing local adaptation. As found previously in other nonhybrid species, climatic differences appear to increase population differentiation. We also report strong population divergence in a gene known to affect beak morphology. Most of the strongly divergent loci among Italian sparrow populations do not seem to be differentiated between its parent species, the house and Spanish sparrows. Unlike in the hybrid, population divergence within each of the parental taxa has occurred mostly at loci with high allele frequency difference between the parental species, suggesting that novel combinations of parental alleles in the hybrid have not necessarily enhanced its evolutionary potential. Rather, our study suggests that constraints linked to incompatibilities may have restricted the evolution of this admixed genome, both during and after hybrid species formation.

Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population

Hybridization between divergent lineages generates new allelic combinations. One mechanism that can hinder the formation of hybrid populations is mitonuclear incompatibility, that is, dysfunctional interactions between proteins encoded in the nuclear and mitochondrial genomes (mitogenomes) of diverged lineages. Theoretically, selective pressure due to mitonuclear incompatibility can affect genotypes in a hybrid population in which nuclear genomes and mitogenomes from divergent lineages admix. To directly and thoroughly observe this key process, we de novo sequenced the 747-Mb genome of the coastal goby, Chaenogobius annularis, and investigated its integrative genomic phylogeographics using RNA-sequencing, RAD-sequencing, genome resequencing, whole mitogenome sequencing, amplicon sequencing, and small RNA-sequencing. Chaenogobius annularis populations have been geographically separated into Pacific Ocean (PO) and Sea of Japan (SJ) lineages by past isolation events around the Japanese archipelago. Despite the divergence history and potential mitonuclear incompatibility between these lineages, the mitogenomes of the PO and SJ lineages have coexisted for generations in a hybrid population on the Sanriku Coast. Our analyses revealed accumulation of nonsynonymous substitutions in the PO-lineage mitogenomes, including two convergent substitutions, as well as signals of mitochondrial lineage-specific selection on mitochondria-related nuclear genes. Finally, our data implied that a microRNA gene was involved in resolving mitonuclear incompatibility. Our integrative genomic phylogeographic approach revealed that mitonuclear incompatibility can affect genome evolution in a natural hybrid population.

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.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

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.

Understanding Admixture: Haplodiploidy to the Rescue

Hybridization has broad evolutionary consequences, from fueling or counteracting speciation to facilitating adaptation to novel environments. Hybridization and subsequent introgression appear widespread along the tree of life. However, our understanding of how distinct evolutionary forces shape admixed genomes and the fate of introgressed genetic variants remains scarce. Most admixture research in animals has focused on diploid organisms. We propose that haplodiploid organisms can help resolve open questions about the genomic consequences of hybridization in natural populations. The ploidy difference between haploid males and diploid females, the availability of genome-wide male haplotypes, and ongoing cases of admixture make haplodiploid organisms promising models to improve our knowledge with regards to the evolution of hybrid genomes.

Evidence for adaptive introgression of exons across a hybrid swarm in deer

Background

Secondary contact between closely related lineages can result in a variety of outcomes, including hybridization, depending upon the strength of reproductive barriers. By examining the extent to which different parts of the genome introgress, it is possible to infer the strength of selection and gain insight into the evolutionary trajectory of lineages. Following secondary contact approximately 8000 years ago in the Pacific Northwest, mule deer (Odocoileus hemionus hemionus) and black-tailed deer (O. h. columbianus) formed a hybrid swarm along the Cascade mountain range despite substantial differences in body size (up to two times) and habitat preference. In this study, we examined genetic population structure, extent of introgression, and selection pressures in freely interbreeding populations of mule deer and black-tailed deer using mitochondrial DNA sequences, 9 microsatellite loci, and 95 SNPs from protein-coding genes.

Results

We observed bi-directional hybridization and classified approximately one third of the 172 individuals as hybrids, almost all of which were beyond the F1 generation. High genetic differentiation between black-tailed deer and mule deer at protein-coding genes suggests that there is positive divergent selection, though selection on these loci is relatively weak. Contrary to predictions, there was not greater selection on protein-coding genes thought to be associated with immune function and mate choice. Geographic cline analyses were consistent across genetic markers, suggesting long-term stability (over hundreds of generations), and indicated that the center of the hybrid swarm is 20-30 km to the east of the Cascades ridgeline, where there is a steep ecological transition from wet, forested habitat to dry, scrub habitat.

Conclusions

Our data are consistent with a genetic boundary between mule deer and black-tailed deer that is porous but maintained by many loci under weak selection having a substantial cumulative effect. The absence of clear reproductive barriers and the consistent centering of geographic clines at a sharp ecotone suggests that ecology is a driver of hybrid swarm dynamics. Adaptive introgression in this study (and others) promotes gene flow and provides valuable insight into selection strength on specific genes and the evolutionary trajectory of hybridizing taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1497-x) contains supplementary material, which is available to authorized users.

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, sex-specific genomic architecture and local adaptation

While gene flow can reduce the potential for local adaptation, hybridization may conversely provide genetic variation that increases the potential for local adaptation. Hybridization may also affect adaptation through altering sexual dimorphism and sexual conflict, but this remains largely unstudied. Here, we discuss how hybridization may affect sexual dimorphism and conflict due to differential effects of hybridization on males and females, and then how this, in turn, may affect local adaptation. First, in species with heterochromatic sexes, the lower viability of the heterogametic sex in hybrids could shift the balance in sexual conflict. Second, sex-specific inheritance of the mitochondrial genome in hybrids may lead to cytonuclear mismatches, for example, in the form of 'mother's curse', with potential consequences for sex ratio and sex-specific expression. Third, sex-biased introgression and recombination may lead to sex-specific consequences of hybridization. Fourth, transgressive segregation of sexually antagonistic alleles could increase sexual dimorphism in hybrid populations. Sexual dimorphism can reduce sexual conflict and enhance intersexual niche partitioning, increasing the fitness of hybrids. Adaptive introgression of alleles reducing sexual conflict or enhancing intersexual niche partitioning may facilitate local adaptation, and could favour the colonization of novel habitats. We review these consequences of hybridization on sex differences and local adaptation, and discuss how their prevalence and importance could be tested empirically.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Genomic and phenotypic consequences of two independent secondary contact zones between allopatric lineages of the anadromous ice goby Leucopsarion petersii

Genetic and phenotypic analyses of independent secondary contact zones between certain pairs of divergent populations offer powerful opportunities to assess whether the consequences vary with different environmental backgrounds. Populations of the ice goby Leucopsarion petersii are distributed throughout the Japanese archipelago and comprise genetically and phenotypically divergent groups in the Japan Sea and the Pacific Ocean. In particular, populations in the Japan Sea have a larger body size and numbers of vertebrae than those in the Pacific Ocean. Herein, we performed integrated analyses of genotypes and phenotypes of two independent secondary contact zones and investigated their consequences. Population genetic analyses revealed asymmetric introgression of the mitochondrial genome of either lineage relative to little admixture of nuclear genomes in both secondary contact zones. On phenotype analyses, vertebral numbers were clearly explained by nuclear genomic ancestry in both secondary contact zones, whereas body size was not, suggesting that a little introgression of nuclear genes regulates body size. Actually, we observed biased introgression of a candidate gene, neuropeptide Y (NPY), which potentially controls body size in the ice goby. Moreover, the body size changes in the introgressed populations possibly affect the introgression patterns of mitochondrial genomes across these zones. Collectively, our results demonstrated that genomic and phenotypic consequences of secondary contact varied in marine variable environments.

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.

Allele-specific expression variation at different ploidy levels in Squalius alburnoides

Allopolyploid plants are long known to be subject to a homoeolog expression bias of varying degree. The same phenomenon was only much later suspected to occur also in animals based on studies of single selected genes in an allopolyploid vertebrate, the Iberian fish Squalius alburnoides. Consequently, this species became a good model for understanding the evolution of gene expression regulation in polyploid vertebrates. Here, we analyzed for the first time genome-wide allele-specific expression data from diploid and triploid hybrids of S. alburnoides and compared homoeolog expression profiles of adult livers and of juveniles. Co-expression of alleles from both parental genomic types was observed for the majority of genes, but with marked homoeolog expression bias, suggesting homoeolog specific reshaping of expression level patterns in hybrids. Complete silencing of one allele was also observed irrespective of ploidy level, but not transcriptome wide as previously speculated. Instead, it was found only in a restricted number of genes, particularly ones with functions related to mitochondria and ribosomes. This leads us to hypothesize that allelic silencing may be a way to overcome intergenomic gene expression interaction conflicts, and that homoeolog expression bias may be an important mechanism in the achievement of sustainable genomic interactions, mandatory to the success of allopolyploid systems, as in S. alburnoides.

Analysis of 11 trace elements in flight feathers of Italian Sparrows in southern Italy: A study of bioaccumulation through age classes, variability in three years of sampling, and relations with body condition

Trace elements have been acknowledged as one of the subtlest environmental hazards in all compartments of the total environment. Enhanced by activities in the anthroposphere, they accumulate in the atmosphere, lithosphere, and hydrosphere. Eventually, trace elements can bioaccumulate or biomagnifiy in the biosphere, with harmful effects on animals occupying higher trophic levels, including humans. Accordingly, there is great interest in assessing and monitoring trace element concentrations in the biosphere, and birds, especially passerines, have been commonly chosen as biomonitors. In this study, the concentration of 11 trace elements was measured (i.e. aluminum, chromium, manganese, iron, nickel, copper, zinc, arsenic, cadmium, barium, and lead) in flight feathers of Italian Sparrows, a common bird species hitherto not analysed in this respect. Samples were collected in an agricultural area in southern Italy, where a mosaic of natural environments, urbanized areas and industrial facilities can be found. Linear mixed modelling was applied to the analysis of flight feathers in juveniles, juvenile birds moulting to adulthood, and adults in three sampling years on 184 birds. Results are timely as they add new data to the scarce available information on Ba and As in bird feathers and showed clear bioaccumulation patterns from juveniles to adults for As, Cr, and Cd. Moreover, the modelling approach showed that the concentration of elements such as As, Cd and Cr can be variable across the years and that some elements, notably Cd and Ba, were inversely correlated with body mass and wing length, respectively, suggesting potential negative effects on bird health. Finally, when modelling bird body condition and trace elements, results showed that Cd and Ba negatively affect birds regardless of age or sampling year. Thus, the Italian Sparrow could be considered as a valuable biomonitor for trace elements in the total environment, especially for Cd and Ba.

Morphologically cryptic Amazonian bird species pairs exhibit strong postzygotic reproductive isolation

We possess limited understanding of how speciation unfolds in the most species-rich region of the planet-the Amazon basin. Hybrid zones provide valuable information on the evolution of reproductive isolation, but few studies of Amazonian vertebrate hybrid zones have rigorously examined the genome-wide underpinnings of reproductive isolation. We used genome-wide genetic datasets to show that two deeply diverged, but morphologically cryptic sister species of forest understorey birds show little evidence for prezygotic reproductive isolation, but substantial postzygotic isolation. Patterns of heterozygosity and hybrid index revealed that hybrid classes with heavily recombined genomes are rare and closely match simulations with high levels of selection against hybrids. Genomic and geographical clines exhibit a remarkable similarity across loci in cline centres, and have exceptionally narrow cline widths, suggesting that postzygotic isolation is driven by genetic incompatibilities at many loci, rather than a few loci of strong effect. We propose Amazonian understorey forest birds speciate slowly via gradual accumulation of postzygotic genetic incompatibilities, with prezygotic barriers playing a less important role. Our results suggest old, cryptic Amazonian taxa classified as subspecies could have substantial postzygotic isolation deserving species recognition and that species richness is likely to be substantially underestimated in Amazonia.

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.

Differential sperm-egg interactions in experimental pairings between two subspecies and their hybrids in a passerine bird

Speciation research has largely overlooked reproductive barriers acting between copulation and the formation of the zygote (i.e., postmating, prezygotic [PMPZ] barriers), especially in internally fertilizing vertebrates. Nonetheless, it is becoming clear that PMPZ reproductive barriers can play a role in the formation and maintenance of species boundaries. We investigated sperm-egg interactions in the recently diverged subspecies pairs of the long-tailed finch, Poephila acuticauda acuticauda and P. a. hecki, to explore potential PMPZ barriers. Specifically, we compared the number of sperm reaching the perivitelline layer (PVL) of the ova, and hence the site of fertilization, in both intra- and inter-subspecies pairings and pairings of F1 hybrid adults with one parental subspecies. Although we found no difference in PVL sperm number among intra- and inter-subspecific pairs, a significantly lower number of sperm reached the site of fertilization in a backcross pairing. As low numbers of PVL sperm appear to be associated with low fertilization success in birds, our findings offer insight into the potential role of postcopulatory processes in limiting gene flow between the subspecies and may help explain the relatively narrow hybrid zone that exists in the wild in this species. Though further work is needed to gain a comprehensive understanding of the morphological, physiological, and molecular mechanisms underlying our results, our study supports the role of PMPZ reproductive barriers in avian speciation, even in recently diverged taxa, that may not yet be fully genetically incompatible.

The Biogeography of Deep Time Phylogenetic Reticulation

Most phylogenies are typically represented as purely bifurcating. However, as genomic data have become more common in phylogenetic studies, it is not unusual to find reticulation among terminal lineages or among internal nodes (deep time reticulation; DTR). In these situations, gene flow must have happened in the same or adjacent geographic areas for these DTRs to have occurred and therefore biogeographic reconstruction should provide similar area estimates for parental nodes, provided extinction or dispersal has not eroded these patterns. We examine the phylogeny of the widely distributed New World kingsnakes (Lampropeltis), determine if DTR is present in this group, and estimate the ancestral area for reticulation. Importantly, we develop a new method that uses coalescent simulations in a machine learning framework to show conclusively that this phylogeny is best represented as reticulating at deeper time. Using joint probabilities of ancestral area reconstructions on the bifurcating parental lineages from the reticulating node, we show that this reticulation likely occurred in northwestern Mexico/southwestern US, and subsequently, led to the diversification of the Mexican kingsnakes. This region has been previously identified as an area important for understanding speciation and secondary contact with gene flow in snakes and other squamates. This research shows that phylogenetic reticulation is common, even in well-studied groups, and that the geographic scope of ancient hybridization is recoverable.

Postcopulatory sexual selection reduces Z-linked genetic variation and might contribute to the large Z effect in passerine birds

The X and Z sex chromosomes play a disproportionately large role in intrinsic postzygotic isolation. The underlying mechanisms of this large X/Z effect are, however, still poorly understood. Here we tested whether faster rates of molecular evolution caused by more intense positive selection or genetic drift on the Z chromosome could contribute to the large Z effect in two closely related passerine birds, the Common Nightingale (Luscinia megarhynchos) and the Thrush Nightingale (L. luscinia). We found that the two species differ in patterns of molecular evolution on the Z chromosome. The Z chromosome of L. megarhynchos showed lower levels of within-species polymorphism and an excess of non-synonymous polymorphisms relative to non-synonymous substitutions. This is consistent with increased levels of genetic drift on this chromosome and may be attributed to more intense postcopulatory sexual selection acting on L. megarhynchos males as was indicated by significantly longer sperm and higher between-male variation in sperm length in L. megarhynchos compared to L. luscinia. Interestingly, analysis of interspecific gene flow on the Z chromosome revealed relatively lower levels of introgression from L. megarhynchos to L. luscinia than vice versa, indicating that the Z chromosome of L. megarhynchos accumulated more hybrid incompatibilities. Our results are consistent with the view that postcopulatory sexual selection may reduce the effective population size of the Z chromosome and thus lead to stronger genetic drift on this chromosome in birds. This can result in relatively faster accumulation of hybrid incompatibilities on the Z and thus contribute to the large Z effect.

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 Mitonuclear Dimension of Neanderthal and Denisovan Ancestry in Modern Human Genomes

Some human populations interbred with Neanderthals and Denisovans, resulting in substantial contributions to modern-human genomes. Therefore, it is now possible to use genomic data to investigate mechanisms that shaped historical gene flow between humans and our closest hominin relatives. More generally, in eukaryotes, mitonuclear interactions have been argued to play a disproportionate role in generating reproductive isolation. There is no evidence of mtDNA introgression into modern human populations, which means that all introgressed nuclear alleles from archaic hominins must function on a modern-human mitochondrial background. Therefore, mitonuclear interactions are also potentially relevant to hominin evolution. We performed a detailed accounting of mtDNA divergence among hominin lineages and used population-genomic data to test the hypothesis that mitonuclear incompatibilities have preferentially restricted the introgression of nuclear genes with mitochondrial functions. We found a small but significant underrepresentation of introgressed Neanderthal alleles at such nuclear loci. Structural analyses of mitochondrial enzyme complexes revealed that these effects are unlikely to be mediated by physically interacting sites in mitochondrial and nuclear gene products. We did not detect any underrepresentation of introgressed Denisovan alleles at mitochondrial-targeted loci, but this may reflect reduced power because locus-specific estimates of Denisovan introgression are more conservative. Overall, we conclude that genes involved in mitochondrial function may have been subject to distinct selection pressures during the history of introgression from archaic hominins but that mitonuclear incompatibilities have had, at most, a small role in shaping genome-wide introgression patterns, perhaps because of limited functional divergence in mtDNA and interacting nuclear genes.

Hybrid chickadees are deficient in learning and memory

Identifying the phenotypes underlying postzygotic reproductive isolation is crucial for fully understanding the evolution and maintenance of species. One potential postzygotic isolating barrier that has rarely been examined is learning and memory ability in hybrids. Learning and memory are important fitness-related traits, especially in scatter-hoarding species, where accurate retrieval of hoarded food is vital for winter survival. Here, we test the hypothesis that learning and memory ability can act as a postzygotic isolating barrier by comparing these traits among two scatter-hoarding songbird species, black-capped (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis), and their naturally occurring hybrids. In an outdoor aviary setting, we find that hybrid chickadees perform significantly worse on an associative learning spatial task and are worse at solving a novel problem compared to both parental species. Deficiencies in learning and memory abilities could therefore contribute to postzygotic reproductive isolation between chickadee species. Given the importance of learning and memory for fitness, our results suggest that these traits may play an important, but as yet overlooked, role in postzygotic reproductive isolation.

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.