# Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

André Carneiro Muniz, Ricardo José Gonzaga Pimenta, Mariana Vargas Cruz, Jacqueline Gomes Rodrigues, Renata Santiago de Oliveira Buzatti, Myriam Heuertz, José Lemos-filho, Maria Bernadete Lovato

► To cite this version:

André Carneiro Muniz, Ricardo José Gonzaga Pimenta, Mariana Vargas Cruz, Jacqueline Gomes Rodrigues, Renata Santiago de Oliveira Buzatti, et al. Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation. Ecology and Evolution, 2022, 12 (1), pp.1-20. 10.1002/ece3.8540. hal-03603921

HAL Id: hal-03603921

https://hal.inrae.fr/hal-03603921

Submitted on 10 Mar 2022

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
RESEARCH ARTICLE

Ecology and Evolution
Open Access

WILEY

# Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

André Carneiro Muniz<fnref n="1" /> | Ricardo José Gonzaga Pimenta<fnref n="1" /> | Mariana Vargas Cruz<fnref n="1" /> ID |

Jacqueline Gomes Rodrigues<fnref n="1" /> | Renata Santiago de Oliveira Buzatti<fnref n="1" /> | Myriam Heuertz<fnref n="2" /> |

José P. Lemos-Filho<fnref n="3" /> | Maria Bernadete Lovato<fnref n="1" /> ID

## Correspondence

Maria Bernadete Lovato, Departamento de Genética, Ecologia e Evolução, Universidade Federal de Minas Gerais, CP 486, Belo Horizonte, MG, 31270-901, Brazil.

Email: lovatomb@icb.ufmg.br

## Funding information

Investissement d'Avenir grant of the ANR, Grant/Award Number: ANR-10-LABEX-0025; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior; Fundação de Amparo a Pesquisa do Estado de Minas Gerais; Conselho Nacional de Desenvolvimento Científico e Tecnológico

## Abstract

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree *Plathymenia reticulata* (n = 233 individuals). Genetic structure analysis of *P. reticulata* was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F<sub>CT</sub> = 0.222 and F<sub>ST</sub> = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F<sup>2</sup>) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in *P. reticulata*. The low genetic diversity of the *P. reticulata* forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

## KEYWORDS

conservation, ecotypes, hybrid zone, Neotropical tree, *Plathymenia*, savanna-forest ecotone

## TAXONOMY CLASSIFICATION

Conservation genetics

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2022 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.

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

<footnote n="1">¹Departamento de Genética, Ecologia e Evolução, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil</footnote>

<footnote n="2">²Biogeco, INRAE, Univ. Bordeaux, Cestas, France</footnote>

<footnote n="3">³Departamento de Botânica, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil</footnote>
# 1 | INTRODUCTION

Understanding the processes that underlie genetic divergence between lineages is one of the major goals of evolutionary biology. The speciation process occurs along a continuum of genetic divergence where the incipient lineages often undergo a period of partial reproductive isolation until they become fully isolated (Christie & Strauss, 2018; Dufresnes et al., 2019; Kopp & Frank, 2005; Wu, 2001). It is frequently recognized that most genetic divergence during speciation occurs owing to allopatric divergence, where there is little or no gene flow between geographically isolated populations (Coyne, 1992). However, a growing body of evidence suggests that speciation with gene flow is more common in nature than previously thought (Choi et al., 2019; Feder et al., 2012; Hanson et al., 2016; Rifkin et al., 2019). The genomic divergence under speciation with gene flow can be very different from that of allopatric divergence because gene flow can homogenize the genomic background, while islands of divergence may result for loci under divergent natural selection (Feder et al., 2012). In speciation with gene flow models, divergence between incipient lineages is mainly driven by two processes, divergent selection that favors the divergence of groups of individuals with extreme traits, and nonrandom or assortative mating (Servedio, 2016).

Differential selection and nonrandom mating may lead to an isolation by environment (IBE) pattern, in which the association between genetic and environmental distance is stronger than that between genetic and spatial distance between populations (Wang & Bradburd, 2014). IBE is agnostic to the cause of divergence and rather a description of the pattern since many processes such as selection against immigrants, reduced hybrid fitness, and biased effective dispersal can generate a pattern of IBE (Wang & Bradburd, 2014). Under IBE, even neutral markers are expected to show some level of correlation with the environmental distance due to the impact of differential gene flow among ecologically similar or dissimilar habitats (Nosil et al., 2008, 2009). Disentangling IBE from isolation by distance (IBD), where lineages diverge due to limited dispersal of individuals, may be statistically difficult since the environmental distance and the spatial distance are often correlated (Nadeau et al., 2016). However, evaluating the relative importance of IBE versus IBD in the genetic divergence can give some insights into the possible causes of evolutionary divergence of lineages.

In the speciation process, some degree of hybridization may be commonly found. Hybridization, estimated to occur in 25% of plant taxa (Mallet, 2005), can create new gene combinations in later generation hybrids, lead to new species through hybrid speciation (Abbott, 2017; Abbott et al., 2013; Maguilla & Escudero, 2016; Paun et al., 2009), and influence the patterns of genetic variability of species through genomic introgression (Francisco et al., 2018; Gompert et al., 2017; Riesberg & Carney, 1998). Furthermore, hybridization can lead to the disruption of reproductive barriers between incipient lineages reverting the speciation process (Abbott et al., 2013; Soltis & Soltis, 2009) or can lead to extinction through genetic or demographic swamping (Rhymer & Simberloff, 1996; Todesco et al.,

2016). Therefore, the study of gene flow between lineages and the genotypic composition of hybrid zones can help us investigate evolutionary questions related to the formation of species and the development of reproductive barriers between lineages.

Neotropical ecosystems harbor high levels of biodiversity and species endemism, but the causes of such remarkable diversification are still poorly understood. The Brazilian Atlantic Forest, constituted by evergreen and semideciduous forest (Oliveira-Filho & Fontes, 2000), and the Cerrado, considered to be the most species-rich savanna of the world, are recognized hotspots for conservation (Myers et al., 2000). The Cerrado features poorer and more acid soils and its climate is characterized by a long dry season with less annual precipitation in comparison with the Atlantic Forest (Eiten, 1972; Klink & Machado, 2005). In contrast with the Atlantic Forest, natural fires are common in the Cerrado (Hoffmann et al., 2003).

The Cerrado and Atlantic Forest are in contact in southeastern Brazil where they form large ecotonal areas where savanna and forest habitats occur in adjacent patches with differences in tree density and species composition (Durigan & Ratter, 2006). In these ecotonal areas, pairs of closely related species from the Cerrado and Atlantic Forest are found side by side (Hoffmann & Franco, 2003; Hoffmann et al., 2003). This spatial proximity can provide opportunities for gene flow between closely related species if reproductive isolation between these taxa is not completed. In fact, a few studies using nuclear markers with sampling in contact zones between the Atlantic Forest and the Cerrado have found gene flow between closely related tree taxa (Cavallari et al., 2010; Lacerda et al., 2002; Muniz et al., 2020; Resende-Moreira et al., 2017).

*Plathymenia reticulata* Benth (Fabaceae, Caesalpinioideae - LPWG, 2017) is the single species of the *Plathymenia* genus (Warwick & Lewis, 2003) and features a widespread distribution in both biomes, the Cerrado and the Atlantic Forest. The species presents two ecotypes, one associated with forest habitats and the other with savanna habitats (Lemos-Filho et al., 2008). The *P. reticulata* forest ecotype usually shows straight trunks reaching up to 30 m in height and up to 70 cm in diameter, whereas the savanna ecotype reaches 6-12 m in height and 30-50 cm in trunk diameter and displays a shorter and tortuously branched trunk (Lorenzi, 1992, Figure S1). The ecotypes also show differences in ecophysiological traits, which are associated to adaptive responses between savanna and forest habitats (Lemos-Filho et al., 2008), such as in vegetative phenology, seed and fruit morphology, seed dormancy and germination, phenotypic plasticity in response to light availability, and stem radial growth (Goulart, Lemos-Filho, & Lovato, 2005, 2006; Goulart et al., 2011; Toledo et al., 2012). In addition, these studies showed that the populations occurring in the ecotone between the Atlantic Forest and the Cerrado are generally morphophysologically intermediate relative to populations from the Atlantic Forest and the Cerrado.

A common garden experiment with *P. reticulata* progenies of populations from the Cerrado, the Atlantic Forest, and ecotonal areas showed that individuals developed the morphological and physiological characteristics of their respective habitats (Goulart et al., 2011), indicating that their differences are genetically determined. Genetic
studies of *Plathymenia* showed admixture in an ecotonal population suggesting gene flow between the two ecotypes (Lacerda et al., 2002) and a pattern of plastid DNA variation congruent with a diverse, persistent central-north Cerrado population as a potential source of population expansion after the Last Glacial Maximum (Novaes et al., 2010). Therefore, *P. reticulata* likely shows a complex pattern of population genetic structure influenced by range dynamics due to Quaternary climate oscillations, the divergence of ecotypes, and admixture in areas where both ecotypes come into contact.

In this study, we sampled *P. reticulata* in the Cerrado and the Atlantic Forest, and in ecotonal areas between the two biomes. We used nuclear microsatellite markers to investigate the patterns and the putative drivers of genetic divergence between savanna and forest ecotypes and to characterize a putative hybrid zone in the ecotonal areas. Specifically, we addressed the following questions: (1) What level of genetic diversity characterizes *P. reticulata* and how is genetic diversity structured geographically, between and within ecotypes? (2) How much gene flow occurs in ecotonal populations, and is gene flow asymmetric between ecotypes? (3) What is the genotypic composition of the hybrid populations, are they mostly F1s or later generation hybrids? (4) What is the effect of environmental variation on the genetic structure? We discuss the evolutionary consequences of divergence in *P. reticulata* and the implications for conservation of genetic diversity of the ecotonal populations.

## 2 | MATERIALS AND METHODS

### 2.1 | Sampling strategy

We sampled 233 individuals of *P. reticulata*, comprising the savanna and forest ecotypes in 10 localities from three different habitats with three localities sampled in the Atlantic Forest (SJF, AJF, and IPF) and the Cerrado (PRS, PTS, and VZS), and four localities sampled in the ecotone between the two biomes (NEE, COE, SUE, and FEE) (Figure 1). In the SUE and FEE localities of the ecotone, both savanna and forest ecotypes occur in sympatry/parapatry, whereas in COE and NEE, only the savanna or forest ecotypes were found, respectively. In the ecotonal localities, we used the letters s (savanna) or f (forest) in the locality code to identify the ecotype we were referring to. For example, in the SUE locality, we sampled s-SUE and f-SUE populations of each ecotype. The classification of individuals in ecotypes was based on differences in size (savanna individuals did not exceed 10 m in height of the canopy, whereas forest individuals measured from 15 to 30 m in height), trunk (savanna individuals had a tortuous and twisted trunk, whereas forest individuals had a straight trunk), and bark (which was more suberous in the savanna ecotype) (Figure S1).

### 2.2 | DNA isolation, amplification, and genotyping

Leaf or bark tissue samples were dried on silica gel and kept frozen at -20°C. Genomic DNA was extracted from leaves using the

sorbitol buffer-based method according to the protocol described by Souza et al. (2012). For bark tissue, we used the protocol developed by Novaes et al. (2009). We genotyped individuals for 11 microsatellite markers previously developed for *P. reticulata*: Pre05, Pre10, Pre11, Pre14, Pre15, Pre16, and Pre23 (Cruz et al., 2012) and Pr18, Pr30, Pr71, and Pr80 (Oliveira et al., 2012). Polymerase chain reactions (PCRs) were performed following the conditions established by Oliveira et al. (2012). Genotyping was performed using the automated sequencers MegaBACE 1000 (GE Healthcare, Chicago, Illinois, EUA) or ABI 3730xl DNA Analyzer (Applied Biosystems, Foster City, California, EUA) through capillary electrophoresis. Allele size was determined by comparison with either GeneScan ROX-500 or LIZ-500 Size Standards (Applied Biosystems) using MegaBACE Fragment Profiler version 1.2 (GE Healthcare) and Geneious version 7.1.3 (Biomatters Ltd., Auckland, New Zealand). The presence of null alleles at each locus was examined by ecotype by population with 1000 Monte Carlo simulations using Micro-Checker version 2.3.3 (Van Oosterhout et al., 2004) and the linkage equilibrium among pairs of loci was evaluated using FSTAT version 2.9.3.2 (Goudet, 2002). Micro-Checker analysis found significant evidence for the presence of null alleles in all *P. reticulata* localities, except for SJF and s-SUE (Table S1). However, the loci did not show significant null alleles regularly across localities, which can indicate that these loci are not in Hardy-Weinberg equilibrium in these populations; therefore, all subsequent analyses were performed with the full, noncorrected dataset. No significant linkage disequilibrium was found among loci in any pairwise comparison at the population level.

### 2.3 | Genetic diversity and structure of *Plathymenia reticulata*, and gene flow between ecotypes

Arlequin version 3.5 (Excoffier & Lischer, 2010) was used to estimate the mean number of alleles per locus (A) and the mean observed ($H_O$) and expected ($H_E$) heterozygosities. The allelic richness ($A_R$) and the inbreeding coefficient ($F IS$) were estimated in FSTAT version 2.9.3.2 (Goudet, 2002). The hierarchical partitioning of genetic variance was investigated using the standard analysis of molecular variance (AMOVA) according to the method of Weir and Cockerham (1984), implemented in Arlequin version 3.5 (Excoffier & Lischer, 2010). We considered variation between ecotypes as the highest hierarchical level of genetic divergence, followed by variation among populations within ecotypes and variation within populations, and computed the corresponding *F* statistics. The pairwise genetic differentiation between localities was evaluated by $F_{ST}$ estimated in Arlequin version 3.5 (Excoffier & Lischer, 2010).

We used POPTREEW (Takezaki et al., 2014) to build a tree from allele frequency data using the neighbor-joining (NJ) method (Saitou & Nei, 1987) based on Nei's standard genetic distance ($D_{ST}$) (Nei, 1972) with sample size bias correction. We used the *P. reticulata* populations without the admixed individuals in order to evaluate the phylogenetic relationships among populations as well as between ecotypes. Bootstrap tests (Felsenstein, 1985) were performed for
<figure><img src="image_3.png" /><figcaption>FIGURE 1 Map showing the sampling localities of *Plathymenia reticulata* in the Brazilian Atlantic Forest (dark gray) and the Cerrado (light gray) biomes, and in the ecotone between the two</figcaption></figure>

the evaluation of confidence intervals of the tree's branches in order to estimate the robustness of the tree.

We investigated the genetic structure as well as estimated admixture in *P. reticulata* using the following algorithms: (1) discriminant analysis of principal components (DAPC) implemented in the *adegenet* package version 2.1.4 in R 4.1.0 (Jombart, 2008; Jombart et al., 2010; R Core Team, 2020) and (2) the Bayesian clustering method implemented in STRUCTURE version 2.3.4 (Hubisz et al., 2009; Pritchard et al., 2000). DAPC has the advantage of not relying on assumptions of Hardy-Weinberg equilibrium or linkage equilibrium. In the DAPC, we used the sequential *k*-means clustering method implemented in the *adegenet* package (Jombart et al., 2010) to choose the number of genetic groups that best explains the genetic variability in *P. reticulata* and then used these groups in DAPC to estimate the posterior assignment probability of individuals to groups. After the evaluation of the population structure based on the two ecotypes, we applied the DAPC algorithm in each ecotype separately using only the individuals that showed a posterior assignment probability higher than 0.85, which we considered to be pure individuals. In STRUCTURE version 2.3.4, we set the number of clusters from *K* = 1 to *K* = 13 and ran 10 independent iterations per *K* assuming an admixture

model with correlated allelic frequencies between localities. We ran 100,000 steps as burn-in followed by 1,000,000 Markov chain Monte Carlo (MCMC) iterations. The most probable number of genetic clusters was determined by the ad hoc statistic ΔK (Evanno et al., 2005), computed in Structure Harvester version 0.6.94 (Earl & vonHoldt, 2012). Besides the Evanno's ΔK method, we also used "MedMeaK" (median of means), "MaxMeaK" (maximum of means), "MedMedK" (median of medians), and "MaxMedK" (maximum of medians) statistics in StructureSelector (Li & Liu, 2018) based on the method of Puechmaille (2016). In this method, a subpopulation is assigned to a cluster if its arithmetic mean (for MedMeaK and MaxMeaK) or its median (for MedMedK and MaxMedK) membership coefficient to that cluster is greater than a threshold value (set to 0.5), therefore ensuring that a subpopulation cannot belong to more than one cluster. We averaged the admixture coefficients using CLUMPP version 1.1.2 (Jakobsson & Rosenberg, 2007) and plotted the values using the ggplot2 package version 3.3.5 in the R software (Wickam, 2009). We also tested for genetic structure within each ecotype, using all individuals identified as pure (<15% or >85% ancestry threshold) in the STRUCTURE analysis of the full dataset. For each ecotype, we ran a separated analysis setting *K* = 1 to *K* = 7.
To estimate rates of recent migration between the two ecotypes (based on phenotype classification), we used BayesAss 3.0 (Wilson & Rannala, 2003). This software uses a Bayesian framework and MCMC iterations to estimate the posterior mean of recent migration rates (mainly the last two generations) and the standard deviation of the marginal posterior distribution from multilocus genotypes for a group of populations in a sample. The software provides estimates using the gametic disequilibrium among immigrants and their descendants. This method relaxes the assumption of previous Hardy-Weinberg equilibrium within populations to assign individuals to populations and identify migrants. We estimated the migration rates between ecotypes in the following datasets: (1) the whole sample of individuals of each ecotype, (2) the ecotonal individuals of the two ecotypes, (3) FEE locality, (4) SUE locality, and (5) forest ecotype of NEE and savanna ecotype of COE, as each of these localities has only one ecotype, but they are spatially close to each other.

## 2.4 | Genotypic composition of *Plathymenia* hybrids in the ecotone

We investigated the genotypic class composition of *P. reticulata* individuals using NewHybrids version 1.1 (Anderson & Thompson, 2002). The implemented method calculates the posterior probability of individuals to belong to predetermined genotypic classes (Anderson & Thompson, 2002). We applied the method using two different sets of genotypic classes. First, we used the six default predetermined genotypic classes of NewHybrids, two purebred parental classes P1 (savanna) and P2 (forest), F1 (savanna × forest) and F2 (F1 × F1) hybrids, and one backcross of F1 with each parental class (B1 = F1 × savanna and B2 = F1 × forest). Secondly, we investigated more advanced generation hybrids using the following additional six classes of backcrosses: B3 = F2 × savanna and B4 = F2 × forest, B5 = savanna × (F1 × savanna) and B6 = forest × (F1 × forest), B7 = savanna × (F2 × savanna) and B8 = forest × (F2 × forest), according to the genotypic classes in Chhatre et al. (2018) and Buck et al. (2020) (Table S2). In both analyses, we ran 1,000,000 MCMC iterations with a burn-in period of 100,000 iterations using Jeffrey's-like priors for mixing proportions and allele frequencies and assuming equal prior probabilities for all genotypic classes.

We further investigated the genotypic composition of *P. reticulata* individuals using the hybrid index estimated in the *introgress* package version 1.2.3 in R (Gompert & Buerkle, 2010). The computation of the hybrid index is based on the maximum likelihood method and estimates the genetic contribution of hybridizing "parental" populations to individuals of unknown ancestry (Buerkle, 2005). Here, we used as pure parental individuals those from the Atlantic Forest and Cerrado localities, respectively, to estimate the hybrid index of those individuals from the ecotone localities. In addition, we calculated the interspecific (interecotypic) heterozygosity ($Q_{12}$), which is the fraction of loci that combine the ancestry of the two parental taxa (Fitzpatrick & Irwin, 2012). The combination of $Q_{12}$ with the hybrid index can give insights into the history of hybridizing

taxa through the summary of the genomic composition of individuals (Gompert & Buerkle, 2016). For example, an F2 individual may show a hybrid index of 0.5 just like F1s. However, F1s show $Q_{12}$ = 1.0, whereas F2 and later generation hybrids have lower values of $Q_{12}$ owing to the breakup of linkage disequilibria because of genome re-shuffling in later generations of hybridization.

## 2.5 | Evaluation of hybrid assignment in the NewHybrids

The efficiency of NewHybrids to identify and classify hybrids was evaluated using purebred parental and hybrid individuals simulated using the function *hybridize* in the *adegenet* package (Jombart, 2008). We simulated five datasets based on the allele frequencies of pure (Q higher than 0.85 or lower than 0.15) *P. reticulata* individuals of the Atlantic Forest and the Cerrado areas, respectively. The simulated datasets comprised 100 individuals of each parental and 100 individuals of each hybrid class, F1, F2, and backcrosses of F1 with each parental (B1 and B2). We also simulated more advanced generation hybrids including F3 (F2 × F2) and F4 (F3 × F3), $B3 = F2 \times savanna$ and $B4 = F2 \times forest$, $B5 = savanna \times (F1 \times savanna)$ and $B6 = forest \times (F1 \times forest)$, $B7 = savanna \times (F2 \times savanna)$ and $B8 = forest \times (F2 \times forest)$, based on Chhatre et al. (2018) and Buck et al. (2020) (Table S2). To evaluate the performance of NewHybrids to distinguish hybrid classes, we ran the simulated dataset of six genotypic classes using six genotypic classes in NewHybrids, the 14 simulated genotypic classes using six genotypic classes in NewHybrids, and the 14 simulated genotypic classes using 12 genotypic classes in NewHybrids (Table S2). Finally, we set four thresholds for posterior probabilities of assignment to determine genotypic classes (0.5, 0.75, 0.85, and 0.9) to evaluate the proportion of correct and incorrect assignments. We also assessed the efficiency of STRUCTURE and *introgress* to detect and classify hybrids based on simulated data (Appendix S1).

## 2.6 | Drivers of divergence of *Plathymenia* *reticulata* ecotypes

We examined the effect of edaphic and bioclimatic variables on genetic differentiation between *P. reticulata* populations using two approaches: (1) Mantel tests as a first approach to assess effects of isolation by distance and isolation by environment and (2) distance-based redundancy analysis, dbRDA. For Mantel tests, we calculated the pairwise spatial distance between populations using the coordinates of the localities separated by ecotypes with the *geodist package* version 0.0.7 in R (Padgham & Sumner, 2021). Thus, we had six ecotonal populations since the populations of SUE and FEE were separated in f-SUE and s-SUE and in f-FEE and s-FEE, respectively. As environmental predictors, we used the edaphic variables *bulk density* of the fine earth fraction (BDOD), *soil pH* (pH), *soil organic carbon* (SOC), and *cation exchange capacity* of the soil (CEC), which
were downloaded from https://files.isric.org/soilgrids/latest/data/
and the 19 bioclimatic variables from the worldclim2 dataset https://
www.worldclim.org/data/bioclim.html (Fick & Hijmans, 2017). We
downloaded a raster file of edaphic data for the soil layer between
15 and 30 cm depth using a spatial resolution of 250 m. The 19
bioclimatic variables were downloaded using a spatial resolution of
30 m. Both datasets were cropped to the extent of the sampling area
for further analysis. We assessed the correlation between edaphic
variables and bioclimatic variables separately through principal com-
ponent analysis using the function rasterPCA of the RStoolbox ver-
sion 0.2.6 (Leutner & Horning, 2019) in order to avoid collinearity
between variables. As the edaphic variables showed low correlation
(<0.7), we used the original dataset for posterior analysis. Bioclimatic
variables showed high correlation; therefore, we based our analyses
on the first five PCA components which together explained 95.5% of
the variance in the bioclimatic variables (Table S3). We calculated the
Mahalanobis distance of the environmental data using the distance
function of the ecodist 2.0.7 package (Goslee & Urban, 2007).

We performed Mantel tests in the R package *ncf* 1.2-9 (Bjornstad, 2020) using the pairwise *F*<sub>ST</sub> between populations as genetic distance. We performed simple Mantel tests between genetic and spatial distance, or between genetic and environmental distance, separately for edaphic and bioclimatic variables. Then, we performed partial Mantel tests between genetic and environmental distances in which the spatial distance was used as a controlling factor, and between genetic and spatial distances in which environment was used as a controlling factor.

To further evaluate the association between environmental
variables and the genetic divergence among populations, we per-
formed a distance-based redundancy analysis (dbRDA) (Legendre
& Legendre, 2014). This multivariate method combines ordination
with multivariate regression approach to find linear relationships
between a matrix of predictor variables and a matrix of response
variables. The method is relatively robust to collinearity among
predictor variables (Wagner & Fortin, 2012) and in the detection of
environmental gradients in landscapes (Kierepka & Latch, 2015). In
these analyses, we used the environmental dataset as predictors and
the scores of the first three axes of a spatial principal components
analysis (sPCA) obtained using the *spca* function from *adegenet*
package as the genetic variables (Jombart et al., 2008). The sPCA
summarizes the genetic variability of a group of individuals while it
controls for spatial autocorrelation between them using Moran's I
statistic (Jombart et al., 2008). We also used the *pcnm* function of
the *vegan* package version 2.5-7 to calculate spatial predictors for
the dbRDA based on the transformation of latitude and longitude of
populations (Oksanen et al., 2020). To find the best dbRDA model,
we performed a stepwise variable selection based on the Akaike in-
formation criterion (AIC) using *ordiR2step* function in *vegan* (Oksanen
et al., 2020) and used the variance inflation factor (VIF) to further as-
sess the collinearity of the variables within the selected model. The
significance of the models, of the variables in the models, and of the
axis of dbRDA was assessed with a permutation test based on 1000
permutations using the *Anova* function in R (R Core Team, 2020).

Finally, we estimated the proportion of variance explained by each
variable (partial R²) using the function varpart of the vegan package
(Oksanen et al., 2020).

3 | RESULTS

3.1 | Genetic diversity and structure of *Plathymenia reticulata*

Plathymenia reticulata localities showed A and A<sub>R</sub> values ranging
from 3.9 to 10.0 and from 3.72 to 6.61, respectively, and H<sub>O</sub> and
H<sub>E</sub> values ranging from 0.361 to 0.667 and 0.456 to 0.749, respec-
tively (Table 1). Populations of the Cerrado harbored higher diver-
sity than those of the Atlantic Forest for all the diversity indices.
In addition, the ecotone localities showed higher diversity than the
Atlantic Forest and Cerrado localities; however, the savanna and for-
est ecotypes in the ecotone exhibited similar overall genetic diver-
sity. The inbreeding coefficients (F<sub>IS</sub>) ranged from -0.005 to 0.330
with three localities from the Cerrado, one from the Atlantic Forest,
and two in the ecotone showing significant positive F<sub>IS</sub> (p < .05 after
sequential Bonferroni correction) (Table 1).

AMOVA showed high divergence among localities (F<sub>ST</sub> = 0.187; Table S4). The hierarchical AMOVA showed that most of the variation among localities is explained by the divergence between ecotypes (F<sub>CT</sub> = 0.222 and F<sub>ST</sub> = 0.307). When the hierarchical AMOVA was performed considering localities of the Atlantic Forest and the Cerrado only, the part of total variation attributed to divergence between ecotypes was even higher (F<sub>CT</sub> = 0.337 and F<sub>ST</sub> = 0.396) (Table S4). Pairwise F<sub>ST</sub> values (p < .05 for all values) were higher between localities from different biomes (F<sub>ST</sub> = 0.327-0.402) than within biomes (F<sub>ST</sub> = 0.072-0.099) and showed large variation in the ecotone (F<sub>ST</sub> = 0.053-0.316, Table S5).

The neighbor-joining (NJ) tree showed that the savanna and for-
est ecotypes are highly divergent from each other (Figure 2). Both
groups showed bootstrap values above 90% indicating a high sup-
port for the grouping of these populations. Furthermore, the NJ tree
showed that the forest populations of the ecotone were divergent
from the populations of the Atlantic Forest (Figure 2).

The number of clusters that best explained the genetic struc-
ture in the *P. reticulata* dataset was *K* = 2 based on the sequential *k*-means clustering method. The DAPC analysis based on this grouping indicated that most of the genetic divergence in *P. reticulata* occurs between forest and savanna ecotypes (Figure 3a). Individuals with intermediate posterior membership probabilities between forest and savanna ecotypes were observed in all ecotone localities, suggesting admixture in the ecotone (Figure 3b). The Bayesian clustering method STRUCTURE showed a similar pattern of genetic divergence in *P. reticulata*. The ΔK method of Evanno et al. (2005) clearly indicated *K* = 2 as the highest hierarchical level of genetic structure with a separation of *P. reticulata* in two genetic groups representing the forest and savanna ecotypes (Figure S2a; Figure 4a). Individuals from localities in the Cerrado and Atlantic Forest did not show admixture,
<table><caption>TABLE 1 Genetic diversity parameters of *Plathymenia reticulata*</caption><thead><tr><th>Locality</th><th>Population code</th><th>State</th><th>N</th><th>A</th><th>A<sub>R</sub></th><th>H<sub>O</sub></th><th>H<sub>E</sub></th><th>F<sub>IS</sub></th></tr></thead><tbody><tr><td colspan="9"><strong>Cerrado</strong></td></tr><tr><td>Pirenópolis</td><td>PRC</td><td>GO</td><td>28</td><td>6.8</td><td>5.40</td><td>0.562</td><td>0.659</td><td><strong>0.149</strong></td></tr><tr><td>Paracatu</td><td>PTC</td><td>MG</td><td>23</td><td>6.7</td><td>5.40</td><td>0.486</td><td>0.607</td><td><strong>0.203</strong></td></tr><tr><td>Várzea da Palma</td><td>VZC</td><td>MG</td><td>17</td><td>6.0</td><td>5.17</td><td>0.495</td><td>0.602</td><td><strong>0.183</strong></td></tr><tr><td>Overall Cerrado</td><td>-</td><td>-</td><td>68</td><td>11.0</td><td>6.36</td><td>0.519</td><td>0.670</td><td><strong>0.226</strong></td></tr><tr><td colspan="9"><strong>Atlantic Forest</strong></td></tr><tr><td>Alto Jequitibá</td><td>AJF</td><td>MG</td><td>22</td><td>5.2</td><td>4.14</td><td>0.361</td><td>0.456</td><td><strong>0.214</strong></td></tr><tr><td>Itapebi</td><td>IPF</td><td>BA</td><td>14</td><td>3.9</td><td>3.72</td><td>0.481</td><td>0.479</td><td>-0.005</td></tr><tr><td>Silva Jardim</td><td>SJF</td><td>RJ</td><td>26</td><td>5.2</td><td>3.96</td><td>0.428</td><td>0.468</td><td>0.087</td></tr><tr><td>Overall Atlantic Forest</td><td>-</td><td>-</td><td>62</td><td>7.9</td><td>4.79</td><td>0.416</td><td>0.494</td><td><strong>0.159</strong></td></tr><tr><td colspan="9"><strong>Ecotone</strong></td></tr><tr><td>Contagem</td><td>COE</td><td>MG</td><td>12</td><td>6.5</td><td>6.19</td><td>0.630</td><td>0.742</td><td>0.157</td></tr><tr><td>Felício dos Santos</td><td>FEE</td><td>MG</td><td>38</td><td>10.0</td><td>6.33</td><td>0.504</td><td>0.749</td><td><strong>0.330</strong></td></tr><tr><td>Ribeirão das Neves</td><td>NEE</td><td>MG</td><td>24</td><td>5.4</td><td>4.45</td><td>0.435</td><td>0.567</td><td><strong>0.237</strong></td></tr><tr><td>Parque Estadual do Sumidouro</td><td>SUE</td><td>MG</td><td>29</td><td>8.8</td><td>6.61</td><td>0.667</td><td>0.747</td><td>0.108</td></tr><tr><td>Overall ecotone</td><td>-</td><td>-</td><td>103</td><td>13.2</td><td>7.13</td><td>0.549</td><td>0.768</td><td><strong>0.286</strong></td></tr><tr><td colspan="9"><strong>Savanna ecotype in ecotone</strong></td></tr><tr><td>Contagem</td><td>s-COE</td><td>MG</td><td>12</td><td>6.5</td><td>6.19</td><td>0.630</td><td>0.742</td><td>0.157</td></tr><tr><td>Felício dos Santos</td><td>s-FEE</td><td>MG</td><td>16</td><td>7.4</td><td>6.07</td><td>0.436</td><td>0.626</td><td><strong>0.312</strong></td></tr><tr><td>Parque Estadual do Sumidouro</td><td>s-SUE</td><td>MG</td><td>15</td><td>7.6</td><td>6.54</td><td>0.614</td><td>0.660</td><td>0.072</td></tr><tr><td>Overall ecotone savanna ecotype</td><td>-</td><td>-</td><td>43</td><td>11.8</td><td>7.09</td><td>0.551</td><td>0.701</td><td><strong>0.216</strong></td></tr><tr><td colspan="9"><strong>Forest ecotype in ecotone</strong></td></tr><tr><td>Felício dos Santos</td><td>f-FEE</td><td>MG</td><td>22</td><td>6.8</td><td>5.17</td><td>0.554</td><td>0.669</td><td><strong>0.176</strong></td></tr><tr><td>Ribeirão das Neves</td><td>f-NEE</td><td>MG</td><td>24</td><td>5.4</td><td>4.45</td><td>0.435</td><td>0.567</td><td><strong>0.237</strong></td></tr><tr><td>Parque Estadual do Sumidouro</td><td>f-SUE</td><td>MG</td><td>14</td><td>6.5</td><td>5.91</td><td>0.718</td><td>0.769</td><td>0.069</td></tr><tr><td>Overall ecotone forest ecotype</td><td>-</td><td>-</td><td>60</td><td>10.0</td><td>6.32</td><td>0.549</td><td>0.706</td><td><strong>0.224</strong></td></tr></tbody></table>

Note: A, mean number of alleles per locus; A<sub>R</sub>, allelic richness; H<sub>O</sub>, mean observed heterozygosity; H<sub>E</sub>, mean expected heterozygosity; F<sub>IS</sub>, inbreeding coefficient. Boldface numbers represent significant F<sub>IS</sub> values (p < .05).

whereas admixture was frequent in individuals from the ecotone (Figure 4a), considering a Q value higher than 0.850 (or lower than 0.150) as a threshold for a pure individual. The highest proportion of admixed individuals was found in f-SUE (86% of individuals) and f-FEE (41% of individuals) in the forest ecotype and in s-COE (67% of individuals) in the savanna ecotype.

In agreement with the STRUCTURE and DAPC, the recent migration rates between ecotypes estimated in BayesAss were high in the ecotonal region (Table 2). The migration rate estimates were higher into the forest ecotype in two of three ecotonal population pairs sampled, with values ranging from 0.024 to 0.119 into the forest ecotype and values ranging from 0.030 to 0.061 into the savanna ecotype (Table 2). When considering the whole sample, which included ecotone, Cerrado, and Atlantic Forest localities, the

migration rate into the forest ecotype was also higher than that into the savanna ecotype with values of 0.040 (SD = 0.01) and 0.024 (SD = 0.09), respectively, and the same was true for the analysis restricted to ecotone localities (Table 2).

In the evaluation of the structure using only the pure individuals of each ecotype in DAPC, we found K = 5 in the forest and K = 7 in the savanna ecotype. In the forest ecotype, all the localities showed a prevailing genetic cluster despite showing high proportions of shared ancestry mainly among the three populations of the Atlantic Forest and among the three ecotonal localities (Figure 3c). In the savanna ecotype, individual assignment probabilities to the clusters were high, but clusters were shared among localities (Figure 3d). The ΔK method of Evanno et al. (2005) on the STRUCTURE analysis of the pure individuals of each ecotype showed K = 3 within both forest
<figure><img src="image_4.png" /><figcaption>FIGURE 2 Neighbor-joining tree based on D<sub>ST</sub> genetic distance among *Plathymenia reticulata* populations with bootstrap support of nodes</figcaption></figure>

and savanna ecotypes (Figure 4, Figure S2b,c). Using the criteria of StructureSelector, K = 5 was the best number of clusters within both ecotypes (Figure 4, Figure S3a,b). Both ecotypes displayed a genetic structure that was partially congruent with the geographic localities, though most localities featured more than one cluster as well as admixed individuals, especially those of the forest ecotype (Figure 4).

## 3.2 | Genotypic composition of *Plathymenia reticulata* in the hybrid zone

The classification of *P. reticulata* individuals in the NewHybrids analysis supported, with high posterior probability (PP > 0.90), a pure parental composition for most individuals in the areas of the Cerrado and Atlantic Forest (Figure 5a). In the ecotone, 18 individuals were classified as F2, with 12 of them resembling the forest ecotype and 6 the savanna ecotype (Figure 5a). Twenty additional individuals were identified as hybrids with posterior probability values between 0.4 and 0.8 for different hybrid genotypic classes involving F1, F2, and backcross to the forest ecotype (B2), but the posterior probability of backcross to the savanna ecotype (B1) never exceeding 10% (Figure 5a). The localities with the highest proportions of F2 individuals were s-COE, f-FEE, and f-SUE with 41%, 27%, and 21% of F2 individuals, respectively (Figure 5a). The classification in 12 genotypic classes assigned as pure, a proportion of 0.636–0.857 of individuals from Atlantic Forest and Cerrado localities, respectively (Figure 5b). Three individuals showed high posterior probabilities of being classified as second-generation backcrosses (B7 or B8) in AJF, PRS, and PTS localities (Figure 5b). Furthermore, in f-SUE, 14% of individuals were classified as F1 when using 12 genotypic classes. The high proportion of individuals classified as later generation hybrids in this

analysis may indicate that the analysis using six genotypic classes misclassified later generation hybrids as pure.

Consistent with the NewHybrids results, the mean hybrid index (Figure 5c) of the ecotone localities indicated high admixture levels with mean values per locality ranging from 0.083 to 0.771 and lower bounds of the 95% confidence (CI) interval ranging from 0.024 to 0.540 and upper bounds of the 95% CI ranging from 0.279 to 0.914. The proportion of individuals per population with a hybrid index between 0.200 and 0.800 ranged from 0.133 to 0.929 in the ecotone localities, with the highest values in COE, NEE, and f-SUE populations (Figure 5c). In comparison, the localities of both biomes showed mean hybrid index values ranging from 0.990 to 0.999 in the Atlantic Forest and from 0.008 to 0.010 in the Cerrado (Figure 5c). The individual interecotypic heterozygosity ($Q_{12}$) ranged from 0.00 to 0.818, with the ecotonal localities showing higher mean values ($Q_{12}$ ranging from 0.261 to 0.403) than the localities of each biome ($Q_{12}$ ranging from 0.071 to 0.104) (Figure 6a). In addition, the interecotypic heterozygosity × hybrid index plots showed a very dispersed distribution for individuals, which indicates that the ecotonal localities mostly comprise later generation hybrids and backcrossed individuals (Figure 6b). A population comprising mostly F1s would show a mean hybrid index closer to 0.5 with interecotypic heterozygosity values closer to the maximum value of 1.0. Finally, some individuals showed higher values than expected for interecotypic heterozygosity with values above the triangle bounds. This occurred mainly in the individuals of localities from Cerrado and Atlantic Forest suggesting that homoplasy may be found in the microsatellite dataset.

## 3.3 | Evaluation of hybrid assignment in NewHybrids

The model evaluation using NewHybrids with six simulated genotypic classes showed an average proportion of correct assignments of pure parental individuals of almost 91% when setting a threshold of 0.85 (Table S6; Figure S4). Only the F2 category showed a proportion of correct assignments lower than 50%, while F1, B1, and B2 displayed 79%, 66%, and 65% correct assignments, respectively. When using 14 simulated crosses for the assignment in 12 genotypic classes, a very low proportion of correct assignment of pure individuals was found, with most of them showing values of posterior probability of assignment below 0.85 (Table S6; Figure S4). The classification was better for the hybrid classes, especially for F1 and F2 with 57% and 59% of individuals correctly assigned to their respective categories (Table S6). The classification of backcrosses never exceeded 10% of correct assignments suggesting that our data could not resolve later generation hybrids correctly, especially for second-generation backcrosses (Table S6). These results indicate a relatively high capacity in the differentiation between admixed and nonadmixed individuals, but a poorer capacity to differentiate among advanced generation hybrids (F2 and backcrosses). Finally, these simulations let us assume that some individuals classified as F2 in the NewHybrids analysis with
FIGURE 3 (a) Distribution of the score values of individuals of each ecotype of *Plathymenia reticulata* based k-means grouping and classification with Discriminant Analysis of Principal Components (DAPC). (b) Bar plot showing individual membership probabilities based on DAPC. (c) Bar plot showing individual membership probabilities based on DAPC for pure individuals from the forest ecotype. (d) Bar plot showing individual membership probabilities based on DAPC for pure individuals from the savanna ecotype

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

six genotypic classes could in fact be backcrosses (Figure S5). The evaluation of the performance of STRUCTURE and *introgress* to identify and classify hybrids based on simulated data is presented in Appendix S1.

## 3.4 | Drivers of divergence of *Plathymenia reticulata* ecotypes

Simple Mantel tests indicated that the dataset displayed both a signature of isolation by distance and of isolation by environment (Table 3). Genetic distance showed a higher correlation with bioclimatic distance (r = .381, *p* = .006) than with spatial distance (r = .307, *p* = .009) or soil distance (r = .308, *p* = .023). The partial Mantel test showed that the correlation of the bioclimatic distance with the genetic distance remained significant when controlling for geographical distance (r = .260, *p* = .035, Table 3). Other partial Mantel tests were not significant.

The sPCA test indicated a significant global structure (*p* = .042) with the first three axes showing significant divergence between the *P. reticulata* populations (Figure S5). The dbRDA analysis identified an association of genetic divergence based on sPCA with Cation Exchange Capacity (CEC), the fourth synthetic variable of the climatic PCA (CLIM4), and with spatial distance between populations (PCNM2) (Table 4). CLIM4 was mainly associated with Precipitation of the Warmest Quarter (bio18, loading = -0.645), Mean Diurnal Range Temperature (bio2, loading = 0.398), and Precipitation of the Wettest Quarter (bio16, loading = 0.309) (Table S3). The CEC and PCNM2 were associated with the first dbRDA axis with loadings of -0.958 and -0.815, respectively. The second dbRDA axis was associated with the CLIM4 variable with a loading value of 0.813 (Table S8). The estimated partial adjusted R² showed that CEC was the variable with the highest explanatory power, with CLIM4 being the second (Table 4). The first axis of dbRDA separated the populations of the Cerrado from those of the Atlantic Forest with the ecotonal populations being intermediary (Figure 7).
<figure><img src="image_6.png" /></figure>

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

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

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

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

Individuals
FIGURE 4 Population structure based on the Bayesian clustering method of STRUCTURE. (a) Bar plot showing individual admixture proportions in K = 2 genetic groups for all individuals of *Plathymenia reticulata*. Bar plots showing individual admixture proportions for (b) K = 3 and (c) K = 5 genetic groups for pure individuals of the forest ecotype. Bar plots showing individuals admixture proportions for (d) K = 3 and (e) K = 5 genetic groups for pure individuals of the savanna ecotype

## 4 | DISCUSSION

Our study reveals genetic divergence between savanna and forest ecotypes of *P. reticulata* and gene flow between ecotypes across ecotonal areas between the Brazilian Cerrado and Atlantic Forest, characterizing a hybrid zone. Hybrid individuals mostly corresponded to second-generation hybrids or backcrosses toward the forest ecotype. We observed asymmetric gene flow between ecotypes, which was higher into the forest ecotype. Our data also suggest that besides the geographical distance, climatic and edaphic variables partially explain the genetic divergence between the ecotypes. The ecotonal areas harbor high genetic diversity relative to the Atlantic Forest and Cerrado biomes for *P. reticulata* and can be an important source of variation for the species, mainly for the forest ecotype.

### 4.1 | Genetic divergence between savanna and forest ecotypes of *Plathymenia reticulata*

The analyses of the genetic structure and divergence of *P. reticulata* showed a clear genetic distinction between savanna and forest ecotypes, consistent with our a priori classification of individuals based on morphological characters of the tree trunk, bark, height, and environment of the trees. The genetic distinction between ecotypes was reinforced by the fact that even in the ecotonal localities this division was recovered using NJ trees, DAPC, and STRUCTURE analyses. Populations from the forest ecotype in the ecotone are more

<table><caption>TABLE 2 Recent immigration rates (m) between forest and savanna ecotypes of *Plathymenia reticulata* estimated by BayeAss</caption><thead><tr><th>Parameter</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td>All localities</td><td></td><td></td></tr><tr><td>m to forest</td><td>0.040</td><td>0.010</td></tr><tr><td>m to savanna</td><td>0.024</td><td>0.009</td></tr><tr><td>All ecotonal localities</td><td></td><td></td></tr><tr><td>m to forest</td><td>0.072</td><td>0.019</td></tr><tr><td>m to savanna</td><td>0.057</td><td>0.021</td></tr><tr><td>COE/NEE localities</td><td></td><td></td></tr><tr><td>m to forest</td><td>0.024</td><td>0.019</td></tr><tr><td>m to savanna</td><td>0.054</td><td>0.035</td></tr><tr><td>FEE locality</td><td></td><td></td></tr><tr><td>m to forest</td><td>0.039</td><td>0.022</td></tr><tr><td>m to savanna</td><td>0.030</td><td>0.023</td></tr><tr><td>SUE locality</td><td></td><td></td></tr><tr><td>m to forest</td><td>0.119</td><td>0.076</td></tr><tr><td>m to savanna</td><td>0.061</td><td>0.031</td></tr></tbody></table>

genetically related to populations from the Atlantic Forest, whereas populations of the savanna ecotype in the ecotone are more related to populations from the Cerrado, despite high genetic admixture in the ecotone. These genetic data are consistent with ecological studies that showed that the two ecotypes maintain their morphophysiological traits in the ecotonal areas, though ecotonal individuals are less dissimilar to the alternative ecotype than when comparing distinct ecotypes in their specific habitats (Goulart et al., 2005, 2006, 2011; Toledo et al., 2012). Differences between ecotypes in several of these traits were demonstrated to be genetically determined when progenies were cultivated in a common garden experiment with different light conditions (Goulart et al., 2011).

We found an association of edaphic and climatic variables with *P. reticulata* genetic divergence that was in some instances stronger than the divergence explained by spatial distance between localities, suggesting a pattern of isolation by environment (IBE). The genetic divergence was associated with Cation Exchange Capacity (CEC) of the soil and with Precipitation of the Warmest Quarter, Precipitation of the Wettest Quarter, and Mean Diurnal Range Temperature (summarized in the synthetic variable CLIM4). We detected IBE using SSR markers that are considered neutral, congruent with a scenario of genome-wide hitchhiking following long-lasting selection that may have affected a small fraction of the genome, in combination with assortative mating (Andrew et al., 2012; Sexton et al., 2014, 2016). The high proportion of variance explained by CCE suggests that the differences in soil between savanna and forest can be important in the divergence in *P. reticulata*. Differences in soil can affect plant growth strategy and be an important driver of plant functional traits and adaptive divergence even at small spatial scales (Schmitt et al., 2021), as is the case in our study where patches of savanna and forest form a mosaic of abrupt habitat shifts within meters in ecotonal areas (Lambers et al., 2011; Reich et al., 2003). Previous studies have shown the adaptive significance of intraspecific variation for population persistence in forest-savanna ecotones, where different resource-use strategies are required in each habitat (Maracahipes et al., 2018; Silva et al., 2019). Soils of Atlantic Forest localities are more fertile and have higher CEC than the Cerrado localities, with ecotonal localities showing intermediary values. CEC can be used as a proxy of soil fertility and nutrient retention, as it relates to the soil's capacity to provide calcium, magnesium, and potassium for plants (https://www.fao.org/; https://www.soilquality.org.au/). Therefore, the high genetic variation we observed across very close-by sites in the ecotone could suggest high adaptive potential, assortative mating of adapted individuals to local soil, and selection against maladapted individuals into the alternative habitat, leading to a pattern of IBE.

Regarding bioclimatic variables, the Mean Diurnal Range in temperature, the Precipitation of the Wettest Quarter, and the
<figure><img src="image_11.png" /></figure>

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

<figure><img src="image_13.png" /><figcaption>FIGURE 5 Genotypic composition of *Plathymenia reticulata*. (a) Genotypic assignment of each individual in two parental populations and four hybrid classes based on NewHybrids. (b) Genotypic assignment of each individual in the six mentioned genotypic classes and six advanced generation hybrids based on NewHybrids. (c) Hybrid index estimated for each individual. P1 = savanna ecotype, P2 = forest ecotype, F1 = P1 × P2, F2 = F1 × F1, B1 = P1 × F11, B2 = P2 × F1, B3 = P1 × F2, B4 = P2 × F2, B5 = P1 × (F1 × P1), B6 = P2 × (F1 × P2), B7 = P1 × (P1 × F2), B8 = P2 × (P2 × F2)</figcaption></figure>

Precipitation of Warmest Quarter were associated with genetic diver-
gence and are related with extremes of temperature within a day and
the seasonality in rainfall along the seasons, respectively. Associations
between genetic divergence and the Mean Diurnal Range (Buzatti
et al., 2019) or between genetic diversity and Precipitation of the
Wettest Quarter (Ribeiro et al., 2016) were already found for *Qualea
grandiflora* and *Annona coriacea*, respectively, two of the most common
tree species of the Cerrado. A similar association was also found in an
annual herb, *Brasilianthus carajensis*, in a mountain savanna in *Canga*

plateaus within the Amazon Forest (Silva et al., 2020). The Cerrado is
a more stressful environment than the Atlantic Forest in many ways,
as it is, for example, characterized by a stronger and more prolonged
drought season. Accordingly, in *P. reticulata*, previous studies on some
of the sampling locations of the present study (COE and NEE in Goulart
et al., 2006, and SUE in Toledo et al., 2012) showed that savanna and
forest ecotypes are differentiated for drought response traits even
in ecotonal areas. Differences between savanna and forest ecotypes
concern also the vegetative phenology, with leaf shed beginning earlier
<figure><img src="image_14.png" /><figcaption>FIGURE 6 (a) Boxplot showing the distribution of interecotype heterozygosity values by sampling localities of *Plathymenia reticulata*. (b) Hybrid index × interecotype heterozygosity values showing the composition of hybrids in the hybrid zone</figcaption></figure>

<table><caption>TABLE 3 Mantel tests showing the correlation between genetic distance and environmental distances (soil and climate) between populations of *Plathymenia reticulata*</caption><thead><tr><th>Soil</th><th>Gen x Soil</th><th>Gen x Geo</th><th>Gen x Soil (geo)</th><th>Gen x Geo (soil)</th></tr></thead><tbody><tr><td>Correlation</td><td>0.308</td><td>0.307</td><td>0.163</td><td>0.163</td></tr><tr><td>p</td><td>.023</td><td>.009</td><td>.137</td><td>.104</td></tr><tr><td>Climate</td><td>Gen x Clim</td><td>Gen x Geo</td><td>Gen x Clim (geo)</td><td>Gen x Geo (clim)</td></tr><tr><td>Correlation</td><td>0.381</td><td>0.307</td><td>0.260</td><td>0.111</td></tr><tr><td>p</td><td>.006</td><td>.008</td><td>.035</td><td>.196</td></tr></tbody></table>

Note: Control variables are given within brackets for partial Mantel tests.
<table><thead><tr><th>Variable</th><th>df</th><th>Sum of Sqs</th><th>F</th><th>Pr(&gt;F)</th><th>Ajusted.R<sup>2</sup></th><th>VIF</th></tr></thead><tbody><tr><td>CEC</td><td>1</td><td>10.225</td><td>7.814</td><td>0.0001</td><td>0.301</td><td>2.228</td></tr><tr><td>CLIM4</td><td>1</td><td>6.793</td><td>5.191</td><td>0.0007</td><td>0.227</td><td>1.368</td></tr><tr><td>PCNM2</td><td>1</td><td>5.514</td><td>4.214</td><td>0.0029</td><td>0.102</td><td>2.326</td></tr><tr><td>Residual</td><td>8</td><td>10.469</td><td></td><td></td><td></td><td></td></tr></tbody></table>

Abbreviations: CEC, Cation Exchange Capacity; CLIM4, Climatic Axis obtained from PCA; PCNM2, Spatial axis obtained from Principal Coordinates Neighboring Matrix; VIF, Variance Inflation Factor.

<caption>TABLE 4 Results of model selection and redundancy analysis testing the association between allelic variation of *Plathymenia reticulata*, spatial distance, and environmental variables</caption>

<figure><img src="image_15.png" /><figcaption>FIGURE 7 Scatter plot showing the distribution of redundancy analysis (RDA) scores of *Plathymenia reticulata* populations and the correlation of environmental variables with the RDA axis</figcaption></figure>

in the Cerrado than in ecotonal areas, and later in the forest (Goulart et al., 2005). In addition, the forest ecotype shows higher stem radial increment, positively related to increasing precipitation, relative to the savanna ecotype, which shows lower midday water potential during the dry season relative to the forest ecotype (Toledo et al., 2012). Trait divergence thus appears in accordance with environmental variation and genetic divergence in *P. reticulata*. Conclusive evidence on the genetic signature of natural selection in the species could be obtained from genome-wide data that might reveal contrasted differentiation of neutral versus potentially adaptive alleles in ecotype divergence.

The high genetic divergence between ecotypes accounting for 22.2–33.7% of the total genetic variance (Table S2) raises questions about the taxonomy of the genus. In comparison, divergence between two *Hymenaea* (Fabaceae) tree species from the Atlantic Forest and Cerrado accounted for 26% of the total genetic variance based on seven nuclear microsatellites (mean $H_E$ = 0.600, Resende-Moreira et al., 2017). In *Eucalyptus globulus*, conversely, differentiation between ecotypes was lower, with pairwise $F_{ST}$ values ranging from 0.03 to 0.21 based on using 12 nuclear microsatellites (mean $H_E$ = 0.84, Foster et al., 2007). Pairwise $F_{ST}$ values between *P. reticulata* localities of the Atlantic Forest and the Cerrado ranged from 0.327 to 0.402. This relatively high level of genetic divergence between *P. reticulata* ecotypes associated to morphophysiological differences found in several previous studies suggests that they could

represent incipient species. These observations are in accordance with the older taxonomic status that considered the *Plathymenia* genus as represented by two species, *P. reticulata* Benth in the Cerrado and *P. foliolosa* Benth in the Atlantic Forest (Heringer, 1956).

## 4.2 | Characterization of *Plathymenia reticulata* in the hybrid zone

The high levels of genetic admixture in several ecotonal areas, with a mean of 40% of hybrids per locality, characterize these areas as a hybrid zone. Most hybrid individuals were classified as later generation hybrids with none showing high probabilities to be F1s and very few showing posterior probability between 0.5 and 0.9 to be classified as F1. Classifying hybrids in genotypic classes based on a few microsatellite markers is challenging, as shown by our simulation study which illustrated the difficulty of distinguishing later generation hybrids from pure individuals depending on the chosen assignment threshold. It should be noted that some studies have shown that individuals previously classified as F2 based on microsatellite data were in fact classified as F1s when using genomic data (Gramlich et al., 2018; Lindtke et al., 2014). In the microsatellite dataset, these misclassified individuals showed levels of interspecies ancestry lower than those expected for F1 hybrids, but unusually higher than those
predicted for later generation hybrids. In our dataset, some individu-
als showed this unusual pattern of high interecotype heterozygosity,
which may suggest that our analysis underestimated the number of
F1s. Therefore, a study based on genomic data may better estimate
the number of early generation hybrids along the *P. reticulata* hybrid
zone. Such a study will also increase the power of classifying later
generation hybrids and help shed light on which hybridization events
and backcross scenarios provide opportunities for introgression be-
tween the ecotypes.

We observed higher gene flow in direction of the forest eco-
type than toward the savanna ecotype. Several causes such as
demographic processes related to range expansion may cause
asymmetric hybridization (Currat et al., 2008; Field et al., 2010).
The savanna ecotype populations could be expanding into forested
habitats and the interbreeding between the savanna and forest
ecotypes in these areas would lead to the invasion of genes of
forest populations by genes of the savanna ecotype. Differences
in dispersal capacity may also cause differential introgression
with the populations with higher dispersal capacity promoting the
spread of alleles to the other populations (Yamamoto et al., 2019).
The *P. reticulata* forest ecotype showed seed structures associ-
ated with higher dispersal capacity (Goulart et al., 2006), which
can indicate that this factor is not the better explanation for the
differential introgression. However, the occurrence of the savanna
ecotype in open habitats may provide opportunities for dispersal
across longer distances which may lead to the capacity of savanna
ecotypes to colonize surrounding forested areas. Natural selec-
tion can also act as a driver of differential introgression (Lexer
et al., 2005; Whitney et al., 2006). Ecological barriers, that is, a
harsh environment, may effectively reduce the gene flow into sa-
vanna populations. Another hypothesis to explain the differences
in gene flow could be related to the longer juvenile stage of forest
individuals compared to the savanna individuals in which frequent
fires appear to select for earlier reproductive maturity (Hoffmann
& Franco, 2003). Furthermore, the Cerrado biome besides typi-
cal savanna vegetation presents more dense wood formations
(Cerradão) and riparian forests where *P. reticulata* savanna indi-
viduals can be found with more similar characteristics to forest
individuals. Thus, the higher effective gene flow toward the forest
ecotype may be due to the fact that some Cerrado individuals are
more able to colonize and survive in forest habitat.

Novaes et al. (2010) suggested that *P. reticulata* lineages di-
 diversified from its center of distribution in the Cerrado shaped by
climatic changes in the Pleistocene. The lower genetic diversity of
the forest ecotype in the Atlantic Forest, compared to the savanna
ecotype, suggests recent expansion of *P. reticulata* toward the
Atlantic Forest possibly involving founder effects that led to the
loss of genetic diversity. The spatial dynamics of Atlantic Forest
and Cerrado biomes during glacial cycles in the Pleistocene, with
the Cerrado progressing eastward into the forests during drier and
colder periods (Behling, 1998; Behling & Negrelle, 2001; Meneses
et al., 2013), may have provided opportunities for the expansion
of *P. reticulata* from the savanna toward previously forested areas.

Phylogeographical studies suggest similar dynamics of range
expansions and contractions in tree species of the Cerrado and
Atlantic Forest during the Pleistocene (Buzatti et al., 2017; Buzatti
et al., 2019; Leal et al., 2016; Novaes et al., 2013; Ribeiro et al.,
2011).

Gene flow in vicariant species such as *P. reticulata* with subspe-
cific taxa typical from forests and savannas was also identified in a
few other tree species in savanna-forest boundaries including eco-
tonal and riparian forests in the Cerrado biome (Cavallari et al., 2010;
Muniz et al., 2020; Resende-Moreira et al., 2017). This can suggest
that hybridization in the ecotonal region between Atlantic Forest
and Cerrado can be a relevant evolutionary process increasing the
genetic diversity of the species of these biomes. The study of other
closely related species of other groups can show how common and
phylogenetically widespread is the gene flow between lineages of
Cerrado and Atlantic Forest.

## 4.3 | Implications for conservation

The highest genetic diversity found here in *Plathymenia* popula-
tions of the savanna-forest ecotone may be due to the hybridiza-
tion between ecotypes/species across this area. This is a pattern
commonly found in hybrid zones where the admixture between
lineages increases genetic diversity (Carlson et al., 2017; Eaton
et al., 2015; Y. Li et al., 2016; Zalapa et al., 2010). Therefore, these
ecotonal areas are an important source of genetic diversity for
*Plathymenia*, mainly for Atlantic Forest populations that show low
genetic diversity. The Atlantic Forest has a history of extensive
habitat loss and fragmentation and currently remains mostly in
small and isolated patches (Ribeiro et al., 2009). In comparison,
the deforestation in the Cerrado, though very quick, was only
significant from the middle of the twentieth century (Silva et al.,
2006). The low genetic diversity of the Atlantic Forest populations
indicates that they may have lower adaptive potential and could
be severely impacted by ongoing habitat loss and fragmentation
in the Atlantic Forest. The forest ecotype, previously considered
as *P. foliolosa*, is considered vulnerable in the IUCN red list (World
Conservation Monitoring Centre, 1998). The conservation ef-
forts for the forest ecotype should be increased and this taxon
should be considered a different management unit from the sa-
vanna ecotype. Thus, in addition to the conservation efforts in the
Atlantic Forest and Cerrado areas, our data point to the impor-
tance of including ecotonal areas between these biomes, which
could carry alleles to augment the genetic diversity in the Atlantic
Forest.

5 | CONCLUSION

Based on multiple independent approaches, with distinct assump-
tions, we showed that the ecotonal areas between the Cerrado and
the Atlantic Forest constitute a hybrid zone between the genetically
differentiated forest and savanna ecotypes of *P. reticulata*. Our results suggest the suitability of the ecotonal areas for studies of evolutionary divergence and speciation. The two *P. reticulata* ecotypes should be considered as different species or at least two different management units with special attention to the forest ecotype. The dynamics of ecotonal areas, as evidenced here, can have a significant role in the origin of the high biodiversity in the Cerrado and Atlantic Forest, highlighting the importance of these areas as hotspot of genetic diversity and conservation.

## ACKNOWLEDGMENTS

We thank the Instituto Chico Mendes de Conservação da Biodiversidade (ICMBio) and the Sistema de Autorização e Informação em Biodiversidade (SISBIO) for providing the license to collect genetic samples. We thank Rosangela L. Brandão for technical assistance in the laboratory. This work was supported by the Fundação de Amparo a Pesquisa do Estado de Minas Gerais (FAPEMIG). We thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for their support to the first author's PhD scholarship and research fellowships to JP Lemos-Filho and MB Lovato, respectively. MH acknowledges support from an Investissement d'Avenir grant of the ANR: CEBA (ANR-10-LABEX-0025).

## CONFLICT OF INTEREST

The authors declare no conflict of interest.

## AUTHOR CONTRIBUTIONS

**André Carneiro Muniz:** Formal analysis (lead); investigation (equal); methodology (equal); writing – original draft (lead); writing – review and editing (equal). **Ricardo José Gonzaga Pimenta:** Formal analysis (supporting); methodology (equal). **Mariana Vargas Cruz:** Formal analysis (supporting); investigation (equal); methodology (equal). **Jacqueline Gomes Rodrigues:** Investigation (supporting); methodology (equal). **Renata Santiago de Oliveira Buzatti:** Investigation (supporting); methodology (equal). **Myriam Heuertz:** Formal analysis (lead); investigation (equal); methodology (equal); writing – review and editing (equal). **Jose P. Lemos-Filho:** Conceptualization (equal); funding acquisition (equal); investigation (equal); resources (equal); writing – review and editing (equal). **Maria Bernadete Lovato:** Conceptualization (equal); formal analysis (equal); funding acquisition (lead); investigation (lead); methodology (equal); project administration (lead); supervision (lead); writing – review and editing (equal).

## DATA AVAILABILITY STATEMENT

Genotype data, geographical coordinates, and scripts to simulate individuals in distinct genotypic classes will be archived in Dryad under https://doi.org/10.5061/dryad.905qfttn0.

## ORCID

Mariana Vargas Cruz https://orcid.org/0000-0001-8569-2234

Maria Bernadete Lovato https://orcid.org/0000-0001-7009-3079

## REFERENCES

Abbott, R. J. (2017). Plant speciation across environmental gradients and the occurrence and nature of hybrid zones. *Journal of Systematics and Evolution*, 55(4), 238–258. https://doi.org/10.1111/jse.12267

Abbott, R., Albach, D., Ansell, S., Arntzen, J. W., Baird, S. J. E., Bierne, N., Boughman, J., Brelsford, A., Buerkle, C. A., Buggs, R., Butlin, R. K., Dieckmann, U., Eroukhmanoff, F., Grill, A., Cahan, S. H., Hermansen, J. S., Hewitt, G., Hudson, A. G., Jiggins, C., ... Zinner, D. (2013). Hybridization and speciation. *Journal of Evolutionary Biology*, 26(2), 229–246. https://doi.org/10.1111/j.1420-9101.2012.02599.x

Anderson, E. C., & Thompson, E. A. (2002). A model-based method for identifying species hybrids using multilocus genetic data. *Genetics*, *160*(3), 1217–1229. https://doi.org/10.1093/genetics/160.3.1217

Andrew, R. L., Ostevik, K. L., Ebert, D. P., & Rieseberg, L. H. (2012). Adaptation with gene flow across the landscape in a dune sunflower. *Molecular Ecology*, 21(9), 2078–2091. https://doi.org/10.1111/j.1365-294X.2012.05454.x

Azani, N., Babineau, M., Bailey, C. D., Banks, H., Barbosa, A. R., Pinto, R. B., Boatwright, J. S., Borges, L. M., Brown, G. K., Bruneau, A., Candido, E., Cardoso, D., Chung, K.-F., Clark, R. P., Conceição, A. D. S., Crisp, M., Cubas, P., Delgado-Salinas, A., Dexter, K. G., ... Zimmerman, E. (2017). A new subfamily classification of the Leguminosae based on a taxonomically comprehensive phylogeny – The Legume Phylogeny Working Group (LPWG). *Taxon*, 66(1), 44–77. https://doi.org/10.12705/661.3

Behling, H. (1998). Late Quaternary vegetational and climatic changes in Brazil. *Palaeobotany and Palynology*, 99(2), 143–156. https://doi.org/10.1016/S0034-6667(97)00044-4

Behling, H., & Negrelle, R. R. B. (2001). Tropical rain forest and climate dynamics of the Atlantic Lowland, Southern Brazil, during the Late Quaternary. *Quaternary Research*, 56(3), 383–389. https://doi.org/10.1006/qres.2001.2264

Bjornstad, O. N. (2020). Package 'ncf': Spatial nonparametric covariance functions. R package version 1.2-9. https://CRAN.R-project.org/package=ncf

Buck, R., Hyasat, S., Hossfeld, A., & Flores-Rentería, L. (2020). Patterns of hybridization and cryptic introgression among one- and four-needled pinyon pines. *Annals of Botany*, 126(3), 401–411. https://doi.org/10.1093/aob/mca045

Buerkle, C. A. (2005). Maximum-likelihood estimation of a hybrid index based on molecular markers. *Molecular Ecology Notes*, 5(3), 684–687. https://doi.org/10.1111/j.1471-8286.2005.01011.x

Buzatti, R. S. O., Lemos-Filho, J. P., Bueno, M. L., & Lovato, M. B. (2017). Multiple *Pleistocene refugia* in the Brazilian cerrado: Evidence from phylogeography and climatic nichemodelling of two Qualea species (Vochysiaceae). *Botanical Journal of the Linnean Society*, 185(3), 307–320. https://doi.org/10.1093/botlinnean/box062

Buzatti, R. S. D. O., Pfeilsticker, T. R., Muniz, A. C., Ellis, V. A., de Souza, R. P., Lemos-Filho, J. P., & Lovato, M. B. (2019). Disentangling the environmental factors that shape genetic and phenotypic leaf trait variation in the tree *Qualea grandiflora* across the Brazilian Savanna. *Frontiers in Plant Science*, 10. https://doi.org/10.3389/fpls.2019.01580

Carlson, M. O., Gazave, E., Gore, M. A., & Smart, C. D. (2017). Temporal genetic dynamics of an experimental, biparental field population of *Phytophthora capsici*. *Frontiers in Genetics*, 8, 26. https://doi.org/10.3389/fgene.2017.00026

Cavallari, M. M., Gimenes, M. A., Billot, C., Torres, R. B., Zucchi, M. I., Cavalheiro, A. J., & Bouvet, J. M. (2010). Population genetic relationships between *Casearia sylvestris* (Salicaceae) varieties occurring sympatrically and allopatrically in different ecosystems in south-east Brazil. *Annals of Botany*, 106(4), 627–636. https://doi.org/10.1093/aob/mcq151

Chhatre, V. E., Evans, L. M., DiFazio, S. P., & Keller, S. R. (2018). Adaptive introgression and maintenance of a trispecies hybrid complex in
range-edge populations of Populus. *Molecular Ecology*, **27**(23),
4820-4838. https://doi.org/10.1111/mec.14820

Choi, J. Y., Purugganan, M. D., & Stacy, E. A. (2019). Divergent selection and primary gene flow shape incipient speciation of a riparian tree on Hawaii Island. *Molecular Biology and Evolution*, **37**(3), 695-710.
https://doi.org/10.1093/molbev/msz259

Christie, K., & Strauss, S. Y. (2018). Along the speciation continuum: Quantifying intrinsic and extrinsic isolating barriers across five million years of evolutionary divergence in California jewelflowers. *Evolution*, **72**(5), 1063-1079. https://doi.org/10.1111/evo.13477

Coyne, J. A. (1992). Genetics and speciation. *Nature*, **355**(6), 511-515.
https://doi.org/10.1038/355511a0

Cruz, M. V., Rodrigues, J. G., Souza, H. A. V. E., & Lovato, M. B. (2012). Isolation and characterization of microsatellite markers for *Plathymenia reticulata* (Fabaceae). *American Journal of Botany*, **99**(5), 210-212.

Currat, M., Ruedi, M., Petit, R. J., & Excoffier, L. (2008). The hidden side of invasions: Massive introgression by local genes. *Evolution*, **62**(8), 1908-1920. https://doi.org/10.1111/j.1558-5646.2008.00413.x

De Lemos-Filho, J. P., Goulart, M. F., & Lovato, M. B. (2008). Populational approach in ecophysiological studies: The case of *Plathymenia reticulata*, a tree from Cerrado and Atlantic Forest. *Brazilian Journal of Plant Physiology*, **20**(3), 205-216. https://doi.org/10.1590/S1677-04202008000300005

Dufresnes, C., Strachinis, I., Suriadna, N., Mykytynets, G., Cogălniceanu, D., Székely, P., & Denoël, M. (2019). Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates). *Molecular Ecology*, **28**(13), 3257-3270. https://doi.org/10.1111/mec.15133

Durigan, G., & Ratter, J. A. (2006). Successional changes in Cerrado and Cerrado/forest ecotonal vegetation in western São Paulo State, Brazil, 1962-2000. *Edinburgh Journal of Botany*, **63**(1), 119-130.
https://doi.org/10.1017/S0960428606000357

Earl, D. A., & vonHoldt, B. M. (2012). STRUCTURE HARVESTER: A web-site and program for visualizing STRUCTURE output and implementing the Evanno method. *Conservation Genetics Resources*, **4**(2), 359-361. https://doi.org/10.1007/s12686-011-9548-7

Eaton, D. A. R., Hipp, A. L., González-Rodríguez, A., & Cavender-Bares, J. (2015). Historical introgression among the American live oaks and the comparative nature of tests for introgression. *Evolution*, **69**(10), 2587-2601. https://doi.org/10.1111/evo.12758

Eiten, G. (1972). The Cerrado vegetational of Brazil. *Botânical Review*, **38**, 201-341.

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

Excoffier, L., & Lischer, H. E. L. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. *Molecular Ecology Resources*, **10**(3), 564-567.
https://doi.org/10.1111/j.1755-0998.2010.02847.x

Feder, J. L., Egan, S. P., & Nosil, P. (2012). The genomics of speciation-with-gene-flow. *Trends in Genetics*, **28**, 342-350. https://doi.org/10.1016/j.tig.2012.03.009

Felsenstein, J. (1985). Confidence limits on phylogenies: An approach using the bootstrap. *Evolution*, **39**(4), 783-791. https://doi.org/10.1111/j.1558-5646.1985.tb00420.x

Fick, S. E., & Hijmans, R. J. (2017). WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. *International Journal of Climatology*, **37**(12), 4302-4315. https://doi.org/10.1002/joc.5086

Field, D. L., Ayre, D. J., Whelan, R. J., & Young, A. G. (2010). Patterns of hybridization and asymmetrical gene flow in hybrid zones of the rare *Eucalyptus aggregata* and common *E. rubida*. *Heredity*, **106**(5), 841-853. https://doi.org/10.1038/hdy.2010.127

Fitzpatrick, B. M., & Irwin (2012). Estimating ancestry and heterozygosity of hybrids using molecular markers. BMC Evolutionary Biology, 12(1), 1-14. https://doi.org/10.1186/1471-2148-12-131

Foster, S. A., McKinnon, G. E., Steane, D. A., Potts, B. M., & Vaillancourt, R. E. (2007). Parallel evolution of dwarf ecotypes in the forest tree *Eucalyptus globulus*. *New Phytologist*, **175**(2), 370-380. https://doi.org/10.1111/j.1469-8137.2007.02077.x

Francisco, P. M., Mori, G. M., Alves, F. M., Tambarussi, E. V., & de Souza, A. P. (2018). Population genetic structure, introgression, and hybridization in the genus Rhizophora along the Brazilian coast. *Ecology and Evolution*, **8**, 3491-3504. https://doi.org/10.1002/ece3.3900

Gompert, Z., & Buerkle, C. A. (2010). introgress: A software package for mapping components of isolation in hybrids. *Molecular Ecology Resources*, **10**(2), 378-384. https://doi.org/10.1111/j.1755-0998.2009.02733.x

Gompert, Z., & Buerkle, C. A. (2016). What, if anything, are hybrids: Enduring truths and challenges associated with population structure and gene flow. *Evolutionary Applications*, **9**(7), 909-923.
https://doi.org/10.1111/eva.12380

Gompert, Z., Mandeville, E. G., & Buerkle, C. A. (2017). Analysis of population genomic data from hybrid zones. *Annual Review of Ecology, Evolution, and Systematics*, **48**(1), 207-229. https://doi.org/10.1146/annurev-ecolsys-110316-022652

Goslee, S., & Urban, D. (2007). The ecodist package for dissimilarity-based analysis of ecological data. *Journal of Statistical Software*, **22**(7), 1-19.

Goudet, J. (2002). FSTAT (version 2.9. 3.2): A program to Estimate and Test gene Diversities and Fixation Indices. http://www.UnilCh/izea/softw are/fstat.html

Goulart, M. F., De Lemos-Filho, J. P., & Lovato, M. B. (2005). Phenological variation within and among populations of *Plathymenia reticulata* in Brazilian Cerrado, the Atlantic Forest and transitional sites. *Annals of Botany*, **96**(3), 445-455. https://doi.org/10.1093/aob/mci193

Goulart, M. F., De Lemos-Filho, J. P., & Lovato, M. B. (2006). Variability in fruit and seed morphology among and within populations of *Plathymenia* (Leguminosae - Mimosoideae) in areas of the Cerrado, the Atlantic Forest, and transitional sites. *Plant Biology*, **8**(1), 112-119. https://doi.org/10.1055/s-2005-865964

Goulart, M. F., Lovato, M. B., de Vasconcellos Barros, F., Valladares, F., & De Lemos-Filho, J. P. (2011). Which extent is plasticity to light involved in the ecotypic differentiation of a tree species from savanna and forest? *Biotropica*, **43**(6), 695-703. https://doi.org/10.1111/j.1744-7429.2011.00760.x

Gramlich, S., Wagner, N. D., & Hörandl, E. (2018). RAD-seq reveals genetic structure of the F2-generation of natural willow hybrids (Salix L.) and a great potential for interspecific introgression. BMC Plant Biology, 18(1). https://doi.org/10.1186/s12870-018-1552-6

Hanson, D., Moore, J. S., Taylor, E. B., Barrett, R. D. H., & Hendry, A. P. (2016). Assessing reproductive isolation using a contact zone between parapatric lake-stream stickleback ecotypes. *Journal of Evolutionary Biology*, **29**(12), 2491-2501. https://doi.org/10.1111/jeb.12978

Heringer, E. P. (1956). O gênero *Plathymenia*. *Anais Da Sociedade Botânica do Brasil*, **7**, 55-64.

Hoffmann, W. A., & Franco, A. C. (2003). Comparative growth analysis of tropical forest and savanna woody plants using phylogenetically independent contrasts. *Journal of Ecology*, **91**(3), 475-484. https://doi.org/10.1046/j.1365-2745.2003.00777.x

Hoffmann, W. A., Orthen, B., & Vargas Do Nascimento, P. K. (2003). Comparative fire ecology of tropical savanna and forest trees. *Functional Ecology*, **17**(6), 720-726. https://doi.org/10.1111/j.1365-2435.2003.00796.x

Hubisz, M. J., Falush, D., Stephens, M., & Pritchard, J. K. (2009). Inferring weak population structure with the assistance of sample group
information. *Molecular Ecology Resources*, 9(5), 1322-1332. https://doi.org/10.1111/j.1755-0998.2009.02591.x

Jakobsson, M., & Rosenberg, N. A. (2007). CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. *Bioinformatics*, 23(14), 1801-1806. https://doi.org/10.1093/bioinformatics/btm233

Jombart, T. (2008). Adegenet: A R package for the multivariate analysis of genetic markers. *Bioinformatics*, 24(11), 1403-1405. https://doi.org/10.1093/bioinformatics/btn129

Jombart, T., Devillard, S., & Balloux, F. (2010). Discriminant analysis of principal components: A new method for the analysis of genetically structured populations. *BMC Genetics*, 11, 94. https://doi.org/10.1186/1471-2156-11-94

Jombart, T., Devillard, S., Dufour, A. B., & Pontier, D. (2008). Revealing cryptic spatial patterns in genetic variability by a new multivariate method. *Heredity*, 101(1), 92-103. https://doi.org/10.1038/hdy.2008.34

Kierepka, E. M., & Latch, E. K. (2015). Performance of partial statistics in individual-based landscape genetics. *Molecular Ecology Resources*, 15(3), 512-525. https://doi.org/10.1111/1755-0998.12332

Klink, C. A., & Machado, R. B. (2005). Conservation of the Brazilian Cerrado. *Conservation Biology*, 19(3), 707-713. https://doi.org/10.1111/j.1523-1739.2005.00702.x

Kopp, A., & Frank, A. K. (2005). Speciation in progress? A continuum of reproductive isolation in *Drosophila bipectinata*. *Genetica*, 125(1), 55-68. https://doi.org/10.1007/s10709-005-4787-8

Lacerda, D. R., de Lemos-Filho, J. P., Acedo, M. D. P., & Lovato, M. B. (2002). Molecular differentiation of two vicariant neotropical tree species, *Plathymenia foliolosa* and *P. reticulata* (Mimosoideae), inferred using RAPD markers. *Plant Systematics and Evolution*, 235, 67-77. https://doi.org/10.1007/s00606-002-0227-8

Lambers, H., Brundrett, M. C., Raven, J. A., & Hopper, S. D. (2011). Plant mineral nutrition in ancient landscapes: High plant species diversity on infertile soils is linked to functional diversity for nutritional strategies. *Plant and Soil*, 348(1-2), 7-27. https://doi.org/10.1007/s11104-011-0977-6

Leal, B. S. S., Palma da Silva, C., & Pinheiro, F. (2016). Phylogeographic studies depict the role of space and time scales of plant speciation in a highly diverse neotropical region. *Critical Reviews in Plant Sciences*, 35(4), 215-230. https://doi.org/10.1080/07352689.2016.1254494

Legendre, P., & Legendre, L. (2014). *Numerical Ecology*. Developments in Environmental Modelling (3 ed., Vol. 24). Elsevier.

Leutner, B., & Horning, N. (2019). *RStoolbox: Tools for Remote Sensing Data Analysis*. R package version 0.2.6. https://CRAN.R-project.org/package=RStoolbox

Lexer, C., Fay, M. F., Joseph, J. A., Nica, M. S., & Heinze, B. (2005). Barrier to gene flow between two ecologically divergent *Populus* species, *P. alba* (white poplar) and *P. tremula* (European aspen): The role of ecology and life history in gene introgression. *Molecular Ecology*, 14(4), 1045-1057. https://doi.org/10.1111/j.1365-294X.2005.02469.x

Li, Y. L., & Liu, J. X. (2018). StructureSelector: A web-based software to select and visualize the optimal number of clusters using multiple methods. *Molecular Ecology Resources*, 18, 176-177. https://doi.org/10.1111/1755-0998.12719

Li, Y., Tada, F., Yamashiro, T., & Maki, M. (2016). Long-term persisting hybrid swarm and geographic difference in hybridization pattern: Genetic consequences of secondary contact between two *Vincetoxicum* species (Apocynaceae-Asclepiadoideae). *BMC Evolutionary Biology*, 16(1), 1-13. https://doi.org/10.1186/s12862-016-0587-2

Lindtke, D., Gompert, Z., lexer, C., & Buerkle, C. A. (2014). Unexpected ancestry of Populus seedlings from a hybrid zone implies a large role for postzygotic selection in the maintenance of species. *Molecular Ecology*, 23(17), 4316-4330. https://doi.org/10.1111/mec.12759

Lorenzi, H. (1992). Árvores Brasileiras. Manual de identificação e cultivo de plantas arbóreas nativas do Brasil (pp. 348). Plantarum Ltda.

Maguilla, E., & Escudero, M. (2016). Cryptic species due to hybridization: A combined approach to describe a new species (carex: Cyperaceae). *PLoS One*, 11(12), 1-24. https://doi.org/10.1371/journal.pone.0166949

Mallet, J. (2005). Hybridization as an invasion of the genome. *Trends in Ecology and Evolution*, 20(5), 229-237. https://doi.org/10.1016/j.tree.2005.02.010

Maracahipes, L., Carlucci, M. B., Lenza, E., Marimon, B. S., Marimon, B. H., Guimarães, F. A. G., & Cianciaruso, M. V. (2018). How to live in contrasting habitats? Acquisitive and conservative strategies emerge at inter- and intraspecific levels in savanna and forest woody plants. *Perspectives in Plant Ecology, Evolution and Systematics*, 34(July), 17-25. https://doi.org/10.1016/j.ppees.2018.07.006

Meneses, M. E. N. S., Costa, M. L., & Behling, H. (2013). Late Holocene vegetation and fire dynamics from a savanna-forest ecotone in Roraima state, northern Brazilian Amazon. *Journal of South American Earth Sciences*, 42, 17-26. https://doi.org/10.1016/j.jsames.2012.10.007

Muniz, A. C., Lemos-Filho, J. P., Souza, H. A., Marinho, R. C., Buzatti, R. S., Heuertz, M., & Lovato, M. B. (2020). The protected tree *Dimorphandra wilsonii* (Fabaceae) is a population of inter-specific hybrids: Recommendations for conservation in the Brazilian Cerrado/Atlantic Forest ecotone. *Annals of Botany*, 126, 191-203. https://doi.org/10.1093/aob/mcaa066

Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., & Kent, J. (2000). Biodiversity hotspots for conservation priorities. *Nature*, 403(6772), 853-858. https://doi.org/10.1038/35002501

Nadeau, S., Meirmans, P. G., Aitken, S. N., Ritland, K., & Isabel, N. (2016). The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines. *Ecology and Evolution*, 6(24), 8649-8664. https://doi.org/10.1002/ece3.2550

Nei, M. (1972). Genetic distance between populations. *The American Naturalist*, 106(949), 283-292. http://www.jstor.org/stable/e/2459777

Nosil, P., Egan, S. P., & Funk, D. J. (2008). Heterogeneous genomic differentiation between walking-stick ecotypes: "Isolation by adaptation" and multiple roles for divergent selection. *Evolution*, 62(2), 316-336. https://doi.org/10.1111/j.1558-5646.2007.00299.x

Nosil, P., Funk, D. J., & Ortiz-Barrientos, D. (2009). Divergent selection and heterogeneous genomic divergence. *Molecular Ecology*, 18(3), 375-402. https://doi.org/10.1111/j.1365-294X.2008.03946.x

Novaes, R. M. L., De Lemos-Filho, J. P., Ribeiro, R. A., & Lovato, M. B. (2010). Phylogeography of *Plathymenia reticulata* (Leguminosae) reveals patterns of recent range expansion towards northeastern Brazil and southern Cerrados in Eastern Tropical South America. *Molecular Ecology*, 19(5), 985-998. https://doi.org/10.1111/j.1365-294X.2010.04530.x

Novaes, R. M. L., Ribeiro, R. A., Lemos-Filho, J. P., & Lovato, M. B. (2013). Concordance between phylogeographical and biogeographical patterns in the Brazilian Cerrado: Diversification of the endemic tree *Dalbergia miscolobium* (Fabaceae). *PLoS One*, 8(12), e82198. https://doi.org/10.1371/journal.pone.0082198

Novaes, R. M. L., Rodrigues, J. G., & Lovato, M. B. (2009). An efficient protocol for tissue sampling and DNA isolation from the stem bark of Leguminosae trees. *Genetics and Molecular Research*, 8(1), 86-96. https://doi.org/10.4238/vol8-1gmr542

Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., & Mcglinn, D. (2020). vegan: Community Ecology Package. R package version 2.5-7. https://CRAN.R-project.org/package=vegan

Oliveira, F. A., Tarazi, R., Menezes, I. P. P., van den Berg, C., Tsai, S. M., & Gaiotto, F. A. (2012). Microsatellite markers for *Plathymenia reticulata* (Leguminosae). *American Journal of Botany*. https://doi.org/10.3732/ajb.1200051

Oliveira-Filho, A. T., Fontes, M. A. L. (2000). Patterns of floristic differentiation among Atlantic Forests in Southeastern Brazil and
the influence of climate. *Biotropica*, 32(4b), 793-810. https://doi.org/10.1111/j.1744-7429.2000.tb00619.x

Padgham, M., & Sumner, M. D. (2021). geodist: Fast, Dependency-Free Geodesic Distance Calculations. R package version 0.0.7. https://CRAN.R-project.org/package=geodist

Paun, O., Forest, F., Fay, M. F., & Chase, M. W. (2009). Hybrid speciation in angiosperms: Parental divergence drives ploidy. *New Phytologist*, 182, 507-518. https://doi.org/10.1111/j.1469-8137.2009.02767.x

Pritchard, J. K., Stephens, M., & Donnelly, P. (2000). Inference of population structure using multilocus genotype data. *Genetics*, 155, 945-959. https://doi.org/10.1111/j.1471-8286.2007.01758.x

Puechmaille, S. J. (2016). The program structure does not reliably recover the correct population structure when sampling is uneven: Subsampling and new estimators alleviate the problem. *Molecular Ecology Resources*, *16*, 608-627. https://doi.org/10.1111/1755-0998.12512

R Core Team. (2020). *R: A Language and Environment for Statistical Computing*. R Foundation for Statistical Computing. https://www.R-project.org/

Reich, P. B., Wright, I. J., Cavender-Bares, J., Craine, J. M., Oleksyn, J., Westoby, M., & Walters, M. B. (2003). The evolution of plant functional variation: Traits, spectra, and strategies. *International Journal of Plant Sciences*, *164*(S3), S143-S164. https://doi.org/10.1086/374368

Resende-Moreira, L. C., Ramos, A. C. S., Scliar, M. O., Silva, R. M., Azevedo, V. C. R., Ciampi, A. Y., Lemos-Filho, J. P., & Lovato, M. B. (2017). Gene flow between vicariant tree species: Insights into savanna-forest evolutionary relationships. *Tree Genetics and Genomes*, *13*(2), 36. https://doi.org/10.1007/s11295-017-1117-0

Rhymer, J. M., & Simberloff, D. (1996). Extinction by hybridization and introgression. *Annual Review of Ecology and Systematics*, *27*(1), 83-109. https://doi.org/10.1146/annurev.ecolsys.27.1.83

Ribeiro, R. A., Lemos-Filho, J. P., Ramos, A. C. S., & Lovato, M. B. (2011). Phylogeography of the endangered rosewood *Dalbergia nigra* (Fabaceae): Insights into the evolutionary history and conservation of the Brazilian Atlantic Forest. *Heredity*, *106*(1), 46-57. https://doi.org/10.1038/hdy.2010.64

Ribeiro, M. C., Metzger, J. P., Martensen, A. C., Ponzoni, F. J., & Hirota, M. M. (2009). The Brazilian Atlantic Forest: How much is left, and how is the remaining forest distributed? Implications for conservation. *Biological Conservation*, *142*(6), 1141-1153. https://doi.org/10.1016/j.biocon.2009.02.021

Ribeiro, P. C. C., Souza, M. L., Muller, L. A. C., Ellis, V. A., Heuertz, M., De Lemos-Filho, J. P., & Lovato, M. B. (2016). Climatic drivers of leaf traits and genetic divergence in the tree *Annona crassiflora*: A broad spatial survey in the Brazilian savannas. *Global Change Biology*, *22*(11), 3789-3803. https://doi.org/10.1111/gcb.13312

Rieseberg, L. H., & Carney, S. E. (1998). Plant hybridization. *New Phytologist*, *140*(4), 599-624. https://doi.org/10.1046/j.1469-8137.1998.00315.x

Rifkin, J. L., Castillo, A. S., Liao, I. T., & Rausher, M. D. (2019). Gene flow, divergent selection and resistance to introgression in two species of morning glories (Ipomoea). *Molecular Ecology*, *28*(7), 1709-1729. https://doi.org/10.1111/mec.14945

Saitou, N., & Nei, M. (1987). The neighbor-joining method: A new method for reconstructing phylogenetic trees. *Molecular Biology and Evolution*, *4*(4), 406-425. https://doi.org/10.1093/oxfordjournals.molbev.a040454

Schmitt, S., Tysklind, N., Hérault, B., & Heuertz, M. (2021). Topography drives microgeographic adaptations of closely related species in two tropical tree species complexes. *Molecular Ecology*, *30*(20), 5080-5093. https://doi.org/10.1111/mec.16116

Servedio, M. R. (2016). Geography, assortative mating, and the effects of sexual selection on speciation with gene flow. *Evolutionary Applications*, *9*(1), 91-102. https://doi.org/10.1111/eva.12296

Sexton, J. P., Hangartner, S. B., & Hoffmann, A. A. (2014). Genetic isolation by environment or distance: Which pattern of gene flow is most common? *Evolution*, *68*(1), 1-15. https://doi.org/10.1111/evo.12258

Sexton, J. P., Hufford, M. B., Bateman, A. C., Lowry, D. B., Meimberg, H., Strauss, S. Y., & Rice, K. J. (2016). Climate structures genetic variation across a species' elevation range: A test of range limits hypotheses. *Molecular Ecology*, *25*(4), 911-928. https://doi.org/10.1111/mec.13528

Silva, A. R., Resende-Moreira, L. C., Carvalho, C. S., Lanes, E. C. M., Ortiz-Vera, M. P., Viana, P. L., & Jaffé, R. (2020). Range-wide neutral and adaptive genetic structure of an endemic herb from Amazonian Savannas. *AoB PLANTS*, *12*(1), 1-11. https://doi.org/10.1093/aobpl/a pla003

Silva, J. F., Farinas, M. R., Felfili, J. M., & Klink, C. A. (2006). Spatial heterogeneity, land use and conservation in the Cerrado region of Brazil. *Journal of Biogeography*, *33*(3), 536-548. https://doi.org/10.1111/j.1365-2699.2005.01422.x

Silva, M. C., Teodoro, G. S., Bragion, E. F. A., & van den Berg, E. (2019). The role of intraspecific trait variation in the occupation of sharp forest-savanna ecotypes. *Flora*, *253*(March), 35-42. https://doi.org/10.1016/j.flora.2019.03.003

Soltis, P. S., & Soltis, D. E. (2009). The role of hybridization in plant speciation. *Annual Review of Plant Biology*, *60*(1), 561-588. https://doi.org/10.1146/annurev.arplant.043008.092039

Souza, H., Muller, L., Brandão, R. L., & Lovato, M. B. (2012). Isolation of high quality and polysaccharide-free DNA from leaves of *Dimorphandra mollis* (Leguminosae), a tree from the Brazilian Cerrado. *Genetics and Molecular Research: GMR*, *11*(1), 756-764. https://doi.org/10.4238/2012.March.22.6

Takezaki, N., Nei, M., & Tamura, K. (2014). POPTREEW: Web version of POPTREEW for constructing population trees from allele frequency data and computing some other quantities. *Molecular Biology and Evolution*, *31*(6), 1622-1624. https://doi.org/10.1093/molbev/msu093

Todesco, M., Pascual, M. A., Owens, G. L., Ostevik, K. L., Moyers, B. T., Hübner, S., Heredia, S. M., Hahn, M. A., Caseys, C., Bock, D. G., & Rieseberg, L. H. (2016). Hybridization and extinction. *Evolutionary Applications*, *9*(7), 892-908. https://doi.org/10.1111/eva.12367

Toledo, M. M., Paiva, E. A. S., Lovato, M. B., & De Lemos-Filho, J. P. (2012). Stem radial increment of forest and savanna ecotypes of a Neotropical tree: Relationships with climate, phenology, and water potential. *Trees - Structure and Function*, *26*(4), 1137-1144. https://doi.org/10.1007/s00468-012-0690-y

Van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M., & Shipley, P. (2004). MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. *Molecular Ecology Notes*, *4*(3), 535-538. https://doi.org/10.1111/j.1471-8286.2004.00684.x

Wagner, H. H., & Fortin, M.-J. (2012). A conceptual framework for the spatial analysis of landscape genetic data. *Conservation Genetics*, *14*(2), 253-261. https://doi.org/10.1007/s10592-012-0391-5

Wang, I. J., & Bradburd, G. S. (2014). Isolation by environment. *Molecular Ecology*, *23*(23), 5649-5662. https://doi.org/10.1111/mec.12938

Warwick, M. C., & Lewis, G. P. (2003). Revision of *Plathymenia* (LEGUMINOSAE - MIMOSOIDEAE). *Edinburgh Journal of Botany*, *60*(2), 111-119. https://doi.org/10.1017/S09604286030000106

Weir, B. S., & Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. *Evolution*, *36*(6), 1358-1370.

Whitney, K. D., Randell, R. A., & Rieseberg, L. H. (2006). Adaptive introgression of herbivore resistance traits in the weeds sunflower *Helianthus annuus*. *American Naturalist*, *167*(6), 794-807. https://doi.org/10.1086/504606

Wickam, H. (2009). *ggplot2: Elegant graphics for data analysis*. Springer. http://HadCoNz/Ggplot2/Book. Accessed: August, 5, 2014.
Wilson, G. A., & Rannala, B. (2003). Bayesian inference of recent migration rates using multilocus genotypes. *Genetics*, *163*(3), 1177-1191.
https://doi.org/10.1093/genetics/163.3.1177

World Conservation Monitoring Centre (1998). *Plathymenia foliolosa*
(Vol. 8235). Vol. 8235

Wu, C.-I.-I. (2001). The genic view of the process of speciation. *Journal of Evolutionary Biology*, *14*(6), 851-865. https://doi.org/10.1046/j.1420-9101.2001.00335.x

Yamamoto, M., Takahashi, D., Horita, K., & Setoguchi, H. (2019). Speciation and subsequent secondary contact in two edaphic endemic primroses driven by Pleistocene climatic oscillation. *Heredity*, *124*(1), 93-107. https://doi.org/10.1038/s41437-019-0245-8

Zalapa, J. E., Brunet, J., & Guries, R. P. (2010). The extent of hybridization and its impact on the genetic diversity and population structure of an invasive tree, *Ulmus pumila* (Ulmaceae). *Evolutionary Applications*, *3*(2), 157-168. https://doi.org/10.1111/j.1752-4571.2009.00106.x

SUPPORTING INFORMATION

Additional supporting information may be found in the online
version of the article at the publisher's website.

**How to cite this article:** Muniz, A. C., Pimenta, R. J. G., Cruz, M. V., Rodrigues, J. G., Buzatti, R. S. D. O., Heuertz, M., Lemos-Filho, J. P., & Lovato, M. B. (2022). Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation. *Ecology and Evolution*, *12*, e8540. https://doi.org/10.1002/ece3.8540