- 1Laboratório de Genética & Biodiversidade, ICB, Universidade Federal de Goiás, Goiânia, Brazil
- 2Laboratório de Ecologia Espacial e Conservação (LEEC), Departamento de Biodiversidade, Universidade Estadual Paulista Júlio de Mesquita Filho, Rio Claro, Brazil
- 3Escola de Agronomia, Universidade Federal de Goiás, Goiânia, Brazil
Changes in landscape structure can affect essential population ecological features, such as dispersal and recruitment, and thus genetic processes. Here, we analyze the effects of landscape metrics on adaptive quantitative traits variation, evolutionary potential, and on neutral genetic diversity in populations of the Neotropical savanna tree Tabebuia aurea. Using a multi-scale approach, we sampled five landscapes with two sites of savanna in each. To obtain neutral genetic variation, we genotyped 60 adult individuals from each site using 10 microsatellite loci. We measured seed size and mass. Seeds were grown in nursery in completely randomized experimental design and 17 traits were measured in seedlings to obtain the average, additive genetic variance (Va) and coefficient of variation (CVa%), which measures evolvability, for each trait. We found that habitat loss increased genetic diversity (He) and allelic richness (AR), and decreased genetic differentiation among populations (FST), most likely due to longer dispersal distance of pollen in landscapes with lower density of flowering individuals. Habitat amount positively influenced seed size. Seeds of T. aurea are wind-dispersed and larger seeds may be dispersed to short distance, increasing genetic differentiation and decreasing genetic diversity and allelic richness. Evolvability (CVa%) in root length decreased with habitat amount. Savanna trees have higher root than shoot growth rate in the initial stages, allowing seedlings to obtain water from water tables. Landscapes with lower habitat amount may be more stressful for plant species, due to the lower plant density, edge effects and the negative impacts of agroecosystems. In these landscapes, larger roots may provide higher ability to obtain water, increasing survival and avoiding dying back because of fire. Despite the very recent agriculture expansion in Central Brazil, landscape changes are affecting neutral and adaptive variation in T. aurea. Several populations have low additive genetic variation for some traits and thus, may have limited evolvability, which may jeopardize species long-term persistence. The effect of habitat loss on highly variable neutral loci may only be detected after a certain threshold of population size is attained, that could become dangerously small masking important losses of heterozygosity endangering species conservation.
Introduction
Agricultural landscapes are now occupying most of Neotropical savannas. The Brazilian Cerrado biome is the largest Neotropical savanna and one of the world’s biodiversity hotspots because of its high level of endemism and threatening (Myers et al., 2000). Over the last 50 years, more than 50% of its vegetation cover has been cleared or transformed into agriculture, pasture or urban area (Sano et al., 2010). Landscapes in the Cerrado biome are now comprised by islands of savanna surrounded by ocean of crops and pastures (see Latrubesse et al., 2019), thus jeopardizing species long-term conservation and ecosystem services maintenance.
Changes in landscape structure, such as habitat amount, edge effect and heterogeneity, can affect essential ecological processes such as dispersal (Corlett, 2017; Bovo et al., 2018), seed predation (Mendes et al., 2016), and pollination (Hadley and Betts, 2011; Duarte et al., 2018), and consequently compromise population persistence (Santos et al., 2016; Regolin et al., 2017). Landscape changes may also affect population connectivity and genetic diversity in different taxa (e.g., Dixo et al., 2009; Carvalho et al., 2015; Jackson and Fahrig, 2016; Moraes et al., 2018; González-Fernández et al., 2019). This occurs mainly due to reduction in population effective size and isolation (Gilpin and Soulé, 1986; Frankham et al., 1999). In this way, understanding how the still-ongoing landscape change affects savanna species may help designing sound conservation and management strategies.
In agricultural landscapes, anthropogenic matrix characteristics may be as important as habitat size and amount due to the potential restriction the agroecosystems may impose to species movement, dispersal and gene flow (Holderegger et al., 2010; Driscoll et al., 2013; Silva et al., 2015; Barros et al., 2019). Matrix type may affect functional connectivity (i.e., the degree to which the matrix facilitates or impedes species movement among habitat patches) because it may be filters or physical barriers to species’ movement (Ricketts, 2001; Shamoon et al., 2018). In plants, landscape structure may affect connectivity due to the influence on pollen dispersal, seed dispersal distance, germination and establishment (Damschen et al., 2014; Trakhtenbrot et al., 2014; Auffret et al., 2017). For instance, population size and productivity may decrease seed germinability, thus decreasing final seed dispersal distance due to habitat quality (Soons and Heil, 2002; Soons et al., 2004, 2005). Thus, although wind-dispersed species may be favored to some degree by habitat fragmentation – due to opened matrix increasing long-distance dispersal (Tackenberg et al., 2003) – changes in landscape composition and configuration may disrupt recruitment (Soons and Heil, 2002).
Most studies addressing the impacts of landscape composition and configuration in plant genetic diversity focus on genetic variation at neutral loci (e.g., Schmidt et al., 2009; Carvalho et al., 2015). Adaptive genetic variation can determine population persistence in the face of environmental changes (Levins, 1968; Meyers and Bull, 2002; Holderegger et al., 2006), such as the rapid landscape changes in the Anthropocene. Genetic variability at adaptive loci depends on the strength of selection and neutral molecular markers might be poor indicators of selection pressures at such loci (McKay and Latta, 2002). The loss of genetic variability at adaptive loci may lead to the loss of individual fitness and thus in population evolutionary potential (Lande, 1988; Reed and Frankham, 2002; Bijlsma and Loeschcke, 2012). Measurement of adaptive genetic variation often requires controlled experiments due to quantitative variation, that are not attainable for many species and thus is often overlooked in plant landscape genetics studies, hindering our understanding of the effects of land use in functional traits and adaption (Holderegger et al., 2006, 2010).
Here, we analyze the effects of landscape composition and configuration in variation at adaptive quantitative traits and at neutral genetic variability in populations of the Neotropical savanna tree Tabebuia aurea (Bignoniaceae); for this we used a multi scale approach. This species is broadly distributed across seasonal and wet savannas of Brazil, Bolivia and Paraguay, but locally it is distributed in well-delimited patches, usually with high density. It is hermaphroditic with mixed-mating reproductive system (Braga and Collevatti, 2011), pollinated by large-sized bees, such as carpenter bees and bumblebees, and seeds are wind-dispersed. Wind may promote long-distance seed dispersal in savanna trees with winged seeds, because long-distance seed dispersal by wind may be associated with plant height and seed morphology (Tackenberg et al., 2003). Because landscape structure may affect dispersal distance, seed germination and establishment we hypothesize that the relationships between our response variables (genetic variability, inbreeding and genetic differentiation) versus explanatory variables (habitat amount, compositional heterogeneity, matrix quality and connectivity) will have a vertex, i.e., a turning point. Away from the turning point, in landscapes with very high or very low habitat amount, compositional heterogeneity, matrix quality and connectivity we expect that genetic variability at neutral and adaptive loci will be lower and genetic differentiation and inbreeding will be higher. Conversely, in landscapes with average habitat amount, compositional heterogeneity, matrix quality and less isolated patches genetic variability at neutral and adaptive loci will be higher and genetic differentiation and inbreeding will be lower. We also hypothesize that population evolvability will be reduced due to habitat and connectivity losses following the trend expected for genetic variability. Because genetic variability is affected by effective populations size and changes in adaptive traits may be a balance between drift and selection (Wright, 1931; Kimura, 1983), we also analyze the effect of effective population size in adaptive and neutral genetic variation.
Materials and Methods
Study Sites and Landscape Sampling Design
The study was carried out in five landscapes of savanna in the Cerrado biome in Goiás state, Central-West Brazil (Figure 1 and Appendix Table S1), one of the most important Brazilian agribusiness regions. Each landscape was defined by a buffer of 6 km radius, where two sites of savanna with populations of T. aurea were sampled (total of 10 sites). The sites were either structurally isolated in different savanna fragments (P1, P2, and P4, Figure 1 and Appendix Table S1) for many landscapes, or in the same large savanna protected area (P3 and P5, Figure 1 and Appendix Table S1). In any case, we selected sampling sites spaced apart by at least 1.4 km (Appendix Table S1).
Figure 1. Geographic distribution of the five landscapes and the 10 sampling sites of Tabebuia aurea in the Brazilian Cerrado. (A) The distribution of the Cerrado biome in Brazil and the landscapes sampled. (B) The landscapes represented by a buffer of 6 km and the sampling sites (black dots). Land use categories are in legends. For details of sampling sites see Appendix Table S1.
We mapped land cover in these five landscapes using high-resolution images available at the map database of the Geographic Information System in ArcGis v.9.3 environment (Esri®). Mapping was performed using visual digitalization and manual classification at a scale of 1:5,000 m, and we made exhaustive field inspection to assure high map accuracy. Our final map comprised 11 different land cover classes: (i) water courses; (ii) savanna; (iii) riparian forest; (iv) seasonal forest; (v) wetland; (vi) pasture; (vii) agriculture (corn or soybean); (viii) rural building; (ix) urban area; (x) road and train rail, and (xi) Eucalyptus spp. plantation (Figure 1).
Genetic Data
Genetic Analysis for Neutral Loci
To obtain genetic parameters for neutral loci, we sampled expanded leaves of ∼60 individuals of T. aurea in each site (Appendix Table S2). In sites with less than 60 we sampled all individuals. Total genomic DNA extraction followed CTAB procedure and all individuals were genotyped using 10 microsatellite loci (Braga et al., 2007) following PCR protocol described in Collevatti et al. (2014). DNA fragments were sized in GS 3500 Genetic Analyzer (Applied Biosystems, CA) using GeneScan ROX 500 size standard (Applied Biosystems, CA), and were scored using GeneMaper v5.0 software (Applied Biosystems, CA). Genotyping errors (allele dropout and null allele) were analyzed using Micro-Checker 2.2.3 software (van Oosterhout et al., 2004).
For each site (10), we estimated genetic diversity (expected heterozygosity under Hardy-Weinberg equilibrium, He; Nei, 1978), allelic richness based on rarefaction analysis (AR; Mousadik and Petit, 1996) and inbreeding coefficient (f; Wright, 1951). All analyses and randomization-based tests for deviation from Hardy-Weinberg equilibrium were performed with the software FSTAT 2.9.3.2 (Goudet, 2002). We estimated genetic differentiation among all pairs of sites nested within landscapes (5) using Wright’s FST, and fixation index (inbreeding) FIS obtained from an analysis of variance of allele frequencies (Weir and Cockerham, 1984). We also estimated GST’ (Hedrick, 2005) that is based on FST, but takes into account observed diversity within population and number of subpopulations, and Jost’ D (Jost, 2008) that is based on effective number of alleles instead of expected heterozygosity. To estimate the contribution of stepwise mutation model to genetic differentiation, we estimated Slatkin’s RST (Slatkin, 1995) and tested the hypothesis that FST = RST using the software Spagedi (Hardy and Vekemans, 2002).
We used microsatellite markers to estimate effective population size (Ne). We used the molecular co-ancestry method (Nomura, 2008) implemented in NeEstimator V2.1 (Do et al., 2014) and tested the hypothesis that genetic variability and differentiation depend on Ne (see below).
Adaptive Quantitative Trait Variation
We obtained quantitative genetic variation for traits related to the fitness of seeds and seedlings, i.e., adaptions related to the initial establishment of seedlings. We focused on seeds and seedlings because these stages are related to plant dispersal and establishment, and seedling is the stage at which plants are most vulnerable to the effects of environment variables (Tomlinson et al., 2012).
We measured 20 quantitative traits related to seed size (longitudinal and transversal diameter) and mass and to seedling growth (height and diameter growth), size (height, leaf size and number, shoot and root length) and mass (root and aboveground). Seed size and mass are strongly related to seedling successful establishment and also with other adaptive traits such as plant height and specific leaf area (Westoby et al., 1996). Height and leaf traits are frequently used to assess the variation in plant ecological strategies because they represent plant response to stress and resource use (Westoby, 1998; Reich, 2000; Westoby et al., 2002). Plant height is related to competitive capacity or disturbance avoidance (Westoby, 1998), which is important in savannas where fire events are frequent and taller plants may avoid top kill due to the death of canopy meristem (Hoffmann et al., 2012). Leaf area represents a functional strategy related to resource retention and photosynthesis rate (Reich et al., 1992; Westoby, 1998; Ackerly et al., 2002; Perez-Harguindeguy et al., 2013) and dry mass is related to the regulation of water loss through evapotranspiration (Poorter et al., 2009), increasing water use efficiency in low moisture soils (Ackerly, 2004; Poorter et al., 2009).
To obtain quantitative genetic variation, we sampled seeds from randomly chosen trees in each site, obtaining a total of 1070 seeds (Appendix Table S2). The number of trees sampled differed due to the number of adults siring fruits. Seeds were measured and weighted to obtain seed traits, SLD (seed longitudinal diameter, mm), STD (seed transversal diameter, mm) and SM (seed mass, mg) (Appendix Table S3). Then, seeds were grown in nursery in a completely randomized experimental design. The number of seedlings analyzed per tree per site differed due to variation in germination (Appendix Table S2). We also obtained the number of days to shoot (time to seed germination, TG) and the proportion of seeds that germinated (PG).
Plant aboveground height and stem diameter were measured 76, 116, 133, and 145 days after seed germination, and growth rate was obtained from the regression coefficient (β) of height (HGR, cm/day) and diameter (DGR, mm/day). We also registered the initial and final seedling height (IH and FH, cm) and diameter (ID and FD, mm). We count the total number of leaves (NL), and leaf length (LL, mm) and width (LW, mm) were obtained from the mean value among three leaves per seedling, measured up to 145 days after germination (Appendix Table S3). After 145 days, seedlings were taken from the nursery pots to measure (Appendix Table S3) aboveground shoot length (ASL, cm), root length (RL, cm). We also obtained aboveground green mass (AGM, g) and dry mass (ADM, g), and root green (RGM, g) and dry masses (RDM, g).
For statistical analyses, we first performed Pearson correlation analyses among quantitative variables for sites and landscapes (Appendix Tables S4, S5), to minimize the correlation in our set of quantitative traits. Therefore, for the subsequent analyses, we excluded the quantitative traits with correlation coefficient r > 0.5. We kept for seed traits SLD, STD and SM, and for seedlings, NL, LW, RL, RDM, and ASL.
We estimated the variance (Va) and the additive genetic coefficient of variation (CVa%), for each selected trait in each site, as a surrogate of population evolutionary potential or evolvability (Houle, 1992; Hansen et al., 2011). We also obtained narrow-sense heritability (h2) of each trait in each site. These parameters were estimated using restricted maximum likelihood (REML) analysis implemented in the software SELEGEN–REML/BLUP (Resende, 2016). We used model 82, implemented for open-pollinated sib families and mixed mating system. In this model, additive genetic variance is corrected by the coefficient . For this, we estimated θ, the coancestry coefficient within family (Vencovsky et al., 2001; Tambarussi et al., 2018), , and the selfing rate (Vencovsky and Crossa, 2003), f is the inbreeding coefficient within population. Inbreeding was obtained from neutral loci (see above). To estimate additive genetic variance and the other parameters (CVa% and h2) we used , where is the genetic variance among families. Because seed traits were measured in seeds sampled in the field, we cannot estimate additive genetic variation. Thus, these parameters were estimated only for seedling traits obtained from controlled experiment in nursery.
To obtain the additive quantitative differentiation between pair of sites nested within landscapes, we estimated QST (Spitze, 1993; Prout and Barker, 1993) and the analogous PST (Leinonen et al., 2006) for each selected trait. QST estimates quantitative genetic differentiation based on population additive genetic variation (Spitze, 1993), and was estimated for seedling traits measured in nursery under experimental conditions. PST is analogous to QST (Leinonen et al., 2006), but estimates phenotypic differentiation when genetic and environmental effects cannot be detached, due to uncontrolled environmental conditions. We estimated PST for seed traits because they were sampled in the field, under uncontrolled environmental conditions. QST and PST parameters were estimated using model 5 implemented in SELEGEN–REML/BLUP software (Resende, 2016), also correcting the additive genetic variance for mixed mating system (see above).
Landscape Metrics
We quantified landscape variables at node and link levels (Figure 2), following Wagner and Fortin (2013) framework. To calculate landscape metrics at the node level, we generated buffers of 0.5 km around each focal site (Figure 2A). We chose this minimum distance to avoid overlap between buffers. For the sites in protected areas (landscapes P3 and P5, Figure 1), we considered the same values for landscape metrics (Figure 2B).
Figure 2. Sampling design to quantify landscape metrics at node and link levels for Tabebuia aurea in the Brazilian Cerrado. (A) Sampling design at node level in landscapes with sites in different patches of savanna. (B) Sampling design at node level in landscapes within the same patch of savanna, which corresponded to protected areas. (C) Sampling design at link level showing the buffers with different radii around the midpoint between the two sampling sites. For details of sampling sites see Figure 1, and Appendix Table S1.
To perform the analysis at the link level we identified the midpoint between the sampling sites, and performed multi-scale analysis, generating different buffers of 2, 4, and 6 km radius around each midpoint (Figure 2C). We chose the minimum buffer size because of pollen dispersal distance (at least 2 km, Braga and Collevatti, 2011). Additionally, buffer of 6 km of radius corresponded to the maximum size to avoid overlap at the link level (Figure 1).
To characterize landscape structure we quantified different landscape metrics related to the composition and configuration of landscape elements. As savanna is considered an essential vegetation type to the establishment of T. aurea, we estimated the total habitat amount calculating the percentage of savanna at both node and link levels at different spatial scales (Appendix Tables S6, S7). To avoid correlation between habitat amount at different spatial scales at the link level, we subtracted the percentage of savanna in the additional scales and considered these new values as an additional metric. For instance, the metric “% savanna 4_2km” corresponded to the percentage of savanna in the scale of 4 km minus the percentage of savanna in the scale of 2 km (Appendix Table S7). To account for the effect of pollinators’ habitat, we summed up the amount of open vegetation where T. aurea can eventually establish, besides savanna (pasture, wetlands) and pollinators’ habitat (seasonal and riparian forests) to generate the variables: (i) % savanna + pasture + wetlands and (ii) savanna + pasture + seasonal forest + riparian forest + wetlands. We also avoided correlation subtracting the values of each metric in additional scales (Appendix Table S7).
The compositional heterogeneity was calculated using the Shannon index (SHDI) available in Fragstats (McGarigal et al., 2012), which included all land cover types within the buffers (Appendix Tables S6, S7). To estimate the functional landscape connectivity we selected all patches of savanna at different spatial scales and made individual buffers of 500 m of radius around each patch. We chose this radius due to foraging distance of potential pollinators of T. aurea such as Bombus spp. and Centris spp. (Braga and Collevatti, 2011). We considered as functional connectivity the sum of the area in hectares of all patches of savanna plus the area connected by the buffers (Appendix Tables S6, S7). The functional landscape connectivity was calculated using GRASS GIS 7.5 (GRASS Development Team, 2018).
To address matrix quality we used the Buffer Matrix Quality (BMQ) metric (Lion et al., 2014; see also Silva et al., 2015), BMQ = ∑PiQi/∑Pi, where Pi is the percentage of each land cover class in each landscape at different spatial scales, and Qi is the quality score of each land cover class considering the requirements and resources that each land cover can provide to focal species. We scored each land cover based in the characteristics of T. aurea’s pollinators, since T. aurea is wind-dispersed and seedlings can eventually establish in wetlands and grasslands. For this, we considered our personal experience, the literature and the opinion of other experts in bees’ behavior. Our scores values ranged from the 0.2 to 1.0, where low values corresponded to areas more impermeable to species (Appendix Tables S6, S7).
For landscape variables at the link level, we first assessed the scale of effect calculating Pearson correlation between the same landscape metrics at different spatial scales. Because all correlations were high (i.e., r > 0.90), we kept only metrics for 2 km spatial scale. We then assessed collinearity between landscape metrics within 2 km spatial scale (Appendix Table S8). At the node level we calculated correlation coefficients between all landscape metrics at 0.5 km spatial scale (Appendix Table S9). We excluded from the subsequent analyses the landscape metrics with correlation coefficient r > 0.5. At node level, we kept habitat amount (%) and compositional heterogeneity (Appendix Table S9). At link level, we used habitat amount at 2 km spatial scale (Appendix Table S8).
Data Analysis
We analyzed the effects of landscape metrics using Generalized Linear Model (GLM) for both neutral genetic parameters and adaptive quantitative traits (Table 1). Because genetic variability may be affected by effective population size (Ne), we modeled genetic parameters using both landscape metrics and Ne (Table 1).
Table 1. Models performed at node and link levels for both neutral and adaptive quantitative traits measured in seeds and seedlings of Tabebuia aurea from 10 sites and 5 landscapes in the Brazilian Cerrado.
We evaluated the distribution of the residuals of all models and used the Gaussian distribution that fitted better to our data. To find the best predictive model at node level we considered both the significance and the Akaike Information Criteria (AIC). We estimated AIC corrected for small sample sizes (AICc) and the difference of each model and the best model (ΔAICci). We also estimated Akaike’s Weight of Evidence (wAICc) as the relative contribution of each model to explain the observed pattern, given a set of competing models (Burnhan and Anderson, 2002). Models with ΔAICc < 2 were considered as equally plausible to explain the observed pattern (Zuur et al., 2009).
For link level, we used the significance to select the best predictive model, because of the small sample size (five landscapes). Additionally, the sampling sites in P5 (F9 and F10, Figure 1 and Appendix Table S1) had outlier behavior in the analysis performed at the link level for neutral loci, hindering model fitting. To minimize this effect, we log transformed (log10) predictor variables. All statistical analyses were performed with R version 3.6.1. (R Core Team, 2019) using the following packages bblme (Bolkerand et al., 2017), visreg (Breheny and Burchett, 2017), and ggplot2 (Wickham, 2016).
Results
Neutral Genetic Diversity and Differentiation
All loci had high genetic diversity and allelic richness (Appendix Table S10). Genetic diversity (He) was high in all populations ranging from 0.874 to 0.921 (Appendix Table S10). Allelic richness (AR) was also high, ranging from 14.2 to 18.5 (Appendix Table S10). Inbreeding coefficient (f) was significant for all populations, ranging from 0.121 to 0.231. FST across landscapes ranged from 0.008 to 0.240, GST’ ranged from 0.115 to 0.304, and Jost’ D ranged from 0.104 to 0.277 (Appendix Table S11). Slatkin’s RST (RST = 0.077, SE = 0.048, p < 0.001) was not significantly different from FST (p = 0.852). Inbreeding coefficients within landscapes (FIS) was also high (Appendix Table S11). All populations showed low effective population sizes (Appendix Table S10), ranging from 19.7 to 56.0.
Adaptive Quantitative Trait Variation
Overall, populations had high additive genetic variance (Va) in quantitative traits, but some had low variance in several traits (e.g., P2F4, Appendix Table S12), leading to low values of additive genetic coefficient of variation and narrow-sense heritability (Appendix Table S12). For instance, Va ranged from 0.0479 to 25.9293 for leaf width, and 0.0025 to 0.6116 for root dry mass (Appendix Table S12). Narrow-sense heritability for seedling traits showed a wide range among sites (Appendix Table S12), ranging from 0.0045 to 0.9419 for leaf width, 0.0024 to 0.6877 and for root dry mass, for example. Genetic differentiation in quantitative traits (QST) was low in most of the landscapes (Appendix Table S13), ranging from 0.005 (RDM) to 0.8140 (RL). Phenotypic differentiation in seed traits (PST) was also low for most pairwise comparisons, ranging from 0.0003 to 0.0835 (Appendix Table S13).
Landscape Effect on Neutral Genetic Diversity and Differentiation
Habitat amount was the model that best explained the variation in neutral genetic diversity (He, wAICc = 0.94; Table 2) and allelic richness (AR, wAICc = 0.87; Table 2). Genetic diversity (Figure 3A) and allelic richness (Figure 3B) decreased with habitat amount. None of the models tested was significant for inbreeding (f, all p > 0.10).
Table 2. Model selection of the competing hypotheses to explain the patterns of variation in neutral genetic variability and adaptive quantitative traits, in populations of Tabebuia aurea in landscapes of the Brazilian Cerrado.
Figure 3. Relationships of neutral genetic variability and habitat amount in five landscapes and 10 sampling sites of Tabebuia aurea in the Brazilian Cerrado. (A) Genetic diversity (He). (B) Allelic richness (AR). Black line is the linear regression fit and shaded area is the 95% confidence interval. All fits were significant (p < 0.05).
At link level, genetic differentiation based on FST was significantly related to habitat amount (β = 0.021; p = 0.031). Populations within landscapes with higher amount of habitat had higher genetic differentiation (Figure 4A). G’ST and Jost’s D were not explained by our model (p > 0.10). Inbreeding coefficient (FIS) was also not significantly related to habitat amount.
Figure 4. Relationship of (A) genetic differentiation at neutral loci (FST) and at (B) adaptive quantitative trait (QST) root length (RL), and habitat amount at 2 km spatial scale in five landscapes of Tabebuia aurea in the Brazilian Cerrado. Black line is the linear regression fit and shaded area is the 95% confidence interval. Fits were significant (p < 0.05).
Landscape Effect on Adaptive Quantitative Traits
Habitat amount explained variation in seed size (SLD) between sites within landscapes, although marginally significant (wAICc = 0.65, p = 0.099; Table 2). Sites in landscapes with higher amount of habitat tended to have larger seeds (Figure 5A). Number of leaves (NL) was better explained by effective population size (Ne, wAICc = 0.74) than by landscape structure, but was marginally significant (p = 0.077; Table 2). Number of leaves was negatively related to Ne (Figure 5B). Landscape structure and Ne did not explain variation in mean STD, SM, RL, RDM, and ASL between sites within landscapes (all p > 0.10).
Figure 5. Relationships of adaptive quantitative traits and habitat amount and effective population size (Ne) in five landscapes of Tabebuia aurea in the Brazilian Cerrado. (A) Mean seed longitudinal diameter (SLD) and habitat amount. (B) Mean number of leaves (NL) and effective population size (Ne). Black line is the linear regression fit and shaded area the 95% confidence interval. All fits were marginally significant (p < 0.10).
Landscape structure explained evolvability of seedling traits (Table 1), measured by the additive genetic coefficient of variation (CVa%). NL (wAICc = 0.70, p = 0.084) and RL (wAICc = 0.48, p = 0.064) coefficients of variation were better explained by habitat amount (Table 2). Coefficient of variation in NL increased with habitat amount (Figure 6D), while RL coefficient of variation decreased with habitat amount (Figure 6A). RL (wAICc = 0.37, p = 0.082) and RDM (wAICc = 0.83, p = 0.062) coefficients of variation were explained by compositional heterogeneity, although marginally significant (Table 2), and both tended to increase with compositional heterogeneity (Figures 6B,C).
Figure 6. Relationships of additive genetic coefficient of variation (CVa%) of adaptive quantitative traits and landscape features in five landscapes and 10 sampling sites of Tabebuia aurea in the Brazilian Cerrado. (A) Additive genetic coefficient of variation (CVa%) in number of leaves (NL) and habitat amount. (B) Additive genetic coefficient of variation (CVa%) in root length (RL) and habitat amount. (C) Additive genetic coefficient of variation (CVa%) in root length (RL) and compositional heterogeneity. (D) Additive genetic coefficient of variation (CVa%) in dry root mass (RDM) and compositional heterogeneity. Black line is the linear regression fit and shaded area is the 95% confidence interval. Fits were significant (p < 0.05) or marginally significant (p < 0.10).
At link level, genetic differentiation at root length (RL) based on QST was significantly related to habitat amount (b = 1.041; p = 0.039). Populations within landscapes with higher amount of habitat had higher quantitative genetic differentiation in RL (Figure 4B). Our models did not explain (p > 0.10) differentiation at other seedling or seed adaptive traits.
Discussion
Pollen Dispersal May Explain the Relationships of Habitat Amount and Neutral Genetic Diversity and Differentiation
Our findings show that contemporary changes in landscape affect genetic diversity and differentiation at neutral microsatellite loci among populations of T. aurea from the Cerrado biome. However, unexpectedly, habitat loss increased genetic diversity and allelic richness and decreased genetic differentiation among populations. This result is most likely due to pollinator behavior and pollen dispersal patterns. Tabebuia aurea has a mixed mating system, with potential long-distance pollen dispersal of at least 2.6 km (Braga and Collevatti, 2011). Pollinators of this species, mainly bumblebees and carpenter bees, have long distance flight capacity, and may carry pollen over long distances (Goulson and Stout, 2001). However, in mass-flowering species, such as T. aurea, the residence time of pollinators may be greater decreasing pollen carry-over and thus effective population size (Frankie et al., 1976). Therefore, the mass-flowering pattern added by the clumped distribution of individuals within populations may affect pollinator foraging behavior, increasing pollinator residence time and leading to a high proportion of low-distance pollen dispersal (Braga and Collevatti, 2011). Habitat loss may decrease population size and density in T. aurea, leading to lower number of individuals flowering within populations, which may decrease the frequency of pollination between nearest neighbors and increase pollen dispersal distance (Kitamoto et al., 2006). In fact, populations in protected areas, such as the landscapes P3 and P5, with the highest habitat amount and structural connectivity, showed the lowest values of genetic diversity and allelic richness, but the highest genetic differentiation for the three parameters (FST, G’ST or Jost’s D).
Our results contrast to other studies with T. aurea (Collevatti et al., 2014) and other savanna tree species (e.g., Caryocar brasiliense, Collevatti et al., 2001) that found no correlation between genetic diversity and habitat loss. This is most likely due to the spatial scale, since the former studies were performed in biome scale, where phylogeographic and biogeographic factors may be more important than local process of landscape changes.
Our findings also evinced high inbreeding within populations and high fixation index (FIS) within landscapes, but not related to landscape structure. In fact, T. aurea have high proportions of self-pollination and biparental inbreeding (Braga and Collevatti, 2011) and high inbreeding in several populations throughout its geographical distribution (Collevatti et al., 2014), which may explain the high values found in this work. Gene dispersal based on seedling fine scale genetic structure is low (nearly 300 m, Braga and Collevatti, 2011), which may also explain the lower differentiation in populations with low habitat amount. Although the method based on fine scale genetic structure may underestimate seed dispersal distance, long-distance dispersal is mainly due to pollen dispersal in many Neotropical plants (Hamilton and Miller, 2002; Ballesteros-Mejia et al., 2016). Although wind may potentially promote long-distance seed dispersal in opened areas (Tackenberg et al., 2003), landscape features may affect seed dispersal distance, seed germination and seedling recruitment (Damschen et al., 2014; Trakhtenbrot et al., 2014). In fact, sites with higher habitat amount had larger seeds (SLD, see below). Although larger seeds may increase seedling survival and thus recruitment, they may disperse to shorter distances (Ganeshaiah and Shaanker, 1991; Greene and Johnson, 1993), increasing genetic relationship within habitat fragments, and thus genetic differentiation, and decreasing genetic diversity and allelic richness.
The Trade-Off Between Seedling Survival and Seed Dispersal May Explain the Relationships Between Adaptive Quantitative Traits and Landscape Features
Our findings show that landscape structure explained the adaptive quantitative variation of traits related to seed size and evolvability of several traits. Populations in sites with higher habitat amount tended to have larger seeds (SLD, seed longitudinal diameter). Seeds of T. aurea are wind dispersed and large seeds may be dispersed to short distance. There is a trade-off between seed size improving germination and seedling survival in savannas (e.g., Lahoreau et al., 2006) and dispersal (Ganeshaiah and Shaanker, 1991; Greene and Johnson, 1993). In fact, SLD is positively correlated with percentage of germination (PG, r = 0.767), time to germinate (TG, r = 0.397), plant height (IH, r = 0.327), and plant growth rate (HGR, r = 0.439). Thus, although seedlings from large seeds may growth faster, they may be dispersed to shorter distances. In the same manner, sites with higher habitat amount may provide more resources for larger seeds (Westoby et al., 1996) that may disperse to shorter distances.
On the other hand, root dry mass (RDM) and root length (RL) evolvability increased with compositional heterogeneity and RL decreased with habitat amount. In addition, RL quantitative differentiation (QST) was also higher in landscapes with higher habitat amount. In fact, root mass and length are important adaptive traits in savanna trees, which have higher root growth rates than shoot growth rates in the initial stages, allowing seedlings to obtain water from water tables (see Goldstein et al., 2008 for a review). Landscapes with higher compositional heterogeneity in agroecosystem-dominated regions may be more stressful for plant species, due to the lower habitat amount, edge effects and the presence of urban areas, roads and agriculture. In these areas, larger roots may provide higher ability to obtain water and nutrients, increasing survival, and avoid dying back because of fire (Hoffmann and Franco, 2003; Hoffmann et al., 2004). Fire is frequent in Cerrado during the dry season, and the frequency of fire is higher in agroecosystem landscapes (Pivello, 2011) that have higher spatial and temporal heterogeneity. In addition, in agroecosystems variation in soil nutrient in heterogeneous landscapes may also contribute to high variation in root size, because of selection to fast growth of roots in nutrient poor and dry areas, and relaxed selection in areas with higher nutrient inputs (for instance in edges of savanna and agriculture areas that are periodically fertilized) and moisture. The positive relationship of number of leaves (NL) additive genetic coefficient of variation and habitat amount is most likely due to the stressful environment in landscapes with lower habitat amount, similarly to RL, but in opposite direction due to negative relationship between root growth rates than shoot growth rates in the initial stages. Traits related to leaf area are correlated to plant water use efficiency, photosynthesis rate and resources retention (Westoby et al., 2002). Landscapes with lower habitat amount may be more instable leading to loss of variation in leaf traits and thus lower evolutionary potential.
Conclusion
We found that landscape changes are affecting genetic diversity and adaptive variation in the savanna tree species T. aurea, despite the very recent agriculture expansion in Central Brazil (∼60 years). On one hand, the loss of habitat increases neutral genetic diversity and decreases genetic differentiation, most likely due to long distance pollen dispersal. The lower genetic differentiation between sites with lower habitat amount may also indicate that agroecosystem matrix is somehow permeable to pollinators. On the other hand, habitat amount decreases adaptive quantitative traits and increases evolvability. It is important to note that, as habitat fragmentation process (sensu Fahrig, 2003) is a recent event when considering the life cycle of the species, the effect of habitat loss on genetic diversity at highly variable neutral loci, such as T. aurea’s microsatellite loci (see Collevatti et al., 2014), may only be detected after a certain threshold of population size is attained due to Allee effect (Berec et al., 2007). As suggested for other savanna species (C. brasiliense; Collevatti et al., 2001), such a threshold could become dangerously small, because the high allelic variation at microsatellite loci could mask important losses of heterozygosity or very drastic events. Indeed, habitat loss had already affected some important quantitative adaptive traits, which evolve faster than neutral loci (McKay and Latta, 2002) and though may respond faster to landscape changes.
Several adaptive quantitative traits show evolutionary potential, meaning that populations still have enough variation to respond to selection, and thus to environmental changes. However, it is important to note that several populations have already low additive genetic variation for some traits at seedling stage. Evolutionary potential at short-term depends on population additive genetic variance (Hansen et al., 2011) and thus, those populations may have no evolvability for some traits, which may jeopardize species long-term persistence in these landscapes. Although mortality is higher at initial stage in savanna tree species, the low additive genetic variation for some traits may be also due to the time length of our study, because variation may be detected in other life-stages, such as juveniles. Despite the limitations of the experiment, which analyzed only seedlings and not other later stages, this result highlights the importance of the maintenance of functional connectivity among sites and landscapes, allowing the migration of advantageous alleles and thus the recovery of evolutionary potential.
Data Availability Statement
The data and additional supporting information may be found in the online version of this article as supporting information.
Author Contributions
RC and MR conceived and funded the work. FR and TA obtained the data. RC, JS, MR, and LC designed the experiment, field sampling, and statistical analyses. JS carried out the analyses. RC, JS, and MR wrote the original draft. All authors contributed to manuscript and approved the final version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by grants to the research network PELD COFA supported by MCT/CNPq/CAPES/(Project No. 441278/2016-7), FAPEG (Project No. 201710267000331) and CAPES/PROCAD (Project No. 88881.068425/2014-01). MR has been funded by FAPESP (Project No. 2013/50421-2). JS received a CAPES/PROCAD postdoctoral fellowship and a FAPESP postdoctoral grant (Project No. 2019/09713-6). FR received a fellowship from FAPEG and TSA from CAPES. We thank Carolina Novaes for helping with fieldwork. RC, MR, and LC have continuously been supported by productivity grants from CNPq, which we gratefully acknowledge.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00259/full#supplementary-material
References
Ackerly, D., Knight, C., Weiss, S., Barton, K., and Starmer, K. (2002). Leaf size, specific leaf area and microhabitat distribution of chaparral woody plants: contrasting patterns in species level and community level analyses. Oecologia 130, 449–457. doi: 10.1007/s004420100805
Ackerly, D. (2004). Functional strategies of chaparral shrubs in relation to seasonal water deficit and disturbance. Ecol. Monogr. 74, 25–44. doi: 10.1890/03-4022
Auffret, A. G., Rico, Y., Bullock, J. M., Hooftman, D. A. P., Pakeman, R. J., Soons, M. B., Suarez-Esteban, A., Traveset, A., Wagner, H. H., and Cousins, S. A. O. (2017). Plant functional connectivity – integrating landscape structure and effective dispersal. J. Ecol. 105, 1648–1656. doi: 10.1111/1365-2745.12742
Ballesteros-Mejia, L., Lima, N. E., Lima-Ribeiro, M. S., and Collevatti, R. G. (2016). Pollination mode and mating system explain patterns in genetic differentiation in Neotropical plants. PLoS One 11:e0158660. doi: 10.1371/journal.pone.0158660
Barros, F. M., Martello, F., Peres, C. A., Pizo, M. A., and Ribeiro, M. C. (2019). Matrix type and landscape attributes modulate avian taxonomic and functional spillover across habitat boundaries in the Brazilian Atlantic Forest. Oikos 128, 1600–1612. doi: 10.1111/oik.05910
Berec, L., Angulo, E., and Courchamp, F. (2007). Multiple Allee effects and population management. Trends Ecol. Evol. 22, 185–191. doi: 10.1016/j.tree.2006.12.002
Bolker, B and R Development Core Team, (2017). bbmle: Tools for General Maximum Likelihood Estimation. R package version 1.0.20. Available online at: https://CRAN.R-project.org/package=bbmle(accessed February 3, 2020).
Bovo, A. A. A., Ferraza, K. M. P.M. B., Magiolia, M., Alexandrino, E. R., Hasui, E., Ribeiro, M. C., and Tobias, J. A. (2018). Habitat fragmentation narrows the distribution of avian functional traits associated with seed dispersal in tropical forest. Perspect. Ecol. Conserv. 16, 90–96. doi: 10.1016/j.pecon.2018.03.004
Braga, A. C., and Collevatti, R. G. (2011). Temporal variation in pollen dispersal and breeding structure in a bee-pollinated Neotropical tree. Heredity 106:911. doi: 10.1038/hdy.2010.134
Braga, A. C., Reis, A. M. M., Leoi, L. T., Pereira, R. W., and Collevatti, R. G. (2007). Development and characterization of microsatellite markers for the tropical tree species Tabebuia aurea (Bignoniaceae). Mol. Ecol. Notes 7, 53–56. doi: 10.1111/j.1471-8286.2006.01521.x
Breheny, P., and Burchett, W. (2017). Visualization of regression models using visreg. R J. 9, 56–71.
Bijlsma, R., and Loeschcke, V. (2012). Genetic erosion impedes adaptive responses to stressful environments. Evol. Appl. 5, 117-129. doi: 10.1111/j.1752-4571.2011.00214.x
Burnhan, K. P., and Anderson, D. R. (2002). Model Selection and Multimodel Inference: An Information-Theoretic Approach. New York, NY: Springer
Carvalho, C. S., Ribeiro, M. C., Côrtes, M. C., Galetti, M., and Collevatti, R. G. (2015). Contemporary and historic factors influence differently genetic differentiation and diversity in a tropical palm. Heredity 115, 216–224. doi: 10.1038/hdy.2015.30
Collevatti, R. G., Grattapaglia, D., and Hay, J. D. (2001). Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci. Mol. Ecol. 10, 349–356. doi: 10.1046/j.1365-294x.2001.01226.x
Collevatti, R. G., Lima-Ribeiro, M. S., Terribile, L. C., Guedes, L. B., Rosa, F. F., and Telles, M. P. (2014). Recovering species demographic history from multi-model inference: the case of a Neotropical savanna tree species. BMC Evol. Biol. 14:213. doi: 10.1186/s12862-014-0213-0
Corlett, R. T. (2017). Frugivory and seed dispersal by vertebrates in tropical and subtropical Asia: an update. Glob. Ecol. Conserv. 11:1e22.
Damschen, E. I., Baker, D. V., Bohrer, G., Nathan, R., Orrock, J. L., Turner, J. R., Brudvig, L. A., Haddad, N. M., Levey, D. J., and Tewksbury, J. J. (2014). How fragmentation and corridors affect wind dynamics and seed dispersal in open habitats. Proc. Natl. Acad. Sci. U.S.A. 111, 3484–3489. doi: 10.1073/pnas.1308968111
Dixo, M., Metzger, J. P., Morgante, J. S., and Zamudio, K. R. (2009). Habitat fragmentation reduces genetic diversity and connectivity among toad populations in the Brazilian Atlantic Coastal Forest. Biol. Conserv. 142, 1560–1569. doi: 10.1016/j.biocon.2008.11.016
Do, C., Waples, R. S., Peel, D., Macbeth, G. M., Tillett, B. J., and Ovenden, J. R. (2014). NeEstimator V2: re-implementation of software for the estimation of contemporary effective population size (Ne) from genetic data. Mol. Ecol. Resour. 14, 209–214. doi: 10.1111/1755-0998.12157
Driscoll, D. A., Banks, S. C., Barton, P. S., Lindenmayer, D. B., and Smith, A. L. (2013). Conceptual domain of the matrix in fragmented landscapes. Trends Ecol. Evol. 28, 605–613. doi: 10.1016/j.tree.2013.06.010
Duarte, G. T., Santos, P. M., Cornelissen, T. G., Ribeiro, M. C., Paglia, A. P. (2018). The effects of landscape patterns on ecosystem services: meta-analyses of landscape services. Lands. Ecol. 33, 1247–1257. doi: 10.1111/j.1461-0248.2008.01157.x
Fahrig, L. (2003). Effects of habitat fragmentation on biodiversity. Ann. Rev. Ecol. Evol. Syst. 34, 487–515.
Frankie, G. W., Opler, P. A., and Bawa, K. S. (1976). Foraging behavior of solitary bees: implications for outcrossing of a Neotropical forest tree species. J. Ecol. 64, 1049–1057.
Frankham, R., Lees, K., Montgomery, M. E., England, P. R., Lowe, E. H., and Briscoe, D. A. (1999). Do population size bottlenecks reduce evolutionary potential? Anim. Conserv. 2, 255–260. doi: 10.1111/j.1469-1795.1999.tb00071.x
Ganeshaiah, K. N., and Shaanker, R. U. (1991). Seed size optimization in a wind dispersed tree Butea monosperma: a trade-off between seedling establishment and pod dispersal efficiency. Oikos 60, 3–6.
Gilpin, M. E., and Soulé, M. E. (1986). “Minimum viable populations: processes of species extinction,” in: Conservation Biology: The Science of Scarcity and Diversity. ed. Soule, [[Au Query:]] ìM. E. (Sunderland, MA: Sinauer Associates Inc), 19–34.
Goldstein, G., Meinzer, F. C., Bucci, S. J., Scholz, F. G., Franco, A. C., and Hoffmann, W. A. (2008). Water economy of Neotropical savanna trees: six paradigms revisited. Tree Physiol. 28, 395–404. doi: 10.1093/treephys/28.3.395
Goulson, D., and Stout, J. C. (2001). Homing ability of the bumblebee Bombus terrestris (Hymenoptera: Apidae). Apidologie 32, 105–111. doi: 10.1051/apido:2001115
González-Fernández, A., Arroyo-Rodríguez, V., Ramírez-Corona, F., Manjarrez, J., Aguilera-Hernández, A., and Sunny, A. (2019). Local and landscape drivers of the number of individuals and genetic diversity of a microendemic and critically endangered salamander. Lands. Ecol. 34, 1989–2000. doi: 10.1007/s10980-019-00871-2
Goudet, J. (2002). FSTAT, a Program to Estimate and Test Gene Diversities and Fixation Indices (Version 2.9.3.2). https://www2.unil.ch/popgen/softwares/fstat.htm
GRASS Development Team, (2018). Geographic Resources Analysis Support System (GRASS) Software, Version 7.4. Beaverton, OR: Open Source Geospatial Foundation.
Greene, D. F., and Johnson, E. A. (1993). Seed mass and dispersal capacity in wind-dispersed diaspores. Oikos 67, 69–74.
Hadley, A. S., and Betts, M. G (2011). The effects of landscape fragmentation on pollination dynamics: absence of evidence not evidence of absence. Biol. Rev. 87, 526–544. doi: 10.1111/j.1469-185X.2011.00205.x
Hamilton, M. B., and Miller, J. R. (2002). Comparing relative rates of pollen and seed gene flow in the island model using nuclear and organelle measures of population structure. Genetics 162, 1897–1909.
Hansen, T. F., Pélabon, C., and Houle, D. (2011). Heritability is not evolvability. Evol. Biol. 38:258. doi: 10.1007/s11692-011-9127-6
Hardy, O. J., and Vekemans, X. (2002). SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol. Ecol. Notes 2, 618–620. doi: 10.1046/j.1471-8286.2002.00305.x
Hedrick, P. W. (2005). A standardized genetic differentiation measure. Evolution 59, 1633–1638. doi: 10.1111/j.0014-3820.2005.tb01814.x
Hoffmann, W., and Franco, A. C. (2003). Comparative growth analysis of tropical forest and savanna woody plants using phylogenetically independent contrasts. J. Ecol. 91, 475–484. doi: 10.1046/j.1365-2745.2003.00777.x
Hoffmann, W. A., Orthen, B., and Franco, A. C. (2004). Constraints to seedling success of savanna and forest trees across the savanna-forest boundary. Oecologia 140, 252–260. doi: 10.1007/s00442-004-1595-2
Hoffmann, W. A., Geiger, E. L., Gotsch, S. G., Rossatto, D. R., Silva, L. C., and Lau, O. L. et al (2012). Ecological thresholds at the savanna-forest boundary: how plant traits, resources and fire govern the distribution of tropical biomes. Ecol. Lett. 15, 759–768. doi: 10.1111/j.1461-0248.2012.01789.x
Holderegger, R., Kamm, U., and Gugerli, F. (2006). Adaptive versus neutral genetic diversity: implications for landscape genetics. Landsc. Ecol. 21, 797–807. doi: 10.1101/gr.104331.109
Holderegger, R., Buehler, D., Gugerli, F., and Manel, S. (2010). Landscape genetics of plants. Trends Plant Sci. 15, 675–683. doi: 10.1016/j.tplants.2010.09.002
Jackson, N. D., and Fahrig, L. (2016). Habitat amount, not habitat configuration, best predicts population genetic structure in fragmented landscapes. Landsc. Ecol. 31, 951–968. doi: 10.1007/s10980-015-0313-2
Jost, L. (2008). GST and its relatives not measure differentiation. Mol. Ecol. 17, 4015–4026. doi: 10.1111/j.1365-294x.2008.03887.x
Kimura, M. (1983). The Neutral Theory of Molecular Evolution. Cambridge: Cambridge University Press.
Kitamoto, N., Ueno, S., Takenaka, A., Tsumura, Y., Washitani, I., and Ohsawa, R. (2006). Effect of flowering phenology on pollen flow distance and the consequences for spatial genetic structure within a population of Primula sieboldii (Primulaceae). Am. J. Bot. 93, 226–233. doi: 10.3732/ajb.93.2.226
Lande, R. (1988). Genetics and demography in biological conservation. Science 241, 1455–1460 doi: 10.1126/science.3420403
Lahoreau, G., Barot, S., Gignoux, J., Hoffmann, W. A., Setterfield, S. A., and Williams, P. R. (2006). Positive effect of seed size on seedling survival in fire-prone savannas of Australia, Brazil and West Africa. J. Trop. Ecol. 22, 719–722. doi: 10.1017/s026646740600349x
Latrubesse, E. M., Arima, E., Ferreira, M. E., Nogueira, S. H., Wittmann, F., Dias, M. S., Dagosta, F. C. P., and Bayer, M. (2019). Fostering water resource governance and conservation in the Brazilian Cerrado biome. Conserv. Sci. Pract. 1:e77.
Leinonen, T., Cano, J. M., Makinen, H., and Merila, J. (2006). Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks. J. Evol. Biol. 19, 1803–1812. doi: 10.1111/j.1420-9101.2006.01182.x
Lion, M. B., Garda, A. A., and Fonseca, C. R. (2014). Split distance: a key landscape metric shaping amphibian populations and communities in forest fragment. Divers. Distrib. 20, 1245–1257. doi: 10.1111/ddi.12228
McGarigal, K., Cushman, S. A., and Ene, E. (2012). FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical and Continuous Maps. Computer Software Program. Amherst, MA: University of Massachusetts.
McKay, J. K., and Latta, R. G. (2002). Adaptive population divergence: markers, QTL and traits. Trends Ecol. Evol. 17, 285–291. doi: 10.1016/s0169-5347(02)02478-3
Mendes, C. P., Ribeiro, M. C., and Galetti, M. (2016). Patch size, shape and edge distance influence seed predation on a palm species in the Atlantic forest. Ecography 39, 465–475. doi: 10.1111/ecog.01592
Meyers, L. A., and Bull, J. J. (2002). Fighting change with change: adaptive variation in an uncertain world. Trends Ecol. Evol. 17, 551–557. doi: 10.1016/s0169-5347(02)02633-2
Moraes, A. M., Ruiz-Miranda, C. R., Galetti, P. M., Niebuhr, B. B., Alexandre, B. R., Muylaerte, R. L., et al. (2018). Landscape resistance influences effective dispersal of endangered golden lion tamarins within the Atlantic Forest. Biol. Conserv. 224, 178–187. doi: 10.1016/j.biocon.2018.05.023
Mousadik, A., and Petit, R. J. (1996). High level of genetic differentiation for allelic richness among populations of the argan tree [Argania spinosa (L.) Skeels] endemic to Morocco. Theor. Appl. Genet. 92, 832–839. doi: 10.1007/BF00221895
Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., and Kent, J. (2000). Biodiversity hotspots for conservation priorities. Nature 403, 853–858. doi: 10.1038/35002501
Nei, M. (1978). Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89, 583–590.
Nomura, T. (2008). Estimation of effective number of breeders from molecular coancestry of single cohort sample. Evol. Appl. 1, 462–474. doi: 10.1111/j.1752-4571.2008.00015.x
van Oosterhout, C. V., Hutchinson, W. F., Wills, D. P. M., and Shipley, P. (2004). MICRO- CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes 4, 535–538. doi: 10.1111/j.1471-8286.2004.00684.x
Perez-Harguindeguy, N., Dıaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234.
Pivello, V. R. (2011). The use of fire in the Cerrado and Amazonian rainforests of Brazil. Fire Ecol. 7, 24–39. doi: 10.4996/fireecology.0701024
Poorter, H., Niinemets, U., Poorter, L., Wright, I. J., and Villar, R. (2009). Causes and consequences of variation in leaf mass per area (LMA): a meta-analysis. New Phytol. 182, 565–588. doi: 10.1111/j.1469-8137.2009.02830.x
Prout, T., and Barker, J. F. S. (1993). F statistics in Drosophila buzzatii: selection, population size and inbreeding. Genetics 134, 369–375.
R Core Team, (2019). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Reed, D. H., and Frankham, R. (2002). Correlation between fitness and genetic diversity. Conserv. Biol. 17, 230–237. doi: 10.1046/j.1523-1739.2003.01236.x
Regolin, A. L., Cherem, J. J., Graipel, M. E., Bogoni, J. A., Ribeiro, J. W., Vancine, M. H., et al. (2017). Forest cover influences occurrence of mammalian carnivores within Brazilian Atlantic Forest. J. Mammal. 98, 1721–1731. doi: 10.1093/jmammal/gyx103
Reich, P. B., Walters, M. B., and Ellsworth, D. S. (1992). Leaf life-span in relation to leaf, plant, and stand characteristics among diverse ecosystems. Ecol. Monogr. 62, 365–392. doi: 10.2307/2937116
Reich, P. B. (2000). Do tall trees scale physiological heights? Trends Ecol. Evol. 15, 41–42. doi: 10.1016/s0169-5347(99)01734-6
Resende, M. D. V. (2016). Software Selegen-REML/BLUP: a useful tool for plant breeding. Crop Breed. Appl. Biotechnol. 16, 330–339. doi: 10.1590/1984-70332016v16n4a49
Ricketts, T. H. (2001). The matrix matters: effective isolation in fragmented landscapes. Am. Natural. 158, 87–99. doi: 10.1086/320863
Sano, E. E., Rosa, R., Brito, J. L. S., and Ferreira, L. G. (2010). Land cover mapping of the tropical savanna region in Brazil. Environ. Monitor. Assess. 166, 113–124. doi: 10.1007/s10661-009-0988-4
Santos, P. M., Chiarello, A. G., Ribeiro, M. C., RIbeiro, J. W., and Paglia, A. P. (2016)Local and landscape influences on the habitat occupancy of the endangered maned sloth Bradypus torquatus within fragmented landscapes. Mammal. Biol. 81, 447–454. doi: 10.1016/j.mambio.2016.06.003
Schmidt, T., Arens, P., Smulders, M. J. M., Billeter, R., Liira, J., Augenstein, I., and Durk, W. (2009). Effects of landscape structure on genetic diversity of Geum urbanum L populations in agricultural landscapes. Flora 204, 549–559. doi: 10.1016/j.flora.2008.07.005
Shamoon, H., Maor, R., Saltz, D., and Dayan, T. (2018). Increased mammal nocturnality in agricultural landscapes results in fragmentation due to cascading effects. Biol. Conserv. 226, 32–41. doi: 10.1016/j.biocon.2018.07.028
Silva, L. G., Ribeiro, M. C., Hasui, E., Costa, C. A., and Cunha, R. G. T. (2015). Patch size, functional isolation, visibility and matrix permeability influences Neotropical Primate occurrence within highly fragmented landscapes. PLoS One 10:e0114025. doi: 10.1371/journal.pone.0114025
Slatkin, M. (1995). A measure of population subdivision based on microsatellite allele frequencies. Genetics 139, 457–462.
Soons, M. B., Heil, G. W., Nathan, R., and Katul, G. G. (2004). Determinants of long-distance seed dispersal by wind in grasslands. Ecology 85, 3056–3068. doi: 10.1890/03-0522
Soons, M. B., and Heil, G. W. (2002). Reduced colonization capacity in fragmented populations of wind-dispersed grassland forbs. J. Ecol. 90, 1033–1043. doi: 10.1046/j.1365-2745.2002.00729.x
Soons, M. B., Messelink, J. H., Jongejans, E., and Heil, G. W. (2005). Habitat fragmentation reduces grassland connectivity for both short-distance and long-distance wind-dispersed forbs. J. Ecol. 93, 1214–1225. doi: 10.1111/j.1365-2745.2005.01064.x
Spitze, K. (1993). Population structure in Daphnia obtusa: quantitative genetic and allozymic variation. Genetics 135, 367–374.
Tackenberg, O., Poschlod, P., and Bonn, S. (2003). Assessment of wind dispersal potential in plant species. Ecol. Monogr. 73, 191–205. doi: 10.1890/0012-9615(2003)073%5B0191:aowdpi%5D2.0.co;2
Trakhtenbrot, A., Katulb, G. G., and Nathan, R. (2014). Mechanistic modeling of seed dispersal by wind over hilly terrain. Ecol. Model. 274, 29–40. doi: 10.1016/j.ecolmodel.2013.11.029
Tambarussi, E. V., Pereira, F. B., Silva, P. H. M., Lee, D., and Bush, D. (2018). Are tree breeders properly predicting genetic gain? A case study involving Corymbia species. Euphytica 214:150.
Tomlinson, K. W., Sterck, F. J., Bongers, F., Silva, D. A., Barbosa, E. R. M., Ward, D., Bakker, F. T., et al. (2012). Biomass partitioning and root morphology of savanna trees across a water gradient. J. Ecol. 100, 1113–1121. doi: 10.1111/j.1365-2745.2012.01975.x
Vencovsky, R., Pereira, M. B., Crisóstomo, J. R., and Ferrerira, M. A. J. F. (2001). Genética e melhoramento de populações mistas. in Recursos Genéticos e Melhoramento: Plantas, eds Nass, LL, Valois, ACC, Melo, IS, Valadares-Inglis, MC (Fundação MT: Rondonópolis), 231–281.
Vencovsky, R., and Crossa, J. (2003). Measurements of representativeness used in genetic resources conservation and plant breeding. Crop Sci. 43,1912–1921. doi: 10.2135/cropsci2003.1912
Wagner, H., and Fortin, M. J. (2013). A conceptual framework for the spatial analysis of landscape genetic data. Conserv. Genet. 14, 253–261. doi: 10.1007/s10592-012-0391-5
Weir, B. S., and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. Evolution 38, 1358–1370. doi: 10.1111/j.1558-5646.1984.tb05657.x
Westoby, M., Leishman, M., and Lord, J. (1996). Comparative ecology of seed size and dispersal. Philos. Transact. R. Soc. Lond. B 351, 1309–1318. doi: 10.1098/rstb.1996.0114
Westoby, M. (1998). A leaf-height-seed (LHS) plant ecology strategy scheme. Plant Soil 199, 213–227.
Westoby, M., Falster, D., Moles, A., Vesk, P., and Wright, I. (2002). Plant ecological strategies: some leading dimensions of variation between species. Ann. Rev. Ecol. Syst. 33, 125–159. doi: 10.1146/annurev.ecolsys.33.010802.150452
Keywords: agroecosystem, Bignoniaceae, Cerrado, fragmentation, genetic diversity, landscape genetics, quantitative genetics, Tabebuia aurea
Citation: Collevatti RG, dos Santos JS, Rosa FF, Amaral TS, Chaves LJ and Ribeiro MC (2020) Multi-Scale Landscape Influences on Genetic Diversity and Adaptive Traits in a Neotropical Savanna Tree. Front. Genet. 11:259. doi: 10.3389/fgene.2020.00259
Received: 30 October 2019; Accepted: 04 March 2020;
Published: 25 March 2020.
Edited by:
Andrew Shirk, University of Washington, United StatesReviewed by:
Evandro Vagner Tambarussi, Universidade Estadual do Centro-Oeste do Paraná, BrazilEduardo Rogério Moribe Barbosa, University of Brasilia, Brazil
Copyright © 2020 Collevatti, dos Santos, Rosa, Amaral, Chaves and Ribeiro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rosane Garcia Collevatti, rosanegc68@hotmail.com