Natural Hybrid Zones Research Articles

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group).

BackgroundHybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.ResultsWe found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.ConclusionsOur findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-015-0385-2) contains supplementary material, which is available to authorized users.

Patterns of genetic diversity and differentiation in resistance gene clusters of two hybridizing European Populus species

Resistance genes (R-genes) are essential for long-lived organisms such as forest trees, which are exposed to diverse herbivores and pathogens. In short-lived model species, R-genes have been shown to be involved in species isolation. Here, we studied more than 400 trees from two natural hybrid zones of the European Populus species Populus alba and Populus tremula for microsatellite markers located in three R-gene clusters, including one cluster situated in the incipient sex chromosome region. The results indicate that all three R-gene clusters present extensive linkage disequilibrium (LD). Outlier tests suggest balancing selection as a potential driver of R-gene diversity, but divergent selection was also detected. Nucleotide-binding site leucine-rich repeat (NBS-LRR) type R-gene clusters exhibit low species differentiation and appear to be affected by shared selection pressures between species, recurrent interspecific gene flow, or both. Our results are consistent with a role for R-gene clusters in the evolution of chromosome XIX, recently put forward as an incipient sex chromosome in Populus.

Mating system as a barrier to gene flow.

Understanding mating system as one of reproductive isolating barriers remains important although this barrier is classified in a different sense from behavioral, ecological, and mechanical isolating barriers. Selfing enhances incipient speciation while outcrossing facilitates species integrity. Here, I study how mating system affects gene exchanges between genetically diverging species in a hybrid zone. Results show that a predominant selfing species has a greater barrier to selective gene flow than does a predominant outcrossing species. Barrier to neutral gene flow convexly changes with the selfing rate due to linkage disequilibrium, with a maximum at around intermediate selfing rate. Asymmetric transient or steady-state barriers to neutral gene flow occur between two sides of a hybrid zone when the neutral gene is affected by its linked selective gene whose alternative alleles are adaptive to heterogeneous habitats. Selfing interacts with both a physical barrier and a density-dependent ecological regulation (a logarithmic model) to strengthen the barriers to neutral and selective gene flow. This theory helps to interpret incipient speciation driven by selfing or to explain the asymmetric gene flow or unequal genomic mixtures between closely related species caused by their asymmetric mating systems in natural hybrid zones.

Genome-wide mapping in a house mouse hybrid zone reveals hybrid sterility loci and Dobzhansky-Muller interactions.

Mapping hybrid defects in contact zones between incipient species can identify genomic regions contributing to reproductive isolation and reveal genetic mechanisms of speciation. The house mouse features a rare combination of sophisticated genetic tools and natural hybrid zones between subspecies. Male hybrids often show reduced fertility, a common reproductive barrier between incipient species. Laboratory crosses have identified sterility loci, but each encompasses hundreds of genes. We map genetic determinants of testis weight and testis gene expression using offspring of mice captured in a hybrid zone between M. musculus musculus and M. m. domesticus. Many generations of admixture enables high-resolution mapping of loci contributing to these sterility-related phenotypes. We identify complex interactions among sterility loci, suggesting multiple, non-independent genetic incompatibilities contribute to barriers to gene flow in the hybrid zone.

Use of a natural hybrid zone for genomewide association mapping of craniofacial traits in the house mouse.

The identification of the genes involved in morphological variation in nature is still a major challenge. Here, we explore a new approach: we combine 178 samples from a natural hybrid zone between two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus), and high coverage of the genome (~ 145K SNPs) to identify loci underlying craniofacial shape variation. Due to the long history of recombination in the hybrid zone, high mapping resolution is anticipated. The combination of genomes from subspecies allows the mapping of both, variation within subspecies and inter-subspecific differences, thereby increasing the overall amount of causal genetic variation that can be detected. Skull and mandible shape were measured using 3D landmarks and geometric morphometrics. Using principal component axes as phenotypes, and a linear mixed model accounting for genetic relatedness in the mapping populations, we identified nine genomic regions associated with skull shape and 10 with mandible shape. High mapping resolution (median size of significant regions = 148 kb) enabled identification of single or few candidate genes in most cases. Some of the genes act as regulators or modifiers of signalling pathways relevant for morphological development and bone formation, including several with known craniofacial phenotypes in mice and humans. The significant associations combined explain 13% and 7% of the skull and mandible shape variation, respectively. In addition, a positive correlation was found between chromosomal length and proportion of variation explained. Our results suggest a complex genetic architecture for shape traits and support a polygenic model.

Hybridization effects and genetic diversity of the common and black-tufted marmoset (Callithrix jacchus and Callithrix penicillata) mitochondrial control region.

Hybridization is continually documented in primates, but effects of natural and anthropogenic hybridization on biodiversity are still unclear and differentiating between these contexts remains challenging in regards to primate evolution and conservation. Here, we examine hybridization effects on the mitochondrial DNA (mtDNA) control region of Callithrix marmosets, which provide a unique glimpse into interspecific mating under distinct anthropogenic and natural conditions. DNA was sampled from 40 marmosets along a 50-km transect from a previously uncharacterized hybrid zone in NE Brazil between the ranges of Callithrix jacchus and Callithrix penicillata. DNA was also collected from 46 marmosets along a 30-km transect in a hybrid zone in Rio de Janeiro state, Brazil, where exotic marmosets appeared in the 1980s. Combining Callithrix DNA sampled inside and outside of these hybrid zones, phylogenetic and network analyses show C. jacchus and C. penicillata being parental species to sampled hybrids. We expand limited Callithrix population genetics work by describing mtDNA diversity and demographic history of these parental species. We show ancient population expansion in C. jacchus and historically constant population size in C. penicillata, with the latter being more genetically diverse than the former. The natural hybrid zone contained higher genetic diversity relative to the anthropogenic zone. While our data suggest hybrid swarm formation within the anthropogenic zone due to removed physical reproductive barriers, this pattern is not seen in the natural hybrid zone. These results suggest different genetic dynamics within natural and anthropogenic hybridization contexts that carry important implications for primate evolution and conservation.

Hybridization, introgression, and the nature of species boundaries.

Species can be defined as populations that are diagnosably distinct, reproductively isolated, cohesive, or exclusive groups of organisms. Boundaries between species in sympatry are maintained by intrinsic barriers to gene exchange; these boundaries may not be uniform in space, in time, or across the genome. Here, we explore the nature of the species boundary, defined as the phenotypes/genes/genome regions that remain differentiated in the face of potential hybridization and introgression. We emphasize that species boundaries are semipermeable, with permeability (gene exchange) being a function of genome region. The early evidence for semipermeable species boundaries came from data on differential introgression in hybrid zones. This "genic view" of species was common in the hybrid zone literature even when few molecular markers were available to characterize genome-wide patterns of variation. Now, molecular tools allow detailed characterization of differentiation between diverging lineages and patterns of variation across natural hybrid zones, but the questions being asked by evolutionary biologists have remained much the same. Recent data (from DNA sequences and genotypes) reinforce earlier conclusions about the semipermeable nature of most species boundaries. However, debate persists over the nature and extent of genome divergence that accompanies speciation.

Natural hybridization, introgression breeding, and cultivar improvement in the genus Actinidia

China is the original home of kiwifruit which are derived from the species complex Actinidia chinensis. A short domestication history of a little more than 100 years characterizes kiwifruit as a unique fruit crop whose cultivar improvement is heavily dependent on exploitation and selection of wild resources. Actinidia species are widespread geographically; their distributions often overlap, and hybridization between taxa is common naturally. Actinidia species that vary in ploidy and complex mixtures of cytotypes are prevalent in natural populations. Here, we review the recent emerging knowledge of natural distribution, biogeography, and population genetics in Actinidia with a particular focus on the pattern of overlapping distribution and natural hybridization among Actinidia species. Based on the comparison of the original geographical localities where most kiwifruit commercial cultivars were selected and the whole geographical range of the species complex A. chinensis, we propose that introgression breeding for new selections and cultivars of kiwifruit should utilize targeted exploration in natural hybrid zones as well as the genomic tools and related genetic resources that are becoming available to an unprecedented extent.

Assortative mating and the maintenance of population structure in a natural hybrid zone.

Understanding the factors that give rise to natural hybrid zones and govern their dynamics and structure is important to predicting the evolutionary consequences of hybridization. Here we use a combination of multigenerational population genetic data, mating patterns from a natural population, behavioral assays, and mark-recapture data within clinal hybrid zones of the genus Xiphophorus to test the role of assortative mating in maintaining population structure and the potential for ongoing genetic exchange between heterospecifics. Our data demonstrate that population structure is temporally robust and driven largely by assortative mating stemming from precopulatory isolation between pure species. Furthermore, mark-recapture data revealed that rates of migration within the same stream reach are far below the level needed to support population structure. In contrast to many empirical studies of natural hybrid zones, there appeared to be no hybrid male dysfunction or discrimination against hybrid males by pure parental females, and hybrid females mated and associated with pure species and hybrid males at random. Despite strong isolation between pure parentals, hybrids therefore can act as a conduit for genetic exchange between heterospecifics, which has been shown to increase the tempo of evolutionary change. Additionally, our findings highlight the complexity of natural hybrid zone dynamics, demonstrating that sexual and ecological selection together can give rise to patterns that do not fit classical models of hybrid zone evolution.

Wolbachia effects in natural populations of Chorthippus parallelus from the Pyrenean hybrid zone.

We evaluate for the first time the effect of Wolbachia infection, involving two different supergroups, on the structure and dynamics of the hybrid zone between two subspecies of Chorthippus parallelus (Orthoptera) in the Pyrenees. Wolbachia infection showed no effects on female fecundity or a slight increment in females infected by F supergroup, although in the last case it has to be well established. Cytoplasmic incompatibility (CI) is confirmed in crosses carried out in the field between individuals from a natural hybrid population. This CI, registered as the relative reduction in embryo production (sh ), was of sh = 0.355 and sh = 0.286 in unidirectional crosses involving B and F supergroups, respectively. CI also occurred in bidirectional crosses (sh = 0.147) but with a weaker intensity. The transmission rates of the two Wolbachia strains (B and F) were estimated by the optimization of a theoretical model to reach the infection frequencies observed in certain population. To fit this scenario, both supergroups should present transmission rates close to 1. Further, we have simulated the infection dynamics, and hence, the capacity of Wolbachia to structure the population of the host insects and to affect to reproduction and genetic introgression in the hybrid zone. This represents a first example of the influence of Wolbachia in an insect natural hybrid zone.

Evidence for natural hybridization betweenPrimula beesianaandP. bulleyana, two heterostylous primroses in NW Yunnan, China

AbstractNatural hybridization was assumed to play an essential role for the diversification ofPrimula; however, only one study of hybridization in the region of the Himalayas has been undertaken. In the present study, we examined another natural hybrid zone where morphologically putative hybrids as well asP. beesianaForrest,P. bulleyanaForrest, andP. poissoniiFranch. co‐occurred. We used molecular data to confirm the parental species of putative hybrids and the unidirectional hybridization pattern betweenP. beesianaandP. bulleyana. Moreover, with reference to synthetic F1s, most hybrids examined are possibly advanced generations, although the possibility of F1hybrids currently examined could not be completely excluded. In addition, pollinator observations on experimental arrays of transplanted parental species showed interspecific pollen flows during visitations of shared pollinators, indicating an incomplete pre‐zygotic barrier betweenP. beesianaandP. bulleyana. Seed productions from both flower morphs of putative hybrids were significantly lower than parental species, suggesting lower reproductive success in these hybrids. Combined with the evidence of recent habitat disturbance in the study area, we might witness the early process of hybridization betweenP. beesianaandP. bulleyana.

Reproductive phenology and physiological traits in the red mangrove hybrid complex (Rhizophora mangle and R. racemosa) across a natural gradient of nutrients and salinity

Factors modulating introgressive hybridization between the red mangrove species Rhizophora mangle and R. racemosa in spatially defined sites are poorly understood. To investigate this, we evaluated the reproductive phenology and the nutrient and physiological traits in those two species and their F1 hybrids genotyped with microsatellite data across a natural hybrid zone from the Pacific coast of Panama. We found no evidence that reproductive phenology represents a barrier to gene flow, because R. mangle and the F1 hybrids produced flowers and propagules throughout the annual cycle, while R. racemosa flowered only in the dry season. Soil nutrient concentrations decreased landward, while soil salinity varied only slightly. Foliar nutrients and δ15N signatures varied according to the soil nutrient gradient, but only foliar phosphorus and carbon varied among species. In contrast, two structural variables (height and trunk diameter) and leaf variables related to salinity tolerance (Na, Cl:Na, K:Na, cation:anion) and water-use efficiency (i.e., δ13C) differed among species, suggesting higher salinity tolerance for R. mangle and F1 hybrids compared with R. racemosa. We conclude that parental species and F1 hybrids differ in salinity tolerance and water-use efficiency, which could be associated with adaptive evolution of the red mangrove hybrid complex.

Interspecific crossing and genetic mapping reveal intrinsic genomic incompatibility between two Senecio species that form a hybrid zone on Mount Etna, Sicily

Studies of hybridizing species can reveal much about the genetic basis and maintenance of species divergence in the face of gene flow. Here we report a genetic segregation and linkage analysis conducted on F2 progeny of a reciprocal cross between Senecio aethnensis and S. chrysanthemifolius that form a hybrid zone on Mount Etna, Sicily, aimed at determining the genetic basis of intrinsic hybrid barriers between them. Significant transmission ratio distortion (TRD) was detected at 34 (∼27%) of 127 marker loci located in nine distinct clusters across seven of the ten linkage groups detected, indicating genomic incompatibility between the species. TRD at these loci could not be attributed entirely to post-zygotic selective loss of F2 individuals that failed to germinate or flower (16.7%). At four loci tests indicated that pre-zygotic events, such as meiotic drive in F1 parents or gametophytic selection, contributed to TRD. Additional tests revealed that cytonuclear incompatibility contributed to TRD at five loci, Bateson-Dobzhansky-Muller (BDM) incompatibilities involving epistatic interactions between loci contributed to TRD at four loci, and underdominance (heterozygote disadvantage) was a possible cause of TRD at one locus. Major chromosomal rearrangements were probably not a cause of interspecific incompatibility at the scale that could be examined with current map marker density. Intrinsic genomic incompatibility between S. aethnensis and S. chrysanthemifolius revealed by TRD across multiple genomic regions in early-generation hybrids is likely to impact the genetic structure of the natural hybrid zone on Mount Etna by limiting introgression and promoting divergence across the genome.

Unidirectional hybridization and reproductive barriers between two heterostylous primrose species in north-west Yunnan, China.

Heteromorphy in flowers has a profound effect on breeding patterns within a species, but little is known about how it affects reproductive barriers between species. The heterostylous genus Primula is very diverse in the Himalaya region, but hybrids there have been little researched. This study examines in detail a natural hybrid zone between P. beesiana and P. bulleyana. Chloroplast sequencing, AFLP (amplified fragment length polymorphism) markers and morphological comparisons were employed to characterize putative hybrids in the field, using synthetic F1s from hand pollination as controls. Pollinator visits to parent species and hybrids were observed in the field. Hand pollinations were conducted to compare pollen tube growth, seed production and seed viability for crosses involving different morphs, species and directions of crossing. Molecular data revealed all hybrid derivatives examined to be backcrosses of first or later generations towards P. bulleyana: all had the chloroplast DNA (cpDNA) of this species. Some individuals had morphological traits suggesting they were hybrids, but they were genetically similar to P. bulleyana; they might have been advanced generation backcrosses. Viable F1s could not be produced with P. bulleyana pollen on P. beesiana females, irrespective of the flower morphs used. Within-morph crosses for each species had very low (<10 %) seed viability, whereas crosses between pin P. bulleyana (female) and pin P. beesiana had a higher seed viability of 30 %. Thus genetic incompatibility mechanisms back up mechanical barriers to within-morph crosses in each species, but are not the same between the two species. The two species share their main pollinators, and pollinators were observed to fly between P. bulleyana and hybrids, suggesting that pollinator behaviour may not be an important isolating factor. Hybridization is strongly asymmetric, with P. bulleyana the only possible mother and all detected hybrids being backcrosses in this direction. Partial ecological isolation and inhibition of heterospecific pollen, and possibly complete barriers to F1 formation on P. beesiana, may be enough to make F1 formation very rare in these species. Therefore, with no F1 detected, this hybrid zone may have a finite life span as successive generations become more similar to P. bulleyana.

Genetic and morphological structure of a spruce hybrid (Picea sitchensis×P. glauca) zone along a climatic gradient

Historic colonization and contemporary evolutionary processes contribute to patterns of genetic variation and differentiation among populations. However, separating the respective influences of these processes remains a challenge, particularly for natural hybrid zones, where standing genetic variation may result from evolutionary processes both preceding and following contact, influencing the evolutionary trajectory of hybrid populations. Where adaptation to novel environments may be facilitated by interspecific hybridization, teasing apart these processes will have practical implications for forest management in changing environments. We evaluated the neutral genetic architecture of the Picea sitchensis (Sitka spruce) × P. glauca (white spruce) hybrid zone along the Nass and Skeena river valleys in northwestern British Columbia using chloroplast, mitochondrial, and nuclear microsatellite markers, in combination with cone morphological traits. Sitka spruce mitotype "capture", evidenced by this species dominating the maternal lineage, is consistent with earlier colonization of the region by Sitka spruce. This "capture" differs from the spatial distribution of chloroplast haplotypes, indicating pollen dispersal and its contribution to geographic structure. Genetic ancestry, based on nuclear markers, was strongly influenced by climate and geography. Highly parallel results for replicate transects along environmental gradients provide support for the bounded hybrid superiority model of hybrid zone maintenance. • This broad-scale analysis of neutral genetic structure indicates the importance of historic and contemporary gene flow, environmental selection, and their interaction in shaping neutral genetic variation within this hybrid zone, informative to seed transfer development and reforestation for future climates.

Low fertility of wild hybrid male flycatchers despite recent divergence

Postzygotic isolation may be important for maintaining species boundaries, particularly when premating barriers are incomplete. Little is known about the course of events leading from minor environmental mismatches affecting hybrid fitness to severe genetic incompatibilities causing sterility or inviability. We investigated whether reduced reproductive success of hybrid males was caused by suboptimal sperm traits or by more severe genetic incompatibilities in a hybrid zone of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) on the island of Öland, Sweden. About 4 per cent hybridization is observed in this population and all female hybrids are sterile. We found no sperm in the ejaculates of most sampled hybrid males, and sperm with abnormal morphology in two hybrids. Furthermore, none of the hybrids sired any offspring because of high levels of hatching failure and extra-pair paternity in their nests. These results from a natural hybrid zone suggest that the spermatogenesis of hybrid males may become disrupted despite little genetic divergence between the parental species.

Epidendrum(Orchidaceae) as a model system for ecological and evolutionary studies in the Neotropics

AbstractThe orchid genus Epidendrum, with 1500 species occurring within the Neotropical region, represents a very promising model system for evolutionary and ecological studies offering an expanded repertoire of research opportunities in the breadth of modern plant biology. Epidendrum displays a significant degree of morphological variation, chromosome number diversity and ecological interactions, which challenges widely held views on reproductive barriers and habitat selection. The widespread geographical distribution of many species and populations offers interesting opportunities to investigate how climatic changes and historic demographic processes shaped the current patterns of genetic and species diversity across different biomes and landscapes. Questions involving chromosome barriers to gene exchange and the role of postzygotic genetic barriers in species cohesion (e.g., the contributions of habitat selection and niche divergence on species cohesion) could be easily addressed when using the variety of natural hybrid zones found across Epidendrum. Several key evolutionary questions could be addressed with this model system, such as the identification of the first stages of adaptive radiation, the evolution of pollination strate gies, the adaptive ecological significance of trait variation and hybridisation, the influence of historical demographic events on lineage diversification and speciation. With the advance of cost–effective molecular techniques and by combining ecological and phenotypic data, researchers can now tackle these questions and foster significant progress in the field of Neotropical plant diversification and evolution.

Sperm-related phenotypes implicated in both maintenance and breakdown of a natural species barrier in the house mouse

The house mouse hybrid zone (HMHZ) is a species barrier thought to be maintained by a balance between dispersal and natural selection against hybrids. While the HMHZ is characterized by frequency discontinuities for some sex chromosome markers, there is an unexpected large-scale regional introgression of a Y chromosome across the barrier, in defiance of Haldane's rule. Recent work suggests that a major force maintaining the species barrier acts through sperm traits. Here, we test whether the Y chromosome penetration of the species barrier acts through sperm traits by assessing sperm characteristics of wild-caught males directly in a field laboratory set up in a Y introgression region of the HMHZ, later calculating the hybrid index of each male using 1401 diagnostic single nucleotide polymorphisms (SNPs). We found that both sperm count (SC) and sperm velocity were significantly reduced across the natural spectrum of hybrids. However, SC was more than rescued in the presence of the invading Y. Our results imply an asymmetric advantage for Y chromosome introgression consistent with the observed large-scale introgression. We suggest that selection on sperm-related traits probably explains a large component of patterns observed in the natural hybrid zone, including the Y chromosome penetration.

Populus hybrid hosts drive divergence in the herbivorous mite, Aceria parapopuli: implications for conservation of plant hybrid zones as essential habitat

Plant hybrid zones receive little conservation attention, yet they may be centers of diversity and evolutionary opportunity for dependent species. In previous studies, cottonwood hybrid zones have been shown to be important drivers of biological diversity and herbivore evolution. Despite these findings, no studies have examined whether hybrid host use drives herbivore genetic divergence across a broad geographic range. Here, we examined the role of Populus hybridization on the evolution of the eriophyid mite, Aceria parapopuli, using ITS1 sequence differentiation. We found support for the hypothesis that Populus hybridization has driven genetic divergence in mites in multiple hybrid zones. Furthermore, our data suggest that hybrid host use has followed at least two instances of mite genetic divergence. Our findings have several important conservation implications. First, they suggest that cottonwood hybrid zones can be important drivers of evolutionary divergence in a dependent herbivore. Second, different hybrid zones represent different ecological environments, and provide independent opportunities for local adaptation and divergence. Although hybrid plants are not considered a high priority for conservation management, and in some cases viewed as “evolutionary dead ends”, our results suggest that new consideration ought to be given to plant hybrid zones. As shown here, natural hybrid zones provide unique ecological and evolutionary opportunities, and essential habitat for dependent species, all of which deserve conservation attention and increased protection.

Genomic scan for single nucleotide polymorphisms reveals patterns of divergence and gene flow between ecologically divergent species

Recent advances in population genomics have triggered great interest in the genomic landscape of divergence in taxa with 'porous' species boundaries. One important obstable of previous studies of this topic was the low genomic coverage achieved. This issue can now be overcome by the use of 'next generation' or short-read DNA-sequencing approaches capable of assaying many thousands of single nucleotide polymorphisms (SNPs) in divergent species. We have scanned the 'porous' genomes of Populus alba and Populus tremula, two ecologically divergent hybridizing forest trees, using >38,000 SNPs assayed by restriction site associated DNA (RAD) sequencing. Windowed analyses indicate great variation in genetic divergence (e.g. the proportion of fixed SNPs) between species, and these results are unlikely to be strongly biased by genomic features of the Populus trichocarpa reference genome used for SNP calling. Divergence estimates were significantly autocorrelated (P < 0.01; Moran's I up to 0.6) along 11 of 19 chromosomes. Many of these autocorrelations involved low divergence blocks, thus suggesting that allele sharing was caused by recurrent gene flow rather than shared ancestral polymorphism. A conspicuous low divergence block of three megabases was detected on chromosome XIX, recently put forward as an incipient sex chromosome in Populus, and was largely congruent with introgression of mapped microsatellites in two natural hybrid zones (N > 400). Our results help explain the origin of the 'genomic mosaic' seen in these taxa with 'porous' genomes and suggest rampant introgression or extensive among-species conservation of an incipient plant sex chromosome. RAD sequencing holds great promise for detecting patterns of divergence and gene flow in highly divergent hybridizing species.