Origin and cross-century dynamics of an avian hybrid zone

Spatiotemporal dynamics in butterfly hybrid zones

Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.

Temporal stability of the hybrid zone between Calocitta magpie-jays revealed through comparison of museum specimens and iNaturalist photos

Hybrid zones are natural experiments for the study of avian evolution. Hybrid zones can be dynamic, moving as species adjust to new climates and habitats, with unknown implications for species and speciation. There are relatively few studies that have comparable modern and historic sampling to assess change in hybrid zone location and width over time, and those studies have generally found mixed results, with many hybrid zones showing change over time, but others showing stability. The white-throated magpie-jay (Calocitta formosa) and black-throated magpie-jay (Calocitta colliei) occur along the western coast of Mexico and Central America. The two species differ markedly in throat color and tail length, and prior observation suggests a narrow hybrid zone in southern Jalisco where individuals have mixed throat color. This study aims to assess the existence and temporal stability of this putative hybrid zone by comparing throat color between georeferenced historical museum specimens and modern photos from iNaturalist with precise locality information. Our results confirm the existence of a narrow hybrid zone in Jalisco, with modern throat scores gradually increasing from the parental ends of the cline toward the cline center in a sigmoidal curve characteristic of hybrid zones. Our temporal comparison suggests that the hybrid zone has not shifted its position between historical (pre-1973) and modern (post-2005) time periods-a surprising result given the grand scale of habitat change to the western Mexican lowlands during this time. An anomalous pocket of white-throated individuals in the northern range of the black-throated magpie-jay hints at the possibility of prehistorical long-distance introduction. Future genomic data will help disentangle the evolutionary history of these lineages and better characterize how secondary contact is affecting both the DNA and the phenotype of these species.

Systematics of Thraupis (Aves, Passeriformes) reveals an extensive hybrid zone between T. episcopus (Blue-gray Tanager) and T. sayaca (Sayaca Tanager)

The Neotropical avian genus Thraupis (Passeriformes, Thraupidae) currently comprises seven species that are widespread and abundant throughout their ranges. However, no phylogenetic hypothesis with comprehensive intraspecific sampling is available for the group and, therefore, currently accepted species limits remain untested. We obtained sequence data for two mitochondrial (ND2, cyt-b) and three non-coding nuclear (TGFB2, MUSK, and βF5) markers from 118 vouchered museum specimens. We conducted population structure and coalescent-based species-tree analyses using a molecular clock calibration. We integrated these results with morphometric and coloration analyses of 1,003 museum specimens to assess species limits within Thraupis. Our results confirm that Thraupis is a monophyletic group and support its origin in the late Miocene and subsequent diversification during the Pleistocene. However, we found conflicts with previous phylogenies. We recovered Thraupis glaucocolpa to be sister to all other species in the genus, and T. cyanoptera to the remaining five species. Our phylogenetic trees and population structure analyses uncovered phylogeographic structure within Thraupis episcopus that is congruent with geographic patterns of phenotypic variation and distributions of some named taxa. The first genetic and phenotypic cluster in T. episcopus occurs east of the Andes and is diagnosed by the white patch on the lesser and median wing coverts, whereas the second group has a blue patch on the wing and distributes to the west of Colombia's eastern Andes. Finally, we present evidence of hybridization and ongoing gene flow between several taxa at different taxonomic levels and discuss its taxonomic implications.

Molecular evidence and ecological niche modeling reveal an extensive hybrid zone among three Bursera species (section Bullockia)

The genus Bursera, includes ~100 shrub and trees species in tropical dry forests with its center of diversification and endemism in Mexico. Morphologically intermediate individuals have commonly been observed in Mexican Bursera in areas where closely related species coexist. These individuals are assumed to result from interspecific hybridization, but no molecular evidence has supported their hybrid origins. This study aimed to investigate the existence of interspecific hybridization among three Mexican Bursera species (Bullockia section: B. cuneata, B. palmeri and B. bipinnata) from nine populations based on DNA sequences (three nuclear and four chloroplast regions) and ecological niche modeling for three past and two future scenario projections. Results from the only two polymorphic nuclear regions (PEPC, ETS) supported the hybrid origin of morphologically intermediate individuals and revealed that B. cuneata and B. bipinnata are the parental species that are genetically closer to the putative hybrids. Ecological niche modeling accurately predicted the occurrence of putative hybrid populations and showed a potential hybrid zone extending in a larger area (74,000 km2) than previously thought. Paleo-reconstructions showed a potential hybrid zone existing from the Last Glacial Maximum (~ 21 kya) that has increased since the late Holocene to the present. Future ecological niche projections show an increment of suitability of the potential hybrid zone for 2050 and 2070 relative to the present. Hybrid zone changes responded mostly to an increase in elevational ranges. Our study provides the first insight of an extensive hybrid zone among three Mexican Bursera species based on molecular data and ecological niche modeling.

Patterns of morphological differentiation within Manacus manacus (Aves: Pipridae) in Colombia: revisiting hypotheses of isolation and secondary contact

Determining the factors responsible for phenotypic variation within species is a fundamental task in evolutionary ecology. Isolation by distance, isolation by environment and clines along secondary contact zones between formerly isolated populations are common patterns of morphological variation. In this study, we evaluated whether patterns of morphological variation exhibit association with isolation by distance, isolation by environment or secondary contact in populations of Manacus manacus with contrasting plumage colour. We used field and museum measurements of five morphological characters from 311 individuals from 88 localities distributed across a highly heterogeneous environmental range within Colombia. Climatic variables explained a higher proportion of the variation than geographic distance among localities. We found differences in wing and culmen length between white and yellow phenotypes of M. manacus in Colombia. Overall morphology did not change in a clinal fashion in relation to the contact zone and the few traits that exhibit clinal variation suggest a distant locality for the contact zone, not in line with our expectations by secondary contact between differentiated populations, suggesting alternative mechanisms for differentiation between subspecies. Our results suggest that this species exhibits marked morphological variation associated with changes in temperature and precipitation. These associations are consistent throughout the geographic range analysed, leading us to postulate that dispersal and local adaptation jointly shape the distribution of phenotypic variation.

A holocenic and dynamic hybrid zone between two cactophilic Drosophila species in a coastal lowland plain of the Brazilian Atlantic Forest

Hybridization and introgression are processes that contribute to shaping biological diversity. The factors promoting the formation of these processes are multiples but poorly explored in a biogeographical and ecological context. In the southeast coastal plain of the Brazilian Atlantic Forest, a hybrid zone was described between two closely related cactophilic species, Drosophila antonietae and D. serido. Here, we revisited and analysed specimens from this hybrid zone to evaluate its temporal and spatial dynamic. We examined allopatric and sympatric populations of the flies using independent sources of data such as mitochondrial and nuclear sequences, microsatellite loci, morphometrics of wings and male genitalia, and climatic niche models. We also verified the emergence of the flies from necrotic tissues of collected cacti to verify the role of host association for the population dynamics. Our results support the existence of a hybrid zone due to secondary contact and limited to the localities where the two species are currently in contact. Furthermore, we detected asymmetric bidirectional introgression and the maintenance of the species integrity, ecological association and morphological characters, suggesting selection and limited introgression. Considering our paleomodels, probably this hybrid zone is recent and the contact occurred during the Holocene to the present day, favoured by range expansion of their populations due to expansion of open and dry areas in eastern South America during palaeoclimatic and geomorphological events.

Distinguishing genomic homogenization from parapatric speciation in an elevationally replacing pair of Ramphocelus tanagers

Geographically connected species pairs with weakly differentiated genomes could either represent cases of genomic homogenization in progress or of incipient parapatric speciation. Discriminating between these processes is difficult because intermediate stages of either may produce weakly differentiated genomes that diverge at few locations. We used coalescent modelling applied to a genome-wide sample of SNPs to discriminate between speciation with gene flow and genomic homogenization in two phenotypically distinct but genomically weakly diverged species of elevationally replacing Ramphocelus tanagers, forming a hybrid zone in the Andean foothills. We found overwhelming support for a model of genomic homogenization following secondary contact. Simulating under this model suggested that our species pair was differentiated (F_ST  = 0.30) at secondary contact but that most of the genome has rapidly homogenized during 254 Ky of high gene flow towards the present (F_ST  = 0.02). Despite extensive genome-wide homogenization, plumage remains distinctive with a narrower than expected geographic cline width, indicating divergent selection on colour. We found two SNPs significantly associated with plumage colour, which retain moderately high F_ST . We conclude that the majority of the genome has fused, but that divergent selection on select loci probably maintains the geographically structured colour differences between these incipient species.

Microstructures amplify carotenoid plumage signals in tanagers

Brilliantly-colored birds are a model system for research into evolution and sexual selection. Red, orange, and yellow carotenoid-colored plumages have been considered honest signals of condition; however, sex differences in feather pigments and microstructures are not well understood. Here, we show that microstructures, rather than carotenoid pigments, seem to be a major driver of male-female color differences in the social, sexually-dimorphic tanager genus Ramphocelus. We comprehensively quantified feather (i) color (using spectrophotometry), (ii) pigments (using liquid chromatography-mass spectrometry (LC-MS)), and (iii) microstructures (using scanning electron microscopy (SEM) and finite-difference time-domain (FDTD) optical modeling). Males have significantly more saturated color patches than females. However, our exploratory analysis of pigments suggested that males and females have concordant carotenoid pigment profiles across all species (MCMCglmm model, female:male ratio = 0.95). Male, but not female, feathers have elaborate microstructures which amplify color appearance. Oblong, expanded feather barbs in males enhance color saturation (for the same amount of pigment) by increasing the transmission of optical power through the feather. Dihedral barbules (vertically-angled, strap-shaped barbules) in males reduce total reflectance to generate "super black" and "velvet red" plumage. Melanin in females explains some, but not all, of the male-female plumage differences. Our results suggest that a widely cited index of honesty, carotenoid pigments, cannot fully explain male appearance. We propose that males are selected to evolve amplifiers-in this case, microstructures that enhance appearance-that are not necessarily themselves linked to quality.

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.

Speciation and gene flow across an elevational gradient in New Guinea kingfishers

Closely related species with parapatric elevational ranges are ubiquitous in tropical mountains worldwide. The gradient speciation hypothesis proposes that these series are the result of in situ ecological speciation driven by divergent selection across elevation. Direct tests of this scenario have been hampered by the difficulty inferring the geographic arrangement of populations at the time of divergence. In cichlids, sticklebacks and Timema stick insects, support for ecological speciation driven by other selective pressures has come from demonstrating parallel speciation, where divergence proceeds independently across replicated environmental gradients. Here, we take advantage of the unique geography of the island of New Guinea to test for parallel gradient speciation in replicated populations of Syma kingfishers that show extremely subtle differentiation across elevation and between historically isolated mountain ranges. We find that currently described high-elevation and low-elevation species have reciprocally monophyletic gene trees and form nuclear DNA clusters, rejecting this hypothesis. However, demographic modelling suggests selection has likely maintained species boundaries in the face of gene flow following secondary contact. We compile evidence from the published literature to show that although in situ gradient speciation in labile organisms such as birds appears rare, divergent selection and post-speciation gene flow may be an underappreciated force in the origin of elevational series and tropical beta diversity along mountain slopes.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Historical climate changes and hybridization shaped the evolution of Atlantic Forest spinetails (Aves: Furnariidae)

Combining phylogeographic approaches and hybrid zone inference in a single framework is a robust way to depict respectively the biogeographic history of lineages and the evolutionary processes responsible for speciation. Here, we studied the spatiotemporal patterns of diversification and characterize the hybrid zone between two Atlantic Forest spinetails (Synallaxis ruficapilla and Synallaxis cinerea) using mitochondrial DNA and nuclear (autosomal and Z-linked) genes. We consistently recovered divergence between and within the two species during the late Pliocene and Pleistocene using an isolation with migration model. Also, our results indicate distinct levels of introgression among lineages. Ecological niche models and demographic inferences, used to infer range distributions throughout the late Quaternary, were not consistent with the hypothesis of a large river as a primary barrier responsible for the divergence of the two species. Instead, a scenario of isolation and divergence followed by geographic expansion and admixture as a consequence of Quaternary climatic oscillations was supported. Paleomodels also were not consistent with the idea that the hybrid zone originated in primary differentiation and favor a secondary contact scenario. Model fitting indicated that clines of different loci spanning the hybrid zone are coincident and concordant. The narrow cline for one Z-linked locus could be indicative of some form of post-zygotic selection hindering genetic homogenization between the two species.

Dissection by genomic and plumage variation of a geographically complex hybrid zone between two Australian non-sister parrot species, Platycercus adscitus and Platycercus eximius

The study of hybrid zones advances understanding of the speciation process, and approaches incorporating genomic data are increasingly used to draw significant conclusions about the impact of hybridisation. Despite the progress made, the complex interplay of factors that can lead to substantially variable hybridisation outcomes are still not well understood, and many systems and/or groups remain comparatively poorly studied. Our study aims to broaden the literature on avian hybrid zones, investigating a potentially geographically and temporally complex putative hybrid zone between two native Australian non-sister parrot species, the pale-headed and eastern rosellas (Platycercus adscitus and Platycercus eximius, respectively). We analysed six plumage traits and >1400 RADseq loci and detected hybrid individuals and an unexpectedly complex geographic structure. The hybrid zone is larger than previously described due to either observer bias or its movement over recent decades. It comprises different subregions where genetic and plumage signals of admixture vary markedly in their concordance. Evidence of contemporary hybridisation (later generation and backcrossed individuals) both within and beyond the previously defined zone, when coupled with a lack of F1 hybrids and differential patterns of introgression among potentially diagnostic loci, indicates a lack of post-zygotic barriers to gene flow between species. Despite ongoing gene flow, species boundaries are likely maintained largely by strong pre-mating barriers. These findings are discussed in detail and future avenues for research into this system are proposed, which would be of benefit to the speciation and hybrid zone literature.