Natural hybridisation between Quercus petraea (Matt.) Liebl. and Quercus pubescens Willd. within an Italian stand as revealed by microsatellite fingerprinting

Interspecific gene flow is frequently reported in the genus Quercus . However, interfertile oak species often seem to remain distinct, even within areas of sympatry. This study employed molecular markers to verify, at a fine scale, the presence of interspecific gene flow in a natural population of Quercus petraea and Quercus pubescens . Within a delimited area of 6 ha, all adult trees belonging to the studied oak complex and seeds from a subsample of such trees were collected and analysed using molecular microsatellite markers. A low interspecific genetic differentiation and a high level of interspecific genetic admixture suggested past hybridisation. Paternity inference of seeds allowed the estimation of pollination frequencies from the three groups of pollen donors ( Q. petraea, Q. pubescens , intermediate). We also assayed pollen viability and germinability of each species group. We observed natural hybridisation between Q. petraea and Q. pubescens, with a predominant component in the direction Q. petraea → Q. pubescens : Q. pubescens displayed a higher level of heterospecific pollination by Q. petraea (25.8%) and intermediate morphotypes (14.7%), compared to Q. petraea acting as pollen receptor (with less than 5% heterospecific pollinations). Intermediate 'mother trees' were pollinated in similar proportions by Q. petraea (23.1%), Q. pubescens (37.8%) and intermediate morphotypes (39.1%). The asymmetrical introgression observed for the studied generation may be caused, among other factors, by the relative abundance of trees from each species group in the studied area.     

Genetic and morphological differentiation patterns between sister species: the case of Onthophagus taurus and Onthophagus illyricus (Coleoptera, Scarabaeidae)

The present study investigated the morphological and genetic differentiation pattern between two sympatric dung beetle sister species, Onthophagus taurus and Onthophagus illyricus . The geometric morphometric approach coupled with the use of molecular markers allowed examination of the nature of interspecific relationships and an outline of the evolutionary and geographical processes that might have led to interspecific differentiation and the present-day partial sympatric and syntopic pattern of distribution. Geometric morphometrics failed to discrimininate the two species on the ground of external morphological traits, but revealed interspecific differences when the shape of the primary sexual traits was taken into account. The use of two different molecular markers (cytochrome oxidase subunit I gene and amplified fragment-length polymorphism) demonstrates that the two species are genetically well differentiated, forming two distinguishable lineages probably diverged during the Pliocene by allopatric speciation. No evidence of past or recent introgression and no support for hybridization were found, suggesting that sympatry was established subsequently, when speciation was accomplished. Both molecular markers and genital shape indicate that phenotypically intermediate individuals, despite their ambiguous appearance, are members of O. illyricus species rather than hybrids.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 197–211.     

Natural hybridization and hybrid zones between Quercus crassifolia and Quercus crassipes (Fagaceae) in Mexico: morphological and molecular evidence

Hybrid zones provide interesting systems to study genetic and ecological interaction between different species. The correct identification of hybrids is necessary to understand the evolutionary process involved in hybridization. An oak species complex occurring in Mexico formed by two parental species, Quercus crassifolia H. & B. and Q. crassipes H. & B., and their putative hybrid species, Q. dysophylla, was analyzed with molecular markers (random amplified polymorphic DNA [RAPDs]) and morphological tools in seven hybrid zones (10 trees per taxa in each hybrid zone) and two pure sites for each parental species (20 trees per site). We tested whether geographic proximity of hybrid plants to the allopatric site of a parental species increases its morphological and genetic similarity with its parent. Seventeen morphological traits were measured in 8700 leaves from 290 trees. Total DNA of 250 individuals was analyzed with six diagnostic RAPD primers. Quercus crassifolia differed significantly from Q. crassipes in all the examined characters. Molecular markers and morphological characters were highly coincident and support the hypothesis of hybridization in this complex, although both species remain distinct in mixed stands. Clusters and a hybrid index (for molecular and morphological data) showed that individuals from the same parental species were more similar among themselves than to putative hybrids, indicating occasional hybridization with segregation in hybrid types or backcrossing to parents. Evidence does not indicate a unidirectional pattern of gene flow.     

Maintenance of genetic and morphological identity in two sibling Syrrhopodon species (Calymperaceae,Bryopsida) despite extensive introgression

Bryophytes are a group of land plants in which the role of hybridization has long been challenged. Using genotyping by sequencing to circumvent the lack of molecular variation at selected loci previously used for phylogeny and morphology, we determine the level of genetic and morphological divergence and reproductive isolation between the sibling Syrrhopodon annotinus and S. simmondsii (Calymperaceae, Bryopsida) that occur in sympatry but in different habitats in lowland Amazonian rainforests. A clear morphological differentiation and a low (0.06), but significant Fst derived from the analysis of 183 single nucleotide polymorphisms were observed between the two species. Conspecific pairs of individuals consistently exhibited higher average kinship coefficients along a gradient of geographic isolation than interspecific pairs. The weak, but significant genetic divergence observed is consistent with growing evidence that ecological specialization can lead to genetic differentiation among bryophyte species. Nevertheless, the spatial genetic structures of the two species were significantly correlated, as evidenced by the significant slope of the Mantel test based on kinship coefficients between pairs of interspecific individuals and the geographic distance separating them. Interspecific pairs of individuals are thus more closely related when they are geographically closer, suggesting that isolation‐by‐distance is stronger than the interspecific reproductive barrier and pointing to interspecific gene flow. We conclude that interspecific introgression, whose role has long been questioned in bryophytes, may take place even in species wherein sporophyte production is scarce due to dioicy, raising the question as to what mechanisms maintain differentiation despite weak reproductive isolation.     

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late‐successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.     

Barriers to sympatry between avian sibling species (Paridae: Baeolophus) in local secondary contact

Range limits and secondary contact zones often occur at ecotones between major associations of habitat and climate. Therefore, understanding processes that limit sympatry between species in such areas provides an important framework for testing biogeographic and evolutionary hypotheses. Theoretical and empirical work has shown that the evolution of species borders is influenced by a complexity of factors, including gene flow from central to peripheral populations and the ability of species to adapt locally to environmental conditions. However, few studies have used bioclimatic models, combined with molecular and morphological data, to predict geographic range limits in the context of gene flow across a secondary contact zone. In this study, I applied these methods to test specific hypotheses about barriers to sympatry between closely related species where they approach and contact each other. Specifically, I examined the importance of historical isolation, local adaptation, and symmetry of gene flow in limiting sympatry and range expansion of ecologically distinct species across environmental gradients. Molecular (mitochondrial DNA, allozymes), morphological, and bioclimatic data were obtained for two avian sibling species (Baeolophus inornatus and B. ridgwayi) that exist in recent, narrow secondary contact in northern California. These species are broadly allopatric and occupy rangewide associations of oak and pinyon-juniper woodlands, respectively, although B. inornatus also inhabits mixed or juniper woodlands locally. Patterns of molecular variation generally were congruent with morphological and bioclimatic data, and support prior evidence for a history of isolation, adaptation, and divergence in distinctive, species-specific vegetation-climate associations. However, molecular and morphological clines fall east of the limit of oaks, and individuals of B. inornatus in this juniper-associated contact zone experience bioclimates that are more similar to B. ridgwayi than to B. inornatus in oak habitat. Thus, B. inornatus is able to adapt and expand locally into the range of its close relative, but not vice versa. These data support the hypothesis that gene flow is asymmetrical where peripheral populations meet at range boundaries. Physiological differences between species may play an important role in influencing these patterns. Empirical studies that highlight the importance of local adaptation and patterns of gene flow in which closely related species contact across ecotones are central to understanding limits on geographic ranges, sympatry, and introgression-a cornerstone of biogeographic and speciation theory.     

Selection and dispersal in a multispecies oak hybrid zone

The four western North American red oak species (Quercus wislizeni, Q. parvula, Q. agrifolia, and Q. kelloggii) are known to produce hybrid products in all interspecific combinations. However, it is unknown whether hybrids are transitory resulting from interspecific gene flow or whether they are maintained through extrinsic selection. Here, we examine cryptic hybrid structure in Q. wislizeni through a broad region including contact and isolation from three other western North American red oaks using amplified fragment length polymorphism molecular markers. All four species were simultaneously detected in the genetic background of individuals morphologically assigned to Q. wislizeni, although the contribution of Q. kelloggii was minor. In some cases, introgression was detected well outside the region of sympatry with one of the parental species. Molecular structure at the individual level indicated this was due to long-distance pollen dispersal and not to local extinction of parental species. Species admixture proportions were correlated with climatic variables and greater proportions of Q. agrifolia and Q. parvula were present in the genetic background of Q. wislizeni in sites with cooler and more humid summers, corresponding with habitat preferences of the parental species. Partial Mantel tests indicated that climate was more important than distance from pollen source in this association. Despite high levels of introgression, species integrity was maintained in some populations in close proximity to the other species, providing further support to environmental selection in determining population genetic structure. Thus, the contribution of species mixtures to population genetic structure varies across the landscape according to availability of pollen, but more importantly to varying environmental selection pressures that produce a complex pattern of hybrid and pure gene pools.     

MATING SYSTEM AND ASYMMETRIC HYBRIDIZATION IN A MIXED STAND OF EUROPEAN OAKS

The sessile (Quercus petraea [Matt.] Liebl.) and pedunculate (Quercus robur L.) oaks are two closely related species having a wide sympatric distribution over Europe. Under natural conditions, they frequently form mixed forests, where hybridization is suspected to occur. In this paper, two different approaches have been applied to the study of the mating system and the interspecific gene flow in a mixed stand formed by the two species. The mating systems of both species have been studied separately by means of the mixed-mating model. The relative contribution of the parental species to the progenies have been estimated with two different methods. The first uses the admixture model. The second is an extension of the mixed-mating model and subdivides the outcrossing rate into intra- and interspecific components. The two species were almost completely outcrossing. This high level of outcrossing and interspecific gene flow could play an important role in the maintenance of the genetic diversity in these long-lived forest tree species. The contribution of the sessile oak to the pedunculate oak progenies varied from 17% to 48%. In contrast, ovules of sessile oak trees appear to be preferentially fertilized by other extreme sessile genotypes. We suggest that interspecific and directional gene flow was responsible for such patterns. Pedunculate oak is considered as a pioneer species and is progressively replaced by sessile oak. Our present findings add a further genetic component to this succession scheme, suggesting that unidirectional gene flow reinforces succession between the two species.     

Natural hybridisation among Quercus glabrescens,Q. rugosa and Q. obtusata (Fagaceae): Microsatellites and secondary metabolites markers

• Natural hybridisation has significant ecological, genetic and evolutionary consequences altering morphological and chemical characters of individuals. Quercus glabrescens, Q. rugosa and Q. obtusata are white oak species well separated by their morphological characters when they occur in allopatry in Mexican temperate forests. However, in sympatry, individuals with atypical morphology have been observed, suggesting hybridisation events.
• In this study, we determined, with microsatellites and secondary metabolites, if interspecific gene flow occurs when these three oak species coexist in sympatry. In total, 180 individuals belonging to seven populations [three allopatric (one for each parental species) and four sympatric sites] were analysed.
• Allopatric populations represent well‐defined genetic groups and the sympatric populations showed genetic evidence of hybridisation between Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata. The hybridisation percentage varied between sites and combination of involved species. We registered the presence of unique flavonoid compounds for Q. glabrescens (caffeic acid and flavonol 2), Q. rugosa (flavonol 5) and Q. obtusata (flavonol 1). Three compounds (quercetin rhamnoside, flavonol 3 and alkyl coumarate) were expressed in all taxa. Finally, the hybrid genotypes identified in this study (Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata) showed specific chemical profiles, resulting from a combination of those of their parental species.
• These results show that hybridisation events between these oak species alter chemical expression of secondary metabolites, creating a mosaic of resources and conditions that provide the substrate for different combinations of foliar‐associated species such as herbivores, endophytic fungi or epiphyte plants.

     

Species distinction in Irish populations of Quercus petraea and Q. robur: morphological versus molecular analyses

BACKGROUND AND AIMS: Populations of oak (Quercus petraea and Q. robur) were investigated using morphological and molecular (AFLP) analyses to assess species distinction. The study aimed to describe species distinction in Irish oak populations and to situate this in a European context. METHODS: Populations were sampled from across the range of the island of Ireland. Leaf morphological characters were analysed through clustering and ordination methods. Putative neutral molecular markers (AFLPs) were used to analyse the molecular variation. Cluster and ordination analyses were also performed on the AFLP markers in addition to calculations of genetic diversity and F-statisitcs. KEY RESULTS: A notable divergence was uncovered between the morphological and molecular analyses. The morphological analysis clearly differentiated individuals into their respective species, whereas the molecular analysis did not. Twenty species-specific AFLP markers were observed from 123 plants in 24 populations but none of these was species-diagnostic. Principal Coordinate Analysis of the AFLP data revealed a clustering, across the first two axes, of individuals according to population rather than according to species. High F(ST) values calculated from AFLP markers also indicated population differentiation (F(ST) = 0.271). Species differentiation accounted for only 13 % of the variation in diversity compared with population differentiation, which accounted for 27 %. CONCLUSIONS: The results show that neutral molecular variation is partitioned more strongly between populations than between species. Although this could indicate that the populations of Q. petraea and Q. robur studied may not be distinct species at a molecular level, it is proposed that the difficulty in distinguishing the species in Irish oak populations using AFLP markers is due to population differentiation masking species differences. This could result from non-random mating in small, fragmented woodland populations. Hybridization and introgression between the species could also have a significant role.     

Effect of a red oak species gradient on genetic structure and diversity of Quercus castanea (Fagaceae) in Mexico

Incipient reproductive barriers are a common characteristic of oak species. Disruption of these barriers promotes changes in diversity and genetic structure of the species involved. Quercus castanea is a red oak with wide geographic distribution in Mexico, which presents atypically high morphological variability when it occurs in sympatry with other red oak species, suggesting that hybridization may explain the observed variation. We tested if the genetic structure and diversity levels of Q. castanea are related to the number of red oak species growing in sympatry. In total, 14 microsatellite (SSRs) primers (six nSSRs and eight cpSSRs) were used in 120 Q. castanea individuals (20/site) belonging to six populations, where the number of red oak species associated varied from zero to five. Results showed a positive and significant relationship between the genetic diversity of Q. castanea and the number of red oak species growing in sympatry, regardless of the marker type or the parameter of genetic diversity analyzed. Also, we found a higher genetic differentiation of Q. castanea populations using cpSSRs in comparison with nSSRs. Our results suggest that temperate forests with high red oaks species richness co-dominated by Q. castanea promote the increase in this species genetic diversity. From a conservation perspective, high genetic diversity levels of foundation species such as Q. castanea may have positive cascade effects extending to other species in the community.     

Genetic variation and differentiation in population of Japanese emperor oak (Quercus dentata Thunb.) and Mongolian oak (quercus mongolica fisch. ex ledeb.) in the south of the Russian far east

Allozyme variation of Japanese emperor oak Quercus dentata Thunb. and Mongolian oak Quercus mongolica Fisch. ex Lebed.) was examined in 11 populations of these species from southern Primorye. Analysis of 18 loci showed that in these populations, 66.7% of genes of Mongolian oak are polymorphic, the number of alleles per locus being 2.28. In Japanese emperor oak, these parameters were respectively 59.8 and 2.36%. The observed heterozygosity in both species did not differ from the expected values, constituting 0.198 in Japanese emperor oak and 0.161 in Maongolian oak. The Nei's distance between Quercus mongolica and Q. dentate was 0.053. The Japanese emperor oak and Mongolian oak populations formed separate clusters in an UPGMA dendrogram. Taking into account morphological and ecological isolation of the species, we suggest that the differences in gene frequencies may reflect selection acting in favor of adaptation to different growth conditions and counteracting between-species hybridization.     

Exploring species limits in two closely related Chinese oaks

BACKGROUND: The species status of two closely related Chinese oaks, Quercus liaotungensis and Q. mongolica, has been called into question. The objective of this study was to investigate the species status and to estimate the degree of introgression between the two taxa using different approaches. METHODOLOGY/PRINCIPAL FINDINGS: Using SSR (simple sequence repeat) and AFLP (amplified fragment length polymorphism) markers, we found that interspecific genetic differentiation is significant and higher than the differentiation among populations within taxa. Bayesian clusters, principal coordinate analysis and population genetic distance trees all classified the oaks into two main groups consistent with the morphological differentiation of the two taxa rather than with geographic locations using both types of markers. Nevertheless, a few individuals in Northeast China and many individuals in North China have hybrid ancestry according to Bayesian assignment. One SSR locus and five AFLPs are significant outliers against neutral expectations in the interspecific FST simulation analysis, suggesting a role for divergent selection in differentiating species. MAIN CONCLUSIONS/SIGNIFICANCE: All results based on SSRS and AFLPS reached the same conclusion: Q. liaotungensis and Q. mongolica maintain distinct gene pools in most areas of sympatry. They should therefore be considered as discrete taxonomic units. Yet, the degree of introgression varies between the two species in different contact zones, which might be caused by different population history or by local environmental factors.     

Evidence for hybridization and introgression within a species-rich oak (<Emphasis Type="Italic">Quercus</Emphasis> spp.) community

Background  

Analysis of interspecific gene flow is crucial for the understanding of speciation processes and maintenance of species integrity. Oaks (genus Quercus, Fagaceae) are among the model species for the study of hybridization. Natural co-occurrence of four closely related oak species is a very rare case in the temperate forests of Europe. We used both morphological characters and genetic markers to characterize hybridization in a natural community situated in west-central Romania and which consists of Quercus robur, Q. petraea, Q. pubescen s, and Q. frainetto, respectively.     

A molecular phylogeny of the neotropical butterfly genus Anartia (Lepidoptera: Nymphalidae)

While Anartia butterflies have served as model organisms for research on the genetics of speciation, no phylogeny has been published to describe interspecific relationships. Here, we present a molecular phylogenetic analysis of Anartia species relationships, using both mitochondrial and nuclear genes. Analyses of both data sets confirm earlier predictions of sister species pairings based primarily on genital morphology. Yet both the mitochondrial and nuclear gene phylogenies demonstrate that Anartia jatrophae is not sister to all other Anartia species, but rather that it is sister to the Anartia fatima-Anartia amathea lineage. Traditional biogeographic explanations for speciation across the genus relied on A. jatrophae being sister to its congeners. These explanations invoked allopatric divergence of sister species pairs and multiple sympatric speciation events to explain why A. jatrophae flies alongside all its congeners. The molecular phylogenies are more consistent with lineage divergence due to vicariance, and range expansion of A. jatrophae to explain its sympatry with congeners. Further interpretations of the tree topologies also suggest how morphological evolution and eco-geographic adaptation may have set species range boundaries.     

Genetic variability and structure of Quercus brantii assessed by ISSR,IRAP and SCoT markers

Persian oak (Quercus brantii Lindl.) is one of the most important woody species of the Zagros forests in Iran. Three molecular marker techniques: start codon targeted (SCoT), inter-simple sequence repeat (ISSR) and inter-retrotransposon amplified polymorphism (IRAP) markers were compared for fingerprinting of 125 individuals of this species collected from different geographical locations of north-west of Iran. A total of 233 bands were amplified by 18 ISSR primers, of which 224 (96.10%) were polymorphic, and 126 polymorphic bands (97.65%) were observed in 129 bands amplified by 10 IRAP primers. Besides, 118 bands were observed for all 10 SCoT primers, of which 113 were polymorphic (95.71%). Average polymorphism information content (PIC) for ISSR, IRAP and SCoT markers was 0.30, 0.32 and 0.38, respectively, and this revealed that SCoT markers were more informative than IRAP and ISSR for the assessment of diversity among individuals. Based on the three different molecular types, cluster analysis revealed that 125 individuals taken for the analysis can be divided into three distinct clusters. The Jaccard's genetic similarity based on the combined data ranged from 0.23 to 0.76. These results suggest that efficiency of SCoT, IRAP and ISSR markers was relatively the same in fingerprinting of individuals. All molecular marker types revealed a low genetic differentiation among populations, indicating the possibility of gene flow between the studied populations. These results have an important implication for Persian oak (Q. brantii) germplasm characterization, improvement, and conservation.     

Co-occurrence and reproductive synchrony do not ensure hybridization between an alien tunicate and its interfertile native congener

Biological invasions can promote secondary contacts between related species and thus provide excellent case studies for investigating the joint ecological and evolutionary trajectories of interfertile taxa. Here, we studied two tunicates living in sympatry, and sometimes in syntopy, in the English Channel, Ciona intestinalis species A (presumed native to the NW Pacific) and species B (native to the N Atlantic). In addition to monitoring their co-existence over time, we examined the level of interspecific gene flow, a process that may increase the invasiveness of the non-native species. The sampling scheme was repeated twice a year for 3 years (six distinct generations) to determine the relative abundance of the two species in 11 localities along the coasts of the English Channel and the Iroise Sea (covering 1600 km) in Brittany, France. We examined 23,000 individuals, including 5315 specimens for which reproductive status was determined. Four species-diagnostic molecular markers traced interspecific gene flow on a random subset of 3048 individuals. Regardless of the sampling date, the two species co-occurred in most of the study sites, with species A showing higher frequency in the autumn. The regional pattern of seasonal variation in relative abundance of the two congeners appears to correspond to different thermal growth optima. Abrupt variations in environmental parameters can act synergistically and may favor the non-native species locally. Despite syntopy, synchronous gamete production and successful in vitro interspecific crosses, only 4.3 % individuals showed an admixed genome (i.e. footprint of present-day or past introgression events), most of them with a species A maternal lineage, of which only one was a putative first generation hybrid. Altogether, efficient barriers seem to prevent interspecific crosses between the two interfertile congeners in the wild: present-day hybridization may have less impact than competitive interactions on the fate of the two study species over their sympatric range.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

Nuclear markers,mitochondrial DNA and male secondary sexual traits variation in a newt hybrid zone (Triturus vulgaris x T. montandoni)

The smooth and the Montandon's newts (Triturus vulgaris and T. montandoni) are genetically similar sister species with highly divergent male secondary sexual traits involved in complex courtship behaviour. Their parapatric ranges overlap at moderate elevations in the Carpathian Mountains where they hybridize readily. Here we present a detailed study of genetic and morphological variation in populations from the area of sympatry. Analysis of variation at seven nuclear markers, mtDNA and male sexual secondary traits was complemented with an ecological survey of breeding sites characteristics. Extensive hybridization was revealed with back-cross individuals similar to either parental species predominating among hybrids. The hybrid zone exhibited a mosaic pattern: the genetic composition of the populations was correlated only weakly with their geographical position. No association with habitat type was found. Departures from Hardy-Weinberg proportions, significant linkage disequilibria and bimodal distribution of genotypes suggest strongly that assortative mating is an important factor shaping the genetic composition of hybrid populations. The pattern of cytonuclear disequilibria did not indicate much asymmetry in interspecific matings. Changes in the frequency of nuclear markers were highly concordant, whereas mtDNA showed much wider bidirectional introgression with 14% excess of T. montandoni haplotype. We argue that the mosaic structure of the newt hybrid zone results mainly from stochastic processes related to extinction and recolonization. Microgeographical differences in mtDNA introgression are explained by historical range shifts. Since morphologically intermediate males were underrepresented when compared to hybrid males identified by genetic markers, sexual selection acting against the morphological intermediates is implied. We discuss the implications of these findings in the context of reinforcement of prezygotic isolation in newts.