Admixture in European Populus hybrid zones makes feasible the mapping of loci that contribute to reproductive isolation and trait differences

The use of admixed human populations to scan the genome for chromosomal segments affecting complex phenotypic traits has proved a powerful analytical tool. However, its potential in other organisms has not yet been evaluated. Here, we use DNA microsatellites to assess the feasibility of this approach in hybrid zones between two members of the 'model tree' genus Populus: Populus alba (white poplar) and Populus tremula (European aspen). We analyzed samples of both species and a Central European hybrid zone (N=544 chromosomes) for a genome-wide set of 19 polymorphic DNA microsatellites. Our results indicate that allele frequency differentials between the two species are substantial (mean delta=0.619+/-0.067). Background linkage disequilibrium (LD) in samples of the parental gene pools is moderate and should respond to sampling schemes that minimize drift and account for rare alleles. LD in hybrids decays with increasing number of backcross generations as expected from theory and approaches background levels of the parental gene pools in advanced generation backcrosses. Introgression from P. tremula into P. alba varies strongly across marker loci. For several markers, alleles from P. tremula are slightly over-represented relative to neutral expectations, whereas a single locus exhibits evidence of selection against P. tremula genotypes. We interpret our results in terms of the potential for admixture mapping in these two ecologically divergent Populus species, and we validate a modified approach of studying genotypic clines in 'mosaic' hybrid zones.     

Clonality and spatial genetic structure in Populus x canescens and its sympatric backcross parent P. alba in a Central European hybrid zone.

Spatial genetic structure (SGS) holds the key to understanding the role of clonality in hybrid persistence, but multilocus SGS in hybrid zones has rarely been quantified. Here, the aim was to fill this gap for natural hybrids between two diploid, ecologically divergent European tree species with mixed sexual/asexual reproduction, Populus alba and P. tremula. Nuclear microsatellites were used to quantify clonality, SGS, and historical gene dispersal distances in up to 407 trees from an extensive Central European hybrid zone including three subpopulation replicates. The focus was on P. x canescens and its backcross parent P. alba, as these two genotypic classes co-occur and interact directly. Sexual recombination in both taxa was more prominent than previously thought, but P. x canescens hybrids tended to build larger clones extending over larger areas than P. alba. The 3.4 times stronger SGS in the P. x canescens genet population was best explained by a combination of interspecific gene flow, assortative mating, and increased clonality in hybrids. Clonality potentially contributes to the maintenance of hybrid zones of P. alba and P. tremula in time and space. Both clonality and SGS need to be taken into account explicitly when designing population genomics studies of locus-specific effects in hybrid zones.     

Genetic analysis of post-mating reproductive barriers in hybridizing European Populus species

Molecular genetic analyses of experimental crosses provide important information on the strength and nature of post-mating barriers to gene exchange between divergent populations, which are topics of great interest to evolutionary geneticists and breeders. Although not a trivial task in long-lived organisms such as trees, experimental interspecific recombinants can sometimes be created through controlled crosses involving natural F(1)'s. Here, we used this approach to understand the genetics of post-mating isolation and barriers to introgression in Populus alba and Populus tremula, two ecologically divergent, hybridizing forest trees. We studied 86 interspecific backcross (BC(1)) progeny and >350 individuals from natural populations of these species for up to 98 nuclear genetic markers, including microsatellites, indels and single nucleotide polymorphisms, and inferred the origin of the cytoplasm of the cross with plastid DNA. Genetic analysis of the BC(1) revealed extensive segregation distortions on six chromosomes, and >90% of these (12 out of 13) favored P. tremula donor alleles in the heterospecific genomic background. Since selection was documented during early diploid stages of the progeny, this surprising result was attributed to epistasis, cyto-nuclear coadaptation, heterozygote advantage at nuclear loci experiencing introgression or a combination of these. Our results indicate that gene flow across 'porous' species barriers affects these poplars and aspens beyond neutral, Mendelian expectations and suggests the mechanisms responsible. Contrary to expectations, the Populus sex determination region is not protected from introgression. Understanding the population dynamics of the Populus sex determination region will require tests based on natural interspecific hybrid zones.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Genetic structure of hybrid zones between Pinus pumila and P. parviflora var. pentaphylla (Pinaceae) revealed by molecular hybrid index analysis

Pinus species have three differently inherited genomes: paternal chloroplast, maternal mitochondrial, and biparental nuclear. Our previous study on the hybrid zones between alpine Pinus pumila and montane to subalpine P. parviflora var. pentaphylla demonstrated contrasting patterns of introgression of two cytoplasmic genomes, i.e., the paternal cpDNA flowed from P. parviflora var. pentaphylla to P. pumila, and the maternal mtDNA flowed in the reverse direction. In the present study, we developed codominant nuclear DNA markers diagnostic or mostly diagnostic for each parental species by single-strand conformation polymorphism (SSCP) of polymerase chain reaction (PCR) products, using expressed sequence tag (EST) primers of Pinus taeda. To describe the introgressive patterns of the nuclear genes, the molecular hybrid index (MHI) showing the overall proportion of alleles inferred to be derived from P. pumila was determined for each plant collected in hybrid zones on Mt. Asahidake and Mt. Higashiazuma, Japan. At Mt. Asahidake, the MHI values changed clinally according to the altitudes at which the plants were collected. However, at Mt. Higashiazuma, there was a gap in the MHI values between the plants above and below the Abies and Tsuga forest zone (alt. 1800-1900 m). This suggested that the zone plays a role in creating an effective barrier to gene flow in the hybrid zone.     

Fitness dynamics within a poplar hybrid zone: I. Prezygotic and postzygotic barriers impacting a native poplar hybrid stand

Hybridization and introgression are pervasive evolutionary phenomena that provide insight into the selective forces that maintain species boundaries, permit gene flow, and control the direction of evolutionary change. Poplar trees (Populus L.) are well known for their ability to form viable hybrids and maintain their distinct species boundaries despite this interspecific gene flow. We sought to quantify the hybridization dynamics and postzygotic fitness within a hybrid stand of balsam poplar (Populus balsamifera L.), eastern cottonwood (P. deltoides Marsh.), and their natural hybrids to gain insight into the barriers maintaining this stable hybrid zone. We observed asymmetrical hybrid formation with P. deltoides acting as the seed parent, but with subsequent introgression biased toward P. balsamifera. Native hybrids expressed fitness traits intermediate to the parental species and were not universally unfit. That said, native hybrid seedlings were absent from the seedling population, which may indicate additional selective pressures controlling their recruitment. It is imperative that we understand the selective forces maintaining this native hybrid zone in order to quantify the impact of exotic poplar hybrids on this native system.     

Matrix attachment regions (MARs) enhance transformation frequency and transgene expression in poplar

We tested the value of a matrix attachment region (MAR) fragment derived from a tobacco gene for increasing the frequency of Agrobacterium-mediated transformation. A binary vector that carried a GUS reporter gene containing an intron and an nptII gene was modified to contain flanking MAR elements within the T-DNA borders. Vectors containing or lacking MARs were then used to transform tobacco, a readily transformabl e poplar clone (Populus tremula × P. alba), and a recalcitrant poplar clone (Populus trichocarpa × P. deltoides). MARs increased GUS gene expression approximately 10-fold in the two hybrid poplar clones and twofold in tobacco one month after cocultivation with Agrobacterium; MARs also increased the frequency of kanamycin-resistant poplar shoots recovered     

Genetic analysis of admixture and patterns of introgression in foundation cottonwood trees (Salicaceae) in southwestern Colorado,USA

Cottonwoods are well known as foundation riparian trees that support diverse communities and drive ecosystem processes. Although hybridization naturally occurs when the distributions of two or more cottonwood species overlap, few cottonwood hybrid zones have been genetically characterized. We use genetic and genomic analyses to characterize patterns of admixture and introgression for a newly described hybrid zone at the intersection of three species (Populus L. Salicaceae—Populus deltoides, Populus fremontii, and Populus angustifolia) in southwestern Colorado, USA. Analysis of nuclear and chloroplast microsatellite marker data detected substantial genetic variation among individuals, revealing that (1) hybridization is occurring between two, not three, species (P. deltoides and P. angustifolia); (2) gene flow is bidirectional; (3) hybrids are not abundant (admixture detected in only 34 of 270 trees), with most being early-generation F₁ hybrids; (4) cytonuclear disequilibria exists and F₁ hybrids tend to retain P. deltoides—like chloroplasts; and (5) roughly 30 % of the nuclear markers deviated from a neutral pattern of introgression, suggesting that selection may play a role in shaping the genetic structure of the hybrid zone in this region. Overall, our results show that despite strong selection maintaining species divergence, transfer of allelic variants across species boundaries can occur. Our study assesses the fine-scale genetic structure of hybridization between P. angustifolia and P. deltoides and lays the foundation for examining how geographic differences in hybrid zone dynamics arise and may influence subsequent ecological and evolutionary processes.     

Spatial distribution of cpDNA and mtDNA haplotypes in a hybrid zone betweenPinus pumila andP. parviflora var.pentaphylla (Pinaceae)

Pinus species exhibit paternal chloroplast inheritance and maternal mitochondrial inheritance. This independent inheritance of two cytoplasmic genomes provides an exceptional environment for discriminating female (seeds) and male (pollen) components of gene flow across hybridizing species. We obtained mitochondrial genetic markers diagnostic toP. parviflora var.pentaphylla andP. pumila by PCR amplification of the intron ofnad1 on mtDNA, and examined the spatial-distribution pattern of the mtDNA haplotypes in a hybrid zone betweenP. parviflora var.pentaphylla andP. pumila in the Tanigawa Mountains of Japan. These data, in conjunction with previous information on cpDNA haplotypes and needle morphology, revealed contrastive patterns of introgression of two cytoplasmic genomes. CpDNA introgression has occurred uni-directionally fromP. parviflora var.pentaphylla toP. pumila. Conversely, mtDNA introgression has occurred in the opposite direction, fromP. pumila toP. parviflora var.pentaphylla. Levels of introgression are roughly equivalent for cpDNA and mtDNA. The contrastive spatial distribution pattern of cpDNA and mtDNA haplotypes could be caused by differential movement of seeds and pollen for interspecific genetic exchange.     

Molecular and morphological patterns of introgression between two large white-headed gull species in a zone of recent secondary contact

Incomplete reproductive isolation promotes gene flow between diverging taxa. However, any gene encoding for traits involved in the reproductive barriers will be less prone to introgression than neutral markers. Comparing introgression rates among loci is thus informative of the number and functions of loci involved in the reproductive barriers. This study aimed at identifying possible mechanisms of restriction to gene flow across a zone of recent secondary contact between Larus argentatus and Larus cachinnans by comparing introgression patterns for nine microsatellite loci, a fragment of mitochondrial DNA and a set of phenotypic traits. The low linkage disequilibrium between neutral nuclear markers indicated introgression without any barrier to gene flow. However, asymmetric introgression of mitochondrial DNA suggested that interspecific crosses may be more successful in one direction. The introgression rate for phenotypic traits was variable and low compared to neutral molecular markers. This was particularly evident in colouration of bare parts: individuals with intermediate colouration were scarcer in sympatry than expected if the genomes recombined freely. We hypothesized that one of these variables, the orbital ring colour, may play a role in mate choice, acting as an incomplete premating barrier through assortative mating. This study emphasizes that multilocus approaches are useful to discriminate among possible mechanisms responsible for the maintenance of hybrid zones.     

Hybrid bridges to gene flow: a case study in milkweeds (Asclepias)

Natural hybridization occurs throughout areas of sympatry for the North American milkweeds Asclepias exaltata and A. syriaca (Asclepiadaceae), even though the formation of F1 hybrid seed is a rare event. For introgressive hybridization to proceed, F1 and advanced hybrids must be released from reproductive barriers and successfully mate with one or both parental species. I investigated the mating system of natural hybrids between A. exaltata and A. syriaca in three populations in Shenandoah National Park, Virginia. Allozyme data and a maximum-likelihood procedure were used to estimate the frequency of six genotypic classes (parentals, F1, F2, and backcrosses) of the hybridizing populations, the pollinia received by hybrid plants, and the paternal parents of seeds produced by hybrids. F1 hybrids, backcross A. syriaca, and parental A. syriaca individuals were common in three hybrid populations. Even though self-pollinations and interhybrid pollinations were common, F2 seed production and the occurrence of F2 individuals were rare in hybrid populations. Hybrid plants received more pollen from A. syriaca than A. exaltata, which resulted in the production of more backcross-A. syriaca seed than backcross-A. exaltata seed. Asclepias exaltata was rare in the hybrid populations, but A. exaltata pollinia were received by hybrids and this species sired between 15% and 36% of the seeds produced on hybrids. The potential for introgression with A. exaltata populations is lower because this species is unsuccessful as the maternal parent in interspecific and backcross hand-pollinations. The asymetry of hybridization with A. syriaca as the maternal parent is further supported by the incorporation of maternally inherited chloroplast DNA markers in hybrids. Hybrid milkweeds frequently backcross with both parental species and may be released from the reproductive barriers that limit the formation of F1 hybrids in natural populations. The direction of interspecific gene flow and introgression in milkweeds is influenced by the reproductive biology of hybrids, the constituency of the surrounding population, and failure of some crosses to produce seeds. Finally, introgressive hybridization remains an important evolutionary force even when the initial formation of F1 hybrids in natural populations is rare.     

Regulation of glutathione synthesis in leaves of transgenic poplar (Populus tremula X P. alba) overexpressing glutathione synthetase

The poplar hybrid Populus tremula X P. alba was transformed with the Escherichia coli gene for glutathione synthetase ( gsh II ) targetted to the cytosol. Leaves of five lines of transgenic plants exhibited glutathione synthetase activities 15- to 60-fold higher than those of wild-type plants. Total glutathione levels and GSH/GSSG ratios were similar in transgenic and wild-type plants. Precursor feeding experiments with cysteine and γ-glutamylcysteine suggest that glutathione synthesis in the cytoplasm is controlled by a multistep procedure that includes (i) the availability of cysteine, (ii) the availability of γ-glutamylcysteine, and (iii) regulation of the activities of both γ-glutamylcysteine synthetase and glutathione synthetase. However step (ii) may set an upper limit for the cellular glutathione content.     

Analysis of the clinal variation pattern inSilene latifolia (Caryophyllaceae) pollen morphology

Reticulate and microechinate pollen forms withinSilene latifolia (S. alba, S. pratensis) were first described from N. America, where the species is an introduced weed. A previous study showed that the two forms also exist in Europe and intergrade along a zone roughly congruent with intergradation zones in seed morphology and flavone glycosylation genotypes. The present survey of pollen from herbarium specimens is more extensive and covers localities from nearly the whole of the species' native Eurasian and North African range. The first axis of a detrended correspondence analysis (DCA) of qualitative characterstate data for 11 pollen morphological characters shows a gradient between the two extreme exine types, reticulate and microechinate. When mapped, the co-ordinates of pollen samples on this axis indicate a broad zone of intergradation which coincides approximately with the somewhat sharper transition between low- and high-tubercle seeds in C. and N. Europe, but becomes diffuse in Italy and the Mediterranean region and diverges from the seed transition zone in SW. Asia and the Middle East. The biological significance of the pollen morphs is unknown; the variation pattern in pollen morphology, unlike that in seeds, is not consistently correlated with macroclimate.     

The strength of reproductive isolation in two hybridizing food-deceptive orchid species

Reproductive isolation is of fundamental importance for maintaining species boundaries in sympatry. In orchids, the wide variety of pollination systems and highly diverse floral traits have traditionally suggested a prominent role for pollinator isolation, and thus for prezygotic isolation, as an effective barrier to gene flow among species. Here, we examined the nature of reproductive isolation between Anacamptis morio and Anacamptis papilionacea, two sister species of Mediterranean food-deceptive orchids, in two natural hybrid zones. Comparative analyses of the two hybrid zones that are located on soils with volcanic origin and have different and well-dated ages consistently revealed that all hybrid individuals were morphologically and genetically intermediate between the parental species, but had strongly reduced fitness. Molecular analyses based on nuclear ITS1 and (amplified fragment length polymorphism) AFLP markers clearly showed that all examined hybrids were F1 hybrids, and that no introgression occurred between parental species. The maternally inherited plastid DNA markers indicated that hybridization between A. morio and A. papilionacea was bidirectional, as confirmed by the molecular analysis of seed families. The genetic architecture of the two hybrid zones suggests that the two parental species easily and frequently hybridize in sympatry as a consequence of partial pollinator overlap but that strong postzygotic barriers reduce hybrid fitness and prevent gene introgression. These results corroborate that chromosomal divergence is instrumental for reproductive isolation between these food-deceptive orchids and suggest that hybridization is of limited importance for their diversification.     

Direction of movements in Hungarian Barn Owls (Tyto alba): gene flow and barriers

Abstract. An analysis of dispersal directions of the barn owl showed that all individuals immigrating to Hungary came from W‐NW‐N. It was shown that immigrating owls breed in Hungary. There is no prevailing direction in emigration from Hungary. The time of fledging does not influence the direction of movement and there is no difference between sexes concerning dispersal direction. The percentages of emigrating owls is greater than that of immigrating ones. These percentages did not differ in relation to most of the analysed countries (Germany, Italy, Switzerland, Poland and countries of the former Yugoslavia and Czechoslovakia) but it differed in relation to Austria. The degree and direction of introgression into and from the transition zone and the recent distribution of the phenotypes are discussed based on the comparative analysis of published European data. These suggest that the subspecies Tyto alba alba and Tyto alba guttata disappear by introgression, to form a phenotypically very variable species.     

GENOMIC AND MORPHOLOGICAL ANALYSIS OF A SEMIPERMEABLE AVIAN HYBRID ZONE SUGGESTS ASYMMETRICAL INTROGRESSION OF A SEXUAL SIGNAL

Hybrid zones are geographic regions where differentiated taxa meet and potentially exchange genes. Increasingly, genomic analyses have demonstrated that many hybrid zones are semipermeable boundaries across which introgression is highly variable. In some cases, certain alleles penetrate across the hybrid zone in only one direction, recombining into the alternate genome. We investigated this phenomenon using genomic (genotyping‐by‐sequencing) and morphological (plumage reflectance spectrophotometry) analyses of the hybrid zone between two subspecies of the red‐backed fairy‐wren (Malurus melanocephalus) that differ conspicuously in a sexual signal, male back plumage color. Geographic cline analyses revealed a highly variable pattern of differential introgression, with many narrow coincident clines combined with several significantly wider clines, suggesting that the hybrid zone is a semipermeable tension zone. The plumage cline was shifted significantly into the genomic background of the orange subspecies, consistent with sexual selection driving asymmetrical introgression of red plumage alleles across the hybrid zone. This interpretation is supported by previous experimental work demonstrating an extra‐pair mating advantage for red males, but the role of genetic dominance in driving this pattern remains unclear. This study highlights the potential for sexual selection to erode taxonomic boundaries and promote gene flow, particularly at an intermediate stage of divergence.     

Interspecific mating in the Piriqueta caroliniana (Turneraceae) complex: effects of pollen load size and composition

Two taxa of Piriqueta (P. caroliniana and P. viridis) form a broad hybrid zone that extends over much of the central Florida peninsula. We used genetic markers to examine the strength of the isolation barriers between these taxa and the patterns of mating at the initial stages of hybridization. Regression models were employed to analyze the effects of pollen load size and the proportions of intra- and interspecific pollen on the frequency of first-generation (F1) hybrid formation. Overall, the postpollination mating barriers between these two taxa were relatively weak. However, there were significant effects of pollen load composition and size on the patterns of hybridization in both taxa with frequency-dependent responses to composition in both taxa. The frequency of F1 hybrid formation was generally lower than expected based on the frequency of each pollen type on the stigma for P. caroliniana recipients. The lower frequencies of F1 seeds in this taxon were apparently due to a greater competitive ability for intraspecific pollen, since hybrid seed formation decreased with increasing pollen load size. Pollen from P. caroliniana donors was also competitively superior on P. viridis recipients, leading to higher than expected frequencies of hybrid seed formation. Pollen from P. caroliniana did suffer higher rates of pollen-tube attrition than intraspecific pollen on P. viridis recipients, so the frequency of hybrid seed formation would be lower when pellen load sizes were small. In general, reproductive isolation mechanisms were stronger in P. caroliniana, suggesting that introgression should occur into P. viridis when these taxa come into close contact. Comparison of these expected patterns of mating to the distribution of hybrid genotypes in Florida provides insights into the relative roles of mating and selection in the evolution of hybrid populations of Piriqueta.     

Comparative genome mapping among Populus adenopoda, P. alba, P. deltoides, P. euramericana and P. trichocarpa

Among the genus Populus, the sections Populus (white poplar), Aigeiros Duby (black poplar) and Tacamahaca Spach contain many tree species of economical and ecological important properties. Two parental maps for the inter-specific hybrid population of Populus adenopoda × P. alba (two species of Populus section) were constructed based on SSR and SRAP markers by means of a two-way pseudo-test cross mapping strategy. The same set of SSR markers developed from the P. trichocarpa (belonging to Tacamahaca section) genome which were used to construct the maps of P. deltoides and P. euramericana (two species of Aigeiros section) was chosen to analyze the genotype of the experimental population of P. adenopoda × P. alba. Using the mapped SSR markers as allelic bridges, the alignment of the white and black poplar maps to each other and to the P. trichocarpa physical map was conducted. The alignment showed high degree of marker synteny and colinearity and the closer relationship between Aigeiros and Tacamahaca sections than that of Populus and Tacamahaca. Moreover, there was evidence for the chromosomal duplication and inter-chromosomal reorganization involving some poplar linkage groups, suggesting a complicated course of fission or fusion in one of the lineages. A poplar consensus map based on the comparisons could be constructed will be useful in practical applications including marker assisted selection.     

A set of novel DNA polymorphisms within candidate genes potentially involved in ecological divergence between Populus alba and P. tremula, two hybridizing European forest trees

We have identified 53 DNA (single nucleotide, microsatellite, and insertion-deletion) polymorphisms within 12 candidate genes potentially involved in ecological differences between Populus alba and P. tremula, two hybridizing European forest trees. The genes represent candidates for functional roles associated with abiotic or biotic stress response, cross-talk, phenology, and leaf development. Distributions within sequences, intraspecific levels of diversity, and genetic divergence (F(ST) ) between species are reported for each polymorphism, as are haplotype frequencies for each gene. The markers will be used for population genomic studies of the barrier to gene flow between these two ecologically divergent forest trees.     

Scale and direction of adaptive introgression between black cottonwood (Populus trichocarpa) and balsam poplar (P. balsamifera)

Introgression can introduce novel genetic variation at a faster rate than mutation alone and result in adaptive introgression when adaptive alleles are maintained in the recipient genome over time by natural selection. A previous study from our group demonstrated adaptive introgression from Populus balsamifera into P. trichocarpa in a target genomic region. Here we expand our local ancestry analysis to the whole genome of both parents to provide a comprehensive view of introgression patterns and to identify additional candidate regions for adaptive introgression genomewide. Populus trichocarpa is a large, fast‐growing tree of mild coastal regions of the Pacific Northwest, whereas P. balsamifera is a smaller stature tree of continental and boreal regions with intense winter cold. The species hybridize where they are parapatric. We detected asymmetric patterns of introgression across the whole genome of these two poplar species adapted to contrasting environments, with stronger introgression from P. balsamifera to P. trichocarpa than vice versa. Admixed P. trichocarpa individuals contained more genomic regions with unusually high levels of introgression (19 regions) and also the largest introgressed genome fragment (1.02 Mb) compared with admixed P. balsamifera (nine regions). Our analysis also revealed numerous candidate regions for adaptive introgression with strong signals of selection, notably related to disease resistance, and enriched for genes that may play crucial roles in survival and adaptation. Furthermore, we detected a potential overrepresentation of subtelomeric regions in P. balsamifera introgressed into P. trichocarpa and possible protection of sex‐determining regions from interspecific gene flow.