Genetic and reproductive evidence of incomplete isolations barriers between Pterodon emarginatus and P. pubescens (Leguminosae, Papilionoideae)Abstract
Pterodon emarginatus Vogel and P. pubescens (Benth.) Benth. are phylogenetically related trees that compose a clade of sister species abundant in the Brazilian Savanna. Despite their morphological differences, some individuals with intermediate morphological characteristics have been reported, indicating the formation of interspecific hybrids. This study proposed to evaluate the genetic structure and controlled pollination of individuals of P. pubescens and P. emarginatus in areas of sympatry with the presence of putative hybrids. For this purpose, we genotyped seven microsatellite loci from 61 individuals collected from four apparent contact zones between P. pubescens and P. emarginatus. Controlled pollination experiments were performed on 4,133 flowers from six trees of P. emarginatus and five of P. pubescens. We observed two genetic clusters (k= 2) that corroborate the divergence between P. pubescens and P. emarginatus. The individual genetic assignment showed evidence of natural hybridization between P. pubescens and P. emarginatus. The genetic assignment did not fully support the visual description of the diagnoses of individual morphological characteristics. Controlled interspecific pollination generated fruit and seed production, possibly indicating the absence of a reproductive barrier at the pollination and ovule fertilization levels between these species. Our results enlarge the understanding of the diversification process of Pterodon species.
Keywords: Controlled pollination; Microsatellite markers; Neotropical tree; Speciation; sucupira-branca
Introduction
The evolutionary consequences of hybridization range from the termination of reproductive isolation between species (Servedio & Noor 2003) to the homogenization of both species’ gene pools (Currat et al. 2008; Marques et al. 2014). The outcome of these hybrids depends on their genetic combinations, demography, population structure, and reproductive viability (Yan et al. 2017). It is not clear which ecological and evolutionary factors drive one outcome of hybridization and not the other (Turchetto et al. 2022).
The presence of both biotic and abiotic factors suggests that the biomes of Savannas and grasslands might shelter a relatively abundant. Since Savannas are relatively young Neotropical biomes, it is likely that their hybrids reflect a minimal divergence among the endemic species that resulted from recent radiation (Schley et al. 2022; Turchetto et al. 2022). Moreover, these Neotropical biomes have widely distributed groups that present clear genetic evidence of hybridization (e.g., Ribeiro et al. 2018). However, there are only a few studies examining the impacts of hybridization on these biomes (Schley et al. 2022; Turchetto et al. 2022).
The Brazilian Savanna biome (also known as Cerrado) is considered one of the greatest biodiversity hotspots in the world (Myers et al. 2000), sheltering two non-endemic tree species belonging to the genus Pterodon (Leguminosae): P. emarginatus and P. pubescens. In addition to sharing the same popular name (“sucupira-branca”), the two species have resembling genetic variations (Lima 2019) and chromosomal structure similarities (Bandel 1974; Coleman & Demenezes 1980; Albernaz 2020), thus being phylogenetically recognized as closely related sibling species (Cardoso et al. 2013). In morphological terms, P. emarginatus (Fig. 1AD) has a glabrous leaf, glabrous or glabrescent rachis containing 4-10 leaflets, truncated to the strongly emarginate apex, and violaceous flowers (Rocha 2006). Pterodon pubescens (Fig. 1BE) has leaves containing 6-19 elliptic leaflets (usually 11-13), slightly ovate, pubescent on both sides, retuse to rounded apex, densely pubescent rachis, and pale pinkish almost white-to-dark pink flowers (Rocha 2006). Both mating systems are allogamous, and the flowers are pollinated by Bombus atratus and Apis mellifera bees, reaching flowering peaks in September, and the dispersing agent is the wind (da Silva Júnior 2012).
Trees of Pterodon emarginatus (A), P. pubescens (B), and potential interspecific hybrid (C) in the natural environment, followed by their respective flowers (D, E, and F). Elaborated by the authors, DMS Rocha & VFM Lima.
Our field observations revealed that both Pterodon species have a disjunct distribution. In addition to occurring in the Cerrado areas further South, the Pterodon pubescens is also present in the states of São Paulo, Mato Grosso do Sul, the south of the states of Mato Grosso, Minas Gerais, Goiás, and in the south of the Federal District. In turn, P. emarginatus occupies the following areas further North: north of the states of Mato Grosso, Federal District, Minas Gerais, and the north center of the states of Goiás, Tocantins, Maranhão, Bahia. Individuals showing intermediate characteristics appeared (SL Lima unpubl. res.) in the regions with overlapping distribution.
In a narrow range in the Distrito Federal (Brazil), Rocha (2006) detected individuals showing intermediate morphological features between both species (Fig. 1CF), in addition to some genetic hybrids (via RAPD). Such a finding reveals the possibility of interspecific permeability of species borders. In turn, a phylogeography study indicated that P. pubescens and P. emarginatus underwent recent divergence; therefore, the barriers of reproductive isolation between them might not be absolute (Lima 2019). Herein, our field observations detected four apparent contact zones between the two species, with flowering adults. Thus, this species pair represents a great study model for reproductive isolation, taxonomic cohesion, and natural hybridization. The reproductive isolation of related species (e.g., Johnson et al. 2015; Kuligowska et al. 2015; Wang et al. 2015) is often associated with post-pollination mechanisms, and might occur in the following situations: if interspecific pollen grains do not germinate on the stigma; pollen tubes do not grow properly in the style, absence of fertilization, and if the embryo and/or endosperm does not develop well (Rieseberg & Carney 1998; Kuligowska et al. 2015). However, such barriers might have occurred if seeds are produced after interspecific cross-pollination. In this sense, controlled crossings might help understand the hybridization between these species.
Melo et al. (2022) used P. pubescens sequence data to develop a set of microsatellite markers (SSR) for studying the genetics of both populations. The author’s proposal was unprecedented in this research line and indicated a high genetic diversity (higher in P. pubescens) with significant intrapopulation inbreeding.
More recently, methods of hybridization detection have estimated the degree of genetic variation shared between adjacent or co-occurring species, including calculating “hybrid indices” to assess the proportion of ancestry inherited from each parent (Schley et al. 2022). In this regard, population-level sampling strategies and molecular markers approaches (e.g., microsatellites, AFLPs, and SNPs from RADSeq) are used (Schley et al. 2022). Studies using molecular markers that can differentiate between closely related species and their hybrids help enlarge the knowledge on the speciation genetic basis and the effects of interspecific crossing on species integrity. Additionally, these studies provide further information on the evolutionary history of the species involved.
Despite the studies addressing these trees of wide economic and cultural potential (e.g., Mors et al. 1967; Bustamante et al. 2010; Basting et al. 2019; Kleinubing et al. 2022), the understanding of their intra- and interspecific genetic variations is yet to be investigated. Therefore, this study aimed to evaluate the genetic composition and controlled pollination of individuals from apparent contact zones of P. pubescens and P. emarginatus containing putative hybrids. Thereby, we sought to answer the three following main questions: (I) Are there genetic groups that support the divergence between P. pubescens and P. emarginatus? (II) Does the genetic composition of individuals characterized as putative hybrids in the field corroborate the hypothesis of natural hybridization between P. pubescens and P. emarginatus? and (III) Does interspecific crossing produce fruits with well-formed seeds?
Materials and methods
Material sampling and identification
We collected young leaves from 61 individuals in four natural sites presenting an apparent sympatry between Pterodon pubescens Benth. (Benth.) and Pterodon emarginatus Vogel (Fig. 2). Each individual collected was identified in the field (Table 1) following the visual description of diagnoses morphological characters (Rocha 2006): P. pubescens (pale pinkish flowers, pubescent leaves, and leaflets with retuse to rounded apex), P. emarginatus (violaceous flowers, glabrous leaves, and leaflets with truncate to strongly emarginate apex), or potentials hybrids (individuals with intermediate or discordant characteristics of Rocha’s (2006) classification).
Representative map of the four natural sites presenting an apparent sympatry of P. pubescens and P. emarginatus. Table 1 shows the details of geographic coordinates and the number of individuals.
Description of the four natural sites in Brazil where apparent sympatry of Pterodon pubescens and P. emarginatus was reported, as well as the presence of potential hybrid (except in Itacaiú/GO). Field recording based on Rocha (2006).
Vouchers were taken to record each phenotype/species and were deposited at Unidade de conservação/PRPI - Herbário UFG (in Goiânia-GO, Brazil) -, as follows: P. pubescens, Herb Nº: 68.409 and 68.413; P. emarginatus, Herb Nº: 68.411, 68.414, 68.415, and 68.416, and the potential hybrid, Herb Nº: 68.410. The collected plant material was identified and packed in plastic bags containing spherical silica gel. Subsequently, these packages were transported to the Laboratório de Genética e Biodiversidade (LGBio) of the Universidade Federal de Goiás (Goiânia-GO, Brazil). After separating the material for DNA extraction, the remaining dried leaves were stored in a deep freezer (-80 ºC).
Molecular analysis
We have obtained total genomic DNA from leaf tissue using CTAB (cetyl-trimethyl ammonium bromide) extraction buffer, according to the extraction protocol by Doyle and Doyle (1987), modified by Ferreira and Grattapaglia (1996). DNA quality and concentration were assessed through electrophoresis on an agarose gel of 1%. The genotypes of the individuals were identified through seven nuclear microsatellite markers, developed for P. pubescens and transferred to P. emarginatus (Melo et al. 2022) (Table 2). PCR reactions were performed in a Veriti™ 96-Well Fast Thermal Cycler (Applied Biosystems®) under the PCR protocol and the thermal cycling conditions described by Melo et al. (2022): one cycle at 94 °C for five min; 30 cycles of 94 °C for one min, annealing temperature for one min (depending on the locus, Tab. 2), and 72 °C for one min; and 72 °C for 30 min to enforce 3′ Taq adenylation. The annealing temperature between P. emarginatus and P. pubescens presents a difference (of 2 ºC) in only one marker (PEM25). Therefore, for the potential hybrid individuals, we tested the PCR reaction with the annealing temperature optimized to P. pubescens. The PCR products were multiplexed, denatured, and size-fractionated using capillary electrophoresis on an ABI 3500 DNA Analyzer (Thermo Fisher Scientific®) with a LIZ (600) molecular size standard (Thermo Fisher Scientific®). The molecular data obtained are available in Table S2 of the Supplementary Material Table S2 - Description of the molecular data used in the genetic analyses of the 61 individuals of Pterodon in this study. Field record based on Rocha (2006). .
List of application multiplexes used for genotyping of the studied individuals using seven microsatellite markers developed by Melo et al. (2022). Forward primers were labeled with specific DYE (fluorochromes) for detection in capillary electrophoresis (genotyping). Potential hybrids were tested at the annealing temperatures of P. pubescens.
Statistical analyses
We estimated the genetic diversity parameters on the FSTAT 2.9.3.2 (Goudet 2002) and GDA Genetic Data Analysis 1.0 (Lewis & Zaykin 2001) programs. The following parameters were estimated: the observed heterozygosity (Ho), the expected heterozygosity under Hardy-Weinberg equilibrium (He), the number of alleles (A), and the average allelic richness based on the minimum sample size (A R ). In addition, we also estimated the probabilities of identity (PI) and paternity exclusion (Q) on the Identity v.1.0 software (Wagner & Sefc 1999) to verify the quality of the microsatellite set in individual discrimination.
The genetic structure of the individuals was analyzed based on the Bayesian clustering method implemented on the STRUCTURE 2.3.4 software (Pritchard et al. 2000). We performed 30 independent runs for K values ranging from one (no genetic structure) to four (assuming four natural sites) genetic clusters. Each run covered 1,000.000 Monte Carlo simulations through Markov chains (MCMC), with a burn-in of 10%, assuming the mixture model with correlated alleles between individuals (admixture model). We then inferred the appropriate K values to explain the data, as proposed by Evanno et al. (2005) and Puechmaille (2016) in the STRUCTURE SELECTOR platform (Li & Liu 2018). The Structure Selector platform allowed us to obtain the coancestry values of each individual belonging to a given cluster and the assignment probability values were assigned through multiple replicates for the same K value (Kopelman et al. 2015) on the CLUMPAK package, provided by the same platform.
Through simulations on the NEWHYBRIDS 1.1 software (Anderson & Thompson 2002), we also evaluated the posterior probability (PP) of each individual belonging to one of the following probability classes of Mendel’s law: Pure Parental A (P. pubescens), Pure Parental B (P. emarginatus), F1 hybrid, F2 hybrid, backcross A (with P. pubescens), and backcross B (with P. emarginatus). The assignment to different hybrid categories was systematized upon the individual belonging to any of the six classes considering a posterior probability ≥ 0.70. Individuals with PP ≥ 0.90 were considered to belong to pure parental lineages, whereas those without probability PP ≥ 0.70 for any of the six Mendel’s laws classes were considered of an uncategorized hybrid origin. This test was performed with no prior information on allele frequencies, “Jeffery’s like priors”, and was based on 100,000 Monte Carlo simulations through Markov chains (MCMC) after a burn-in period of 100,000 to ensure the convergence of chains and homogeneity of runs.
Controlled pollination experiments
We conducted controlled pollination experiments (in 2019 and 2021) in six trees of P. emarginatus from a fragment of Cerrado, at the Planaltina Campus of the Universidade de Brasília (Planaltina-DF, Brazil). For P. pubescens, the experiments used five trees from a fragment of Cerrado at the Campus of Universidade de Brasília (Brasília-DF, Brazil). We performed hand pollinations during the co-flowering period to verify whether intra- and heterospecific pollen depositions produced fruits with well-formed seeds, in addition to comparing the seed set rates between intraspecific and interspecific cross-pollinations. We also performed manual self-pollination and checked the occurrence of self-pollination and natural production of fruit and seed (control). The following individuals were chosen for our analyses: a monospecific population of P. pubescens with previously genotyped trees (Rocha 2006) and a monospecific population of P. emarginatus with trees from a contiguous area to that studied by Rocha (2006), identified as P. emarginatus (both by morphological characteristics and RAPD).
The following four treatments were applied to the bagged flowers of each species (see Table 5 for sample size): (1) heterospecific pollination (hybridization): stigmas of several flowers of one of the species were pollinated with pollen from the other species; (2) conspecific cross-pollination: stigmas of several flowers were pollinated with pollen from flowers of another distant plant of the same species (at least 30 m); (3) manual self-pollination: stigmas of several flowers were pollinated with pollen from flowers of the same tree; (4) spontaneous self-pollination: flowers in pre-anthesis and bottom of several inflorescences were counted and the inflorescences were bagged to check whether self-pollination occurs spontaneously; (5) Control: counting of the flowers in several inflorescences of different plants of both species for tagging and fruit counting. From experiments 1 to 4, the inflorescences of the flowers examined were bagged before and after hand pollination. In experiments 1 and 2, the flowers were emasculated before anthers opened, manually pollinated, and bagged again until the fruits were set. Subsequently, the fruits were counted, and the reproductive success rates were compared between treatments. The resulting well-formed fruits were collected and counted during the fruiting period.
We adopted general linear mixed models (GLMM) assuming a binomial error distribution (individual as a random factor) to test whether the fruiting rate differed between species and treatment; also in addition, the interaction between species and treatment was considered. Posthoc comparisons between treatment levels were established through the Tukey's method. All analyses were carried out in R version 4.1.2 (R Core Team 2021), using the packages lme4 v4.1-28 (Bates et al. 2014) and Emmeans v1.7.2 (Lenth 2022) for model fit and posthoc tests.
Results
Descriptive estimates of genetic variability
High estimate values of combined paternity exclusion probability (Q ~ 0.99985) were shown by the microsatellite set, thus demonstrating a strong ability of the markers to exclude potentially false paternity. We also found low values of combined probability of identity (I ~ 8.441*10-11), indicating that the microsatellite set presents a strong power of individual discrimination since the estimated PI is almost zero.
The PEM25 locus is the only one that has a different annealing temperature between species. The first PCR test amplified the locus PEM25 in the hybrid individuals at the annealing temperature of P. pubescens (54 ºC). Moreover, all loci were polymorphic, and 96 alleles (Table 3) were obtained, with an average allele richness per locus of 13.446. Higher values of allele number and heterozygosity were observed, thus indicating a greater content of marker information (Table 3).
The estimates of genetic variability parameters were analyzed from seven microsatellite loci for the sample of individuals registered as P. pubescens, P. emarginatus, and their potential hybrids. A: number of alleles; A R: average allelic richness based on minimum sample size (n= 53); He: expected heterozygosity under HWE; Ho: observed heterozygosity; PI: identity probability; Q: paternity exclusion probability; IC: combined probability of genetic identity; and QC: combined probability of paternity exclusion. All values are rounded to three decimal places.
Genetic structure
The Bayesian analysis using STRUCTURE suggested the following numbers of clusters: k= 2 by the Puechmaille (2016) method and k= 3 by the Evanno et al. (2005) method (Fig. 3). The k= 2 scenario summarized the information better than the k= 3 (Fig. 3). The k= 2 comprises the two parental species as contributors to the genetic structure of our samples, showing distinct cluster assignments between the P. emarginatus (PEM) and P. pubescens (PPU) individuals (coancestry values > 0.82). The k= 3 scenarios present a substructure of P. pubescens into two genetic clusters.
Bayesian clustering based on coancestry values, as follows: A) the formation of two genetic groups (k= 2) highlighted in dark pink and purple by the Puechmaille (2016) method and B) the formation of three genetic groups (k= 3) highlighted in dark pink, light pink, and purple by the Evanno et al. (2005) method. All data were grouped according to the field records, based on Rocha (2006), as follows: PEM - P. emarginatus, PPU - P. pubescens, and PSP - potential hybrids.
We identified the two parental species in k= 2 as contributors to the total gene pool sampled; therefore, individuals whose coancestry value was above or equal to 0.800 in any of the clusters were considered canonical species (for detailed coancestry values, see Table S1 in Supplementary Material Table S1 - Table with data compiled from the field record based on Rocha (2006) and the genetic assignment by coancestry of STRUCTURE and posterior probability (PP) in NEWHYBRIDS. For STRUCTURE (k= 2 obtained by Puechmaille (2016) method): Cluster 1 - P. pubescens, and Cluster 2 - P. emarginatus. For NEWHYBRIDS - the posterior probability of each individual belonging to one of the following probability classes of Mendel’s law: PPU - Pure Parental A (P. pubescens), PEM - Pure Parental B (P. emarginatus), F1HYB - F1 Hybrid, F2HYB - F2 Hybrid, BPPU - Backcross with P. pubescens, and BPEM - Backcross with P. emarginatus. ). Thus, we assigned 40 individuals as P. pubescens, 15 as P. emarginatus, and six as hybrids.
The genetic assignment by STRUCTURE showed differences (compared with field record) in the number of individuals assigned, as follows: from 34 to 40 for P. pubescens, from 10 to 15 for P. emarginatus, and from 17 to 6 for the hybrids. Based on six individuals with a clear hybrid identity, this initial evaluation revealed a genetic admixture between the species.
Hybrid assignment
The Bayesian analysis by NEWHYBRIDS provided a clear genetic distinction of the individuals (for detailed posterior probability values, see Table S1 in Supplementary Material Table S1 - Table with data compiled from the field record based on Rocha (2006) and the genetic assignment by coancestry of STRUCTURE and posterior probability (PP) in NEWHYBRIDS. For STRUCTURE (k= 2 obtained by Puechmaille (2016) method): Cluster 1 - P. pubescens, and Cluster 2 - P. emarginatus. For NEWHYBRIDS - the posterior probability of each individual belonging to one of the following probability classes of Mendel’s law: PPU - Pure Parental A (P. pubescens), PEM - Pure Parental B (P. emarginatus), F1HYB - F1 Hybrid, F2HYB - F2 Hybrid, BPPU - Backcross with P. pubescens, and BPEM - Backcross with P. emarginatus. ). In total, 87% of the individuals (n= 53) had significant values (≥ 0.700) of posterior probability (PP), thus belonging to either of the classes.
The genetic identification of the 61 individuals analyzed showed 13 P. emarginatus pure lineage (PP ≥ 0.930), 33 P. pubescens (PP ≥ 0.760), out of which only 28 were pure lineage (PP ≥ 0.900), and seven F2 hybrids (PP ≥ 0.773) (Fig. 4). The eight individuals wigth no significant probability (< 0.700) of belonging to any of the six classes were considered uncategorized hybrids, however with high genetic admixture values. Two of these uncategorized hybrids showed a higher posterior probability (PP) of being P. pubescens (0.452 < PP < 0.465), two P. emarginatus (0.520 < PP < 0.697), and the other four F2 hybrids (0.543 < PP < 0.633).
Posterior probability (PP) according to the category analyses on NEWHYBRIDS to verify hybridization between “sucupira-branca” species. Pure parent A (P. pubescens) is in pink; Pure parent B (P. emarginatus) is in purple. The size proportion of the categories in the individual bars indicates a higher assignment probability per class. All individuals were grouped according to the field records, based on Rocha (2006).
This fine genetic assignment provided 33 P. pubescens individuals, 13 P. emarginatus, and 15 hybrids. The NEWHYBRIDS results show the following differences in the number of individuals (compared with the field record): from 34 to 33 for P. pubescens, from 10 to 13 for P. emarginatus, and from 17 to 15 for the hybrids.
Agreement between hybrid assignment methods
We performed a combined evaluation of the field observation and genetic identification assignments (STRUCTURE and NEWHYBRIDS) (Table 4 and Table S1 in Supplementary Material Table S1 - Table with data compiled from the field record based on Rocha (2006) and the genetic assignment by coancestry of STRUCTURE and posterior probability (PP) in NEWHYBRIDS. For STRUCTURE (k= 2 obtained by Puechmaille (2016) method): Cluster 1 - P. pubescens, and Cluster 2 - P. emarginatus. For NEWHYBRIDS - the posterior probability of each individual belonging to one of the following probability classes of Mendel’s law: PPU - Pure Parental A (P. pubescens), PEM - Pure Parental B (P. emarginatus), F1HYB - F1 Hybrid, F2HYB - F2 Hybrid, BPPU - Backcross with P. pubescens, and BPEM - Backcross with P. emarginatus. ). Approximately 67% of the individuals (Six P. emarginatus, 31 P. pubescens, and four hybrids) were assigned to the same group in both the field and genetic identifications. About ~18% (n= 11) of the individuals showed agreement in the genetic identification but differing in the field records. These individuals presented the following behavior: six were genetically P. emarginatus and potentially hybrid according to the field records, two were genetically P. pubescens and potentially hybrid according to the field records, two were genetically hybrids and P. emarginatus according to the field records, and one was genetically P. emarginatus and P. pubescens according to the field record. The agreement between the field records and STRUCTURE reached ~70% of the individuals, while the field record and NEWHYBRIDS agreed in ~75% of the cases.
The classification was compiled from the field record according to the morphological diagnosis by Rocha (2006) and the genetic identification obtained by STRUCTURE’s coancestry and posterior probability (PP) in NEWHYBRIDS. Discrepancies in genetic identification between STRUCTURE and NEWHYBRIDS were classified as “Questionable hybrid” individuals. *Discrepancies between field records and genetic identification.
For both methods, only ~15% of the individuals (n= 9) showed different genetic assignments (STRUCTURE and NEWHYBRIDS), thus being considered “questionable hybrid” individuals (see Table 4 for detailed data). 10% (n= 6) of the total sampled individuals (n= 61) were genetically identified as hybrids by both software programs.
Controlled pollination experiment
All treatments of P. emarginatus generated ripe fruits (Table 5). The rate of fruit and seed set was higher for cross and interspecific pollinations comparing with self-pollination, which produced one manually self-pollinated fruit and one spontaneous fruit. Interspecific and cross-pollinations also produced fruits with seeds on Pterodon pubescens (Table 6). Self-pollination could not set any fruits, thus indicating that the species is completely allogamous and self-incompatible. The control treatment had very low rates of fruit set. The final fruiting rate was the same for both species (X2= 3.42, p= 0.18) but different between the treatments (X2= 44.44, p < 0.0001). Both species had the same outcome in the interspecific cross-pollination treatment (z= 1.33, p= 0.85) and the intraspecific cross-pollination treatment did not generate different results either for the P. pubescens (z= 2.12, p= 9.94) or P. emarginatus (z= 0.63, p= 0.99) hybrids.
Results of controlled pollination experiments conducted during the flowering of P. emarginatus in 2019 and 2021.
Results of controlled pollination experiments conducted during the flowering of P. pubescens in 2019 and 2021.
Discussion
Hybridization plays an important role in plant evolution (Abbott et al. 2016) and varied according to the taxa (Taylor & Larson 2019), presumably due to the biological differences between species pairs, such as timing and speciation mode (Payseur & Rieseberg 2016). Moreover, interspecific gene flow via natural hybrid zones might be a source of genetic variability for adaptation to environmental changes (Janes & Hamilton 2017; Turchetto et al. 2022).
Species that have diverged recently and occur in contact zones provide a unique opportunity to study the evolutionary process involved in speciation. Our study identified sympatric genetic individuals of P. emarginatus and P. pubescens and presenting potential hybrids. Our analysis was based on highly variable and informative microsatellite loci (Melo et al. 2022) with great potential for individual and species discrimination (paternity exclusion probability > 0.99 and almost no identity probability). Overall, the nuclear microsatellite loci demonstrated that the two “sucupira-branca” species are genetically distinct.
Based on our field observations and the literature, P. emarginatus seems to be distributed toward north and P. pubescens toward the south of the Cerrado, with a latitudinal contact range between them, thus forming contact zones (Rocha 2006). Rocha (2006) found individuals living in such a latitudinal range - our sampling site (Figure 1 and 2) - presenting morphological features and genetic composition that shown to be intermediary of the two species.
Contact zones opened space for the study of hybrid production and interspecific gene flow. Our results demonstrated that interspecific pollination experiments generated fruits and seeds, thus indicating that the two species have no reproductive barrier for pollen germination on the stigma, pollen tube development in the style, or ovule fertilization levels. Despite the absence of such barriers, the numbers of backcrosses and hybrids are relatively low; therefore, the barriers might be of an ecological nature. Our field observations revealed that some barriers might occur at the floral biology level since flowers of P. pubescens open earlier (around 6:00h) than those of P. emarginatus (around 8:00h), despite the time overlap between 8:00 to 12:00h. In addition, the two species might share pollinators, such as Bombus sp. Furthermore, we did not test whether there is a difference between the seed viability produced during interspecific and intraspecific cross-pollinations. There might be variations in seed germination, as well as seedling and plant development. In this sense, further studies should investigate their floral and pollination biology, seed germination, and seedling development to better understand how the two species avoid intercrossing.
The control test produced fewer fruits and seeds than manual cross-pollinations, which may be explained by a natural pollination deficit. The two cases of seed set after self-pollination of 206 flowers do not indicate the species are self-fertile because they can be a product of cross-pollen contamination. But further study might be conducted to better clarify these species’ self-compatibility systems.
Our results corroborate the hypothesis of interspecific genetic mixing between P. emarginatus and P. pubescens at the contact zone, which is often reported in the literature for other phylogenetically closely related species (Lorenz-Lemke et al. 2006; Taylor et al. 2014; Ley & Hardy 2017). Additionally, our analyses confirm a hybrid genetic status in six individuals at the contact zones, out of which four also have intermediate morphology, according to the field records (Table 4: ID 53, 54, 58, and 59).
We also detected that the morphological and genetic records indicated different assignments, as reported in other studies involving phylogenetically close plants (Teixeira et al. 2019; Schnitzler et al. 2020). In fact, the unpredictable phenotype expression in hybrids hampers the morphological hybrid diagnosis (Rieseberg et al. 1993; Teixeira et al. 2019). Since several genetic and environmental factors influence the hybrid phenotypic expression, genetic analyses using molecular markers should provide more robust diagnoses for hybrids presenting the phenotypic expression of a parental species, as our results show (López-Caamal & Tovar-Sánchez 2014).
We found a uniform hybridization pattern at the four sampling sites, given that they were limited to F2 hybrids (see NEWHYBRIDS results), indicating a cross between F1 hybrids. However, NEWHYBRIDS analyses may underestimate backcrosses when evaluating parental species that have recently diverged, thus classifying individuals as belonging to the pure lineage of either parent (Vähä & Primmer 2006). Such a scenario could indicate a lower frequency of backcrosses in our observations.
The genetic composition of each studied individual, with clear genetic distinction, confirm that P. emarginatus and P. pubescens should be considered independent species. However, there are no definite reproductive barriers between the two species, thus promoting typical hybrid genetic mixing. We suggest that the small genetic admixture in the canonical species (see STRUCTURE results) might be associated with size homoplasy since the alleles have the same size but the sequence is different (Estoup et al. 2002). Nonetheless, such a hypothesis should be verified through genotyping by sequencing the microsatellites.
An easy and successful interspecific crossing depends directly on the phylogenetic relationship between the species involved in hybridization, in addition to chromosomal homology. Otherwise, there might be incongruities/incompatibilities between the species’ genomes. Pterodon emarginatus and P. pubescens have the same karyotypic number 2n= 16, with small and morphologically similar chromosomes (Bandel 1974; Coleman & Demenezes 1980; Albernaz 2020). More recently, Albernaz (2020) studied the genomes of P. pubescens and P. emarginatus and found highly similar repetitive fractions. Thus, there are several elements in P. pubescens and P. emarginatus that facilitate interspecific hybridization, such as phylogenetic proximity, recent divergence, and vast cytomolecular similarity (both structural and numerical).
It is difficult to detect the processes that lead to species hybridization, even more so for those with a history of introgression and ancestral polymorphisms persisting in the speciation process, such as “sucupira-branca” (Lima 2019). The function of these hybrids in diversifying the studied species are unknown, as well as whether they reinforce reproductive isolation (due to the low viability of hybrids). However, hybrid formation does not necessarily have a specific function, it might simply be result from incomplete species barriers, such as in genetic material exchange (e.g., Zhang et al. 2016). In this sense, further studies should detail the reproductive biology and viability of hybrids of P. emarginatus and P. pubescens to better understand the role of this phenomenon in their evolutionary history.
Our study enlarges the knowledge on the diversification process of the “sucupira-branca” species. The microsatellite set indicated genetic groups that distinguish the P. pubescens and P. emarginatus individuals. Although approximately 10% of the studied individuals were genetically identified as hybrids, the intermediate morphological features per se do not allow us to identify an individual as hybrid. These contrasts found between the genetic and morphological attributions of P. pubescens, P. emarginatus, and their potential hybrids, highlight that integrative studies must investigate the ecological and evolutionary scenarios of these species deeper.
Although Rocha (2006) and Sonsin-Oliveira et al. (2022) has suggested the existence of hybridization between P. emarginatus and P. pubescens, our study is the first to reveal the reproductive viability of hybrid formation, in addition to genetic evidence of natural hybrids between these species. On the other hand, considering their evolutionary proximity, studies adopting more specific approaches, including the genetic identification of hybrid seeds or broad genomic analyses, should differentiate recent hybridization from the retention of ancestral polymorphism/ ancient hybridization.
Acknowledgments
We acknowledge MCTIC/CNPq (proc.431320/2018-7) and PRONEM/FAPEG (call Nº.07/2016, proc.201710267000539) for financial support. Current research is developed in the context of the INCT em Ecologia, Evolução e Conservação da Biodiversidade, supported by MCTIC/CNPq (proc.465610/2014-5) and FAPEG. VFML was supported by a fellowship from CNPq and CAPES. ARG has been supported by a fellowship from CNPq. TNS has been supported by productivity fellowships from CNPq. We are grateful to the collection team: LS Santos, MPC Telles, LJ Chaves, and RB Pinto. We thank LLima Bergamini for the statistical analysis done in the controlled pollination experiments.
References
-
Abbott RJ, Barton NH, Good JM. 2016. Genomics of hybridization and its evolutionary consequences. Molecular Ecology 25: 2325-2332. doi: 10.1111/mec.13685
» https://doi.org/10.1111/mec.13685 - Albernaz VB. 2020. Citogenômica de Pterodon pubescens e citogenética comparativa com P. emarginatus (Leguminosae). MSc Thesis, Universidade Federal de Goiás, Brasil.
-
Anderson EC, Thompson EA. 2002. A model-based method for identifying species hybrids using multilocus genetic data. Genetics 160: 1217-1229. doi: 10.1093/genetics/160.3.1217
» https://doi.org/10.1093/genetics/160.3.1217 -
Bandel G. 1974. Chromosome numbers and evolution in the Leguminosae. Caryologia 27: 17-32. doi: 10.1080/00087114.1974.10796558
» https://doi.org/10.1080/00087114.1974.10796558 -
Basting RT, Spindola HM, Sousa IMO et al 2019. Pterodon pubescens and Cordia verbenacea association promotes a synergistic response in antinociceptive model and improves the anti-inflammatory results in animal models. Biomedicine and Pharmacotherapy 112: 108693. doi: 10.1016/j.biopha.2019.108693
» https://doi.org/10.1016/j.biopha.2019.108693 -
Bates D, Mächler M, Bolker B, Walker S. 2014. Fitting linear mixed-effects models using lme4. Journal of Statistical Software: 1406.5823. doi: 10.48550/arXiv.1406.5823
» https://doi.org/10.48550/arXiv.1406.5823 -
Bustamante KGL, Lima ADF, Soares ML et al 2010. Avaliação da atividade antimicrobiana do extrato etanólico bruto da casca da sucupira branca (Pterodon emarginatus Vogel)-Fabaceae. Revista Brasileira de Plantas Medicinais 12: 341-345. doi: 10.1590/S1516-05722010000300012
» https://doi.org/10.1590/S1516-05722010000300012 -
Cardoso D, Pennington RT, Queiroz LP et al 2013. Reconstructing the deep-branching relationships of the papilionoid legumes. South African Journal of Botany 89: 58-75. doi: 10.1016/j.sajb.2013.05.001
» https://doi.org/10.1016/j.sajb.2013.05.001 -
Coleman JR, Demenezes EM. 1980. Chromosome numbers in Leguminosae from the State of Sao Paulo. Rhodora 82: 475-481. https://www.jstor.org/stable/23311939
» https://www.jstor.org/stable/23311939 -
Currat M, Ruedi M, Petit RJ, Excoffier L. 2008. The hidden side of invasions: Massive introgression by local genes. Evolution: International Journal of Organic Evolution 62: 1908-1920. doi: 10.1111/j.1558-5646.2008.00413.x
» https://doi.org/10.1111/j.1558-5646.2008.00413.x - da Silva Júnior MC. 2012. 100 Árvores do Cerrado - sentido restrito: Guia de campo. Brasília, Ed. Rede de Sementes do Cerrado.
- Doyle JJ, Doyle JL. 1987. A rapid DNA isolation method for small quantities of fresh tissues. Phytochem Bull 19: 11-15.
-
Estoup A, Jarne P, Cornuet JM. 2002. Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Molecular Ecology 11: 1591-1604. doi: 10.1046/j.1365-294X.2002.01576.x
» https://doi.org/10.1046/j.1365-294X.2002.01576.x -
Evanno G, Regnaut S, Goudet J. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Molecular Ecology 14: 2611-2620. doi: 10.1111/j.1365-294X.2005.02553.x
» https://doi.org/10.1111/j.1365-294X.2005.02553.x - Ferreira ME, Grattapaglia D. 1996. Introdução ao uso de marcadores moleculares em análise genética. Documento 20. Brası́lia, Embrapa Cenargen.
-
Goudet J. 2002. FSTAT: A program to estimate and test gene diversities and fixation indices (version 2.9.3.2). http://www2.unil.ch/popgen/softwares/fstat.htm 16 Sep. 2021.
» http://www2.unil.ch/popgen/softwares/fstat.htm -
Janes JK, Hamilton JA. 2017. Mixing it up: The role of hybridization in forest management and conservation under climate change. Forests 8: 237. doi: 10.3390/f8070237
» https://doi.org/10.3390/f8070237 -
Johnson MA, Price DK, Price JP, Stacy EA. 2015. Postzygotic barriers isolate sympatric species of Cyrtandra (Gesneriaceae) in Hawaiian montane forest understories. American Journal of Botany 102: 1870-1882. doi: 10.3732/ajb.1500288
» https://doi.org/10.3732/ajb.1500288 -
Kleinubing SA, Outuki PM, Santos ÉDS et al 2022. Stability Studies and the In Vitro Leishmanicidal Activity of Hyaluronic Acid-Based Nanoemulsion Containing Pterodon pubescens Benth. Oil. Colloids and Interfaces 6: 64. doi: 10.3390/colloids6040064
» https://doi.org/10.3390/colloids6040064 -
Kopelman NM, Mayzel J, Mattias J, Rosenberg NA, Mayrose I. 2015. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Molecular Ecology Resources 15: 1179-1191. doi: 10.1111/1755-0998.12387
» https://doi.org/10.1111/1755-0998.12387 -
Kuligowska K, Lütken H, Christensen B, Skovgaard I, Linde M, Winkelmann T, Müller R. 2015. Evaluation of reproductive barriers contributes to the development of novel interspecific hybrids in the Kalanchoë genus. BMC Plant Biology 15: 15. doi: 10.1186/s12870-014-0394-0
» https://doi.org/10.1186/s12870-014-0394-0 -
Lenth RV. 2022. emmeans: Estimated marginal means, aka least-squares means. R package version 1.7.2. https://CRAN.R-project.org/package=emmeans 3 Sep. 2022.
» https://CRAN.R-project.org/package=emmeans -
Lewis PO, Zaykin D. 2001. GDA (Genetic Data Analysis): Computer Program for the Analysis of Allelic Data (version 1.1). http://phylogeny.uconn.edu/software/ 16 Sep. 2021.
» http://phylogeny.uconn.edu/software/ -
Ley AC, Hardy OJ. 2017. Hybridization and asymmetric introgression after secondary contact in two tropical African climber species, Haumania danckelmaniana and Haumania liebrechtsiana (Marantaceae). International Journal of Plant Sciences 178: 421-430. doi: 10.1086/691628
» https://doi.org/10.1086/691628 -
Li YL, Liu JX. 2018. StructureSelector: A web-based software to select and visualize the optimal number of clusters using multiple methods. Molecular Ecology Resources 18: 176-177. doi: 10.1111/1755-0998.12719
» https://doi.org/10.1111/1755-0998.12719 - Lima SL. 2019. Filogeografia de Pterodon emarginatus e Pterodon pubescens (LEGUMINOSAE). PhD Thesis, Universidade Federal de Goiás, Brasil.
-
López-Caamal A, Tovar-Sánchez E. 2014. Genetic, morphological, and chemical patterns of plant hybridization. Revista Chilena de Historia Natural 87: 16. doi: 10.1186/s40693-014-0016-0
» https://doi.org/10.1186/s40693-014-0016-0 -
Lorenz-Lemke AP, Mäder G, Muschner VC et al 2006. Diversity and natural hybridization in a highly endemic species of Petunia (Solanaceae): A molecular and ecological analysis. Molecular Ecology 15: 4487-4497. doi: 10.1111/j.1365-294X.2006.03100.x
» https://doi.org/10.1111/j.1365-294X.2006.03100.x -
Marques I, Draper D, Riofrío L, Naranjo C. 2014. Multiple hybridization events, polyploidy and low postmating isolation entangle the evolution of neotropical species of Epidendrum (Orchidaceae). BMC Evolutionary Biology 14: 20. doi: 10.1186/1471-2148-14-20
» https://doi.org/10.1186/1471-2148-14-20 -
Melo PZ, Antunes AM, Fernandes JG et al 2022. Development of microsatellite markers in Pterodon pubescens and transferability to Pterodon emarginatus, two Brazilian plant species with medicinal potential. Rodriguésia 73: e01112021. doi: 10.1590/2175-7860202273091
» https://doi.org/10.1590/2175-7860202273091 -
Mors WB, Fo MFS, Monteiro HJ, Gilbert B, Pellegrino J. 1967. Chemoprophylactic agent in schistosomiasis: 14, 15-epoxygeranylgeraniol. Science 157: 950-951. doi: 10.1126/science.157.3791.950
» https://doi.org/10.1126/science.157.3791.950 -
Myers N, Mittermeier RA, Mittermeier CG, Fonseca GAB, Kent J. 2000. Biodiversity hotspots for conservation priorities. Nature 403: 853-858. doi: 10.1038/35002501
» https://doi.org/10.1038/35002501 -
Payseur BA, Rieseberg LH. 2016. A genomic perspective on hybridization and speciation. Molecular Ecology 25: 2337-2360. doi: 10.1111/mec.13557
» https://doi.org/10.1111/mec.13557 -
Pritchard JK, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945-959. doi: 10.1093/genetics/155.2.945
» https://doi.org/10.1093/genetics/155.2.945 -
Puechmaille SJ. 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. doi: 10.1111/1755-0998.12512
» https://doi.org/10.1111/1755-0998.12512 -
R Core Team. 2021. R: A language and environment for statistical computing. Vienna, Austria, R Foundation for Statistical Computing. https://www.R-project.org/ 3 Sep. 2022.
» https://www.R-project.org/ -
Ribeiro PL, Pereira ACS, Borba EL, Giulietti AM. 2018. Genetic and morphological diversity and evidence of hybridization in the “sempre-vivas” (Comanthera, Eriocaulaceae) endemic to the Chapada Diamantina, Bahia, Brazil. Flora 238: 60-71. doi: 10.1016/j.flora.2017.09.013
» https://doi.org/10.1016/j.flora.2017.09.013 -
Rieseberg LH, Ellstrand NC, Arnold M. 1993. What can molecular and morphological markers tell us about plant hybridization? Critical Reviews in Plant Sciences 12: 213-241. doi: 10.1080/07352689309701902
» https://doi.org/10.1080/07352689309701902 - Rieseberg LH, Carney SE. 1998. Plant hybridization. New Phytologist 140: 599-624.
- Rocha DMS. 2006. Aspectos taxonômicos, genéticos e reprodutivos de Pterodon pubescens Benth. e P. emarginatus Vog. (Leguminosae, Dipteryxeae). PhD Thesis, Universidade Estadual de Campinas, Brasil.
-
Schley RJ, Twyford AD, Pennington RT. 2022. Hybridization: A ‘double-edged sword’for Neotropical plant diversity. Botanical Journal of the Linnean Society 199: 331-356. doi: 10.1093/botlinnean/boab070
» https://doi.org/10.1093/botlinnean/boab070 -
Schnitzler CK, Turchetto C, Teixeira MC, Freitas LB. 2020. What could be the fate of secondary contact zones between closely related plant species? Genetics and Molecular Biology 43: e20190271. doi: 10.1590/1678-4685-GMB-2019-0271
» https://doi.org/10.1590/1678-4685-GMB-2019-0271 -
Servedio MR, Noor MAF. 2003. The Role of Reinforcement in Speciation: Theory and Data. Annual Review of Ecology, Evolution, and Systematics 34: 339-364. doi: 10.1146/annurev.ecolsys.34.011802.132412
» https://doi.org/10.1146/annurev.ecolsys.34.011802.132412 -
Sonsin-Oliveira J, Santiago VHD, Fagg CW, Oliveira RCD, Rocha DMS. 2022. Sucupira-branca (Pterodon-Fabaceae): Does wood anatomy support the distinction among species and hybrid? Acta Botanica Brasilica 36: e2021abb0124. https://doi.org/10.1590/0102-33062021abb0124
» https://doi.org/10.1590/0102-33062021abb0124 -
Taylor SA, Larson EL. 2019. Insights from genomes into the evolutionary importance and prevalence of hybridization in nature. Nature Ecology & Evolution 3: 170-177. doi: 10.1038/s41559-018-0777-y
» https://doi.org/10.1038/s41559-018-0777-y -
Taylor SA, Curry RL, White TA, Ferretti V, Lovette I. 2014. Spatiotemporally consistent genomic signatures of reproductive isolation in a moving hybrid zone. Evolution 68: 3066-3081. doi: 10.1111/evo.12510
» https://doi.org/10.1111/evo.12510 -
Teixeira MC, Turchetto C, Hartke C, Schnitzler CK, Freitas LB. 2019. Morphological and genetic perspectives of hybridization in two contact zones of closely related species of Petunia (Solanaceae) in southern Brazil. Acta Botanica Brasilica 33: 734-740. doi: 10.1590/0102-33062019abb0143
» https://doi.org/10.1590/0102-33062019abb0143 -
Turchetto C, Segatto ALA, Turchetto-Zolet AC. 2022. Biotic and abiotic factors in promoting the starting point of hybridization in the Neotropical flora: Implications for conservation in a changing world. Botanical Journal of the Linnean Society 200: 285-302. doi: 10.1093/botlinnean/boac042
» https://doi.org/10.1093/botlinnean/boac042 -
Vähä JP, Primmer CR. 2006. Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Molecular Ecology 15: 63-72. doi: 10.1111/j.1365-294X.2005.02773.x
» https://doi.org/10.1111/j.1365-294X.2005.02773.x -
Wagner HW, Sefc KM. 1999. IDENTITY 1.0 Centre for Applied Genetics. University of Agricultural Sciences, Vienna. http:/boku. ac. at/zag/forsgh/identity.htm 16 Sep. 2021.
» http:/boku. ac. at/zag/forsgh/identity.htm -
Wang LL, Zhang C, Tian B et al 2015. Reproductive isolation is mediated by pollen incompatibility in sympatric populations of two Arnebia species. Ecology and Evolution 5: 5838-5846. doi: 10.1002/ece3.1849
» https://doi.org/10.1002/ece3.1849 -
Yan LJ, Burgess KS, Milne R, Fu CN, Li DZ, Gao LM. 2017. Asymmetrical natural hybridization varies among hybrid swarms between two diploid Rhododendron species. Annals of Botany 120: 51-61. doi: 10.1093/aob/mcx039
» https://doi.org/10.1093/aob/mcx039 -
Zhang W, Dasmahapatra KK, Mallet J, Moreira GR, Kronforst MR. 2016. Genome-wide introgression among distantly related Heliconius butterfly species. Genome biology 17: 25. doi: 10.1186/s13059-016-0889-0
» https://doi.org/10.1186/s13059-016-0889-0
Supplementary Material
The following online material is available for this article:
Publication Dates
-
Publication in this collection
11 Aug 2023 -
Date of issue
2023
History
-
Received
14 Oct 2022 -
Accepted
30 June 2023
