- Research article
- Open Access
Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end
© Gamisch et al. 2015
- Received: 6 July 2015
- Accepted: 28 August 2015
- Published: 16 September 2015
The transition from outcrossing to selfing has long been portrayed as an ‘evolutionary dead end’ because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a ‘rostellum’) co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.
The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4–3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.
The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the ‘rostellum’) is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.
- Internal Transcribe Spacer
- Markov Chain Monte Carlo
- Mating Type
- Extinction Rate
- Bayesian Inference Analysis
Flowering plants (angiosperms) display a bewildering diversity and, in many cases, complexity of mating systems compared to most animal groups [1, 2]. The overwhelming majority of angiosperms are characterized by cross-pollination of their flowers, whereby mechanisms preventing self-fertilization range from the separation of sex functions within flowers (herkogamy, dichogamy) and plants (monoecy), to separate sexes (dioecy), and self-incompatibility (SI) systems [2, 3]. Still, about 20–25 % of all flowering plant species are predominantly selfing , which is commonly thought to reflect two main advantages, which are not necessarily mutually exclusive. First, selfing genes should enjoy a two-fold transmission advantage over genes for outcrossing [5–7] and, second, the loss or paucity of pollinators or mates should selectively favour selfing individuals/populations because it provides reproductive assurance [8–10]. On the other hand, inbreeding depression has long been recognized to counteract the evolution of selfing and explaining the maintenance of outcrossing [11–13]. In fact, selfing lineages or taxa are typically assumed to be an evolutionary ‘blind alley’, or ‘dead end’ ([14–17]; selfing as an evolutionary dead end or SEDE hypothesis), presumably because of their limited capacity to adapt to changing environment, or their susceptibility to the accumulation of deleterious mutations [2, 7, 14]. Taking a macroevolutionary perspective, however, our body of knowledge is still limited regarding the historical evolutionary dynamics of mating system transitions among closely related taxa differing in outcrossing vs. selfing strategies and over extended time scales.
Several underlying assumptions of the SEDE hypothesis, as originally proposed by Stebbins [15, 16], have turned out to be more complex than initially thought. Nonetheless, there is wide consensus about two major predictions: (1) evolutionary change is unidirectional and irreversible from outcrossing to selfing; and (2) selfing lineages suffer from an increased risk of extinction owing to a reduced potential for adaptation [7, 14, 17]. Consistent with both connotations, several molecular phylogenies have inferred outcrossing systems to be ancestral, with selfing species commonly located at terminal branches [18–23]. However, the historical dynamics of outcrossing-to-selfing transitions are often difficult to reconstruct and a few studies have even pointed out the possibility of reversals [21, 24–29]. Moreover, potential departures from theory are generally suspected to result from methodological difficulties associated with, for example, phylogenetic estimation per se, incomplete or missing taxonomic accounts, ambiguous character state coding, and/or subjectively chosen character transformation weights [17, 24, 25, 30, 31]. Other challenges are posed by the presence of intraspecific mating type polymorphisms (e.g., ), and possible confounding effects of diversification and character state evolution [28, 32, 33]. Perhaps most strikingly, the second claim of the SEDE hypothesis, concerning the assumed higher extinction [viz. negative net diversification (speciation minus extinction) rate] associated with selfing, has largely been ignored in the literature , possibly because of the difficulties of obtaining accurate and unbiased estimates of extinction rates from molecular phylogenies [33–35].
However, recent developments of model-based methods of phylogenetic inference, molecular dating, and ancestral character state reconstructions [36, 37] should overcome at least some of the past methodological obstacles encountered in phylogenetic studies of mating type evolution. More reliable tests of the SEDE hypothesis are also expected from a recently developed, yet sparingly employed tree-based model of character evolution (Binary State Speciation and Extinction [BiSSE]; [33, 38]) that allows the joint estimation of rates of binary character change as well as state-dependent rates of extinction and speciation [39–41]. By capitalizing on these recent advances, the present study investigates mating type shifts (outcrossed vs. self-fertilized) during the diversification of a well-defined clade of epidendroid orchids of the genus Bulbophyllum Thouars from Madagascar and adjacent islands. No study to date has subjected members of Orchidaceae to such kind of macroevolutionary analysis, and there are several reasons detailed below why this group of tropical orchids provides a particular interesting system for the study of selfing evolution.
Observing a relatively high incidence of auto-pollination in Bulbophyllum is rather unexpected because this genus has long been considered to be highly pollinator-dependent (see above). Even though estimates of auto-pollination in Orchidaceae as a whole (ca. 25,000 spp.) are relatively high with about 31 % of species [51, 52], auto-pollination is generally considered to be rare in species from the tropics [45, 53], excepting small islands (e.g., La Réunion; ). When compared to other earlier phylogenetic studies of selfing evolution [7, 20, 21, 26–28, 55, 56], our study system is notable in at least two respects: (1) selfing is not achieved through the breakdown of a floral multi-organ (e.g., heterostyly) or genetic (e.g., SI) mechanism associated with outcrossing but rather due to a structural, discrete-state character transition affecting floral micro-morphology (i.e., rostellum presence/absence); and (2) selfing morphs (‘selfers’) apparently only coexist with their pollinator-dependent conspecifics (‘outcrossers’), i.e., selfing is never fixed within species. Because conspecific individuals differing in morph-type are virtually identical in overall (floral, vegetative) phenotype [47, 48], selfing may have evolved only recently and multiple times within clade C. This hypothesis, however, has yet to be tested since phylogenetic relationships within clade C still remain obscure due to incomplete taxon and gene sampling in an earlier phylogenetic analysis .
The overall aim of this study was to examine the polarity, frequencies, time scales, and consequences of mating type transitions during the diversification of clade C in order to test major predictions raised by the SEDE hypothesis, i.e., irreversibility of transitions to selfing and high extinction rate of selfing species. Our first aim was to infer a robust, fossil-calibrated phylogeny of clade C using DNA sequence data (five plastid and three nuclear loci) to explore the timing and tempo of species diversification. Second, we reconstructed ancestral character states onto the phylogenetic tree to evaluate when and how often transitions occurred between outcrossing and ‘selfing’ (viz. dimorphic) mating types, and whether uni- or bi-directionally. Finally, we used the state-dependent BiSSE model of character evolution to recover the processes proposed by the SEDE hypothesis, and to test in particular whether the evolution of selfing in these tropical orchids was associated with higher levels of extinction. Overall, this study is the first to demonstrate within a macroevolutionary framework that, contrary to prevailing views, selfing not always needs to be an endpoint of mating type evolution in flowering plants, suggesting instead more frequent origins of this character than previously thought.
Study system and mating system categorization
The most recent taxonomic treatment of sects. Calamaria, Humblotiorchis, and Bifalcula recognizes 32 species (; Fischer et al., unpubl. data). Fifteen of those species were included in the study of Madagascan Bulbophyllum conducted by , using a combined sequence matrix of one nuclear ribosomal (ITS) and four plastid gene regions, and were identified as a monophyletic group (‘clade C’; Bayesian posterior probability [PP] of 0.93), albeit with poorly resolved interspecies relationships.
Recent micro-morphological and experimental studies [47, 48] revealed that mating type in clade C is directly related to column (gynostemium) type, most often involving the presence or absence of the rostellum. In detail, of the 29 species analysed by  (1–25 individuals per species; mean ± SD: 7.17 ± 6.11; 208 individuals in total, 21 species were found to possess only pollinator-dependent individuals (‘outcrossers’) with the conventional (erect) rostellum (Type I); by contrast, seven species proved to be dimorphic for mating type, comprising both Type I individuals and others (Type II) capable of vector-less, autonomous fruit set (‘selfers’) due to the lack of the rostellum (Fig. 1). Moreover, a single species (B. bicoloratum) was found to have three mating types by additionally possessing selfing individuals with a receptive (sub-erect) rostellum (Type III; see also ).
Coding of mating types for 30 Bulbophyllum clade C species included in this study, following the mating type classification of . Species fixed for outcrossing were coded as state ‘0’, and dimorphic (outcrossing/selfing) species as state ‘1’ (‘any-instance’ coding; see text for details). For two species (B. spp. nov. ‘B’ and ‘D’) mating type was coded as missing (‘NA’)
B. capuronii Bosser
B. complanatum H.Perrier
B. implexum Jum. & H.Perrier.
B. minutum Thouars
B. sp. nov. ‘A’
B. sp. nov. ‘B’
B. bicoloratum Schltr.
B. cirrhoglossum H.Perrier
B. elliottii Rolfe
B. erectum Thouars
B. hildebrandtii Rchb.f.
B. histrionicum G.A.Fischer & P.J.Cribb
B. incurvum Thouars
B. lecouflei Bosser
B. luteobracteatum Jum. & H.Perrier
B. malawiense B.Morris
B. obtusatum (Jum. & H.Perrier) Schltr.
B. occultum Thouars
B. pervillei Rolfe
B. pusillum (H.Perrier) G.A.Fischer & P.J.Cribb
B. quadrifarium Rolfe
B. rubrum Jum. & H.Perrier
B. ruginosum H.Perrier
B. senghasii G.A.Fischer & A.Sieder
B. trifarium Rolfe
B. sp. nov. ‘C’
B. sp. nov. ‘D’
B. sp. nov. ‘E1’
B. sp. nov. ‘E2’
B. humblotii Rolfe
The BiSSE method used in the present study (detailed below) requires that the states that are compared must be coded as a binary character [21, 39, 40]. If species exhibit multiple (‘polymorphic’) character states, the most common state is usually taken as the only one present (e.g., [21, 40, 56]). However, such a frequency-dependent coding strategy is prone to sampling bias if sample sizes of character statement assessment are low , as in the present study. To accommodate this latter problem, we used the ‘any instance’ method (i.e., fixation of the ancestral state = 0; and polymorphism or fixation of the derived state = 1; [59, 60]) and coded the mating types of clade C species (; Table 1) as being either outcrossing (0) or selfing (1), and thus regardless of the fact that selfing species also harbour outcrossing morphs. Following [59–61], we justify this ‘any-instance’ coding, as (1) population-level estimates of morph frequencies are mostly lacking for these tropical orchids (see above); and (2) some assessment of selfing as a derived condition in clade C can be made a priori by comparison with other Madagascan Bulbophyllum species used as outgroups [47, 48].
Taxon sampling and DNA sequencing
Thirty species (65 accessions) out of the 32 species ascribed to sects. Calamaria (23/25 spp.), Humblotiorchis (1/1 sp.), and Bifalcula (6/6 spp.) were included in our study [Additional file 1). This taxon sampling comprises 15 species that were not sampled in , six yet undescribed species (‘spp. nov.’), and 28 species with known mating type, including eight selfers (Table 1). We were unable to obtain material for DNA extraction from the two remaining species of sect. Calamaria, i.e., B. sp. ‘F’ (with unknown mating type) and B. cryptostachium (Type I; ). We also included, as outgroups, one species each from sects. Kainochilus, Inversiflorum, and Alcistachys, whose close relationship to clade C is supported on molecular and morphological grounds [49, 50]. Like all extra-clade C species of Madagascan Bulbophyllum, these outgroups have so far proved to be monomorphic for the outcrossing mating type .
Most Bulbophyllum samples used in this study were field-collected in Madagascar, with all necessary permits obtained by the Département des Eaux et Fôrets (Madagascar), or received from botanic gardens (Parc Botanique et Zoologique de Tsimbazaza; Parc National de La Réunion; Vienna; Salzburg) and herbaria (REU, SZU, TAN, WU). Additional DNA samples were obtained from the Jodrell Laboratories, Kew. Whole genomic DNA was extracted from fresh, silica-dried leaf material using the 2x CTAB method . A few accessions, preserved as leaf extract on FTA® Mini Cards, were processed using the FTA® Starter Pack (Whatman, Florham Park, USA) according to the manufacturer’s specifications.
We PCR amplified and sequenced five plastid and three nuclear gene regions. The former comprised four intergenic spacers (atpI–atpH, psbA–trnH, trnD–trnE, trnT–trnS [including the rps4 gene]) and a hypothetical open reading frame gene region (ycf1). Nuclear regions were the internal transcribed spacer (ITS) region of ribosomal DNA (including the 5.8S gene) and portions of two low-copy genes, i.e., the structural gene encoding phosphoenolpyruvate carboxylase (PEPC; ), and the floral homoeotic gene pistillata/globosa (PI; ). For amplification we used protocols and primers from  for ITS, psbA–trnH, and trnD–trnE, [65, 66]) for atpI–atpH and trnT–trnS, and  for ycf1. Details of the newly designed primers for PEPC and PI are provided in [Additional file 2]. 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 as implemented in Geneious version 5.4 . For all plastid and nuclear regions, 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 1). We included outgroups and multiple intra-specific samples in the phylogenetic analyses and removed outgroups for the dating analysis in Beast version 1.6.1  (described further below), but removed each of those in all subsequent Beast, diversification, and mating-type transition analyses, resulting in a total of 30 terminal species. For these latter analyses, multiple accessions of the same species were pruned by keeping the most early-branching accession.
Four data partitions were defined, corresponding to each nuclear dataset (ITS, PEPC, PI) and the five plastid regions combined. We identified the best-fit models of nucleotide substitution for each partition using the Akaike Information Criterion (AIC) in jModeltest version 0.1.1 , i.e., the GTR + G model for the ITS and plastid datasets and the HKY + G model for PEPC and PI. These models were then applied to the phylogenetic analyses of each partition using Bayesian Inference (BI) as implemented in MrBayes version 3.1.2 , with B. variegatum (sect. Alcistachys) assigned as outgroup (note, MrBayes allows for only one outgroup sample). The Markov chain Monte Carlo (MCMC) algorithm was run with two independent runs of four Metropolis-coupled chains with 20 million generations each, with trees sampled every 1,000 generations. After assessing convergence between runs by monitoring the standard deviation of split frequencies in MrBayes and using the effective sampling size (ESS > 200) criterion in Tracer version 1.5 , the first 25 % of sampled trees were discarded as burn-in (based on stabilization of the log-likelihood of the trees). A 50 % majority-rule consensus tree was constructed from the remaining 15,000 trees to estimate posterior probability values for each node. For comparison, node support was also estimated by bootstrapping under maximum parsimony (MP) in Paup* version 4.0b10  using 1000 heuristic search replicates and the following settings: 100 random addition sequence replicates with three trees held at each step, tree bisection reconnection (TBR) branch swapping, and saving no more than 10 trees per replicate. The degree of incongruence among the four data partitions was checked by looking for conflicting clades that were strongly supported in terms of both Bayesian posterior probability (PP ≥ 0.95) and parsimony bootstrap percentage (BP ≥ 85) (see also [73, 74]). Since no such incongruence was detected see [Additional files 3, 4, 5 and 6], we concatenated all partitions into a single alignment, and re-analysed it under the above BI and MP settings.
Bayesian hypothesis testing of monophyletic selfers
Results of the monophyly tests for selfing Bulbophyllum clade C species using logarithmic Bayes Factor (BF) comparisons
Marginal tree likelihood
All selfers of clade C
All selfers of sects. Calamaria + Humblotiorchis
All selfers of sect. Calamaria
Evolutionary transitions in mating types
Based on the data of mating type (Table 1), we reconstructed transitions between outcrossing and selfing mating types in clade C using a two rate (Mk2) model of character evolution in the Bayesian (MCMC) framework of the submodule Multistate of BayesTraits version 1.0 (; http://www.evolution.rdg.ac.uk). Following [59–61], we designated the species as outcrossing or selfing using the ‘any-instance’ coding method (see above). Mating types were traced over the last 5,000 post-burn-in trees of the BI analysis to incorporate branch length and topology uncertainty. The MCMC analyses were run during 5,050,000 generations, with a ratedev parameter of 110 (obtained from initial test runs), a reversible-jump hyperprior with an exponential prior (mean seeded from a uniform distribution on the interval 0 to 30), and a burn-in of 50,000 generations.
Estimating times of divergence
Estimation of divergence times of clade C was performed in Beast with outgroup taxa removed and using the same unlinked gene partitions and substitution models as in the BI analysis (see above). A birth-death (BD) speciation model was specified as tree prior and an uncorrelated lognormal relaxed clock was assumed . Five independent MCMC runs were performed with 20 million generations each, sampling every 2,000 generations. Following the removal of a conservative burn-in of 15 %, the MCMC samples were combined using LogCombiner version 1.6.1  and inspected in Tracer to confirm convergence of the chain to stationary and assess sampling adequacy (ESS > 200). The results of the Beast analysis were summarized with TreeAnnotator version 1.6.1 . The remaining 42,500 trees were summarized onto the topology of 50% majority-rule consensus tree of the BI analysis with a threshold of zero for the posterior probability of clades and with node heights estimated using median values. The resulting chronogram (Fig. 3c) was visualized in Figtree version 1.3.1 (http://tree.bio.ed.ac.uk/software/figtree/). As there are no fossils of Bulbophyllum that can be used to assess the timing of divergence of clade C, we used a multi-step secondary calibration approach based on three fossils of other orchid genera and two geological calibration points taken from island ages. See [Additional file 7] for full details.
Course and mode of diversification
Results of fitting rate-constant and rate-variable diversification models to the maximum clade credibility chronogram of Bulbophyllum clade C (Fig. 3c) using birth-death likelihood (BDL) analysis
Estimates of state-dependent diversification/transitions rates
Median maximum likelihood (ML) estimates of parameters, log-likelihoods and AIC values of alternative BiSSE models estimated across 1,000 Beast trees of Bulbophyllum clade C. Speciation rates are λ, extinction rates are μ, and character transition rates are q. Outcrossing is coded as ‘O’ and selfing as ‘S’
% trees1 rejecting H0
q OS = 0
λ O = 0
λ O = λ S
q OS = q SO
q SO = 0
μ O = μ S
μ S = 0
μ O = 0
λ S = 0
λ S = 0, μ O = 0
λ S = 0, μ O = 0, q SO = 0
λ S = 0, μ O = 0, μ S = 0, q SO = 0
All models were fitted by ML nonlinear optimization across a sample of 1,000 of our Beast trees, using a heuristic starting point based (by default) on the state-independent BD model. For each model, ML parameter estimates, log-likelihoods and AIC values were recorded as median values (across the 1,000 trees), with the frequency distribution of full-model parameter values plotted along with their 95 % credible intervals. Model results were evaluated and compared using both median AIC values as well as likelihood ratio tests (LRTs). When differing by less than two ΔAIC units, models were considered being of essentially equivalent fit . For each tree, LRTs were calculated as twice the difference in log-likelihoods between the full (more complex) model and the constrained (simpler) model (H0), which should follow a χ 2 distribution with degrees of freedom equal to the difference in the number of free parameters allowed by the two models . We then recorded the percentage of trees where the full model provided a significant improvement (P < 0.05) over the constrained model.
Phylogenetics and hypothesis testing of monophyletic selfers
Our final dataset consisted of 670 bp from ITS, 765 bp of PEPC, 358 bp of PI, and 4712 bp of five plastid regions (atpI–atpH, psbA–trnH, trnD–trnE, trnT–trnS, ycf1). The combined nuclear/plastid dataset included 6,505 nucleotide positions, 477 of which were parsimony informative. The BI analysis of the four individual partitions [Additional files 3, 4, 5 and 6] yielded largely congruent tree topologies albeit with different levels of resolution. In our focal BI tree of the combined analysis (Fig. 2), which produced highest node support values, members of clade C grouped according to their sectional affiliation (all PP = 1; BP ≥ 99). The species-poor sect. Bifalcula and the monotypic sect. Humblotiorchis (comprising one selfing species each) formed successive sister taxa to the species-rich sect. Calamaria (including six selfers). While these sectional groupings and relationships are in agreement with previous phylogenetic studies , our increased taxon and gene sampling resulted in better-resolved interspecies relationships, with most intra-specific accessions recovered as monophyletic (PP > 0.95; BP 55–100), except for the paraphyletic B. implexum and the likely polyphyletic B. sp. nov. ‘E1’. However, especially within sect. Calamaria, several internal branches were not supported well (PP < 0.83; BP ≤ 50), indicating some phylogenetic uncertainty for character state reconstructions. Nonetheless, the scattered phylogenetic distribution of the eight selfing species strongly suggested that they do not form a monophyletic group (Fig. 2). In support of this, tree topology tests (Table 2) decisively rejected models in which all selfing species of clade C, or only those of its sublineages (sects. Calamaria + Humblotiorchis; or sect. Calamaria) were constrained to be monophyletic (all BF values > 10), pointing at multiple independent shifts in mating type across clade C (but see below).
Evolutionary transitions of mating type in clade C
Our Mk2 (MCMC) ancestral state reconstructions of mating type over 5,000 BI trees yielded mostly similar results when compared to those obtained with ML reconstructions over 1,000 Beast trees following the best fitting BiSSE model (Fig. 3a, b; Table 4). Both analyses identified outcrossing as the most probable ancestral state of clade C, and revealed that transitions to selfing were in general concentrated in terminal branches (marked with asterisks in Fig. 3a, b). However, under the Mk2 reconstruction, two subterminal nodes had relatively high probabilities of being selfing (i.e., the ancestral nodes of, respectively, B. occultum/B. bicoloratum/B. sp. nov. ‘D’, and B. pusillum/B.quadrifarium; Fig. 3a), resulting in six transitions. By contrast, the best fitting BiSSE model reconstructed all internal nodes as being outcrossing, resulting in eight independent transitions to selfing along terminal branches leading to extant species (Fig. 3b). In any event, these results qualitatively indicate that mating type transitions within clade C are a strikingly recurrent evolutionary phenomenon, but which requires further temporal delimitation and quantitative evidence as provided in the following sections.
Our relaxed molecular clock analyses (Fig. 3c, [Additional file 8]) dated the onset of diversification (crown age) of clade C to the Late Tertiary/Early Pliocene, at ca. 5.3 Ma (95 % highest posterior density [HPD] of 8.5–2.7 Ma). By contrast, the most recent common ancestor (MRCA) of sects. Humblotiorchis/Calamaria was placed in the mid-Pliocene, at ca. 4.4 Ma (95 % HPD of 5.3–3.4 Ma), and the crown age of sect. Bifalcula in the Late Pliocene, at ca. 3.1 Ma (95 % HPD of 4.2–2.0 Ma). The origins of most extant clade C species were dated to the Quaternary (≤ 2.6 Ma), excepting few species inferred as Late Tertiary/Pliocene relics (i.e., B. erectum, B. humblotii, B. capuronii). Notably, five of the eight selfing species originated within a relatively recent time frame starting at about the mid-Pleistocene (ca. 1.3 Ma; see node marked by an asterisk in Fig. 3c; 95 % HPD of 1.9–0.8 Ma). By contrast, of the three remaining selfers, one originated in the Early Pleistocene (B. complanatum [ca. 2.4 Ma]) and two in the Pliocence (B. erectum [ca. 3.1 Ma] and B. humblotii [ca. 4.37 Ma]; Fig. 3c). All these latter taxa evolved selfing independently (Fig. 3), leaving the possibility of a more complex history of selfing (in terms of speciation, extinction and/or transitions) over extended time scales.
Diversification and transition rates
The LTT plots for the phylogeny of clade C (inset Fig. 3c), derived from the chronogram (Fig. 3c) and 1,000 dated Beast trees, suggested an almost constant lineage accumulation over time, with no apparent upturn near the present as typical for either high (background) extinction rates or recent and rapid speciation rates [82, 85]. Concomitantly, the ∆AICRC test in Laser (Table 3) indicated no significant departure of the best-fit rate-constant model, i.e., pure-birth (PB) without extinction (AIC = -48.97), from the best-fit rate-variable model (YR3) with likewise no extinction but three distinct rates of diversification (AIC = -52.22). The overall speciation rate (λ) of clade C calculated under the null (PB) model was 0.53 species myr-1. This value is within the range observed for other plant taxa that have been identified as undergoing rapid rates of relatively recent speciation (e.g., [27, 88, 89]).
According to recent reformulations [7, 14], Stebbins’ SEDE hypothesis conflates two distinct claims: (1) transition rates from selfing to outcrossing are zero; and (2) the net diversification rate (speciation rate minus extinction rate) is negative for selfing taxa, implying a higher risk of extinction (i.e., shorter evolutionary lifespan) for selfing relative to outcrossing taxa [2, 17, 20, 25]. Our results for Madagascan Bulbophyllum clade C support the first claim of the SEDE hypothesis (irreversibility of transitions to selfing) and we find some, albeit no definite evidence for the second (higher extinction rates for selfers). In principle, therefore, our data are largely congruent with predictions raised by the SEDE hypothesis , but some of the patterns observed suggest that selfing represents a starting-point rather than a terminus of mating-system evolution in these tropical orchids. Below we discuss the evidence and possible causes for the strikingly recurrent and probably very recent transitions to selfing in our study system.
Multiple origins and irreversibility of selfing
We support the first assumption of the SEDE hypothesis in that outcrossing (rostellum presence, Type I) represents the ancestral state of Bulbophyllum clade C (Fig. 3) and that transitions exclusively occurred from outcrossing to selfing (Types II/III), as also confirmed by the state-dependent BiSSE analyses (Table 4). Although our ancestral state reconstructions using the Mk2 model indicate two subterminal nodes where transitions could have occurred (Fig. 3a), this is not supported by the ML reconstructions under the best fitting BiSSE model (Fig. 3b, c). Accordingly, selfing via rostellum loss (Type II) most likely evolved independently at eight (rather than only six) times within this clade of 30 species. In addition, selfing via rostellum receptivity (Type III) evolved uniquely in the trimorphic B. bicoloratum [47, 48]. These results also imply that the sharing of selfing/outcrossing morph variation among the eight focal species is unlikely due to the retention of an ancestral polymorphism, and the consistent rejection of their monophyly (Table 2) further supports this view.
Unidirectional transitions from outcrossing to selfing are commonly reported in flowering plants [3, 18–20, 24, 26, 28] but see [25, 27]. However, this is one of few studies to estimate a time-calibrated transition rate to selfing based on a discrete-state character. The only comparable studies we are aware of, focussing on the breakdown of SI in Solanaceae  and the loss of heterostyly in Primula , reported higher and lower rates than observed here (i.e., q OS 0.56 and 0.04–0.10, respectively, vs. 0.21 [this study], all in units myr-1). We note that BiSSE is currently the best method to test character state irreversibility, as it is less prone to falsely infer reversions [7, 14], and has been used previously for clades of similarly small size (e.g., [37, 78, 90–93]). However, it has relatively low power to detect rate differences among states of a binary trait when dealing with topologies of less than about 300 terminals and a high degree of asymmetry in tip states [33–35]. Our clade C phylogeny has 30 tips of which 8 (26.7 %) have state one (i.e., selfing). The results of our BiSSE analysis must therefore be treated cautiously (see also below).
The fact that selfing has arisen multiple, independent times within clade C may imply a selective advantage of this trait ([55, 92, 94]; see further below), and that inbreeding depression, as the principal counterweight to this advantage (e.g., [13, 14]), has not impeded the evolution of selfing in these tropical orchids. In support of the latter hypothesis, controlled mating experiments have found no evidence of inbreeding depression at the level of seed viability and fruit set in, respectively, the outcrossing morph of Madagascan B. occultum (U. Jaros, unpubl. data) and neotropical Bulbophyllum species . There is also little evidence of such early-acting (pre-dispersal) inbreeding depression in other self-compatible orchid species [96–98]. At least in part, this may be explained by the fact that selfed offspring usually does not depend on parental (i.e., endosperm) resources in Orchidaceae [43, 99].
But what intrinsic factor(s) then could promote the parallel evolution of selfing via rostellum loss within clade C? One possible genetic explanation for this parallelism is that it ultimately reflects recurrent, perhaps loss-of-function mutations at only one or a few genes affecting the timing of rostellum development (see also ). Morphological divergence within and between plant species involving the presence versus absence of structures, as observed here for rostellum type, has often been demonstrated to have such a simple genetic basis [100–103], and selection on such traits should lead to rapid evolutionary response [101, 103, 104]. On the other hand, given the morphological and developmental complexity of the rostellum [47, 48, 105], it is not unexpected that we find no evidence for the regain of this organ. In general, such kind of irreversibility is expected from ‘Dollo’s Law’, which posits that complex characters, once lost, are unlikely to re-evolve [2, 106], as also demonstrated for other morphological structures associated with outcrossing (e.g., heterostyly: [2, 20] but see ). However, at least for the orchid species studied here, a non-mutually exclusive explanation is that transitions to selfing evolved too recently to detect reversals (see below).
Explanations for the ‘twiggy’ phylogenetic distribution of selfing
It has long been noted that selfing taxa often occur in scattered positions at the tips of plant phylogenies (e.g., Amsinckia: ; Triticeae: ; Leavenworthia: [22, 23]. Such’twiggy’ (or ‘tippy’) distributions have typically been interpreted as a hallmark signature of the elevated extinction risk of selfers, as proposed by the second part of the SEDE hypothesis [15, 17]. However, these’twiggy’ distributions can also arise if transitions from the ancestral state are rare or asymmetrical, but this has rarely been tested using sophisticated analysis methods, such as the state-dependent BiSSE model [40, 41, 108, 109], and with partly contrasting results. For example, in Primula , this latter model proved useful in explaining the phylogenetic ‘twigginess’ of the selfing state by highly asymmetrical transitions to selfing from outcrossing, albeit coupled with increased speciation and zero extinction (positive r). Interestingly, these results demonstrate that selfing lineages, perhaps freed from selective constraints, can evolve at higher rates than outcrossing ones, but also suggests that sufficient time is needed for the negative effects of selfing to accumulate (; J. de Vos, pers. comm.).
For Bulbophyllum, our data best fit two BiSSE models in which the outcrossing state is generally associated with high rates of speciation and zero extinction (λ O = 0.68/071; μ O = 0), and the selfing state with zero speciation (λ S = 0) (Table 4); however, we did not find definite evidence for the prediction of increased extinction of selfing species, even though the second-best model with moderate extinction (μ S = 0.25; and thus negative r) fit reasonably well compared to the best-fit model without extinction (μ S = 0) (ΔAIC < 2; Table 4). This lack of ‘measurable’ extinction may simply reflect the above caveat about the statistical power of BiSSE, coupled with the relatively short branch lengths subtending the selfing species of Bulbophyllum (Fig. 3c). In sum, the results of our BiSSE analyses must be regarded as inconclusive on whether or not extinction contributes to the twiggy phylogenetic distribution of selfing in clade C (Fig. 3), and thus in addition to the repeated (and irreversible) transition to this character, coupled with strongly reduced (zero) speciation in selfing compared to outcrossing taxa. That said, if selfing had evolved only recently in this group, as discussed below, then time might have been insufficient to detect elevated extinction and/or speciation rates associated with this state ([40, 56]; see also ).
Relatively recent shifts in mating system
There is mounting evidence from population genetic studies indicating that outcrossing-to-selfing transitions are often rather recent in flowering plants [7, 14], with most origins dated to Quaternary (≤ 2.6 Ma) or even millennial time scales [22, 110–115]. Our results lend support to this view in that selfing evolved independently, in all probability at five times, within clade C since about the mid-Pleistocene (ca. ≤ 1.3 Ma; see node marked by an asterisk in Fig. 3c). However, given that all eight transitions are most likely associated with terminal branches leading to extant species (Fig. 3c), it remains unclear when selfing evolved within a particular species. This could have happened at any point along such a branch, with the possibility in three species (i.e., B. complanatum, B. erectum, B. humblotii) of tracing back to even Early Pleistocene/Pliocene times. Nonetheless, for each of the eight focal species, there is circumstantial evidence to suggest that selfing evolved only recently, possibly over Late Quaternary, if not historic time scales.
First, molecular genetic data in particular from Brassicaceae indicate that selfing tends to evolve much more recently within a species than the time of separation from its closest extant relative [13, 113]. For example, the common ancestor of Arabidopsis thaliana and A. lyrata (which likely diverged ca. 5 Ma  or 8.0–17.9 Ma ) was presumably outcrossing and self-incompatible, yet A. thaliana is thought to have acquired self-compatibility only within the last 1.8 million years [113–115]. An even more recent origin of selfing seems plausible for each of the eight Bulbophyllum species, where – unlike in A. thaliana – selfing individuals still co-exist with ancestrally outcrossing ones at the intraspecific level. In addition, apart from their diagnostic rostellum traits, these conspecific morphs hardly differ in other morphological or functional changes to the flowers [47, 48]. This further suggests that these morphs only recently diverged, whereby insufficient time may have elapsed for such a ‘floral selfing syndrome’  to develop, as has been demonstrated in a large number of genera where selfing lineages have been derived from outcrossing ancestors (reviewed in ; but see ). On the other hand, if selfing had evolved in the more distant past in Bulbophyllum, we would expect to observe at least some of the older species (e.g., B. humblotii) to be entirely fixed for selfing, especially given the evidence for low inbreeding depression in this group (see above). Accordingly, the eight clade C species that currently vary for mating type may well be in a transitory phase towards increased selfing, as postulated for other species polymorphic for selfing/outcrossing individuals [55, 112].
Potential environmental catalysts for the evolution of selfing
Observations on cultivated material suggest that the selfing morphs (Types II/III) of Madagascan Bulbophyllum species are unlikely to reproduce through outcrossed pollinia, and thus most likely behave as exclusive rather than partial selfers . If this holds true in nature, then the origin and maintenance of selfing within clade C is unlikely due to a mere ‘selective transmission advantage’ of selfing [5–7], because this advantage disappears when pollen discounting (i.e., a reduction in male outcross success that results from selfing) is maximised . It is feasible, therefore, that recurrent transitions to selfing within this clade mostly, if not exclusively, reflect an adaptive ‘reproductive assurance strategy’ [8–10] in response to similar selective environments that favour reproduction when pollinators or mates are rare or absent, conditions that might accrue from island colonization or habitat fragmentation [7, 119, 120]. In our system, however, we find no support for the oft-invoked idea that the evolution of selfing promotes long-distance colonization [2, 8, 15, 121]. Although in some species (B. humblotii, B. occultum, B. pusillium) selfing may have increased colonizing ability of neighbouring islands, such as La Réunion, all these taxa vary for mating type in Madagascar [48, 54]. On the other hand, we suspect that potential catalysts that may have initiated fragmented habitat conditions in Madagascar include (1) climate-induced vegetation shifts during the (Late) Quaternary [122–126]; and (2) human-mediated degradation of primary forest over the last 1,800 years [127–129]. Effects of both Quaternary climate change and human-modified landscapes on the evolution of selfing in flowering plants are well documented [22, 104, 111, 114, 130]. Our phylogenetic data cannot determine the relative importance of these ‘shallow-time’ effects to mating system evolution in Bulbophyllum clade C; however, they emphasize the need for comprehensive field and population-level data (e.g., on morph frequencies, pollinators, demography, niche requirements, phylogeography, and genomics) to elucidate the past and/or present selective environments under which selfing evolved within a particular Bulbophyllum species .
In this study, we examined the evolution of selfing (auto-pollination) in Madagascan Bulbophylum orchids (clade C) in the light of two predictions raised by the classical SEDE hypothesis: (1) irreversible transitions from outcrossing to selfing; and (2) elevated extinction rates associated with the selfing state. Although we find convincing support for the first prediction, results remain inconclusive on whether or not extinction contributes to the twiggy phylogenetic distribution of selfing in clade C. Despite this latter limitation, we suspect that selfing has evolved too recently in this group to detect elevated extinction and/or speciation rates associated with this state. We further hypothesize that a simple genetic basis of rostellum loss may explain the strikingly recurrent transitions to selfing in our study system (i.e., exclusively within species), perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales. Future demographic and genomic studies of these and other tropical orchids with an apparent lability towards rostellum abortion have thus the potential to unravel novel perspectives on the ability of some natural systems to rapidly evolve selfing as a means of avoiding extinction in the face of ever changing environments.
The data sets supporting the results of this article are included within the article (and its Additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 and 13). Nucleotide sequence data supporting the results of this article are available in GenBank [see Additional file 1]. The phylogenetic data sets supporting the results of this article are available in the Dryad Digital Repository [doi:10.5061/dryad.35935].
We thank members of the Botanical Garden of Vienna University (M. Kiehn, A. Sieder), the Parc Botanique et Zoologique de Tsimbazaza (S. Rapanarivo, J. Andriantiana), and the Université de la Réunion (T. Pailler) for general support and help with fieldwork in Madagascar and La Réunion, the Département des Eaux et Fôrets (Madagascar) and the Parc National de La Réunion for issuing collecting permits, the curators of REU, SZU, TAN, and WU for the loan of Bulbophyllum specimens, the Jodrell Laboratories (Kew) for providing DNA samples, and G. Hochschartner for making available unpublished sequences. We also thank G. Brunauer and M. Affenzeller for lab support and E. Egger and her team for Bulbophyllum cultivation at Salzburg Botanical Garden. We extend special thanks to R.J. Abbott (University of St. Andrews), H.P. Linder (University of Zürich.) and H. Schäfer (Technical University Munich) for stimulating discussions on this research topic, and to J. de Vos (Brown University) and three anonymous reviewers for very insightful comments on an earlier version of this manuscript. This study was funded by the FWF (Austrian Science Fund) grants P20726-B03 and P17124-B0 to H.P.C. and G.A.F., respectively. All authors disclose any potential sources of conflict of interest.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Takhtajan A. Evolutionary trends in flowering plants. New York: Columbia Univ Press; 1991.Google Scholar
- Barrett SCH. The evolution of plant reproductive systems: how often are transitions irreversible? Proc R Soc B. 2013;280:20130913.PubMed CentralPubMedView ArticleGoogle Scholar
- Barrett SCH. The evolution of plant sexual diversity. Nat Rev Genet. 2002;3:274–84.PubMedView ArticleGoogle Scholar
- Barrett SCH, Eckert CG. Variation and evolution of mating systems in seed plants. In: Kawano S, editor. Biological approaches and evolutionary trends in plants. New York: Academic; 1990. p. 229–54.Google Scholar
- Fisher RA. Average excess and average effect of a gene substitution. Ann Eugenics. 1941;11:53–63.View ArticleGoogle Scholar
- Holsinger KE. Reproductive systems and evolution in vascular plants. Proc Natl Acad Sci U S A. 2000;97:7037–42.PubMed CentralPubMedView ArticleGoogle Scholar
- Wright SI, Kalisz S, Slotte T. Evolutionary consequences of self-fertilization in plants. Proc R Soc B. 2013;280:20130133.PubMed CentralPubMedView ArticleGoogle Scholar
- Baker HG. Self-compatibility and establishment after ‘long distance’ dispersal. Evolution. 1955;9:347–8.View ArticleGoogle Scholar
- Baker HG. Support for Baker’s law as a rule. Evolution. 1967;21:853–6.View ArticleGoogle Scholar
- Cheptou P-O. Allee effect and self-fertilization in hermaphrodites: reproductive assurance in demographically stable populations. Evolution. 2004;58:2613–21.PubMedView ArticleGoogle Scholar
- Knight T. Experiments on the fecundation of vegetables. Phil Trans R Soc B. 1799;89:195–204.View ArticleGoogle Scholar
- Darwin C. The effects of cross and self fertilisation in the vegetable kingdom. London: John Murray; 1876.Google Scholar
- Charlesworth D. Evolution of plant breeding systems. Curr Biol. 2006;16:R726–35.PubMedView ArticleGoogle Scholar
- Igic B, Busch JW. Is self-fertilization an evolutionary dead end? New Phytol. 2013;198:386–97.PubMedView ArticleGoogle Scholar
- Stebbins GL. Self-fertilization and population variability in higher plants. Am Nat. 1957;41:337–54.View ArticleGoogle Scholar
- Stebbins GL. Flowering plants: evolution above the species level. Cambridge: Belknap; 1974.View ArticleGoogle Scholar
- Takebayashi N, Morrell PL. Is self-fertilization an evolutionary dead end? Revisiting an old hypothesis with genetic theories and a macroevolutionary approach. Am J Bot. 2001;88:1143–50.PubMedView ArticleGoogle Scholar
- Barrett SCH, Harder LD, Worley AC. The comparative biology of pollination and mating in flowering plants. Phil Trans R Soc B. 1996;351:1271–80.View ArticleGoogle Scholar
- Kohn JR, Graham SW, Morton B, Doyle JJ, Barrett SCH. Reconstruction of the evolution of reproductive characters in Pontederiaceae using phylogenetic evidence from chloroplast DNA restriction-site variation. Evolution. 1996;50:1454–69.View ArticleGoogle Scholar
- Schoen DJ, Johnston MO, L’Heureux A-M, Marsolais JV. Evolutionary history of the mating system in Amsinckia (Boraginaceae). Evolution. 1997;51:1090–9.View ArticleGoogle Scholar
- Escobar JS, Cenci A, Bolognini J, Haudry A, Laurent S, David J, et al. An integrative test of the dead-end hypothesis of selfing evolution in Triticeae (Poaceae). Evolution. 2010;64:2855–72.PubMedGoogle Scholar
- Busch JW, Joly S, Schoen DJ. Demographic signatures accompanying the evolution of selfing in Leavenworthia alabamica Jeremiah. Mol Biol Evol. 2011;28:1717–29.PubMedView ArticleGoogle Scholar
- Herman AC, Busch JW, Schoen DJ. Phylogeny of Leavenworthia S-alleles suggests unidirectional mating system evolution and enhanced positive selection following an ancient population bottleneck. Evolution. 2012;66:1849–61.PubMedView ArticleGoogle Scholar
- Armbruster WS. Evolution of plant pollination systems: hypotheses and tests with the neotropical vine Dalechampia. Evolution. 1993;47:1480–505.View ArticleGoogle Scholar
- Bena G, Lejeune B, Prosperi J-M, Olivieri I. Molecular phylogenetic approach for studying life-history evolution: the ambiguous example of the genus Medicago L. Proc R Soc Lond B. 1998;265:1141–51.View ArticleGoogle Scholar
- Goodwillie C. Multiple origins of self-compatibility in Linanthus section Leptosiphon (Polemoniaceae): phylogenetic evidence from internal transcribed-spacer sequence data. Evolution. 1999;53:1387–95.View ArticleGoogle Scholar
- Ferrer MM, Good-Avila SV. Macrophylogenetic analyses of the gain and loss of self-incompatibility in the Asteraceae. New Phytol. 2006;173:401–14.View ArticleGoogle Scholar
- Igic B, Bohs L, Kohn JR. Ancient polymorphism reveals unidirectional breeding system shifts. Proc Natl Acad Sci U S A. 2006;103:1359–63.PubMed CentralPubMedView ArticleGoogle Scholar
- Domes K, Norton RA, Maraun M, Scheu S. Reevolution of sexuality breaks Dollo’s law. Proc Natl Acad Sci U S A. 2007;104:7139–44.PubMed CentralPubMedView ArticleGoogle Scholar
- Pagel M. Inferring the historical patterns of biological evolution. Nature. 1999;401:877–84.PubMedView ArticleGoogle Scholar
- Barraclough TG, Birky CW, Burt A. Diversification in sexual and asexual organisms. Evolution. 2003;57:2166–72.PubMedView ArticleGoogle Scholar
- Maddison WP. Confounding asymmetries in evolutionary diversification and character change. Evolution. 2006;60:1743–6.PubMedView ArticleGoogle Scholar
- Maddison WP, Midford PE, Otto SP. Estimating a binary character’s effect on speciation and extinction. Syst Biol. 2007;56:701–10.PubMedView ArticleGoogle Scholar
- Rabosky DL. Extinction rates should not be estimated from molecular phylogenies. Evolution. 2010;64:1816–24.PubMedView ArticleGoogle Scholar
- Davis M, Midford P, Maddison W. Exploring power and parameter estimation of the BiSSE method for analyzing species diversification. BMC Evol Biol. 2013;13:38.PubMed CentralPubMedView ArticleGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.PubMed CentralPubMedView ArticleGoogle Scholar
- Goldberg EE, Igic B. On phylogenetic tests of irreversible evolution. Evolution. 2008;62:2727–41.PubMedView ArticleGoogle Scholar
- FitzJohn RG, Maddison WP, Otto SP. Estimating trait-dependent speciation and extinction rates from incompletely resolved phylogenies. Syst Biol. 2009;58:595–611.PubMedView ArticleGoogle Scholar
- Goldberg EE, Kohn JR, Lande R, Robertson KA, Smith SA, Igic B. Species selection maintains self-incompatibility. Science. 2010;330:493–5.PubMedView ArticleGoogle Scholar
- Johnson MTJ, FitzJohn RG, Smith SD, Rausher MD, Otto SP. Loss of sexual recombination and segregation is associated with increased diversification in evening primroses. Evolution. 2011;65:3230–40.PubMedView ArticleGoogle Scholar
- Goldberg EE, Igic B. Tempo and mode in plant breeding system evolution. Evolution. 2012;66:3701–9.PubMedView ArticleGoogle Scholar
- Sieder A, Rainer H, Kiehn M: 2007 CITES checklist for Bulbophyllum and allied taxa (Orchidaceae). [www.cites.org/common/com/nc/tax_ref/Bulbophyllum.pdf]
- Dressler RL. The orchids: natural history and classification. Cambridge: Harvard Univ Press; 1981.Google Scholar
- van der Cingel NA. An atlas of orchid pollination America, Africa, Asia and Australia. Rotterdam: Balkema; 2001.Google Scholar
- Catling PM. Auto-pollination in the Orchidaceae. In: Arditti J, editor. Orchid biology: reviews and perspectives v. Portland: Timber Press; 1990. p. 121–58.Google Scholar
- Borba EL, Semir J. Temporal variation in pollinarium size after its removal in species of Bulbophyllum: a different mechanism preventing self-pollination in Orchidaceae. Plant Syst Evol. 1999;217:197–204.View ArticleGoogle Scholar
- Gamisch A, Staedler YM, Schönenberger J, Fischer GA, Comes HP. Histological and micro-CT evidence of stigmatic rostellum receptivity promoting auto-pollination in the Madagascan orchid Bulbophyllum bicoloratum. PLoS One. 2013;8, e72688.PubMed CentralPubMedView ArticleGoogle Scholar
- Gamisch A, Fischer GA, Comes HP. Recurrent polymorphic mating type variation in Madagascan Bulbophyllum species exemplifies high incidence of auto-pollination in tropical orchids. Bot J Linn Soc. 2014;175:242–58.PubMed CentralPubMedView ArticleGoogle Scholar
- Fischer GA. Evolution of the orchid genus Bulbophyllum in Madagascar. PhD thesis. University of Vienna: Austria; 2007.Google Scholar
- Fischer GA, Gravendeel B, Sieder A, Andriantiana J, Heiselmayer P, Cribb PJ, et al. Evolution of resupination of Madagascan species of Bulbophyllum (Orchidaceae). Mol Phylogenet Evol. 2007;45:358–76.PubMedView ArticleGoogle Scholar
- Peter CI. Pollinators floral deception and evolutionary processes in Eulophia (Orchidaceae) and its allies. PhD thesis. University of KwaZulu-Natal: Republic of South Africa; 2009.Google Scholar
- Peter CI, Johnson SD. Autonomous self-pollination and pseudo-fruit set in South African species of Eulophia (Orchidaceae). S Afr J Bot. 2009;75:791–7.View ArticleGoogle Scholar
- Micheneau C, Johnson SD, Fay MF. Orchid pollination: from Darwin to the present day. Bot J Linn Soc. 2009;161:1–19.View ArticleGoogle Scholar
- Jacquemyn H, Micheneau C, Roberts DL, Pailler T. Elevational gradients of species diversity breeding system and floral traits of orchid species on Réunion Island. J Biogeogr. 2005;32:1751–61.View ArticleGoogle Scholar
- Busch JW, Urban L. Insights gained from 50 years of studying the evolution of self-compatibility in Leavenworthia (Brassicaceae). Evol Biol. 2011;38:15–27.View ArticleGoogle Scholar
- de Vos JM, Hughes CE, Schneeweiss GM, Moore BR, Conti E. Heterostyly accelerates diversification via reduced extinction in primroses. Proc R Soc B. 2014;281:20140075.PubMed CentralPubMedView ArticleGoogle Scholar
- Pridgeon AM, Cribb PJ, Chase MW, Rasmussen FN. Genera Orchidacearum Vol 6 Epidendroideae (Part 3). Oxford: Oxford University Press; 2014.Google Scholar
- Cribb P, Hermans J. Field guide to the orchids of Madagascar. Royal Botanic Gardens Kew: Richmond; 2009.Google Scholar
- Wiens JJ. Polymorphism in systematics and comparative biology. Annu Rev Ecol Syst. 1999;30:327–62.View ArticleGoogle Scholar
- Wiens JJ. Coding morphological variation within species and higher taxa for phylogenetic analysis. In: Wiens JJ, editor. Phylogenetic analysis of morphological data. Washington: Smithsonian Institute Press; 2000. p. 115–45.Google Scholar
- Dávalos LM, Cirranello AL, Geisler JH, Simmons NB. Understanding phylogenetic incongruence: lessons from phyllostomid bats. Biol Rev. 2012;87:991–1024.PubMed CentralPubMedView ArticleGoogle Scholar
- Doyle JJ, Doyle JL. A rapid DNA isolation procedure for small quantities of fresh leaf material. Phytochem Bull. 1987;19:11–5.Google Scholar
- Gehrig H, Heute V, Kluge M. New partial sequences of phosphoenolpyruvate carboxylase as molecular phylogenetic markers. Mol Phylogenet Evol. 2001;20:262–74.PubMedView ArticleGoogle Scholar
- Schlüter PM, Kohl G, Stuessy TF, Paulus HF. A screen of low-copy nuclear genes reveals the LFY gene as phylogenetically informative in closely related species of orchids (Ophrys). Taxon. 2007;56:493–504.Google Scholar
- Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, et al. The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am J Bot. 2005;92:142–66.PubMedView ArticleGoogle Scholar
- Shaw J, Lickey EB, Schilling EE, Small RL. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. Am J Bot. 2007;94:275–88.PubMedView ArticleGoogle Scholar
- Neubig KM, Whitten WM, Carlsward BS, Blanco MA, Endara L, Williams NH, et al. Phylogenetic utility of ycf1 in orchids: a plastid gene more variable than matK. Plant Syst Evol. 2009;277:75–84.View ArticleGoogle Scholar
- Drummond AJB, Ashton S, Buxton M, Cheung A, Cooper C, Duran M, et al. 2011 Geneious v.5.4. [http://www.geneious.com/]
- Posada D. JModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.PubMedView ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.PubMedView ArticleGoogle Scholar
- Rambaut A, Suchard MA, Xie D, Drummond AJ. 2013 Tracer v.1.5. [http://treebioedacuk/software/tracer/]
- Swofford DL. PAUP*: phylogenetic analysis using parsimony (* and other methods) version 4.0b10. Sunderland: Sinauer Associates; 2001.Google Scholar
- Norup MV, Dransfield J, Chase MW, Barfod AS, Fernando ES, Baker WJ. Homoplasious character combinations and generic delimitation: a case study from the Indo-Pacific arecoid palms (Arecaceae: Areceae). Am J Bot. 2006;93:1065–80.PubMedView ArticleGoogle Scholar
- Salvo G, Bacchetta G, Ghahremaninejad F, Conti E. Phylogenetic relationships of Ruteae (Rutaceae): new evidence from the chloroplast genome and comparisons with non-molecular data. Mol Phylogenet Evol. 2008;49:736–48.PubMedView ArticleGoogle Scholar
- Kass RE, Raftery AE. Bayes factors. J Am Stat Assoc. 1995;90:773–95.View ArticleGoogle Scholar
- Suchard MA, Weiss RE, Sinsheimer JS. Bayesian selection of continuous time Markov chain evolutionary models. Mol Biol Evol. 2001;18:1001–13.PubMedView ArticleGoogle Scholar
- Pagel M, Meade A, Barker D. Bayesian estimation of ancestral character states on phylogenies. Syst Biol. 2004;53:673–84.PubMedView ArticleGoogle Scholar
- Lynch VJ, Wagner GP. Did egg-laying boas break Dollo’s law? Phylogenetic evidence for reversal to oviparity in sand boas (Eryx: Boidae). Evolution. 2010;64:207–16.PubMedView ArticleGoogle Scholar
- Silvestro D, Zizka G, Schulte K. Disentangling the effects of key innovations on the diversification of Bromelioideae (Bromeliaceae). Evolution. 2014;68:163–75.PubMedView ArticleGoogle Scholar
- Drummond AJ, Nicholls GK, Rodrigo AG, Solomon W. Estimating mutation parameters population history and genealogy simultaneously from temporally spaced sequence data. Genetics. 2002;161:1307–20.PubMed CentralPubMedGoogle Scholar
- Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.PubMedView ArticleGoogle Scholar
- Nee S, May RM, Harvey PH. The reconstructed evolutionary process. Phil Trans R Soc Lond B. 1994;344:305–11.View ArticleGoogle Scholar
- Rabosky DL. LASER: a maximum likelihood toolkit for detecting temporal shifts in diversification rates from molecular phylogenies. Evol Bioinform Online. 2006;2:247–50.PubMed CentralGoogle Scholar
- Rabosky DL, Lovette IJ. Problems detecting density-dependent diversification on phylogenies: reply to Bokma. Proc R Soc B. 2009;276:995–7.PubMed CentralView ArticleGoogle Scholar
- Rabosky DL. Likelihood methods for inferring temporal shifts in diversification rates. Evolution. 2006;60:1152–64.PubMedView ArticleGoogle Scholar
- Cabezas P, Sanmartín I, Paulay G, Macpherson E, Machordom A. Deep under the sea: unraveling the evolutionary history of the deep-sea squat lobster Paramunida (Decapoda Munididae). Evolution. 2012;66:1878–96.PubMedView ArticleGoogle Scholar
- Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.Google Scholar
- Baldwin BG, Sanderson MJ. Age and rate of diversification of the Hawaiian silversword alliance (Compositae). Proc Natl Acad Sci U S A. 1998;95:9402–6.PubMed CentralPubMedView ArticleGoogle Scholar
- Egan AN, Crandall KA. Divergence and diversification in North American Psoraleeae (Fabaceae) due to climate change. BMC Biol. 2008;6:55.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith SDW, Miller RE, Otto SP, FitzJohn RG, Rausher MD. The effects of flower color transitions on diversification rates in morning glories (Ipomoea subg. Quamoclit Convolvulaceae). In: Long M, Gu H, Zhou Z, editors. Darwin's Heritage Today. Bejing: Higher Education Press; 2010. p. 202-226 Google Scholar
- Mayrose I, Zhan SH, Rothfels CJ, Magnuson-Ford K, Barker MS, Rieseberg LH, et al. Recently formed polyploid plants diversify at lower rates. Science. 2011;333:1257.PubMedView ArticleGoogle Scholar
- Fernández-Mazuecos M, Blanco-Pastor JL, Gómez JM, Vargas P. Corolla morphology influences diversification rates in bifid toadflaxes (Linaria sect, Versicolores). Ann Bot. 2013;112:1705–22.PubMed CentralPubMedView ArticleGoogle Scholar
- Breitkopf H, Onstein RE, Cafasso D, Schluter PM, Cozzolino S. Multiple shifts to different pollinators fuelled rapid diversification in sexually deceptive Ophrys orchids. New Phytol. 2014;207:377–89.PubMedView ArticleGoogle Scholar
- Wiens JJ, Brandley MC, Reeder TW. Why does a trait evolve multiple times within a clade? Repeated evolution of snakelike body form in squamate reptiles. Evolution. 2006;60:123–41.PubMedGoogle Scholar
- Borba EL, Sheppard GJ, Semir J. Reproductive systems and crossing potential in three species of Bulbophyllum (Orchidaceae) occurring in Brazilian ‘campo rupestre’ vegetation. Plant Syst Evol. 1999;217:205–14.View ArticleGoogle Scholar
- Ortiz-Barney E, Ackerman JD. The cost of selfing in Encyclia cochleata (Orchidaceae). Plant Syst Evol. 1999;219:55–64.View ArticleGoogle Scholar
- Tremblay RL, Ackerman JD, Zimmerman JK, Calvo RN. Variation in sexual reproduction in orchids and its evolutionary consequences: a spasmodic journey to diversification. Biol J Linn Soc. 2005;84:1–54.View ArticleGoogle Scholar
- Kuehnle AR. Orchids. In: Anderson NO, editor. Flower breeding and genetics. Dordrecht: Springer; 2007. p. 539–60.Google Scholar
- Harder LD, Richards SA, Routley MB. Effects of reproductive compensation gamete discounting and reproductive assurance on mating-system diversity in hermaphrodites. Evolution. 2007;62:167–72.View ArticleGoogle Scholar
- Gottlieb LD. Genetics and morphological evolution in plants. Am Nat. 1984;123:681–709.View ArticleGoogle Scholar
- Coyne JA, Lande R. The genetic basis of species differences in plants. Am Nat. 1985;126:141–5.View ArticleGoogle Scholar
- Kim M, Cui ML, Cubas P, Gillies A, Lee K, Chapman MA, et al. Regulatory genes control a key morphological and ecological trait transferred between species. Science. 2008;322:1116–9.PubMedView ArticleGoogle Scholar
- Lankinen A. Upper petal lip colour polymorphism in Collinsia heterophylla (Plantaginaceae): genetic basis within a population and its use as a genetic marker. J Genet. 2009;88:205–15.PubMedView ArticleGoogle Scholar
- Foxe JP, Slotte T, Stahl EA, Neuffer B, Hurka H, Wright SW. Recent speciation associated with the evolution of selfing in Capsella. Proc Natl Acad Sci U S A. 2009;106:5241–5.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurzweil H, Kocyan A. Ontogeny of orchid flowers. In: Arditti J, editor. Orchid biology: reviews and perspectives. Dordrecht: Kluwer; 2002. p. 2–138.Google Scholar
- Bull JJ, Charnov EL. On irreversible evolution. Evolution. 1985;62:2727–41.Google Scholar
- McDill J, Repplinger M, Simpson BB, Kadereit JW. The phylogeny of Linum subfamily Linoideae with implications for their systematics biogeography and evolution of heterostyly. Syst Bot. 2009;34:386–405.View ArticleGoogle Scholar
- Schwander T, Crespi BJ. Twigs on the tree of life? Neutral and selective models for integrating macroevolutionary patterns with microevolutionary processes in the analysis of asexuality. Mol Ecol. 2009;18:28–42.PubMedView ArticleGoogle Scholar
- Ng J, Smith SD. How traits shape trees: new approaches for detecting character state-dependent lineage diversification. J Evol Biol 2014, doi:101111/jeb12460.Google Scholar
- Tang C, Toomajian C, Sherman-Broyles S, Plagnol V, Guo Y-L, Hu TT, et al. The evolution of selfing in Arabidopsis thaliana. Science. 2007;317:1070–2.PubMedView ArticleGoogle Scholar
- Guo Y-L, Bechsgaard JS, Slotte T, Neuffer B, Lascoux M, Weigel D, et al. Recent speciation of Capsella rubella from Capsella grandiflora associated with loss of self-incompatibility and an extreme bottleneck. Proc Natl Acad Sci U S A. 2009;106:5246–51.PubMed CentralPubMedView ArticleGoogle Scholar
- Pettengill JB, Moeller DA. Tempo and mode of mating system evolution between incipient Clarkia species. Evolution. 2012;66:1210–25.PubMedView ArticleGoogle Scholar
- Bechsgaard JS, Castric V, Charlesworth D, Vekemans X, Schierup MH. The transition to self‐compatibility in Arabidopsis thaliana and evolution within S‐haplotypes over 10 Myr. Mol Biol Evol. 2006;23:1741–50.PubMedView ArticleGoogle Scholar
- Shimizu KK, Kudoh H, Kobayashi MJ. Plant sexual reproduction during climate change: gene function in natura studied by ecological and evolutionary systems biology. Ann Bot. 2011;108:777–87.PubMed CentralPubMedView ArticleGoogle Scholar
- Tsuchimatsu T, Suwabe K, Shimizu-Inatsugi R, Isokawa S, Pavlidis P, Städler T, et al. Evolution of self-compatibility in Arabidopsis by a mutation in the male specificity gene. Nature. 2010;464:1342–6.PubMedView ArticleGoogle Scholar
- Ornduff R. Reproductive biology in relation to systematics. Taxon. 1969;18:121–33.View ArticleGoogle Scholar
- Sicard A, Lenhard M. The selfing syndrome: a model for studying the genetic and evolutionary basis of morphological adaptation in plants. Ann Bot. 2011;107:1433–43.PubMed CentralPubMedView ArticleGoogle Scholar
- de Vos JM, Wüest RO, Conti E. Small and ugly? Phylogenetic analyses of the ‘selfing syndrome’ reveal complex evolutionary fates of monomorphic primrose flowers. Evolution. 2014;68:1042–57.PubMedView ArticleGoogle Scholar
- Eckert CG, Kalisz S, Geber MA, Sargent R, Elle E, Cheptou P-O, et al. Plant mating systems in a changing world. Trends Ecol Evol. 2010;25:35–43.PubMedView ArticleGoogle Scholar
- Schueller SK. Self-pollination in island and mainland populations of the introduced hummingbird-pollinated plant, Nicotiana glauca (Solanaceae). Am J Bot. 2004;91:672–81.PubMedView ArticleGoogle Scholar
- Micheneau C, Fournel J, Gauvin-Bialecki A, Pailler T. Auto-pollination in a long-spurred endemic orchid (Jumellea stenophylla) on Reunion Island (Mascarene Archipelago Indian Ocean). Plant Syst Evol. 2008;272:11–22.View ArticleGoogle Scholar
- Burney DA. Change and fire ecology as factors in the Quaternary biogeography of Madagascar. In: Lourenço WR, editor. Biogéographie de Madagascar. Paris: Editions de l’ORSTOM; 1996. p. 49–58.Google Scholar
- Gasse F, Van Campo E. Late Quaternary environmental changes from a pollen and diatom record in the southern tropics (Lake Tritrivakely Madagascar). Palaeogeogr Palaeoclimatol Palaeoecol. 2001;167:287–308.View ArticleGoogle Scholar
- Strijk JS, Noyes RD, Strasberg D, Cruaud C, Gavory F, Chase MW, et al. In and out of Madagascar: dispersal to peripheral islands insular speciation and diversification of Indian Ocean daisy trees (Psiadia Asteraceae). PLoS One. 2012;7, e42932.PubMed CentralPubMedView ArticleGoogle Scholar
- Burney DA, Burney LP, Godfrey LR, Jungers WL, Goodman SM, Wright HT, et al. A chronology for late prehistoric Madagascar. J Hum Evol. 2004;47:25–63.PubMedView ArticleGoogle Scholar
- Rakotoarisoa JE, Raheriarisena M, Goodman SM. Late Quaternary climatic vegetational shifts in an ecological transition zone of northern Madagascar: insights from genetic analyses of two endemic rodent species. J Evol Biol. 2013;26:1019–34.PubMedView ArticleGoogle Scholar
- Jarosz L. Defining and explaining tropical deforestation: shifting cultivation and population growth in colonial Madagascar (1896-1940). Econ Geogr. 1993;69:366–79.PubMedView ArticleGoogle Scholar
- Harper GJ, Steininger MK, Tucker CJ, Juhn D, Hawkins F. Fifty years of deforestation and forest fragmentation in Madagascar. Environ Conserv. 2007;34:1–9.View ArticleGoogle Scholar
- Cable S. New directions and challenges for the conservation of the flora of Madagascar. In: Bramwell D, Caujapé-Castells J, editors. The biology of Island Floras. Cambridge: Univ Press Cambridge; 2011. p. 425–42.View ArticleGoogle Scholar
- Catling PM, Bennett BA. Discovery of a possibly relict outbreeding morphotype of Sparrow’s-egg Lady’s-slipper Orchid, Cypripedium passerinum, in southwestern Yukon. Can Field Nat. 2007;121:295–8.Google Scholar
- Barrett SCH, Arunkumar A, Wright SI. The demography and population genomics of evolutionary transitions to self-fertilization in plants. Phil Trans R Soc B. 2014;369:20130344.PubMed CentralPubMedView ArticleGoogle Scholar