Search

The concept of a “speciation continuum” has gained popularity in recent decades. It emphasizes speciation as a continuous process that may be studied by comparing contemporary population pairs that show differing levels of divergence. In their recent perspective article in Evolution, Stankowski and Ravinet provided a valuable service by formally defining the speciation continuum as a continuum of reproductive isolation, based on opinions gathered from a survey of speciation researchers. While we agree that the speciation continuum has been a useful concept to advance the understanding of the speciation process, some intrinsic limitations exist. Here, we advocate for a multivariate extension, the speciation hypercube, first proposed by Dieckmann et al. in 2004, but rarely used since. We extend the idea of the speciation cube and suggest it has strong conceptual and practical advantages over a one-dimensional model. We illustrate how the speciation hypercube can be used to visualize and compare different speciation trajectories, providing new insights into the processes and mechanisms of speciation. A key strength of the speciation hypercube is that it provides a unifying framework for speciation research, as it allows questions from apparently disparate subfields to be addressed in a single conceptual model.

Hybridization is thought to reactivate transposable elements (TEs) that were efficiently suppressed in the genomes of the parental hosts. Here, we provide evidence for this ‘genomic shock hypothesis’ in the fission yeast Schizosaccharomyces pombe. In this species, two divergent lineages (Sp and Sk) have experienced recent, likely human induced, hybridization. We used long-read sequencing data to assemble genomes of 37 samples derived from 31 S. pombe strains spanning a wide range of ancestral admixture proportions. A comprehensive TE inventory revealed exclusive presence of long terminal repeat (LTR) retrotransposons. Sequence analysis of active full-length elements, as well as solo-LTRs, revealed a complex history of homologous recombination. Population genetic analyses of syntenic sequences placed insertion of many solo-LTRs prior to the split of the Sp and Sk lineages. Most full-length elements were inserted more recently after hybridization. With the exception of a single full-length element with signs of positive selection, both solo-LTRs, and in particular, full-length elements carried signatures of purifying selection indicating effective removal by the host. Consistent with reactivation upon hybridization, the number of full-length LTR retrotransposons, varying extensively from zero to 87 among strains, significantly increased with the degree of genomic admixture. This study gives a detailed account of global TE diversity in S. pombe, documents complex recombination histories within TE elements and provides evidence for the ‘genomic shock hypothesis’.

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal ( Arctocephalus galapagoensis ) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

With the advent of chromatin-interaction maps, chromosome-level genome assemblies have become a reality for a wide range of organisms. Scaffolding quality is, however, difficult to judge. To explore this gap, we generated multiple chromosome-scale genome assemblies of an emerging wild animal model for carcinogenesis, the California sea lion (Zalophus californianus). Short-read assemblies were scaffolded with two independent chromatin interaction mapping data sets (Hi-C and Chicago), and long-read assemblies with three data types (Hi-C, optical maps and 10X linked reads) following the "Vertebrate Genomes Project (VGP)" pipeline. In both approaches, 18 major scaffolds recovered the karyotype (2n = 36), with scaffold N50s of 138 and 147 Mb, respectively. Synteny relationships at the chromosome level with other pinniped genomes (2n = 32-36), ferret (2n = 34), red panda (2n = 36) and domestic dog (2n = 78) were consistent across approaches and recovered known fissions and fusions. Comparative chromosome painting and multicolour chromosome tiling with a panel of 264 genome-integrated single-locus canine bacterial artificial chromosome probes provided independent evaluation of genome organization. Broad-scale discrepancies between the approaches were observed within chromosomes, most commonly in translocations centred around centromeres and telomeres, which were better resolved in the VGP assembly. Genomic and cytological approaches agreed on near-perfect synteny of the X chromosome, and in combination allowed detailed investigation of autosomal rearrangements between dog and sea lion. This study presents high-quality genomes of an emerging cancer model and highlights that even highly fragmented short-read assemblies scaffolded with Hi-C can yield reliable chromosome-level scaffolds suitable for comparative genomic analyses.

Adaptive divergence is the key evolutionary process generating biodiversity by means of natural selection. Yet, the conditions under which it can arise in the presence of gene flow remain contentious. To address this question, we subjected 132 sexually reproducing fission yeast populations, sourced from two independent genetic backgrounds, to disruptive ecological selection and manipulated the level of migration between environments. Contrary to theoretical expectations, adaptive divergence was most pronounced when migration was either absent (allopatry) or maximal (sympatry), but was much reduced at intermediate rates (parapatry and local mating). This effect was apparent across central life-history components (survival, asexual growth and mating) but differed in magnitude between ancestral genetic backgrounds. The evolution of some fitness components was constrained by pervasive negative correlations (trade-off between asexual growth and mating), while others changed direction under the influence of migration (for example, survival and mating). In allopatry, adaptive divergence was mainly conferred by standing genetic variation and resulted in ecological specialization. In sympatry, divergence was mainly mediated by novel mutations enriched in a subset of genes and was characterized by the repeated emergence of two strategies: an ecological generalist and an asexual growth specialist. Multiple loci showed consistent evidence for antagonistic pleiotropy across migration treatments providing a conceptual link between adaptation and divergence. This evolve-and-resequence experiment shows that rapid ecological differentiation can arise even under high rates of gene flow. It further highlights that adaptive trajectories are governed by complex interactions of gene flow, ancestral variation and genetic correlations.

The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here, we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genus Arctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 and 2.5 Ma. Otaria diverged first, followed by Phocarctos and then four major lineages within Arctocephalus. However, we found Zalophus to be nonmonophyletic, with California (Zalophus californianus) and Steller sea lions (Eumetopias jubatus) grouping closer than the Galapagos sea lion (Zalophus wollebaeki) with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family. [Hybridization; ILS; phylogenomics; Pleistocene; Pliocene; monophyly.]

Recombination is a central biological process with implications for many areas in the life sciences. Yet we are only beginning to appreciate variation in the recombination rate along the genome and among individuals, populations and species. Spurred by technological advances, we are now able to bring variation in this key biological parameter to centre stage. Here, we review the conceptual implications of recombination rate variation and guide the reader through the assumptions, strengths and weaknesses of genomic inference methods, including population-based, pedigree-based and gamete-based approaches. Appreciation of the differences and commonalities of these approaches is a prerequisite to formulate a unifying and comparative framework for understanding the molecular and evolutionary mechanisms shaping, and being shaped by, recombination.

Within hybrid zones of socially monogamous species, the number of mating opportunities with a conspecific can be limited. As a consequence, individuals may mate with a heterospecific (social) partner despite possible fitness costs to their hybrid offspring. Extra‐pair copulations with a conspecific may thus arise as a possible post hoc strategy to reduce the costs of hybridization. We here assessed the rate of extra‐pair paternity in the hybrid zone between all‐black carrion crows (Corvus (corone) corone ) and grey hooded crows (C. (c.) cornix ) and tested whether extra‐pair paternity (EPP) was more likely in broods where parents differed in plumage colour. The proportion of broods with at least one extra‐pair offspring and the proportion of extra‐pair offspring were low overall (6.98% and 2.90%, respectively) with no evidence of hybrid broods having higher EPP rates than purebred nests.

The major histocompatibility complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assembled and annotated the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5 Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologues among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

The evolution of genetic sex determination is often accompanied by degradation of one of the proto sex chromosomes. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAfZWf, AAmZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf< 1) and dosage balance (0.5f:ZZmCorvus (corone) spp.) raised under common garden conditions. In addition, we characterized aspects of the regulatory landscape quantifying genome-wide ATAC-seq and 5mC methylation profiles. Partial dosage compensation was explained by female upregulation of Z-linked genes accompanied by increased chromatin accessibility on the female Z chromosome. 5mC methylation was strongly reduced in open chromatin-regions and GC islands and showed chromosome-, but no sex-specific variation. With the exception of the pseudo-autosomal region (PAR), female upregulation of gene expression was evenly spread across the Z chromosome without evidence for regional epigenetic regulation, as has for example been suggested for the male hypermethylated region (MHM) in chicken. Our results support the hypothesis that partial dosage compensation in female heterogametic systems is subject to chromosome-wide, epigenetic control mediated by differential chromatin accessibility between the sexes.

International audience; Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically, the probability of inbreeding mediated by mating system and/or population demography. Here, we investigated whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global data set of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstructed demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We found a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (1.5 Mb) were found in low latitude populations, and populations of known conservation concern. These include a Scottish killer whale, for which 37.8% of the autosomes were comprised of ROH >1.5 Mb in length. The fate of this population, in which only two adult males have been sighted in the past five years, and zero fecundity over the last two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.

The evolution of genetic barriers opposing interspecific gene flow is key to the origin of new species. Drawing from information on over 400 admixed genomes sourced from replicate transects across the European hybrid zone between all-black carrion crows and grey-coated hooded crows, we decipher the interplay between phenotypic divergence and selection at the molecular level. Over 68% of plumage variation was explained by epistasis between the gene NDP and a ~2.8-megabase region on chromosome 18 with suppressed recombination. Both pigmentation loci showed evidence for divergent selection resisting introgression. This study reveals how few, large-effect loci can govern prezygotic isolation and shield phenotypic divergence from gene flow.

Hybrid zones provide a window into the evolutionary processes governing species divergence. Yet, the contribution of mate choice to the temporal and spatial stability of hybrid zones remains poorly explored. Here, we investigate the effects of assortative mating on hybrid-zone dynamics by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and gray-coated hooded crows. In the best-fit model, narrow clines of the two mating-trait loci were maintained by a moderate degree of assortative mating inducing pre- and postzygotic isolation via positive frequency-dependent selection. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favoring gray-coated phenotypes ∼1 200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of the genetic mating-trait architecture for spatio-temporal hybrid-zone dynamics.

Adaptive divergence is the key evolutionary process generating biodiversity by means of natural selection. Yet, the conditions under which it can arise in the presence of gene flow remain contentious. To address this question, we subjected 132 sexually reproducing fission yeast populations sourced from two independent genetic backgrounds to disruptive ecological selection and manipulated the level of migration between environments. Contrary to theoretical expectations, adaptive divergence was most pronounced when migration was either absent (‘allopatry’) or maximal (‘sympatry’), but was much reduced at intermediate rates (‘parapatry’, ‘local mating’). This effect was apparent across central life history components (survival, asexual growth, and mating), but differed in magnitude between ancestral genetic backgrounds. The evolution of some fitness components was constrained by pervasive negative correlations (trade-off between asexual growth and mating), while others changed direction under the influence of migration (e.g. survival and mating). In allopatry, adaptive divergence was mainly conferred by standing genetic variation and resulted in ecological specialization. In sympatry, divergence was mainly mediated by novel mutations enriched in a subset of genes and was characterized by the repeated emergence of two strategies: an ecological generalist and an asexual growth specialist. Multiple loci showed consistent evidence for antagonistic pleiotropy across migration treatments and provide a conceptual link between adaptation and divergence. This evolve-and-resequence experiment demonstrates that rapid ecological differentiation can arise even under high rates of gene flow. It further highlights that adaptive trajectories are governed by complex interactions of gene flow, ancestral variation and genetic correlations.

The phylogeny and systematics of fur seals and sea lions (Otariidae) have long been studied with diverse data types, including an increasing amount of molecular data. However, only a few phylogenetic relationships have reached acceptance because of strong gene-tree species tree discordance. Divergence times estimates in the group also vary largely between studies. These uncertainties impeded the understanding of the biogeographical history of the group, such as when and how trans-equatorial dispersal and subsequent speciation events occurred. Here we used high-coverage genome-wide sequencing for 14 of the 15 species of Otariidae to elucidate the phylogeny of the family and its bearing on the taxonomy and biogeographical history. Despite extreme topological discordance among gene trees, we found a fully supported species tree that agrees with the few well-accepted relationships and establishes monophyly of the genusArctocephalus. Our data support a relatively recent trans-hemispheric dispersal at the base of a southern clade, which rapidly diversified into six major lineages between 3 to 2.5 Ma.Otariadiverged first, followed byPhocarctosand then four major lineages withinArctocephalus. However, we foundZalophusto be non-monophyletic, with California(Z. californianus)and Steller sea lions(Eumetopias jubatus)grouping closer than the Galapagos sea lion (Z. wollebaeki)with evidence for introgression between the two genera. Overall, the high degree of genealogical discordance was best explained by incomplete lineage sorting resulting from quasi-simultaneous speciation within the southern clade with introgresssion playing a subordinate role in explaining the incongruence among and within prior phylogenetic studies of the family.

Structural variation (SV) constitutes an important type of genetic mutations providing the raw material for evolution. Here, we uncover the genome-wide spectrum of intra- and interspecific SV segregating in natural populations of seven songbird species in the genus Corvus. Combining short-read (N = 127) and long-read re-sequencing (N = 31), as well as optical mapping (N = 16), we apply both assembly- and read mapping approaches to detect SV and characterize a total of 220,452 insertions, deletions and inversions. We exploit sampling across wide phylogenetic timescales to validate SV genotypes and assess the contribution of SV to evolutionary processes in an avian model of incipient speciation. We reveal an evolutionary young (~530,000 years) cis-acting 2.25-kb LTR retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth and evolutionary significance of SV segregating in natural populations and highlight the need for reliable SV genotyping., Structural genomic variation can fuel evolution. Here, authors present genome data from seven Corvus species and unearth structural variants that vary between incipient crow species in Europe, with implications for premating isolation involving plumage patterning.

Meiotic drivers have been proposed as a potent evolutionary force underlying genetic and phenotypic variation, genome structure, and also speciation. Due to their strong selective advantage, they are expected to rapidly spread through a population despite potentially detrimental effects on organismal fitness. Once fixed, autosomal drivers are cryptic within populations and only become visible in between‐population crosses lacking the driver or corresponding suppressor. However, the assumed ubiquity of meiotic drivers has rarely been assessed in crosses between populations or species. Here we test for meiotic drive in hybrid embryos and offspring of Timor and Australian zebra finches—subspecies that have evolved in isolation for about two million years—using 38,541 informative transmissions of 56 markers linked to either centromeres or distal chromosome ends. We did not find evidence for meiotic driver loci on specific chromosomes. However, we observed a weak overall transmission bias toward Timor alleles at centromeres in females (transmission probability of Australian alleles of 47%, nominal p = 6 × 10–5). While this is in line with the centromere drive theory, it goes against the expectation that the subspecies with the larger effective population size (i.e., the Australian zebra finch) should have evolved the more potent meiotic drivers. We thus caution against interpreting our finding as definite evidence for centromeric drive. Yet, weak centromeric meiotic drivers may be more common than generally anticipated and we encourage further studies that are designed to detect also small effect meiotic drivers., We test for meiotic drive in hybrid embryos and offspring of Timor and Australian zebra finches. We did not find evidence for meiotic driver loci on specific chromosomes but observed an overall transmission bias towards Timor alleles at centromeres in females.

Hybrid zones provide a window into the evolutionary processes governing species divergence. While the role of postzygotic isolation has been extensively characterized in the context of hybrid zones, the contribution of prezygotic isolation is less well explored. Here, we investigate the effects of assortative mate choice, the underlying preference function and mating-trait architecture, and the strength of sexual selection on hybrid zone dynamics. We explore this question by means of a mathematical model parameterized with phenotype and genotype data from the hybrid zone between all-black carrion and grey-coated hooded crows. The best-fit model resulted in narrow clines for two mating-trait loci coding for colour phenotype maintained by a moderate degree of assortative mating. Epistasis between the two loci induced hybrid-zone movement in favor of alleles conveying dark plumage followed by a shift in the opposite direction favouring grey-coated phenotypes ∼1,200 generations after secondary contact. Unlinked neutral loci diffused near-unimpeded across the zone. These results were generally robust to the choice of matching rule (self-referencing or parental imprinting) and effects of genetic drift. Overall, this study illustrates under which conditions assortative mating can maintain steep clines in mating-trait loci without generalizing to genome-wide reproductive isolation. It further emphasizes the importance of mating-trait architecture for spatio–temporal hybrid-zone dynamics.

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (Ne ) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage-specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in Ne . Utilizing population-level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in Ne would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population-scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (FST , PBS, dxy ) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.

Accurate and contiguous genome assembly is key to a comprehensive understanding of the processes shaping genomic diversity and evolution. Yet, it is frequently constrained by constitutive heterochromatin, usually characterized by highly repetitive DNA. As a key feature of genome architecture associated with centromeric and subtelomeric regions, it locally influences meiotic recombination. In this study, we assess the impact of large tandem repeat arrays on the recombination rate landscape in an avian speciation model, the Eurasian crow. We assembled two high-quality genome references using single-molecule real-time sequencing (long-read assembly [LR]) and single-molecule optical maps (optical map assembly [OM]). A three-way comparison including the published short-read assembly (SR) constructed for the same individual allowed assessing assembly properties and pinpointing misassemblies. By combining information from all three assemblies, we characterized 36 previously unidentified large repetitive regions in the proximity of sequence assembly breakpoints, the majority of which contained complex arrays of a 14-kb satellite repeat or its 1.2-kb subunit. Using whole-genome population resequencing data, we estimated the population-scaled recombination rate (ρ) and found it to be significantly reduced in these regions. These findings are consistent with an effect of low recombination in regions adjacent to centromeric or subtelomeric heterochromatin and add to our understanding of the processes generating widespread heterogeneity in genetic diversity and differentiation along the genome. By combining three different technologies, our results highlight the importance of adding a layer of information on genome structure that is inaccessible to each approach independently.

The ratio of nonsynonymous over synonymous sequence divergence, dN/dS, is a widely used estimate of the nonsynonymous over synonymous fixation rate ratio ω, which measures the extent to which natural selection modulates protein sequence evolution. Its computation is based on a phylogenetic approach and computes sequence divergence of protein-coding DNA between species, traditionally using a single representative DNA sequence per species. This approach ignores the presence of polymorphisms and relies on the indirect assumption that new mutations fix instantaneously, an assumption which is generally violated and reasonable only for distantly related species. The violation of the underlying assumption leads to a time-dependence of sequence divergence, and biased estimates of ω in particular for closely related species, where the contribution of ancestral and lineage-specific polymorphisms to sequence divergence is substantial. We here use a time-dependent Poisson random field model to derive an analytical expression of dN/dS as a function of divergence time and sample size. We then extend our framework to the estimation of the proportion of adaptive protein evolution α. This mathematical treatment enables us to show that the joint usage of polymorphism and divergence data can assist the inference of selection for closely related species. Moreover, our analytical results provide the basis for a protocol for the estimation of ω and α for closely related species. We illustrate the performance of this protocol by studying a population data set of four corvid species, which involves the estimation of ω and α at different time-scales and for several choices of sample sizes.

Over the past decade, the spotted wing Drosophila, Drosophila suzukii, has invaded Europe and America and has become a major agricultural pest in these areas, thereby prompting intense research activities to better understand its biology. Two draft genome assemblies based on short-read sequencing were released in 2013 for this species. Although valuable, these resources contain pervasive assembly errors and are highly fragmented, two features limiting their values. Our purpose here was to improve the assembly of the D. suzukii genome. For this, we generated PacBio long-read sequencing data at 160X sequence coverage and assembled a novel, contiguous D. suzukii genome. We obtained a high-quality assembly of 270 Mb (with 546 contigs, a N50 of 2.6Mb, a L50 of 15, and a BUSCO score of 95%) that we called WT3-2.0. We found that despite 16 rounds of full-sib crossings the D. suzukii strain that we sequenced has maintained high levels of polymorphism in some regions of its genome (ca. 19Mb). As a consequence, the quality of the assembly of these regions was reduced. We explored possible origins of this high residual diversity, including the presence of structural variants and a possible heterogeneous admixture pattern of North American and Asian ancestry. Overall, our WT3-2.0 assembly provides a higher quality genomic resource compared to the previous one in terms of general assembly statistics, sequence quality and gene annotation. This new D. suzukii genome assembly is therefore an improved resource for high-throughput sequencing approaches, as well as manipulative genetic technologies to study D. suzukii.

Funding was provided by the European Research Council (ERCStG-336536 FuncSpecGen to J.B.W.W.), the Swedish Research Council Vetenskapsrådet (621-2013-4510 to J.B.W.W.), the Knut and Alice Wallenberg Foundation (to J.B.W.W.), the Lawski foundation (to V.E.K. and J.B.W.W.), the German Research Foundation (KU 3402/1-1 to V.E.K.), the UK’s Biotechnology and Biological Sciences Research Council (BB/G023913/2 to C.R.), and the New Zealand Marsden Fund (to G.R.H.). Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species. Publisher PDF

Structural variation (SV) accounts for a substantial part of genetic mutations segregating across eukaryotic genomes with important medical and evolutionary implications. Here, we characterized SV across evolutionary time scales in the songbird genus Corvus using de novo assembly and read mapping approaches. Combining information from short-read (N = 127) and long-read re-sequencing data (N = 31) as well as from optical maps (N = 16) revealed a total of 201,738 insertions, deletions and inversions. Population genetic analysis of SV in the Eurasian crow speciation model revealed an evolutionary young (~530,000 years) cis-acting 2.25-kb retrotransposon insertion reducing expression of the NDP gene with consequences for premating isolation. Our results attest to the wealth of SV segregating in natural populations and demonstrate its evolutionary significance.

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

Funding: European Research Council (ERCStG-336536 FuncSpecGen to J.W.), the Swedish Research Council Vetenskapsrådet (621-2013-4510 to J.W.), Knut and Alice Wallenberg Foundation (to J.W.) and Tovetorp fieldstation through Stockholm University. Cooperation is a prevailing feature of many animal systems. Coalitionary aggression, where a group of individuals engages in coordinated behaviour to the detriment of conspecific targets, is a form of cooperation involving complex social interactions. To date, evidence has been dominated by studies in humans and other primates with a clear bias towards studies of male-male coalitions. We here characterize coalitionary aggression behaviour in a group of female carrion crows consisting of recruitment, coordinated chase, and attack. The individual of highest social rank liaised with the second most dominant individual to engage in coordinated chase and attack of a lower ranked crow on several occasions. Despite active intervention by the third most highly ranked individual opposing the offenders, the attack finally resulted in the death of the victim. All individuals were unrelated, of the same sex, and naive to the behaviour excluding kinship, reproduction, and social learning as possible drivers. Instead, the coalition may reflect a strategy of the dominant individual to secure long-term social benefits. Overall, the study provides evidence that members of the crow family engage in coordinated alliances directed against conspecifics as a possible means to manipulate their social environment. Publisher PDF

The effective size of a population (Ne), which determines its level of neutral variability, is a key evolutionary parameter. Ne can substantially depart from census sizes of present-day breeding populations (NC) as a result of past demographic changes, variation in life-history traits and selection at linked sites. Using genome-wide data we estimated the long-term coalescent Ne for 17 pinniped species represented by 36 population samples (total n = 458 individuals). Ne estimates ranged from 8,936 to 91,178, were highly consistent within (sub)species and showed a strong positive correlation with NC ( $${R}_{\mathrm{adj}}^2$$ = 0.59; P = 0.0002). Ne/NC ratios were low (mean, 0.31; median, 0.13) and co-varied strongly with demographic history and, to a lesser degree, with species’ ecological and life-history variables such as breeding habitat. Residual variation in Ne/NC, after controlling for past demographic fluctuations, contained information about recent population size changes during the Anthropocene. Specifically, species of conservation concern typically had positive residuals indicative of a smaller contemporary NC than would be expected from their long-term Ne. This study highlights the value of comparative population genomic analyses for gauging the evolutionary processes governing genetic variation in natural populations, and provides a framework for identifying populations deserving closer conservation attention. Using comparative population genomics across pinnipeds, this study explores how demographic change and life-history traits are correlated to the effective size of a population and conservation status.

The fission yeastSchizosaccharomyces pombeis an important model organism, but its natural diversity and evolutionary history remain under-studied. In particular, the population genomics of theS. pombemitochondrial genome (mitogenome) has not been thoroughly investigated. Here, we assembled the complete circular-mapping mitogenomes of 192S. pombeisolatesde novo, and found that these mitogenomes belong to 69 non-identical sequence types ranging from 17618 bp to 26910 bp in length. Using the assembled mitogenomes, we identified 20 errors in the reference mitogenome and discovered two previously unknown mitochondrial introns. Analysing sequence diversity of these 69 types of mitogenomes revealed two highly distinct clades, with only three mitogenomes exhibiting signs of inter-clade recombination. This diversity pattern suggests that currently availableS. pombeisolates descend from two long-separated ancestral lineages. This conclusion is corroborated by the diversity pattern of the recombination-repressedK-region located between donor mating-type locimat2andmat3in the nuclear genome. We estimated that the two ancestralS. pombelineages diverged about 31 million generations ago. These findings shed new light on the evolution ofS. pombeand the datasets generated in this study will facilitate future research on genome evolution.

Uncovering the genetic basis of species diversification is a central goal in evolutionary biology. Yet, the link between the accumulation of genomic changes during population divergence and the evolutionary forces promoting reproductive isolation is poorly understood. Here, we analysed 124 genomes of crow populations with various degrees of genome-wide differentiation, with parallelism of a sexually selected plumage phenotype, and ongoing hybridization. Overall, heterogeneity in genetic differentiation along the genome was best explained by linked selection exposed on a shared genome architecture. Superimposed on this common background, we identified genomic regions with signatures of selection specific to independent phenotypic contact zones. Candidate pigmentation genes with evidence for divergent selection were only partly shared, suggesting context-dependent selection on a multigenic trait architecture and parallelism by pathway rather than by repeated single-gene effects. This study provides insight into how various forms of selection shape genome-wide patterns of genomic differentiation as populations diverge., Genetic changes accumulate as populations diverge and new species emerge. Here, Vijay et al. resequence the genomes of crow populations at various stages of genetic differentiation, and shed light on the evolutionary processes acting during the origin and hybridization of crows in Europe and Asia.

Parallelism in genomic landscapes of differentiation, conserved genomic features and the role of linked selection

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global dataset of killer whale genomes in a rare attempt to elucidate global population structure in a non-human species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species’ range, likely associated with founder effects and allelic surfing during post-glacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

Mutation and recombination are key evolutionary processes governing phenotypic variation and reproductive isolation. We here demonstrate that biodiversity within all globally known strains of Schizosaccharomyces pombe arose through admixture between two divergent ancestral lineages. Initial hybridization was inferred to have occurred ∼20–60 sexual outcrossing generations ago consistent with recent, human-induced migration at the onset of intensified transcontinental trade. Species-wide heritable phenotypic variation was explained near-exclusively by strain-specific arrangements of alternating ancestry components with evidence for transgressive segregation. Reproductive compatibility between strains was likewise predicted by the degree of shared ancestry. To assess the genetic determinants of ancestry block distribution across the genome, we characterized the type, frequency, and position of structural genomic variation using nanopore and single-molecule real-time sequencing. Despite being associated with double-strand break initiation points, over 800 segregating structural variants exerted overall little influence on the introgression landscape or on reproductive compatibility between strains. In contrast, we found strong ancestry disequilibrium consistent with negative epistatic selection shaping genomic ancestry combinations during the course of hybridization. This study provides a detailed, experimentally tractable example that genomes of natural populations are mosaics reflecting different evolutionary histories. Exploiting genome-wide heterogeneity in the history of ancestral recombination and lineage-specific mutations sheds new light on the population history of S. pombe and highlights the importance of hybridization as a creative force in generating biodiversity.

Despite an increasing appreciation of the importance of host-microbe interactions in ecological and evolutionary processes, the factors shaping microbial communities in wild populations remain poorly understood. We therefore exploited a natural experiment provided by two adjacent Antarctic fur seal (Arctocephalus gazella) colonies of high and low social density and combined 16S rRNA metabarcoding with microsatellite profiling of mother-offspring pairs to investigate environmental and genetic influences on skin microbial communities. Seal-associated bacterial communities differed profoundly between the two colonies, despite the host populations themselves being genetically undifferentiated. Consistent with the hypothesis that social stress depresses bacterial diversity, we found that microbial alpha diversity was significantly lower in the high-density colony. Seals from one of the colonies that contained a stream also carried a subset of freshwater-associated bacteria, indicative of an influence of the physical environment. Furthermore, mothers and their offspring shared similar microbial communities, in support of the notion that microbes may facilitate mother-offspring recognition. Finally, a significant negative association was found between bacterial diversity and heterozygosity, a measure of host genetic quality. Our study thus uncovers a complex interplay between environmental and host genetic effects, while also providing empirical support for the leash model of host control, which posits that bacterial communities are driven not only by bottom-up species interactions, but also by top-down host regulation. Taken together, our findings have broad implications for understanding host-microbe interactions as well as prokaryotic diversity in general. This article is protected by copyright. All rights reserved. This article is protected by copyright. All rights reserved.

Functional validation of candidate genes for adaptation and speciation remains challenging. We here exemplify the utility of a method quantifying individual mRNA transcripts in revealing the molecular basis of divergence in feather pigment synthesis during early-stage speciation in crows. Using a padlock probe assay combined with rolling circle amplification we quantified cell-type specific gene expression in the native, histological context of growing feather follicles. Expression of Tyrosinase related protein 1 (TYRP1), Solute carrier family 45 member 2 (SLC45A2) and Hematopoietic prostaglandin D synthase (HPGDS) was melanocyte-limited and significantly reduced in follicles from hooded crow explaining the substantially lower melanin content in grey vs. black feathers. The central upstream transcription factor Microphthalmia-associated transcription factor (MITF) only showed differential expression specific to melanocytes - a feature not captured by bulk RNA-seq. Overall, this study provides insight into the molecular basis of an evolutionary young transition in pigment synthesis, and demonstrates the power of histologically explicit, statistically substantiated single-cell gene expression quantification for functional genetic inference in natural populations.

Recent advances in high throughput sequencing have transformed the study of wild organisms by facilitating the generation of high quality genome assemblies and dense genetic marker datasets. These resources have the potential to significantly advance our understanding of diverse phenomena at the level of species, populations and individuals, ranging from patterns of synteny through rates of linkage disequilibrium (LD) decay and population structure to individual inbreeding. Consequently, we used PacBio sequencing to refine an existing Antarctic fur seal (Arctocephalus gazella) genome assembly and genotyped 83 individuals from six populations using restriction site associated DNA (RAD) sequencing. The resulting hybrid genome comprised 6,169 scaffolds with an N50 of 6.21 Mb and provided clear evidence for the conservation of large chromosomal segments between the fur seal and dog (Canis lupus familiaris). Focusing on the most extensively sampled population of South Georgia, we found that LD decayed rapidly, reaching the background level ofr2= 0.09 by around 26 kb, consistent with other vertebrates but at odds with the notion that fur seals experienced a strong historical bottleneck. We also found evidence for population structuring, with four main Antarctic island groups being resolved. Finally, appreciable variance in individual inbreeding could be detected, reflecting the strong polygyny and site fidelity of the species. Overall, our study contributes important resources for future genomic studies of fur seals and other pinnipeds while also providing a clear example of how high throughput sequencing can generate diverse biological insights at multiple levels of organisation.

Animal coloration is one of the most conspicuous phenotypic traits in natural populations and has important implications for adaptation and speciation. Changes in coloration can occur over surprisingly short evolutionary timescales, while recurrence of similar colour patterns across large phylogenetic distances is also common. Even though the genetic basis of pigment production is well understood, little is known about the mechanisms regulating colour patterning. In this study, we shed light on the molecular elements regulating regional pigment production in two genetically near-identical crow taxa with striking differences in a eumelanin-based phenotype: black carrion and grey-coated hooded crows. We produced a high-quality genome annotation and analysed transcriptome data from a 2 × 2 design of active melanogenic feather follicles from head (black in both taxa) and torso (black in carrion and grey in hooded crow). Extensive, parallel expression differences between body regions in both taxa, enriched for melanogenesis genes (e.g. ASIP, CORIN, and ALDH6), indicated the presence of cryptic prepatterning also in all-black carrion crows. Meanwhile, colour-specific expression (grey vs. black) was limited to a small number of melanogenesis genes in close association with the central transcription factor MITF (most notably HPGDS, NDP and RASGRF1). We conclude that colour pattern differences between the taxa likely result from an interaction between divergence in upstream elements of the melanogenesis pathway and genes that provide an underlying prepattern across the body through positional information. A model of evolutionary stable prepatterns that can be exposed and masked through simple regulatory changes may explain the phylogenetically independent recurrence of colour patterns that is observed across corvids and many other vertebrate groups.

Females of many pinniped species generally exhibit strong fine-scale philopatry, but it is unclear over what spatial scale this behavior may translate into genetic population structure. We conducted a population genetic survey in the Galapagos fur seal, Arctocephalus galapagoensis, an endangered pinniped endemic to a small geographic range in the northwest of the Galapagos archipelago. To assess patterns of genetic diversity levels and population differentiation, we analyzed part of the mitochondrial control region (mtDNA) and 18 microsatellites DNA markers. We detected similar levels of genetic diversity to many other pinniped species (h = 0.86, pi = 0.012, A = 7.44) despite severe anthropogenic exploitation in the nineteenth century and recurrent population crashes due to recent climatic perturbations associated with El Nio Southern Oscillation events. We further found remarkably strong fine-scale matrilineal population structure, with 33.9 % of the mtDNA variation being partitioned among colonies separated by as little as 70 km swimming distance. In contrast, population structure inferred from nuclear markers was weak. Our findings provide further evidence that natal philopatry can translate into fine-scale genetic population structure in highly mobile species. We discuss the relevance of our results for the fine-scale conservation management of this species with a very restricted geographic range.

Parasites exert a major impact on the eco-evolutionary dynamics of their hosts and the associated biotic environment. Migration constitutes an effective means for long-distance invasions of vector-borne parasites and promotes their rapid spread. Yet, ecological and spatial information on population-specific host-parasite connectivity is essentially lacking. Here, we address this question in a system consisting of a transcontinental migrant species, the European barn swallow (Hirundo rustica) which serves as a vector for avian endoparasites in the genera Plasmodium, Haemoproteus and Leucocytozoon. Using feather stable isotope ratios as geographically informative markers, we first assessed migratory connectivity in the host: Northern European breeding populations predominantly overwintered in dry, savannah-like habitats in Southern Africa, whereas Southern European populations were associated with wetland habitats in Western Central Africa. Wintering areas of swallows breeding in Central Europe indicated a migratory divide with both migratory programmes occurring within the same breeding population. Subsequent genetic screens of parasites in the breeding populations revealed a link between the host's migratory programme and its parasitic repertoire: controlling for effects of local breeding location, prevalence of Africa-transmitted Plasmodium lineages was significantly higher in individuals overwintering in the moist habitats of Western Central Africa, even among sympatrically breeding individuals with different overwintering locations. For the rarer Haemoproteus parasites, prevalence was best explained by breeding location alone, whereas no clear pattern emerged for the least abundant parasite Leucocytozoon. These results have implications for our understanding of spatio-temporal host-parasite dynamics in migratory species and the spread of avian borne diseases.

A central paradigm in conservation biology is that population bottlenecks reduce genetic diversity and population viability. In an era of biodiversity loss and climate change, understanding the determinants and consequences of bottlenecks is therefore an important challenge. However, as most studies focus on single species, the multitude of potential drivers and the consequences of bottlenecks remain elusive. Here, we combined genetic data from over 11,000 individuals of 30 pinniped species with demographic, ecological and life history data to evaluate the consequences of commercial exploitation by 18th and 19th century sealers. We show that around one third of these species exhibit strong signatures of recent population declines. Bottleneck strength is associated with breeding habitat and mating system variation, and together with global abundance explains much of the variation in genetic diversity across species. Overall, bottleneck intensity is unrelated to IUCN status, although the three most heavily bottlenecked species are endangered. Our study reveals an unforeseen interplay between human exploitation, animal biology, demographic declines and genetic diversity., Historical hunting has caused documented declines in pinnipeds, but the extent to which hunting caused genetic bottlenecks among species was unknown. Here, the authors show evidence of severe bottlenecks in several pinniped species, particularly those that breed on land.