Complex patterns of differentiation and gene flow underly the divergence of aposematic phenotypes in Oophaga poison frogs

Hybridization and introgression can have complex consequences for both species evolution and conservation. Here, we investigated the origin and characteristics of a putative hybrid zone between two South American poison dart frog species, Oophaga anchicayensis and the critically endangered Oophaga lehmanni, which are heavily sought after on the illegal pet market. Using a combination of phenotypic (49 traits) and genomic (ddRADseq) data, we found that the putative hybrids are morphologically distinct from their parental species and confirmed genomic signatures of admixture in these populations. Several lines of evidence (hybrid indices, interspecific hybrid heterozygosity, genomic clines, comparisons with simulated hybrids and demographic modelling) support the conclusion that these populations are not comprised of early‐generation hybrids and thus, they probably did not arise as a result of illegal translocations associated with wildlife trafficking. Instead, they probably represent an independent lineage which has persisted through isolation and has only relatively recently re‐established gene flow with both parental species. Furthermore, we detected signals of differential introgression from parental species into these hybrid populations which suggest relaxed stabilizing selection on these aposematic colour morphs, potentially via context‐dependent female choice. These populations thus provide a fascinating window into the role of hybridization, isolation and female choice in the diversification of South American poison dart frogs. In addition, our results underline the importance of landscape conservation measures to protect, not only known localities of nominal species, but also the phenotypic and genomic variation harbored by admixed lineages which represent crucial repositories for the impressive diversity in this system.     

Admixture and the organization of genetic diversity in a butterfly species complex revealed through common and rare genetic variants

Detailed information about the geographic distribution of genetic and genomic variation is necessary to better understand the organization and structure of biological diversity. In particular, spatial isolation within species and hybridization between them can blur species boundaries and create evolutionary relationships that are inconsistent with a strictly bifurcating tree model. Here, we analyse genome‐wide DNA sequence and genetic ancestry variation in Lycaeides butterflies to quantify the effects of admixture and spatial isolation on how biological diversity is organized in this group. We document geographically widespread and pervasive historical admixture, with more restricted recent hybridization. This includes evidence supporting previously known and unknown instances of admixture. The genome composition of admixed individuals varies much more among than within populations, and tree‐ and genetic ancestry‐based analyses indicate that multiple distinct admixed lineages or populations exist. We find that most genetic variants in Lycaeides are rare (minor allele frequency <0.5%). Because the spatial and taxonomic distributions of alleles reflect demographic and selective processes since mutation, rare alleles, which are presumably younger than common alleles, were spatially and taxonomically restricted compared with common variants. Thus, we show patterns of genetic variation in this group are multifaceted, and we argue that this complexity challenges simplistic notions concerning the organization of biological diversity into discrete, easily delineated and hierarchically structured entities.     

Genomic consequences of a recent three‐way admixture in supplemented wild brown trout populations revealed by local ancestry tracts

Understanding the evolutionary consequences of human‐mediated introductions of domesticated strains into the wild and their subsequent admixture with natural populations is of major concern in conservation biology. However, the genomic impacts of stocking from distinct sources (locally derived vs. divergent) on the genetic integrity of wild populations remain poorly understood. We designed an approach based on estimating local ancestry along individual chromosomes to provide a detailed picture of genomic admixture in supplemented populations. We used this approach to document admixture consequences in the brown trout Salmo trutta, for which decades of stocking practices have profoundly impacted the genetic make‐up of wild populations. In southern France, small local Mediterranean populations have been subject to successive introductions of domestic strains derived from the Atlantic and Mediterranean lineages. To address the impact of stocking, we evaluate the extent of admixture from both domestic strains within populations, using 75,684 mapped SNPs obtained from double‐digested restriction site‐associated DNA sequencing. Then, the chromosomal ancestry profiles of admixed individuals reveal a wider diversity of hybrid and introgressed genotypes than estimated using classical methods for inferring ancestry and hybrid pedigrees. In addition, the length distribution of introgressed tracts retained different timings of introgression between the two domestic strains. We finally reveal opposite consequences of admixture on the level of polymorphism of the recipient populations between domestic strains. Our study illustrates the potential of using the information contained in the genomic mosaic of ancestry tracts in combination with classical methods based on allele frequencies for analysing multiple‐way admixture with population genomic data.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Mixed signals from hybrid genomes

Admixture results from interbreeding between individuals from different populations or species that were previously genetically isolated from each other (Fig.  1 ). Identifying admixture events in the genome is not always a straightforward task, because the genetic signature left behind fades with time as recombination events fragment the genomic segments introduced during the interbreeding event. Additionally, when the genetic architecture of populations or species that admix is not very different (e.g. they coalesce to a common ancestor recently), admixture signatures may be difficult to detect. Ignoring the effects of admixture can, however, pose severe problems for population genetic analyses that rely on the distribution of polymorphic markers across the genome. In this issue of Molecular Ecology, Bosse et al. ( 2014 ) analyse genomic data from modern pigs to understand hybridization processes that occurred between domestic pigs from European and Asiatic origin, and between pigs and wild boars. Their results are interesting regarding the fine‐scale distribution of admixture across the pig genome, and the way in which this admixture biases estimates of the effective population size in European domestic pigs. The implications of these results are significant, as they serve as a cautionary note on genomic analyses that depend on the distribution of polymorphic variants in potentially admixed populations.     

Genomic sequencing reveals historical,demographic and selective factors associated with the diversification of the fire‐associated fungus Neurospora discreta

Delineating microbial populations, discovering ecologically relevant phenotypes and identifying migrants, hybrids or admixed individuals have long proved notoriously difficult, thereby limiting our understanding of the evolutionary forces at play during the diversification of microbial species. However, recent advances in sequencing and computational methods have enabled an unbiased approach whereby incipient species and the genetic correlates of speciation can be identified by examining patterns of genomic variation within and between lineages. We present here a population genomic study of a phylogenetic species in the Neurospora discreta species complex, based on the resequencing of full genomes (~37 Mb) for 52 fungal isolates from nine sites in three continents. Population structure analyses revealed two distinct lineages in South–East Asia, and three lineages in North America/Europe with a broad longitudinal and latitudinal range and limited admixture between lineages. Genome scans for selective sweeps and comparisons of the genomic landscapes of diversity and recombination provided no support for a role of selection at linked sites on genomic heterogeneity in levels of divergence between lineages. However, demographic inference indicated that the observed genomic heterogeneity in divergence was generated by varying rates of gene flow between lineages following a period of isolation. Many putative cases of exchange of genetic material between phylogenetically divergent fungal lineages have been discovered, and our work highlights the quantitative importance of genetic exchanges between more closely related taxa to the evolution of fungal genomes. Our study also supports the role of allopatric isolation as a driver of diversification in saprobic microbes.     

Estimating Individual Admixture Proportions from Next Generation Sequencing Data

Inference of population structure and individual ancestry is important both for population genetics and for association studies. With next generation sequencing technologies it is possible to obtain genetic data for all accessible genetic variations in the genome. Existing methods for admixture analysis rely on known genotypes. However, individual genotypes cannot be inferred from low-depth sequencing data without introducing errors. This article presents a new method for inferring an individual’s ancestry that takes the uncertainty introduced in next generation sequencing data into account. This is achieved by working directly with genotype likelihoods that contain all relevant information of the unobserved genotypes. Using simulations as well as publicly available sequencing data, we demonstrate that the presented method has great accuracy even for very low-depth data. At the same time, we demonstrate that applying existing methods to genotypes called from the same data can introduce severe biases. The presented method is implemented in the NGSadmix software available at http://www.popgen.dk/software.     

Adaptive admixture in the West African bovine hybrid zone: insight from the Borgou population

Understanding the adaptive response to environmental fluctuations represents a central issue in evolutionary biology. Population admixture between divergent ancestries has often been considered as an efficient short‐term adaptation strategy. Cattle populations from the West African Bos taurus × Bos indicus hybrid zone represent a valuable resource to characterize the effect of such adaptive admixture at the genome level. We here provide a detailed assessment of the global and local genome ancestries of the Borgou breed, one of the most representative cattle of this hybrid zone. We analysed a large data set consisting of 38 100 SNPs genotyped on 203 Borgou and 591 individuals representative of all the different cattle ancestries. At the global genomic level, we show that Borgou is a stabilized admixed breed whose origin (c. 130 years ago) traces back to the great African rinderpest pandemic, several centuries after the last admixture events, the West African zebus originate from (c. 500 years ago). To identify footprints of adaptive admixture, we combined the identification of signatures of selection and the functional annotation of the underlying genes using systems biology tools. The detection of the SILV coat coloration gene likely under artificial selection may be viewed as a validation of our approach. Overall, our results suggest that the long‐term presence of pathogens and the intermediate environmental conditions are the main acting selective pressures. Our analytical framework can be extended to other model or nonmodel species to understand the process that shapes the patterns of genetic variability in hybrid zones.     

Insight into the roles of selection in speciation from genomic patterns of divergence and introgression in secondary contact in venomous rattlesnakes

Investigating secondary contact of historically isolated lineages can provide insight into how selection and drift influence genomic divergence and admixture. Here, we studied the genomic landscape of divergence and introgression following secondary contact between lineages of the Western Diamondback Rattlesnake (Crotalus atrox) to determine whether genomic regions under selection in allopatry also contribute to reproductive isolation during introgression. We used thousands of nuclear loci to study genomic differentiation between two lineages that have experienced recent secondary contact following isolation, and incorporated sampling from a zone of secondary contact to identify loci that are resistant to gene flow in hybrids. Comparisons of patterns of divergence and introgression revealed a positive relationship between allelic differentiation and resistance to introgression across the genome, and greater‐than‐expected overlap between genes linked to lineage‐specific divergence and loci that resist introgression. Genes linked to putatively selected markers were related to prominent aspects of rattlesnake biology that differ between populations of Western Diamondback rattlesnakes (i.e., venom and reproductive phenotypes). We also found evidence for selection against introgression of genes that may contribute to cytonuclear incompatibility, consistent with previously observed biased patterns of nuclear and mitochondrial alleles suggestive of partial reproductive isolation due to cytonuclear incompatibilities. Our results provide a genome‐scale perspective on the relationships between divergence and introgression in secondary contact that is relevant for understanding the roles of selection in maintaining partial isolation of lineages, causing admixing lineages to not completely homogenize.     

Locus‐specific ancestry to detect recent response to selection in admixed Swiss Fleckvieh cattle

Identification of selection signatures is one of the current endeavors of evolutionary genetics. Admixed populations may be used to infer post‐admixture selection. We calculated local ancestry for Swiss Fleckvieh, a composite of Simmental (SI) and Red Holstein Friesian (RHF), to infer such signals. Illumina Bovine SNP50 BeadChip data for 300 admixed, 88 SI and 97 RHF bulls were used. The average RHF ancestry across the whole genome was 0.70. To identify regions with high deviation from average, we considered two significance thresholds, based on a permutation test and extreme deviation from normal distribution. Regions on chromosomes 13 (46.3–47.3 Mb) and 18 (18.7–25.9 Mb) passed both thresholds in the direction of increased SI. Extended haplotype homozygosity within (iHS) and between (Rsb) populations was calculated to explore additional patterns of pre‐ and post‐admixture selection signals. The Rsb score of admixed and SI was significant in a wide region of chromosome 18 (6.6–24.6 Mb) overlapped with one area of strong local ancestry deviation. FTO, with pleiotropic effect on milk and fertility, NOD2 on dairy and NKD1 and SALL1 on fertility traits are located there. Genetic differentiation of RHF and SI (F_st), an alternative indicator of pre‐admixture selection in pure populations, was calculated. No considerable overlap of peaks of local ancestry deviations and F_st was observed. We found two regions with significant signatures of post‐admixture selection in this very young composite, applying comparatively stringent significance thresholds. The signals cover relatively large genomic areas and did not allow pinpointing of the gene(s) responsible for the apparent shift in ancestry proportions.     

Patterns of Admixture and Population Structure in Native Populations of Northwest North America

The initial contact of European populations with indigenous populations of the Americas produced diverse admixture processes across North, Central, and South America. Recent studies have examined the genetic structure of indigenous populations of Latin America and the Caribbean and their admixed descendants, reporting on the genomic impact of the history of admixture with colonizing populations of European and African ancestry. However, relatively little genomic research has been conducted on admixture in indigenous North American populations. In this study, we analyze genomic data at 475,109 single-nucleotide polymorphisms sampled in indigenous peoples of the Pacific Northwest in British Columbia and Southeast Alaska, populations with a well-documented history of contact with European and Asian traders, fishermen, and contract laborers. We find that the indigenous populations of the Pacific Northwest have higher gene diversity than Latin American indigenous populations. Among the Pacific Northwest populations, interior groups provide more evidence for East Asian admixture, whereas coastal groups have higher levels of European admixture. In contrast with many Latin American indigenous populations, the variance of admixture is high in each of the Pacific Northwest indigenous populations, as expected for recent and ongoing admixture processes. The results reveal some similarities but notable differences between admixture patterns in the Pacific Northwest and those in Latin America, contributing to a more detailed understanding of the genomic consequences of European colonization events throughout the Americas.     

The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm

Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.     

AN EVALUATION OF THE HYBRID SPECIATION HYPOTHESIS FOR XIPHOPHORUS CLEMENCIAE BASED ON WHOLE GENOME SEQUENCES

Once thought rare in animal taxa, hybridization has been increasingly recognized as an important and common force in animal evolution. In the past decade, a number of studies have suggested that hybridization has driven speciation in some animal groups. We investigate the signature of hybridization in the genome of a putative hybrid species, Xiphophorus clemenciae, through whole genome sequencing of this species and its hypothesized progenitors. Based on analysis of this data, we find that X. clemenciae is unlikely to have been derived from admixture between its proposed parental species. However, we find significant evidence for recent gene flow between Xiphophorus species. Although we detect genetic exchange in two pairs of species analyzed, the proportion of genomic regions that can be attributed to hybrid origin is small, suggesting that strong behavioral premating isolation prevents frequent hybridization in Xiphophorus. The direction of gene flow between species is potentially consistent with a role for sexual selection in mediating hybridization.     

Becoming pure: identifying generational classes of admixed individuals within lesser and greater scaup populations

Estimating the frequency of hybridization is important to understand its evolutionary consequences and its effects on conservation efforts. In this study, we examined the extent of hybridization in two sister species of ducks that hybridize. We used mitochondrial control region sequences and 3589 double‐digest restriction‐associated DNA sequences (ddRADseq) to identify admixture between wild lesser scaup (Aythya affinis) and greater scaup (A. marila). Among 111 individuals, we found one introgressed mitochondrial DNA haplotype in lesser scaup and four in greater scaup. Likewise, based on the site‐frequency spectrum from autosomal DNA, gene flow was asymmetrical, with higher rates from lesser into greater scaup. However, using ddRADseq nuclear DNA, all individuals were assigned to their respective species with >0.95 posterior assignment probability. To examine the power for detecting admixture , we simulated a breeding experiment in which empirical data were used to create F1 hybrids and nine generations (F2–F10) of backcrossing. F1 hybrids and F2, F3 and most F4 backcrosses were clearly distinguishable from pure individuals, but evidence of admixed histories was effectively lost after the fourth generation. Thus, we conclude that low interspecific assignment probabilities (0.011–0.043) for two lesser and nineteen greater scaup were consistent with admixed histories beyond the F3 generation. These results indicate that the propensity of these species to hybridize in the wild is low and largely asymmetric. When applied to species‐specific cases, our approach offers powerful utility for examining concerns of hybridization in conservation efforts, especially for determining the generational time until admixed histories are effectively lost through backcrossing.     

Seascape genetics along environmental gradients in the Arabian Peninsula: insights from ddRAD sequencing of anemonefishes

Understanding the processes that shape patterns of genetic structure across space is a central aim of landscape genetics. However, it remains unclear how geographical features and environmental variables shape gene flow, particularly for marine species in large complex seascapes. Here, we evaluated the genomic composition of the two‐band anemonefish Amphiprion bicinctus across its entire geographical range in the Red Sea and Gulf of Aden, as well as its close relative, Amphiprion omanensis endemic to the southern coast of Oman. Both the Red Sea and the Arabian Sea are complex and environmentally heterogeneous marine systems that provide an ideal scenario to address these questions. Our findings confirm the presence of two genetic clusters previously reported for A. bicinctus in the Red Sea. Genetic structure analyses suggest a complex seascape configuration, with evidence of both isolation by distance (IBD) and isolation by environment (IBE). In addition to IBD and IBE, genetic structure among sites was best explained when two barriers to gene flow were also accounted for. One of these coincides with a strong oligotrophic–eutrophic gradient at around 16–20?N in the Red Sea. The other agrees with a historical bathymetric barrier at the straight of Bab al Mandab. Finally, these data support the presence of interspecific hybrids at an intermediate suture zone at Socotra and indicate complex patterns of genomic admixture in the Gulf of Aden with evidence of introgression between species. Our findings highlight the power of recent genomic approaches to resolve subtle patterns of gene flow in marine seascapes.     

海三棱藨草及其近缘种基因组大小的测定

基因组大小是物种基因组的重要特征,通常用DNA C值来衡量,能够用于快速判断基因组倍性,并为分类学与进化生物学提供重要依据。海三棱藨草(Scirpus mariqueter)是长江口和杭州湾具有重要生态意义的标志性物种,被认为是扁秆藨草(S. planiculmis)和藨草(S. triqueter)的杂交种,因染色体小而难以准确确定倍性。近年来,部分研究者指出该物种的分类和命名存在疑点。该研究通过基因组Survey分析检测海三棱藨草样本CJ1的基因组特征,测序深度约为120 ×,并以绿豆(Vigna radiata)为参考标准,利用流式细胞术测定了海三棱藨草及其同域近缘种扁秆藨草和藨草以及海三棱藨草和扁秆藨草的杂交F1共13个样本的DNA C值和相对倍性。结果表明:(1)基因组Survey分析测得CJ1的基因组大小为244.12 Mbp,杂合率为0.68%,重复序列比例为42.38%,GC含量为37.25%。(2)流式细胞术测得来自不同区域的海三棱藨草各样本的基因组倍性相同,1C值在234.87 ～ 242.5 Mbp之间,其中CJ1的基因组大小与基因组Survey检测结果高度一致。(3)扁秆藨草的1C值在251.77 ～ 264.13 Mbp之间,藨草1C值为537.33 Mbp。根据上述基因组大小,认为海三棱藨草不可能是这两者的杂交种。该研究补充了海三棱藨草及其近缘种的基因组特征,为后续全基因组测序奠定基础,同时也否定了海三棱藨草起源于扁杆藨草和藨草杂交的假说。     

Detecting introgressive hybridization between free-ranging domestic dogs and wild wolves (Canis lupus) by admixture linkage disequilibrium analysis

Occasional crossbreeding between free-ranging domestic dogs and wild wolves (Canis lupus) has been detected in some European countries by mitochondrial DNA sequencing and genotyping unlinked microsatellite loci. Maternal and unlinked genomic markers, however, might underestimate the extent of introgressive hybridization, and their impacts on the preservation of wild wolf gene pools. In this study, we genotyped 220 presumed Italian wolves, 85 dogs and 7 known hybrids at 16 microsatellites belonging to four different linkage groups (plus four unlinked microsatellites). Population clustering and individual assignments were performed using a Bayesian procedure implemented in structure 2.1, which models the gametic disequilibrium arising between linked loci during admixtures, aiming to trace hybridization events further back in time and infer the population of origin of chromosomal blocks. Results indicate that (i) linkage disequilibrium was higher in wolves than in dogs; (ii) 11 out of 220 wolves (5.0%) were likely admixed, a proportion that is significantly higher than one admixed genotype in 107 wolves found previously in a study using unlinked markers; (iii) posterior maximum-likelihood estimates of the recombination parameter r revealed that introgression in Italian wolves is not recent, but could have continued for the last 70 (+/- 20) generations, corresponding to approximately 140-210 years. Bayesian clustering showed that, despite some admixture, wolf and dog gene pools remain sharply distinct (the average proportions of membership to wolf and dog clusters were Q(w) = 0.95 and Q(d) = 0.98, respectively), suggesting that hybridization was not frequent, and that introgression in nature is counteracted by behavioural or selective constraints.     

Excessive and asymmetrical removal of heterozygous sites by maxSH biases downstream population genetic inference: Implications for hybridization between two primroses

Techniques of reduced-representation sequencing (RRS) have revolutionized ecological and evolutionary genomics studies. Precise establishment of orthologs is a critical challenge for RRS, especially when a reference genome is absent. The proportion of shared heterozygous sites across samples is an alternative criterion for filtering paralogs. In the prevailing pipeline for variant calling of RRS data – PYRAD/IPYRAD, maxSH is an often overlooked parameter with implications to detecting and filtering paralogs according to shared heterozygosity. Using empirical genotyping by sequencing data of two primroses (Primula alpicola Stapf and Primula florindae Ward) and their putative hybrids, and extra data sets of Californian golden cup oaks, we explore the impact of maxSH on filtering paralogs and further downstream analyses. Our study sheds light on the simultaneous validity and risk of filtering paralogs using maxSH, and its significant effects on downstream analyses of outlier detection, population assignment, and demographic modeling, emphasizing the importance of attention to detail during bioinformatic processes. The mutual confirmation between results of population assignment and demographic modeling in this study suggested maxSH = 0.10 has a potentially excessive and asymmetrical effect on the removal of truly shared heterozygous sites as paralogs. These results indicate that hybridization origin hypotheses of putative hybrids represented by results with maxSH = 0.25 and 0.50 are more credible. In conclusion, we revealed the critical hazard of paralogs filtration according to sharing heterozygosity at first, so that we propose to use specific protocols, rather than maxSH, to filter potential paralogs for closely related lineages.     

Likelihood‐based inference of population history from low‐coverage de novo genome assemblies

Short‐read sequencing technologies have in principle made it feasible to draw detailed inferences about the recent history of any organism. In practice, however, this remains challenging due to the difficulty of genome assembly in most organisms and the lack of statistical methods powerful enough to discriminate between recent, nonequilibrium histories. We address both the assembly and inference challenges. We develop a bioinformatic pipeline for generating outgroup‐rooted alignments of orthologous sequence blocks from de novo low‐coverage short‐read data for a small number of genomes, and show how such sequence blocks can be used to fit explicit models of population divergence and admixture in a likelihood framework. To illustrate our approach, we reconstruct the Pleistocene history of an oak‐feeding insect (the oak gallwasp Biorhiza pallida), which, in common with many other taxa, was restricted during Pleistocene ice ages to a longitudinal series of southern refugia spanning the Western Palaearctic. Our analysis of sequence blocks sampled from a single genome from each of three major glacial refugia reveals support for an unexpected history dominated by recent admixture. Despite the fact that 80% of the genome is affected by admixture during the last glacial cycle, we are able to infer the deeper divergence history of these populations. These inferences are robust to variation in block length, mutation model and the sampling location of individual genomes within refugia. This combination of de novo assembly and numerical likelihood calculation provides a powerful framework for estimating recent population history that can be applied to any organism without the need for prior genetic resources.     

Recombinant hybrids retain heterozygosity at many loci: new insights into the genomics of reproductive isolation in Populus

The maintenance of species barriers in the face of gene flow is often thought to result from strong selection against intermediate genotypes, thereby preserving genetic differentiation. Most speciation genomic studies thus aim to identify exceptionally divergent loci between populations, but divergence will be affected by many processes other than reproductive isolation (RI) and speciation. Through genomic studies of recombinant hybrids sampled in the wild, genetic variation associated with RI can be observed in situ, because selection against incompatible genotypes will leave detectable patterns of variation in the hybrid genomes. To better understand the mechanisms directly involved in RI, we investigated three natural ‘replicate’ hybrid zones between two divergent Populus species via locus‐specific patterns of ancestry across recombinant hybrid genomes. As expected, genomic patterns in hybrids and their parental species were consistent with the presence of underdominant selection at several genomic regions. Surprisingly, many loci displayed greatly increased between‐species heterozygosity in recombinant hybrids despite striking genetic differentiation between the parental genomes, the opposite of what would be expected with selection against intermediate genotypes. Only a limited, reproducible set of genotypic combinations was present in hybrid genomes across localities. In the absence of clearly delimited ‘hybrid habitats’, our results suggest that complex epistatic interactions within genomes play an important role in advanced stages of RI between these ecologically divergent forest trees. This calls for more genomic studies that test for unusual patterns of genomic ancestry in hybridizing species.