Demographically explicit scans for barriers to gene flow using gIMble.

Identifying regions of the genome that act as barriers to gene flow between recently diverged taxa has remained challenging given the many evolutionary forces that generate variation in genetic diversity and divergence along the genome, and the stochastic nature of this variation. Progress has been impeded by a conceptual and methodological divide between analyses that infer the demographic history of speciation and genome scans aimed at identifying locally maladaptive alleles i.e. genomic barriers to gene flow. Here we implement genomewide IM blockwise likelihood estimation (gIMble), a composite likelihood approach for the quantification of barriers, that bridges this divide. This analytic framework captures background selection and selection against barriers in a model of isolation with migration (IM) as heterogeneity in effective population size (Ne) and effective migration rate (me), respectively. Variation in both effective demographic parameters is estimated in sliding windows via pre-computed likelihood grids. gIMble includes modules for pre-processing/filtering of genomic data and performing parametric bootstraps using coalescent simulations. To demonstrate the new approach, we analyse data from a well-studied pair of sister species of tropical butterflies with a known history of post-divergence gene flow: Heliconius melpomene and H. cydno. Our analyses uncover both large-effect barrier loci (including well-known wing-pattern genes) and a genome-wide signal of a polygenic barrier architecture.

Chromosome Fissions and Fusions Act as Barriers to Gene Flow between Brenthis Fritillary Butterflies.

Chromosome rearrangements are thought to promote reproductive isolation between incipient species. However, it is unclear how often, and under what conditions, fission and fusion rearrangements act as barriers to gene flow. Here we investigate speciation between two largely sympatric fritillary butterflies, Brenthis daphne and Brenthis ino. We use a composite likelihood approach to infer the demographic history of these species from whole-genome sequence data. We then compare chromosome-level genome assemblies of individuals from each species and identify a total of nine chromosome fissions and fusions. Finally, we fit a demographic model where effective population sizes and effective migration rate vary across the genome, allowing us to quantify the effects of chromosome rearrangements on reproductive isolation. We show that chromosomes involved in rearrangements experienced less effective migration since the onset of species divergence and that genomic regions near rearrangement points have a further reduction in effective migration rate. Our results suggest that the evolution of multiple rearrangements in the B. daphne and B. ino populations, including alternative fusions of the same chromosomes, have resulted in a reduction in gene flow. Although fission and fusion of chromosomes are unlikely to be the only processes that have led to speciation between these butterflies, this study shows that these rearrangements can directly promote reproductive isolation and may be involved in speciation when karyotypes evolve quickly.

Experimental introgression in Drosophila: Asymmetric postzygotic isolation associated with chromosomal inversions and an incompatibility locus on the X chromosome.

Interspecific gene flow (introgression) is an important source of new genetic variation, but selection against it can reinforce reproductive barriers between interbreeding species. We used an experimental approach to trace the role of chromosomal inversions and incompatibility genes in preventing introgression between two partly sympatric Drosophila virilis group species, D. flavomontana and D. montana. We backcrossed F1 hybrid females from a cross between D. flavomontana female and D. montana male with the males of the parental species for two generations and sequenced pools of parental strains and their reciprocal second generation backcross (BC2 mon and BC2 fla) females. Contrasting the observed amount of introgression (mean hybrid index, HI) in BC2 female pools along the genome to simulations under different scenarios allowed us to identify chromosomal regions of restricted and increased introgression. We found no deviation from the HI expected under a neutral null model for any chromosome for the BC2 mon pool, suggesting no evidence for genetic incompatibilities in backcrosses towards D. montana. In contrast, the BC2 fla pool showed high variation in the observed HI between different chromosomes, and massive reduction of introgression on the X chromosome (large X-effect). This observation is compatible with reduced recombination combined with at least one dominant incompatibility locus residing within the X inversion(s). Overall, our study suggests that genetic incompatibilities arising within chromosomal inversions can play an important role in speciation.

Do chromosome rearrangements fix by genetic drift or natural selection? Insights from Brenthis butterflies.

Large-scale chromosome rearrangements, such as fissions and fusions, are a common feature of eukaryote evolution. They can have considerable influence on the evolution of populations, yet it remains unclear exactly how rearrangements become established and eventually fix. Rearrangements could fix by genetic drift if they are weakly deleterious or neutral, or they may instead be favoured by positive natural selection. Here, we compare genome assemblies of three closely related Brenthis butterfly species and characterize a complex history of fission and fusion rearrangements. An inferred demographic history of these species suggests that rearrangements became fixed in populations with large long-term effective size (Ne ), consistent with rearrangements being selectively neutral or only very weakly underdominant. Using a recently developed analytic framework for characterizing hard selective sweeps, we find that chromosome fusions are not enriched for evidence of past sweeps compared to other regions of the genome. Nonetheless, we do infer a strong and recent selective sweep around one chromosome fusion in the B. daphne genome. Our results suggest that rearrangements in these species likely have weak absolute fitness effects and fix by genetic drift. However, one putative selective sweep raises the possibility that natural selection may sometimes play a role in the fixation of chromosome fusions.

Discordant Pleistocene population size histories in a guild of hymenopteran parasitoids.

Signatures of past changes in population size have been detected in genome-wide variation in many species. However, the causes of such demographic changes and the extent to which they are shared across co-distributed species remain poorly understood. During Pleistocene glacial maxima, many temperate European species were confined to southern refugia. While vicariance and range expansion processes associated with glacial cycles have been widely documented, it is unclear whether refugial populations of co-distributed species have experienced shared histories of population size change. We analyse whole-genome sequence data to reconstruct and compare demographic histories during the Quaternary for Iberian refuge populations in a single ecological guild (seven species of chalcid parasitoid wasps associated with oak cynipid galls). For four of these species, we find support for large changes in effective population size (Ne ) through the Pleistocene that coincide with major climate events. However, there is little evidence that the timing, direction and magnitude of demographic change are shared across species, suggesting that demographic histories in this guild are largely idiosyncratic, even at the scale of a single glacial refugium.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species.

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

Whole-genome data reveal the complex history of a diverse ecological community.

How widespread ecological communities assemble remains a key question in ecology. Trophic interactions between widespread species may reflect a shared population history or ecological fitting of local pools of species with very different population histories. Which scenario applies is central to the stability of trophic associations and the potential for coevolution between species. Here we show how alternative community assembly hypotheses can be discriminated using whole-genome data for component species and provide a likelihood framework that overcomes current limitations in formal comparison of multispecies histories. We illustrate our approach by inferring the assembly history of a Western Palearctic community of insect herbivores and parasitoid natural enemies, trophic groups that together comprise 50% of terrestrial species. We reject models of codispersal from a shared origin and of delayed enemy pursuit of their herbivore hosts, arguing against herbivore attainment of "enemy-free space." The community-wide distribution of species expansion times is also incompatible with a random, neutral model of assembly. Instead, we reveal a complex assembly history of single- and multispecies range expansions through the Pleistocene from different directions and over a range of timescales. Our results suggest substantial turnover in species associations and argue against tight coevolution in this system. The approach we illustrate is widely applicable to natural communities of nonmodel species and makes it possible to reveal the historical backdrop against which natural selection acts.

Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp.

Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems.

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd-RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2-3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster- and species-specific marker analyses revealed that, apart from four early-generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1-3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.

What haplodiploids can teach us about hybridization and speciation.

Most evolutionary theory focuses on species that reproduce through sexual reproduction where both sexes have a diploid chromosome count. Yet a substantial proportion of multicellular species display complex life cycles, with both haploid and diploid life stages. A classic example is haplodiploidy, where females develop from fertilized eggs and are diploid, while males develop from unfertilized eggs and are haploid. Although haplodiploids make up about 15% of all animals (de la Filia et al. ), this type of reproduction is rarely considered in evolutionary theory. In this issue of Molecular Ecology, Patten et al. () develop a theoretical model to compare the rate of nuclear and mitochondrial introgression in haplodiploid and diploid species. They show that when two haplodiploid species hybridize, nuclear genes are much less likely to cross the species barrier than if both species were to be diploids. The reason for this is that only half of the offspring resulting from matings between haplodiploid species are true hybrids: sons from such mating only inherit their mother genes and therefore only contain genes of the maternal species. Truly, hybrid males can only occur through backcrossing of a hybrid female to a male of one of the parental species. While this twist of haplodiploid transmission genetics limits nuclear introgression, mitochondrial genes, which are maternally inherited, are unaffected by the scarcity of hybrid males. In other words, the rate of mitochondrial introgression is the same for haplodiploid and diploid species. As a result, haplodiploid species on average show a bias of mitochondrial compared to nuclear introgression.

Testing models of speciation from genome sequences: divergence and asymmetric admixture in Island South-East Asian Sus species during the Plio-Pleistocene climatic fluctuations.

In many temperate regions, ice ages promoted range contractions into refugia resulting in divergence (and potentially speciation), while warmer periods led to range expansions and hybridization. However, the impact these climatic oscillations had in many parts of the tropics remains elusive. Here, we investigate this issue using genome sequences of three pig (Sus) species, two of which are found on islands of the Sunda-shelf shallow seas in Island South-East Asia (ISEA). A previous study revealed signatures of interspecific admixture between these Sus species (Genome biology, 14, 2013, R107). However, the timing, directionality and extent of this admixture remain unknown. Here, we use a likelihood-based model comparison to more finely resolve this admixture history and test whether it was mediated by humans or occurred naturally. Our analyses suggest that interspecific admixture between Sunda-shelf species was most likely asymmetric and occurred long before the arrival of humans in the region. More precisely, we show that these species diverged during the late Pliocene but around 23% of their genomes have been affected by admixture during the later Pleistocene climatic transition. In addition, we show that our method provides a significant improvement over D-statistics which are uninformative about the direction of admixture.

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.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana.

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

A new genus of oak gallwasp, Cyclocynips Melika, Tang & Sinclair (Hymenoptera: Cynipidae: Cynipini), with descriptions of two new species from Taiwan.

A new genus of cynipid oak gallwasp-Cyclocynips Melika, Tang, & Sinclair (Hymenoptera: Cynipidae: Cynipini), with two new species--C. uberis and C. tumorvirgae--reared from galls on oaks of the Quercus subgenus Cyclobalanopsis is described from Taiwan. Descriptions of asexual generation adults and their diagnostic characters are presented. The likelihood of yet undiscovered sexual generations and the evolution of host-plant associations in these species are discussed.

The genome sequence of the Chalkhill Blue, Lysandra coridon (Poda, 1761).

We present a genome assembly from an individual male Lysandra coridon (the Chalkhill Blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 541 megabases in span. Most of the assembly is scaffolded into 90 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,334 protein coding genes.

The genome sequence of the Large Skipper, Ochlodes sylvanus, (Esper, 1777).

We present a genome assembly from an individual female Ochlodes sylvanus, the Large Skipper (Arthropoda; Insecta; Lepidoptera; Hesperiidae). The genome sequence is 380 megabases in span. Most of the assembly (99.97%) is scaffolded into 30 chromosomal pseudomolecules, including the assembled W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 17.1 kilobases in length. Gene annotation of this assembly on Ensembl identified 13,451 protein coding genes.

The genome sequence of the Mazarine Blue, Cyaniris semiargus (Rottemburg, 1775).

We present a genome assembly from an individual male Cyaniris semiargus (the Mazarine Blue; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequence is 441.5 megabases in span. Most of the assembly is scaffolded into 24 chromosomal pseudomolecules, including the assembled Z sex chromosome. The mitochondrial genome has also been assembled and is 15.4 kilobases in length. Gene annotation of this assembly on Ensembl identified 16,408 protein coding genes.

The genome sequence of the Brown Argus, Aricia agestis (Denis & Schiffermüller, 1775).

We present genome assemblies from two male Aricia agestis specimens (the Brown Argus; Arthropoda; Insecta; Lepidoptera; Lycaenidae). The genome sequences are 435.3 and 437.4 megabases in span. Each assembly is scaffolded into 23 chromosomal pseudomolecules, including the Z sex chromosome. The mitochondrial genomes were assembled and are 15.47 and 15.45 kilobases in length. Gene annotation of these assemblies on Ensembl identified 12,688 and 12,654 protein coding genes.

A likelihood-based comparison of population histories in a parasitoid guild.

Little is known about the stability of trophic relationships in complex natural communities over evolutionary timescales. Here, we use sequence data from 18 nuclear loci to reconstruct and compare the intraspecific histories of major Pleistocene refugial populations in the Middle East, the Balkans and Iberia in a guild of four Chalcid parasitoids (Cecidostiba fungosa, Cecidostiba semifascia, Hobbya stenonota and Mesopolobus amaenus) all attacking Cynipid oak galls. We develop a likelihood method to numerically estimate models of divergence between three populations from multilocus data. We investigate the power of this framework on simulated data, and-using triplet alignments of intronic loci-quantify the support for all possible divergence relationships between refugial populations in the four parasitoids. Although an East to West order of population divergence has highest support in all but one species, we cannot rule out alternative population tree topologies. Comparing the estimated times of population splits between species, we find that one species, M. amaenus, has a significantly older history than the rest of the guild and must have arrived in central Europe at least one glacial cycle prior to other guild members. This suggests that although all four species may share a common origin in the East, they expanded westwards into Europe at different times.

Model-based comparisons of phylogeographic scenarios resolve the intraspecific divergence of cactophilic Drosophila mojavensis.

The cactophilic fly Drosophila mojavensis exhibits considerable intraspecific genetic structure across allopatric geographic regions and shows associations with different host cactus species across its range. The divergence between these populations has been studied for more than 60years, yet their exact historical relationships have not been resolved. We analysed sequence data from 15 intronic X-linked loci across populations from Baja California, mainland Sonora-Arizona and Mojave Desert regions under an isolation-with-migration model to assess multiple scenarios of divergence. We also compared the results with a pre-existing sequence data set of eight autosomal loci. We derived a population tree with Baja California placed at its base and link their isolation to Pleistocene climatic oscillations. Our estimates suggest the Baja California population diverged from an ancestral Mojave Desert/mainland Sonora-Arizona group around 230,000-270,000years ago, while the split between the Mojave Desert and mainland Sonora-Arizona populations occurred one glacial cycle later, 117,000-135,000years ago. Although we found these three populations to be effectively allopatric, model ranking could not rule out the possibility of a low level of gene flow between two of them. Finally, the Mojave Desert population showed a small effective population size, consistent with a historical population bottleneck. We show that model-based inference from multiple loci can provide accurate information on the historical relationships of closely related groups allowing us to set into historical context a classic system of incipient ecological speciation.