Original Research ARTICLE
Anthropogenic Disturbances Eroding the Genetic Diversity of a Threatened Palm Tree: A Multiscale Approach
- 1Applied Ecology and Conservation Lab, Programa de Pós-Graduação em Ecologia e Conservação da Biodiversidade, Universidade Estadual de Santa Cruz, Ilhéus, Brazil
- 2Laboratório de Marcadores Moleculares, Centro de Biotecnologia e Genética, Universidade Estadual de Santa Cruz, Ilhéus, Brazil
Habitat loss and the illegal exploitation of natural resources are among the main drivers of species extinction around the world. These disturbances act at different scales, once changes in the landscape composition and configuration operate at large scales and exploitation of natural resources at local scales. Evidence suggests that both scales are capable of triggering genetic erosion in the remaining populations. However, most of the studies so far did not evaluate simultaneously the effects of these disturbances on genetic diversity and structure of plants. In this study, we used a multiple scale approach to empirically evaluate the impacts caused by local and landscape scale disturbances in the genetic diversity and structure of an endangered palm tree, Euterpe edulis. We sampled and genotyped with microsatellite markers 544 juveniles of E. edulis in 17 fragments of Atlantic Forest in Brazil. In addition, we estimated the local logging rate and the forest cover and isolation at landscape scale. We found that the palm populations have not undergone any recent bottleneck events and that only logging intensification had affected the fixation index and the number of private alleles. Additionally, we did not detect any evidence of spatial genetic structure or genetic divergence associated with environmental disturbance variables at different scales. However, we identified distinct genetic clusters, which may indicate a reduction of gene flow between fragments that were previously a continuous habitat. Our results show that local disturbances, which act directly on population size reduction, such as logging, modified the genetic diversity more rapidly, whereas genetic structure is probably more influenced by large-scale modifications. In this way, to maximize the conservation efforts of economically exploited species, we recommend to increase the inspection to reduce the illegal exploitation, and reforestation of degraded areas, in order to increase the gene flow in Atlantic Forest fragments.
The conversion of natural environments into anthropogenic landscapes is one of the main causes of biodiversity loss worldwide (Laurance et al., 2014; Newbold et al., 2015; Miraldo et al., 2016). Mounting evidence shows the pervasive effects of habitat loss on several taxonomic groups (Mortelliti et al., 2010; Ferreira et al., 2015; Muylaert et al., 2016; Rocha-Santos et al., 2017). The remaining forest patches are being reduced and isolated with negative consequences for species diversity. In this context, the remaining populations are prone to gene flow decrease, alleles loss due to genetic drift, and inbreeding depression. The consequence of these genetic changes might be a lower ability to adapt to further environmental changes (Bijlsma and Loeschcke, 2012; Fountain et al., 2016; Rhoads et al., 2017; Browne and Karubian, 2018). Several studies have corroborated this prediction (Dixo et al., 2009; Zhang et al., 2012; Wood et al., 2017), however, the intensity and the velocity of the response to impacts vary according to life history and environmental characteristics (Vranckx et al., 2012; Lino et al., 2019).
In human-modified landscapes, species of economic value are also locally vulnerable due to (i) the direct removal of individuals or (ii) the indirect exploitation of resources, such as seeds, impairing the regeneration capacity of these species (Homma, 1992; Peres et al., 2003). Acute disturbances, such as logging, are capable of causing severe demographic depletion over generations (Richardson and Peres, 2016) and increase the risk of genetic erosion in the remaining populations (Shaanker et al., 2003; Shaanker et al., 2004; Dai et al., 2018). Theoretical and empirical studies have observed reductions in allele number and heterozygosity and an increase in spatial genetic distance in exploited populations (Cruse-Sanders et al., 2005;Lacerda et al., 2008; Sebbenn et al., 2008). However, investigations that simultaneously assessed the individual impact caused by landscape and local scale disturbances are scarce (but see González-Fernández et al., 2019).
The combination of anthropogenic pressures at the landscape and local scale is striking in the tropics. The tropics stand out worldwide due to the high rates of species biodiversity and endemism, but the increasing human population size and forest loss is a growing threat (Hansen et al., 2013; Lewis et al., 2015). In this scenario, understanding the individual impact of different types of disturbance on the genetic diversity of exploited species could better direct conservation efforts. For this reason, in this study, we used a multi-level approach to investigate the effects of anthropogenic disturbances at the landscape and local scale on the genetic diversity and structure of Euterpe edulis populations. This species is native to the Atlantic Rainforest (Reis et al., 2000a), a biome that currently conserves only 12% of its original coverage (Ribeiro et al., 2009). Moreover, this palm is in the list of the Brazilian species threatened with extinction because of the population decline recorded in the last decades (Martinelli and Moraes, 2013). The species has great economic importance due to the illegal harvest to the commercialization of the apical meristem (Galetti and Fernandez, 1998; Matos and Bovi, 2002). The harvest culminates in the death of individuals, as there is no regrowth after cutting (Ferri and Cavalcante, 1997). In the tropics, changes in the landscape composition and configuration, such as forest cover loss and isolation provide more accessibility to forest resources, increasing illegal hunt and logging (Tabarelli et al., 2004). These activities might in turn, negatively affect seed dispersal by vertebrates and seedling recruitment (Peres and Palacios, 2007; Gutiérrez-Granados, 2011). Using a multiscale inference approach, we evaluated how landscape composition and configuration, measured at several spatial scales, and local disturbances, affect the genetic diversity of E. edulis. We predicted that estimates of genetic diversity and structure would be strongly influenced by the synergetic effects of disturbances acting at different scales. Thus, models containing a combination of landscape and local variables would better explain the genetic diversity of the species.
Materials and Methods
Study Area and Sites Selection
The study region was located in the Atlantic Rainforest in Southern Bahia State, Brazil. Deforestation in the region started in the mid-1980’s and was intensified in the 1990s during the cocoa crisis, the main economic product at that time (Rocha, 2006). The region presents some of the last remnants of the northeastern Atlantic Forest (Araújo et al., 1998) and still harbors a large number of flora and fauna species, including several endemics (Thomas et al., 1998). The land use history resulted in a mosaic of forest patches in different successional stages immersed in a matrix of pasture, rubber and eucalyptus plantations, and cacao agroforests (Landau, 2003; Sambuichi, 2003; Landau et al., 2008). The native vegetation is tropical lowland rainforest and the climate is classified as Af, according to Köppen, hot and moist, without a distinct dry season (Gouvêa, 1969).
After an intense process of ground-truthing, we mapped the land-use of 3,470 km2 using high-resolution satellite images (i.e. RapidEye from 2009-2010, QuickBird and WorldView both from 2011; with resolutions of 5, 0.5, and 0.6 m, respectively). We avoided to sample in Montane and Restinga Forests, which resulted in a subdivision into two regions (north and south, Figure 1A) mainly due to a sandy stripe between them. The history of deforestation and land use are different between the two regions. The matrix of the north region is more heterogeneous and in the south region the matrix is mainly dominated by pastures. However, the two regions show similar soil, topography and vegetation types (Thomas et al., 1998). The main land-use classes were classified as: forests (mature and secondary native forest), cattle pasture, and plantations of cacao, rubber tree, and Eucalyptus sp.) From the map, we selected 58 possible forest fragments adopting as exclusion criteria highly demanding access areas, and indigenous lands. Then, we used previous studies in the region to select a subset of 17 forest fragments, at least 2 km apart, which spanned a large range of landscape-scale forest cover (Table S1), and for which the occurrence of E. edulis was known (Santos et al., 2015; Soares et al., 2015; Benchimol et al., 2017).
Figure 1 Location of the 17 forest fragments in which the populations of Euterpe edulis were sampled in the Atlantic Forest of southern Bahia. (A) Distribution of fragment center points (circle) and characterization of the main forms of land use in the region; (B) Example of the seven buffers created around the fragment’s center point for landscape metric calculations.
Our sampling occurred between 2014 and 2016, and we used Google Earth to evaluate substantial changes in land use cover during the period between mapping (2009–2011) and data collection. After finding stability in land use cover during this period, we performed all analyzes with forest cover calculations obtained by the high-resolution mapping.
Focal Species and Sampling
Euterpe edulis is a monoecious palm tree, with annual and predominantly crossed reproduction (Mantovani and Morellato, 2000; Gaiotto et al., 2003; Castro et al., 2007). The flowers are abundant in nectar and pollen and attract a wide variety of insects (Reis et al., 2000b). Despite this, its pollination is mainly performed by bees belonging to different taxonomic groups such as Meliponini, Euglossini, and Honeybees (Reis et al., 2000b; Dorneles et al., 2013; Santos et al., 2018a; Santos et al., 2018b). E. edulis fruits have a pulp rich in lipid and fiber (Silva et al., 2013; Da Silva et al., 2014) and are used as food by 58 bird species and 20 mammal species (Galetti et al., 2013). However, the main seed dispersers are large frugivores such as toucans and cotingas. Small birds such as thrushes are especially important in defaunated areas (Reis and Kageyama, 2000; Galetti et al., 2013; Santos et al., 2018a; Santos et al., 2018b).
We randomly established three 50 × 10 forest-plots in each of the 17 forest fragments. Plots were located at least 50 m apart from each other and from the nearest edge. In each forest plot, we sampled all E. edulis juvenile (i.e. individuals with pinnate leaves, height ≤0.15 and >1.00 m) based on the categories proposed by Silva et al. (2009). Subsequently, all individuals were georeferenced and had a leaf tissue sampled. We chose the juvenile ontogenetic stage because the genetic parameters of the early developmental stages better represent the effects of recent environmental disturbances (Vranckx et al., 2012). In addition, juveniles might show more consistent responses to anthropogenic disturbances because they exhibit higher numerical stability in terms of population fluctuation than seedlings (Conte et al., 2000). All collected individuals were numbered and later, 32 individuals were drawn in each of the 17 population for genotyping. This procedure resulted in a total sampling of 544 individuals. The determination of the number of individuals per population was based on a previous study that detected little benefit on the accuracy of allele frequencies and diversity estimates for microsatellite loci above 30 individuals (Hale et al., 2012).
DNA Extraction and Genotyping
DNA extraction followed the CTAB protocol (Doyle and Doyle, 1987). We genotyped all E. edulis individuals using 17 nuclear microsatellite markers developed for the species (Gaiotto et al., 2001). Multiplex PCR reactions were performed in a Veriti™ Thermal Cycler (Applied Biosystems, Foster City, CA) with two triplex combinations (EE43, EE45, EE52 and EE47, EE59, EE63) and two duplex combinations (EE2, EE32 and EE8, EE23). We also performed single locus PCRs to the following markers: EE3; EE5; EE9; EE15; EE25; EE48; and EE54. Subsequently, we submitted only the EE5 marker product to the individual electrophoresis system, while the PCR products of the other markers were submitted to the multiload electrophoresis system adapted from Gaiotto et al. (2003) and were organized with the following combinations: (i) Pentaload I, (EE2, EE32, EE8, EE23, EE3); (ii) Pentaload II (EE43, EE45, EE52, EE48, EE54); (iii) Tetraload (EE47, EE59, EE63, EE25); and (iv) Biload (EE9 and EE15). The mix submitted to the electrophoresis was composed of 2 µl of the PCR product (or of the mix of the PCR products for the cases of the multiload system), 0.2 µl of GeneScan™ 500 Liz™ (Applied Biosystems, Thermo-Fisher Scientific, Inc., Waltham, MA, USA) and 7.8 µl of deionized formamide (Applied Biosystems). Later, we performed genotyping on the ABI 3500 Genetic Analyzer (Applied Biosystems, Foster City, CA). The sizing of fragments was obtained with GeneMarker® software 2.2 (SoftGenetics, State College, PA, USA).
We evaluated the Hardy–Weinberg equilibrium within each sample population with the package DiveRsity (Keenan et al., 2013) on the R environment 3.5.2 (R Core Team, 2018) and tested the linkage disequilibrium between all pairs of loci through FSTAT 2.9.3 (Goudet, 2001). The null alleles were corrected using the method of Oosterhout (van Oosterhout et al., 2006), in MICRO-CHECKER software v. 2.2.3 (van Oosterhout et al., 2004) after the results of 999 bootstrap simulations with a confidence interval of 95%. After null alleles correction, we perform the estimations of all genetic parameters. The DiveRsity package was used to calculate observed and expected heterozygosity (HO and HE, respectively), allelic richness, and fixation index with a 95% confidence interval estimated by 999 bootstrap simulations. In addition, we used the GENALEX 6.5 program (Peakall and Smouse, 2012) to calculate the number of private alleles in each population. Considering that detection of private alleles may be biased by sample insufficiency, we acknowledge that this estimate might include alleles that occur exclusively in a single population but also those presented in low frequencies (<0.05) that could not be detected. Finally, we used the software BOTTLENECK v.1.2.02 (Cornuet and Luikart, 1996; Piry et al., 1999) to investigate the occurrence of recent genetic bottlenecks in the populations sampled. We chose the two-phase mutation model, due to it is the most recommended for microsatellite data (Piry et al., 1999; Williamson-Natesan, 2005). We fitted this model with 95% of single-step mutations and 5% of mutations of multiple steps, as suggested by Piry et al. (1999). Subsequently we used the Wilcoxon test with 5,000 iterations to evaluate the occurrence of excess heterozygosity.
Landscape Metrics and Local Variable
To evaluate whether the genetic diversity of E. edulis is influenced by anthropogenic disturbances, we related the genetic estimates with landscape and local attributes. At the landscape scale, we evaluated metrics of landscape composition and configuration that might affect genetic diversity and structure (Balkenhol et al., 2013; Jackson and Fahrig, 2016). These landscape attributes were calculated using the program FRAGSTATS v126.96.36.1993 (McGarigal et al., 2002) for seven buffers of different radii sizes, ranging from 0.5 to 2 km, each 250 m, from the central plot in each forest fragment (Figure 1B). The highest radius was chosen because it covers the foraging distances of important large seed dispersers of E. edulis, such as toucans, which have a medium-sized seed dispersal distance ranging from 269 to 449 m (Holbrook, 2011). In addition it also includes different foraging distances reported for potential pollinators of E. edulis, such as Plebeia doryana (maximum 540 m) and Apis mellifera that can reach distances of hundreds of meters of foraging (Zurbuchen et al., 2010). We evaluated (i) forest cover (a proxy of habitat amount) as landscape composition descriptor and (ii) mean patch size, proximity index, and edge density as landscape configuration descriptors. For all scales, forest cover was calculated as the area occupied by the sum of mature and secondary native forests divided by the total landscape area. We used mature and secondary forests as a proxy of habitat because natural populations of E. edulis occurred in both types of forests. The calculation of all landscape configuration metrics was performed at the class level and was based only on the areas occupied by forest remnants within each landscape. Then, we evaluated the Pearson correlation among all landscape variables. We excluded from the subsequent analyses mean patch size and edge density due to the high correlation with forest cover amount in more than one spatial scale (r > 0.4) (Table S2).
At the local scale, we estimated logging activity as a disturbance variable related to E. edulis harvest. We used logging as a proxy of harvesting because both activities are intensified in anthropogenic landscapes and close to urban centers where they become one of the main income sources for local populations (Shaanker et al., 2003; Tabarelli et al., 2004). We also opt to estimate logging because the record of the illegal palm harvest is very scarce due to the rapid stipe decomposition when compared with hardwood trees (personal observation). The logging estimative in each fragment was made in 2014 when we evaluated a 100 × 8 m plot, distant at least 50 m from the closest edge. In all plots, we counted the wood stumps with the diameter at ground level ≥20 cm.
To evaluate the spatial scale to which the environmental context influenced the genetic diversity (scale of effect) we used the multifit function (Huais, 2018) and generalized additive models (GAM: Hastie and Tibshirani, 1990). We related simultaneously the genetic attributes to the forest cover and proximity index quantified on seven spatial scales (buffers with 0.5, 0.75, 1.0, 1.5, 1.25, 1.75, 2.0 km radius). In addition, we evaluated the spatial autocorrelation of the residual models using the selected scale by applying a Moran I index (Fortin and Dale, 2005) (Table S3). We did not find any spatial autocorrelation in our data (Supplementary Material), hence it was not necessary to consider the spatial structure in further analyses.
We performed GAM to capture linear and nonlinear relationships between variables. All possible additive combination of independent variables (landscape and local attributes) and the genetic parameters were tested, totaling seven models. We tested for the concurvity among predictors i.e. the nonparametric counterpart of multicolinearity in linear regressions. We also included a null model containing only the intercept. Finally, we performed a model selection approach using the Akaike Information Criterion corrected for small sample sizes (AICc; Burnham and Anderson, 2002) for each one of the genetic parameters. Models with ΔAIC values ≤ 2.00 were considered to be equally plausible. But we chose the most parsimonious model, using the criterion of at least weight among plausible models.
The genetic differentiation among populations of E. edulis was investigated through three distinct approaches: we used the DiveRsity package to calculate the Gst-statistics (Hedrick, 2005) and Wright’s FST The statistical significance of these two analysis was evaluated by 1,000 randomizations with a confidence interval of 95%. The genetic differentiation measures estimated by Gst values are more suitable for microsatellite locus than those generated by Fst because they standardize the bias generated by the polymorphic character of these markers (Hedrick, 2005). However, Meirmans and Hedrick (2011)recommend that Fst values be presented in current studies to enable comparisons with older studies. The third approach of the genetic differentiation used in our study was the discriminant analysis of principal components (DAPC: Jombart et al., 2010) performed by the Adegenet package in R (Jombart, 2008). This analysis creates a model in which the genetic variation is partitioned between groups and within groups with the objective of maximizing variation between groups (Jombart and Collins, 2015). We initially used the K-means algorithm of the “find cluster” function and retained all the major components to detect evidence of genetic clusters in our populations (Figure S1A). We then applied the Bayesian Information criterion (BIC) to infer the best number of clusters by identifying the lowest BIC value (Figure S1B). Subsequently, we used the cross-validation function (Xval.Dapc) to identify the ideal number of principal component (PCs) to be retained based on the least mean square error after 1,000 stratified random simulations (Figure S1C). After obtaining this result, we performed the DAPC, which in the first moment transformed the genetic data using principal component analysis and then submitted the number of retained PCs to linear DAPC.
Genetic Diversity Within Populations
Our results did not show evidence of genotypic linkage disequilibrium in any of the pairs of loci. In addition, we detected the presence of null alleles in all loci with frequency varying from 0.5 to 16%. The average of genetic parameters estimated with original or corrected data for null alleles lead to similar results (data not shown). However, we opted to show all results based on the corrected data for null alleles (Table 1). We detected a total number of 347 alleles in the 17 loci by the 544 juvenile individuals from all 17 E. edulis populations evaluated. We used highly polymorphic microsatellites and the number of alleles generated for each marker ranged from 7 (EE43) to 33 (EE52), with a mean of 20.4 alleles for loci.
We did not detect signal of bottleneck in any population, since the unicaudal probability for the excess of heterozygotes between microsatellite markers was not significant (P > 0.7) and the proportion of alleles in the frequency intervals followed an L distribution format (Table 1).
Environmental Landscape and Local Variable
The scale of effect that strongly influences the genetic parameters varied in relation to the response variable and also to the landscape metric (Table S1). The observed and expected allelic richness and heterozygosity were not affected by any landscape or local variable (null model selected; Table 2). In contrast, the number of private alleles and the coefficient of inbreeding were significantly affected by logging intensity within the fragment (Figure 2). The concurvity measures were very small in the selected models, suggesting negligible concurvity (estimate < 0.0005).
Table 2 Result of the multimodel inference using generalized additive models (GAM) for the different parameters of genetic diversity. In all cases, we include the null model, containing only the intercept and error parameters.
Figure 2 Effects of logging (number of stumps) on the genetic diversity of E. edulis populations. (A) private alleles; (B) fixation index (f). The gray area corresponds to the confidence interval of the generalized additive models.
Paired GST ranged from 0.077 to 0.563 (0.292 ± 0.095). All estimates of genetic differentiation given by this parameter were significantly different from zero (Figure 3). Estimated values for paired FST ranged from 0.022 to 0.19 (0.089 ± 0.032) and although genetic differentiation between populations was less pronounced than that presented by GST (Table S4), we also found a significant genetic structure pattern through this analyze. The ideal number of genetic groupings indicated by the lowest BIC in the DAPC analysis was 9 (K9 = 1,048.768) (Figure S1B). However, only three clusters were formed by a clear genetic differentiation (populations 6, 15 e 16) (Figure 4). The other groups were composed of different proportion of individuals that belonged to distinct geographic populations, which indicates a great mixture of the gene pool of those populations.
Figure 3 The GST values paired among the 17 populations of E. edulis located in the Atlantic Forest of Southern Bahia.
Figure 4 Discriminant analysis of principal components (DAPC) for 17 populations of E. edulis from the Atlantic Forest of Southern Bahia, Brazil. The dots symbolize individuals within the population and the circles represent the genetic clusters defined by the DAPC.
Multilevel approach studies have been increasing in landscape ecology literature in recent years (Graham and Blake, 2001; Buskirk, 2005; Calamari et al., 2018), but are still scarce in the genetic field (González-Fernández et al., 2019). Our study represent a step forward in filling this gap, since it simultaneously evaluated the impacts of human activities at landscape and local scales on the genetic diversity of E. edulis populations in the Atlantic Forest. Despite the demographic decline recorded in many localities of E. edulis occurrence (Galetti and Fernandez, 1998; Matos and Bovi, 2002), we did not find evidence of a genetic bottleneck in any of the populations sampled. However, we found that the number of effective alleles and the fixation index were influenced by logging activity. In addition, we found evidence of genetic structure among the populations investigated suggests a limitation on gene flow that operates at larger scales.
The geographic distribution of E. edulis was originally recorded in almost all Atlantic Forest biome, with high densities in subcanopy of dense ombrophilous forests (Reis et al., 2000a). The current threat status of E. edulis is vulnerable to extinction, mainly due to the 30% population decline registered in the last six decades (Martinelli and Moraes, 2013). In spite of that, our study revealed that the populations investigated showed no signs of recent bottlenecks and still retain high genetic diversity. This result is compatible with previous studies (Conte et al., 2008; Carvalho et al., 2015; Santos et al., 2015), and also corroborates the hypothesis that recent population declines are not sufficient to produce changes in genetic diversity levels (Klank et al., 2012; Montes et al., 2016). This is because although genetic bottleneck events result commonly in the loss of rare alleles—frequency < 0.05—(Allendorf, 1986), heterozygosity levels are poorly influenced since these alleles contribute little to the formation of heterozygous individuals (Hartl and Clark, 1997). Therefore, our results suggest that E. edulis populations reflect a historical condition of the original species distribution since it corresponds to the expected diversity pattern for populations with a broad geographic distribution (Hamrick et al., 1992).
Contrary to our initial expectations, we found that only local scale disturbances affect some estimates of genetic diversity. Previous studies suggested the importance of landscape scale effects on the genetic diversity. For instance, a simulation study by Jackson and Fahrig (2014) found that while species abundance is influenced by local scale disturbances, genetic diversity tends to be more affected by landscape composition and configuration at large spatial scales. In addition, Carvalho et al. (2015) assessing the effects of recent landscape changes on genetic diversity of the same species, E. edulis, found the resistance of the matrix as one of the main factors affecting the allelic richness of these populations. These authors conducted the study in the southeastern Atlantic Forest, thus a possible explanation for the lack of landscape effects on our study, may be attributed to differences in land-use history between regions. The fragmentation of the southeastern Atlantic Forest began with coffee plantations in the 19th century (Dean, 1996). In contrast, intensive deforestation in our region began only during the 1980s (Rocha, 2006). A meta-analysis by Schlaepfer et al. (2018) identified that the age of anthropogenic fragmentation is a determining factor to detect genetic diversity loss. Thus, we believe that the recent landscape modification in the southern Bahia still did not trigger a genetic diversity erosion in E. edulis populations.
On the other hand, we observed a faster negative effect of logging activity, on the number of private alleles and fixation index. The inverted U-shaped relationship between logging and the number of private alleles showed that both areas with low and high logging activity had few alleles, which are unshared among populations. One possible explanation for this pattern is that the small number of private alleles in the poorly explored areas reflected the historical condition of the geographic distribution of E. edulis, with great populations and high level of gene flow (Gaiotto et al., 2003, Santos et al., 2015). In contrast, the low number of private alleles in intensively exploited areas was probably the result of genetic drift, which is a major cause of random loss of alleles in reduced populations (Frankham et al., 2002).
The relationship between logging and fixation index was expressed by a non-linear relationship, evidencing an increasing inbreeding with logging intensification. The results obtained are in agreement with the expected pattern for the traditional harvest system of this species that compromises the demographic structure of the populations by leaving only a few reproductive individuals in the fragments. (Reis et al., 2000a). In these scenarios an inbred pattern is common because the natural regeneration of populations is associated with the presence of few matrices (Murawski et al., 1994). On the other hand, in efficient forest restoration practices in which several matrices of distinct geograpthical origins is used the inbred patterns can be diluted (Zucchi et al., 2018). The genetic structure of E. edulis populations showed a significant difference in the gene pool of all populations compared. In addition, we found that 3 out of 17 populations investigated had completely distinct gene pools due to the low probability of sharing common alleles among them. This result is in agreement with Santos et al. (2015; 2016) who investigated the effect of habitat loss on the genetic structure of E. edulis, and recorded that deforestation caused a decline in the distance of gene flow of juvenile individuals. In addition, a recent study with the Oenocarpus bataua palm tree found that habitat loss and fragmentation were mainly responsible for increasing structure and reducing genetic diversity in female gametes of this species (Browne and Karubian, 2018). Thus, the limitation of pollen exchange and especially in processes involving long-distance seed dispersal should be the main reason for the strong genetic structure found in our populations, since currently there are only 12% of the original habitat for this species (Ribeiro et al., 2009). Finally, we found no evidence of spatial genetic structure on a fine scale distance. This result was already expected, because E. edulis is a predominantly allogenic, pollinated and dispersed by a large number of organisms that can move at long distances (Gaiotto et al., 2003; Dorneles et al., 2013; Galetti et al., 2013).
Conclusions and Conservation Implications
To conclude, our study revealed that the assessment of local and landscape-scale anthropogenic pressures may provide different information about the genetic vulnerability of a species. Our findings showed that intensification of logging activities affected E. edulis genetic diversity more rapidly than recent landscape modifications. However, the strong genetic structure found in all populations in our study suggests a limitation of gene flow resulting from habitat loss on a regional scale. Although current populations of E. edulis still maintain high levels of genetic diversity, this situation tends not to be maintained if local and landscape-scale anthropogenic pressures are not attenuated. This is because the loss of private alleles and increased inbreeding, as found in our results, can affect the suitability and adaptability of E. edulis populations if kept small and isolated for many generations. For this reason, to maximize conservation efforts for species of long-term economic interest, we recommend increasing inspection to reduce illegal exploitation of these species. In addition, we also suggest investments in reforestation of degraded areas in order to increase gene flow among Atlantic Forest fragments and rebalance inbreeding levels in other populations.
Data Availability Statement
The SSR genotypes are available on the Dryad Digital Repository under https://doi.org/10.5061/dryad.ghx3ffbjh.
All authors contributed to the conception of the work. Conceived the study design: LAS, EC, FG. Data collection for the manuscript: LAS, DF, LRS. Wrote the manuscript: LAS. Provided input and revised the manuscript: LAS, EC, LRS, DF, FG. All authors approved the final manuscript.
The authors thank Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (Grant N° 426828/2016-0) and Universidade Estadual de Santa Cruz (UESC) (Grant N° 00220.1100.1681) for financial support. We also thank for grants to LAS (Fapesb—Ph.D.), EC and FG (CNPq-PQ), and LRS (Capes-PNPD).
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.
We are thankful to the team of REDE SISBIOTA and the Instituto Nacional de Ciência e Tecnologia em Estudos Interdisciplinares e Transdisciplinares em Ecologia e Evolução (INCT IN-TREE) for their helpful scientific discussions, Flora Bittencourt for help with the laboratory analysis and Pavel Dodonov for help in statistical analysis.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.01090/full#supplementary-material
Araújo, M., Alger, K., Rocha, R., Mesquita, E. C. A. B. (1998). A Mata Atlântica do sul da Bahia: situação atual, ações e perspectivas. Conselho Nacional da Reserva da Biosfera da Mata Atlântica. São Paulo caderno 8.
Balkenhol, N., Pardini, R., Cornelius, C., Fernandes, F., Sommer, S. (2013). Landscape-level comparison of genetic diversity and differentiation in a small mammal inhabiting different fragmented landscapes of the Brazilian Atlantic Forest. Conserv. Genet. 14, 355–367. doi: 10.1007/s10592-013-0454-2
Benchimol, M., Talora, D. C., Mariano-Neto, E., Oliveira, T. L. S., Leal, A., Mielke, M. S. (2017). Losing our palms: The influence of landscape-scale deforestation on Arecaceae diversity in the Atlantic forest. For. Ecol. Manage. 384, 314–322. doi: 10.1016/j.foreco.2016.11.014
Calamari, N., Vilella, F. J., Sica, Y. V., Mercuri, P. A. (2018). Patch and landscape responses of bird abundance to fragmentation in agroecosystems of east-central Argentina. Avian Conserv. Ecol. 13, 3. doi: 10.5751/ACE-01222-130203
Carvalho, C. S., Ribeiro, M. C., Côrtes, M. C., Galetti, M., 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
Castro, E. R., Galetti, M., Morellato, P. C. (2007). Reproductive phenology of Euterpe edulis (Arecaceae) along a gradient in the Atlantic rainforest of Brazil. Aust. J. Bot. 55, 725–735. doi: 10.1071/BT07029
Conte, R., Reis, M. S., Reis, A., Mantovani, A., Mariot, A., Fantini, A. C., et al. (2000). “Dinâmica da regeneração natural de Euterpe edulis,” in Euterpe edulis Martius (Palmiteiro) biologia, conservação e manejo. Eds. Reis, M. S., Reis, A. (Itajaí;: Herbário Barbosa Rodrigues), 106–130.
Conte, R., Reis, M. S., Mantovani, A., Vencovsky, R. (2008). Genetic structure and mating system of Euterpe edulis Mart. populations: a comparative analysis using microsatellite and alloenzyme markers. J. Hered. 99, 476–482. doi: 10.1093/jhered/esn055
Cruse-Sanders, J. M., Hamrick, J., Ahumada, J. A. (2005). Consequences of harvesting for genetic diversity in American ginseng (Panax quinquefolius L.): a simulation study. Biodivers. Conserv. 14, 493–504. doi: 10.1007/s10531-004-7308-7
Da Silva, N. A., Rodrigues, E., Mercadante, A. Z., De Rosso, V. V. (2014). Phenolic compounds and carotenoids from four fruits native from the Brazilian Atlantic forest. J. Agric. Food Chem. 62, 5072–5084. doi: 10.1021/jf501211p 11
Dai, Z., Si, C., Zhai, D., Huang, P., Qi, S., Lin, Y., et al. (2018). Genetic effects of historical anthropogenic disturbance on a long-lived endangered tropical tree Vatica mangachapoi. J. For. Res. 29, 291–299. doi: 10.1007/s11676-017-0470-z
Dixo, M., Metzger, J. P., Morgante, J. S., 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
Dorneles, L. L., Zillikens, A., Steiner, J., Padilha, M. T. S. (2013). Biologia da polinizaça~o do açaizeiro-do-sul Euterpe edulis Mart. (Arecaceae) e associaça~o com abelhas sociais (Apidae: Apini) em sistema agroflorestal na Ilha de Santa Catarina. Rev. Iheringia 68:1, 47–57.
Ferreira, P. A., Boscolo, D., Carvalheiro, L. G., Biesmeijer, J. C., Rocha, P. L. B., Viana, B. F. (2015). Responses of bees to habitat loss in fragmented landscapes of Brazilian Atlantic Rainforest. Landscape Ecol. 30, 2067–2078. doi: 10.1007/s10980-015-0231-3
Fountain, T., Nieminen, M., Sirén, J., Wong, S. C., Lehtonen, R., Hanski, I. (2016). Predictable allele frequency changes due to habitat fragmentation in the Glanville fritillary butterfly. Proc. Natl. Acad. Sci. U. S. A. 11, 413–418. doi: 10.1073/pnas.1600951113
Gaiotto, F. A., Brondani, R. P. V., Grattapaglia, D. (2001). Microsatellite markers for heart of palm—Euterpe edulis and E. oleracea Mart. (Arecaceae). Mol. Ecol. Notes 1, 86–88. doi: 10.1046/j.1471-8278.2001.00036.x
Gaiotto, F. A., Grattapaglia, D., Vencovsky, R. (2003). Genetic structure, mating system, and long-distance gene flow in heart of palm (Euterpe edulis Mart.). Heredity 94, 399–406. doi: 10.1093/jhered/esg087
Galetti, M., Fernandez, J. C. (1998). Palm heart harvesting in the Brazilian Atlantic Forest: changes in industry structure and the illegal trade. J. Appl. Ecol. 35, 294–301. doi: 10.1046/j.1365-2664.1998.00295.x
Galetti, M., Guevara, R., Côrtes, M. C., Fadini, R., Von-Matter, S., Leite, A., et al. (2013). Functional extinction of birds drives rapid evolutionary changes in seed size. Science 340, 1086–1090. doi: 10.1126/science.1233774
González-Fernández, A., Arroyo-Rodríguez, V., Ramírez-Corona, F., Manjarrez, J., Aguilera-Hernández, A., Sunny, A. (2019). Local and landscape drivers of the number of individuals and genetic diversity of a microendemic and critically endangered salamander. Landscape Ecol. 34, 1–12. doi: 10.1007/s10980-019-00871-2
Goudet, J. (2001). FSTAT, a program to estimate and test gene diversities and fixation indices ver. 2.9.3. Available from http://www.unil.ch/izea/softwares/fstat.html. Goudet, J. (1995) FSTAT v-1.2. A computer program to calculate F-statistics. J. Hered. 86, 485–486. doi: 10.1093/oxfordjournals.jhered.a111627
Gouvêa, J. B. S. (1969). Contribuição a geomorfologia do sul da Bahia. Área dos baixos cursos dos rios Pardo e Jequitinhonha. Comun. Técnica - Cent. Pesqui. do Cacau (CEPEC)/Comissão Exec. Planej. da Lavoura Cacau. (CEPLAC) Vol. 35. Itabuna: BA, 1–11.
Gutiérrez-Granados, G. (2011). Effect of logging on rodent scatter-hoarding dynamics in tropical forests: implications for plant recruitment. Integr. Zool. 6, 74–80. doi: 10.1111/j.1749-4877.2011.00234.x
Hale, M. L., Burg, T. M., Steeves, T. E. (2012). Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PloS One 7, e45170. doi: 10.1371/journal.pone.0045170
Hansen, M. C., Potapov, P. V., Moore, R., Hancher, M., Turubanova, S. A. (2013). Tyukavina, AHigh-resolution global maps of 21st century forest cover change. al.Science 342, 850–853. doi: 10.1126/science.1244693
Jackson, N. D., Fahrig, L. (2016). Habitat amount, not habitat configuration, best predicts population genetic structure in fragmented landscapes. Landscape Ecol. 31, 951–968. doi: 10.1007/s10980-015-0313-2
Jombart, T., Collins, C. (2015). A tutorial for discriminant analysis of principal components (DAPC) using adegenet 2.0.0. London: Imperial College London, MRC Centre for Outbreak Analysis and Modelling.
Jombart, T., Devillard, S., Balloux, F. (2010). Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 11, 94. doi: 10.1186/1471-2156-11-94
Keenan, K., McGinnity, P., Cross, T. F., Crozier, W. W., Prodöhl, P. A. (2013). diveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors. Methods Ecol. Evol. 4, 782–788. doi: 10.1111/2041-210X.12067
Klank, C., Ghazoul, J., Pluess, A. R. (2012). Genetic variation and plant performance in fragmented populations of globeflowers (Trollius europaeus) within agricultural landscapes. Conserv. Genet. 13, 873–884. doi: 10.1007/s10592-012-0337-y
Lacerda, A. E. B., Kanashiro, M., Sebbenn, A. M. (2008). Effects of Reduced Impact Logging on genetic diversity and spatial genetic structure of a Hymenaea courbaril population in the Brazilian Amazon Forest. For. Ecol. Manage. 255, 1034–1043. doi: 10.1016/j.foreco.2007.10.009
Landau, E. C. (2003). “Padrões de ocupação espacial da paisagem na Mata Atlântica do sudeste da Bahia, Brasil,” in Corredor de Biodiversidade na Mata Atlântica do Sul da Bahia. Eds. Prado, P. I., Landau, E. C., Moura, R. T., Pinto, L. P. S., Fonseca, G. A. B., Alger, K. (Ilhéus: Publicação em CD-ROM, IESB/CI/CABS/UFMG/UNICAMP).
Landau, E. C., Hirsch, A., Musinsky, J. (2008). “Vegetation cover and land use in the Atlantic forest of southern Bahia, Brazil, based on satellite imagery: a comparison among municipalities,” in The Atlantic Coastal Forest of northeastern Brazil. Eds. Thomas, W. W., Britton, E. G. (New York: The New York Botanical Garden Press).
Lino, A., Fonseca, C., Rojas, D., Fischer, E., Pereira, M. J. R. (2019). A meta-analysis of the effects of habitat loss and fragmentation on genetic diversity in mammals. Mammalian Biol. 94, 69–76. doi: 10.1016/j.mambio.2018.09.006
Mantovani, A., Morellato, P. (2000). “Fenologia da floração, frutificação, mudança foliar e aspectos da biologia floral do palmiteiro,” in Euterpe edulis Martius (Palmiteiro) Biologia, Conservação e Manejo, 23–38. Eds. Reis, M. S., Reis, A. (Itajaí, Brazil: Herbário Barbosa Rodrigues).
McGarigal, K., Cushman, S. A., Neel, M. C., Ene, E. (2002) Fragstats: Spatial pattern analysis program for categorical maps. Computer software program produced by the authors at the University of Massachusetts, Amherst. Available online at http://www.umass.edu/landeco/research/fragstats/fragstats.html
Miraldo, A., Li, S., Borregaard, M. K., Flórez-Rodríguez, A., Gopalakrishnan, S., Rizvanovic, M., et al. (2016). An Anthropocene map of genetic diversity. Science 353, 1532–1535. doi: 10.1126/science.aaf4381
Montes, I., Iriondo, M., Manzano, C., Santos, M., Conklin, D., Carvalho, G. R., et al. (2016). No loss of genetic diversity in the exploited and recently collapsed population of Bay of Biscay anchovy (Engraulis encrasicolus, L.). Mar. Biol. 163, 98. doi: 10.1007/s00227-016-2866-2
Mortelliti, A., Fagiani, S., Battisti, C., Capizzi, D., Boitani, L. (2010). Independent effects of habitat loss, habitat fragmentation and structural connectivity on forest-dependent birds. Divers. Distrib. 16, 941–951. doi: 10.1111/j.1472-4642.2010.00701.x
Murawski, D. A., Gunatilleke, I. A. U. N., Bawa, K. S. (1994). The effects of selective logging on inbreeding in Shorea megistophylla (Dipterocarpaceae) from Sri Lanka. Conserv. Biol. 8, 997–1002. doi: 10.1046/j.1523-1739.1994.08040997.x
Newbold, T., Hudson, L. N., Hill, S. L. L., Contu, S., Lysenko, I., Senior, R. A., et al. (2015). Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50. doi: 10.1038/nature14324
Peakall, R., Smouse, P. F. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research— An update. Bioinformatics 28, 2537–2539. doi: 10.1111/j.1471-8286.2005.01155.x
Peres, C. A., Baider, C., Zuidema, P. A., Wadt, L. H. O., Kainer, K. A., Gomes-Silva, D. A. P., et al. (2003). Demographic Threats to the Sustainability of Brazil Nut Exploitation. Science 302, 2112–2114. doi: 10.1126/science.1091698
Peres, C. A., Palacios, E. (2007). Basin-wide effects of game harvest on vertebrate population densities in Amazonian forests: Implications for animal-mediated seed dispersal. Biotropica 39, 304–315. doi: 10.1111/j.1744-7429.2007.00272.x
Piry, S. G., Luikart, G., Cornuet, J. M. (1999). BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. J. Hered. 90, 502–503. doi: 10.1093/jhered/90.4.502
Reis, A., Kageyama, P. Y. (2000). Dispersão de sementes do palmiteiro (Euterpe edulis Martius – Palmae). In: Reis, M.S., Reis, A. (eds). Euterpe edulis Martius (Palmiteiro): biologia, conservação e manejo. Herbário Barbosa Rodrigues. Sellowia, 45–48, 60–92.
Reis, M. S., Guerra, M. P., Nodari, R. O., Ribeiro, R. J., Reis, A. (2000a). Distribuição geográfica e situação atual das populações na área de ocorrência de Euterpe edulis Martius. Sellowia 49-52, 324–335.
Reis, M. S., Fantini, A. C., Nodari, R. O., Reis, A., Guerra, M. P., Mantovani, A. (2000). Management and conservation of natural populations in Atlantic Rain Forest: the case study of palm heart (Euterpe edulis Martius). Biotropica 32, 894–902. doi: 10.1111/j.1744-7429.2000.tb00627.x
Rhoads, E. A., Williams, P. K., Krane, C. M. (2017). High inbreeding and low connectivity among Ambystoma texanum populations in fragmented Ohio forests. Ecol. Evol. 7, 11135–11147. doi: 10.1002/ece3.3637
Ribeiro, M. C., Metzger, J. P., Martensen, A. C., Ponzoni, F. J., Hirota, M. M. (2009). The Brazilian Atlantic Forest: How much is left, and how is the remaining forest distributed? Implications conservation. Biol. Conserv. 142, 1141–1153. doi: 10.1016/j.biocon.2009.02.021
Rocha-Santos, L., Benchimol, M., Mayfield, M. M., Faria, D., Pessoa, M. S., Talora, D. C., et al. (2017). Functional decay in tree community within tropical fragmented landscapes: Effects of landscape-scale forest cover. PloS One 12, e0175545. doi: 10.1371/journal.pone.0175545
Sambuichi, R. H. R. (2003). Ecologia da vegetação arbórea de cabruca – Mata Atlântica raleada utilizada para cultivo de cacau – na região sul da Bahia. [doctor’s thesis]. [Brasília (DF)]: Universidade de Brasília.
Santos, A. S., Cazetta, E., Morante-Filho, J. C., Baumgarten, J., Faria, D., Gaiotto, F. A. (2015). Lessons from a palm: genetic diversity and structure in anthropogenic landscapes from Atlantic Forest, Brazil. Conserv. Genet. 16, 1295–1302. doi: 10.1007/s10592-015-0740-2
Santos, A. S., Cazetta, E., Dodonov, P., Faria, D., Gaiotto, F. A. (2016). Landscape-scale deforestation decreases gene flow distance of a keystone tropical palm, Euterpe edulis Mart (Arecaceae). Ecol. Evol. 6, 6586–6598. doi: 10.1002/ece3.2341
Santos, J., Varassin, I. G., Muschner, V. C., Ovaskainen, O. (2018b). Estimating seed and pollen dispersal kernels from genetic data demonstrates a high pollen dispersal capacity for an endangered palm species. Am. J. Bot. 105, 1802–1812. doi: 10.1002/ajb2.1176
Schlaepfer, D. R., Braschler, B., Rusterholz, H. P., Baur, B. (2018). Genetic effects of anthropogenic habitat fragmentation on remnant animal and plant populations: A meta-analysis. Ecosphere 9, e02488. doi: 10.1002/ecs2.2488
Sebbenn, A. M., Degen, B., Azevedo, V. C. R., Silva, M. B., Lacerda, A. E. B., Ciampi, A. Y., et al. (2008). Modelling the long-term impacts of selective logging on genetic diversity and demographic structure of four tropical tree species in the Amazon forest. For. Ecol. Manage. 254, 335–349. doi: 10.1016/j.foreco.2007.08.009
Shaanker, R. U., Ganeshaiah, K. N., Rao, M. N., Ravikanth, G. (2003). “Genetic Diversity of NTFP Species: Issues and Implications,” in Proceedings of International Workshop on Policies, Management, Utilization and Conservation of Non-Timber Forest Products in the South Asia Region. Eds. Hiremath, A. J., Joseph, G. C., Shaanker, R. U. (Bangalore: ATREE and Forest Research Support Programme for Asia and the Pacific, FAO), 40–44.
Shaanker, R. U., Ganeshaiah, K. N., Rao, M. N., Aravind, N. A. (2004). Ecological Consequences of Forest Use: From Genes to Ecosystem - A Case Study in the Biligiri Rangaswamy Temple Wildlife Sanctuary, South India. 347–363: Conservat. Soc. 2:2.
Silva, P., Carmo, L., Silva, G., Silveira-Diniz, M. F., Casemiro, R. C., Spoto, M. H. F. (2013). Physical, chemical, and lipid composition of juçara (Euterpe edulis Mart.) Pulp. Braz. J. Food. Nutr. 24, 7–13.
Soares, L. A. S. S., Faria, D., Vélez-Garcia, F., Vieira, E. M., Talora, D. C., Cazetta, E. (2015). Implications of Habitat Loss on Seed Predation and Early Recruitment of a Keystone Palm in Anthropogenic Landscapes in the Brazilian Atlantic Rainforest. PloS One 10, e0133540. doi: 10.1371/journal.pone.0133540.
Thomas, W. W., Carvalho, A. M., Amorim, A. M. A., Garrison, J., Arbeláez, A. L. (1998). Plant endemism in two forests in southern Bahia, Brasil. Biodivers. Conserv. 7, 311–322. doi: 10.1023/A:1008825627656
van Oosterhout, C., Hutchinson, W. F., Wills, D. P., 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
van Oosterhout, C., Weetman, D., Hutchison, W. F. (2006). Estimation and adjustmentof microsatellite null alleles in non-equilibrium populations. Mol. Ecol. Notes 6, 255–256. doi: 10.1111/j.1471-8286.2005.01082.x
Vranckx, G., Jacquemyn, H., Muys, B., Honnay, O. (2012). Meta-analysis of susceptibility of woody plants to loss of genetic diversity through habitat fragmentation. Conserv. Biol. 26, 228–237. doi: 10.1111/j.1523-1739.2011.01778.x
Wood, D. A., Bui, T. D., Overton, C. T., Vandergast, A. G., Casazza, M. L., Hull, J. M., et al. (2017). A century of landscape disturbance and urbanization of the San Francisco Bay region affects the present-day genetic diversity of the California Ridgway’s rail (Rallus obsoletus obsoletus). Conserv. Genet. 18, 131–146. doi: 10.1007/s10592-016-0888-4
Zhang, X., Shi, M., Shen, D., Chen, X. (2012). Habitat Loss other than Fragmentation per se Decreased Nuclear and Chloroplast Genetic Diversity in a Monoecious Tree. PloS One 7, e39146. doi: 10.1371/journal.pone.0039146
Zucchi Maria, I., Sujii, P. S., Mori, G. M., Viana, J. P. G., Grando, C., Silvestre, E. A., et al. (2018). Genetic diversity of reintroduced tree populations in restoration plantations of the Brazilian Atlantic Forest. Restoration Ecol. 26, 694–701. doi: 10.1111/rec.12620
Keywords: landscape genetics, tropical rainforest, conservation, threatened species, molecular marker
Citation: Soares LASS, Cazetta E, Santos LR, França DdS and Gaiotto FA (2019) Anthropogenic Disturbances Eroding the Genetic Diversity of a Threatened Palm Tree: A Multiscale Approach. Front. Genet. 10:1090. doi: 10.3389/fgene.2019.01090
Received: 01 July 2019; Accepted: 09 October 2019;
Published: 07 November 2019.
Edited by:Rosane Garcia Collevatti, Universidade Federal de Goiás, Brazil
Reviewed by:Alexandre Magno Sebbenn, Instituto Florestal, Brazil
Kaiser Dias Schwarcz, Instituto Federal de Educação, Ciência e Tecnologia de Brasília (IFB), Brazil
Juliana Silveira Dos Santos, Universidade Federal de Goiás, Brazil
Andrzej Oleksa, Kazimierz Wielki University of Bydgoszcz, Poland
Copyright © 2019 Soares, Cazetta, Santos, França and Gaiotto. 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: Fernanda Amato Gaiotto, email@example.com