Evolutionarily inspired solutions to the crop pollination crisis.

The global decline in insect diversity threatens pollination services, potentially impacting crop production and food security. Here, we argue that this looming pollination crisis is generally approached from an ecological standpoint, and that consideration of evolutionary principles offers a novel perspective. First, we outline that wild plant species have overcome 'pollination crises' throughout evolutionary history, and show how associated principles can be applied to crop pollination. We then highlight technological advances that can be used to adapt crop flowers for optimal pollination by local wild pollinators, especially by increasing generalization in pollination systems. Thus, synergies among fundamental evolutionary research, genetic engineering, and agro-ecological science provide a promising template for addressing a potential pollination crisis, complementing much-needed strategies focused on pollinator conservation.

Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families.

Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION® platform, with and without the use of a short fragment depletion kit.

Finding New Communities: A Principle of Neuronal Network Reorganization in Alzheimer's Disease.

Background: Graph-theoretical analyses have been previously used to investigate changes in the functional connectome in patients with Alzheimer's disease (AD). However, these analyses generally assume static organizational principles, thereby neglecting a fundamental reconfiguration of functional connections in the face of neurodegeneration. Methods: Here, we focus on differences in the community structure of the functional connectome in young and old individuals and patients with AD. Patients with AD, moreover, underwent molecular imaging positron emission tomography by using [18F]AV1451 to measure tau burden, a major hallmark of AD. Results: Although the overall organizational principles of the community structure of the human functional connectome were preserved even in advanced healthy aging, they were considerably changed in AD. We discovered that the communities in AD are re-organized, with nodes changing their allegiance to communities, thus resulting in an overall less efficient re-organized community structure. We further discovered that nodes with a tendency to leave the communities displayed a relatively higher tau pathology burden. Discussion: Together, this study suggests that local tau pathology in AD is associated to fundamental changes in basic organizational principles of the human connectome. Our results shed new light on previous findings obtained by using the graph theory in AD and imply a general principle of the brain in response to neurodegeneration.

Diversification in evolutionary arenas-Assessment and synthesis.

Understanding how and why rates of evolutionary diversification vary is a key issue in evolutionary biology, ecology, and biogeography. Evolutionary rates are the net result of interacting processes summarized under concepts such as adaptive radiation and evolutionary stasis. Here, we review the central concepts in the evolutionary diversification literature and synthesize these into a simple, general framework for studying rates of diversification and quantifying their underlying dynamics, which can be applied across clades and regions, and across spatial and temporal scales. Our framework describes the diversification rate (d) as a function of the abiotic environment (a), the biotic environment (b), and clade-specific phenotypes or traits (c); thus, d ~ a,b,c. We refer to the four components (a-d) and their interactions collectively as the "Evolutionary Arena." We outline analytical approaches to this framework and present a case study on conifers, for which we parameterize the general model. We also discuss three conceptual examples: the Lupinus radiation in the Andes in the context of emerging ecological opportunity and fluctuating connectivity due to climatic oscillations; oceanic island radiations in the context of island formation and erosion; and biotically driven radiations of the Mediterranean orchid genus Ophrys. The results of the conifer case study are consistent with the long-standing scenario that low competition and high rates of niche evolution promote diversification. The conceptual examples illustrate how using the synthetic Evolutionary Arena framework helps to identify and structure future directions for research on evolutionary radiations. In this way, the Evolutionary Arena framework promotes a more general understanding of variation in evolutionary rates by making quantitative results comparable between case studies, thereby allowing new syntheses of evolutionary and ecological processes to emerge.

A Phylogenomic Analysis of the Floral Transcriptomes of Sexually Deceptive and Rewarding European Orchids, Ophrys and Gymnadenia.

The orchids (Orchidaceae) constitute one of the largest and most diverse families of flowering plants. They have evolved a great variety of adaptations to achieve pollination by a diverse group of pollinators. Many orchids reward their pollinators, typically with nectar, but the family is also well-known for employing deceptive pollination strategies in which there is no reward for the pollinator, in the most extreme case by mimicking sexual signals of pollinators. In the European flora, two examples of these different pollination strategies are the sexually deceptive genus Ophrys and the rewarding genus Gymnadenia, which differ in their level of pollinator specialization; Ophrys is typically pollinated by pseudo-copulation of males of a single insect species, whilst Gymnadenia attracts a broad range of floral visitors. Here, we present and describe the annotated floral transcriptome of Ophrys iricolor, an Andrena-pollinated representative of the genus Ophrys that is widespread throughout the Aegean. Furthermore, we present additional floral transcriptomes of both sexually deceptive and rewarding orchids, specifically the deceptive Ophrys insectifera, Ophrys aymoninii, and an updated floral transcriptome of Ophrys sphegodes, as well as the floral transcriptomes of the rewarding orchids Gymnadenia conopsea, Gymnadenia densiflora, Gymnadenia odoratissima, and Gymnadenia rhellicani (syn. Nigritella rhellicani). Comparisons of these novel floral transcriptomes reveal few annotation differences between deceptive and rewarding orchids. Since together, these transcriptomes provide a representative sample of the genus-wide taxonomic diversity within Ophrys and Gymnadenia (Orchidoideae: Orchidinae), we employ a phylogenomic approach to address open questions of phylogenetic relationships within the genera. Specifically, this includes the controversial placement of O. insectifera within the Ophrys phylogeny and the placement of "Nigritella"-type morphologies within the phylogeny of Gymnadenia. Whereas in Gymnadenia, several conflicting topologies are supported by a similar number of gene trees, a majority of Ophrys gene topologies clearly supports a placement of O. insectifera as sister to a clade containing O. sphegodes.

Emergence of a floral colour polymorphism by pollinator-mediated overdominance.

Maintenance of polymorphism by overdominance (heterozygote advantage) is a fundamental concept in evolutionary biology. In most examples known in nature, overdominance is a result of homozygotes suffering from deleterious effects. Here we show that overdominance maintains a non-deleterious polymorphism with black, red and white floral morphs in the Alpine orchid Gymnadenia rhellicani. Phenotypic, metabolomic and transcriptomic analyses reveal that the morphs differ solely in cyanidin pigments, which are linked to differential expression of an anthocyanidin synthase (ANS) gene. This expression difference is caused by a premature stop codon in an ANS-regulating R2R3-MYB transcription factor, which is heterozygous in the red colour morph. Furthermore, field observations show that bee and fly pollinators have opposite colour preferences; this results in higher fitness (seed set) of the heterozygous morph without deleterious effects in either homozygous morph. Together, these findings demonstrate that genuine overdominance exists in nature.

The complete plastid genomes of Ophrys iricolor and O. sphegodes (Orchidaceae) and comparative analyses with other orchids.

Sexually deceptive orchids of the genus Ophrys may rapidly evolve by adaptation to pollinators. However, understanding of the genetic basis of potential changes and patterns of relationships is hampered by a lack of genomic information. We report the complete plastid genome sequences of Ophrys iricolor and O. sphegodes, representing the two most species-rich lineages of the genus Ophrys. Both plastomes are circular DNA molecules (146754 bp for O. sphegodes and 150177 bp for O. iricolor) with the typical quadripartite structure of plastid genomes and within the average size of photosynthetic orchids. 213 Simple Sequence Repeats (SSRs) (31.5% polymorphic between O. iricolor and O. sphegodes) were identified, with homopolymers and dipolymers as the most common repeat types. SSRs were mainly located in intergenic regions but SSRs located in coding regions were also found, mainly in ycf1 and rpoC2 genes. The Ophrys plastome is predicted to encode 107 distinct genes, 17 of which are completely duplicated in the Inverted Repeat regions. 83 and 87 putative RNA editing sites were detected in 25 plastid genes of the two Ophrys species, all occurring in the first or second codon position. Comparing the rate of nonsynonymous (dN) and synonymous (dS) substitutions, 24 genes (including rbcL and ycf1) display signature consistent with positive selection. When compared with other members of the orchid family, the Ophrys plastome has a complete set of 11 functional ndh plastid genes, with the exception of O. sphegodes that has a truncated ndhF gene. Comparative analysis showed a large co-linearity with other related Orchidinae. However, in contrast to O. iricolor and other Orchidinae, O. sphegodes has a shift of the junction between the Inverted Repeat and Small Single Copy regions associated with the loss of the partial duplicated gene ycf1 and the truncation of the ndhF gene. Data on relative genomic coverage and validation by PCR indicate the presence, with a different ratio, of the two plastome types (i.e. with and without ndhF deletion) in both Ophrys species, with a predominance of the deleted type in O. sphegodes. A search for this deleted plastid region in O. sphegodes nuclear genome shows that the deleted region is inserted in a retrotransposon nuclear sequence. The present study provides useful genomic tools for studying conservation and patterns of relationships of this rapidly radiating orchid genus.

Trans-generational inheritance of herbivory-induced phenotypic changes in Brassica rapa.

Biotic stress can induce plastic changes in fitness-relevant plant traits. Recently, it has been shown that such changes can be transmitted to subsequent generations. However, the occurrence and extent of transmission across different types of traits is still unexplored. Here, we assessed the emergence and transmission of herbivory-induced changes in Brassica rapa and their impact on interactions with insects. We analysed changes in morphology and reproductive traits as well as in flower and leaf volatile emission during two generations with leaf herbivory by Mamestra brassicae and Pieris brassicae and two subsequent generations without herbivory. Herbivory induced changes in all trait types, increasing attractiveness of the plants to the parasitoid wasp Cotesia glomerata and decreasing visitation by the pollinator Bombus terrestris, a potential trade-off. While changes in floral and leaf volatiles disappeared in the first generation after herbivory, some changes in morphology and reproductive traits were still measurable two generations after herbivory. However, neither parasitoids nor pollinators further discriminated between groups with different past treatments. Our results suggest that transmission of herbivore-induced changes occurs preferentially in resource-limited traits connected to plant growth and reproduction. The lack of alterations in plant-insect interactions was likely due to the transient nature of volatile changes.

Molecular mechanisms of adaptation and speciation: why do we need an integrative approach?

Understanding divergent adaptation and ecological speciation requires the synthesis of multiple approaches, including phenotypic characterization, genetics and genomics, realistic assessment of fitness and population genetic modelling. Current research in this field often approaches this problem from one of two directions: either a mechanistic approach-seeking to link phenotype, genotype and fitness, or a genomic approach-searching for signatures of divergence or selection across the genome. In most cases, these two approaches are not synthesized, and as a result, our understanding is incomplete. We argue that research in adaptation and evolutionary genetics needs to integrate these approaches for multiple reasons, including progress towards understanding the architecture and evolutionary history of adaptation and speciation loci, the ability to untangle linkage and pleiotropy, increased knowledge of mechanisms of genomic evolution and insights into parallel evolutionary events. Identifying the genetic underpinnings of adaptation and ecological speciation is not necessarily the end goal of research, but it is an integral part of understanding the evolutionary process. As a result, it is critical to utilize both genetic and genomic approaches. Challenges remain, particularly in nonmodel organisms and in our ability to synthesize results from multiple experimental systems. Nonetheless, advances in genetic and genomic techniques are increasingly available in a diverse array of systems, and the time is ripe to exploit the synthesis of these two approaches to increase our understanding of evolution.

Herbivore-Induced DNA Demethylation Changes Floral Signalling and Attractiveness to Pollinators in Brassica rapa.

Plants have to fine-tune their signals to optimise the trade-off between herbivore deterrence and pollinator attraction. An important mechanism in mediating plant-insect interactions is the regulation of gene expression via DNA methylation. However, the effect of herbivore-induced DNA methylation changes on pollinator-relevant plant signalling has not been systematically investigated. Here, we assessed the impact of foliar herbivory on DNA methylation and floral traits in the model crop plant Brassica rapa. Methylation-sensitive amplified fragment length polymorphism (MSAP) analysis showed that leaf damage by the caterpillar Pieris brassicae was associated with genome-wide methylation changes in both leaves and flowers of B. rapa as well as a downturn in flower number, morphology and scent. A comparison to plants with jasmonic acid-induced defence showed similar demethylation patterns in leaves, but both the floral methylome and phenotype differed significantly from P. brassicae infested plants. Standardised genome-wide demethylation with 5-azacytidine in five different B. rapa full-sib groups further resulted in a genotype-specific downturn of floral morphology and scent, which significantly reduced the attractiveness of the plants to the pollinator bee Bombus terrestris. These results suggest that DNA methylation plays an important role in adjusting plant signalling in response to changing insect communities.

Amino Acid Change in an Orchid Desaturase Enables Mimicry of the Pollinator's Sex Pheromone.

Mimicry illustrates the power of selection to produce phenotypic convergence in biology [1]. A striking example is the imitation of female insects by plants that are pollinated by sexual deception of males of the same insect species [2-4]. This involves mimicry of visual, tactile, and chemical signals of females [2-7], especially their sex pheromones [8-11]. The Mediterranean orchid Ophrys exaltata employs chemical mimicry of cuticular hydrocarbons, particularly the 7-alkenes, in an insect sex pheromone to attract and elicit mating behavior in its pollinators, males of the cellophane bee Colletes cunicularius [11-13]. A difference in alkene double-bond positions is responsible for reproductive isolation between O. exaltata and closely related species, such as O. sphegodes [13-16]. We show that these 7-alkenes are likely determined by the action of the stearoyl-acyl-carrier-protein desaturase (SAD) homolog SAD5. After gene duplication, changes in subcellular localization relative to the ancestral housekeeping desaturase may have allowed proto-SAD5's reaction products to undergo further biosynthesis to both 7- and 9-alkenes. Such ancestral coproduction of two alkene classes may have led to pollinator-mediated deleterious pleiotropy. Despite possible evolutionary intermediates with reduced activity, amino acid changes at the bottom of the substrate-binding cavity have conferred enzyme specificity for 7-alkene biosynthesis by preventing the binding of longer-chained fatty acid (FA) precursors by the enzyme. This change in desaturase function enabled the orchid to perfect its chemical mimicry of pollinator sex pheromones by escape from deleterious pleiotropy, supporting a role of pleiotropy in determining the possible trajectories of adaptive evolution.

Modeling the two-locus architecture of divergent pollinator adaptation: how variation in SAD paralogs affects fitness and evolutionary divergence in sexually deceptive orchids.

Divergent selection by pollinators can bring about strong reproductive isolation via changes at few genes of large effect. This has recently been demonstrated in sexually deceptive orchids, where studies (1) quantified the strength of reproductive isolation in the field; (2) identified genes that appear to be causal for reproductive isolation; and (3) demonstrated selection by analysis of natural variation in gene sequence and expression. In a group of closely related Ophrys orchids, specific floral scent components, namely n-alkenes, are the key floral traits that control specific pollinator attraction by chemical mimicry of insect sex pheromones. The genetic basis of species-specific differences in alkene production mainly lies in two biosynthetic genes encoding stearoyl-acyl carrier protein desaturases (SAD) that are associated with floral scent variation and reproductive isolation between closely related species, and evolve under pollinator-mediated selection. However, the implications of this genetic architecture of key floral traits on the evolutionary processes of pollinator adaptation and speciation in this plant group remain unclear. Here, we expand on these recent findings to model scenarios of adaptive evolutionary change at SAD2 and SAD5, their effects on plant fitness (i.e., offspring number), and the dynamics of speciation. Our model suggests that the two-locus architecture of reproductive isolation allows for rapid sympatric speciation by pollinator shift; however, the likelihood of such pollinator-mediated speciation is asymmetric between the two orchid species O. sphegodes and O. exaltata due to different fitness effects of their predominant SAD2 and SAD5 alleles. Our study not only provides insight into pollinator adaptation and speciation mechanisms of sexually deceptive orchids but also demonstrates the power of applying a modeling approach to the study of pollinator-driven ecological speciation.

Multiple shifts to different pollinators fuelled rapid diversification in sexually deceptive Ophrys orchids.

Episodes of rapid speciation provide unique insights into evolutionary processes underlying species radiations and patterns of biodiversity. Here we investigated the radiation of sexually deceptive bee orchids (Ophrys). Based on a time-calibrated phylogeny and by means of ancestral character reconstruction and divergence time estimation, we estimated the tempo and mode of this radiation within a state-dependent evolutionary framework. It appears that, in the Pleistocene, the evolution of Ophrys was marked by episodes of rapid diversification coinciding with shifts to different pollinator types: from wasps to Eucera bees to Andrena and other bees. An abrupt increase in net diversification rate was detected in three clades. Among these, two phylogenetically distant lineages switched from Eucera to Andrena and other bees in a parallel fashion and at about the same time in their evolutionary history. Lack of early radiation associated with the evolution of the key innovation of sexual deception suggests that Ophrys diversification was mainly driven by subsequent ecological opportunities provided by the exploitation of novel pollinator groups, encompassing many bee species slightly differing in their sex pheromone communication systems, and by spatiotemporal fluctuations in the pollinator mosaic.

Genic rather than genome-wide differences between sexually deceptive Ophrys orchids with different pollinators.

High pollinator specificity and the potential for simple genetic changes to affect pollinator attraction make sexually deceptive orchids an ideal system for the study of ecological speciation, in which change of flower odour is likely important. This study surveys reproductive barriers and differences in floral phenotypes in a group of four closely related, coflowering sympatric Ophrys species and uses a genotyping-by-sequencing (GBS) approach to obtain information on the proportion of the genome that is differentiated between species. Ophrys species were found to effectively lack postpollination barriers, but are strongly isolated by their different pollinators (floral isolation) and, to a smaller extent, by shifts in flowering time (temporal isolation). Although flower morphology and perhaps labellum coloration may contribute to floral isolation, reproductive barriers may largely be due to differences in flower odour chemistry. GBS revealed shared polymorphism throughout the Ophrys genome, with very little population structure between species. Genome scans for FST outliers identified few markers that are highly differentiated between species and repeatable in several populations. These genome scans also revealed highly differentiated polymorphisms in genes with putative involvement in floral odour production, including a previously identified candidate gene thought to be involved in the biosynthesis of pseudo-pheromones by the orchid flowers. Taken together, these data suggest that ecological speciation associated with different pollinators in sexually deceptive orchids has a genic rather than a genomic basis, placing these species at an early phase of genomic divergence within the 'speciation continuum'.

Transcriptome and proteome data reveal candidate genes for pollinator attraction in sexually deceptive orchids.

BACKGROUND: Sexually deceptive orchids of the genus Ophrys mimic the mating signals of their pollinator females to attract males as pollinators. This mode of pollination is highly specific and leads to strong reproductive isolation between species. This study aims to identify candidate genes responsible for pollinator attraction and reproductive isolation between three closely related species, O. exaltata, O. sphegodes and O. garganica. Floral traits such as odour, colour and morphology are necessary for successful pollinator attraction. In particular, different odour hydrocarbon profiles have been linked to differences in specific pollinator attraction among these species. Therefore, the identification of genes involved in these traits is important for understanding the molecular basis of pollinator attraction by sexually deceptive orchids. RESULTS: We have created floral reference transcriptomes and proteomes for these three Ophrys species using a combination of next-generation sequencing (454 and Solexa), Sanger sequencing, and shotgun proteomics (tandem mass spectrometry). In total, 121 917 unique transcripts and 3531 proteins were identified. This represents the first orchid proteome and transcriptome from the orchid subfamily Orchidoideae. Proteome data revealed proteins corresponding to 2644 transcripts and 887 proteins not observed in the transcriptome. Candidate genes for hydrocarbon and anthocyanin biosynthesis were represented by 156 and 61 unique transcripts in 20 and 7 genes classes, respectively. Moreover, transcription factors putatively involved in the regulation of flower odour, colour and morphology were annotated, including Myb, MADS and TCP factors. CONCLUSION: Our comprehensive data set generated by combining transcriptome and proteome technologies allowed identification of candidate genes for pollinator attraction and reproductive isolation among sexually deceptive orchids. This includes genes for hydrocarbon and anthocyanin biosynthesis and regulation, and the development of floral morphology. These data will serve as an invaluable resource for research in orchid floral biology, enabling studies into the molecular mechanisms of pollinator attraction and speciation.

The genetic basis of pollinator adaptation in a sexually deceptive orchid.

In plants, pollinator adaptation is considered to be a major driving force for floral diversification and speciation. However, the genetic basis of pollinator adaptation is poorly understood. The orchid genus Ophrys mimics its pollinators' mating signals and is pollinated by male insects during mating attempts. In many species of this genus, chemical mimicry of the pollinators' pheromones, especially of alkenes with different double-bond positions, plays a key role for specific pollinator attraction. Thus, different alkenes produced in different species are probably a consequence of pollinator adaptation. In this study, we identify genes that are likely involved in alkene biosynthesis, encoding stearoyl-acyl carrier protein (ACP) desaturases (SAD), in three closely related Ophrys species, O. garganica, O. sphegodes, and O. exaltata. Combining floral odor and gene expression analyses, two SAD homologs (SAD1/2) showed significant association with the production of (Z)-9- and (Z)-12-alkenes that were abundant in O. garganica and O. sphegodes, supporting previous biochemical data. In contrast, two other newly identified homologs (SAD5/6) were significantly associated with (Z)-7-alkenes that were highly abundant only in O. exaltata. Both molecular evolutionary analyses and pollinator preference tests suggest that the alkenes associated with SAD1/2 and SAD5/6 are under pollinator-mediated divergent selection among species. The expression patterns of these genes in F1 hybrids indicate that species-specific expression differences in SAD1/2 are likely due to cis-regulation, while changes in SAD5/6 are likely due to trans-regulation. Taken together, we report a genetic mechanism for pollinator-mediated divergent selection that drives adaptive changes in floral alkene biosynthesis involved in reproductive isolation among Ophrys species.

Floral isolation is the main reproductive barrier among closely related sexually deceptive orchids.

Floral isolation is an important component of pollinator-driven speciation. However, up to now, only a few studies have quantified its strength and relative contribution to total reproductive isolation. In this study, we quantified floral isolation among three closely related, sympatric orchid species of the genus Ophrys by directly tracking pollen flow. Ophrys orchids mimic their pollinators' mating signals, and are pollinated by male insects during mating attempts. This pollination system, called sexual deception, is usually highly specific. However, whether pollinator specialization also conveys floral isolation is currently under debate. In this study, we found strong floral isolation: among 46 tracked pollen transfers in two flowering seasons, all occurred within species. Accounting for observation error rate, we estimated a floral isolation index ≥0.98 among each pair of species. Hand pollination experiments suggested that postpollination barriers were effectively absent among our study species. Genetic analysis based on AFLP markers showed a clear species clustering and very few F(1) hybrids in natural populations, providing independent evidence that strong floral isolation prevents significant interspecies gene flow. Our results provide the first direct evidence that floral isolation acts as the main reproductive barrier among closely related plant species with specialized pollination.

Evidence for progenitor-derivative speciation in sexually deceptive orchids.

BACKGROUND AND AIMS: Sexually deceptive orchids of the genus Ophrys use mimicry of pollinator females to attract specific pollinators. Pollinator shifts may drive speciation in Ophrys, since novel pollinators may in principle act as isolating factors immediately. It is thus possible that evolution of novel species occurs rapidly and with a progenitor-derivative pattern. The aims of this study are to compare genetic structure and diversity among widespread and geographically restricted Ophrys taxa, to test whether genetic structure is associated with specific pollinators, and to investigate whether any widespread species may have acted as a progenitor for the evolution of more restricted taxa. METHODS: Genetic differentiation and diversity were investigated in O. leucadica and O. cinereophila, the two taxa of the Ophrys fusca sensu lato complex widespread in the Aegean, and three geographically restricted taxa from Rhodes, O. attaviria, O. parvula and O. persephonae, all differing in their specific pollinators. This was done using amplified fragment length polymorphism (AFLP) DNA fingerprinting, and sequencing of the low-copy nuclear gene LEAFY (LFY). KEY RESULTS: All taxa were found to be separate genetic entities, with O. leucadica forming two geographic groups from the west and east of the Aegean. Genetic structure was significantly shaped by pollinators and geography, and comparison of sequence and AFLP data revealed ancestral polymorphisms shared among several taxa. Among the sampled taxa, O. leucadica harbours the greatest genetic differentiation and geographic structure, and the highest genetic diversity. Part of the genome of O. parvula, endemic to Rhodes, may be derived from O. leucadica. CONCLUSIONS: Pollinators probably influence the genetic structure of the investigated Ophrys species. The genetic pattern identified is consistent with O. leucadica being the oldest of the sampled taxa, making O. leucadica a candidate progenitor species from which more restricted taxa such as O. parvula may have evolved.