Rock outcrop orchids reveal the genetic connectivity and diversity of inselbergs of northeastern Brazil

Background Because of their fragmented nature, inselberg species are interesting biological models for studying the genetic consequences of disjoint populations. Inselbergs are commonly compared with oceanic islands, as most of them display a marked ecological isolation from the surrounding area. The isolation of these rock outcrops is reflected in the high number of recorded endemic species and the strong floristic differences between individual inselbergs and adjacent habitats. We examined the genetic connectivity of orchids Epidendrum cinnabarinum and E. secundum adapted to Neotropical inselbergs of northeastern Brazil. Our goals were to identify major genetic divergences or disjunctions across the range of the species and to investigate potential demographic and evolutionary mechanisms leading to lineage divergence in Neotropical mountain ecosystems. Results Based on plastid markers, high genetic differentiation was found for E. cinnabarinum (FST = 0.644) and E. secundum (FST = 0.636). Haplotypes were not geographically structured in either taxon, suggesting that restricted gene flow and genetic drift may be significant factors influencing the diversification of these inselberg populations. Moreover, strong differentiation was found between populations over short spatial scales, indicating substantial periods of isolation among populations. For E. secundum, nuclear markers indicated higher gene flow by pollen than by seeds. Conclusions The comparative approach adopted in this study contributed to the elucidation of patterns in both species. Our results confirm the ancient and highly isolated nature of inselberg populations. Both species showed similar patterns of genetic diversity and structure, highlighting the importance of seed-restricted gene flow and genetic drift as drivers of plant diversification in terrestrial islands such as inselbergs.


Background
Inselbergs are isolated rock outcrops typically harboring rupicolous vegetation and embedded within a landscape composed of contrasting plant communities. Because of their disconnected geographic nature, inselbergs are frequently compared with oceanic islands, since most of them display a marked ecological isolation from the surrounding area [1]. The isolation and ancient age of these rock outcrops is reflected in the high number of recorded endemic species (reviewed by [2]) and the strong floristic differences observed between individual inselbergs and their surrounding habitats [3,4]. Some authors have also suggested that inselbergs may have acted as refugia for xerophilic or cold-adapted species during glacial/interglacial cycles [4,5]. In addition, most inselbergs maintain their typical attributes irrespective of geographic location, enabling broad-scale comparisons even between different and contrasting biomes, such as mesic forests and seasonally dry plant communities [2]. As biological islands, inselbergs are promising ecosystems for biogeographic and evolutionary studies, comparable to their oceanic counterparts [6].
Phylogeographic studies of inselberg-adapted species have provided insights into evolutionary processes leading to diversification of lineages and species [7][8][9][10][11]. Such studies have uncovered strong phylogeographic structure, high population differentiation and extensive genetic diversity levels, supporting the view of inselbergs as centers of species diversity and endemism. Species occurring in small isolated patches are expected to experience reduced gene flow, significant genetic drift and high levels of population divergence [6]. Indeed, genetic studies of inselberg species have provided strong support for these expectations [7,[12][13][14][15]. In isolated populations, gene flow is constrained, and genetic drift is expected to be the predominant force governing allele frequencies [8,16].
Porembsky et al. [17] have noted that plant community organization on rock outcrops is driven by stochastic colonization events. Thus, phylogeographic studies of multiple co-occurring species with similar biological traits (e.g. pollination and seed dispersion) may be particularly useful for understanding the role of stochasticity (i.e. genetic drift) in the evolution of these naturally fragmented populations. By comparing the phylogeographic structures of different species, one can infer whether the current plant community has been historically stable as evidenced by geographically similar genetic structure. Alternatively, if species distributions are ephemeral over evolutionary time and more influenced by intrinsic species preferences, a mixture of phylogeographic structures is expected [18,19].
Phylogeographic studies focused on organisms associated with mesic forest communities have indicated a strong association between the occurrence of glacial cycles and the fragmentation of forest-dwelling species, supporting the classical model of tropical refugia [20][21][22]. On the other hand, species associated with open biomes, such as grassland, savannas and dry forests, show variable demographic responses to past glacial cycles, suggesting a more complex scenario [21,23]. In this regard, cross validation of environmental envelope models (EEMs according to terminology suggested by [24]) and molecular genetic data have confirmed different past demographic scenarios for a wide array of species. Fragmentation [25], expansion [26] and long-term persistence [27] have been detected in organisms associated with open vegetation communities. By combining EEMs and molecular genetic tools, different phylogeographic hypotheses depicting the role of past climatic fluctuations in lineage diversification can be tested.
The Caatinga is the largest, most diverse dry seasonal tropical forest biome in the Neotropical region [28]. Despite the extensive geographic area of this biome, relatively few phylogeographic studies have been conducted using species from the Caatinga [27,[29][30][31][32]. This paucity of studies precludes broad conclusions regarding the impact of past climate oscillation events on the genetic structure of organisms. The high number of endemic species [28] coupled with old lineage-divergence times [27,33,34] implies an ancient origin for Caatinga communities. In addition, paleoclimate models [35] indicate that dry forest distribution within the Caatinga was stable during the Last Glacial Maximum (LGM). Indeed, high species diversity and floristic differences have been found among Caatinga rock outcrops [3,36], supporting the hypothesis of inselbergs as ancient, stable refuges of diversity [4,5].
In the present study, we used nuclear and plastid microsatellite markers to analyze phylogeographic structure and genetic diversity of Epidendrum cinnabarinum and E. secundum (Orchidaceae), two widespread, cooccurring species commonly found on inselbergs in northeastern Brazil. Epidendrum cinnabarinum, a polyploid species with 2n = 240 chromosomes [37], occurs in Caatinga and Brazilian Atlantic Forest (BAF) inselbergs and on sand dune vegetation along the seashore in northeastern Brazil. Epidendrum secundum exhibits a diploid chromosome count of 2n = 56 in most populations [37]; it has a much broader geographic distribution, occurring in Central America, Guiana and Andean ranges and the Brazilian Central Shield [38]. This species occurs preferentially on rock outcrops, with abundant populations observed on Caatinga and BAF inselbergs. In this study, we also applied EEMs to explore the demography of both species during important climatic oscillation events: the LGM (21,000 ka) and Last Interglacial (LIG; 120,000 ka). The following specific questions were addressed: (1) What is the current extent of genetic structure and the degree of isolation among disjoint populations on inselbergs? (2) Were the current range distributions of both species stable during Late Quaternary climatic oscillations? and (3) Are there phylogeographic breaks separating populations from different biomes (Caatinga and BAF) and ecoregions (Chapada Diamantina and Planalto da Borborema)? We also considered the phylogeographical and genetic structure of E. cinnabarinum and E. secundum in light of results based on climate modeling, paleovegetation reconstructions and the island-like nature of populations distributed on inselbergs.

Plastid and nuclear genetic diversity of sampled populations
Analysis of seven plastid loci recovered a total of 10 haplotypes for E. cinnabarinum and 12 haplotypes for E. secundum ( Figure 1C, 1D, Table 1, Additional file 1:  Table S1), with no haplotype sharing between species. One to five haplotypes were found within E. cinnabarinum populations, while one to six haplotypes were detected within E. secundum populations. Similar diversity levels were detected for both species concerning haplotype richness (E. cinnabarinum 0.000-3.066; E. secundum 0.000-4.873) and haplotype diversity (E. cinnabarinum 0.000-0.813; E. secundum 0.000-0.800) ( Table 1). In two sympatric populations, TO and BZ, genetic diversity parameters were identical for both species (Table 1).
With respect to six nuclear markers genotyped in E. secundum, moderate levels of genetic diversity were observed in most populations ( Table 2). The number of alleles per population ranged from 19 to 23, and allelic richness per population varied from 2.28 to 2.77. Expected heterozygosity per population ranged from 0.442 to 0.608. One to nine private alleles were observed within populations, with values of private allelic richness ranging from 0.11 to 0.72. Inbreeding coefficients, which ranged from −0.032 to 0.156, were not significantly different from zero in any population. The presence of genotyping errors due to stuttering or null alleles (frequencies ranging from 0 to 0.09, estimated using Brookfield Equation one) was ruled out using MICRO-CHECKER tests [39].
Genetic structure of E. cinnabarinum A haplotype network based on plastid markers showed loops (ambiguities), and no apparent geographic structure was observed ( Figure 1C). The most frequent haplotype in both BAF and Caatinga biomes (Hc1) was found in 46.8% of individuals and in 5 out of 11 sampled populations. Hc5, the second most frequent haplotype, was also distributed in both biomes and occupied a central position in the haplotype network ( Figure 1C). Only 3 (Hc1, Hc5 and Hc6) out of 10 haplotypes were found to occur in more than one population. High genetic differentiation was observed across populations, with F ST = 0.644 and G ST = 0.672. According to an analysis of molecular variance (AMOVA) ( Table 3), a high proportion of the genetic variability in the haplotype data resided among populations (64.41%, P < 0.0001), with only 35.59% accounted for within populations. Hierarchical AMOVA did not support a division between Caatinga and BAF populations (Table 3; P = 0.704). Pairwise F ST comparisons among populations (Additional file 2: Table S2) ranged from 1.000 to 0.000, with most values found to be significant (P < 0.005). Haplotype differentiation between populations separated over short spatial scales was observed between populations QE and TO (12 km apart) and between BZ and FE (18 km apart) (Additional file 2: Table S2, Figure 1A). Other than the strong genetic differentiation found among populations, no among-population phylogeographic structure or isolation by distance was detected (P = 0.114 and P = 0.084, respectively).

Genetic structure of E. secundum
Analysis of plastid markers yielded a network lacking clear geographic structuring ( Figure 1D). Haplotype Hs5 was shared between populations from Caatinga and BAF biomes, and also between Chapada Diamantina and Planalto da Borborema ecoregions ( Figure 1B, 1D). Hs5 was also the most frequent haplotype, found in 54.5% of individuals and in six out of nine sampled populations. Seven haplotypes were restricted to Chapada Diamantina populations ( Figure 1B).
The only hypothesis of population differentiation significantly supported by AMOVA was that between populations from Chapada Diamantina and the remaining localities ( Table 4). The partitioning of plastid genetic diversity and structure among E. secundum populations was also very similar to the patterns observed for E. cinnabarinum. AMOVA results indicated that a high proportion of uncovered genetic variability was found among populations (63.64%, P < 0.0001), with only 36.36% attributed to within-population variability (Table 4). Genetic differentiation was high across all populations, with F ST = 0.636 and G ST = 0.632. As observed for E. cinnabarinum, contrasting haplotypes were detected between populations BZ and FE. Pairwise F ST comparisons among populations ranged from 1.000 to 0.000, with  most values significant (P < 0.005; Additional file 3: Table S3). No phylogeographic structure was detected among populations, as R ST was not significantly larger than F ST (P = 0.503). Isolation by distance among populations was also not detected (P = 0.822). Compared with plastid markers, nuclear markers exhibited lower levels of genetic differentiation as estimated by Φ ST (0.114), G ST (0.122) and G' ST (0.255). An AMOVA also generated different results than for plastid markers, with most of the genetic variation partitioned within (95.42%, P < 0.0001) rather than among (4.58%, P = 0.087) populations (Table 4). Differentiation between BAF and Caatinga populations and between Chapada Diamantina and Planalto da Borborema was not significant (Table 4). No sign of phylogeographic structure was detected among populations (P = 0.329). A Mantel test did not support isolation by distance among populations based on Φ ST , G ST or G′ ST (all P > 0.05).
As shown in Additional file 4: Figure S1, simulations performed in STRUCTURE consistently identified K = 2 clusters. Most analyzed individuals showed admixed ancestry, however, and genetic subdivisions among populations were unclear ( Figures 1B and 2). Admixture proportions (Q) for population BZ showed strong assignment to cluster 2 (average assignment proportion = 0.92; Figure 2), whereas specimens from the southernmost population MC showed assignment proportions (0.77; Figure 2) associated with cluster 1. Intermediate assignment proportions were observed for most individuals from remaining populations ( Figure 2).    Table 2). The lowest observed critical value (M C ) obtained from simulations performed with different values of parameter θ, proportion of single-step mutations (pg) and average size of non-one-step mutations (Δg) was 0.589 (Table 5). According to this Mc threshold, none of the populations showed signatures of past genetic bottlenecks (Table 2). Using values of genetic differentiation among populations BZ, DE, TO, FE and MC obtained for nuclear (G ST = 0.122) and plastid (G ST = 0.782) markers, the ratio of pollen flow to seed flow was estimated as 25.85. The value of this ratio suggests that gene flow via pollen in E. secundum is > 20 times higher than that occurring via seeds.

Potential ancient distributions of E. cinnabarinum and E. secundum
Distributions under LIB, LGM and current climatic conditions were successfully estimated for E. cinnabarinum and E. secundum. As indicated by high area-under-thecurve (AUC) values (0.969 for E. secundum and 0.968 for E. cinnabarinum), the EEM analysis performed well. Climate variables contributed to both species models, with the E. cinnabarinum model additionally adjusted using a geological eras variable. Past distribution ranges differed between the two species, mainly at their peripheral limits. From the LIG to the present, E. cinnabarinum exhibited a pronounced decrease in its presence at Chapada Diamantina and other inland mountain ranges, with a more recent expansion towards northern Planalto da Borborema and coastal regions after the LGM ( Figure 3A, 3B and 3C). Populations distributed in the southern portion of the northeastern Brazilian seashore did not show marked demographic oscillations. A decrease in the inland distribution of E. secundum was also observed, but without any further expansion ( Figure 3D Table 1 for population details.  E. secundum has a fragmented distribution within Planalto da Borborema, where most populations are found at the border between Caatinga and BAF biomes. Since the LIG, the distribution of this species has not been continuous between Chapada Diamantina and Planalto da Borborema.

Discussion
Genetic studies of inselberg species can provide insights into the combined effects of genetic drift and restricted gene flow on evolution and diversification of lineages restricted to disjunct populations [1,2,5]. Phylogeographic studies of inselberg species have indeed confirmed many  Table 1 for population details.
In some extreme cases of isolation, within-inselberg interspecific gene flow between congeneric species has been found to be higher than intraspecific gene flow between different inselbergs [7,14]. To understand the effects of long-term disjoint distribution on genetic architecture and evolution of inselberg populations, we carried out a phylogeographic study of two rock-outcrop orchid species showing very similar life traits. The use of plastid markers revealed marked genetic differentiation among populations. In contrast, nuclear genetic differentiation was much lower in E. secundum, highlighting the important role of pollen dispersal, and consequently pollination services, in species cohesion. This result has also been observed in other plants [7,40]. Results from EEM analysis of E. secundum suggest long-term stability of populations from Chapada Diamantina and the southern portion of Planalto da Borborema ( Figure 3D, 3E and 3F); in contrast, results for E. cinnabarinum imply a decrease in distribution within Chapada Diamantina, in agreement with its current rarity in this ecoregion. After the LGM, E. cinnabarinum expanded its distribution to northern sand dune vegetation communities. This expansion most likely occurred from inland inselberg populations, not from southern coastal ones, as extensive haplotype sharing was found among seashore and inland populations. The joint use of phylogeography and EEMs helped to clarify the origin of genetic diversity and population differentiation in these inselberg species, emphasizing the importance of such studies in terrestrial island-like environments.

Marked genetic differentiation among inselberg populations
Genetic differentiation among populations of both E. cinnabarinum and E. secundum was not significantly associated with BAF and Caatinga biomes (Tables 3  and 4). In contrast, extensive haplotype sharing was found among populations distributed in both biomes ( Figure 1A and 1B). Haplotype sharing was also found between populations from inselbergs and sand dune vegetation, in agreement with floristic similarities reported in previous studies [41]. Because of their xeric nature, inselbergs can provide suitable conditions for drought tolerant species even within humid forest patches [2]. Inselbergs embedded in mesic forests often harbor xeric plant communities and can be considered refugia for species adapted to drier climatic conditions, highlighting their importance for comparative studies across different biomes [2,4]. Paleontological evidence [42] indicates that alternating dry and wet periods have influenced the BAF distribution of plant species, particularly along the northeastern coast. Inselberg populations within the BAF (populations JA, CB, TB and JI) may thus represent relicts of xeric plant communities that probably extended further towards the coast during drier periods, thereby allowing gene exchange with populations currently placed in the Caatinga biome. The most outstanding pattern recovered by the use of plastid markers was the deep among-population differentiation within E. cinnabarinum and within E. secundum ( Figure 1A and 1B, Additional files 2 and 3: Table S2 and S3). Even populations separated by a few kilometers exhibited strong differentiation based on the plastidial genome (populations TO and QE, 12 km apart, F ST = 0.985; BZ and FE, 18 km apart, F ST = 1.000; ES and RR, 34 km apart, F ST = 0.727). The significant pairwise genetic differentiation among most populations, the absence of a pattern of isolation by distance, and the lack of significantly phylogeographically structured haplotypes suggest that genetic drift, as observed in other rock outcrop species [7][8][9]43], has been an important force shaping plastidial genetic variation in E. cinnabarinum and E. secundum. In finite populations, such as those inhabiting inselbergs and oceanic islands, genetic drift may be the predominant force governing allele frequencies of neutral loci [6,44]. Because effective population size of maternally inherited organelles is expected to be an order of magnitude lower than that of nuclear genes, genetic drift is expected to play a stronger role in shaping frequencies of organellar alleles than those of nuclear alleles [45]. In addition to the spatial dimension, genetic drift includes an intrapopulation temporal component in which a series of subpopulations diverge from one another over time. Fine-scale genetic analysis [43,46] coupled with long-term studies (e.g. [8]) may help to clarify the impact of drift on the genetic structure of inselberg populations.
The genetic structure observed among E. secundum populations was higher for plastid than for nuclear markers ( Figure 2, Table 4). Most haplotype diversity (63.64%) was found to reside among populations; this result is in contrast to nuclear markers, in which diversity among populations was not significant (4.58%). This difference in genetic structure between genomic compartments was translated into a pollen/seed flow ratio of 25.85, indicating that pollen-mediated gene dispersal is 25-fold more efficient than dispersal through seeds. The inference of extensive pollen flow between E. secundum populations is also supported by observations of low inbreeding coefficients (Table 2) and a high proportion of genetic diversity within populations (Table 4). In addition, Bayesian assignment based on nuclear markers indicated intermediate admixture proportions for many individuals and specimens assigned to different nuclear genetic clusters in the same population ( Figure 1B); these results also suggest higher levels of gene flow by pollen. Efficient gene flow by pollen is normally expected for food-deceptive orchids [40,47], as pollinators avoid plants in the same patch, thereby promoting pollen flow over long distances and reducing the chances of geitonogamous pollination. In general, low levels of nuclear genetic structure are observed in orchids (reviewed by [48]). Pollen flow has proven to be crucial to the maintenance of species cohesion in other fragmented inselberg species, highlighting the importance of pollinators in promoting gene exchange among rock outcrops [7,8,43].
Despite the strong genetic differentiation observed between populations from Chapada Diamantina and the remaining localities (Table 4), phylogeographic structure was not significant in E. secundum. This surprising situation is probably due to ancestral haplotype sharing (Hs5) among most populations, indicating ancient dispersal events. Another explanation for the non-significance of phylogeographic structure is the fact that multiple, phylogenetically unrelated haplotypes are found in the same populations, especially in the Chapada Diamantina ecoregion. Chapada Diamantina is an old and very stable region; demographic changes in peripheral populations, such as those inhabiting the BAF and from mountains south of Chapada Diamantina, may have contributed to the exchange of unrelated haplotypes-a pattern observed in other studies [32,34,49]. Indeed, species distribution modeling shows a broader past geographic distribution for E. secundum, with a former range that includes mountains to the south and north of Chapada Diamantina.

Demographic changes in E. secundum and E. cinnabarinum
For both species, results obtained using EEMs have shed light on whether haplotype sharing among populations is due to regular gene flow, past long-distance dispersal events or past continuous distribution. For E. cinnabarinum, a general pattern of expansion from the inland towards the seashore since the LIG is suggested ( Figure 3A, 3B and 3C). This species appears to have been more commonly distributed in Chapada Diamantina during the LIG, showing a further decrease in distribution until the present. Epidendrum cinnabarinum is in fact currently very rare in Chapada Diamantina, where its populations are composed of only a few individuals. Colonization of the northern portion of the seashore occurred during the last 20,000 years, probably from Planalto da Borborema populations. The haplotype sharing detected among populations and the low diversity of Mamanguape and Cabo de Santo Agostinho populations are consistent with this hypothesis. Because only ancestral haplotypes are shared among populations, recurrent gene flow is unlikely; in this case, long-distance dispersal may have been responsible for the northern seashore colonization. Results from Behling [42] suggest that xeric vegetation replaced mesic plant communities along the northern portion of the northeastern Brazilian seashore in the last 10,000 years, a trend likely favoring the expansion of E. cinnabarinum from inland inselbergs towards dry coastal vegetation zones. Haplotype sharing between inland and coastal populations and the inferred expansion of inland populations towards the seashore have also been observed in other phylogeographic studies along the Brazilian coast, confirming the connection between the BAF and other adjacent biomes such as the Caatinga [32,49], Cerrado [50,51] and Pampas [40,52].
Analysis of EEMs recovered different demographic patterns for E. secundum ( Figure 3D, 3E, 3F). The longterm persistence of Chapada Diamantina populations is in agreement with the high levels of genetic diversity observed in Morro do Chapéu and Palmeiras populations (Tables 1 and 2) and the absence of bottlenecks ( Table 5). The genetic differentiation of Chapada Diamantina and Planalto da Borborema populations was also confirmed by the EEM analysis, which demonstrated a clear discontinuity in E. secundum distribution along these two ecoregions in the last 120,000 years. High levels of plant endemism and biodiversity have been recorded for Chapada Diamantina [53,54], and floristic results indicate different species compositions between rock outcrops from Chapada Diamantina and Planalto da Borborema [3]. In addition, phylogeographic studies have indicated the occurrence of older lineages [32] and genetic differentiation among populations [27,34] occurring within Chapada Diamantina. The patterns of genetic diversity found within Chapada Diamantina populations are in agreement with refuge theory, which predicts higher levels of genetic diversity and older lineages to be present in historically climate-stable regions [5,22].
The relictual nature of inselberg populations is also evident in the structure of the haplotype networks ( Figure 1C and 1D). Several mutational steps were observed in the networks calculated for E. cinnabarinum and E. secundum, suggesting substantial periods of isolation among populations of both species. Deep divisions among haplotypes may be reflective of populations with long divergence histories [16], as demonstrated by phylogeographic studies of Caatingan vertebrate animals having fragmented distributions [27,29,31,32]. Divergence times estimated for Caatingan organisms indicate major crown clade splits to have occurred in the Miocene (between 11,000 ka and 5,000 ka), followed by divergence within these major groups during the Pliocene-Pleistocene [20,27,33,34]. These results are also in agreement with the long-term persistence of populations inhabiting inselbergs [1]. Accordingly, lower levels of genetic differentiation have been observed in other Epidendrum species inhabiting well-connected landscapes (e.g. sand dunes and swamp vegetation) [40,55].

Conclusions
The population genetic structure recovered for both E. cinnabarinum and E. secundum reflects the isolated nature of their inselberg habitats. Similar to the results of studies of other inselberg species [7,9,11,14,15], significant differentiation at plastid loci was observed between populations separated by only a few kilometers (Additional files 2 and 3: Table S2 and S3). Genetic differentiation among populations was not correlated with geographic distance or ecoregions, suggesting that genetic drift may be a significant factor influencing the diversification of inselberg species. According to Porembski et al. [17], inselberg species composition is difficult to predict even over short distances, probably because of stochastic colonization events. Long-term persistence of populations was also supported in our study, suggesting that both species have been restricted to the specific inselberg habitat over long time frames, and have not been reduced to rock outcrops because of climatic instability. The results of EEM analysis indicate that fragmentation and reduction of E. secundum distribution has occurred in mountain ranges south of Chapada Diamantina. Future studies should consequently include a broad sampling of this species over adjacent mountain ranges to explore how such distribution oscillations have affected current levels of genetic diversity in these populations. In northeastern Brazil, inselbergs have different mineral origins and are distributed across different biomes and climatic clines. Thus, selection for divergent habitats may also be of great importance during lineage diversification within species under restricted gene flow. Given the contrasting patterns of genetic structure observed between nuclear and plastid markers, cytonuclear incompatibilities may contribute to the first stages of reproductive isolation among divergent lineages. Future efforts should take advantage of next-generation sequencing approaches [56] coupled with experimental studies [57,58] to clarify the role of drift, selection and cytonuclear interactions in the diversification of inselberg species.

Study species
Epidendrum L. is the largest genus of Orchidaceae in the Neotropics, with about 1,500 species distributed from the southeastern United States to northern Argentina [59]. The genus contains many species with wide distributions and high morphological diversity [38]. Epidendrum cinnabarinum and E. secundum are perennial, terrestrial species commonly found on inselbergs whose locations include sea-level to mountain elevations. Both species belong to subgenus Amphyglottium [60,61], a group of food-deceptive orchids whose flowers are visited by several butterfly species despite the lack of nectar reward (reviewed by [38]). Although the species are self-compatible, pollinators are necessary for seed set (Pinheiro unpub. res., [62]).

Sampling design
Samples of E. cinnabarinum and E. secundum were obtained from 11 and 9 locations, respectively ( Table 6, Figure 1A and B). Sample sites covered most of the geographic distribution of these species in the Caatinga, with additional samples collected in the BAF, an adjacent biome with higher levels of humidity and species diversity (perennial mesic forest). To estimate the degree of isolation of inselberg populations, samples were collected from sites ranging from 12 to 1,000 km apart. Most populations were sampled on rock outcrops; E. cinnabarinum was also collected in sand dune vegetation (populations MA, PI and AB, Table 6). Both species occurred in sympatry in populations TO, BZ and FE. 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. [63].

Genetic diversity of sampled populations
All sampled populations were characterized for levels of diversity based on plastid DNA markers. The number of haplotypes in each population, haplotype diversity and haplotype richness were estimated using RAREFAC v. 3.5 software [67]. Estimates of haplotype richness were corrected for differences in sample size using the rarefaction method.
Nuclear microsatellite diversity of five E. secundum populations was characterized according to number of alleles, number of private alleles, allelic richness, private allelic richness, expected heterozygosity and inbreeding coefficient [68], which were calculated using the programs MSA v. 4.05 [69] and HP-RARE v. 1.0 [70]. Departures from Hardy-Weinberg Equilibrium of withinpopulation inbreeding coefficients were identified using exact tests in GENEPOP v. 4.0 [71]. The microsatellite data set was tested for genotyping errors due to stuttering, short allele dominance and null alleles by means of a Monte Carlo simulation of expected allele-size differences as implemented in MICRO-CHECKER [39].

Plastid genetic structuring
The geographical structure of genetic variation in plastid DNA was investigated for both species through several approaches. A median-joining network [72] based on plastid DNA was constructed using the program NET-WORK v. 4.5.1.0 (www.fluxus-engineering.com). To assess whether observed genetic differentiation was due to drift, we tested the hypothesis that G ST = R ST [73] (where Slatkin's R ST [74] estimates the contribution of stepwiselike mutations to genetic differentiation) following Pons and Petit [75] using the program PERMUT/CpSSR. To test for the presence of phylogeographic structure among populations (i.e. R ST significantly larger than G ST ), 10,000 permutations of R ST values were performed. Partitioning of genetic diversity within and among populations and between populations from BAF (coastal) and Caatinga (inland) regions was assessed by AMOVA. For E. secundum, two additional models were considered, one corresponding to differentiation between populations from Chapada Diamantina and Planalto da Borborema, and the other considering populations from Chapada Diamantina separated from remaining localities. Genetic differentiation between populations was also measured by pairwise comparisons of F ST using the program ARLEQUIN v. 3.5. Finally, the hypothesis that populations Populations sampled with their identification code (ID), elevation above sea level, habitat description and sample size, analysed for nuclear and plastid markers in Epidendrum cinnabarinum (C) and E. secundum (S) from Northeast Brazil, distributed in the Caatinga (CAA) and Brazilian Atlantic Forest (BAF) biomes. Populations are indicated as shown on the maps in Figure 1A and B.
were differentiated because of isolation by distance was tested by assessing the correlation between pairwise geographic distances and pairwise F ST values using a Mantel test in the program GENEPOP v. 4.0. Correlation significance was estimated after performing 10,000 permutations between pairwise geographic distance and pairwise genetic differentiation matrices.

Nuclear genetic structuring
Nuclear genetic differentiation among populations of E. secundum was measured using the unbiased estimator of relative differentiation G ST [76] and the standardized genetic differentiation measure G′ ST [77]. Partitioning of genetic diversity was examined at different hierarchical levels using AMOVA as implemented in ARLEQUIN v. 3.5. The hypothesis that populations were differentiated because of isolation by distance was tested by assessing the correlation between pairwise geographic distance and pairwise G′ ST values using a Mantel test in GENEPOP v. 4.0, similar to the approach used for plastid markers. In addition, to test for the presence of phylogeographic structure among populations, R ST was calculated by permuting allele sizes among alleles for 10,000 permutations using the program SPAGEDI v. 1.3d [78], and was then compared with F ST values. Phylogeographic structure was inferred when R ST was significantly larger than F ST . Bayesian assignment analysis (in STRUCTURE v. 2.3.3 [79]) was used to assign individuals to genetic clusters (K) and to estimate admixture proportions (Q) for each individual. A set of models was chosen in which individuals had admixed ancestries and correlated allele frequencies. Ten replicate runs were completed for each K value for K = 1-10. Markov chain Monte Carlo runs consisted of 600,000 generations after an initial burn-in of 250,000 generations. The most probable number of genetic clusters (K) present in the data were defined following Evanno et al. [80] using the program Structure Harvester v. 6.0 [81].

Demographic analyses
Recent population size reductions (i.e. genetic bottlenecks) were tested in E. secundum with nuclear microsatellitebased M-ratios [82] using ARLEQUIN v. 3.5. Significance for each population was assessed by comparisons of M-ratios and critical values (M C values) obtained by simulating the distribution of M-ratios under specific demographic and mutational conditions using the software CRITICAL_M.EXE (http://swfsc.noaa.gov/textblock.aspx? Division=FED&id=3298). The critical value M C is set at the lower 5% tail of this distribution, and bottlenecks are detected when the M-ratio value is below the calculated M C threshold. Different values of M C were simulated by modifying θ (0.5, 2.0 and 10.0), pg (0.1 and 0.3) and Δg (2.0 and 3.5). Low M C values are more conservative as a bottleneck must be of greater intensity to drop below this level [82]; for this reason, the lowest obtained M C value was used to check bottleneck significance.
To examine the relative contribution of pollen vs. seed flow to total gene flow among E. secundum populations, the parameter G ST for nuclear and plastid loci was compared based on populations genotyped for both plastid and nuclear markers (TO, DE, BZ, FE and MC). Pollen/ seed flow ratio was estimated following Ennos [83] and Petit et al. [84] using Equation 1 presented by the latter authors.
Environmental envelope models for E. cinnabarinum and E. secundum To test the hypothesis that inselbergs were refugia through late Quaternary climate cycles, we projected EEMs from the current situation onto late Quaternary scenarios. We used the maximum entropy algorithm of Maxent v. 3.3.3e [85,86] to obtain current and past distributions. We used current and past bioclimatic variables derived from monthly temperature and rainfall (Bioclim scheme). Current data were obtained from WorldClim [87]. Two past scenarios were considered: LIG (~120,000-140,000 years BP), based on Otto-Bliesner et al. [88], and LGM (~21,000 years BP), based on Paleoclimate Modelling Intercomparison Project Phase II [89] considering the CCSM3 model [90]. To incorporate geological differentiation between outcrops, we also used three variables derived from Schenk et al. [91]. Geological variables were major geological types (sedimentary/igneous and metamorphic rocks), eras and geological age. Major geological types and eras were categorical and age was continuous. All variables were adjusted to a spatial resolution of 30 arcsec. Highly correlated variables were removed, and jackknifing was used to estimate variable importance. The final models were obtained considering only variables with contributions to AUC higher than 75%. Five model replicates were run for each one of the presence-only methods, with 75% of occurrences used for calibration and different subsets (25%) used for validation. Mean AUC was used to assess the performance of the models [92], where 1 was the maximum prediction and 0.5 suggested a random prediction. Past distributions were estimated by projecting the current relationship onto scenarios of past climate, assuming that current relationships between climate and distribution were maintained.
We obtained 226 E. cinnabarinum and 242 E. secundum occurrences from our field GPS records and georeferenced herbarium data extracted from the speciesLink project (http://splink.cria.org.br). All points were verified with Google Earth to ensure localities were not placed in heavily urbanized areas. The localities used for EEM analysis are available from the authors upon request to avoid illegal plant collecting.