Reproductive character displacement allows two sexually deceptive orchids to coexist and attract the same specific pollinator

An increased divergence in characters between species in secondary contact can be shaped by selection against competition for a common resource (ecological character displacement, ECD) or against maladapted hybridization (reproductive character displacement, RCD). These selective pressures can act between incipient species (reinforcement) or well-separated species that already completed the speciation process, but that can still hybridize and produce maladapted hybrids. Here, we investigated two well-separated sexually deceptive orchid species that, unusually, share their specific pollinator. Sympatric individuals of these species are more divergent than allopatric ones in floral characters involved in a mechanical isolating barrier, a pattern suggestive of RCD. To experimentally test this scenario, we built an artificial sympatric population with allopatric individuals. We measured flower characters, genotyped the offspring in natural and artificial sympatry and estimated fertility of hybrids. Different from naturally sympatric individuals, allopatric individuals in artificial sympatry hybridized widely. Hybrids showed lower pollination success and seed viability than parentals. Character displacement did not affect plant pollination success. These findings suggest that RCD evolved between these species to avoid hybridization and that selection on reinforcement may be very strong even in plants with highly specialized pollination.


Introduction
Reproductive isolation can evolve between incipient species as a byproduct of niche divergence and may occur at different degrees of allopatry, parapatry or sympatry (Coyne and Orr 2004). When reproductive barriers have mostly evolved in allopatry, secondary contacts may represent a test for the speciation process, as interspecific gene flow can be completely or partly halted by previously formed reproductive barriers, or extensive hybridization and eventually genetic swamping can occur (Widmer et al. 2009;Zitari et al. 2012;Todesco et al. 2016).
If hybrids are less fit than parentals, natural selection is expected to favour local phenotypes that minimize gamete waste in interspecific crosses, and thus secondary contact zones become places where reproductive barriers may continue to evolve and reinforce, completing the speciation process (e.g., Gerhardt 1994;Noor 1995;Rundle and Schluter 1998;Servedio and Noor 2003;Haavie et al. 2004;Lemmon 2009;Sobel et al. 2010;Roda et al. 2017;Spriggs et al. 2019). The consequence of this so-called sympatric reinforcement of prezygotic isolation (Dobzhansky 1940), is an increased divergence of sympatric populations of the two incipient species (Noor 1999). Similar increased divergence can also occur in secondary contact zones between well-separated species that have already completed the speciation process, but that can still hybridize giving rise to sterile or partly sterile hybrids. In these latter cases, natural selection can decrease interspecific mating opportunities by displacing characters involved in reproduction, hence reducing the gamete waste when gene flow is already precluded by postzygotic isolation (Reproductive character displacement, RCD sensu Armbruster and Muchhala 2009). However, while sympatric reinforcement of prezygotic isolation (sensu Dobzhansky 1940) prevents gene flow that would lead to the genetic swamping of the incipient species, RCD relates to species that are already isolated to a degree sufficient to prevent their genetic amalgamation when in contact (Armbruster and Muchhala 2009).
Character divergence in different species may also result from selection to alleviate competition for a shared ecological resource, independently from their reproductive compatibility. In these cases, a high ecological similarity could preclude stable coexistence in secondary contact zones and local selection could optimize the use of the shared resource, a process that has been termed ecological character displacement (ECD; Brown and Wilson 1956;Pfennig and Pfennig 2009). This process is thus qualitatively different from reinforcement and RCD, though it generates similar patterns (e.g., Geyer and Palumbi 2003;Pfennig and Pfennig 2009) to the point that they have often been treated as a single process (e.g. van der Niet et al. 2006;Smith and Rausher 2008).
Under some circumstances, both reinforcing pressures (RCD and ECD) may occur on the same character when species come into contact secondarily. In these cases, it is difficult to determine whether the divergence is caused by competition for resources, or costly hybridization, as decreased hybridization is likely to be the pleiotropic consequence of selection for optimizing the use of a common resource in contact zones. A key example for this dual action of reinforcing selection are flowering plants and their pollinators. Competition for pollinator services (i.e., competition for a shared resource) may cause pollen limitation and lower seed sets (e.g., Caruso 2000;Bell et al. 2005) and natural selection will favour character divergence (morphological and/or phenological characters) between coflowering species that have met secondarily (i.e., ECD). However, if pre-zygotic isolation is not complete, the sharing of common pollinators may promote the transfer of heterospecific pollen and result in the formation of inviable or sterile hybrids, therefore lowering 1 3 plant fitness (e.g., Fishman and Wyatt 1999;Norton et al. 2015;Hopkins and Rausher 2012), and selection against hybrid mating may enhance divergence in reproductive characters (i.e., RCD). Thus, competition for pollinators represents the best example of convergence between ECD and RCD as it can be considered an indirect type of RCD or ECD as pollinators represent both a form of plant reproductive isolation and a shared ecological resource for co-flowering plant species (Dayan and Simberloff 2005). Furthermore, in order to reduce the amount of maladaptive hybridization, premating isolation is less costly than postmating isolation, as it prevents the waste of gametes (Coyne and Orr 2004;Price 2007) and, because in plants, pollinator isolation is an early premating reproductive barrier, it can be predicted as one of the main targets of reinforcing selection.
Pollinator isolation reduces competition (for pollinators) between coexisting plant species (Waser 1978;Mitchell et al. 2009), and several studies have provided evidence of ECD in which this competition was the most plausible fuel (Armbruster and Herzig 1984;Stone et al. 1998;Aizen and Vázquez 2006;Muchhala and Potts 2007). Indeed, ECD is particularly common in plant species with generalized pollination where diverged lineages, once in secondary contact, often share and compete for the overlapping pollinator sets (Norton et al. 2015 and references therein). Differently, ECD involving pollinator-related characters is less common among species with distinct, specific pollinators as competition for pollinators is unlikely upon secondary contact (Beans 2014). Indeed, pollinator-specialized plant species may coexist more easily than generalist plant species (Kozak and Wiens 2006;van der Niet et al. 2006) once they are not sharing their specialized pollinator. In contrast, generalized pollinated plants can easier displace floral characters in order to avoid pollinator interspecific movements, without the extreme consequence on fitness (Hopkins and Rausher 2012) expected for pollinator-specialised plants.
Among the pollinator-specialised plant species, orchids employing sexually deceptive strategies are considered one of most extreme examples of species-specific pollination. Each Ophrys species attracts its exclusive pollinator by producing a specific floral odour bouquet. This extremely specific pollination strategy assures ethological pollinator-mediated isolation (Ayasse et al. 2011) even among closely related sympatric sexually deceptive orchids (Xu et al. 2011;Whitehead and Peakall 2014) and the fast radiation of sexually deceptive Ophrys is driven by shift in pollinator niches (Breitkopf et al. 2015; but see also Bateman et al. 2018).
Here, we examine an unusual case in which two sexually deceptive orchid species, i.e. Ophrys normanii and O. chestermanii, belonging to different Ophrys clades (Gögler et al. 2009), converged to the same species-specific pollinator, males of the cuckoo bumblebee Bombus vestalis, by producing the same odour bouquets and overlapping in flowering time (Gögler et al. 2015). The two species are thought to have evolved in allopatry (Gögler et al. 2009;Lussu et al. 2018), but now meet in secondary contact, to the extent that the distribution of O. normanii is largely overlapping with that of O. chestermanii. In spite of phylogenetic distance, Gögler et al. (2015) found an insignificant difference in viable seed production from intra and interspecific crosses, as often reported in this orchid genus, whose reproductive barriers largely rely on pre zygotic isolation (Scopece et al. 2007). Furthermore, hybrid seeds germinated and produced viable flowering F1 hybrids suggestive of weak, if any, early postzygotic reproductive incompatibilities. Nevertheless, it remains to be tested whether other late-acting post zygotic barriers (intrinsic, such as hybrid sterility/unviability or extrinsic, such as hybrid attractiveness to pollinators) prevent genetic admixture between O. normanii and O. chestermanii. Indeed, despite the apparent absence of premating and early post zygotic reproductive barriers, the two species maintain their genetic integrity in sympatric population, hybrids being very rare in nature (Gögler et al.

3
2009). Interestingly, flower characters, as width and length of the stigmatic cavity and pollinaria length were found to be more different between sympatric than between allopatric populations of the two species (Gögler et al. 2015), a pattern suggestive of RCD (see Howard 1993). These diverging floral characters have been shown to be related to pollinator isolation in a sort of lock and key pattern shaping a mechanical barrier between the two sympatric species (Gögler et al. 2015;Lussu et al. 2019). Accordingly, sympatric individuals with floral characters impeding interspecific pollen transfer (i.e. more divergent characters) are expected to be positively selected, as opposed to allopatric individuals that experience no selection on these characters. However, if floral characters have already diverged to the point where the displaced species are no longer competing, then the signal of this reinforcing selection may also be undetectable in the contact zone. At the same time, observation of increased character divergence in sympatric populations by itself is not sufficient to demonstrate the role of these characters in species displacement as a direct link between character divergence and reproductive isolation for species pair has to be shown. Finally, character divergence may be to reduce costly hybridization (RCD) or to increase pollinator efficiency thereby reducing competition for services of the shared pollinator without shifting to a different pollinators (ECD).
To resolve these alternative hypotheses we performed a transplantation experiment by creating an artificially sympatric population where we first determined whether distributions of floral variation in our study populations are consistent across years, and consequently tested the strength of reproductive isolation between allopatric individuals and compared it with that between sympatric individuals (Martin and Harding 1981). More specifically, we addressed the following questions: 1. Is the spatial distribution of character displacement consistent with allopatric/sympatric occurrence of the two Ophrys species? 2. Are divergent floral characters correlated with pollination success (i.e., ECD)? 3. Is hybridization a cost and is it alleviated in plants with more divergent floral characters (i.e., RCD)? Is hybrid fitness lower than parental fitness? 4. Is the character divergence maintained in common garden environment and thus likely under genetic control? 5. Are allopatric individuals more likely to hybridize in artificial sympatry than natural sympatric individuals?

Study system and sample collection
Ophrys normanii Wood and O. chestermanii (Wood) Gölz and Reinhard are two rare orchid species endemic to Sardinia. Ophrys normanii distribution consists of only three small populations located in the south-western portion of the island while O. chestermanii occurs in the same area and in two other localities along the Tyrrhenian coast (Scrugli 1990). Ophrys normanii shows a slightly different ecological preference but both species are sympatric in the locality of Perdu Carta. The two species belong to two phylogenetically distant clades, with O. normanii belonging to the O. tenthredinifera clade and O. chestermanii to the O. holoserica clade (Gögler et al. 2009;Breitkopf et al. 2015) and have clearly distinguishable flower shape ( Supplementary Fig. 1).
In 2017 we collected 25 individuals of O. normanii at flowering stage (to allow unambiguous identification) from an allopatric population (Fluminimaggiore) and haphazardly placed them in a small allopatric population of O. chestermanii of approx. 30 blooming individuals (Reigraxius, Domusnovas) ( Supplementary Fig. 2) to gather a new artificial sympatric population with comparable frequencies of each species. The same source populations of O. normanii and O. chestermanii were previously investigated by Gögler et al. (2015). We then labeled all transferred O. normanii individuals, and 27 O. resident chestermanii individuals in this new artificial sympatric population. As in orchids, plant phenotype is also determined by resources allocated in the tuber produced in the previous year, only in 2018, i.e. one year after the transplantation, we collected a flower for morphological measurements from each labeled individual of O. normanii and O. chestermanii in the artificial sympatric population and in allopatric and natural sympatric populations.
At the end of the flowering season, both in 2017 and in 2018, we collected the fruits produced by each labeled individual and we also recorded pollination success in natural (Perdu Carta) and artificial (Reigraxius, Domusnovas) sympatric populations by counting the number of fruits relative to the number of flowers. At the end of the experiment, transplanted individuals were transferred back to their source population.

Hybrid fitness
To assess hybrid fitness, we used F1 hybrids between O. chestermanii and O. normanii produced from manual crosses in Gögler et al. (2015). From these crosses, three years after in vitro germination, three blooming F1 hybrid plants were obtained (two O. normanii × O. chestermanii and one O. chestermanii × O. normanii). Here, we crossed these three F1 plants (six manual crosses, each plant providing and receiving pollen from the other two) as described in Scopece et al. (2007).
In 2018, we also collected leaf material from eight labelled putative O. normanii × O. chestermanii hybrid individuals identified, on morphological bases, in another newly discovered small sympatric zone (Buggerru, loc. San Niccolò) (Fig. S2). DNA was extracted and Internal Transcribed Spacers 2 (ITS2) sequences were used to confirm potential hybrid individuals (as described below). Of the eight labeled individuals, six showed parental genetic profiles (four O. normanii and two O. chestermanii) and two showed a genetic profile compatible with the F1 hybrid (i.e. presence of heterozygote sites as described in Gögler et al. 2015). From these two putative hybrids, we collected the fruits and assessed seed viability as described in Scopece et al. (2007). With this approach we couldn't distinguish the hybrid generation of the two putative hybrids (F1, F2, BCs) nor we could identify the pollen donor of their fruits.

Morphological measurements
Flowers collected from allopatric, natural sympatric and artificial sympatric populations in 2018 were photographed under a binocular microscope (Wild Heerbrugg M8) with a digital camera (Nikon Nikkor1J1) using a millimeter scale as a reference. As previously done by Gögler et al. (2015), from each flower we measured four characters: width and height of stigmatic cavity, pollinaria length, and length of outer central tepal (a flower character potentially not involved in the placement of the pollinator on the flower) ( Supplementary Fig. 3). To measure pollinaria length we first removed pollinaria by using a measuring lens under the binocular microscope. Measurements were performed using ImageJ 1.33 software (Rasband, National Institutes of Health, USA).

Seed pool paternity test
Ripe fruits collected in natural and artificial sympatric populations in 2017 and 2018 from O. normanii and O. chestermanii individuals were harvested and stored at 4 °C until using them for molecular analysis. DNA from pooled seeds from a single fruit was extracted following Doyle and Doyle (1987). The ITS2 region was amplified by polymerase chain reaction (PCR) using primers annealing with the 3′ region of the 5.8 S and the 5′ region of the 25 S rDNA genes, by using the conditions described in Gögler et al. (2015) and purified by exonuclease I and SAP treatments (Fermentas Inc., Hannover, MD). Purified PCR products were sequenced in both directions using a modification of the Sanger dideoxy method as implemented in a double-stranded DNA cycle sequencing system with fluorescent dyes. Sequence reactions were then run on a 3130 Automated DNA Sequencer and the Sequencing Analysis software (Applied Biosystems, USA) was used to check electropherograms. Sequences were assembled with BioEdit and manually edited. ITS2 sequences showed three indels between O. chestermanii and O. normanii; hybrid fruits were thus recognized by the presence of heterozygosity at these specific sites (Gögler et al. 2015).

Data analysis
As our data were not normally distributed, we tested for differences in flower characters between O. normanii and O. chestermanii individuals in allopatry, natural sympatry and artificial sympatry using a Mann Whitney U test. In the artificial sympatric population, we also tested the differences of flower characters between O. normanii and O. chestermanii hybridizing and non-hybridizing individuals. We did not test for differences in pollinaria length between hybridizing individuals, because our experimental approach does not allow the pollen donor of observed fruits to be identified. Pollination success data were compared between natural and artificial sympatry using a Mann Whitney U test. We performed a Pearson's correlation between investigated flower characters and pollination success.
Principal component analysis (PCA) of flower character involved in reproductive isolation (i.e., width and height of stigmatic cavity, pollinaria length) was also conducted as a descriptor of overall morphological differences between species in natural and artificial sympatry. We also compared floral characters of individuals from the natural and artificial sympatric populations sampled in the present study (2018) with those of individuals sampled in the same source populations by Gögler et al. (2015).
Fruit production in manual crosses between artificial F1 hybrids was compared to fruit production in intraspecific crosses reported in Gögler et al. (2015) using a Fisher's exact test. Seed viability of natural hybrids (from Buggerru, loc. San Niccolò) was compared to seeds viability in intraspecific crosses reported in Gögler et al. (2015) using a Fisher's exact test. All analyses were conducted using SPSS (v. 21.0.0.0).

Results
Molecular analyses showed that none of the 11 fruits collected in 2017 in the natural sympatric population contained hybrid seeds, whilst 9 out of 12 (75%) fruits collected in the artificial sympatric population were of hybrid origin (5 out of 7 on O. chestermanii and 4 out of 5 on O. normanii). Similarly, in 2018, molecular analyses showed 8 out of 32 (25%) fruits collected in the artificial sympatric population were of hybrid origin (4 out of 23 fruits collected on O. chestermanii and 4 out of 9 fruits collected on O. normanii).
Overall    Gögler et al. (2015) and in the present study 1 3 Mean pollination success (calculated as number of fruits relative to the number of flowers) ranged from 0.31 to 0.53 in the two sampled years (Table 1) (Table S1).
Hybrid fitness was lower than parental fitness in terms of fruit production after manual crosses. From the 6 artificial pollinations carried out between the three F1 artificial hybrids only one developed into a fruit (i.e., 16.6%). Fruit production in crosses among F1 hybrids is thus significantly lower than that of intraspecific crosses reported in Gögler et al. (2015): 14 out of 15, i.e. 93.3%, in O. normanii and 11 out of 14, i.e. 80%, in O. chestermanii, average fruiting production 86.65%) (Fisher's exact test, p = 0.0021). In the single fruit obtained in crosses between F1 hybrids, we observed no viable seeds.
From the two putative hybrids identified by their additive genetic profile in the newly discovered sympatric zone of Buggerru (loc. San Niccolò), we collected two fruits from natural pollinated flowers. Seed viability in these two fruits was 37.1% on average. Seed viability in these natural hybrids is thus significantly lower than that of intraspecific crosses reported in Gögler et al. (2015;49.51% in O. normanii and 62.67% in O. chestermanii, average intraspecific crosses 50.09%) (Fisher's exact test, p = 0.0105). Thus, also in natural conditions, hybrid fitness was lower than parental fitness, in terms of seed viability.

Discussion
Despite a long-standing debate regarding the terminology, reinforcement may be intended as the evolution of divergence in any reproductive character in sympatric plant species in response to competition, or in response to the costs associated with sharing pollen with another species (i.e., the loss of pollen to other species, the clogging of stigmas with heterospecific pollen, or the potential for unfit hybrid seed; Pfennig and Pfennig 2009;Beans 2014). Natural selection can act in sympatric populations to strengthen reproductive isolation between hybridizing species: this may occur between incipient sister species, but also between distinct species. Here, we investigated two species belonging to different Ophrys lineages (Gögler et al. 2009;Breitkopf et al. 2015), that converged in attracting the same species-specific pollinator. The two species apparently have no premating isolation with the exclusion of morphological divergence among floral characters that produces a mechanical barrier in sympatry. The floral characters potentially involved in this barrier were found to be more differentiated in sympatry than in allopatry, suggesting the presence of reproductive character displacement (Gögler et al. 2015). While the related literature offers ample evidence and theoretical support, experimental tests validating the role of characters in RCD are less common, particularly between distant plant species, as several criteria need to be satisfied to show reinforcement (as explicitly identified in Hopkins 2013). Here, we experimentally tested the role of floral divergence in RCD by building up an artificial sympatric population with allopatric individuals and we found that in contrast to sympatric individuals, allopatric individuals widely hybridize once in contact. This finding supports the role of character displacement in generating the mechanical barrier that halts costly hybridization in these two orchids when in sympatry.
The observation of character displacement between allopatric and sympatric populations does not represent evidence for RCD, as a putative character displacement pattern could be explained by any number of other biotic and abiotic factors. A prerequisite for the evolution of RCD is that natural plant populations contain heritable variation in the displacing characters and that the difference is not (only) due to plasticity. Given the complexity of orchid life cycle and the difficulty in cultivation we cannot perform crosses for true estimates of heritability. However, we found that individuals from the same natural sympatric population sampled by Gögler et al. (2015) were not significantly different from those collected in present study and that character difference was maintained between species (with the exception of tepal length, a neutral character for pollination) (Fig. 2).
If there is a genetic basis to the difference between allopatric and sympatric individuals, then this same difference should be apparent when individuals from allopatric populations are grown together (in a sort of common garden). Accordingly, we found that character differences between populations were maintained when allopatric individuals of O. normanii were transferred in the O. chestermanii population to build the artificial sympatric zone (Fig. 1). Furthermore, by comparing character expression in the same individuals grown in allopatric and sympatric arrays, we also tested for plasticity, as the putative character displacement pattern may be caused by plastic responses to local environmental conditions. This was done across two growing years, as, in terrestrial orchids, the epigean part of the plant is renewed yearly from the small bulb-like tuber generated in the previous growing season. We found that the investigated floral characters still have a plasticity component, while maintaining the difference between allopatric and sympatric arrays across years (Fig. 2). It has been suggested that RCD may occur more readily in plastic species, because plasticity permits survival among competitors long enough for selection to narrow the reaction norm of each species in opposite directions (Beans 2014). However, we did not find that sympatric individuals have a more fixed sympatric phenotype than allopatric ones (Rice and Pfennig, 2007), as character variance was comparable across sites and years ( Figs. 1 and 2), so that gene flow, standing genetic variation, plasticity and character complexity still subtend the actual variation, in spite of the selective pressure on divergent characters in sympatric zones.
Sympatric and allopatric sites investigated here are geographically proximate (approx. 20 km, Fig. S2) and are very likely to have similar resource availability in terms of abundance of the specific pollinator. By comparing pollination success between natural and artificial sympatric areas, we found that it was slightly, but significantly, higher in natural sympatry (Table 1). However, in artificial sympatry, pollination success was not correlated to any of the investigated characters (Supplementary Table 1). This suggests that character variation does not affect plant pollination success by optimizing the pollinator resources (as expected for ECD). Interestingly, and as typically found in sexually deceptive orchids, the plants in all the sites were severely pollinator limited (i.e. a very low number of fruits compared to flower number), a condition that boosts the action of pollinator mediated selection.
A prerequisite for the existence of reproductive character displacement is that mating between diverging taxa and/or hybridization is costly. Gögler et al. (2015) found that artificial F1 hybrids can grow healthily and bloom, thus excluding early postzygotic barriers but without testing any late postzygotic barrier. Here we found that crosses between the same artificial F1 hybrids result in low fruit production and no viable seeds (i.e. late postzygotic barrier). In addition, fruits of the natural hybrids identified through a molecular approach in Buggerru (loc. San Niccolò), contained a lower amount of viable seeds when compared to parental species in same conditions. Late generation hybrids can partially restore fertility, particularly when fecundated by pollen of parental species (i.e. backcrosses). Although our approach impedes the identification of both parental individuals (i.e. the identification of pollen donor), we found that seed viability is still highly reduced when the hybrid is the ovule donor.
Increased production of maladapted hybrid progeny in artificial sympatry compared to a natural sympatric area, where character displacement occurred, confirms that the decrease in interspecific mating increases net fitness of individuals with displaced characters, and that natural selection favours divergence in phenotypes among sympatric species. The allopatric morphs produced 75% and 25% hybrid fruits in 2017 and 2018, while the sympatric morphs produced no hybrid fruits in examined years, indicating that the sympatric morphs had probably evolved to almost completely halt interspecific pollination. This is supported by the fact that no hybrid genotype was detected in the molecular analysis of Gögler et al. (2009) and that only two putative hybrid plantlets were identified in the subsequent analysis of Gögler et al. (2015) in the natural sympatric zone of Perdu Carta. On a cautionary note, however, it must be noted that the high level of interspecific pollination detected in the artificial sympatry could be partially due to the lower average distance between individuals of the two species (around 10 mt) compared to the species distribution in the natural sympatric area (around 50 mt). Shorter distance between the two species in the artificial sympatric array may have favored pollinator flights between species and thus the high interspecific pollen transfer (50% on average, i.e., as expected from random chance). We suspect that the incomplete admixture of the two species on a small scale in the natural sympatric zone allows for the potential contribution of microhabitat isolation in further reducing interspecific pollen flow.
However, in contrast to the expectation that individuals with similar phenotypes should interact more strongly, by comparing the plant phenotype and the seed progeny we found that hybridizing and non-hybridizing individuals of O. chestermanii and of O. normanii were not significantly different in terms of height of stigmatic cavity or in terms of width of stigmatic cavity (i.e. the displaced characters), so that we do not detect a direct link between the displaced character and decrease hybridization or mating between diverging taxa (Fig. 3). In sympatry, the main factor preventing gene flow between these two sympatric species is a combination of the longer pollinaria of O. normanii and the smaller and differently shaped stigmatic cavity in O. chestermanii, a sort of lock and key mechanism as clearly confirmed by video recording (Gögler et al. 2015). We argue that this is an intrinsic limitation of our experimental design as, by analyzing fruits, we can only measure the displaced characters in the ovule donor plant, but not in the pollen donor plant, so that we cannot fully evaluate the effect of the characters combination.
Our study on these rare and endemic orchids has evident limits in the replications of allopatric and sympatric zones and of in the small number of examined artificial and natural hybrids. Still, our results clearly indicate that the character displacement detected in the area of secondary contact can be considered proof of reinforcing selection against maladaptive hybridization (e.g., RCD), rather than ECD resulting from competition for pollinator service (see above). Differently from generalist species, increased floral divergence in highly specialized plants, which eventually come into secondary contact and share their specialized pollinator, is expected to be more strongly constrained by the specific requirements of the pollinator. Indeed, both the odor bouquet and phenology are not displaced between species in sympatry, as they may potentially have detrimental effects on plant pollination success (Smith and Rausher 2008;Gögler et al. 2015).This, coupled with a condition of strong pollinator limitation, as typical of deceptive orchids, suggests that the reinforcing pressure can be very strong, so that we detected phenotypic divergence between sympatric and allopatric populations, even in a very narrow geographic context (less than 100 km 2 ) and, most likely, in the presence of intraspecific gene flow between allopatric and sympatric populations.
Compared to generalist pollinated plants that may vary in reinforcing processes to avoid overlap in pollinators in contact zones (as shifting in phenology, changes in pollinator spectrum, scent, colour, etc.) highly specialist pollinated plants are expected to be severely constrained by the very strict interaction with their unique pollinator. Nevertheless, here we report on a case of reinforcement of reproductive isolation upon secondary contact occurring, in plants, like orchids, with extremely specialized pollination systems.