The Role of Hybridization in Species Formation and Persistence.

Hybridization, or interbreeding between different taxa, was traditionally considered to be rare and to have a largely detrimental impact on biodiversity, sometimes leading to the breakdown of reproductive isolation and even to the reversal of speciation. However, modern genomic and analytical methods have shown that hybridization is common in some of the most diverse clades across the tree of life, sometimes leading to rapid increase of phenotypic variability, to introgression of adaptive alleles, to the formation of hybrid species, and even to entire species radiations. In this review, we identify consensus among diverse research programs to show how the field has progressed. Hybridization is a multifaceted evolutionary process that can strongly influence species formation and facilitate adaptation and persistence of species in a rapidly changing world. Progress on testing this hypothesis will require cooperation among different subdisciplines.

Cycles of fusion and fission enabled rapid parallel adaptive radiations in African cichlids.

Although some lineages of animals and plants have made impressive adaptive radiations when provided with ecological opportunity, the propensities to radiate vary profoundly among lineages for unknown reasons. In Africa's Lake Victoria region, one cichlid lineage radiated in every lake, with the largest radiation taking place in a lake less than 16,000 years old. We show that all of its ecological guilds evolved in situ. Cycles of lineage fusion through admixture and lineage fission through speciation characterize the history of the radiation. It was jump-started when several swamp-dwelling refugial populations, each of which were of older hybrid descent, met in the newly forming lake, where they fused into a single population, resuspending old admixture variation. Each population contributed a different set of ancient alleles from which a new adaptive radiation assembled in record time, involving additional fusion-fission cycles. We argue that repeated fusion-fission cycles in the history of a lineage make adaptive radiation fast and predictable.

Repeated genetic adaptation to altitude in two tropical butterflies.

Repeated evolution can provide insight into the mechanisms that facilitate adaptation to novel or changing environments. Here we study adaptation to altitude in two tropical butterflies, Heliconius erato and H. melpomene, which have repeatedly and independently adapted to montane habitats on either side of the Andes. We sequenced 518 whole genomes from altitudinal transects and found many regions differentiated between highland (~ 1200 m) and lowland (~ 200 m) populations. We show repeated genetic differentiation across replicate populations within species, including allopatric comparisons. In contrast, there is little molecular parallelism between the two species. By sampling five close relatives, we find that a large proportion of divergent regions identified within species have arisen from standing variation and putative adaptive introgression from high-altitude specialist species. Taken together our study supports a role for both standing genetic variation and gene flow from independently adapted species in promoting parallel local adaptation to the environment.

Multispecies colour polymorphisms associated with contrasting microhabitats in two Mediterranean wrasse radiations.

Intraspecific colour polymorphisms (CPs) present unique opportunities to study fundamental evolutionary questions, such as the link between ecology and phenotype, mechanisms maintaining genetic diversity and their putative role in speciation. Wrasses are highly diverse in ecology and morphology and harbour a variety of colour-polymorphic species. In the Mediterranean Sea, wrasses of the tribe Labrini evolved two species radiations each harbouring several species with a brown and a green morph. The colour morphs occur in complete sympatry in mosaic habitats with rocky outcrops and Neptune grass patches. Morph-specific differences had not been characterized yet and the evolutionary forces maintaining them remained unknown. With genome-wide data for almost all Labrini species, we show that species with CPs are distributed across the phylogeny, but show evidence of hybridization. This suggests that the colour morphs are either ancient and have been lost repeatedly, that they have evolved repeatedly or have been shared via hybridization. Focusing on two polymorphic species, we find that each colour morph is more common in the microhabitat providing the best colour match and that the morphs exhibit additional behavioural and morphological differences further improving crypsis in their respective microhabitats. We find little evidence for genetic differentiation between the morphs in either species. Therefore, we propose that these colour morphs represent a multi-niche polymorphism as an adaptation to the highly heterogeneous habitat. Our study highlights how colour polymorphism (CP) can be advantageous in mosaic habitats and that Mediterranean wrasses are an ideal system to study trans-species polymorphisms, i.e. polymorphisms maintained across several species, in adaptive radiations.

Identification of a novel sex determining chromosome in cichlid fishes that acts as XY or ZW in different lineages.

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.

Genomics of altitude-associated wing shape in two tropical butterflies.

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 studied the genomic basis of wing shape in two Neotropical butterflies that inhabit large geographical ranges. Heliconius butterflies at high elevations have been shown to generally have rounder wings than those in the lowlands. We reared over 1,100 butterflies from 71 broods of H. erato and H. melpomene in common-garden conditions and showed that wing aspect ratio, that is, elongatedness, is highly heritable in both species and that elevation-associated wing aspect ratio differences are maintained. Genome-wide associations with a published data set of 666 whole genomes from across a hybrid zone, uncovered a highly polygenic basis to wing aspect ratio variation in the wild. We identified several genes that have roles in wing morphogenesis or wing aspect ratio variation in Drosophila flies, 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 for optix in H. erato. Thus, we present the first insights into the heritability and genomic basis of within-species wing aspect ratio in two Heliconius species, adding to a growing body of evidence that polygenic adaptation may underlie many ecologically relevant traits.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly 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.

Movement of transposable elements contributes to cichlid diversity.

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.

The ecological and genomic basis of explosive adaptive radiation.

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.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning.

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.

Rapid generation of ecologically relevant behavioral novelty in experimental cichlid hybrids.

The East African cichlid radiations are characterized by repeated and rapid diversification into many distinct species with different ecological specializations and by a history of hybridization events between nonsister species. Such hybridization might provide important fuel for adaptive radiation. Interspecific hybrids can have extreme trait values or novel trait combinations and such transgressive phenotypes may allow some hybrids to explore ecological niches neither of the parental species could tap into. Here, we investigate the potential of second-generation (F2) hybrids between two generalist cichlid species from Lake Malawi to exploit a resource neither parental species is specialized on: feeding by sifting sand. Some of the F2 hybrids phenotypically resembled fish of species that are specialized on sand sifting. We combined experimental behavioral and morphometric approaches to test whether the F2 hybrids are transgressive in both morphology and behavior related to sand sifting. We then performed a quantitative trait loci (QTL) analysis using RADseq markers to investigate the genetic architecture of morphological and behavioral traits. We show that transgression is present in several morphological traits, that novel trait combinations occur, and we observe transgressive trait values in sand sifting behavior in some of the F2 hybrids. Moreover, we find QTLs for morphology and for sand sifting behavior, suggesting the existence of some loci with moderate to large effects. We demonstrate that hybridization has the potential to rapidly generate novel and ecologically relevant phenotypes that may be suited to a niche neither of the parental species occupies. Interspecific hybridization may thereby contribute to the rapid generation of ecological diversity in cichlid radiations.

The Persistence of Polymorphisms across Species Radiations.

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.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths.

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to 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 invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. 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.

Microclimate buffering and thermal tolerance across elevations in a tropical butterfly.

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.

The coincidence of ecological opportunity with hybridization explains rapid adaptive radiation in Lake Mweru cichlid fishes.

The process of adaptive radiation was classically hypothesized to require isolation of a lineage from its source (no gene flow) and from related species (no competition). Alternatively, hybridization between species may generate genetic variation that facilitates adaptive radiation. Here we study haplochromine cichlid assemblages in two African Great Lakes to test these hypotheses. Greater biotic isolation (fewer lineages) predicts fewer constraints by competition and hence more ecological opportunity in Lake Bangweulu, whereas opportunity for hybridization predicts increased genetic potential in Lake Mweru. In Lake Bangweulu, we find no evidence for hybridization but also no adaptive radiation. We show that the Bangweulu lineages also colonized Lake Mweru, where they hybridized with Congolese lineages and then underwent multiple adaptive radiations that are strikingly complementary in ecology and morphology. Our data suggest that the presence of several related lineages does not necessarily prevent adaptive radiation, although it constrains the trajectories of morphological diversification. It might instead facilitate adaptive radiation when hybridization generates genetic variation, without which radiation may start much later, progress more slowly or never occur.

Eukaryote hybrid genomes.

Interspecific hybridization is the process where closely related species mate and produce offspring with admixed genomes. The genomic revolution has shown that hybridization is common, and that it may represent an important source of novel variation. Although most interspecific hybrids are sterile or less fit than their parents, some may survive and reproduce, enabling the transfer of adaptive variants across the species boundary, and even result in the formation of novel evolutionary lineages. There are two main variants of hybrid species genomes: allopolyploid, which have one full chromosome set from each parent species, and homoploid, which are a mosaic of the parent species genomes with no increase in chromosome number. The establishment of hybrid species requires the development of reproductive isolation against parental species. Allopolyploid species often have strong intrinsic reproductive barriers due to differences in chromosome number, and homoploid hybrids can become reproductively isolated from the parent species through assortment of genetic incompatibilities. However, both types of hybrids can become further reproductively isolated, gaining extrinsic isolation barriers, by exploiting novel ecological niches, relative to their parents. Hybrids represent the merging of divergent genomes and thus face problems arising from incompatible combinations of genes. Thus hybrid genomes are highly dynamic and undergo rapid evolutionary change, including genome stabilization in which selection against incompatible combinations results in fixation of compatible ancestry block combinations within the hybrid species. The potential for rapid adaptation or speciation makes hybrid genomes a particularly exciting subject of in evolutionary biology. Here we summarize how introgressed alleles or hybrid species can establish and how the resulting hybrid genomes evolve.

Altitude and life-history shape the evolution of Heliconius wings.

Phenotypic divergence between closely related species has long interested biologists. Taxa that inhabit a range of environments and have diverse natural histories can help understand how selection drives phenotypic divergence. In butterflies, wing color patterns have been extensively studied but diversity in wing shape and size is less well understood. Here, we assess the relative importance of phylogenetic relatedness, natural history, and habitat on shaping wing morphology in a large dataset of over 3500 individuals, representing 13 Heliconius species from across the Neotropics. We find that both larval and adult behavioral ecology correlate with patterns of wing sexual dimorphism and adult size. Species with solitary larvae have larger adult males, in contrast to gregarious Heliconius species, and indeed most Lepidoptera, where females are larger. Species in the pupal-mating clade are smaller than those in the adult-mating clade. Interestingly, we find that high-altitude species tend to have rounder wings and, in one of the two major Heliconius clades, are also bigger than their lowland relatives. Furthermore, within two widespread species, we find that high-altitude populations also have rounder wings. Thus, we reveal novel adaptive wing morphological divergence among Heliconius species beyond that imposed by natural selection on aposematic wing coloration.

A Combinatorial View on Speciation and Adaptive Radiation.

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.

Differential introgression of a female competitive trait in a hybrid zone between sex-role reversed species.

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.