# Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis-scenario

Sebastián Escobar, Jean-christophe Pintaud, Henrik Balslev, Rodrigo Bernal, Mónica Moraes Ramírez, Betty Millán, Rommel Montúfar

► To cite this version:

Sebastián Escobar, Jean-christophe Pintaud, Henrik Balslev, Rodrigo Bernal, Mónica Moraes Ramírez, et al.. Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis-scenario. Ecology and Evolution, 2018, 8 (16), pp.8030-8042. 10.1002/ece3.4216. hal-03356574

HAL Id: hal-03356574

https://hal.umontpellier.fr/hal-03356574

Submitted on 28 Sep 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis-scenario

Article in Ecology and Evolution · July 2018

DOI: 10.1002/ece3.4216

CITATIONS

4

READS

350

7 authors, including:

<figure><img src="image_2.png" /></figure>

Sebastián Escobar
Fundación para la Conservación de los Andes Tropicales (FCAT)

8

PUBLICATIONS

6

CITATIONS

SEE PROFILE

<figure><img src="image_3.png" /></figure>

Rodrigo Bernal
Reserva Natural Guadualito

144

PUBLICATIONS

2,466

CITATIONS

SEE PROFILE

<figure><img src="image_4.png" /></figure>

Henrik Balslev
Aarhus University

360

PUBLICATIONS

11,703

CITATIONS

SEE PROFILE

<figure><img src="image_5.png" /></figure>

Mónica Moraes R.
Universidad Mayor de San Andres

121

PUBLICATIONS

1,035

CITATIONS

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

<figure><img src="image_6.png" /></figure>

Impact of palm harvest in northwest South America View project

<figure><img src="image_7.png" /></figure>

Enciclopedia de plantas útiles del Ecuador View project
# Genetic structuring in a Neotropical palm analyzed through an Andean orogenesis-scenario

Sebastián Escobar<fnref n="1" />,<fnref n="2" /> | Jean-Christophe Pintaud<fnref n="3" /> | Henrik Balslev<fnref n="2" /> | Rodrigo Bernal<fnref n="4" /> |
Mónica Moraes Ramírez<fnref n="5" /> | Betty Millán<fnref n="6" /> | Rommel Montúfar<fnref n="1" /> ID

## Correspondence

Rommel Montúfar, Facultad de Ciencias Exactas y Naturales, Pontificia Universidad Católica del Ecuador, Av. 12 de Octubre 1076 y Roca, Quito, Ecuador.

Email: rjmontufar@puce.edu.ec

## Funding information

Seventh Framework Programme, Grant/Award Number: 212631; Pontificia Universidad Católica del Ecuador; Det Frie Forskningsråd, Grant/Award Number: 4181-00158

## Abstract

Andean orogenesis has driven the development of very high plant diversity in the Neotropics through its impact on landscape evolution and climate. The analysis of the intraspecific patterns of genetic structure in plants would permit inferring the effects of Andean uplift on the evolution and diversification of Neotropical flora. In this study, using microsatellite markers and Bayesian clustering analyses, we report the presence of four genetic clusters for the palm *Oenocarpus bataua* var. *bataua* which are located within four biogeographic regions in northwestern South America: (a) Chocó rain forest, (b) Amotape-Huancabamba Zone, (c) northwestern Amazonian rain forest, and (d) southwestern Amazonian rain forest. We hypothesize that these clusters developed following three genetic diversification events mainly promoted by Andean orogenic events. Additionally, the distinct current climate dynamics among northwestern and southwestern Amazonia may maintain the genetic diversification detected in the western Amazon basin. Genetic exchange was identified between the clusters, including across the Andes region, discarding the possibility of any cluster to diversify as a distinct intraspecific variety. We identified a hot spot of genetic diversity in the northern Peruvian Amazon around the locality of Iquitos. We also detected a decrease in diversity with distance from this area in westward and southward direction within the Amazon basin and the eastern Andean foothills. Additionally, we confirmed the existence and divergence of *O. bataua* var. *bataua* from var. *oligo-carpus* in northern South America, possibly expanding the distributional range of the latter variety beyond eastern Venezuela, to the central and eastern Andean cordilleras of Colombia. Based on our results, we suggest that Andean orogenesis is the main driver of genetic structuring and diversification in *O. bataua* within northwestern South America.

## KEYWORDS

genetic divergence, genetic diversity, microsatellite markers, *Oenocarpus*, phylogeography

<footnote n="1">¹Facultad de Ciencias Exactas y Naturales,
Pontificia Universidad Católica del Ecuador,
Quito, Ecuador</footnote>

<footnote n="2">²Department of Bioscience, Ecoinformatics and Biodiversity Group, Aarhus University, Aarhus, Denmark</footnote>

<footnote n="3">³DYNADIV, UMR DIADE, Institut de Recherche pour le Développement, Montpellier, France</footnote>

<footnote n="4">⁴Instituto de Ciencias Naturales, Universidad Nacional de Colombia, Bogotá, Colombia</footnote>

<footnote n="5">⁵Instituto de Ecología, Universidad Mayor de San Andrés, La Paz, Bolivia</footnote>

<footnote n="6">⁶Museo de Historia Natural, Universidad Nacional Mayor de San Marcos (UNMSM), Lima, Perú</footnote>
# 1 | INTRODUCTION

Plant diversification at the regional level is strongly linked to tectonism and subsequent climate modification (Hoorn, Wesselingh, Ter Steege, et al., 2010). The orogeny of the Andes during the Cenozoic is considered one of the most important drivers of the very high plant diversity found in the Neotropics (Antonelli, Nylander, Persson, & Sanmartín, 2009; Gentry, 1982; Hoorn, Wesselingh, Ter Steege, et al., 2010; Luebert & Weigend, 2014; Turchetto-Zolet, Pinheiro, Salgueiro, & Palma-Silva, 2012). Andean uplift promoted diversification in northern South America by increasing habitat heterogeneity and geographic vicariance which may promote genetic isolation and ultimately speciation (Antonelli & Sanmartín, 2011; Mora et al., 2010). Particularly, cross-Andean vicariance (at both sides of the Andes) may have played an important role in Neotropical plant evolution as at least half of the Ecuadorian Flora could have appeared through vicariance after the Andean uplift (Balslev, 1988). Jointly, the uplift also affected regional climate by increasing rainfall along the eastern flanks of the Andes due to the disruption of large-scale atmospheric circulation (Garreaud, Vuille, Compagnucci, & Marengo, 2009; Luebert & Weigend, 2014), resulting in the modification of the continental hydrological network (Hoorn, Wesselingh, Ter Steege, et al., 2010). As plants are responsive to past geological changes (Klimova, Hoffman, Gutierrez-Rivera, Leon de la Luz, & Ortega-Rubio, 2017; Meerow et al., 2015), the imprints of Andean uplift could be detected in the intraspecific genetic patterns of Neotropical plants with cross-Andean distribution. In this sense, analyzing the genetic structure of plant species would allow inferring the effect of Andean orogenesis on their intraspecific diversification and evolution. Furthermore, this information could shed light about the generation of the highly diverse Neotropical flora.

The tropical Andes are divided into two major geographical units: central Andes and northern Andes (Weigend, 2002). The central Andes (central and southern Peru, northern Chile, and southern Bolivia) uplifted during the Paleogene (65–34 Ma), while the northern Andes (Colombia, Ecuador, and northern Peru) did so during the Neogene (23–2.6 Ma). The central and northern Andes are separated by a low transition zone known as Amotape-Huancabamba which harbors high levels of species diversity and endemism (Borchsenius, 1997; Byg, Vormisto, & Balslev, 2006; Struwe, Haag, Heiberg, & Grant, 2009; Weigend, 2002, 2004). The Andean uplift formerly divided the former continuous rain forest into Amazonia and Chocó forests (Turchetto-Zolet et al., 2012) and promoted the reaccommodation of adjacent foreland basins (Hoorn, Wesselingh, Ter Steege, et al., 2010). Intense mountain building occurred during the middle Miocene (12 Ma) allowing the parallel formation of a new aquatic system known as "Pebas" in the western Amazon basin (Gregory-Wodzicki, 2000; Hoorn, Wesselingh, Hovikoski, & Guerrero, 2010; Meerow et al., 2015; Mora et al., 2010). At that time, the Andes may have started restricting gene flow between eastern and western rain forests on each side of the region, promoting vicariance processes through genetic divergence (Dick, Roubik, Gruber, & Bermingham, 2004; Trénel, Hansen, Normand, & Borchsenius, 2008). After

the Pebas system disappeared and was replaced by the Amazon drainage system (10 Ma) and the establishment of the Amazon River (7 Ma), the Amazon foreland basins became overfilled by the Andean influx of water and sediments over millions of years (Hoorn, Wesselingh, Ter Steege, et al., 2010). This change in the continental hydrological system along with new events of intense tectonism during the Pliocene (5 Ma) allowed for the formation of palaeoarches in the Amazonian lowlands which promoted allopatric diversification between populations located on each side of the arches (Hoorn, Wesselingh, Ter Steege, et al., 2010; Hubert et al., 2007). The basins that formed after the uplift of the arches eventually became overfilled with Andean influx, hiding the arches underground and reconnecting Amazonian isolated biota (Espurt et al., 2010; Hubert et al., 2007).

Despite the increasing evidence that supports cross-Andean genetic divergence in both highland and lowland Neotropical plant populations (Dick, Abdul-Salim, & Bermingham, 2003; Dick & Heuertz, 2008; Hardesty et al., 2010; Motamayor et al., 2008) and also animals (Lovette, 2004), genetic connectivity between cross-Andean regions has been identified (Pérez-Escobar et al., 2017; Rymer, Dick, Vendramin, Buonamici, & Boshier, 2012; Trénel et al., 2008). Historically, biotic exchange has occurred between the Chocó and the Amazon through low passes in the Andes that have functioned as dispersal corridors. The Amotape-Huancabamba zone in southwestern Ecuador/northwestern Peru is a region where dispersal corridors, such as the Huancabamba depression or the Girón-Paute deflection, have facilitated historical cross-Andean dispersal in low- and midelevation lineages in an east-west direction and vice versa (Quintana, Peninngton, Ulloa Ulloa, & Balslev, 2017; Weigend, 2002, 2004) during favorable climatic conditions (Haffer, 1967). Another identified dispersal corridor is Las Cruces mountain pass in the eastern cordillera of Colombia (Dick, Bermingham, Lemes, & Gribel, 2007). These dispersal corridors may have hindered diversification processes in the region by maintaining genetic connectivity between cross-Andean regions (Trénel et al., 2008).

Besides Andean uplift, Pleistocenic climatic shifts (2.5 Ma), as explained in the theory of refugia (Haffer, 1969), were believed to be the main drivers of diversification in the region through a continuous fragmentation process of Amazonian forest. Nevertheless, its role has lately been given less emphasis as both highland and lowland organisms already diversified during Andean orogeny before the Pleistocene (Antonelli & Sanmartín, 2011; Hoorn, Wesselingh, Ter Steege, et al., 2010; Hughes & Eastwood, 2006). Modern rainfall and temperature patterns seem to be related to species richness at large time scales, but their effect on short time scales is often less evident (Antonelli & Sanmartín, 2011; Eiserhardt, Svenning, Kissling, & Balslev, 2011; Field et al., 2009). Additionally, nutrient availability also seems to be an abiotic factor that explains biodiversity accumulation (Antonelli & Sanmartín, 2011; Tuomisto, Zuquim, & Cárdenas, 2014).

*Oenocarpus bataua* is a Neotropical palm that provides a good model for exploring the structuring of genetic diversity at the regional level due to its wide distribution in northern South America
(Balick, 1986), growing in several ecoregions such as the Chocó, the Amazon basin, and the Andean slopes. Despite its wide geographical distribution, its intraspecific variability has been poorly studied with only two allopatric varieties described as follows: (a) *O. bataua* var. *bataua* distributed in northwestern South America (Chocó region, Amazonia basin, and Andean slopes) from sea level to 1,400 m.a.s.l. (Henderson, 1995) and (b) *O. bataua* var. *oligocarpus* distributed in northeastern South America (Guianas, Suriname, eastern Venezuela, eastern Amazonia, and Trinidad Island) at lower altitudes up to 580 m.a.s.l. (Balick, 1986; Henderson, 1995). As it was first described by Martius in 1823, *O. bataua* has had a complex taxonomic history (Montúfar & Pintaud, 2008). Both varieties were originally described as the species *O. bataua* and *Jessenia oligocarpus* (Balick, 1986) but later included as intraspecific variation into the single species *O. bataua* (Henderson, 1995).

Here, we describe the intraspecific genetic structure of the Neotropical palm *O. bataua* in northwestern South America, and we hypothesize possible diversification scenarios linked with the evolution of the Neotropical landscape shaped by Andean orogeny and other drivers of diversification. To do so, we analyzed cross-Andean populations of *O. bataua* var. *bataua* in Colombia, Ecuador, Peru, and

Bolivia using microsatellite markers and Bayesian clustering analyses. Populations of *O. bataua* var. *oligocarpus* from French Guiana were studied only in a first step of the analysis in order to explore the genetic relationship between both varieties. We also compared the levels of genetic diversity and the inbreeding coefficients within populations to provide valuable information for its conservation and management. Finally, we propose additional research in order to improve the knowledge about the diversification and genetics of Neotropical plants. To our knowledge, this is the first work exploring the influence of Andean uplift on the regional genetic structure of a wild palm using its intraspecific genetic diversity.

## 2 | MATERIALS AND METHODS

### 2.1 | Study species

*Oenocarpus bataua* is an arborescent, allogamous, monoecious, and dominant Neotropical palm, reaching 26 m height with stems up to 30 cm in diameter (Dransfield et al., 2008; Henderson, 1995). It displays different habitat preferences at the regional scale. For instance, in the western Amazon basin it grows in well-drained soils in

<figure><img src="image_8.png" /><figcaption>**FIGURE 1** Locations of the 644 *Oenocarpus bataua* samples obtained. Individuals collected in Colombia, Ecuador, Peru, and Bolivia correspond to localities of *O. bataua* var. *bataua*, while the individuals collected in French Guiana correspond to *O. bataua* var. *oligocarpus*. The bigger circumferences represent localities with *n* > 15</figcaption></figure>
terra firme forests, in the central Amazon lowlands it grows in poorly drained soils in swamp forests, and in the intermediate zone it does not show any particular habitat preference (Kahn & de Granville, 1992; Montúfar & Pintaud, 2006). *Oenocarpus bataua* flowers all year round, with peaks depending on geographic region (García, 1998; Núñez-Avellaneda & Rojas-Robles, 2008). Its main pollinators are beetles, but also bees, flies, and even bats may pollinate it (Barfod, Hagen, & Borchsenius, 2011). Its fruits are consumed by birds, primates, and rodents, with the principal long-distance seed dispersers varying in relation to geography (Karubian, Ottewell, Link, & Di Fiore, 2015; Karubian, Sork, Roorda, Duraes, & Smith, 2010). Its fruits and the oil extracted from them are an important nutritional resource for local human populations in the Amazon basin (Balick, 1992; Balslev & Barfod, 1987; Montúfar et al., 2010; Peralta Rivero, Zonta, Moraes, & Ríos, 2008).

## 2.2 | Sample collection

Young leaves and root tissues were collected and stored in plastic bags with silica gel from a total of 644 *O. bataua* individuals at 33 localities through several sampling periods performed from August 2006 to September 2012. Five hundred and ninety-one samples were collected from 28 localities in northwestern South America (Colombia, Ecuador, Peru, and Bolivia) and were identified as var. *bataua*; the remaining 53 samples were collected from five localities in French Guiana and identified as var. *oligocarpus* (Figure 1, Supporting Information Table S1). No samples were collected in the Brazilian Amazon region due to logistic limitations. Both adult and juvenile individuals were collected during the survey. Sampling from closely neighboring plants (distance < 15 m) was avoided to prevent temporal and full-sibling sampling bias.

Sampling was done in a variety of ecosystems, including montane (Andes), semideciduous (southwestern Ecuador), and wet tropical lowland (Chocó and the Amazon) forests. Samples were obtained mostly from primary and secondary forests, although some sampling was done in pastures where deforestation had left few individuals standing. A particular case was the Ecuadorian Amazon locality of Chiriap (a Shuar indigenous community) where *O. bataua* samples were collected from a managed agricultural system.

## 2.3 | DNA extraction, primer screening, and genotyping

Laboratory procedures were performed at IRD's (Institut de Recherche pour le Développement) Genetrop facilities in Montpellier, France. DNA extraction was done using the DNeasy plant mini kit (Qiagen) and its concentration was analyzed with a Thermo Scientific NanoDrop™ ND-100 spectrophotometer. A total of 35 microsatellite markers (simple sequence repeats) of *O. bataua* or related (species within the tribe Euterpeae) were screened for amplification within our sampling (Gaiotto, Brondani, & Grattapaglia, 2001; Lepsch-Cunha, Lund, & Hamilton, 2003; Montúfar, Mariac, Pham, & Pintaud, 2006). Six of these nuclear loci provided highly

polymorphic SSR data and were therefore selected for further analysis (Ob02, Ob06, Ob08, Ob16, AC5-3#4, and AG5-5#1). One additional intron nuclear microsatellite (AG1) was used for genotyping (Ludeña et al., 2011).

SSR loci were amplified in two separate multiplex polymerase chain reactions (PCR): the first included four loci (Ob02, Ob06, Ob08, and Ob16), while the second included three (AC5-3#4, AG5-5#1, and AG1). Forward primers were fluorescently labeled using VIC, NED, 6-FAM, and PET dyes. PCR amplifications were performed in a final volume of 10 µl as follows: 4 µl of DNA template (5 ng/µl), 5 µl of 2x Multiplex Mastermix, and 1 µl of 10x primer mix (containing 3 or 4 pairs of primers at 0.2 µM each primer). Thermal cycling conditions consisted of an initial denaturation step for 15 min at 96°C followed by 30 cycles of denaturation for 30 s at 94°C, annealing for 90 s at 56°C, and elongation for 90 s at 72°C. A final 10 min extension step at 72°C was added. For each PCR product, a final dilution factor of ~1/300 was obtained using 0.15 µl of size marker 500 LIZ (GeneScan™) and water. The final product of each PCR was then analyzed in a 3500 Genetic Analyzer sequencer (Applied Biosystems™), and the resulting chromatograms were examined using GeneMapper software v3.7 (Applied Biosystems™).

## 2.4 | Statistical analysis

For all the sampled localities and individuals, we performed a linkage disequilibrium test for each pair of loci using the software Genepop v4.2 with default parameters (Raymond & Rousset, 1995) and also estimated the presence of null alleles with the software FreeNA (Chapuis & Estoup, 2007).

### 2.4.1 | Individual-based analysis

In order to determine a major population structure in *O. bataua*, a Bayesian clustering method included in the software Structure v2.3 (Pritchard, Stephens, & Donnelly, 2000) was used to assign all 644 individuals to different clusters. To determine the optimal number of clusters (K), Structure was run under the default model of ancestry and population intercorrelation (population admixture and correlated allele frequencies) without prior information about the samples' geographic origin. Five independent Markov chain Monte Carlo (MCMC) runs were performed using 10⁴ burn-in generations followed by 10⁵ sampling generations, K ranging = 1-10. From these data, the ΔK statistic, developed by Evanno, Regnaut, and Goudet (2005), was computed to infer the optimal number of clusters (K).

To identify areas of genetic discontinuity within the distribution range of *O. bataua*, a spatial Bayesian clustering analysis was performed using Geneland software (Guillot, Mortier, & Estoup, 2005). Each individual's geo-referenced and genotypic information was used to determine its posterior probability of belonging to a certain cluster. One MCMC run was performed using the resulting K value obtained from the Structure analysis (following Evanno et al., 2005) applying the following parameters: 10⁵ iterations,
thinning = 1,000, allele frequencies correlated, and with uncertainty in the coordinates. Inbreeding coefficient ($F_{IS}$) and differentiation values ($F_{ST}$) among pairwise clusters were also obtained from Geneland. Based on the Bayesian clustering analyses, we obtained a hierarchical AMOVA in order to understand how the genetic variation is partitioned between varieties, between populations within varieties, and within populations, using the software Arlequin v3.5.2.2 (Excoffier, Laval, & Schneider, 2005) with 1,000 permutations.

We then repeated the Bayesian clustering analyses using only the samples previously identified as var. *bataua* (*n* = 566) in order to infer the population structure within this variety. Next, we used the mean *F* value computed by Structure to explore the dispersal history among genetic clusters of var. *bataua*. Under population admixture and correlated allele frequencies, the program returns an *F* value (*F*_{ST} analogue) that describes the degree of genetic differentiation of a certain cluster from a hypothetical ancestral population (Falush, Stephens, & Pritchard, 2003). The dispersal history of the clusters can be inferred as a path from low to high *F* values, assuming constant rates of genetic drift in all the clusters (Trénel et al., 2008). The mean *F* value of each cluster identified within the var. *bataua* was obtained after averaging the *F* values of the five independent MCMC with the determined optimal number of clusters (*K*). Additionally, phylogenetic relationships between clusters were depicted to determine the sequence of divergence between the clusters. A neighbor-joining (NJ) tree was constructed in MEGA7 (Kumar, Stecher, & Tamura, 2016) using a mean matrix of allele frequency divergence among clusters (net nucleotide distance) that resulted from the analysis in Structure. The robustness of the NJ branches

was evaluated using PHYLIP v3.6 (Felsenstein, 2005) through 1,000 bootstrap replications.

## 2.4.2 | Population-based analysis

At this level, we worked with 18 localities (*n* > 15) previously identified as var. *bataua* by the first Bayesian clustering analyses, and each locality was treated as a population. In order to determine how genetic diversity was spatially distributed, allelic richness (*A*) was calculated using the rarefaction procedure implemented in FSTAT v2.9.3 (Goudet, 2001) in order to control for the effect of the different sampling size between populations. A simple linear regression between *A* and altitude from each population was performed to determine its association. Another linear regression was done between *A* and the distance from each population to the population with the highest *A* value (most diverse) to check whether genetic diversity decreased as we moved away from a center of diversity. Expected heterozygosity (*H*<sub>O</sub>), observed heterozygosity (*H*<sub>D</sub>), and the inbreeding coefficient (*F*<sub>IS</sub>) for each population were estimated with Arlequin.

# 3 | RESULTS

Seven of the 21 pairs of loci were detected in linkage disequilibrium (*p* < 0.05; Supporting Information Table S2). All loci showed low average estimates of null allele frequency (AC5-3#4 = 0.024, AG5-5#1 = 0.013, AG1 = 0.022, Ob02 = 0.059, Ob06 = 0.057, Ob08 = 0.071, Ob16 = 0.068). The presence of null alleles did not affect our data as the *F*<sub>ST</sub> automatically generated for all loci by the

<figure><img src="image_9.png" /><figcaption>**FIGURE 2** Genetic clusters identified within *Oenocarpus bataua* (var. *bataua* and var. *oligocarpus*) using 644 samples. These were identified with a spatial Bayesian clustering analysis conducted in Geneland (Guillot et al., 2005) using posterior probabilities to belong to one of *K* = 2 clusters as identified in Structure (Pritchard et al., 2000) with the statistical analysis developed by Evanno et al. (2005). Each point represents a sampled locality, while the lines represent the probability of membership to a determined cluster</figcaption></figure>
software FreeNA was similar before (0.179) and after (0.175) correction for null alleles.

## 3.1 | Individual-based analysis

Based on the ΔK statistic (Evanno et al., 2005), we determined the best K at K = 2 (Supporting Information Figure S1) with a clear separation between eastern (French Guiana, var. *oligocarpus*) and western (Chocó, Andean, and western Amazonian forests, var. *bataua*) populations. An exception to this pattern was the relatedness of the San Francisco population (Colombia), with the individuals identified as var. *oligocarpus* by the Bayesian analyses (Figure 2; Supporting Information Figure S2), even though it was identified in the field as var. *bataua*. Therefore, this population was not included in the following analyses on var. *bataua*. Results produced in Geneland (Figure 2) showed a strong assignment of the individuals to each corresponding cluster (var. *bataua* and var. *oligocarpus*) with posterior probabilities of membership of 1. Despite the clear divergence

between the two varieties, the $F_{ST}$ value obtained was moderate (Table 1). The AMOVA showed that most of the variation (73.23%) was harbored within populations, as is expected for an allogamous and long-lived perennial species (Hamrick & Godt, 1990). Genetic variation between varieties was 13.65% and between populations within varieties was 13.12% (all *p* < 0.01**).

Within var. *bataua*, the Structure analysis identified a peak at K = 4 (Supporting Information Figure S1). The bar plot generated by Structure (Supporting Information Figure S3) and the Geneland analysis (Figure 3) assigned the samples to four distinct groups that represent four regions of particular biodiversity: Chocó rain forest (hereafter CHO), Amotape-Huancabamba Zone (AMO), northwestern Amazonian rain forest plus northwestern Bolivia (NWA), and southwestern Amazonian rain forest (SWA) (Figure 3). Based on the Structure bar plot, historical gene flow has been maintained between clusters as a great number of individuals show some probability of belonging to more than one cluster. The first cluster, CHO, is formed by all samples located within the Chocó rain forest. The AMO cluster is formed by cross-Andean populations located in the Andean slopes of southern Ecuador/northern Peru, suggesting genetic connectivity through the Andes as the Geneland analysis (Figure 3) and the Structure bar plot (Supporting Information Figure S3) showed. The populations from the western Amazon were split among two different clusters: (a) the NWA cluster that included populations along the Napo River basin in Colombia, Ecuador, and Peru, and populations in northwestern Bolivia; (b) the SWA cluster was formed by all samples from southwestern Amazonia along the Ucayali and Madre de Dios River basins in Peru. The average *F* values obtained for the clusters in Structure suggest that the NWA and SWA clusters experienced

<table><caption>**TABLE 1** *F*<sub>IS</sub> (inbreeding coefficient) and pairwise *F*<sub>ST</sub> (fixation index) values obtained from Arlequin (Excoffier et al., 2005) for the *bataua* and *oligocarpus* genetic clusters</caption><thead><tr><th rowspan="2">Variety</th><th rowspan="2">n</th><th rowspan="2">F<sub>IS</sub></th><th colspan="2">F<sub>ST</sub></th></tr><tr><th>bataua</th><th>oligocarpus</th></tr></thead><tbody><tr><td>bataua</td><td>566</td><td>0.144</td><td>-</td><td>-</td></tr><tr><td>oligocarpus</td><td>78</td><td>0.327</td><td>0.167</td><td>-</td></tr></tbody></table>

Note. n = sample size.

<figure><img src="image_10.png" /><figcaption>**FIGURE 3** Four genetic clusters identified within *Oenocarpus bataua* var. *bataua* (AMO: Amotape-Huancabamba zone; CHO: Chocó rain forests; NWA: northwestern Amazonia rain forests + northwestern Bolivia; SWA: southwestern Amazonia rain forests) using 566 samples. These were identified with a spatial Bayesian clustering analysis conducted in Geneland (Guillot et al., 2005) using posterior probabilities to belong to one of K = 4 clusters as identified in Structure (Pritchard et al., 2000) with the statistical analysis developed by Evanno et al. (2005). Each point represents a sampled locality, while the lines represent the probability of membership to a determined cluster</figcaption></figure>
<table><caption>**TABLE 2** Inbreeding coefficients ($F_{IS}$) and pairwise fixation values ($F_{ST}$) obtained from Geneland (Guillot et al., 2005) for the four genetic clusters identified in Structure (Pritchard et al., 2000) within *Oenocarpus bataua* var. *bataua*. Divergence estimates from a hypothetical ancestral population (F) were obtained from Structure (Pritchard et al., 2000)</caption><thead><tr><th rowspan="2">Cluster</th><th rowspan="2">F<sub>IS</sub></th><th rowspan="2">F</th><th colspan="4">F<sub>ST</sub></th></tr><tr><th>CHO</th><th>AMO</th><th>NWA</th><th>SWA</th></tr></thead><tbody><tr><td>CHO</td><td>0.018</td><td>0.191</td><td>-</td><td></td><td></td><td></td></tr><tr><td>AMO</td><td>0.105</td><td>0.163</td><td>0.098</td><td>-</td><td></td><td></td></tr><tr><td>NWA</td><td>0.186</td><td>0.031</td><td>0.087</td><td>0.058</td><td>-</td><td></td></tr><tr><td>SWA</td><td>0.097</td><td>0.053</td><td>0.134</td><td>0.091</td><td>0.050</td><td>-</td></tr></tbody></table>

Note. AMO: Amotape-Huancabamba zone; CHO: Chocó rain forests; NWA: northwestern Amazonia rain forests + northwestern Bolivia; SWA: southwestern Amazonia rain forests.

the lowest genetic drift after the ancestral population subdivided, whereas CHO experienced the highest. The highest genetic differentiation ($F_{ST}$ values) between clusters was found between CHO compared to SWA and CHO compared to AMO, whereas the lowest was between NWA and SWA (Table 2), showing that *O. bataua* populations within the Chocó region have diverged more than the others. These results were consistent with the phylogenetic analysis (Figure 4) in which NWA and SWA were the closest evolutionary units, and CHO was the most distant unit.

## 3.2 | Population-based analysis

The lowland Amazonian populations of Intuto and Jenaro Herrera, located near Iquitos in Peru within the NWA cluster, had the highest genetic diversity (Table 3). Other Amazonian populations such as Pucallpa and Iñapari, located within the SWA cluster, also had high genetic diversity. Populations at higher elevations such as Pacifica and Zamora, located within the cluster AMO, showed less diversity than those in lowland Amazonia. The lowest diversity was found in the westernmost populations within the CHO cluster and within Villaseca (AMO cluster) in southwestern Ecuador. The population of San Buenaventura in Bolivia, located at the southern limit of our sampling, showed low diversity in terms of allelic richness. No significant regression was found between A and altitude ($p = 0.062$) or between A and the distance from each population to the most diverse population (Intuto; $p = 0.058$); however, their $p$ values were in the threshold of statistical significance. As all trans-Andean (west of the Andes) populations (Esmeraldas, El Chontal, Bilsa, Villaseca) showed values of A under the regression line (Supporting Information Figure S4), we repeated the procedures obviating these populations. Then, a significant correlation ($p = 0.045^{**}$) was detected between A and the distance from each population to Intuto, whereas no significant correlation was found between A and altitude ($p = 0.054$).

The inbreeding coefficient tended to be low with few exceptions (Table 3). The highest $F_{IS}$ value was observed in Villaseca. High and significant $F_{IS}$ values were also found in the Andean populations

<figure><img src="image_11.png" /><figcaption>**FIGURE 4** A neighbor-joining analysis showing the phylogenetic relationships between the four genetic clusters identified within *Oenocarpus bataua* var. *bataua* (AMO: Amotape-Huancabamba zone; CHO: Chocó rain forests; NWA: northwestern Amazonia rain forests + northwestern Bolivia; SWA: southwestern Amazonia rain forests) using 566 samples. It was performed in MEGA (Kumar et al., 2016) using a mean matrix of allele frequency divergence among clusters (net nucleotide distance) resulted from the analysis in Structure (Pritchard et al., 2000) that determined genetic clusters within *Oenocarpus bataua* var. *bataua* populations. The robustness of the neighbor-joining branches was evaluated using PHYLIP (Felsenstein, 2005) through 1,000 bootstrap replications</figcaption></figure>

Shuaro and Rioja in Peru and for the Amazonian population of San Buenaventura in Bolivia, where deforestation processes have occurred recently and were observed during sampling. An unexpected moderate inbreeding coefficient was observed in Yasuní in western Ecuador, located within a protected national park. A particular case with excess of heterozygosity came from Chiriap in Ecuador.

# 4 | DISCUSSION

## 4.1 | Andean uplift as possible driver of divergence

The four genetic clusters identified in the Bayesian analysis within var. *bataua* (Figure 3) correlated to major ecoregions recognized within northwestern South America (Dinerstein et al., 1995; Olson & Dinerstein, 2002; Weigend, 2002). We hypothesize that Andean uplift promoted three events of diversification that shaped the genetic structure of var. *bataua* into four clusters. Despite not being able to prove this hypothesis with our data, we will explore possible orogenic scenarios that can explain the divergence observed. We suggest that a first diversification event in var. *bataua* occurred between the Chocó region and the Amazon basin. Cross-Andean divergence has been also reported for several rain forest trees in the Neotropics (Dick & Heuertz, 2008; Dick et al., 2003, 2007; Hardesty et al., 2010; Motamayor et al., 2008; Rymer et al., 2012). Although this seems like a logical explanation, we cannot discard the possibility that the cross-Andean distribution of *O. bataua* may be due to long-distance dispersal processes after the Andes reached its current height just 2.7 Ma (Gregory-Wodzicki, 2000; Mora et al., 2010). In this sense, a dated phylogeny of *O. bataua* populations would help to elucidate whether this diversification event, and the other two detected, shares a time frame with the Andean orogenic
<table><thead><tr><th>Locality</th><th>Cluster</th><th>n</th><th>A</th><th>H<sub>e</sub></th><th>H<sub>o</sub></th><th>F<sub>IS</sub></th></tr></thead><tbody><tr><td>Esmeraldas</td><td>CHO</td><td>16</td><td>6.191</td><td>0.70</td><td>0.81</td><td>-0.158</td></tr><tr><td>El Chontal</td><td>CHO</td><td>32</td><td>5.671</td><td>0.70</td><td>0.75</td><td>-0.074</td></tr><tr><td>Bilsa</td><td>CHO</td><td>31</td><td>5.653</td><td>0.73</td><td>0.84</td><td>-0.172</td></tr><tr><td>Villaseca</td><td>CHO/AMO</td><td>19</td><td>4.823</td><td>0.77</td><td>0.47</td><td>0.377**</td></tr><tr><td>Pachicusa</td><td>AMO</td><td>19</td><td>6.475</td><td>0.73</td><td>0.70</td><td>0.031</td></tr><tr><td>Zamora</td><td>AMO</td><td>57</td><td>6.074</td><td>0.71</td><td>0.74</td><td>-0.046</td></tr><tr><td>Rioja</td><td>AMO/NWA</td><td>18</td><td>7.718</td><td>0.76</td><td>0.61</td><td>0.157**</td></tr><tr><td>Yasuní</td><td>NWA</td><td>30</td><td>6.967</td><td>0.75</td><td>0.60</td><td>0.17**</td></tr><tr><td>Pantoja</td><td>NWA</td><td>26</td><td>7.294</td><td>0.74</td><td>0.65</td><td>0.038</td></tr><tr><td>Chiriap</td><td>AMO/NWA</td><td>30</td><td>6.917</td><td>0.72</td><td>0.88</td><td>-0.314</td></tr><tr><td>Intuto</td><td>NWA</td><td>32</td><td>9.523</td><td>0.86</td><td>0.82</td><td>0.034</td></tr><tr><td>Jenaro Herrera</td><td>NWA</td><td>30</td><td>9.371</td><td>0.81</td><td>0.86</td><td>-0.066</td></tr><tr><td>Pucallpa</td><td>SWA</td><td>43</td><td>8.072</td><td>0.76</td><td>0.75</td><td>0.007</td></tr><tr><td>Ahuaytía</td><td>SWA</td><td>32</td><td>6.385</td><td>0.67</td><td>0.66</td><td>0.01</td></tr><tr><td>Shuaro</td><td>SWA</td><td>15</td><td>7.054</td><td>0.71</td><td>0.53</td><td>0.242**</td></tr><tr><td>Iñapari</td><td>SWA</td><td>30</td><td>7.848</td><td>0.77</td><td>0.82</td><td>-0.064</td></tr><tr><td>Tambopata</td><td>SWA</td><td>35</td><td>6.952</td><td>0.72</td><td>0.68</td><td>0.045</td></tr><tr><td>San Buenaventura</td><td>NWA</td><td>17</td><td>5.806</td><td>0.75</td><td>0.60</td><td>0.23**</td></tr></tbody></table>

Notes. AMO: Amotape-Huancabamba zone; CHO: Chocó rain forests; NWA: northwestern Amazonia rain forests + northwestern Bolivia; SWA: southwestern Amazonia rain forests.

n = sample size.

**Highly significant (p < 0.01).

<caption>TABLE 3 Diversity values and inbreeding coefficients (F<sub>IS</sub>) for the 18 *Oenocarpus bataua* var. *bataua* populations with n > 15. Allelic richness (A) was calculated using the rarefaction procedure implemented in FSTAT (Goudet, 2001), whereas expected (H<sub>e</sub>) and observed (H<sub>o</sub>) heterozygosity, and the inbreeding coefficient (F<sub>IS</sub>) were obtained from Arlequin (Excoffier et al., 2005)</caption>

events discussed or whether they happened more recently during the Pliocene.

Despite the presence of the Andes, genetic connectivity seems to occur between cross-Andean populations in south- ern Ecuador/northern Peru as the Bayesian clustering analyses showed. The disjunctive distribution of the cluster AMO sug- gested that cross-Andean populations within the Amotape- Huancabamba zone have maintained genetic connectivity across the Andes through dispersal corridors. This cluster may have developed due to a second diversification process influenced by the Andes that occurred when populations within the Amotape- Huancabamba zone diverged from those in the Amazon basin and the Chocó region. A similar pattern of genetic structuring was reported for *Theobroma cacao* (Motamayor et al., 2008; Thomas et al., 2012), where cross-Andean populations located in southern Ecuador were grouped into a single cluster. The same trend was reported for a clade within the genus *Macrocarpaea* in the family Gentianaceae (Struwe et al., 2009). Apparently, the Amotape- Huancabamba zone could be an intermediate zone that receives genetic information from both trans- (west of the Andes) and cis- Andean (east of the Andes) populations due to the presence of cross-Andean dispersal corridors.

A more recent diversification process was identified in the western Amazon basin due to the presence of two distinct genetic clusters, NWA and SWA. The split between these two clusters may have occurred when the western Amazon drainage basin divided into two foreland basins (Roddaz, Viers, Brusset, Baby, & Hérail,

2005) due to the uplift of a palaeoarch known as Fitzcarrald Arch in central Peru during the Pliocene (~4 Ma; Espurt et al., 2010). This uplift, associated with Andean tectonics, created a NE-SW- trending barrier for gene flow in the western Amazon basin that could have promoted a third diversification event influenced by Andean uplift. The geographical division of the clusters NWA (excluding the northwestern Bolivia section) and SWA correlated with the historical division of the northern and southern foreland basins. After the basins overfilled, no geographical barriers have been present in the western Amazon basin (Espurt et al., 2010); however, the intraspecific diversification between both clusters has maintained to the present. The location of these two ancient basins presents different climatic conditions currently, being northwestern Amazonia a more humid and less seasonal region than southwestern Amazonia (Silman, 2007). Therefore, current climatic dynamics may be contributing to the maintaining of in- traspecific genetic diversification in *O. bataua* as palms are highly sensitive to climatic conditions (Eiserhardt et al., 2011). Variation in climate can influence flowering phenology among populations (Welt, Litt, & Franks, 2015), which may alter their gene flow pat- terns (Franks & Weis, 2009) and even promote reproductive isola- tion (Martin, Bouck, & Arnold, 2007). It is worth mentioning that the location of the clusters NWA and SWA partially correlates with the location of two Pleistocenic forest refuges (Napo and East Peruvian) proposed by Haffer (1969); however, as the theory of refugia was shown to be based on sampling artifacts (Nelson, Ferreira, da Silva, & Kawasaki, 1990), we strongly support the
hypothesis that these clusters developed after geographical isolation and posterior reproductive isolation. Therefore, the determination of the flowering phenology of the populations within the clusters NWA and SWA would shed light on the influence of climatic dynamics on the genetic structuring of *O. bataua* in the western Amazon basin.

The apparent disjunct distribution of the NWA cluster may be a sampling artifact caused by lack of sampling in more eastern localities. It is possible that the two sections of this cluster are linked by unsampled areas in western Brazil. A similar pattern of disjunct distribution was reported for *T. cacao* (Motamayor et al., 2008), where Amazonian populations in northern Peru were genetically related to populations in southwestern Brazil.

## 4.2 | Is var. *oligocarpus* distributed beyond eastern Venezuela?

The assignment of the San Francisco population (Magdalena River valley, between the central and eastern Andean cordilleras of Colombia) to the *oligocarpus* cluster was a surprising result. Geographically, var. *oligocarpus* is distributed in the eastern Amazon and has not been reported in northwestern South America where San Francisco is located. The particular genetic footprint of the San Francisco population is strong enough to be interpreted as a population closely related to var. *oligocarpus*. Nevertheless, other explanations could be also analyzed in following studies, such as the potential hybridization process between var. *bataua* and *oligocarpus*, or even between *O. bataua* and *O. minor*, whose hybrids have been reported in the Magdalena River valley (Núñez-Avellaneda, 2007).

Our study confirmed the presence of two intraspecific varieties in *O. bataua* that have experienced strong genetic divergence. The geographic distribution of var. *oligocarpus* coincided with the location of the Guiana Shield, which is made up of exposed Precambrian rock that formed about 2 Ga. The geology of this area allowed for the development of a forest with a different floristic composition compared to the forests of western and central Amazonia (Gibbs & Barron, 1993). Therefore, the divergence between var. *bataua* and var. *oligocarpus* was not related to the uplift of the Andes. The molecular differentiation found between these two varieties agreed with previous studies (Montúfar, 2007; Montúfar & Pintaud, 2008). Despite the strong differentiation between the two varieties, our data were insufficient to support the hypothesis of the botanists A. Grisebach and H. Wendland, who originally described *J. oligocarpus* (var. *oligocarpus*) in 1864 as a distinctive species from *O. bataua* (Balick, 1986). The study of gene flow between the two varieties in sympatric zones (probably within the Brazilian and Venezuelan Amazon) would enhance our knowledge about the genetic patterns between them. Furthermore, the implementation of studies that determine the presence or absence of reproductive isolation (floral morphology, phenology) between the two varieties would help to cast light on their biological divergence.

## 4.3 | Genetic diversity hot spot in northwestern Amazonia

Our results added evidence to a pattern of geographic correlation between high genetic diversity and high species diversity (Lankau & Strauss, 2007; Palma-Silva et al., 2009; Vellend & Geber, 2005); however, we did not test this correlation statistically. In this study, we reported a genetic diversity hot spot of *O. bataua* in the Peruvian rain forest around Iquitos (Intuto and Jenaro Herrera populations), which spatially coincides with the high palm and tree species diversity found in these forests (Alvez-Valles et al., 2017; Gentry, 1988; Kristiansen et al., 2011; Ter Steege et al., 2013; Valencia, Balslev, & Paz y Miño, 1994; Vormisto, Svenning, Hall, & Balslev, 2004). High genetic diversity has been also identified in cultivated cacao within this area (Motamayor et al., 2008; Sereno, Albuquerque, Vencovsky, & Figueira, 2006; Thomas et al., 2012). The high genetic diversity of *O. bataua* also was consistent with the high species diversity of the genus *Oenocarpus* in this region. The locality of La Pedrera in southeastern Colombia (~350 km from Iquitos) harbors the highest *Oenocarpus* species diversity, with six described species (Bernal, Galeano, & Henderson, 1991; Galeano & Bernal, 2010), showing that northwestern Amazonian forests and their surroundings are a center of species and genetic diversity for this genus.

Within the Amazon region, the genetic diversity of *O. bataua* peaks in the northern Peruvian Amazon and then decreases westward and southward as we move away from this area. This pattern is statistically weak when trans-Andean populations are included in the analysis as the Andes may interfere with long-distance dispersal events between cross-Andean regions (Dick et al., 2004), altering the genetic patterns between regions. This gene flow constraint may explain why trans-Andean populations harbor lower diversity than expected. This pattern of genetic diversity added evidence on the role of northwestern Amazonian rain forests as a center of diversity in *Oenocarpus*.

The high diversity harbored within northwestern Amazonia may be related to high resource availability expressed as annual rainfall and soil cation concentration (Antonelli & Sanmartín, 2011; Tuomisto et al., 2014). This region presents high rates of water availability and climate stability due to convective rain caused by the Andes even during glacial periods (Kristiansen et al., 2011; Pitman, 2000; Tuomisto et al., 2014). It also harbors higher nutrient soils than central and eastern Amazonia due to the deposition and accumulation of material eroded during Andean orogeny (Higgins et al., 2011; Hoorn, Wesselingh, Ter Steege, et al., 2010; Tuomisto et al., 2014). It is possible that diversity depends on resource availability; however, we have ignored the specific mechanisms for this.

*Oenocarpus bataua* maintained medium levels of genetic diversity within its populations when compared with other Neotropical palms genotyped with microsatellite markers. Average genetic diversity of *O. bataua* populations in terms of allelic richness (A = 6.93) was lower when compared with other lowland palms such as *Euterpe edulis* (A = 14.46; Conte, Sedrez dos Reis, Mantovani, & Vencovsky,
2008) or *Bactris gasipaes* (A = 7.5; Galluzzi et al., 2015). On the other hand, *O. bataua* showed higher diversity when compared with the Andean palm *Ceroxylon echinulatum* (A = 1.82; Trénel et al., 2008). Although there is variation in the experimental designs of the studies mentioned above, they allow inferring the levels of genetic diversity in *O. bataua* compared to other Neotropical palms.

High inbreeding coefficients were encountered in populations growing in forests under high deforestation pressures, such as Shuaro, Rioja, or San Buenaventura. A particular high inbreeding coefficient was determined for the population of Villaseca (southeastern Ecuador), located in an isolated and heavily deforested area. This population also represents the edge of the species' ecological range, growing in semideciduous forest; therefore, its conservation would secure this ecological adaptation of *O. bataua*. In fragmented areas, the inbreeding coefficient increases along with linkage disequilibrium (Zartman, McDaniel, & Shaw, 2006). Thus, the inbreeding coefficients reported here could explain the linkage disequilibrium detected on one-third of pairs of loci. Forest degradation and fragmentation produces high inbreeding within populations (Rhoads, Williams, & Krane, 2017); however, as the generation time of *O. bataua* may be of more than 50 years, we could not relate the observed inbreeding coefficients with recent deforestation processes. Furthermore, the high inbreeding coefficient found in the Yasuní population was also not related to deforestation as it is located within a protected area. It is possible that these populations may be under the Wahlund effect which refers to an excess of homozygotes (high F<sub>IS</sub> values) due to population structuring (Garnier-Géré & Chikhi, 2013); however, a detailed study is needed to confirm this. On the other hand, an excess of heterozygosity was determined in the Chiriap population in the Ecuadorian Amazon, which is managed by local indigenous people who disperse seeds of different provenances of the Amazonian forest surroundings in order to maintain a healthy population for lipid extraction; however, this value was not significant.

## 4.4 | Conclusions and future perspectives

We detected three events of genetic diversification within *O. bataua* var. *bataua* that promoted the structuring of four genetic clusters in northwestern South America. In order to determine the effective role of the Andean orogenesis as a driver of intraspecific genetic diversification in the region, we encourage the realization of dated phylogenies or coalescence analyses in other Neotropical plant species with similar distribution to *O. bataua* in a similar or bigger geographical extent as this study. Based on our results, we can suggest that Andean uplift may be the main driver of genetic diversification for this palm in northwestern South America, showing the influence of this orogenic event on the evolution and diversification of Neotropical flora. As gene flow is present among the identified genetic clusters of var. *bataua*, it seems improbable that one of these clusters would diversify as a different intraspecific entity (variety) unless gene flow becomes restricted. Addressing the role of dispersal corridors in the genetic structure of widely distributed plants would provide

more accurate information about the influence of the Andes in the diversification of Neotropical flora. Additional studies are necessary to understand whether populations with high inbreeding coefficients are under the Wahlund effect. Furthermore, we confirmed the existence of two varieties with strong divergence within the wide distribution of *O. bataua* whose divergence is not related to Andean orogenesis. We strongly recommend a more detailed exploring of genetic diversification in the northern Andes in Colombia and Venezuela in order to have a better understanding of the influence of the Andes within this region. Furthermore, we suggest the study of the flowering phenology of the populations within the NWA and SWA clusters to determine whether differences in climate are altering flowering times and gene flow patterns in these populations.

## ACKNOWLEDGMENTS

We would like to thank three anonymous reviewers for enhancing the quality of a former version of this manuscript. We thank Cédric Mariac, Marie Couderc, Julissa Roncal, and the entire DYNADIV team from the IRD (Montpellier, France) for helping during laboratory work. Special thanks go to Jean-Jacques de Granville (IRD) for the valuable support during the sampling in French Guiana; to José Manuel Barreiro, María José Sanín, Gabriel Rivadeneira, Luis López, and John Colorado for assistance during fieldwork; and to Christine Fram for grammar and style correction of the manuscript. FP7-PALMS-project (Palms Harvest Impacts in Tropical Forests; EU-contract 212631 and the Danish Council for Independent Research-Natural Sciences [grant#4181-00158] to Henrik Balslev) and PUCE (Project "Oenocarpus bataua") supported this research financially. The Ministerio del Ambiente-Ecuador (N°004-IC-FLO-DPPZ/MA, N°06-2012-INVESTIGACIÓN-B-DP/MS/MAE, N°009-IC-FAU/FLO-DPZCH-MA), Dirección General Forestal y de Fauna Silvestre-Peru (0118-2009-AG-DGFFS-DGEFFS), Ministerio de Medio Ambiente y Agua-Bolivia (1088), and Autoridad Nacional de Licencias Ambientales-Colombia (ANLA 0255) gave us permission to collect the needed samples. Samples from French Guiana were obtained with the collaboration of the IRD, France.

## AUTHOR CONTRIBUTION

R. M., J.-C. P., and H. B. developed the study; S. E., R. M., and J.-C. P. made the sampling; J.-C. P., R. B., M. M. R., and B. M. provided logistics for sampling and obtained collection permits; J.-C. P. provided logistics for laboratory work; S. E. performed laboratory work and did the statistical analyses. S. E. and R. M. wrote the manuscript; R. M., H. B., R. B., M. M. R., B. M., and J.-C. P. corrected the manuscript.

## DATA ACCESSIBILITY

Microsatellite and geographic data are available at https://doi.org/10.5061/dryad.1r4p8.
ORCID

Rommel Montúfar ID http://orcid.org/0000-0001-5309-4889

REFERENCES

Alvez-Valles, C. M., Balslev, H., Carvalho, F. A., Garcia-Villacorta, R., Grandez, C., & Neto, L. M. (2017). Endemism and conservation of Amazon palms. *Biodiversity and Conservation*, 27(3), 1-20.

Antonelli, A., Nylander, J. A., Persson, C., & Sanmartín, I. (2009). Tracing the impact of the Andean uplift on Neotropical plant evolution. *Proceedings of the National Academy of Sciences*, 106(24), 9749-9754.
https://doi.org/10.1073/pnas.0811421106

Antonelli, A., & Sanmartín, I. (2011). Why are there so many plant species in the Neotropics? *Taxon*, 60(2), 403-414.

Balick, M. J. (1986). *Systematics and economic botany of the Oenocarpus-Jessenia (Palmae) complex*. Advances in Economic Botany 3. Bronx, NY: New York Botanical Garden.

Balick, M. J. (1992). *Jessenia y Oenocarpus: Palmas aceiteras neotropicales dignas de ser domesticadas*. Rome, Italy: Organización de la Naciones Unidas para la Agricultura y la Alimentación (FAO).

Balslev, H. (1988). Distribution patterns of Ecuadorean plant species. *Taxon*, 37, 567-577. https://doi.org/10.2307/1221100

Balslev, H., & Barfod, A. (1987). Ecuadorian palms - An overview. *Opera Botanica*, 92, 17-35.

Barfod, A. S., Hagen, M., & Borchsenius, F. (2011). Twenty-five years of progress in understanding pollination mechanisms in palms (Arecaceae). *Annals of Botany*, 108, 1503-1516. https://doi.org/10.1093/aob/mcr192

Bernal, R. G., Galeano, C., & Henderson, A. (1991). Notes on *Oenocarpus* (Palmae) in the Colombian Amazon. *Brittonia*, 43(3), 154-164. https://doi.org/10.2307/2807045

Borchsenius, F. (1997). Patterns of plant species endemism in Ecuador. *Biodiversity and Conservation*, 6, 379-399. https://doi.org/10.1023/A:1018312724137

Byg, A., Vormisto, J., & Balslev, H. (2006). Using the useful: Characteristics of used palms in south-eastern Ecuador. *Environment, Development and Sustainability*, 8, 495-506. https://doi.org/10.1007/s10668-006-9051-6

Chapuis, M.-P., & Estoup, A. (2007). Microsatellite null alleles and estimation of population differentiation. *Molecular Biology and Evolution*, 24(3), 621-631. https://doi.org/10.1093/molbev/msl191

Conte, R., Sedrez dos Reis, M., Mantovani, A., & Vencovsky, R. (2008). Genetic structure and mating system of *Euterpe edulis* Mart. populations: A comparative analysis using microsatellite and allozyme markers. *Journal of Heredity*, 99(5), 476-482. https://doi.org/10.1093/jhered/esn055

Dick, C. W., Abdul-Salim, K., & Bermingham, E. (2003). Molecular systematics analysis reveals cryptic Tertiary diversification of a widespread tropical rainforest tree. *American Naturalist*, 162, 691-703. https://doi.org/10.1086/379795

Dick, C. W., Bermingham, E., Lemes, M. R., & Gribel, R. (2007). Extreme long-distance dispersal of the lowland tropical rainforest tree *Ceiba pentandra* L. (Malvaceae) in Africa and the Neotropics. *Molecular Ecology*, 16(14), 3039-3049. https://doi.org/10.1111/j.1365-294X.2007.03341.x

Dick, C. W., & Heuertz, M. (2008). The complex biogeographic history of a widespread tropical tree species. *Evolution*, 62(11), 2760-2774. https://doi.org/10.1111/j.1558-5646.2008.00506.x

Dick, C. W., Roubik, D. W., Gruber, K. F., & Bermingham, E. (2004). Long-distance gene flow and cross-Andean dispersal of lowland rainforest bees (Apidae: Euglossini) revealed by comparative mitochondrial DNA phylogeography. *Molecular Ecology*, 13, 3775-3785. https://doi.org/10.1111/j.1365-294X.2004.02374.x

Dinerstein, E., Olson, D. M., Graham, D. J., Webster, A. L., Primm, S. A., Bookbinder, M. P., & Ledec, G. (1995). *A conservation assessment of the terrestrial ecoregions of Latin America and the Caribbean*. Washington, DC: The World Bank. https://doi.org/10.1596/0-8213-3295-3

Dransfield, J., Uhl, N. W., Asmussen, C. B., Baker, W. J., Harley, M. M., & Lewis, C. E. (2008). *Genera Palmarum*. Kew, UK: Royal Botanic Gardens.

Eiserhardt, W. L., Svenning, J. C., Kissling, W. D., & Balslev, H. (2011). Geographical ecology of the palms (Arecaceae): Determinants of diversity and distributions across spatial scales. *Annals of Botany*, 108(8), 1391-1416. https://doi.org/10.1093/aob/mcr146

Espurt, N., Baby, P., Brusset, S., Roddaz, M., Hermoza, W., & Barbarand, J. (2010). The Nazca Ridge and uplift of the Fitzcarrald Arch: Implications for regional geology in Northern South America. In C. Hoorn, & F. P. Wesselingh (Eds.), *Amazonia, landscape and species evolution: A look into the past* (pp. 89-100). Hoboken, NJ: Blackwell Publishing.

Evanno, G., Regnaut, S., & Goudet, J. (2005). Detecting the number of clusters of individuals using the software structure: A simulation study. *Molecular Ecology*, 14, 2611-2630. https://doi.org/10.1111/j.1365-294X.2005.02553.x

Excoffier, L., Laval, G., & Schneider, S. (2005). ARLEQUIN version 3.0: An integrated software package for population genetics data analysis. *Evolutionary Bioinformatics Online*, 1, 47-50.

Falush, D., Stephens, M., & Pritchard, J. (2003). Inference of Population Structure Using Multilocus Genotype Data: Linked Loci and Correlated Allele Frequencies. *Genetics Society of America*, 164, 1567-1587.

Felsenstein, J. (2005). *PHYLIP (Phylogeny Inference Package) version 3.6*. Seattle, WA: Department of Genome Sciences, University of Washington.

Field, R., Hawkins, B. A., Cornell, H. V., Currie, D. J., Diniz-Filho, J. A. F., Guégan, J. F., ... Turner, J. R. G. (2009). Spatial species-richness gradients across scales: A meta-analysis. *Journal of Biogeography*, 36, 132-147. https://doi.org/10.1111/j.1365-2699.2008.01963.x

Franks, S. J., & Weis, A. E. (2009). Climate change alters reproductive isolation and potential gene flow in an annual plant. *Evolutionary Applications*, 2(4), 481-488. https://doi.org/10.1111/j.1752-4571.2009.00073.x

Gaiotto, F., Brondani, R., & Grattapaglia, D. (2001). Microsatellite markers for heart of palm- *Euterpe edulis* and *E. oleracea* Mart. (Arecaceae). *Molecular Ecology*, 1(1-2), 86-88.

Galeano, G., & Bernal, R. G. (2010). *Palmas de Colombia. Guia de campo*. Bogotá, Colombia: Editorial Universidad Nacional de Colombia.

Galluzzi, G., Dufour, D., Thomas, E., Van Zonneveld, M., Salamanca, A. F. E., Toro, A. G., ... González-Mejía, A. (2015). An integrated hypothesis on the domestication of *Bactris gasipaes*. *PLoS One*, 10(12), e0144644. https://doi.org/10.1371/journal.pone.0144644

García, M. (1998). *Observaciones de polinización en Jessenia bataua* (Arecaceae) (Bachelor thesis). Pontificia Universidad Católica del Ecuador, Quito, Ecuador.

Garnier-Géré, P., & Chikhi, L. (2013). Population subdivision, Hardy-Weinberg equilibrium and the Wahlund effect. *eLS*. https://doi.org/10.1002/9780470015902.a0005446.pub3

Garreaud, R. D., Vuille, M., Compagnucci, R., & Marengo, J. (2009). Present-day South American climate. *Palaeogeography, Palaeoclimatology, Palaeoecology*, 281, 180-195. https://doi.org/10.1016/j.palaeo.2007.10.032

Gentry, A. H. (1982). Neotropical floristic diversity: Phytogeographical connections between Central and South America, Pleistocene climatic fluctuations, or an accident of the Andean orogeny? *Annals of the Missouri Botanical Garden*, 69, 557-593. https://doi.org/10.2307/2399084

Gentry, A. H. (1988). Tree species richness of upper Amazonian forests. *Proceedings of the National Academy of Sciences*, 85, 156-159. https://doi.org/10.1073/pnas.85.1.156
Gibbs, A. K., & Barron, C. N. (1993). *The geology of the Guiana Shield*. Oxford, UK: Oxford University Press.

Goudet, J. (2001). *FSTAT, a program to estimate and test gene diversities and fixation indices* version 2.9.3. https://www2.unil.ch/popgen/softwares/fstat.htm

Gregory-Wodzicki, K. M. (2000). Uplift history of the Central and Northern Andes: A review. *Geological Society of America Bulletin*, **112**(7), 1091–1105. https://doi.org/10.1130/0016-7606(2000)112&It;1091:UHOTCA&gt;gt;2.0.CO;2

Guillot, G., Mortier, F., & Estoup, A. (2005). *Geneland: A computer package for landscape genetics*. *Molecular Ecology Notes*, **5**, 712–715. https://doi.org/10.1111/j.1471-8286.2005.01031.x

Haffer, J. (1967). Speciation in Colombian forest birds west of the Andes. *American Museum Novitates*, **2294**, 1–57.

Haffer, J. (1969). Speciation in Amazonian forest birds. *Science*, **165**(3889), 131–137. https://doi.org/10.1126/science.165.3889.131

Hamrick, J. L., & Godt, M. J. W. (1990). Allozyme diversity in plant species. In A. H. D. Brown, M. T. Clegg, A. L. Kahler, & B. S. Weir (Eds.), *Plant population genetics, breeding and genetic resources* (pp. 43–63). Sunderland, MA: Sinauer.

Hardesty, B. D., Dick, C. W., Hamrick, J. L., Degen, B., Hubbell, S. P., & Bermingham, E. (2010). Geographic Influence on Genetic Structure in the Widespread Neotropical Tree *Simarouba amara* (*Simaroubaceae*). *Tropical Plant Biology*, **3**, 28–39. https://doi.org/10.1007/s12042-010-9044-3

Henderson, A. (1995). *The palms of the Amazon*. Oxford, UK: Oxford University Press.

Higgins, M. A., Ruokolainen, K., Tuomisto, H., Llerena, N., Cardenas, G., Phillips, O. L., ... Räsänen, M. (2011). Geological control of floristic composition in Amazonian forests. *Journal of Biogeography*, **38**(11), 2136–2149. https://doi.org/10.1111/j.1365-2699.2011.02585.x

Hoorn, C., Wesselingh, F. P., Hovikoski, J., & Guerrero, J. (2010). The development of the Amazonian mega-wetland (Miocene; Brazil, Colombia, Peru, Bolivia). In C. Hoorn & F. P. Wesselingh (Eds.), *Amazonia, landscape and species evolution: A look into the past* (pp. 123–142). Hoboken, NJ: Blackwell Publishing.

Hoorn, C., Wesselingh, F. P., Ter Steege, H., Bermudez, M. A., Mora, A., Sevink, J., ... Antonelli, A. (2010). Amazonia through time: Andean uplift, climate change, landscape evolution and biodiversity. *Science*, **330**, 927–931. https://doi.org/10.1126/science.1194585

Hubert, N., Duponchelle, F., Nunez, J., García-Dávila, C., Paugy, D., & Renno, J. F. (2007). Phylogeography of the piranha genera *Serrasalmus* and *Pygocentrus*: Implications for the diversification of the Neotropical ichthyofauna. *Molecular Ecology*, **16**(10), 2115–2136. https://doi.org/10.1111/j.1365-294X.2007.03267.x

Hughes, C., & Eastwood, R. (2006). Island radiation on a continental scale: Exceptional rates of plant diversification after uplift of the Andes. *Proceedings of the National Academy of Sciences*, **103**(27), 10334–10339. https://doi.org/10.1073/pnas.0601928103

Kahn, F., & de Granville, J. J. (1992). *Palms in forest ecosystems of Amazonia*. Berlin, Germany: Springer-Verlag. https://doi.org/10.1007/978-3-642-76852-1

Karubian, J., Ottewell, K., Link, A., & Di Fiore, A. (2015). Genetic consequences of seed dispersal to sleeping trees by white-bellied spider monkeys. *Acta Oecologica*, **68**, 50–58. https://doi.org/10.1016/j.actao.2015.07.005

Karubian, J., Sork, V. L., Roorda, T., Duraes, R., & Smith, T. (2010). Destination-based seed dispersal homogenizes genetic structure of a tropical palm. *Molecular Ecology*, **19**, 1745–1753. https://doi.org/10.1111/j.1365-294X.2010.04600.x

Klimova, A., Hoffman, J. I., Gutierrez-Rivera, J. N., Leon de la Luz, J., & Ortega-Rubio, A. (2017). Molecular genetic analysis of two native desert palm genera, *Washingtonia* and *Brahea*, from the Baja California Peninsula and Guadalupe Island. *Ecology and Evolution*, **7**, 4919–4935. https://doi.org/10.1002/ece3.3036

Kristiansen, T., Svenning, J. C., Pedersen, D., Eiserhardt, W. L., Grández, C., & Balslev, H. (2011). Local and regional palm (Arecaceae) species richness patterns and their cross-scale determinants in the western Amazon. *Journal of Ecology*, **99**, 1001–1015. https://doi.org/10.1111/j.1365-2745.2011.01834.x

Kumar, S., Stecher, G., & Tamura, K. (2016). MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. *Molecular Biology and Evolution*, **37**(7), 1870–1874. https://doi.org/10.1093/molbev/msw054

Lankau, R. A., & Strauss, S. Y. (2007). Mutual feedbacks maintain both genetic and species diversity in a plant community. *Science*, **317**(5844), 1561–1563. https://doi.org/10.1126/science.1147455

Lepsch-Cunha, N., Lund, C. A., & Hamilton, M. B. (2003). Isolation and characterization of nuclear microsatellite loci in the tropical arboreal palm *Oenocarpus bacaba* (Arecaceae). *Molecular Ecology Notes*, **3**, 435–437. https://doi.org/10.1046/j.1471-8286.2003.00473.x

Lovette, I. J. (2004). Molecular phylogeny and plumage signal evolution in a trans Andean and circum Amazonian avian species complex. *Molecular Phylogenetics and Evolution*, **32**, 512–523. https://doi.org/10.1016/j.ympev.2004.01.007

Ludeña, B., Chabrillange, N., Aberlenc-Bertossi, F., Adam, H., Tregear, J., & Pintaud, J.-C. (2011). Phylogenetic utility of the nuclear genes *AGAMOUS 1* and *PHYTOCHROME B* in palms (Arecaceae): An example within Bactridinae. *Annals of Botany*, **108**, 1433–1444. https://doi.org/10.1093/aob/mcr191

Luebert, F., & Weigend, M. (2014). Phylogenetic insights into Andean plant diversification. *Frontiers in Ecology and Evolution*, **2**, 27.

Martin, N. H., Bouck, A. C., & Arnold, M. L. (2007). The genetic architecture of reproductive isolation in Louisiana irises: Flowering phenology. *Genetics*, **175**(4), 1803–1812. https://doi.org/10.1534/genetics.106.068338

Martius, K. P. F. von (1823). *Historia naturalis palmarum*. Leipzig, Germany: T. O. Weigel.

Meerow, A. W., Noblick, L., Salas-Leiva, D. E., Sanchez, V., Francisco-Ortega, J., Jestrow, B., & Nakamura, K. (2015). Phylogeny and historical biogeography of the cocosoid palms (Arecaceae, Arcoideae, Cocoseae) inferred from sequences of six WRKY gene family loci. *Cladistics*, **31**(5), 509–534. https://doi.org/10.1111/cla.12100

Montúfar, R. (2007). *Structure génétique, biochimique, morphologique et écologique de Oenocarpus bataua* Mart. (Arecaceae): *perspectives pour la valorisation durable d'une ressource forestière néotropicale* (PhD thesis). Université de Montpellier, Montpellier, France.

Montúfar, R., Laffargue, A., Pintaud, J.-C., Hamon, S., Avallone, S., & Dussert, S. (2010). *Oenocarpus bataua* Mart. (Arecaceae): Rediscovering a Source of High Oleic Vegetable Oil from Amazonia. *Journal of the American Oil Chemists' Society*, **87**(2), 167–172. https://doi.org/10.1007/s11746-009-1490-4

Montúfar, R., Mariac, C., Pham, J., & Pintaud, J.-C. (2006). Isolation of 23 polymorphic microsatellite loci in the Neotropical palm *Oenocarpus bataua* Martius (Arecaceae). *Molecular Ecology Notes*, **7**(1), 75–78. https://doi.org/10.1111/j.1471-8286.2006.01532.x

Montúfar, R., & Pintaud, J.-C. (2006). Variation in species composition, abundance and microhabitat preferences among western Amazonian terra firme palm communities. *Botanical Journal of the Linnean Society*, **151**, 127–140. https://doi.org/10.1111/j.1095-8339.2006.00528.x

Montúfar, R., & Pintaud, J.-C. (2008). Estatus taxonómico de *Oenocarpus bataua* (Euterpeae, Arecaceae) inferido por secuencias del ADN cloroplástico. *Revista Peruana de Biología*, **15**, 73–78.

Mora, A., Baby, P., Roddaz, M., Parra, M., Brusset, S., Hermoza, W., & Espurt, N. (2010). Tectonic history of the Andes and sub-Andean zones: Implications for the development of the Amazon drainage basin. In C. Hoorn & F. P. Wesselingh (Eds.), *Amazonia, landscape and species evolution: A look into the past* (pp. 38–60). Hoboken, NJ: Blackwell Publishing.

Motamayor, J. C., Lachenaud, P., da Silva e Mota, J. W., Loor, R., Kuhn, D. N., Brown, J. S., & Schnell, R. J. (2008). Geographic and
genetic population differentiation of the Amazonian chocolate tree (*Theobroma cacao* L.). *PLoS One*, 3(10), e33311. https://doi.org/10.1371/journal.pone.0003311

Nelson, B. W., Ferreira, C. A., da Silva, M. F., & Kawasaki, M. L. (1990). Endemism centres, refugia and botanical collection density in Brazilian Amazonia. *Nature*, 345(6277), 714. https://doi.org/10.1038/345714a0

Núñez-Avellaneda, L. A. (2007). *Evaluación del aislamiento reproductivo entre tres especies simpátricas* Oenocarpus bataua, Oenocarpus mapora y Euterpe precatoria (*Palmae:Arecoideae: Euterpeae*) en un bosque de los Andes en Colombia (Master thesis). Universidad Nacional de Colombia, Bogotá, Colombia.

Núñez-Avellaneda, L. A., & Rojas-Robles, R. (2008). Biología reproduc-
tiva y ecología de la polinización de la palma milpesos Oenocarpus
bataua en los andes colombianos. *Caldasia*, 30(1), 101-125.

Olson, D. M., & Dinerstein, E. (2002). The Global 200: Priority ecoregions for global conservation. *Annals of the Missouri Botanical Garden*, *89*, 199-224. https://doi.org/10.2307/3298564

Palma-Silva, C., Lerer, C., Paggi, G. M., Barbará, T., Bered, F., & Bodanese-Zanettini, M. H. (2009). Range-wide patterns of nuclear and chloroplast DNA diversity in *Vriesea gigantea* (Bromeliaceae), a neotropical forest species. *Heredity*, 103, 503–512. https://doi.org/10.1038/hdy.2009.116

Peralta Rivero, C. A., Zonta, M., Moraes, V., & Ríos, R. (2008). *Efecto del aprovechamiento comercial de frutos de majo* (Oenocarpus bataua C. Martius) en estructura y densidad poblacional, en tres comunidades en el Norte Amazónico de Bolivia (Bachelor thesis). Universidad Autónoma del Beni, Trinidad, Bolivia.

Pérez-Escobar, O. A., Gottschling, M., Chomicki, G., Condamine, F. L., Klitgård, B. B., Pansarin, E., & Gerlach, G. (2017). Andean mountain building did not preclude dispersal of lowland epiphytic orchids in the Neotropics. *Scientific Reports*, 7, 4919. https://doi.org/10.1038/s41598-017-04261-z

Pitman, N. C. A. (2000). *A large-scale inventory of two Amazonian tree communities* (PhD thesis). Duke University, Durham, NC.

Pritchard, J. K., Stephens, M., & Donnelly, P. (2000). Inference of pop-
ulation structure using multilocus genotype data. *Genetics*, *155*,
945-959.

Quintana, C., Peninngton, R. T., Ulloa Ulloa, C., & Balslev, H. (2017). Biogeographic barriers in the Andes: Is the Amotape-Huancabamba zone a dispersal barrier for dry forest plants? *Annals of the Missouri Botanical Garden*, 102(3), 542–550. https://doi.org/10.3417/D-17-00003A

Raymond, M., & Rousset, F. (1995). *GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism*. *Journal of Heredity*, 86, 248–249. https://doi.org/10.1093/oxfordjournals.jhered.a111573

Rhoads, E. A., Williams, P. K., & Krane, C. M. (2017). High inbreeding and low connectivity among *Ambystoma texanum* populations in fragmented Ohio forests. *Ecology and Evolution*, 7, 11135-11147. https://doi.org/10.1002/ece3.3637

Roddaz, M., Viers, J., Brusset, S., Baby, P., & Hérail, G. (2005). Sediment provenances and drainage evolution of the Neogene Amazonian foreland basin. *Earth and Planetary Science Letters*, 239, 57-78. https://doi.org/10.1016/j.epsl.2005.08.007

Rymer, P. D., Dick, C. W., Vendramin, G. G., Buonamici, A., & Boshier, D. (2012). Recent phylogeographic structure in a widespread 'weedy' Neotropical tree species, *Cordia alliodora* (Boraginaceae). *Journal of Biogeography*, 40, 693-706.

Sereno, M. L., Albuquerque, P. S. B., Vencovsky, R., & Figueira, A. (2006). Genetic diversity and natural population structure of cacao (*Theobroma cacao* L.) from the Brazilian Amazon evaluated by microsatellite markers. *Conservation Genetics*, 7, 13–24. https://doi.org/10.1007/s10592-005-7568-0

Silman, M. R. (2007). Plant species diversity in Amazonian forests. In M. B. Bush & J. R. Flenley (Eds.), *Tropical rainforest responses to climatic*

change (pp. 269-294). Chichester, UK: Praxis Publishing Ltd. https://doi.org/10.1007/978-3-540-48842-2

Struwe, L., Haag, S., Heiberg, E., & Grant, J. R. (2009). Andean speciation and vicariance in Neotropical *Macrocarpaea* (Gentianaceae-Helieae). *Annals of the Missouri Botanical Garden*, 96, 450–469. https://doi.org/10.3417/2008040

Ter Steege, H., Pitman, N. C. A., Sabatier, D., Baraloto, C., Salomão, R. P., Guevara, J. E., ... Silman, N. (2013). Hyperdominance in the Amazonian Tree Flora. *Science*, 342(6156), 1243092. https://doi.org/10.1126/science.1243092

Thomas, E., Van Zonneveld, M., Loo, J., Hodgkin, T., Galluzzi, G., & van Etten, J. (2012). Present spatial diversity patterns of *Theobroma cacao* L. in the Neotropics reflect genetic differentiation in Pleistocene refugia followed by human-influenced dispersal. *PLoS One*, 7(10), e47676. https://doi.org/10.1371/journal.pone.0047676

Trénel, P., Hansen, M. M., Normand, S., & Borchsenius, F. (2008). Landscape genetics, historical isolation and cross-Andean gene flow in the wax palm, *Ceroxylon echinulatum* (Arecaceae). *Molecular Ecology*, 17(15), 3528–3540.

Tuomisto, H., Zuquim, G., & Cárdenas, G. (2014). Species richness and diversity along edaphic and climatic gradients in Amazonia. *Ecography*, 37(11), 1034–1046.

Turchetto-Zolet, A. C., Pinheiro, F., Salgueiro, F., & Palma-Silva, C. (2012). Phylogeographical patterns shed light on evolutionary process in South America. *Molecular Ecology*, 22, 1193–1213.

Valencia, R., Balslev, H., & Paz y Miño, G. C. (1994). High tree alpha-
diversity in Amazonian Ecuador. *Biodiversity and Conservation*, 3,
21-28. https://doi.org/10.1007/BF00115330

Vellend, M., & Geber, M. A. (2005). Connections between species diversity and genetic diversity. *Ecology Letters*, 8, 767–781. https://doi.org/10.1111/j.1461-0248.2005.00775.x

Vormisto, J., Svenning, J. C., Hall, P., & Balslev, H. (2004). Diversity and dominance in palm (Arecaceae) communities in *terra firme* forests in the western Amazon basin. *Journal of Ecology*, 92, 577–588. https://doi.org/10.1111/j.0022-0477.2004.00904.x

Weigend, M. (2002). Observations on the Biogeography of the Amotape-Huancabamba Zone in Northern Peru. *The Botanical Review*, 68(1), 38-54. https://doi.org/10.1663/0006-8101(2002)068[0038:OOTBOT]2.0.CO;2

Weigend, M. (2004). Additional observations on the biogeography of the Amotape-Huancabamba zone in Northern Peru: Defining the South-Eastern limits. *Revista Peruana de Biologia*, 11(2), 127-134.

Welt, R. S., Litt, A., & Franks, S. J. (2015). Analysis of population genetic structure and gene flow in an annual plant before and after a rapid evolutionary response to drought. *AoB PLANTS*, 7, plv026. https://doi.org/10.1093/aobpla/plv026

Zartman, C. E., McDaniel, S. F., & Shaw, A. J. (2006). Experimental habitat fragmentation increases linkage disequilibrium but does not affect genetic diversity or population structure in the Amazonian liverwort *Radula flaccida*. *Molecular Ecology*, 15(9), 2305–2315. https://doi.org/10.1111/j.1365-294X.2006.02929.x

SUPPORTING INFORMATION

Additional supporting information may be found online in the
Supporting Information section at the end of the article.

How to cite this article: Escobar S, Pintaud J-C, Balslev H,
et al. Genetic structuring in a Neotropical palm analyzed
through an Andean orogenesis-scenario. *Ecol Evol.*
2018;00:1-13. https://doi.org/10.1002/ece3.4216
**Graphical Abstract**

The contents of this page will be used as part of the graphical abstract of html only.

It will not be published as part of main article.

<figure><img src="image_12.png" /></figure>

Andean orogeny explains the structuring of genetic diversity in the palm *Oenocarpus bataua* into four genetic clusters within northwestern South America. These clusters maintain genetic exchange between them, even across the Andes. A hotspot of genetic diversity was determined in the northern Peruvian Amazonia. Furthermore, we confirm the divergence of the intraspecific varieties *bataua* and *oligocarpus*.