- Research article
- Open Access
Unidirectional transitions in nectar gain and loss suggest food deception is a stable evolutionary strategy in Epidendrum (Orchidaceae): insights from anatomical and molecular evidence
BMC Plant Biologyvolume 18, Article number: 179 (2018)
Nectar gain and loss are important flower transitions observed in angiosperms, and are particularly common in orchids. To understand such transitions, the availability of detailed anatomical data and species-level phylogenies are crucial. We investigated the evolution of food deception in Epidendrum, one of the largest orchid genera, using genus phylogeny to map transitions between nectar gain and loss among different clades. Associations between anatomical and histochemical changes and nectar gain and loss were examined using fresh material available from 27 species. The evolution of nectar presence/absence in Epidendrum species was investigated in a phylogenetic framework of 47 species, using one nuclear and five plastid DNA regions available from GenBank and sequenced in this study.
The presence or absence of nectar was strongly associated with changes in the inner epidermal tissues of nectaries. Nectar-secreting species have unornamented epidermal tissue, in contrast to the unicellular trichomes found on the epidermis of food deceptive species. Bayesian tests confirmed that transitions occurred preferentially from nectar presence to nectar absence across the Epidendrum phylogeny. In addition, independent nectar loss events were found across the phylogeny, suggesting a lack of constraint for these transitions.
Ornamented nectaries may play an important role in the deceptive pollination strategy by secreting volatile organic compounds and providing tactile stimuli to pollinators. The recurrent and apparently irreversible pattern of nectar loss in Epidendrum suggests that food deception may constitute an alternative evolutionarily stable strategy, as observed in other orchid groups.
Floral rewards are one of the most conspicuous adaptations of angiosperms, playing a central role in the attraction and maintenance of pollinators . However, a large number of species deceive their pollinators by mimicking signals that are associated with food, sexual partners, or oviposition sites. The loss of flower rewards has several biological implications for plant populations. In general, pollinators avoid rewardless flowers by visiting distant plants, enhancing interpopulation genetic exchange, and consequently reducing the genetic structure of populations  and the levels of geitonogamy [3, 4]. Studies that added artificial nectar to rewardless flowers observed a significant increase in fruit set and self-pollination . According to Johnson et al. , selection toward rewardless species is favored when pollinators are common, and nectar production is selected when pollinators are scarce.
Orchids deceive their pollinators in several ways, but generalized food deception is the most common strategy . Despite the lack of nectar in many species, nectaries are widespread in the family, exhibiting extensive variation in morphology, anatomy and position. Nectaries can be associated with the lip callus, lip spur, column, or with a nectar spur formed by the fusion of lip, column, and sepal margins [7, 8]. In several genera of subtribe Laeliinae, there is a nectary-like structure below the column called a cuniculus. The cuniculus is a floral tube that originates at the fusion of the lip claw and column and penetrates the pedicel between the ovary and the perianth tube (Additional file 1: Figure S1). The presence of a cuniculus is often associated with Lepidoptera pollination, which is commonly found in Epidendrum L. (reviewed by ). Few anatomical studies have investigated this structure in depth (but see ), and whether nectar is present within the structure is unknown in many species. In Epidendrum, nectar has been detected only in species showing swollen cuniculi, in which the presence of copious nectar is evident . Thus, rewarding and rewardless Epidendrum species are distinguishable only through observing free nectar availability in floral nectaries [11, 12].
To understand the evolution of food deception, detailed anatomical data have been mapped into phylogenies, revealing the transitions between nectar gain and loss across different clades. A lack of floral nectar is a potential ancestral trait in the Orchidaceae due to the lack of nectar in the ancestral genera Apostasia and Neuwiedia, which offer pollen as a reward . Within the family, independent and reversible transitions have been observed in different clades. There are groups for which food deception is the ancestral state, as in the tribe Orchideae , and cases in which nectar reward evolved secondarily from food deception, as observed in the genera Disa  and Anacamptis . A gain or loss of nectar may result in different anatomical responses. Transitions may be associated with deep structural changes in nectaries, as observed in the subtribe Orchidinae , or may result in the pattern observed in Disa, in which little epidermal change has been observed between nectar-producing and nectarless species, suggesting that subcellular modifications may control nectar production .
In the present study, we investigated the evolution of nectar reward in the genus Epidendrum, which, with approximately 1500 species, is considered one of the largest plant genera of the neotropics . The presence/absence of nectar in Epidendrum species was inferred based only on visual inspection of nectaries, and detailed histochemical and/or anatomical analyses have been rarely reported in this genus [12, 20]. Due to the apparent lack of nectar, many authors have considered food deception to be widespread within Epidendrum (reviewed by ). Indirect evidence gathered from different data sources has been interpreted as collective evidence indicating the presence of a food deceptive strategy. The short visit times of pollinators that abandon the inflorescence after visiting a single flower , low fruit set , and the absence of nuclear genetic structure among populations  may indicate the existence of food deception in the analyzed species. According to , food deception in orchids tends to reduce the time spent by pollinators in the same inflorescence, decreasing the levels of geitonogamy and fruit set. Despite the putative lack of nectar, no Epidendrum species have been submitted to detailed anatomical analyses aimed at exploring the structure and nature of the epidermal tissue inside the nectary, or cuniculus.
Due to the large size of the genus, we focused our study on a smaller clade, the group Amphyglottium (sensu ), which has a previous phylogenetic hypothesis based on amplified fragment length polymorphism (AFLP) fingerprinting and plastid markers . New nuclear and plastid regions were sequenced to obtain a robust phylogenetic hypothesis for the group and to study the divergence time of main clades within the subgenus. Detailed anatomical observations on nectaries from all species were also included, using stereomicroscopy and scanning electron microscopy (SEM). The primary aim of the present study was to use a dated molecular phylogeny to study the evolution of nectar reward in the group Amphyglottium and selected Epidendrum species to identify nectar-producing and nectarless species. The following specific questions were addressed: (1) What is the timing and tempo of species diversification in the group Amphyglottium? (2) How many times have nectar and specific nectary types evolved in the group? (3) Are anatomical changes associated with the gain and loss of nectar production? Overall, this study supported the recurrent and unidirectional evolution of nectar loss in Epidendrum nectaries, which differ in anatomical features compared with nectaries from nectar-secreting species.
Structure and identification of the cuniculus glands
Epidendrum species have a cuniculus (a glandular chamber) positioned along the stylar canal between the perianth tube and the pedicellate ovary (Additional file 1: Figure S1). The inner epidermis of the cuniculus is always glandular, and its lumen varies among species and is wider in nectar-producing structures. Here we described the glands at secretory phase in flowers at anthesis. Four glandular morphologies (Figs. 1 and 2) were identified: (i) an almost uniform epidermis lacking trichomes, stomata, or cuticle disruption (designated “ordinary epidermis” Figs. 1a–g, 2a, b) and three sizes of glandular unicellular trichomes; (ii) short trichomes, composed of secretory papillae (Figs. 1h–p, 2c); (iii) medium trichomes (Figs. 1q–v, 2d); and (iv) long trichomes (Figs. 1w–z, 2e), which were at least twice as long as the medium ones.
Almost all nectar-producing flowers produced an exudate profuse enough to be tested with glucose strips. The exudates of E. densiflorum, E. mirmecophorum, E. nocturnum, and E. vesicatum tested positive for reducing sugars. The cuniculi of these species, as well as the remaining ones of the ordinary epidermis type, showed positive results for reducing sugars in the vacuoles of the epidermis and parenchyma (usually three cell layers), and were therefore identified as floral nectaries (Fig. 2f).
The glandular trichomes had thick cuticles and cell walls, regardless of their size. In addition, a conspicuous collenchyma tissue composing the one-to-three-layered subepidermal region was noteworthy (Fig. 2c–e). The absence of reducing sugars in all glandular trichomes and collenchyma tissue was also remarkable (Fig. 2g–i). Conversely, the occurrence of lipophilic droplets in the protoplast of the trichomes, observed under TEM (Fig. 2j, k), indicated that these glands may act as osmophores. In summary, all unornamented cuniculi showed signs of nectar secretion, whereas no signs of nectar were observed in ornamented cuniculi. This pattern was used to map the ancestral states in the phylogenetic tree and to infer the presence of the food deception strategy in the genus Epidendrum (see below).
Divergence times of Epidendrum clades and species
The Bayesian inferences for the first (primary calibration) and second (secondary calibration) analyses based on the combined data sets resulted in largely congruent topologies that proved stable in repeated runs. Tracers found stabilized likelihoods, smooth distributions, and high effective sample sizes for all parameters of all Bayesian runs. Bayesian inferences of the partitioned data sets (primary calibration = matK and rbcL; secondary calibration = matK, rbcL, trnL–trnF, trnT–trnL, rpl32–trnL, and ITS) produced the most resolved and well-supported trees (Additional file 2: Figure S2, Fig. 3).
The maximum credibility tree for the fossil-calibrated relaxed molecular clock analysis of the family Orchidaceae is shown in Additional file 1: Figure S2. Support values for the different clades and branches were high, with only a few BPP values below 0.90. Our age estimates indicate that the subfamilies Orchidoideae and Epidendroideae started to diversify in the Eocene, 51.6 Ma and 42.2 Ma, respectively (Additional file 1: Figure S2). The diversification of Higher Epidendroids appears to have started during the Oligocene, approx. 34.5 Ma. Our age estimates also indicate that the subtribe Laeliinae diversification started during the Early Miocene, approx. 20.1 Ma.
The maximum credibility tree for the calibrated relaxed molecular clock analysis of the Higher Epidendroids, focusing on the subtribe Laeliinae and the genus Epidendrum, is shown in Fig. 3. The support values were high for the clade including Encyclia, Cattleya, and Meiracyllium, with BPP > 0.90. All Epidendrum species clustered together, with high support values (PP > 0.90), supporting the monophyly of the genus. However, most clades within the genus showed low resolution, with BPP < 0.90, except group Amphyglottium. High support values were also observed within group Amphyglottium, and only the group formed by E. radicans and E. flexuosum showed support values below 0.90 BPP. The informal clades recovered by Pinheiro et al. (2009a) were also present in our analysis. Four main clades were observed, the Andean clade (E. calanthum, E. ibaguense, and E. macrocarpum), which is sister to the Tuberculata (E. cochlidium, E. secundum, and E. xanthinum) and Atlantic clades (E. cinnabarinum, E. denticulatum, E. flammeus, E. fulgens, and E. puniceoluteum), and the clade formed by E. orchidiflorum and E. purpureum, which is sister to the remaining species (Fig. 3). According to our dating, the genus Epidendrum started its diversification during the Late Miocene, approximately 12 Ma. Group Amphyglottium started to diversify in the early Pliocene, approximately 5.3 Ma.
Nectar presence/absence transitions
Ancestral state reconstruction suggested unidirectional transitions, with nectar presence as the most common ancestral state of most species (Fig. 4). The number of transitions from NA to NP (minimum 0, maximum 1, average 1.48) was lower than the converse situation (minimum 2, maximum 9, average 5.06). Indeed, MCMC analyses confirmed that the log-likelihood of the restricted model, allowing only transitions from NP to NA (q10), was significantly better than the unrestricted model or the restricted model allowing only transitions from NA to NP (q01) (Table 1). Nectar presence was the ancestral character state detected for most Epidendrum species, even in clades in which nectarless species were nested with nectar-producing species. The ancestral character state for the group Amphyglottium is uncertain, as nectar presence was recovered as the ancestral character state for this clade in only 35% of the trees. However, a different picture was observed within group Amphyglottium, in which most species do not produce nectar, and nectar absence was the primary ancestral condition in almost all clades.
We examined a total of 1100 species described by Hágsater et al. in Icones Orchidacearum (Additional file 3: Table S2), and a smooth inner epidermal surface was observed in the cuniculus (lack of epidermal unicellular trichomes) of most species. An ornamented cuniculus was found in 113 species (Additional file 3: Table S2), suggesting that approximately 10.3% of Epidendrum species do not produce nectar.
The evolution of rewardless species may have several biological implications, such as altering pollinator behavior, patterns of gene exchange, type of reproductive isolation, and rates of speciation. These implications are particularly well known for a few orchid groups, such as subfamily Orchidoideae [2, 3]. The availability of well-sampled phylogenies coupled with detailed anatomical analyses has provided crucial information for understanding how nectarless species evolved [17, 18]. The availability of such informative datasets has been a major challenge in understanding the evolution of rewardless species. In this study, datasets of DNA sequences and nectary anatomy were used to investigate the evolution of rewardless species in Epidendrum, mapping the transitions of nectar-producing and nectarless species over time. Epidendrum is a monophyletic group that started its diversification during the Late Miocene, likely from a nectar-secreting ancestor (Figs. 3 and 4). Transitions between nectar-secreting and nectarless species occurred independently several times along different time frames. Transitions were significantly unidirectional, with nectar presence the most probable ancestral state for most clades (Table 1). Despite several gains and losses detected along the phylogeny, almost all species from group Amphyglottium showed rewardless flowers. Comparing the tissue morphologies of closely related nectar-producing and nectarless species revealed that gain or loss of nectar production was followed by major histological changes (Figs. 1 and 2). Nectar-secreting species exhibited a uniform, unornamented inner epidermis of the cuniculus, with signs of reducing sugars, indicating nectar secretion. In contrast, glandular unicellular trichomes were found in the cuniculi of nectarless species. These trichomes were identified as osmophores and may play an important role in attracting pollinators by signaling a putative reward. The recurrent origin of glandular unicellular trichomes in deceptive species suggests convergent functional evolution, achieved through a widespread absence of morphological constraints associated with nectar losses, as observed in other orchid groups [17, 18].
The timing of diversification and nectary transitions in Epidendrum
The genus Epidendrum started its diversification during the Late Miocene (approx. 12 Ma, Fig. 3), coinciding with a period of accelerated Andean uplift . In fact, major shifts in diversification rates occurred in Cymbidieae during the Late Miocene as well, and Pleurothallidinae diversified along altitudinal gradients after major Andean uplifts . Givnish et al.  also highlighted the importance of epiphytic habit, pollination by Lepidoptera, and occurrence in tropical cordilleras as drivers of orchid diversification. Epidendrum species share most of the traits discussed above, and it is probable that multiple ecological drivers have contributed to the extensive diversification observed in this genus. Most Epidendrum species occur along different altitudes in the Andean region as epiphytes in humid forests and are pollinated by Lepidoptera .
Our anatomical data clearly associate nectar-producing species with the presence of an unornamented cuniculus, which is a widespread feature in Epidendrum. Of the 1100 species descriptions analyzed, only 10.3% showed ornamented cuniculi, which may indicate nectar absence. Thus, contrary to previous reports [9, 19], nectar-producing flowers may constitute the rule rather than the exception in this genus. However, this result needs a careful interpretation and additional confirmation by future studies in the genus, since our sample is very limited when considering the size of the genus Epidendrum.
How rewardless flowers contribute to speciation has been a matter of debate in the literature. Inda et al.  and Givnish et al.  discussed a potential association between increased diversification rates and the occurrence of pollination by deceit in orchids. These authors focused their discussion on different genera from subfamily Orchidoideae, for which the evolutionary outcomes of deceit strategies are well known . In Disa, transitions between deceptive and rewarding species did not change speciation rates . Species-level phylogenies have revealed that transitions between deceptive and rewarding pollination systems are not strongly constrained [15, 17]. Thus, rewardless species may evolve in response to certain ecological circumstances, such as pollen limitation, and arise from infrequent pollinator visits. On the other hand, the high levels of species diversity observed in Epidendrum may be a result of the high levels of genetic differentiation found in nectar-producing species, compared with rewardless species . Pollinators restrict their foraging range when visiting nectar-producing species, increasing the levels of genetic divergence among plant populations.
Species from group Amphyglottium are nested in a well-supported clade that started diverging during the Early Pliocene (5.3 Ma, Fig. 3). Nectar was consistently absent in all species except E. orchidiflorum (Fig. 4). Transitions were significantly asymmetric along the phylogeny, with character changes occurring preferentially from nectar presence to nectar absence (Table 1). Thus, nectar-secreting species did not evolve from ancestors lacking reward, suggesting that food deception may constitute an alternative, evolutionarily stable strategy, as observed in other orchid groups . These species are well known for their ability to thrive in harsh environments such as exposed rocky outcrops, sand dunes, and poor soils . Nectar production consumes a significant amount of energy, and nectar loss may represent a selective advantage in stressful habitats, such as those where species from group Amphyglottium are found. A similar pattern was observed by Antonelli et al.  in a group of Cattleya species (formerly Hoffmannseggella). All of these species are food deceptive , and their occurrence is concentrated in quartzite soils and rocks, which are extremely nutrient-poor. There are additional potential benefits of not producing nectar. For example, the increase in pollinator sharing commonly observed in food deceptive orchids may compensate for strong fluctuations in insect communities . Food deceptive species commonly exhibit weak population genetic structures because pollinators abandon the population when nectarless flowers are visited, transporting pollen over long distances and increasing interpopulation gene exchange [2, 3]. Moreover, hybridization is also commonly observed in food deceptive species from group Amphyglottium [4, 27]. Interspecific gene exchange may increase the adaptive responses of species to stressful habitat conditions by exchanging adaptive genes via introgression, as observed in other plant groups . Thus, food deception may constitute an alternative evolutionarily stable strategy in particular ecological conditions , enhancing intra- and interspecific gene exchange and decreasing the extinction risk of hybridizing species.
Anatomical changes in nectar-secreting and nectarless species
The structures of the cuniculus nectaries we observed (Figs. 1 and 2) differed from the descriptions available for the cuniculus of Brassavola flagellaris (Laeliinae; ) in the resemblance of this structure to the osmophores identified in the Epidendrum species lined with unicellular trichomes. The presence of reducing sugars in the exudate, however, undoubtedly identified these glands as nectaries . The nectary type and the absence of trichomes, stomata, or cuticle disruption suggest that nectar is released from nectar-producing cells by cuticular pores . Indeed, the cuticle covering the nectary cells may be completely permeable, allowing the reabsorption of nectar, as observed in Platanthera chlorantha  and Epidendrum vesicatum (Cardoso-Gustavson, pers. obs.).
In Orchidaceae, the osmophores are usually formed by a single layer of epidermal cells or may have secretory trichomes [31, 32]. Here, all trichome-bearing cuniculi, regardless of their size, were identified as osmophores. This identification was based on the presence of many lipophilic droplets in the protoplast, visible by TEM. Lipophilic droplets (also referred to as lipid droplets/bodies or lipophilic inclusions when referring to droplets in the cytoplasm or phenolic compounds in the vacuole, respectively; ) are a typical feature of osmophores, being identified in different orchid groups [31, 32, 34]. These droplets can eventually be observed in some nectaries (see ), although they may be restricted to specific stages of secretion . Nectaries comprising unicellular trichomes resembling the Epidendrum osmophores were described in Disa (Disinae; ). Indeed, glands were involved in the evolution of pollinating systems in orchid groups , constituting the fundamental functional element in the transition from nectarless to nectar-producing flowers .
The occurrence of lipophilic droplets in the protoplast can be directly associated with the low polarity of the volatile organic compounds (VOC)  composing the floral fragrance emitted by flowers. Indeed, among the species we analyzed, the only one with an available floral VOC profile was the nectarless moth-pollinated Epidendrum ciliare , in which the occurrence of terpenoids was remarkable. The occurrence of terpenoids has been associated with lipophilic droplets in the cytoplasm [40, 41], corroborating our identification of osmophores based on the massive amounts of these structures in glandular trichomes. The restricted chamber-like environment close to the trichomes may hold the highest concentration of VOCs, and their diffusion across the cuniculus to the atmosphere may act as a clue for pollinators. In this sense, the cuniculus may act as a VOC holder, in the same sense as the nectar holders (i.e., where nectar is confined before being released) found in Asclepiadaceae flowers . Thus, the cuniculus may provide tactile and olfactory stimuli to pollinators, enhancing their attraction and exploratory behavior, as suggested in other orchid groups [17, 43].
The recurrent and apparently irreversible pattern of nectar loss in Epidendrum suggests that food deception may constitute an alternative evolutionarily stable strategy, as observed in other orchid groups [17, 18]. The evolution of rewardless species may have several biological implications, such as altering pollinator behavior, patterns of gene exchange, type of reproductive isolation, and rates of speciation. The presence and absence of nectar in different Epidendrum species offer an opportunity to explore the ecological forces associated to the evolution of reward and rewardless species. Unfortunately, Epidendrum is one of the orchid genus poorly understood, despite its impressive diversity. Information about pollinators, reproductive biology and seed set in reward and rewardless species are needed, in order to understand the ecological conditions in which both strategies evolve. Population genetic and phylogeographic studies have focused mainly species from group Amphyglottium (reviewed by ) which is composed mainly by nectarless species, reducing our ability to understand the patterns of gene exchange in nectar offering species. For example, our hypothesis that ornamented nectaries may provide tactile stimuli to pollinators need to be investigated by monitoring pollinators behavior, observing the depth of penetration, period of residence and movement of the proboscis within the cuniculus. We are aware that, giving the size of the genus, our sample is limited, and some conclusions may be interpreted cautiously, such as the inference of nectar secreting species using the inner morphology of the cuniculus. Nonetheless, our study offer new insights and hypothesis that need to be tested in future studies, in order to clarify what evolutionary mechanisms are involved in the diversification of hyper diverse genera such as Epidendrum.
Study system and taxon sampling
The genus Epidendrum has been the focus of major taxonomic treatments in the last few years, focused mainly on describing new species. We focused our analysis on the group Amphyglottium due to its monophyly  and the availability of fresh flowers for anatomical analyses. Flowers were sampled from 27 species (31 accessions) maintained in the living orchid collection at the Instituto de Botanica (São Paulo, Brazil), from which 20 species (24 accessions) were also included in the phylogeny. The phylogenetic inference was based on sequence data obtained from GenBank and new data generated in the present study. Details are given in Additional file 4: Table S1.
The cuniculus glands from the flowers of 27 selected species (see Additional file 4: Table S1) were investigated at their secretory stage, when the flowers were at anthesis. The presence of nectar was assessed and established by glucose strips (whenever possible) and in situ by histochemical detection of reducing sugars.
Ovaries were excised from flowers, longitudinally sectioned to observe the epidermis of the cuniculus, and immediately fixed in 2.5% glutaraldehyde in phosphate buffer (pH 7.4, 0.1 M) for 24 h. The samples were washed in phosphate buffer, dehydrated in an ethanol series, and stored in 70% ethanol. The material was further dehydrated to 100% ethanol and rinsed in a hexamethyldisilazane (HMDS) series (33.3, 50.0, and 66.6% v/v in 100% ethanol) and then three times in 100% HMDS for 1 min each  to dry the material. Samples were mounted on stubs, coated with gold in a Leica ACE 200 sputtering system and viewed with a FEI Quanta 250 at 10 kV. Digital images were edited using Adobe Photoshop version 7.0 (Adobe Systems, Inc., San Jose, California, USA).
Fixed ovaries (see above) were embedded using standard methods for Leica historesin (Heraeus Kulzer, Germany) and serially sectioned at 5-μm thicknesses. The sections were stained with toluidine blue O  and mounted in water. Observations and digital images were acquired with an Olympus BX53 compound microscope equipped with an Olympus Q-Color 5 digital camera and Image Pro Express 6.3 software.
The presence of glucose in the exudate was assessed with a glucose strip test (Inlab Diagnostica; Alamar Tecno Cientifica, Buenos Aires, Argentina) whenever the nectar was profuse enough to be tested. Immediately after the excision of ovaries from fresh flowers, their internal liquid content was removed with microcapillary tubes and applied to glucose strips. For the in situ detection of reducing sugars, fresh cross-sectioned ovaries were treated with Fehling’s reagent, equal parts copper (II) sulfate 6.93% w:v, sodium potassium tartrate 34.6, and 12% sodium hydroxide w:w:v solutions heated to a pre-boiling temperature. After treatment, the sections were immediately observed under a light microscope.
Ovarian regions were isolated and immediately fixed at room temperature in 1.5% glutaraldehyde and 0.2 M cacodylate buffer, pH 7.25. After several washes with cacodylate buffer, the material was post-fixed in 1% OsO4 for 2 h, dehydrated to 100% ethanol, and them embedded in LR White embedding medium (Ted Pella Inc., Redding, CA, USA). The sections were mounted on copper grids and viewed with a Jeol JEM 2100 TEM at 80 kV.
For molecular analysis, leaf samples were sliced into small pieces and transferred to silica gel for drying. Total genomic DNA was extracted as described by Pinheiro et al. . We PCR-amplified and sequenced five plastid regions (matK, rbcL, trnL–trnF, trnT–trnL, and rpl32–trnL) and one nuclear region, the internal transcribed spacer (ITS) region (including the 5.8S gene), using the protocols and primers from Pinheiro et al.  for trnL–trnF, trnT–trnL, and rpl32–trnL, and van den Berg et al.  for matK, rbcL, and ITS. We believe more nuclear markers would be helpful but, due to the extensive chromosome variation found in Epidendrum species, it is very difficult to amplify and sequence low copy nuclear regions.
Amplification products were sequenced at Macrogen, Inc. (Seoul, South Korea) using the above PCR primers. Sequences were edited and aligned using the MUSCLE algorithm with default settings implemented in GENEIOUS version 5.4 software (Biomatters, Ltd., Auckland, New Zealand). We also performed a final editing of the sequence alignment by hand. For all plastid and nuclear regions, the sequence length variation among species was low. Hence, gaps were generally treated as missing data. GenBank sequence accession numbers and voucher details are summarized in Additional file 4: Table S1. We included outgroups and multiple intra-specific samples in the phylogenetic analyses, resulting in a total of 58 terminal accessions from 54 species. Two misidentifications in Pinheiro et al.  were corrected: the specimen identified as E. myrmecophorum was re-identified as E. orchidiflorum, and the specimen identified as E. incisum was re-identified as E. macrocarpum.
Phylogeny and time calibration of Epidendrum divergence
The joint posterior distribution of topologies and divergence times was estimated using Bayesian Markov chain Monte Carlo (MCMC) implemented in BEAST 1.8.1 software  in a two-step analysis with different calibration schemes. First, we estimated the divergence time of the main Orchidaceae clades with the following node calibrations: age for subtribe Goodyerinae, 15–20 Ma based on fossil records of orchid pollinia ; and divergence time estimates  for subfamily Epidendroideae (mean: 49 Ma), and the informal group “Higher Epidendroids” (mean: 39 Ma). This data set included the plastid DNA sequences ribulose-bisphosphate carboxylase (rbcL) and maturase K (matK) for 39 orchids from the aligned sequence matrix by Gustafsson et al. , plus ten Epidendrum species, resulting in 49 orchid taxa.
The second-step analysis was performed for Epidendrum using the age estimates of the first analysis as a secondary calibration: 1) age for subtribe Laellinae, normal prior distribution with mean 20.1 Ma and standard deviation 3.8 (95% CI, 13.85–26.35 Ma); and 2) age for clade “Meiracyllium, Cattleya and Epidendrum”: normal prior distribution with mean 12.62 Ma and standard deviation 4.2 (95% CI, 5.72–19.53 Ma). Fifteen species (19 accessions) from group Amphyglottium were analyzed, and 32 Epidendrum species were used as outgroups, plus 7 species belonging to four different Laeliinae genera (Arpophyllum, Cattleya, Encyclia, and Meiracyllium). This data set consisted of five plastid regions (matK, rbcL, trn–trnF, trnT–trnL, and rpl32–trnL) and one nuclear region, the ITS region (including the 5.8S gene). We also included the same rbcL and matK sequences used for the outgroups and the Epidendrum species used in the first-step calibration analysis. The General Time Reversible (GTR) model was selected for the nuclear marker ITS and the plastid markers rpl32–trnL, trnL–trnF, and trnT–trnL. The GTR + Γ + I model was selected for plastid markers matK and rbcL. All partitions were analyzed linked with an uncorrelated lognormal relaxed clock prior and a Yule tree prior. Three separated runs were performed in BEAST with 50,000,000 generations and a sampling frequency of every 1000 generations. Log files were analyzed with Tracer v. 1.6  to assess convergence and determine whether the effective sample sizes were larger than 200 for all parameters. The resulting trees were then combined with LogCombiner v. 1.8.1, with a burn-in of 10%. A maximum credibility tree was produced using TreeAnnotator v. 1.8.1 .
Preliminary analysis in our data revealed a covariance with a mean of 0.12 and 95% confidence intervals ranging from 0.0551–0.3055 (values close to zero indicate low covariance). The low covariance found in our data favor the use of BEAST instead of R8S . As reported in the literature [50, 54], using the penalized likelihood with uncorrelated data would fail to obtain robust age estimates. Thus, it would be difficult to compare age estimates obtained by R8S and Beast, precluding categorical assertion son which of these methods yields the most correct divergence estimate .
The congruence between nuclear and plastid data sets was assessed by comparing the topologies and posterior probabilities of the strict consensus trees of two main data partitions, the nuclear ITS, and all plastid regions together, according to the method of Xiang et al. . The test was performed by visually comparing the support and resolution of each of the clades in the separate analyses that had a higher posterior probability (PP) > 90 . Because most incongruences were restricted to clades showing low support values (Bayesian posterior probability [BPP] < 0.90 – data not shown), only the Bayesian inference based on the combined nuclear and plastid data sets was used for phylogenetic inference, divergence time estimations, and analysis of character evolution.
Evolution of nectar gain and loss
Nectar presence and absence was directly inferred from anatomical observations and histochemical tests. In addition, we interpreted the presence of an ornamented cuniculus as a sign of nectar absence. Conversely, the presence of an unornamented cuniculus was indicative of nectar presence, according to our results (Additional file 4: Table S1). Ancestral character state reconstruction analysis was used to map the transitions between nectar absence (NA) and nectar presence (NP) in Epidendrum species. The NA and NP states were coded as binary data (0 and 1, respectively) and optimized onto the combined (total evidence = nuclear + plastid datasets) Bayesian inference tree under a Maximum Likelihood criterion using the Mesquite v. 3.31 package . The number and directionality of transitions between NA and NP states were quantified using the Summarize State Changes Over Trees function in Mesquite. The number of state transitions was summarized over 10,000 chronograms.
Differences in the rates of transitions between NA and NP were tested following Pinheiro et al. , using 10,000 chronograms. Three models were compared using likelihood ratio tests: (1) the unrestricted model, in which the probability of the two types of transitions, from NA to NP (q01) and the converse (q10), were calculated; (2) a restricted model in which only NA to NP transitions were permitted (i.e., q10 = 0); and (3) an alternative, restricted model in which only NP to NA transitions were permitted (i.e., q01 = 0). The likelihood ratio test was used to compare the two likelihoods derived from the unrestricted and each of the restricted models. In the restricted models, ancestral states were fixed to NA or NP. MCMC analyses were performed using the BayesMultistates program  implemented in BayesTraits v. 3.0.1 (http://www.evolution.rdg.ac.uk/BayesTraits.html). The MCMC analyses were run for 5,050,000 generations, with a uniform prior distribution and a burn-in of 50,000 generations.
The presence and absence of nectar was directly inferred for the 27 species investigated in detail based on anatomical observations. We also used an indirect approach to infer nectar presence/absence in species for which fresh material was not available, using associations between epidermal characters and the presence/absence of nectar. We found a strong association between an ornamented cuniculus (presence of epidermal unicellular trichomes) and lack of nectar. In contrast, an unornamented cuniculus (uniform epidermis lacking trichomes) was always associated with nectar presence (see description in results). Thus, to estimate the proportion of nectarless Epidendrum species, we used the species described by Hágsater et al. in Icones Orchidacearum (Additional file 3: Table S2) as a proxy for the whole genus. We restricted our analysis to the species described in Icones Orchidacearum due to the detailed and standardized nature of its flower descriptions, providing information for almost the entire genus (approximately 1100 species). Detailed distribution maps for all species analyzed here are also provided by Icones Orchidacearum.
Amplified fragment length polymorphism
Bayesian posterior probability
General Time Reversible model
Internal transcribed spacer
Million of years
Markov-chain Monte Carlo
Scanning electron microscopy
Transmission electron microscopy
Volatile organic compounds
Van der Pijl L, Dodson CH: Orchid flowers--their Pollination and Evolution. 1966, University of Miami Press, Coral Gables, Florida, USA.
Cozzolino S, Widmer A. Orchid diversity: an evolutionary consequence of deception? Trends Ecol Evol. 2005;20:487–94.
Scopece G, Cozzolino S, Johnson SD, Schiestl FP. Pollination efficiency and the evolution of specialized deceptive pollination systems. Am Nat. 2010;175:98–105.
Pinheiro F, Gouveia TMZM, Cozzolino S, et al. Strong but permeable barriers to gene exchange between sister species of Epidendrum. Am J Bot. 2016;103:1472–82.
Jersáková J, Johnson SD, Kindlmann P. Mechanisms and evolution of deceptive pollination in orchids. Biol Rev. 2006;81:219–35.
Johnson SD, Peter CI, Agren J. The effects of nectar addition on pollen removal and geitonogamy in the non-rewarding orchid Anacamptis morio. Proc R Soc Lond B Biol Sci. 2004;271:803–9.
Stpiczyńska M. Nectar resorption in the spur of Platanthera chlorantha (Custer) Rchb. – structural and microautoradiographic studies. Plant Syst Evol. 2003;238:119–26.
Stpiczyńska M, Davies KL, Gregg A. Comparative account of nectary structure in Hexisea imbricata (Lindl.) Rchb. F.(Orchidaceae). Ann Bot. 2005;95:749–56.
Pinheiro F, Cozzolino S. Epidendrum (Orchidaceae) as a model system for ecological and evolutionary studies in the Neotropics. Taxon. 2013;62:77–88.
Stpiczyńska M, Davies KL, Kaminska M. Structure of the cuniculus nectary in Brassavola flagellaris barb. Rodr. (Laeliinae Benth., Orchidaceae). Acta Agrobot. 2010;63:3–10.
Almeida AM, Figueiredo RA. Ants visit nectaries of Epidendrum denticulatum (Orchidaceae) in a Brazilian rainforest: effects on herbivory and pollination. Braz J Biol. 2003;63:551–8.
Pansarin ER, Amaral MCE. Reproductive biology and pollination mechanisms of Epidendrum secundum (Orchidaceae). Floral variation: a consequence of natural hybridization? Plant Biol. 2008;10:211–9.
Kocyan A, Endress PK. Floral structure and development of Apostasia and Neuwiedia (Apostasioideae) and their relationships to other Orchidacea. Int J Plant Sci. 2001;162:847–67.
Inda LA, Pimentel M, Chase MW. Phylogenetics of tribe Orchideae (Orchidaceae: Orchidoideae) based on combined DNA matrices: inferences regarding timing of diversification and evolution of pollination syndromes. Ann Bot. 2012;110:71–90.
Johnson SD, Hobbhahn N, Bytebier B. Ancestral deceit and labile evolution of nectar production in the African orchid genus Disa. Biol Lett. 2013;9:20130500.
Cozzolino S, Caputo P, Aceto S, Widmer A, Dafni A. Speciation processes in eastern Mediterranean Orchis sl species: molecular evidence and the role of pollination biology. Isr J Plant Sci. 2001;49:91–103.
Bell AK, Roberts DL, Hawkins JA, Rudall PJ, Box MS, Bateman RM. Comparative micromorphology of nectariferous and nectarless labellar spurs in selected clades of subtribe Orchidinae (Orchidaceae). Bot J Linn Soc. 2009;160:369–87.
Hobbhahn N, Johnson SD, Bytebier B, Yeung EC, Harder LD. The evolution of floral nectaries in Disa (Orchidaceae: Disinae): recapitulation or diversifying innovation? Ann Bot. 2013;112:1303–19.
Hágsater E, Soto-Arenas MA. Epidendrum L. In: Pridgeon AM, Cribb PJ, Chase MW, Rasmussen FN, editors. Genera Orchidacearum, vol. 4. Oxford: Oxford University Press; 2005. p. 236–51.
Pansarin ER, Pansarin LM. Reproductive biology of Epidendrum tridactylum (Orchidaceae: Epidendroideae): a reward-producing species and its deceptive flowers. Plant Syst Evol. 2014;300:321–8.
Pinheiro F, Koehler S, Corrêa AM, Salatino MLF, Salatino A, Barros F. Phylogenetic relationships and infrageneric classification of Epidendrum subgenus Amphiglottium (Laeliinae, Orchidaceae). Plant Syst Evol. 2009a;283:165–77.
Gregory-Wodzicki KM. Uplift history of the central and northern Andes: a review. Geol Soc Am Bull. 2000;112:1091–105.
Pérez-Escobar OA, Chomicki G, Condamine FL, et al. Recent origin and rapid speciation of Neotropical orchids in the world's richest plant biodiversity hotspot. New Phytol. 2017;215:891–905.
Givnish TJ, Spalink D, Ames M, et al. Orchid phylogenomics and multiple drivers of their extraordinary diversification. Proc R Soc Lond B Biol Sci. 2015;282:20151553.
Antonelli A, Verola CF, Parisod C, Gustafsson ALS. Climate cooling promoted the expansion and radiation of a threatened group of south American orchids (Epidendroideae: Laeliinae). Biol J Linn Soc Lond. 2010;100:597–607.
Verola CF. Estudos biossistemáticos em espécies de Hoffmannseggella H.G. Jones (Orchidaceae: Laeliinae) ocorrentes nos complexos rupestres de altitude. 2008, (Doctoral dissertation, [Campinas]).
Pinheiro F, Barros F, Palma-Silva C, et al. Hybridization and introgression across different ploidy levels in the Neotropical orchids Epidendrum fulgens and E. puniceoluteum (Orchidaceae). Mol Ecol. 2010;19:3981–94.
Rieseberg LH, Kim SC, Randell RA, et al. Hybridization and the colonization of novel habitats by annual sunflowers. Genetica. 2007;129:149–65.
Lemos RCC, Silva DC, Melo-de-Pinna GFA. A structural review of foliar glands in Passiflora L.(Passifloraceae). PLoS One. 2017;12:e0187905.
Paiva EAS. How does the nectar of stomata-free nectaries cross the cuticle? Acta Bot Bras. 2017;31:525–30.
Pansarin LM, Castro MM, Sazima M. Osmophore and elaiophores of Grobya amherstiae (Catasetinae, Orchidaceae) and their relation to pollination. Bot J Linn Soc. 2009;159:408–15.
Melo MC, Borba EL, Paiva EAS. Morphological and histological characterization of the osmophores and nectaries of four species of Acianthera (Orchidaceae: Pleurothallidinae). Plant Syst Evol. 2010;286:141–51.
Curry KJ, McDowell LM, Judd WS, Stern WL. Osmophores, floral features, and systematics of Stanhopea (Orchidaceae). Am J Bot. 1991;78:610–23.
Pridgeon AM, Stern WL. Ultrastructure of osmophores in Restrepia (Orchidaceae). Am J Bot. 1983;70:1233–43.
Paiva EA. Anatomy, ultrastructure, and secretory activity of the floral nectaries in Swietenia macrophylla (Meliaceae). Am J Bot. 2012;99:1910–7.
Pansarin ER, Maciel AA. Evolution of pollination systems involving edible trichomes in orchids. AoB PLANTS 2017;9:plx033.
Neubig KM, Carlsward BS, Whitten MW, Williams NH. Nectary structure and nectar in Sobralia and Elleanthus (sobralieae: orchidaceae). Lankesteriana. 2015;15:113–27.
Piechulla B, Pott MB. Plant scents-mediators of inter-and intraorganismic communication. Planta. 2003;217:687–9.
Moya S, Ackerman JD. Variation in the floral fragrance of Epidendrum ciliare (Orchidaceae). Nord J Bot. 1993;13:41–7.
He X, Sun Y, Zhu RL. The oil bodies of liverworts: unique and important organelles in land plants. CRC Crit Rev Plant Sci. 2013;32:293–302.
Widhalm JR, Jaini R, Morgan JA, Dudareva N. Rethinking how volatiles are released from plant cells. Trends Plant Sci. 2015;20:545–50.
Monteiro MM, Demarco D. Corona development and floral nectaries of Asclepiadeae (Asclepiadoideae, Apocynaceae). Acta Bot Bras. 2017;31:420–32.
Davies KL, Stpiczyńska M. Structure and distribution of floral trichomes in Lycaste and Sudamerlycaste (Orchidaceae: Maxillariinae s.L.). Bot J Linn Soc. 2010;164:409–21.
Jeger RN, Lichtenfeld Y, Peretz H, Shany B, Vago R, Baranes D. Visualization of the ultrastructural interface of cells with the outer and inner surface of coral skeletons. J Electron Microsc. 2009;58:47–53.
Sakai WS. Simple method for differential staining of paraffin embedded plant material using toluidine blue O. Stain Technol. 1973;48:247–9.
Pinheiro F, Palma-Silva C, Barros F, et al. Chloroplast microsatellite markers for the Neotropical orchid genus Epidendrum, and cross amplification in other Laeliinae species (Orchidaceae). Conserv Genet Resour. 2009b;1:505–11.
Van den Berg C, Goldman DH, Freudenstein JV, Pridgeon AM, Cameron KM, Chase MW. An overview of the phylogenetic relationships within Epidendroideae inferred from multiple DNA regions and recircumscription of Epidendreae and Arethuseae (Orchidaceae). Am J Bot. 2005;92:613–24.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Ramirez SR, Gravendeel B, Singer RB, Marshall CR, Pierce NE. Dating the origin of the Orchidaceae from a fossil orchid with its pollinator. Nature. 2007;448:1042–5.
Gustafsson ALS, Verola CF, Antonelli A. Reassessing the temporal evolution of orchids with new fossils and a Bayesian relaxed clock, with implications for the diversification of the rare south American genus Hoffmannseggella (Orchidaceae: Epidendroideae). BMC Evol Biol. 2010;10:177.
Rambaut A, Drummond A. Tracer v1.6.0. 2013, [http://beast.bio.ed.ac.uk/Tracer].
Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.
Sanderson MJ. r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics. 2003;19:301–2.
Lohmann LG, Bell CD, Calió MF, Winkworth RC. Pattern and timing of biogeographical history in the Neotropical tribe Bignonieae (Bignoniaceae). Bot J Linn Soc. 2012;171:154–70.
Magallón S. Using fossils to break long branches in molecular dating: a comparison of relaxed clocks applied to the origin of angiosperms. Syst Biol. 2010;59:384–99.
Xiang XG, Schuiteman A, Li DZ, et al. Molecular systematics of Dendrobium (Orchidaceae, Dendrobieae) from mainland Asia based on plastid and nuclear sequences. Mol Phylogenet Evol. 2013;69:950–60.
Wiens JJ. Combining data sets with different phylogenetic histories. Syst Biol. 1998;47:568–81.
Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. 2011, [http://mesquiteproject.org].
Pinheiro F, Cafasso D, Cozzolino S, Scopece G. Transitions between self-compatibility and self-incompatibility and the evolution of reproductive isolation in the large and diverse tropical genus Dendrobium (Orchidaceae). Ann Bot. 2015;116:457–67.
Pagel M, Meade A, Barker D. Bayesian estimation of ancestral character states on phylogenies. Syst Biol. 2004;53:673–84.
Raftery AE. Hypothesis testing and model selection In: Gilks WR, Richardson S, Spiegelhalter DJ, editors. Markov Chain Monte Carlo in Practice. London: Chapman & Hall; 1996. p. 163–87.
We thank Alexandre Antonelli, Maria Fernanda Calió, Thiago André and Thaís Vasconcelos for discussions and comments regarding the early ideas for this manuscript.
This research was supported by a CNPq Universal grant (No. 440367/2014–0) and fellowships from PNPD/CAPES awarded to FP (300927/2016–9) and PCG (0469-13-0), and CNPq awarded to PCG (0588/12–1). FP and CPS also acknowledges the productivity grant from CNPq (300927/2016–9, 300819/2016–1, respectively).
Availability of data and materials
All of the phylogenetic sequence data generated in this study was deposited in GenBank (National Center for Biotechnology Information) under accession numbers: ITS, MH218741 - MH218763; matK, MH218764 - MH218779; trnT-L, MH218803 - MH218807; rpl32–trnL: MH218793 - MH218798; rbcL, MH218780 - MH218792; trnL-F, MH218799 - MH218802; (Additional file 4: Table S1). All additional materials supporting the results of the article are included as additional files.
Ethics approval and consent to participate
Plant materials used in this study were collected in field with necessary permissions from local authorities of Forest Department. All voucher specimens have been deposited in publicly available herbaria (Additional file 4: Table S1).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Details of Epidendrum flowers, showing the cuniculi. (A) E. coronatum flower. The red rectangle indicates the pedicel region dissected by longitudinal sections in the remaining pictures (B–D). (B) Detail of the cuniculus of E. cristatum. The dotted line indicates the area from which most samples were taken for anatomical analyses. (C) Unornamented cuniculus of E. orchidiflorum. (D) Ornamented cuniculus of E. fulgens. Scale bars = 1.0 cm. (PDF 6036 kb)
Figure S2. Time-calibrated tree of the Orchidaceae, focusing on the subfamilies Orchidoideae and Epidendroideae, based on matK and rbcL plastid regions, estimated using BEAST software. Thin branches indicate posterior probabilities below 0.9. Circles indicate age-constrained nodes. The yellow circle indicates the calibration point for subtribe Goodyerinae (Pachyplectron–Dossinia), the blue circle indicates the calibration point for Epidendroideae, and the green circle indicates the calibration point for the Higher Epidendroids. Numbers at nodes represent median ages in millions of years (Ma). (PDF 665 kb)
Table S2. Species described in Icones Orchidacearum (1993–2016) showing ornamented cuniculus, including their informal group and specific reference. (DOCX 20 kb)
Table S1. Epidendrum species and allied genera analyzed in this study, including GenBank access numbers by molecular marker, including the type of analyses performed (PC = primary calibration, SC = secondary calibration, AN = anatomy), and presence or absence of nectar with the corresponding source of information. (DOCX 30 kb)