Mapping seasonal migration in a songbird hybrid zone -- heritability, genetic correlations, and genomic patterns linked to speciation

Seasonal migration is a widespread behavior relevant for adaptation and speciation, yet knowledge of its genetic basis is limited. We leveraged advances in tracking and sequencing technologies to bridge this gap in a well-characterized hybrid zone between songbirds that differ in migratory behavior. Migration requires the coordinated action of many traits, including orientation, timing, and wing morphology. We used genetic mapping to show these traits are highly heritable and genetically correlated, explaining how migration has evolved so rapidly in the past and suggesting future responses to climate change may be possible. Many of these traits mapped to the same genomic regions and small structural variants indicating the same, or tightly linked, genes underlie them. Analyses integrating transcriptomic data indicate cholinergic receptors could control multiple traits. Furthermore, analyses integrating genomic differentiation further suggested genes underlying migratory traits help maintain reproductive isolation in this hybrid zone.

The emerging importance of cross-ploidy hybridisation and introgression

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Integration of natural selection across the life cycle stabilizes a marine mussel hybrid zone

Hybridization among related species is now recognized as common but it remains unclear how hybrid zones persist for prolonged periods. Here, we test the hypothesis that selection in different components of the life cycle may stabilize a hybrid zone. A hybrid zone occurs in southwest England between the marine mussels Mytilus edulis and M. galloprovincialis. Previous studies have found strong directional selection against alleles from M. edulis occurs among hybrids in the adult stage. Traditional hybrid zone models argue that alleles that are selected within the hybrid zone are replaced by migration from neighboring parental population into the hybrid zone. In this system, however, migration occurs out of this hybrid zone into neighboring parental populations. This hybrid zone should therefore be unstable and dissipate, yet this zone has persisted for more than 30 years. We tested and rejected the hypothesis that differences in fecundity may select for M. edulis alleles within this hybrid zone and thus counter the selection observed against these alleles among adults. We also tested the hypothesis that selection during the larval stage may counter selection against M. edulis alleles in the adult stage. We found that selection favors M. edulis alleles during the veliger stage of larval development. The direction and strength of selection during the larval stage are sufficient to counter strong selection during the adult portion of the life cycle. This hybrid zone is stabilized by opposing forms of directional selection operating in different portions of the life cycle.

Tracking hybrid viability across life stages in a natural avian contact zone

Hybrid inviability is an important post-zygotic reproductive barrier between species, but emerging signs of reduced viability can be difficult to study across the lifespan of natural hybrids. We use a combination of long-term monitoring, extra-pair paternity, and mitochondrial DNA identification in a natural hybrid zone of Ficedula flycatchers to detect emerging signs of intrinsic hybrid inviability across their entire lifespan. We evaluate possible evidence of Darwin's corollary to Haldane's rule, predicting asymmetries in inviability between hybrids resulting from reciprocal crosses, due to incompatible genetic factors with sex-specific inheritance patterns. We found higher hatching failure among mixed-species pairs, possibly indicating early developmental impairments associated with specific parental genetic combinations. Adult hybrids had a higher basal mortality rate than both parental species and different age-specific mortality trajectories. There were signs of differences in age-independent mortality rates between the reciprocal hybrid crosses: hybrids with a pied flycatcher mother experienced slightly increased mortality later in life. Using an exceptional dataset with many natural hybrids tracked across life stages, we provide evidence for several emerging signs of reduced hybrid viability. Incompatibilities between alleles located on autosomes and uniparentally inherited factors such as Z-linked and/or mitochondrial genes are strong candidates underlying intrinsic hybrid dysfunction in this system.

Bayesian hybrid index and genomic cline estimation with the R package gghybrid

Admixture, the interbreeding of individuals from differentiated source populations, is now known to be a widespread phenomenon. Genomic studies of natural hybridisation can help to answer many questions on the impacts of admixture on adaptive evolution, reproductive isolation, and speciation. When a large variety of admixture proportions between two source populations exist, both geographic and genomic cline analysis are suitable methods for inferring biased, restricted or excessive gene flow at individual loci into the foreign genomic background, providing evidence for reproductive isolation, selection across an environmental transition, balancing selection, and adaptive introgression. Genomic cline analysis replaces geographic location with genome-wide hybrid index and is therefore useable in circumstances that violate geographic cline assumptions. Here, I introduce gghybrid, an R package for simple and flexible Bayesian estimation of Buerkle's hybrid index and Fitzpatrick's logit-logistic genomic clines using bi-allelic data, suitable for both small and large datasets. gghybrid allows any ploidy and uses Structure input file format. It has separate functions for hybrid index and cline estimation, treating each individual and locus respectively as an independent analysis, making it highly parallelisable. Admixture proportions from other software can alternatively be used in cline analysis, alongside parental allele frequencies. Parameters can be fixed and samples pooled for statistical model comparison with AIC or waic. Here, I describe the functions, pipeline, and statistical properties of gghybrid. Simulations reveal that model comparison with waic is preferred, and use of Bayesian posterior distributions and p values to select candidate non-null loci is problematic and should be avoided.

The effect of host admixture on wild house mouse gut microbiota is weak when accounting for spatial autocorrelation

The question of how interactions between the gut microbiome and vertebrate hosts contribute to host adaptation and speciation is one of the major problems in current evolutionary research. Using bacteriome and mycobiome metabarcoding, we examined how these two components of the gut microbiota vary with the degree of host admixture in secondary contact between two house mouse subspecies (Mus musculus musculus and M. m. domesticus). We used a large data set collected at two replicates of the hybrid zone and model-based statistical analyses to ensure the robustness of our results. Assuming that the microbiota of wild hosts suffers from spatial autocorrelation, we directly compared the results of statistical models that were spatially naive with those that accounted for spatial autocorrelation. We showed that neglecting spatial autocorrelation can strongly affect the results and lead to misleading conclusions. The spatial analyses showed little difference between subspecies, both in microbiome composition and in individual bacterial lineages. Similarly, the degree of admixture had minimal effects on the gut bacteriome and mycobiome and was caused by changes in a few microbial lineages that correspond to the common symbionts of free-living house mice. In contrast to previous studies, these data do not support the hypothesis that the microbiota plays an important role in host reproductive isolation in this particular model system.

Ecological niche contributes to the persistence of the Western x Glaucous-winged Gull hybrid zone

Hybrid zones occur in nature when populations with limited reproductive barriers overlap in space. Many hybrid zones persist over time, and different models have been proposed to explain how selection can maintain hybrid zone stability. More empirical studies are needed to elucidate the role of ecological adaptation in maintaining stable hybrid zones. Here, we investigated the role of exogenous factors in maintaining a hybrid zone between western gulls (Larus occidentalis) and glaucous-winged gulls (L. glaucescens). We used ecological niche models (ENMs) and niche similarity tests to quantify and examine the ecological niches of western gulls, glaucous-winged gulls, and their hybrids. We found evidence of niche divergence between all three groups. Our results best support the bounded superiority model, providing further evidence that exogenous selection favoring hybrids may be an important factor in maintaining this stable hybrid zone.

Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree

BACKGROUND AND AIMS

In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado.

METHODS

The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated.

KEY RESULTS

High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction.

CONCLUSIONS

The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Demography-driven and adaptive introgression in a hybrid zone of the Armeria syngameon

Syngameons represent networks of otherwise distinct species connected by limited gene exchange. Although most studies have focused on how species maintain their cohesiveness despite gene flow, there are additional relevant questions regarding the evolutionary dynamics of syngameons and their drivers, as well as the success of their members and the network as a whole. Using a ddRADseq approach, we analysed the genetic structure, genomic clines and demographic history of a coastal hybrid zone involving two species of the Armeria (Plumbaginaceae) syngameon in southern Spain. We inferred that a peripheral population of the sand dune-adapted A. pungens diverged from the rest of the conspecific populations and subsequently hybridized with a locally more abundant pinewood congener, A. macrophylla. Both species display extensive plastid DNA haplotype sharing. Genomic cline analysis identified bidirectional introgression, but more outlier loci with excess A. pungens than A. macrophylla ancestry, suggesting the possibility of selection for A. pungens alleles. This is consistent with the finding that the A. pungens phenotype is selected for in open habitats, and with the strong correlation found between ancestry and phenotype. Taken together, our analyses suggest an intriguing scenario in which bidirectional introgression may, on the one hand, help to avoid reduced levels of genetic diversity due to the small size and isolated location of the A. pungens range-edge population, thereby minimizing demographic risks of stochastic extinction. On the other hand, the data also suggest that introgression into A. macrophylla may allow individuals to grow in open, highly irradiated, deep sandy, salt-exposed habitats.

Seasonal differences in the timing of flight between the invasive winter moth and native Bruce spanworm promotes reproductive isolation

The European winter moth, Operophtera brumata L. (Lepidoptera: Geometridae), was accidentally introduced to North America on at least 4 separate occasions, where it has been hybridizing with the native Bruce spanworm, O. bruceata Hulst, at rates up to 10% per year. Both species are known to respond to the same sex pheromones and to produce viable offspring, but whether they differ in the seasonal timing of their mating flights is unknown. Therefore, we collected adult male moths weekly along 2 transects in the northeastern United States and genotyped individuals using polymorphic microsatellite markers as males of these 2 species cannot be differentiated morphologically. Along each transect, we then estimated the cumulative proportions (i.e., the number of individuals out of the total collected) of each species on each calendar day. Our results indicate that there are significant differences between the species regarding their seasonal timing of flight, and these allochronic differences likely are acting to promote reproductive isolation between these 2 species. Lastly, our results suggest that the later flight observed by winter moth compared to Bruce spanworm may be limiting its inland spread in the northeastern United States because of increased exposure to extreme winter events.

Speciation patterns of related species under the hybrid zone: A case study of three sclerophyllous oaks in the east Himalaya-Hengduan Mountains

Speciation is often accompanied by frequent gene exchanges, which have been gradually recognized as a common phenomenon in nature. Although gene flow may influence different reproductive isolations, the specific mechanism of the process still requires more experimental evidence, especially in hybrid populations that have no significant differentiation and isolation. To address this challenge, this study aims to elucidate the underlying mechanisms of sympatry and parapatry in related species. Three sclerophyllous oaks (Quercus spinosa, Quercus aquifolioides and Quercus rehderiana), which are mainly distributed in the sympatry/parapatry of the East Himalaya-Hengduan Mountains and adjacent regions, were used to explore the population dynamics and evolutionary history. Based on 12,420 genome-wide single nucleotide polymorphism datasets, gene flow detection showed that the three species did not present any obvious genetic barriers. Evolutionary analysis indicated that the three species diverged during the Tertiary Period, and no migration events occurred in the early stages of species divergence. Combined with the data of 19 ecological factors, geological movements and climatic turbulence caused the rapid radiated differentiation of the three species during the Neocene, and similar selective pressures resulted in the same evolutionary pattern based on demographic history analysis. In addition, the predicted niche occupancy profiles and Generalized Dissimilarity Modelling revealed that the three species occupied distinct niches and had significant differences in ecological adaptation, which may explain the specific morphological characteristics of the different species. Therefore, we believe that the populations of the three related species underwent adaptive evolution in different habitats during the early stages of divergence. This study provides new experimental evidence of the formation patterns of parallel speciation.

Hybridization despite elaborate courtship behavior and female choice in Neotropical tree frogs

The mechanisms of hybridization can be elucidated by analyzing genotypes as well as phenotypes that could act as premating barriers, as the reproductive interactions among heterospecifics can alter the evolutionary history of species. In frogs, hybrids typically occur among species that reproduce explosively (in dense aggregations) with few opportunities for mate selection but are rare in species with elaborate courtship behaviors that may prevent erroneous mating. Using 21 microsatellite markers, we examined hybridization in the prolonged-breeding tree frogs Bokermannohyla ibitiguara and B. sazimai sampled within a contact zone in the Brazilian savanna (72 tadpoles; 74 adults). We also compared acoustic and morphological data. We confirmed both parental species genetically; STRUCTURE results confirmed 14 hybrids, 11 of which were second-generation according to NEWHYBRIDS, all with intermediate values of genetic dissimilarities compared to the parentals. Morphological and acoustic analyses revealed that hybrids showed variable but not necessarily intermediate phenotypes. Moreover, 2 hybrids exhibited call types different from parentals. The reproduction of B. ibitiguara involves territorial and aggressive males, elaborate courtships with acoustic and tactile stimuli, choosy females, and opportunistic strategies. Our study uncovers a rare case of viable hybridization among closely related frogs with such a combination of complex courtship behaviors and mate choice. We discuss the likely directionality and mechanisms behind this phenomenon, and highlight the importance of investigating hybridization even in species that show elaborate reproduction and female choice to advance our understanding of animal diversification.

Experimentally testing mate preference in an avian system with unidirectional bill color introgression

Mating behavior can play a key role in speciation by inhibiting or facilitating gene flow between closely related taxa. Hybrid zones facilitate a direct examination of mating behavior and the traits involved in establishing species barriers. The long-tailed finch (Poephila acuticauda) has two hybridizing subspecies that differ in bill color (red and yellow), and the yellow bill phenotype appears to have introgressed ~350 km eastward following secondary contact. To examine the role of mate choice on bill color introgression, we performed behavioral assays using natural and manipulated bill colors. We found an assortative female mating preference for males of their own subspecies when bill color was not manipulated. However, we did not find this assortative preference in trials based on artificially manipulated bill color. This could suggest that assortative preference is not fixed entirely on bill color and instead may be based on a different trait (e.g., song) or a combination of traits, or alternatively may be due to lower statistical power alongside the bill manipulations being unconvincing to the female choosers. Intriguingly, we find a bias in the inheritance of bill color in captive bred F1 hybrid females. Previous modeling suggests that assortative mate preference and this kind of partial dominance in the underlying genes may together contribute to introgression, making the genetic architecture of bill color in this system a priority for future research.

Tandem Repeat DNA Provides Many Cytological Markers for Hybrid Zone Analysis in Two Subspecies of the Grasshopper Chorthippus parallelus

Recent advances in next generation sequencing (NGS) have greatly increased our understanding of non-coding tandem repeat (TR) DNA. Here we show how TR DNA can be useful for the study of hybrid zones (HZ), as it serves as a marker to identify introgression in areas where two biological entities come in contact. We used Illumina libraries to analyse two subspecies of the grasshopper Chorthippus parallelus, which currently form a HZ in the Pyrenees. We retrieved a total of 152 TR sequences, and used fluorescent in situ hybridization (FISH) to map 77 families in purebred individuals from both subspecies. Our analysis revealed 50 TR families that could serve as markers for analysis of this HZ, using FISH. Differential TR bands were unevenly distributed between chromosomes and subspecies. Some of these TR families yielded FISH bands in only one of the subspecies, suggesting the amplification of these TR families after the geographic separation of the subspecies in the Pleistocene. Our cytological analysis of two TR markers along a transect of the Pyrenean hybrid zone showed asymmetrical introgression of one subspecies into the other, consistent with previous findings using other markers. These results demonstrate the reliability of TR-band markers for hybrid zone studies.

The gut microbiome reflects ancestry despite dietary shifts across a hybrid zone

The microbiome is critical to an organism's phenotype, and its composition is shaped by, and a driver of, eco-evolutionary interactions. We investigated how host ancestry, habitat and diet shape gut microbial composition in a mammalian hybrid zone between Neotoma lepida and N. bryanti that occurs across an ecotone between distinct vegetation communities. We found that habitat is the primary determinant of diet, while host genotype is the primary determinant of the gut microbiome-a finding further supported by intermediate microbiome composition in first-generation hybrids. Despite these distinct primary drivers, microbial richness was correlated with diet richness, and individuals that maintained higher dietary richness had greater gut microbial community stability. Both relationships were stronger in the relative dietary generalist of the two parental species. Our findings show that host ancestry interacts with dietary habits to shape the microbiome, ultimately resulting in the phenotypic plasticity that host-microbial interactions allow.

Species recognition limits mating between hybridizing ant species

Identifying mechanisms limiting hybridization is a central goal of speciation research. Here, we studied premating and postmating barriers to hybridization between two ant species, Formica selysi and Formica cinerea. These species hybridize in the Rhône valley in Switzerland, where they form a mosaic hybrid zone, with limited introgression from F. selysi into F. cinerea. There was no sign of temporal isolation between the two species in the production of queens and males. With choice experiments, we showed that queens and males strongly prefer to mate with conspecifics. Yet, we did not detect postmating barriers caused by genetic incompatibilities. Specifically, hybrids of all sexes and castes were found in the field and F1 hybrid workers did not show reduced viability compared to nonhybrid workers. To gain insights into the cues involved in species recognition, we analyzed the cuticular hydrocarbons (CHCs) of queens, males, and workers and staged dyadic encounters between workers. CHC profiles differed markedly between species, but were similar in F. cinerea and hybrids. Accordingly, workers also discriminated species, but they did not discriminate F. cinerea and hybrids. We discuss how the CHC-based recognition system of ants may facilitate the establishment of premating barriers to hybridization, independent of hybridization costs.

Restricted X chromosome introgression and support for Haldane's rule in hybridizing damselflies

Contemporary hybrid zones act as natural laboratories for the investigation of species boundaries and may shed light on the little understood roles of sex chromosomes in species divergence. Sex chromosomes are considered to function as a hotspot of genetic divergence between species; indicated by less genomic introgression compared to autosomes during hybridization. Moreover, they are thought to contribute to Haldane's rule, which states that hybrids of the heterogametic sex are more likely to be inviable or sterile. To test these hypotheses, we used contemporary hybrid zones of Ischnura elegans, a damselfly species that has been expanding its range into the northern and western regions of Spain, leading to chronic hybridization with its sister species Ischnura graellsii. We analysed genome-wide SNPs in the Spanish I. elegans and I. graellsii hybrid zone and found (i) that the X chromosome shows less genomic introgression compared to autosomes, and (ii) that males are underrepresented among admixed individuals, as predicted by Haldane's rule. This is the first study in Odonata that suggests a role of the X chromosome in reproductive isolation. Moreover, our data add to the few studies on species with X0 sex determination system and contradict the hypothesis that the absence of a Y chromosome causes exceptions to Haldane's rule.

Extensive Sympatry and Frequent Hybridization of Ecologically Divergent Aquatic Plants on the Qinghai-Tibetan Plateau

Hybridization has fascinated biologists in recent centuries for its evolutionary importance, especially in plants. Hybrid zones are commonly located in regions across environmental gradients due to more opportunities to contact and ecological heterogeneity. For aquatic taxa, intrazonal character makes broad overlapping regions in intermediate environments between related species. However, we have limited information on the hybridization pattern of aquatic taxa in alpines, especially submerged macrophytes. In this study, we aimed to test the hypotheses that niche overlap and hybridization might be extensive in related aquatic plants across an altitudinal gradient. We evaluated the niche overlap in three related species pairs on the Qinghai-Tibetan Plateau and assessed the spatial pattern of hybrid populations. Obvious niche overlap and common hybridization were revealed in all three pairs of related aquatic plants. The plateau edge and river basins were broad areas for the sympatry of divergent taxa, where a large proportion of hybrid populations occurred. Hybrids are also discretely distributed in diverse habitats on the plateau. Differences in the extent of niche overlap, genetic incompatibility and phylogeographic history might lead to variation differences in hybridization patterns among the three species pairs. Our results suggested that plateau areas are a hotspot for ecologically divergent aquatic species to contact and mate and implied that hybridization may be important for the freshwater biodiversity of highlands.

Climate Warming Since the Holocene Accelerates West-East Communication for the Eurasian Temperate Water Strider Species Aquarius paludum

Holocene climate warming has dramatically altered biological diversity and distributions. Recent human-induced emissions of greenhouse gases will exacerbate global warming and thus induce threats to cold-adapted taxa. However, the impacts of this major climate change on transcontinental temperate species are still poorly understood. Here, we generated extensive genomic datasets for a water strider, Aquarius paludum, which was sampled across its entire distribution in Eurasia and used these datasets in combination with ecological niche modeling (ENM) to elucidate the influence of the Holocene and future climate warming on its population structure and demographic history. We found that A. paludum consisted of two phylogeographic lineages that diverged in the middle Pleistocene, which resulted in a “west–east component” genetic pattern that was probably triggered by Central Asia-Mongoxin aridification and Pleistocene glaciations. The diverged western and eastern lineages had a second contact in the Holocene, which shaped a temporary hybrid zone located at the boundary of the arid–semiarid regions of China. Future predictions detected a potentially novel northern corridor to connect the western and eastern populations, indicating west–east gene flow would possibly continue to intensify under future warming climate conditions. Further integrating phylogeographic and ENM analyses of multiple Eurasian temperate taxa based on published studies reinforced our findings on the “west–east component” genetic pattern and the predicted future northern corridor for A. paludum. Our study provided a detailed paradigm from a phylogeographic perspective of how transcontinental temperate species differ from cold-adapted taxa in their response to climate warming.

A migratory divide spanning two continents is associated with genomic and ecological divergence

Migratory divides are contact zones between breeding populations with divergent migratory strategies during the nonbreeding season. These locations provide an opportunity to evaluate the role of seasonal migration in the maintenance of reproductive isolation, particularly the relationship between population structure and features associated with distinct migratory strategies. We combine light-level geolocators, genomic sequencing, and stable isotopes to investigate the timing of migration and migratory routes of individuals breeding on either side of a migratory divide coinciding with genomic differentiation across a hybrid zone between barn swallow (Hirundo rustica) subspecies in China. Individuals west of the hybrid zone, with H. r. rustica ancestry, had comparatively enriched stable-carbon and hydrogen isotope values and overwintered in eastern Africa, whereas birds east of the hybrid zone, with H. r. gutturalis ancestry, had depleted isotope values and migrated to southern India. The two subspecies took divergent migratory routes around the high-altitude Karakoram Range and arrived on the breeding grounds over 3 weeks apart. These results indicate that assortative mating by timing of arrival and/or selection against hybrids with intermediate migratory traits may maintain reproductive isolation between the subspecies, and that inhospitable geographic features may have contributed to the diversification of Asian avifauna by influencing migratory patterns.

Discrimination of conspecifics from heterospecifics in a hybrid zone: Behavioral and chemical cues in ants

Species and nestmate recognition in social insects occurs mostly through cuticular hydrocarbons acting as chemical cues. These compounds generate a colony-specific odor profile depending on genetic and environmental factors. Species and nestmate recognition results in specific behavioral responses, regulating the level of aggression toward other individuals during an interaction. Although species discrimination and recognition cues have been poorly studied in the context of interspecific hybridization, such systems offer an opportunity to further investigate the influence of heritable and environmental factors on recognition. We explored the strength of discrimination in a hybrid zone between two ant species-Tetramorium immigrans and T. caespitum-by comparing cuticular hydrocarbon profiles and measuring intra- and interspecific worker aggression in both areas of sympatry and areas of allopatry among species. Species cuticular hydrocarbon profiles were well-differentiated and interspecific aggression was high, revealing highly discriminating species recognition cues. Hybrids' cuticular hydrocarbon profiles consisted of a mixture of the parental bouquets, but also exhibited hybrid-specific patterns. Behavioral assays showed that T. immigrans is as aggressive toward hybrids as toward heterospecifics. Finally, aggression between heterospecific workers was lower when interacting individuals came from areas of sympatry among species than from areas of allopatry. Taken as a whole, these findings paint a particularly complex picture of the recognition system in T. immigrans, T. caespitum, and their hybrids, and highlight that hybrid zones afford a still underexplored opportunity for investigating recognition mechanisms and discrimination between species.

Expansion of the rare Eucalyptus risdonii under climate change through hybridization with a closely related species despite hybrid inferiority

BACKGROUND AND AIMS

Hybridization is increasingly recognized as an integral part of the dynamics of species range expansion and contraction. Thus, it is important to understand the reproductive barriers between co-occurring species. Extending previous studies that argued that the rare Eucalyptus risdonii was expanding into the range of the surrounding E. amygdalina by both seed and pollen dispersal, we here investigate the long-term fitness of both species and their hybrids and whether expansion is continuing.

METHODS

We assessed the survival of phenotypes representing a continuum between the two pure species in a natural hybrid swarm after 29 years, along with seedling recruitment. The performance of pure species as well as of artificial and natural hybrids was also assessed over 28 years in a common garden trial.

KEY RESULTS

In the hybrid zone, E. amygdalina adults showed greater mortality than E. risdonii, and the current seedling cohort is still dominated by E. risdonii phenotypes. Morphologically intermediate individuals appeared to be the least fit. Similar results were observed after growing artificial first-generation and natural hybrids alongside pure species families in a common garden trial. Here, the survival, reproduction, health and growth of the intermediate hybrids were significantly less than those of either pure species, consistent with hybrid inferiority, although this did not manifest until later reproductive ages. Among the variable progeny of natural intermediate hybrids, the most E. risdonii-like phenotypes were the most fit.

CONCLUSIONS

This study contributes to the increasing number of reports of hybrid inferiority in Eucalyptus, suggesting that post-zygotic barriers contribute to the maintenance of species integrity even between closely related species. However, with fitness rapidly recovered following backcrossing, it is argued that hybridization can still be an important evolutionary process, in the present case appearing to contribute to the range expansion of the rare E. risdonii in response to climate change.

Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree Plathymenia reticulata (n = 233 individuals). Genetic structure analysis of P. reticulata was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F CT = 0.222 and F ST = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F2) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in P. reticulata. The low genetic diversity of the P. reticulata forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata)

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7–10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Mito-nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern "Rocky Mountain'' color form with ferruginous mid-abdominal segments (B. m. melanopygus) and a southern "Pacific'' form with black mid-abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color-controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito-nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.

A de novo genome assembly and annotation of the southern flying squirrel (Glaucomys volans)

Northern (Glaucomys sabrinus) and southern (Glaucomys volans) flying squirrels are widespread species distributed across North America. Northern flying squirrels are common inhabitants of the boreal forest, also occurring in coniferous forest remnants farther south, whereas the southern flying squirrel range is centered in eastern temperate woodlands. These two flying squirrel species exhibit a hybrid zone across a latitudinal gradient in an area of recent secondary contact. Glaucomys hybrid offspring are viable and can successfully backcross with either parental species, however, the fitness implications of such events are currently unknown. Some populations of G. sabrinus are endangered, and thus, interspecific hybridization is a key conservation concern in flying squirrels. To provide a resource for future studies to evaluate hybridization and possible introgression, we sequenced and assembled a de novo long-read genome from a G. volans individual sampled in southern Ontario, Canada, while four short-read genomes (two G. sabrinus and two G. volans, all from Ontario) were resequenced on Illumina platforms. The final genome assembly consisted of approximately 2.40 Gb with a scaffold N50 of 455.26 Kb. Benchmarking Universal Single-Copy Orthologs reconstructed 3,742 (91.2%) complete mammalian genes and genome annotation using RNA-Seq identified the locations of 19,124 protein-coding genes. The four short-read individuals were aligned to our reference genome to investigate the demographic history of the two species. A principal component analysis clearly separated resequenced individuals, while inferring population size history using the Pairwise Sequentially Markovian Coalescent model noted an approximate species split 1 million years ago, and a single, possibly recently introgressed individual.

Genetic structure in Louisiana Iris species reveals patterns of recent and historical admixture

PREMISE

When divergent lineages come into secondary contact, reproductive isolation may be incomplete, thus providing an opportunity to investigate how speciation is manifested in the genome. The Louisiana Irises (Iris, series Hexagonae) comprise a group of three or more ecologically and reproductively divergent lineages that can produce hybrids where they come into contact. We estimated standing genetic variation to understand the current distribution of population structure in the Louisiana Irises.

METHODS

We used genotyping-by-sequencing techniques to sample the genomes of Louisiana Iris species across their ranges. We sampled 20 populations (n = 632 individuals) across 11,249 loci and used Entropy and PCA models to assess population genetic data.

RESULTS

We discovered evidence for interspecific gene flow in parts of the range. Our analysis revealed patterns of population structure at odds with widely accepted nominal taxonomy. We discovered undescribed hybrid populations, designated as belonging to the I. brevicaulis lineage. Iris nelsonii shared significant ancestry with only one of the purported parent species, I. fulva, evidence inconsistent with a hybrid origin.

CONCLUSIONS

This study provides several key findings important to the investigation of standing genetic variation in the Louisiana Iris species complex. Compared to the other nominal species, I. brevicaulis contains a large amount of genetic diversity. In addition, we discovered a previously unknown hybrid zone between I. brevicaulis and I. hexagona along the Texas coast. Finally, our results do not support the long-standing hypothesis that I. nelsonii has mixed ancestry from three parental taxa.

Assortative mating and epistatic mating-trait architecture induce complex movement of the crow hybrid zone

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

Host genetic factors associated with the range limit of a European hantavirus

The natural host ranges of many viruses are restricted to very specific taxa. Little is known about the molecular barriers between species that lead to the establishment of this restriction or generally prevent virus emergence in new hosts. Here, we identify genomic polymorphisms in a natural rodent host associated with a strong genetic barrier to the transmission of European Tula orthohantavirus (TULV). We analysed the very abrupt spatial transition between two major phylogenetic clades in TULV across the comparatively much wider natural hybrid zone between evolutionary lineages of their reservoir host, the common vole (Microtus arvalis). Genomic scans of 79,225 single nucleotide polymorphisms (SNPs) in 323 TULV‐infected host individuals detected 30 SNPs that were consistently associated with the TULV clades CEN.S or EST.S in two replicate sampling transects. Focusing the analysis on 199 voles with evidence of genomic admixture at the individual level (0.1–0.9) supported statistical significance for all 30 loci. Host genomic variation at these SNPs explained up to 37.6% of clade‐specific TULV infections. Genes in the vicinity of associated SNPs include SAHH, ITCH and two members of the Syngr gene family, which are involved in functions related to immune response or membrane transport. This study demonstrates the relevance of natural hybrid zones as systems not only for studying processes of evolutionary divergence and speciation, but also for the detection of evolving genetic barriers for specialized parasites.

Genetic ancestry predicts male-female affiliation in a natural baboon hybrid zone

Opposite-sex social relationships are important predictors of fitness in many animals, including several group-living mammals. Consequently, understanding sources of variance in the tendency to form opposite-sex relationships is important for understanding social evolution. Genetic contributions are of particular interest due to their importance in long-term evolutionary change, but little is known about genetic effects on male-female relationships in social mammals, especially outside of the mating context. Here, we investigate the effects of genetic ancestry on male-female affiliative behaviour in a hybrid zone between the yellow baboon, Papio cynocephalus, and the anubis baboon, Papio anubis, in a population in which male-female social bonds are known predictors of life span. We place our analysis within the context of other social and demographic predictors of affiliative behaviour in baboons. Genetic ancestry was the most consistent predictor of opposite-sex affiliative behaviour we observed, with the exception of strong effects of dominance rank. Our results show that increased anubis genetic ancestry is associated with a subtle, but significantly higher, probability of opposite-sex affiliative behaviour, in both males and females. Additionally, pairs of anubis-like males and anubis-like females were the most likely to socially affiliate, resulting in moderate assortativity in grooming and proximity behaviour as a function of genetic ancestry. Our findings indicate that opposite-sex affiliative behaviour partially diverged during baboon evolution to differentiate yellow and anubis baboons, despite overall similarities in their social structures and mating systems. Furthermore, they suggest that affiliative behaviour may simultaneously promote and constrain baboon admixture, through additive and assortative effects of ancestry, respectively.

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones

PREMISE

The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data.

METHODS

We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape.

RESULTS

Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations.

CONCLUSIONS

The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

Whole-genome sequencing reveals asymmetric introgression between two sister species of cold-resistant leaf beetles

A large number of genetic variation studies have identified cases of mitochondrial genome introgression in animals, indicating that reproductive barriers among closely related species are often permeable. Because of its sheer size, the impact of hybridization on the evolution of the nuclear genome is more difficult to apprehend. Only a few studies have explored it recently thanks to recent progress in DNA sequencing and genome assembly. Here, we analysed whole-genome sequence variation among multiple individuals of two sister species of leaf beetles inside their hybrid zone, in which asymmetric mitochondrial genome introgression had previously been established. We used a machine learning approach based on computer simulations for training to identify regions of the nuclear genome that were introgressed. We inferred asymmetric introgression of ≈2% of the genome, in the same direction that was observed for the mitochondrial genome. Because a previous study based on a reduced-representation sequencing approach was not able to detect this introgression, we conclude that whole-genome sequencing is necessary when the fraction of the introgressed genome is small. We also analysed the whole-genome sequence of a hybrid individual, demonstrating that hybrids have the capacity to backcross with the species for which virtually no introgression was observed. Our data suggest that one species has recently invaded the range of the other and/or some alleles that where transferred from the invaded into the invading species could be under positive selection and may have favoured the adaptation of the invading species to the Alpine environment.

Haplotype tagging reveals parallel formation of hybrid races in two butterfly species

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.

Lack of assortative mating might explain reduced phenotypic differentiation where two grasshopper species meet

Hybridization is an evolutionary process with wide‐ranging potential outcomes, from providing populations with important genetic variation for adaptation to being a substantial fitness cost leading to extinction. Here, we focussed on putative hybridization between two morphologically distinct species of New Zealand grasshopper. We collected Phaulacridium marginale and Phaulacridium otagoense specimens from a region where mitochondrial introgression had been detected and where their habitat has been modified by introduced mammals eating the natural vegetation and by the colonization of many non‐native plant species. In contrast to observations in the 1970s, our sampling of wild pairs of grasshoppers in copula provided no evidence of assortative mating with respect to species. Geometric morphometrics on pronotum shape of individuals from areas of sympatry detected phenotypically intermediate specimens (putative hybrids), and the distribution of phenotypes in most areas of sympatry was found to be unimodal. These results suggest that hybridization associated with anthropogenic habitat changes has led to these closely related species forming a hybrid swarm, with random mating. Without evidence of hybrid disadvantage, we suggest a novel hybrid lineage might eventually result from the merging of these two species.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Primary hybrid zone formation in Tephroseris helenitis (Asteraceae), following postglacial range expansion along the central Northern Alps

Distinguishing between secondary versus primary hybrid zone formation remains a challenging task as, for instance, the time window in which these historical (vicariant) versus contemporary (environmental-selective) processes are distinguishable may be relatively narrow. Here, we examine the origin and structure of a transition zone between two subspecies of Tephroseris helenitis along the central Northern Alps, using molecular (AFLP) and morphological (achene type) data in combination with ecological niche models (ENMs) to hindcast ranges at the Last Glacial Maximum (LGM) and mid-Holocene. Samples were collected over a c. 350 km long transect, largely covered by ice during the LGM. Genetically nonadmixed individuals of subspp. helenitis versus salisburgensis dominated the westernmost versus eastern transect areas, with admixed individuals occurring in between. Clines for achene morphology and outlier loci potentially under climate-driven selection were steep, largely noncoincidental, and displaced to the east of the cline centre for neutral AFLPs. During the LGM, ssp. helenitis should have been able to persist in a refugium southwest of the transect, while suitable habitat for ssp. salisburgensis was apparently absent at this time. Together with patterns of genetic and clinal variation, our ENM data are suggestive of a primary hybrid zone that originated after the species' postglacial, eastward expansion. The observed clinal changes may thus reflect random/nonadaptive processes during expansion and selection on particular loci, and possibly achene type, in response to a long-term, west-to-east climate gradient in the direction of more stressful (e.g., wetter/cooler) conditions. Overall, this study adds to the vast hybrid zone literature a rare example of a hybrid zone caused by primary differentiation within a plant species, underlaid by historical range expansion.

Characterization of a pericentric inversion in plateau fence lizards (Sceloporus tristichus): evidence from chromosome-scale genomes

Spiny lizards in the genus Sceloporus are a model system among squamate reptiles for studies of chromosomal evolution. While most pleurodont iguanians retain an ancestral karyotype formula of 2n = 36 chromosomes, Sceloporus exhibits substantial karyotype variation ranging from 2n =  22 to 46 chromosomes. We present two annotated chromosome-scale genome assemblies for the Plateau Fence Lizard (Sceloporus tristichus) to facilitate research on the role of pericentric inversion polymorphisms on adaptation and speciation. Based on previous karyotype work using conventional staining, the S. tristichus genome is characterized as 2n =  22 with six pairs of macrochromosomes and five pairs of microchromosomes and a pericentric inversion polymorphism on chromosome 7 that is geographically variable. We provide annotated, chromosome-scale genomes for two lizards located at opposite ends of a dynamic hybrid zone that are each fixed for different inversion polymorphisms. The assembled genomes are 1.84-1.87 Gb (1.72 Gb for scaffolds mapping to chromosomes) with a scaffold N50 of 267.5 Mb. Functional annotation of the genomes resulted in ∼15K predicted gene models. Our assemblies confirmed the presence of a 4.62-Mb pericentric inversion on chromosome 7, which contains 62 annotated coding genes with known functions. In addition, we collected population genomics data using double digest RAD-sequencing for 44 S. tristichus to estimate population structure and phylogeny across the Colorado Plateau. These new genomic resources provide opportunities to perform genomic scans and investigate the formation and spread of pericentric inversions in a naturally occurring hybrid zone.

Extensive hybridization reveals multiple coloration genes underlying a complex plumage phenotype

Coloration is an important target of both natural and sexual selection. Discovering the genetic basis of colour differences can help us to understand how this visually striking phenotype evolves. Hybridizing taxa with both clear colour differences and shallow genomic divergences are unusually tractable for associating coloration phenotypes with their causal genotypes. Here, we leverage the extensive admixture between two common North American woodpeckers-yellow-shafted and red-shafted flickers-to identify the genomic bases of six distinct plumage patches involving both melanin and carotenoid pigments. Comparisons between flickers across approximately 7.25 million genome-wide SNPs show that these two forms differ at only a small proportion of the genome (mean F_ST = 0.008). Within the few highly differentiated genomic regions, we identify 368 SNPs significantly associated with four of the six plumage patches. These SNPs are linked to multiple genes known to be involved in melanin and carotenoid pigmentation. For example, a gene (CYP2J19) known to cause yellow to red colour transitions in other birds is strongly associated with the yellow versus red differences in the wing and tail feathers of these flickers. Additionally, our analyses suggest novel links between known melanin genes and carotenoid coloration. Our finding of patch-specific control of plumage coloration adds to the growing body of literature suggesting colour diversity in animals could be created through selection acting on novel combinations of coloration genes.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

Landscape and taxon age are associated with differing patterns of hybridization in two Eucalyptus (Myrtaceae) subgenera

BACKGROUND AND AIMS

Hybridization is an important evolutionary process that can have a significant impact on natural plant populations. Eucalyptus species are well known for weak reproductive barriers and extensive hybridization within subgenera, but there is little knowledge of whether patterns of hybridization differ among subgenera. Here, we examine eucalypts of Western Australia's Stirling Range to investigate how patterns of hybridization are associated with landscape and taxon age between the two largest Eucalyptus subgenera: Eucalyptus and Symphyomyrtus. In doing so, we tested a hypothesis of OCBIL (old, climatically buffered, infertile landscape) theory that predicts reduced hybridization on older landscapes.

METHODS

Single nucleotide polymorphism markers were applied to confirm the hybrid status, parentage and genetic structure of five suspected hybrid combinations for subg. Eucalyptus and three combinations for subg. Symphyomyrtus.

KEY RESULTS

Evidence of hybridization was found in all combinations, and parental taxa were identified for most combinations. The older parental taxa assessed within subg. Eucalyptus, which are widespread on old landscapes, were identified as well-defined genetic entities and all hybrids were exclusively F1 hybrids. In addition, many combinations showed evidence of clonality, suggesting that the large number of hybrids recorded from some combinations is the result of long-term clonal spread following a few hybridization events rather than frequent hybridization. In contrast, the species in subg. Symphyomyrtus, which typically occur on younger landscapes and are more recently evolved, showed less distinction among parental taxa, and where hybridization was detected, there were high levels of introgression.

CONCLUSIONS

Reduced hybridization in subg. Eucalyptus relative to extensive hybridization in subg. Symphyomyrtus affirmed the hypothesis of reduced hybridization on OCBILs and demonstrate that clade divergence times, landscape age and clonality are important drivers of differing patterns of speciation and hybridization in Eucalyptus.

Asymmetry in fitness-related traits of later-generation hybrids between two invasive species

PREMISE

The importance of hybridization to invasion has been frequently discussed, with most studies focusing on the comparison of fitness-related traits between F1 hybrids and their parents and the consequences of such fitness differences. However, relatively little attention has been given to late-generation hybrids. Different fitness landscapes could emerge in later generations after hybrids cross with each other or backcross with their parents, which may play an important role in plant invasion and subsequent speciation.

METHODS

In this study, artificial crosses were conducted to generate multiple generations, including F1, F2, and backcrosses between two invasive species: Cakile edentula (self-compatible) and C. maritima (self-incompatible). Putative hybrids were also collected in the sympatric zone and compared with their co-occurring parents for phenotypic and genetic differences.

RESULTS

Genetic data provided evidence of hybridization happening in the wild, and phenotypic comparisons showed that natural hybrids had intermediate traits between the two species but showed more similarity to C. maritima than to C. edentula. The asymmetry was further identified in artificial generations for several phenotypic characters. Furthermore, backcrosses exhibited different patterns of variation, with backcrosses to C. maritima having higher reproductive output than their counterparts.

CONCLUSIONS

Our results suggest that hybridization and introgression (backcrossing) in Cakile species is asymmetric and most likely to favor the proliferation of C. maritima genes in the mixed population and thus help its establishment, a finding that could not be predicted by characterizing F1 hybrids.

Natural selection maintains species despite frequent hybridization in the desert shrub Encelia

Natural selection is an important driver of genetic and phenotypic differentiation between species. For species in which potential gene flow is high but realized gene flow is low, adaptation via natural selection may be a particularly important force maintaining species. For a recent radiation of New World desert shrubs (Encelia: Asteraceae), we use fine-scale geographic sampling and population genomics to determine patterns of gene flow across two hybrid zones formed between two independent pairs of species with parapatric distributions. After finding evidence for extremely strong selection at both hybrid zones, we use a combination of field experiments, high-resolution imaging, and physiological measurements to determine the ecological basis for selection at one of the hybrid zones. Our results identify multiple ecological mechanisms of selection (drought, salinity, herbivory, and burial) that together are sufficient to maintain species boundaries despite high rates of hybridization. Given that multiple pairs of Encelia species hybridize at ecologically divergent parapatric boundaries, such mechanisms may maintain species boundaries throughout Encelia.

Pre-introduction introgression contributes to parallel differentiation and contrasting hybridization outcomes between invasive and native marine mussels

Non-native species experience novel selection pressures in introduced environments and may interbreed with native lineages. Species introductions therefore provide opportunities to investigate repeated patterns of adaptation and introgression across replicated contact zones. Here, we investigate genetic parallelism between multiple introduced populations of the invasive marine mussel, Mytilus galloprovincialis, in the absence (South Africa and California) and presence of hybridization with a native congener (Mytilus planulatus in Batemans Bay and Sydney Harbour, Australia). Repeatability in post-introduction differentiation from native-range populations varied between genetically distinct Atlantic and Mediterranean lineages, with Atlantic-derived introductions displaying high differentiation (maxF_ST  > 0.4) and parallelism at outlier loci. Identification of long noncoding RNA transcripts (lncRNA) additionally allowed us to clarify that parallel responses are largely limited to protein-coding loci, with lncRNAs likely evolving under evolutionary constraints. Comparisons of independent hybrid zones revealed differential introgression most strongly in Batemans Bay, with an excess of M. galloprovincialis ancestry and resistance to introgression at loci differentiating parental lineages (M. planulatus and Atlantic M. galloprovincialis). Additionally, contigs putatively introgressed with divergent alleles from a closely related species, Mytilus edulis, showed stronger introgression asymmetries compared with genome-wide trends and also diverged in parallel in both Atlantic-derived introductions. These results suggest that divergent demographic histories experienced by introduced lineages, including pre-introduction introgression, influence contemporary admixture dynamics. Our findings build on previous investigations reporting contributions of historical introgression to intrinsic reproductive architectures shared between marine lineages and illustrate that interspecific introgression history can shape differentiation between colonizing populations and their hybridization with native congeners.

Population genetic structure underlying the geographic variation in beetle structural colour with multiple transition zones

We studied the population genetic structure underlying the geographic variation in the structural colour of the geotrupid dung beetle, Phelotrupes auratus, which exhibits metallic body colours of different reflectance wavelengths perceived as red, green and indigo. These forms occur parapatrically in an area of Japan. The colour variation was not related to variation in climatic factors. Using single nucleotide polymorphisms (SNPs) from restriction-site-associated DNA sequences, we discriminated five groups of populations (west/red, south/green, south/indigo, south/red and east/red) by a combination of genetic clusters (west, south and east) and three colour forms. There were three transition zones for the colour forms: two between the red and green forms were hybrid zones with steep genetic clines, which implies the existence of barriers to gene flow between regions with different colours. The remaining transition zone between the green and indigo forms lacked genetic differentiation, despite the evident colour changes, which implies regionally specific selection on the different colours. In a genomewide association study, we identified four SNPs that were associated with the red/green or indigo colour and were not linked with one another, which implies that the coloration was controlled by multiple loci, each affecting the expression of a different colour range. These loci may have controlled the transitions between different combinations of colours. Our study demonstrates that geographic colour variation within a species can be maintained by nonuniform interactions among barriers to gene flow, locally specific selection on different colours, and the effects of different colour loci.

Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone

Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the “normal” integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which “purge” extreme variants in hybrids (e.g. reduced fitness of hybrid shape).

Selection on a small genomic region underpins differentiation in multiple color traits between two warbler species

Speciation is one of the most important processes in biology, yet the study of the genomic changes underlying this process is in its infancy. North American warbler species Setophaga townsendi and Setophaga occidentalis hybridize in a stable hybrid zone, following a period of geographic separation. Genomic differentiation accumulated during geographic isolation can be homogenized by introgression at secondary contact, whereas genetic regions that cause low hybrid fitness can be shielded from such introgression. Here, we examined the genomic underpinning of speciation by investigating (1) the genetic basis of divergent pigmentation traits between species, (2) variation in differentiation across the genome, and (3) the evidence for selection maintaining differentiation in the pigmentation genes. Using tens of thousands of single nucleotide polymorphisms (SNPs) genotyped in hundreds of individuals within and near the hybrid zone, genome-wide association mapping revealed a single SNP associated with cheek, crown, breast coloration, and flank streaking, reflecting pleiotropy (one gene affecting multiple traits) or close physical linkage of different genes affecting different traits. This SNP is within an intron of the RALY gene, hence we refer to it as the RALY SNP. We then examined between-species genomic differentiation, using both genotyping-by-sequencing and whole genome sequencing. We found that the RALY SNP is within one of the highest peaks of differentiation, which contains three genes known to influence pigmentation: ASIP, EIF2S2, and RALY (the ASIP-RALY gene block). Heterozygotes at this gene block are likely of reduced fitness, as the geographic cline of the RALY SNP has been narrow over two decades. Together, these results reflect at least one barrier to gene flow within this narrow (∼200 kb) genomic region that modulates plumage difference between species. Despite extensive gene flow between species across the genome, this study provides evidence that selection on a phenotype-associated genomic region maintains a stable species boundary.

Biased assessment of ongoing admixture using STRUCTURE in the absence of reference samples

Detection of hybridization and introgression is important in ecological research as in conservation and evolutionary biology. STRUCTURE is one of the most popular software to study introgression and allows estimating what proportion of the genome of each individual belongs to each ancestral population, even in cases where no reference sample from the ancestral nonadmixed populations is previously identified. In spite of its frequent use, some studies have indicated that ancestry estimates may not always be reliable. We simulated population data under different conditions with regard to the genetic differentiation between ancestral populations, number of loci considered, number of alleles per marker and hybridization rate, and analysed data with STRUCTURE. When reference samples were not included, the comparison of the known degree of admixture for each simulated individual and the value estimated with STRUCTURE revealed a strong underestimation of the level of introgression, classifying many admixed individuals as nonadmixed. This derives from an inaccurate estimation of the ancestral allele frequencies. When samples from the nonadmixed ancestral population were included as reference in the analyses, the bias in the estimations was reduced. The most accurate estimates were obtained when potentially admixed samples were few in relation to reference samples. Thus, whenever possible, a very large proportion of nonadmixed reference samples should be included in admixture assessments and different approaches should be combined. The misestimate of the amount of introgression can impair our understanding of the evolutionary history of species and misguide conservation efforts.

A reassessment of the diversity of green toads (Bufotes) in the circum-Aegean region

The dynamic biogeography of glacial refugia may cause complex patterns of genetic admixture between parapatric taxa, which in turn can mislead their systematics, diversity, and distributions. We investigated this issue for green toads (Bufotes) inhabiting the circum-Aegean region, a biodiversity hotspot of the Eastern Mediterranean. A previous phylogeographic study based on mitochondrial and microsatellite loci identified the hybrid zone between the European (viridis) and Anatolian (sitibundus) lineages of B. viridis all over the Balkan Peninsula, but subsequent range-wide genomic analyses (>1000 SNPs) located this transition in Turkey, a thousand kilometers eastwards. In order to clarify the diversity and taxonomy of the circum-Aegean populations, we reconciled these conflicting findings by integrating previous data with pure sitibundus individuals. Our results confirmed that the viridis/sitibundus hybrid zone extends in Western Anatolia, but that southeastern European populations feature cytonuclear discordances and a high and structured microsatellite diversity. This remarkable situation may stem from a massive geographic displacement of the hybrid zone during the last glacial fluctuations, an underappreciated yet seemingly common feature among the herpetofauna of the region. Our study thus contributes to the rising view that mitochondrial DNA can be a poor predictor of current phylogeographic structure, hence the need for genomic data, especially for narrowly distributed taxa. Finally, the analyses unambiguously support the distinction of a micro-endemic clade of green toads unique to some Cyclades islands, for which we provide a formal taxonomic description.