Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

Past human expansions shaped the spatial pattern of Neanderthal ancestry

The worldwide expansion of modern humans (Homo sapiens) started before the extinction of Neanderthals (Homo neanderthalensis). Both species coexisted and interbred, leading to slightly higher introgression in East Asians than in Europeans. This distinct ancestry level has been argued to result from selection, but range expansions of modern humans could provide an alternative explanation. This hypothesis would lead to spatial introgression gradients, increasing with distance from the expansion source. We investigate the presence of Neanderthal introgression gradients after past human expansions by analyzing Eurasian paleogenomes. We show that the out-of-Africa expansion resulted in spatial gradients of Neanderthal ancestry that persisted through time. While keeping the same gradient orientation, the expansion of early Neolithic farmers contributed decisively to reducing the Neanderthal introgression in European populations compared to Asian populations. This is because Neolithic farmers carried less Neanderthal DNA than preceding Paleolithic hunter-gatherers. This study shows that inferences about past human population dynamics can be made from the spatiotemporal variation in archaic introgression.

Multiple introgressions shape mitochondrial evolutionary history in Drosophila paulistorum and the Drosophila willistoni group

Hybridization and the consequent introgression of genomic elements is an important source of genetic diversity for biological lineages. This is particularly evident in young clades in which hybrid incompatibilities are still incomplete and mixing between species is more likely to occur. Drosophila paulistorum, a representative of the Neotropical Drosophila willistoni subgroup, is a classic model of incipient speciation. The species is divided into six semispecies that show varying degrees of pre- and post-mating incompatibility with each other. In the present study, we investigate the mitochondrial evolutionary history of D. paulistorum and the willistoni subgroup. For that, we perform phylogenetic and comparative analyses of the complete mitochondrial genomes and draft nuclear assemblies of 25 Drosophila lines of the willistoni and saltans species groups. Our results show that the mitochondria of D. paulistorum are polyphyletic and form two non-sister clades that we name α and β. Identification and analyses of nuclear mitochondrial insertions further reveal that the willistoni subgroup has an α-like mitochondrial ancestor and strongly suggest that both the α and β mitochondria of D. paulistorum were acquired through introgression from unknown fly lineages of the willistoni subgroup. We also uncover multiple mitochondrial introgressions across D. paulistorum semispecies and generate novel insight into the evolution of the species.

Homoploid hybrids are common but evolutionary dead ends, whereas polyploidy is not linked to hybridization in a group of Pyrenean saxifrages

Hybridization and polyploidy are major forces in plant evolution. Homoploid hybridization can generate new species via hybrid speciation, or modify extant evolutionary lineages through introgression. Polyploidy enables instantaneous reproductive isolation from the parental lineage(s) and is often coupled with evolutionary innovations, especially when linked to hybridization. While allopolyploidy is a well-known and common mechanism of plant speciation, the evolutionary role of autopolyploidy might have been underestimated. Here, we studied the saxifrages of Saxifraga subsection Saxifraga in the Pyrenees, which easily hybridise and include polyploid populations of uncertain origin, as a model to unravel evolutionary consequences and origin of hybridization and polyploidy. Additionally, we investigate the phylogenetic relationship between the two subspecies of the endemic S. pubescens to ascertain whether they should rather be treated as different species. For these purposes, we combined ploidy-informed restriction associated DNA analyses, plastid DNA sequences and morphological data on a comprehensive population sample of seven species. Our results unravel multiple homoploid hybridization events at the diploid level between different species pairs, but with limited evolutionary impact. The ploidy-informed analyses reveal that all tetraploid populations detected in the present study belong to the widespread alpine species S. moschata. Although of autopolyploid origin, they are to some extent morphologically differentiated and underwent a different evolutionary pathway than their diploid parent. However, the high plastid DNA diversity and the internal structure within eastern and western population groups suggest multiple origins of the polyploids. Finally, our phylogenetic analyses show that S. pubescens and S. iratiana are clearly not sister lineages, and should consequently be considered as independent species.

Midfacial Morphology and Neandertal-Modern Human Interbreeding

Ancient DNA from, Neandertal and modern human fossils, and comparative morphological analyses of them, reveal a complex history of interbreeding between these lineages and the introgression of Neandertal genes into modern human genomes. Despite substantial increases in our knowledge of these events, the timing and geographic location of hybridization events remain unclear. Six measures of facial size and shape, from regional samples of Neandertals and early modern humans, were used in a multivariate exploratory analysis to try to identify regions in which early modern human facial morphology was more similar to that of Neandertals, which might thus represent regions of greater introgression of Neandertal genes. The results of canonical variates analysis and hierarchical cluster analysis suggest important affinities in facial morphology between both Middle and Upper Paleolithic early modern humans of the Near East with Neandertals, highlighting the importance of this region for interbreeding between the two lineages.

Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations

Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.

Assessing conservation status with extensive but low-resolution data: Application of frequentist and Bayesian models to endangered Athabasca River rainbow trout

Use of extensive but low-resolution abundance data is common in the assessment of species at-risk status based on quantitative decline criteria under International Union for Conservation of Nature (IUCN) and national endangered species legislation. Such data can be problematic for 3 reasons. First, statistical power to reject the null hypothesis of no change is often low because of small sample size and high sampling uncertainty leading to a high frequency of type II errors. Second, range-wide assessments composed of multiple site-specific observations do not effectively weight site-specific trends into global trends. Third, uncertainty in site-specific temporal trends and relative abundance are not propagated at the appropriate spatial scale. A common result is the propensity to underestimate the magnitude of declines and therefore fail to identify the appropriate at-risk status for a species. We used 3 statistical approaches, from simple to more complex, to estimate temporal decline rates for a designatable unit (DU) of rainbow trout in the Athabasca River watershed in western Canada. This DU is considered a native species for purposes of listing because of its genetic composition characterized as >0.95 indigenous origin in the face of continuing introgressive hybridization with introduced populations in the watershed. Analysis of abundance trends from 57 time series with a fixed-effects model identified 33 sites with negative trends, but only 2 were statistically significant. By contrast, a hierarchical linear mixed model weighted by site-specific abundance provided a DU-wide decline estimate of 16.4% per year and a 3-generation decline of 93.2%. A hierarchical Bayesian mixed model yielded a similar 3-generation decline trend of 91.3% and the posterior distribution showed that the estimate had a >99% probability of exceeding thresholds for an endangered listing. We conclude that the Bayesian approach was the most useful because it provided a probabilistic statement of threshold exceedance in support of an at-risk status recommendation.

Genetic traces of hybrid zone movement across a fragmented habitat

Theoretical and empirical studies suggest that the structure and position of hybrid zones can change over time. Evidence for moving hybrid zones has been directly inferred by repeated sampling over time, or indirectly through the detection of genetic footprints left by the receding species and the resulting asymmetric patterns of introgression across markers. We here investigate a hybrid zone formed by two subspecies of the Iberian golden-striped salamander, Chioglossa lusitanica, using a panel of 35 nuclear loci (31 SNPs and 4 allozymes) and one mitochondrial locus in a transect in central Portugal. We found concordant and coincident clines for most of the nuclear loci (n=22, 63%), defining a narrow hybrid zone of ca. 6 km wide, with the centre positioned ca. 15 km south of the Mondego river. Asymmetric introgression was observed at another 14 loci. Their clines are displaced towards the north, with positions located either close to the Mondego river (n=6), or further northwards (n=8). We interpret these profiles as genetic traces of the southward displacement of C. lusitanica lusitanica by C. l. longipes over the wider Mondego river valley. We noted the absence of significant linkage disequilibrium and we inferred low levels of effective selection per locus against hybrids, suggesting that introgression in the area of species replacement occurred under a neutral diffusion process. A species distribution model suggests that the C. lusitanica hybrid zone coincides with a narrow corridor of fragmented habitat. From the position of the displaced clines, we infer that patches of locally suitable habitat trapped some genetic variants that became disassociated from the southward moving hybrid zone. This study highlights the influence of habitat availability on hybrid zone movement.

Simulated patterns of mitochondrial diversity are consistent with partial population turnover in Bronze Age Central Europe

OBJECTIVES

The analysis of ancient mitochondrial DNA from osteological remains has challenged previous conclusions drawn from the analysis of mitochondrial DNA from present populations, notably by revealing an absence of genetic continuity between the Neolithic and modern populations in Central Europe. Our study investigates how to reconcile these contradictions at the mitochondrial level using a modeling approach.

MATERIALS AND METHODS

We used a spatially explicit computational framework to simulate ancient and modern DNA sequences under various evolutionary scenarios of post Neolithic demographic events and compared the genetic diversity of the simulated and observed mitochondrial sequences. We investigated which-if any-scenarios were able to reproduce statistics of genetic diversity similar to those observed, with a focus on the haplogroup N1a, associated with the spread of early Neolithic farmers.

RESULTS

Demographic fluctuations during the Neolithic transition or subsequent demographic collapses after this period, that is, due to epidemics such as plague, are not sufficient to explain the signal of population discontinuity detected on the mitochondrial DNA in Central Europe. Only a scenario involving a substantial genetic input due to the arrival of migrants after the Neolithic transition, possibly during the Bronze Age, is compatible with observed patterns of genetic diversity.

DISCUSSION

Our results corroborate paleogenomic studies, since out of the alternative hypotheses tested, the best one that was able to recover observed patterns of mitochondrial diversity in modern and ancient Central European populations was one were immigration of populations from the Pontic steppes during the Bronze Age was explicitly simulated.

Genetic and morphological data demonstrate hybridization and backcrossing in a pair of salamanders at the far end of the speciation continuum

Deeply diverged yet hybridizing species provide a system to investigate the final stages of the speciation process. We study a hybridizing pair of salamander species-the morphologically and genetically drastically different newts Triturus cristatus and T. marmoratus-with a panel of 32 nuclear and mitochondrial genetic markers. Morphologically identified hybrids are mostly of the F1 generation and mothered by T. cristatus. The sex ratio of the F1 hybrid class is reciprocally skewed, with a preponderance of females in T. cristatus-mothered hybrids and males in T. marmoratus-mothered hybrids. This amounts to the Haldane effect operating in one direction of the cross. Deeper generation hybrids are occasionally produced, possibly including F1 hybrid × backcross hybrid offspring. Interspecific gene flow is low, yet skewed toward T. cristatus. This asymmetry may be caused by hybrid zone movement, with the superseding species being predisposed to introgression. The persisting gene flow between deeply differentiated species supports the notion that full genetic isolation may be selected against. Conversely, published morphological data suggest that introgressive hybridization is detrimental, with digital malformations occurring more frequently in the area of sympatry. Finally, to assist field identification, both within the area of natural range overlap and concerning anthropogenic introductions elsewhere, we document the phenotypical variation of two generations of hybrids compared with both parental species. We suggest that fluctuating range boundaries, ecological segregation, cytonuclear incompatibilities and hybrid breakdown through Bateson-Dobzhansky-Muller incompatibilities all contribute to species integrity, despite incomplete isolation during secondary contact.

Spatially explicit paleogenomic simulations support cohabitation with limited admixture between Bronze Age Central European populations

The Bronze Age is a complex period of social, cultural and economic changes. Recent paleogenomic studies have documented a large and rapid genetic change in early Bronze Age populations from Central Europe. However, the detailed demographic and genetic processes involved in this change are still debated. Here we have used spatially explicit simulations of genomic components to better characterize the demographic and migratory conditions that may have led to this change. We investigated various scenarios representing the expansion of pastoralists from the Pontic steppe, potentially linked to the Yamnaya cultural complex, and their interactions with local populations in Central Europe, considering various eco-evolutionary factors, such as population admixture, competition and long-distance dispersal. Our results do not support direct competition but rather the cohabitation of pastoralists and farmers in Central Europe, with limited gene flow between populations. They also suggest occasional long-distance migrations accompanying the expansion of pastoralists and a demographic decline in both populations following their initial contact. These results link recent archaeological and paleogenomic observations and move further the debate of genomic changes during the early Bronze Age.

The Genomic Processes of Biological Invasions: From Invasive Species to Cancer Metastases and Back Again

The concept of invasion is useful across a broad range of contexts, spanning from the fine scale landscape of cancer tumors up to the broader landscape of ecosystems. Invasion biology provides extraordinary opportunities for studying the mechanistic basis of contemporary evolution at the molecular level. Although the field of invasion genetics was established in ecology and evolution more than 50 years ago, there is still a limited understanding of how genomic level processes translate into invasive phenotypes across different taxa in response to complex environmental conditions. This is largely because the study of most invasive species is limited by information about complex genome level processes. We lack good reference genomes for most species. Rigorous studies to examine genomic processes are generally too costly. On the contrary, cancer studies are fortified with extensive resources for studying genome level dynamics and the interactions among genetic and non-genetic mechanisms. Extensive analysis of primary tumors and metastatic samples have revealed the importance of several genomic mechanisms including higher mutation rates, specific types of mutations, aneuploidy or whole genome doubling and non-genetic effects. Metastatic sites can be directly compared to primary tumor cell counterparts. At the same time, clonal dynamics shape the genomics and evolution of metastatic cancers. Clonal diversity varies by cancer type, and the tumors’ donor and recipient tissues. Still, the cancer research community has been unable to identify any common events that provide a universal predictor of “metastatic potential” which parallels findings in evolutionary ecology. Instead, invasion in cancer studies depends strongly on context, including order of events and clonal composition. The detailed studies of the behavior of a variety of human cancers promises to inform our understanding of genome level dynamics in the diversity of invasive species and provide novel insights for management.

Introgression dynamics from invasive pigs into wild boar following the March 2011 natural and anthropogenic disasters at Fukushima

Natural and anthropogenic disasters have the capability to cause sudden extrinsic environmental changes and long-lasting perturbations including invasive species, species expansion and influence evolution as selective pressures force adaption. Such disasters occurred on 11 March 2011, in Fukushima, Japan, when an earthquake, tsunami and meltdown of a nuclear power plant all drastically reformed anthropogenic land use. Using genetic data, we demonstrate how wild boar (Sus scrofa leucomystax) have persevered against these environmental changes, including an invasion of escaped domestic pigs (Sus scrofa domesticus). Concurrently, we show evidence of successful hybridization between pigs and native wild boar in this area; however in future offspring, the pig legacy has been diluted through time. We speculate that the range expansion dynamics inhibit long-term introgression and introgressed alleles will continue to decrease at each generation while only maternally inherited organelles will persist. Using the gene flow data among wild boar, we assume that offspring from hybrid lineages will continue dispersal north at low frequencies as climates warm. We conclude that future risks for wild boar in this area include intraspecies competition, revitalization of human-related disruptions and disease outbreaks.

Insights into invasive species from whole-genome resequencing

Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.

Evaluating Causes of Current Genetic Gradients of Modern Humans of the Iberian Peninsula

The history of modern humans in the Iberian Peninsula includes a variety of population arrivals sometimes presenting admixture with resident populations. Genetic data from current Iberian populations revealed an overall east–west genetic gradient that some authors interpreted as a direct consequence of the Reconquista, where Catholic Kingdoms expanded their territories toward the south while displacing Muslims. However, this interpretation has not been formally evaluated. Here, we present a qualitative analysis of the causes of the current genetic gradient observed in the Iberian Peninsula using extensive spatially explicit computer simulations based on a variety of evolutionary scenarios. Our results indicate that the Neolithic range expansion clearly produces the orientation of the observed genetic gradient. Concerning the Reconquista (including political borders among Catholic Kingdoms and regions with different languages), if modeled upon a previous Neolithic expansion, it effectively favored the orientation of the observed genetic gradient and shows local isolation of certain regions (i.e., Basques and Galicia). Despite additional evolutionary scenarios could be evaluated to more accurately decipher the causes of the Iberian genetic gradient, here we show that this gradient has a more complex explanation than that previously hypothesized.