Reproductive Isolation Research Articles

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse.

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Beyond color and pattern: elucidating the factors associated with intraspecific aggression in the mimic poison frog (Ranitomeya imitator)

Assortative or disassortative aggression – when individuals display more aggression towards conspecifics with similar or different phenotypic characteristics – can either maintain polymorphisms or facilitate gene flow between populations depending on which direction the aggression is aimed. Deciphering which factors elicit or prevent aggression is crucial to improving our knowledge of the origin and maintenance of reproductive barriers and subsequent speciation. The Peruvian mimic poison frog, Ranitomeya imitator, is a monogamous and territorial species that has evolved into four distinct color-pattern morphs in a mimetic radiation. Here we use historical landscape genetic data and competition trials between male individuals sourced from different populations and color-pattern morphs to show that the level of aggression between individuals is not associated with color morph or body size but rather with source population. Individuals spent more time in combat with individuals from their own deme (genetically homogeneous population), irrespective of color morph or size. These findings indicate that genotypic similarity is correlated with increased aggression in R. imitator, though the mechanism by which R. imitator males identify conspecifics as territorial threats remains unclear. As body size and color morph were not significantly associated with aggression levels, this study emphasizes the necessity of further research to identify whether other phenotypic traits are influencing territorial behavior between male frogs, and if these factors play a role in increasing gene flow, or conversely, the formation of reproductive barriers between populations.

Genetic drift promotes and recombination hinders speciation on holey fitness landscapes.

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

Divergent warning patterns influence male and female mating behaviours in a tropical butterfly.

Traits under divergent ecological selection that also function during mating can be important in maintaining species boundaries. Few studies have considered mutual mate choice, where both males and females base mating decisions on the same trait. Wing colouration in Heliconius butterflies evolved as a warning signal but also functions as a mating cue. We investigated the contribution of visual preference to assortative mating in an aposematic butterfly Heliconius cydno in the context of reproductive isolation with its sympatric, visually distinct relative Heliconius melpomene. Heliconius cydno have conspicuous white bands on their forewings, whereas those of H. melpomene are red in colour. We predicted that both sexes of H. cydno contributed to assortative mating by exhibiting visual preference towards conspecific wing colouration. We analysed published and new data from preference experiments, in which males were presented with conspecific and H. melpomene females. We also recorded female responses and mating outcomes in choice experiments, involving conspecific males with either the original white or artificially painted red forewing bands. Both sexes of H. cydno responded more positively towards the conspecific colouration, and males strongly preferred females of its own colours. In contrast, male colouration did not predict mating outcomes in female choice experiments. As courtships are initiated by males in butterflies, our findings suggest that female visual preference might be of secondary importance in H. cydno. Our data also suggest that the contribution of visual preference to reproductive isolation might be unequal between H. cydno and its sympatric relative H. melpomene.

Genetic Cause of Hybrid Lethality Observed in Reciprocal Interspecific Crosses between Nicotiana simulans and N. tabacum.

Hybrid lethality, a type of postzygotic reproductive isolation, is an obstacle to wide hybridization breeding. Here, we report the hybrid lethality that was observed in crosses between the cultivated tobacco, Nicotiana tabacum (section Nicotiana), and the wild tobacco species, Nicotiana simulans (section Suaveolentes). Reciprocal hybrid seedlings were inviable at 28 °C, and the lethality was characterized by browning of the hypocotyl and roots, suggesting that hybrid lethality is due to the interaction of nuclear genomes derived from each parental species, and not to a cytoplasmic effect. Hybrid lethality was temperature-sensitive and suppressed at 36 °C. However, when hybrid seedlings cultured at 36 °C were transferred to 28 °C, all of them showed hybrid lethality. After crossing between an N. tabacum monosomic line missing one copy of the Q chromosome and N. simulans, hybrid seedlings with or without the Q chromosome were inviable and viable, respectively. These results indicated that gene(s) on the Q chromosome are responsible for hybrid lethality and also suggested that N. simulans has the same allele at the Hybrid Lethality A1 (HLA1) locus responsible for hybrid lethality as other species in the section Suaveolentes. Haplotype analysis around the HLA1 locus suggested that there are at least six and two haplotypes containing Hla1-1 and hla1-2 alleles, respectively, in the section Suaveolentes.

Reproductive isolation in a three-way contact zone.

Contact zones between divergent forms within a species provide insight into the role of gene flow in adaptation and speciation. Previous work has focused on contact zones involving only two divergent forms, but in nature, many more than two populations may overlap simultaneously and experience gene flow. Patterns of introgression in wild populations are, therefore, likely much more complicated than is often assumed. We begin to address this gap in current knowledge by investigating patterns of divergence and introgression across a complex natural contact zone. We use phenotypic and genomic data to confirm the existence of a three-way contact zone among divergent freshwater resident, saltwater resident and saltwater migratory three-spined stickleback (Gasterosteus aculeatus) on the island of North Uist, Scottish Western Isles. We find evidence for hybridization, mostly between saltwater resident and saltwater migratory forms. Despite hybridization, genomic analyses reveal pairwise islands of divergence between all forms that are maintained across the contact zone. Genomic cline analyses also provide evidence for selection and/or hybrid incompatibilities in divergent regions. Divergent genomic regions occur across multiple chromosomes and involve many known adaptive loci and several chromosomal inversions. We also identify distinct immune gene expression profiles between forms, but no evidence for transgressive expression in hybrids. Our results suggest that reproductive isolation is maintained in this three-way contact zone, despite some hybridization, and that reduced recombination in chromosomal inversions may play an important role in maintaining this isolation.

Natural variation in the zinc-finger-encoding exon of Prdm9 affects hybrid sterility phenotypes in mice.

PRDM9-mediated reproductive isolation was first described in the progeny of Mus musculus musculus (MUS) PWD/Ph and Mus musculus domesticus (DOM) C57BL/6J inbred strains. These male F1 hybrids fail to complete chromosome synapsis and arrest meiosis at prophase I, due to incompatibilities between the Prdm9 gene and hybrid sterility locus Hstx2. We identified 14 alleles of Prdm9 in exon 12, encoding the DNA-binding domain of the PRDM9 protein in outcrossed wild mouse populations from Europe, Asia, and the Middle East, 8 of which are novel. The same allele was found in all mice bearing introgressed t-haplotypes encompassing Prdm9. We asked whether 7 novel Prdm9 alleles in MUS populations and the t-haplotype allele in 1 MUS and 3 DOM populations induce Prdm9-mediated reproductive isolation. The results show that only combinations of the dom2 allele of DOM origin and the MUS msc1 allele ensure complete infertility of intersubspecific hybrids in outcrossed wild populations and inbred mouse strains examined so far. The results further indicate that MUS mice may share the erasure of PRDM9msc1 binding motifs in populations with different Prdm9 alleles, which implies that erased PRDM9 binding motifs may be uncoupled from their corresponding Prdm9 alleles at the population level. Our data corroborate the model of Prdm9-mediated hybrid sterility beyond inbred strains of mice and suggest that sterility alleles of Prdm9 may be rare.

A lethal mitonuclear incompatibility in complex I of natural hybrids

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1–3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4–6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Evidence that genetic drift not adaptation drives fast-Z and large-Z effects in Ficedula flycatchers.

The sex chromosomes have been hypothesized to play a key role in driving adaptation and speciation across many taxa. The reason for this is thought to be the hemizygosity of the heteromorphic part of sex chromosomes in the heterogametic sex, which exposes recessive mutations to natural and sexual selection. The exposure of recessive beneficial mutations increases their rate of fixation on the sex chromosomes, which results in a faster rate of evolution. In addition, genetic incompatibilities between sex-linked loci are exposed faster in the genomic background of hybrids of divergent lineages, which makes sex chromosomes contribute disproportionately to reproductive isolation. However, in birds, which show a Z/W sex determination system, the role of adaptation versus genetic drift as the driving force of the faster differentiation of the Z chromosome (fast-Z effect) and the disproportionate role of the Z chromosome in reproductive isolation (large-Z effect) are still debated. Here, we address this debate in the bird genus Ficedula flycatchers based on population-level whole-genome sequencing data of six species. Our analysis provides evidence for both faster lineage sorting and reduced gene flow on the Z chromosome than the autosomes. However, these patterns appear to be driven primarily by the increased role of genetic drift on the Z chromosome, rather than an increased rate of adaptive evolution. Genomic scans of selective sweeps and fixed differences in fact suggest a reduced action of positive selection on the Z chromosome.

The role of temporal reproductive isolation, trophic polymorphism and growth rate fluctuations in the diversification of Arctic charr Salvelinus alpinus (L.) in Lake Kalarskii Davatchan, Transbaikalia, Russia

The role of temporal reproductive isolation, trophic polymorphism and growth rate fluctuations in the diversification of Arctic charr Salvelinus alpinus (L.) in Lake Kalarskii Davatchan, Transbaikalia, Russia

Coupling of Barriers to Gene Exchange: Causes and Consequences.

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three ‘pure’ species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an ‘evolutionary syngameon’ where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Pronounced differentiation on the Z chromosome and parts of the autosomes in crowned sparrows contrasts with mitochondrial paraphyly: implications for speciation.

When a single species evolves into multiple descendent species, some parts of the genome can play a key role in the evolution of reproductive isolation while other parts flow between the evolving species via interbreeding. Genomic evolution during the speciation process is particularly interesting when major components of the genome-for instance, sex chromosomes vs. autosomes vs. mitochondrial DNA-show widely differing patterns of relationships between three diverging populations. The golden-crowned sparrow (Zonotrichia atricapilla) and the white-crowned sparrow (Zonotrichia leucophrys) are phenotypically differentiated sister species that are largely reproductively isolated despite possessing similar mitochondrial genomes, likely due to recent introgression. We assessed variation in more than 45,000 single nucleotide polymorphisms to determine the structure of nuclear genomic differentiation between these species and between two hybridizing subspecies of Z. leucophrys. The two Z. leucophrys subspecies show moderate levels of relative differentiation and patterns consistent with a history of recurrent selection in both ancestral and daughter populations, with much of the sex chromosome Z and a large region on the autosome 1A showing increased differentiation compared to the rest of the genome. The two species Z. leucophrys and Z. atricapilla show high relative differentiation and strong heterogeneity in the level of differentiation among various chromosomal regions, with a large portion of the sex chromosome (Z) showing highly divergent haplotypes between these species. Studies of speciation often emphasize mitochondrial DNA differentiation, but speciation between Z. atricapilla and Z. leucophrys appears primarily associated with Z chromosome divergence and more moderately associated with autosomal differentiation, whereas mitochondria are highly similar due apparently to recent introgression. These results add to the growing body of evidence for highly heterogeneous patterns of genomic differentiation during speciation, with some genomic regions showing a lack of gene flow between populations many hundreds of thousands of years before other genomic regions.

The role of divergent host use and geography in the evolution of habitat isolation and sexual isolation among sister species of Belonocnema gall wasps.

Ecology and geography can play important roles in the evolution of reproductive isolation across the speciation continuum, but few studies address both at the later stages of speciation. This notable gap in knowledge arises from the fact that traditional ecological speciation studies have predominantly focused on the role of ecology in initiating the speciation process, while many studies exploring the effect of geography (e.g., reinforcement) concentrate on species pairs that lack divergent ecological characteristics. We simultaneously examine the strength of habitat isolation and sexual isolation among three closely related species of Belonocnema gall-forming wasps on two species of live oaks, Quercus virginiana and Q. geminata, that experience divergent selection from their host plants and variable rates of migration due to their geographic context. We find that the strength of both habitat isolation and sexual isolation is lowest among allopatric species pairs with the same host plant association, followed by allopatric species with different host plant associations, and highest between sympatric species with different host-plant associations. This pattern suggests that divergent selection due to different host use interacts with geography in the evolution of habitat isolation and sexual isolation during the later stages of speciation of Belonocnema wasps.

Postzygotic barriers persist despite ongoing introgression in hybridizing Mimulus species.

The evolution of postzygotic isolation is thought to be a key step in maintaining species boundaries upon secondary contact, yet the dynamics and persistence of hybrid incompatibilities in naturally hybridizing species are not well understood. Here, we explore these issues using genetic mapping in three independent populations of recombinant inbred lines between naturally hybridizing monkeyflowers, Mimulus guttatus and Mimulus nasutus, from the sympatric Catherine Creek population. We discover that the three M. guttatus founders differ dramatically in admixture history, with nearly a quarter of one founder's genome introgressed from M. nasutus. Comparative genetic mapping in the three RIL populations reveals three new putative inversions, each one segregating among the M. guttatus founders, two due to admixture. We find strong, genome-wide transmission ratio distortion in all RILs, but patterns are highly variable among the three populations. At least some of this distortion appears to be explained by epistatic selection favouring parental genotypes, but tests of inter-chromosomal linkage disequilibrium also reveal multiple candidate Dobzhansky-Muller incompatibilities. We also map several genetic loci for hybrid pollen viability, including two interacting pairs that coincide with peaks of distortion. Remarkably, even with this limited sample of three M. guttatus lines, we discover abundant segregating variation for hybrid incompatibilities with M. nasutus, suggesting this population harbours diverse contributors to postzygotic isolation. Moreover, even with substantial admixture, hybrid incompatibilities between Mimulus species persist, suggesting postzygotic isolation might be a potent force in maintaining species barriers in this system.

Phylogenomic analyses and chromosome ploidy identification reveal multiple cryptic species in Allium sikkimense complex (Amaryllidaceae).

Polyploidization is a process that typically leads to instantaneous reproductive isolation and has, therefore, been considered as one of the major evolutionary forces in the species-rich Hengduan Mountains (HM), yet this topic remains poorly studied in the region. Allium sikkimense and its relatives (about eight species) compose a natural diploid-polyploid complex with the highest diversity in the HM and adjacent areas. A combination of nuclear ribosomal DNA (nrDNA), plastome, transcriptome, and ploidy identification through chromosome counting and flow cytometry is employed to reconstruct the phylogenetic relationships in this complex and to investigate the frequency and the evolutionary significance of polyploidy in the complex. The plastome failed to resolve the phylogenetic relationships of the different species in the A. sikkimense complex, and the phylogenetic tree based on nrDNA also has limited resolution. However, our study reveals a well-resolved phylogenetic framework for species in the A. sikkimense complex using more than 1,000 orthologous genes from the transcriptome data. Previously recognized morphospecies A. sikkimense are non-monophyletic and comprise at least two independently evolved lineages (i.e., cryptic species), each forming a clade with different diploid species in this complex. The embedded pattern of octoploid A. jichouense and tetraploid A. sp. nov. within different polyploid samples of A. sikkimense supports a possible scenario of budding speciation (via niche divergence). Furthermore, our results reveal that co-occurring species in the A. sikkimense complex usually have different ploidy levels, suggesting that polyploidy is an important process for reproductive isolation of sympatric Allium species. Phylogenetic network analyses suggested that the phylogenetic relationships of the A. sikkimense complex, allowing for reticulation events, always fit the dataset better than a simple bifurcating tree. In addition, the included or exserted filaments, which have long been used to delimit species, are highly unreliable taxonomically due to their extensive parallel and convergent evolution.

Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data.

While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown. To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites. Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study.

Cross compatibility and endogenous phytohormone profiles in interspecific hybridization between Iris tectorum and Iris germanica

In order to introduce new valuable trait to other species, the interspecific hybridization is used in ornamental breeding. The success of interspecific hybridization is hindered by incompatibility occurring at various stages of hybridization. To character the barriers of interspecific hybridization between Iris tectorum and Iris germanica cultivars, this study investigated the compatibility of pre-fertilization and post-fertilization, and profiles of endogenous phytohormones in ovary of self-pollination and interspecific hybridization. After pollination, pollen tubes germinated, grew through style, and entered ovules in both self-pollination and interspecific hybridization. Thus, the interspecific incompatibility between I. tectorum and I. germanica cultivars were not caused by pre-fertilization barrier. The observations of embryo and endosperm development revealed that incompatibility of interspecific-hybridization occurred at post-fertilization stage during the development of the embryo. The degeneration of endosperms was observed at 15 days after pollination and no embryo was observed in the interspecific hybridization. The abnormal development of embryos and endosperms indicates that reproductive barriers between I. tectorum and I. germanica cultivars happened at post-fertilization stage. To character the profiles of endogenous phytohormones in the ovary of self-pollination and interspecific hybridization, concentration of auxins, cytokinins (CKs), gibberellins (GAs), abscisic acid (ABA), Jasmonate (JA), salicylic acid (SA), and their derivatives were analyzed by UPLC-ESI-MS/MS. The results showed concentrations of indole-3-acetic acid (IAA) and bioactive free based of CKs were highest in ovary of I. tectorum self-pollination (IT-Self), which has normally developed embryos and endosperms. Four bioactive forms of GAs were identified in this study, while the patterns of their concentration changes were not constant among different materials. Besides, the concentration of ABA was highest in unfertilized ovary of I. tectorum (IT-CK). According to the profiles of endogenous phytohormones, the abnormal development of embryos in interspecific hybridization might be caused by the low concentration of IAA and CKs, proper concentration of these phytohormone could be added into in vitro culture medium to rescue embryos of interspecific hybridization. The investigation of cross compatibility and endogenous phytohormone profiles in this study will provide insight for overcoming the barrier of interspecific hybridization and embryo rescue in Iris genus.

Exploring the role of polymorphic interspecies structural variants in reproductive isolation and adaptive divergence in Eucalyptus.

Structural variations (SVs) play a significant role in speciation and adaptation in many species, yet few studies have explored the prevalence and impact of different categories of SVs. We conducted a comparative analysis of long-read assembled reference genomes of closely related Eucalyptus species to identify candidate SVs potentially influencing speciation and adaptation. Interspecies SVs can be either fixed differences or polymorphic in one or both species. To describe SV patterns, we employed short-read whole-genome sequencing on over 600 individuals of Eucalyptus melliodora and Eucalyptus sideroxylon, along with recent high-quality genome assemblies. We aligned reads and genotyped interspecies SVs predicted between species reference genomes. Our results revealed that 49,756 of 58,025 and 39,536 of 47,064 interspecies SVs could be typed with short reads in E. melliodora and E. sideroxylon, respectively. Focusing on inversions and translocations, symmetric SVs that are readily genotyped within both populations, 24 were found to be structural divergences, 2,623 structural polymorphisms, and 928 shared structural polymorphisms. We assessed the functional significance of fixed interspecies SVs by examining differences in estimated recombination rates and genetic differentiation between species, revealing a complex history of natural selection. Shared structural polymorphisms displayed enrichment of potentially adaptive genes. Understanding how different classes of genetic mutations contribute to genetic diversity and reproductive barriers is essential for understanding how organisms enhance fitness, adapt to changing environments, and diversify. Our findings reveal the prevalence of interspecies SVs and elucidate their role in genetic differentiation, adaptive evolution, and species divergence within and between populations.

Gene regulation and speciation in a migratory divide between songbirds

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration – an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns – no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.