Search

Sex chromosomes are evolutionarily labile in many animals, and fusion with an autosome is one of the mechanisms by which they can evolve. Sex chromosome-autosome (SA) fusions can reduce sexual conflict and promote adaptation and reproductive isolation among species. Recently, advances in genomics have fuelled the discovery of SA fusions in multiple vertebrates and in some invertebrates such as Lepidoptera. Here, we discovered multiple SA fusions in thesara/saphoclade of the classical adaptive radiation ofHeliconiusbutterflies. While the vast majority ofHeliconiusspecies have 21 chromosomes, species of the particularly rapidly diversifyingsaphoclade have up to 60 chromosomes. We found a sex chromosome fusion with chromosome 4 at the base of thesaphoclade and two additional sex chromosome fusions, each shared by two species. These sequential fusions between autosomes and sex chromosomes make theHeliconius saphoclade an ideal system to study the role of neo-sex chromosomes in adaptive radiations and the degeneration of sex chromosomes over time. Our study adds to a small but growing number of examples in butterflies with sex chromosome fusions that will help to unravel the importance of such rearrangements in the evolution of Lepidoptera and eukaryotes in general.Author SummaryChromosome number and structure are fundamental characteristics of a species. However, chromosomal rearrangements can occur spontaneously within a species, which can have significant consequences for recombination, segregation, speciation, and adaptation. Despite the importance of chromosomal rearrangements in evolution, we still have limited knowledge of the drivers and consequences of chromosomal fusions and fissions in natural populations, especially in adaptive radiations such asHeliconiusbutterflies. For the first time, we document chromosomal fusions between sex chromosomes and autosomes inHeliconius, particularly in species from thesara/saphoclade. Our research provides evidence for sex-autosome fusions involving autosomes 4, 9, and 14. All of these fusions seem to be associated with speciation events in this clade, with the sex-autosome 4 fusion being the oldest one. Although we do not yet understand the role or evolutionary consequences of these fusions, our study shows that chromosomal structure can evolve rapidly within a clade and generate chromosomal diversity. Overall, our findings contribute to a better understanding of chromosome evolution in highly diverse species.

Funder: The Branco Weiss Fellowship – Society in Science; doi: http://dx.doi.org/10.13039/501100001710, UNLABELLED: Sex determination systems are highly conserved among most vertebrates with genetic sex determination, but can be variable and evolve rapidly in some. Here, we study sex determination in a clade with exceptionally high sex chromosome turnover rates. We identify the sex determining chromosomes in three interspecific crosses of haplochromine cichlid fishes from Lakes Victoria and Malawi. We find evidence for different sex determiners in each cross. In the Malawi cross and one Victoria cross the same chromosome is sex-linked but while females are the heterogametic sex in the Malawi species, males are the heterogametic sex in the Victoria species. This chromosome has not previously been reported to be sex determining in cichlids, increasing the number of different chromosomes shown to be sex determining in cichlids to 12. All Lake Victoria species of our crosses are less than 15,000 years divergent, and we identified different sex determiners among them. Our study provides further evidence for the diversity and evolutionary flexibility of sex determination in cichlids, factors which might contribute to their rapid adaptive radiations. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10750-021-04560-7.

Speciation rates vary considerably among lineages, and our understanding of what drives the rapid succession of speciation events within young adaptive radiations remains incomplete1–11. The cichlid fish family provides a notable example of such variation, with many slowly speciating lineages as well as several exceptionally large and rapid radiations12. Here, by reconstructing a large phylogeny of all currently described cichlid species, we show that explosive speciation is solely concentrated in species flocks of several large young lakes. Increases in the speciation rate are associated with the absence of top predators; however, this does not sufficiently explain explosive speciation. Across lake radiations, we observe a positive relationship between the speciation rate and enrichment of large insertion or deletion polymorphisms. Assembly of 100 cichlid genomes within the most rapidly speciating cichlid radiation, which is found in Lake Victoria, reveals exceptional ‘genomic potential’—hundreds of ancient haplotypes bear insertion or deletion polymorphisms, many of which are associated with specific ecologies and shared with ecologically similar species from other older radiations elsewhere in Africa. Network analysis reveals fundamentally non-treelike evolution through recombining old haplotypes, and the origins of ecological guilds are concentrated early in the radiation. Our results suggest that the combination of ecological opportunity, sexual selection and exceptional genomic potential is the key to understanding explosive adaptive radiation. Analyses of the genomes of cichlid species reveal that the combination of ecological opportunity, sexual selection and exceptional genomic potential is the key to understanding explosive adaptive radiation in cichlids.

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to high altitude in two divergent tropical butterflies, H. erato and H. melpomene, which have repeatedly and independently adapted to high elevations on either side of the Andean mountains. We sequenced 518 whole genomes from elevational transects and found many regions under selection at high altitude, with repeated genetic differentiation across multiple replicates, including allopatric comparisons. In contrast, there is little ‘molecular parallelism’ between H. erato and H. melpomene. With a further 85 whole genomes of five close relatives, we find that a large proportion divergent regions have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a key role of standing genetic variation and gene flow from pre-adapted species in promoting parallel genetic local adaptation to the environment.

Significance A defining goal in genetics is linking variation in DNA sequence to trait evolution between populations and, ultimately, species. Genome sequencing efficiently captures such variation but typically in millions of tiny fragments that omit haplotype or linkage information. We present “haplotagging,” a simple, rapid linked-read sequencing technique that allows high-throughput sequencing without sacrificing haplotype information. We validated this affordable approach for whole-genome haplotyping in large populations. We used haplotagging to investigate the rise of a novel hybrid morph in parallel hybrid zones of two comimetic Heliconius butterfly species in Ecuador. Our results reveal that strikingly parallel divergences in their genomes produced coordinated shifts in haplotype frequencies across the hybrid zone, giving rise to comimetic hybrid morphs in each species., Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present “haplotagging,” a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.

Speciation is often thought of as a slow process due to the waiting times for mutations that cause incompatibilities, and permit ecological differentiation or assortative mating. Cases of rapid speciation and particularly cases of rapid adaptive radiation into multiple sympatric species have remained somewhat mysterious. We review recent findings from speciation genomics that reveal an emerging commonality among such cases: reassembly of old genetic variation into new combinations facilitating rapid speciation and adaptive radiation. The polymorphisms in old variants frequently originated from hybridization at some point in the past. We discuss why old variants are particularly good fuel for rapid speciation, and hypothesize that variation in access to such old variants might contribute to the large variation in speciation rates observed in nature.

Understanding how organisms adapt to their local environment is central to evolution. With new whole-genome sequencing technologies and the explosion of data, deciphering the genomic basis of complex traits that are ecologically relevant is becoming increasingly feasible. Here we study the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges.Heliconiusbutterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1100 butterflies from 71 broods ofH. eratoandH. melpomenein common-garden conditions and show that wing aspect ratio, i.e. elongatedness, is highly heritable in both species and elevation-associated wing shape differences are maintained. Genome-wide associations with a published dataset of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing shape variation in the wild. We identify several genes that have roles in wing morphogenesis or wing shape variation inDrosophilaflies, making them promising candidates for future studies. There was little evidence for molecular parallelism in the two species, with only one shared candidate gene, nor for a role of the four known colour pattern loci, except foroptixinH. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing shape in twoHeliconiusspecies, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Background Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism. Findings We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity. Conclusions We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.

Studies on polymorphisms have been foundational to our understanding of evolution. The presence of different phenotypic morphs is sometimes considered a precursor to speciation in which morphs evolve into different species. While speciation should initially reduce genetic variation in daughter versus parental species, a common pattern is the recurrence of the same phenotypic polymorphism in many species of a clade. Despite the ubiquity of these persistent polymorphisms, there is little discussion of their evolutionary origins. How does the genetic variation underpinning such polymorphisms cross speciation boundaries? What selection pressures maintain the morphs in multiple daughter species? Using diverse case studies, we highlight the characteristics of polymorphisms and selection regimes influencing the likelihood of polymorphism retention across species radiations.

Microclimatic variability in tropical forests plays a key role in shaping species distributions and their ability to cope with environmental change, especially for ectotherms. Nonetheless, currently available climatic datasets lack data from the forest interior and, furthermore, our knowledge of thermal tolerance among tropical ectotherms is limited. We therefore studied natural variation in the microclimate experienced by tropical butterflies in the genus Heliconius across their Andean range in a single year. We found that the forest strongly buffers temperature and humidity in the understorey, especially in the lowlands, where temperatures are more extreme. There were systematic differences between our yearly records and macroclimate databases (WorldClim2), with lower interpolated minimum temperatures and maximum temperatures higher than expected. We then assessed thermal tolerance of 10 Heliconius butterfly species in the wild and found that populations at high elevations had significantly lower heat tolerance than those at lower elevations. However, when we reared populations of the widespread H. erato from high and low elevations in a common-garden environment, the difference in heat tolerance across elevations was reduced, indicating plasticity in this trait. Microclimate buffering is not currently captured in publicly available datasets, but could be crucial for enabling upland shifting of species sensitive to heat such as highland Heliconius. Plasticity in thermal tolerance may alleviate the effects of global warming on some widespread ectotherm species, but more research is needed to understand the long-term consequences of plasticity on populations and species., Highlighted Article: Tropical forests along the Andes can greatly buffer climate. Butterflies inhabiting high elevations are less thermally tolerant but common-garden experiments suggest this is largely due to phenotypic plasticity.

Hybridization between invasive and native species has raised global concern given the dramatic increase in species range shifts and pest outbreaks due to climate change, development of suitable agroecosystems, and anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasiveHelicoverpa armigeraandH. zeain Brazil. We carried out whole-genome resequencing ofHelicoverpamoths from Brazil in two temporal samples: during the outbreak ofH. armigerain 2013, and more recent populations from 2017. There is evidence for a burst of hybridization and widespread introgression from localH. zeainto invasiveH. armigeracoinciding withH. armigeraexpansion in 2013. However, inH. armigera, admixture proportions were reduced between 2013 and 2017, indicating a decline in hybridization rates. Recent populations also showed shorter introgressed tracks suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for introgression of a single region including an insecticide-resistance allele from the invasiveH. armigerainto localH. zea,which increased in frequency over time but was localized within the genome. In summary, despite extensive gene-flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression and hybridization, with significant implications for pest management in future population expansions and introductions of novel resistance genes from new invasiveH. armigerapopulations.Author summarySecondary contact occurs when related species with non-overlapping ranges are geographically reunited. Scenarios of secondary contact have increased due to anthropogenic movement of species outside of their native range, often resulting in invasive species that successfully spread and stabilised in the new environment. This is the case forHelicoverpa armigera, a major agricultural pest in the Old World that has recently invaded the Americas, where it reunited with its closest relative,H. zea. While some authors reported hybridisation, and hypothesised about the potential emergence of novel ecotypes and the exchange of pesticide-resistant genes, these outcomes have not been tested yet. We examine these outcomes by sequencing individuals from both species in Brazil, collected in 2013 after outbreaks ofH. armigerawere reported, and individuals collected during 2017. We discovered that despite hybridisation, these moths have not collapsed into a single species nor formed new ecotypes, and that the species distinctiveness is maintained through selection against most of the foreign genotypes that cross species boundaries. However, we found that hybridisation mediated the rapid acquisition of aH. armigeragene conferring resistance to pyrethroids byH. zea. The overall decline in populations of both species during the interval covered by this study means that our results are likely to reflect the consequences of hybridization events early after invasion, despite the likely ongoing introduction ofH. armigeragenetic diversity through trade across the South American continent. Our results provide a rare example of adaptive transferral of variation right after invasion and elucidate the dynamics of insecticide resistance evolution inH. zea.

Phenotypic divergence between closely related species has long interested biologists. Taxa that inhabit a range of environments and have known and diverse, natural histories, can help understand how different selection pressures shape diverging traits. In butterflies, wing colour patterns have been extensively studied, whereas wing shape diversity is less well understood despite its importance for flight. Here we study a measure of wing shape, aspect ratio, and wing size in a large dataset of over 3200 individuals, representing 13Heliconiusspecies, from across the Neotropics. We assess the relative importance of phylogenetic relatedness, natural history and habitat in determining wing shape and size. We find that both larval and adult behavioural ecology affect patterns of adult size dimorphisms. On one hand, species with solitary larvae have larger adult males, in contrast to gregariousHeliconiusspecies, and indeed most Lepidoptera. On the other hand, species in the pupal-mating clade are smaller overall than those in the adult-mating clade. Interestingly, while controlling for phylogeny, sex ratios and allometry, we find that species inhabiting higher altitudes have rounder wings and, in one of the two majorHeliconiusclades, are also bigger in size than their lowland relatives. Thus, we reveal novel adaptive wing morphological divergence amongHeliconiusspecies beyond that imposed by natural selection on aposematic wing colouration. Our study highlights the value of phylogenetic comparative methods in study systems that have diverse and well-studied natural histories to disentangle the selection pressures shaping adaptive phenotypes.

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male-male competition in maintaining or eroding species barriers. Although female-female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex-role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female-female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex-role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.

Adaptive radiation research typically relies on the study of evolution in retrospective, leaving the predictive value of the concept hard to evaluate. Several radiations, including the cichlid fishes in the East African Great Lakes, have been studied extensively, yet no study has investigated the onset of the intraspecific processes of niche expansion and differentiation shortly after colonization of an adaptive zone by cichlids. Haplochromine cichlids of one of the two lineages that seeded the Lake Victoria radiation recently arrived in Lake Chala, a lake perfectly suited for within-lake cichlid speciation. Here, we infer the colonization and demographic history, quantify phenotypic, ecological and genomic diversity and diversification, and investigate the selection regime to ask if the population shows signs of diversification resembling the onset of adaptive radiation. We find that since their arrival in the lake, haplochromines have colonized a wide range of depth habitats associated with ecological and morphological expansion and the beginning of phenotypic differentiation and potentially nascent speciation, consistent with the very early onset of an adaptive radiation process. Moreover, we demonstrate evidence of rugged phenotypic fitness surfaces, indicating that current ecological selection may contribute to the phenotypic diversification.

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

The genetic basis of parallel evolution of similar species is of great interest in evolutionary biology. In the adaptive radiation of Lake Victoria cichlid fishes, sister species with either blue or red-back male nuptial coloration have evolved repeatedly, often associated with shallower and deeper water, respectively. One such case is blue and red-backed Pundamilia species, for which we recently showed that a young species pair may have evolved through “hybrid parallel speciation”. Coalescent simulations suggested that the older species P. pundamilia (blue) and P. nyererei (red-back) admixed in the Mwanza Gulf and that new “nyererei-like” and “pundamilia-like” species evolved from the admixed population. Here, we use genome scans to study the genomic architecture of differentiation, and assess the influence of hybridization on the evolution of the younger species pair. For each of the two species pairs, we find over 300 genomic regions, widespread across the genome, which are highly differentiated. A subset of the most strongly differentiated regions of the older pair are also differentiated in the younger pair. These shared differentiated regions often show parallel allele frequency differences, consistent with the hypothesis that admixture-derived alleles were targeted by divergent selection in the hybrid population. However, two-thirds of the genomic regions that are highly differentiated between the younger species are not highly differentiated between the older species, suggesting independent evolutionary responses to selection pressures. Our analyses reveal how divergent selection on admixture-derived genetic variation can facilitate new speciation events.

Adaptive radiation research typically relies on the study of evolution in retrospective, leaving the predictive value of the concept hard to evaluate. Several radiations, including the cichlid fishes in the East African Great Lakes, have been studied extensively, yet no study has investigated the onset of the intraspecific processes of niche expansion and differentiation shortly after colonization of an adaptive zone by cichlids. Haplochromine cichlids of one of the two lineages that seeded the Lake Victoria radiation recently arrived in Lake Chala, a lake perfectly suited for within-lake cichlid speciation. Here, we infer the colonization and demographic history, quantify phenotypic, ecological and genomic diversity and diversification, and investigate the selection regime to ask if the population shows signs of diversification resembling the onset of adaptive radiation. We find that since their arrival in the lake, haplochromines have colonized a wide range of depth habitats associated with ecological and morphological expansion and the beginning of phenotypic differentiation and potentially nascent speciation, consistent with the very early onset of an adaptive radiation process. Moreover, we demonstrate evidence of rugged phenotypic fitness surfaces, indicating that current ecological selection may contribute to the phenotypic diversification.

Genetic linkage maps are essential for comparative genomics, high quality genome sequence assembly and fine scale quantitative trait locus (QTL) mapping. In the present study we identified and genotyped markers via restriction-site associated DNA (RAD) sequencing and constructed a genetic linkage map based on 1,597 SNP markers of an interspecific F2 cross of two closely related Lake Victoria cichlids (Pundamilia pundamilia and P. sp. “red head”). The SNP markers were distributed on 22 linkage groups and the total map size was 1,594 cM with an average marker distance of 1.01 cM. This high-resolution genetic linkage map was used to anchor the scaffolds of the Pundamilia genome and estimate recombination rates along the genome. Via QTL mapping we identified a major QTL for sex in a ∼1.9 Mb region on Pun-LG10, which is homologous to Oreochromis niloticus LG 23 (Ore-LG23) and includes a well-known vertebrate sex-determination gene (amh).

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

Modes and mechanisms of speciation are best studied in young species pairs. In older taxa, it is increasingly difficult to distinguish what happened during speciation from what happened after speciation. Lake Victoria cichlids in the genus Pundamilia encompass a complex of young species and polymorphic populations. One Pundamilia species pair, P. pundamilia and P. nyererei, is particularly well suited to study speciation because sympatric population pairs occur with different levels of phenotypic differentiation and reproductive isolation at different rocky islands within the lake. Genetic distances between allopatric island populations of the same nominal species often exceed those between the sympatric species. It thus remained unresolved whether speciation into P. nyererei and P. pundamilia occurred once, followed by geographical range expansion and interspecific gene flow in local sympatry, or if the species pair arose repeatedly by parallel speciation. Here, we use genomic data and demographic modelling to test these alternative evolutionary scenarios. We demonstrate that gene flow plays a strong role in shaping the observed patterns of genetic similarity, including both gene flow between sympatric species and gene flow between allopatric populations, as well as recent and early gene flow. The best supported model for the origin of P. pundamilia and P. nyererei population pairs at two different islands is one where speciation happened twice, whereby the second speciation event follows shortly after introgression from an allopatric P. nyererei population that arose earlier. Our findings support the hypothesis that very similar species may arise repeatedly, potentially facilitated by introgressed genetic variation.

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this., Author Summary Ecological speciation can be defined as the evolution of new, reproductively isolated, species driven by natural selection and ecologically-mediated sexual selection. Its genomic signature has mainly been studied in ecotypes and emerging species that started diverging hundreds to thousands of generations ago, while little is known about the very early stages of species divergence. To fill this knowledge gap, we studied whether and how threespine stickleback, which have adapted either to lake or to stream environments in less than 150 years, differ across their genomes. We found several segments of the genome to be clearly divergent between lake and stream ecotypes, even when both forms breed side by side in the same area. Strikingly, this genomic differentiation was mainly concentrated in one region with low to intermediate recombination rates and clustered around genes controlling ecotype-specific phenotypic traits. Our findings suggest that genomic differentiation can arise despite gene flow already very early at the onset of speciation, and that its occurrence may be facilitated by the genomic organization of genes that control traits involved in adaptation and reproductive isolation.

Ecological speciation is the evolution of reproductive isolation as a consequence of direct divergent natural selection or ecologically mediated divergent sexual selection. While the genomic signature of the former has been extensively studied in recent years, only few examples exist for genomic differentiation where environment-dependent sexual selection has played an important role. Here, we describe a very young (~90 years old) population of threespine sticklebacks exhibiting phenotypic and genomic differentiation between two habitats within the same pond. We show that differentiation among habitats is limited to male throat colour and nest type, traits known to be subject to sexual selection. Divergence in these traits mirrors divergence in much older benthic and limnetic stickleback species pairs from North American west coast lakes, which also occur in sympatry but are strongly reproductively isolated from each other. We demonstrate that in our population, differences in throat colour and breeding have been stable over a decade, but in contrast to North American benthic and limnetic stickleback species, these mating trait differences are not accompanied by divergence in morphology related to feeding, predator defence or swimming performance. Using genomewide SNP data, we find multiple genomic islands with moderate differentiation spread across several chromosomes, whereas the rest of the genome is undifferentiated. The islands contain potential candidate genes involved in visual perception of colour. Our results suggest that phenotypic and multichromosome genomic divergence of these morphs was driven by environment-dependent sexual selection, demonstrating incipient speciation after only a few decades of divergence in sympatry.