Herbivory and leaf traits of two tree species from different successional stages in a tropical dry forest

ABSTARCT

Leaf traits are analyzed as essential drivers for the evolutionary and ecological role of plant defense mechanisms. Plants show leaf trait variation as a response to the diverse environmental conditions, like different successional stages. Those variations can impact leaf herbivory rates and drive changes in the allocation of plant resources. In this study, we aimed at comparing the expression of leaf defenses in established plants over different successional stages in a tropical dry forest to understand how these defenses modify the plant–herbivore interactions based on herbivory rates. We analyzed physical leaf traits (specific leaf area, thickness, and density), nutrient content (N, P, and K), total phenolic compounds, and leaf herbivory, of the native tree species Aspidosperma pyrifolium and Cenostigma pyramidale, in early- and late-successional stage areas. Results showed that the plant investment in defenses varies according to the successional stage and that both species have similar defense strategies, confirming the resource availability hypothesis. Individuals from the early stage adopt a strategy to lower sclerophylly, higher nutrient content and less phenolic compounds, while the late-stage individuals showed an opposite behaviour. For both species in this study, the average percentage of leaf herbivory observed was 40% higher in the early compared to the late-stage area. Our data indicate that plant defenses are tightly coupled to sclerophylly and investments in secondary metabolites, and the environmental conditions of different successional stage drive that plasticity in such leaf traits.

KEYWORDS:

• Leaf nutritional content
• natural regeneration
• phenolic compounds
• sclerophylly
• semiarid

Introduction

Plant–insect interactions are essential for maintaining the biodiversity and the function in ecosystems [1]. The production of phenolic compounds and a higher leaf sclerophylly are examples of the many traits that plants exhibit as defenses, aimed to reduce the loss of photosynthetic tissue due to possible herbivore attacks [2,3]. These characteristics are mainly influenced by environmental conditions, such as nutrient and soil water availability [4–6].

Many hypotheses that describe the plant resource allocation and defense patterns have been postulated, with the “Resource Availability” being one of the most well known and mediated by the leaf economics spectrum [7–9]. In pioneer plant species, plants invest more resources for growth, despite defensive leaf traits, making the replacement of tissue lost to herbivory more beneficial [4,7]. On the other hand, plant species in late successional stages often adopt a conservative (or slow maturation) growth approach by allocating more resources to the formation of tissues that have a higher longevity and providing additional protection against herbivores [7]. Distinct ecological strategies are considered positive in the same successional stage, but usually the establishment of a plant species in different succession stages is limited by trade-offs [10–12].

In the tropics, huge forest areas are used for farming and cattle raising, while other already-deforested areas are abandoned after usage, creating patches of forests under different successional stages [13]. This anthropic disturbance has been widely observed in tropical dry forests, especially in South America [14]. In Brazil, 80% of dry forests are in regeneration and 40% of these are maintained in the initial stages by constant land use [15]. In areas undergoing natural regeneration, the plasticity displayed on leaf traits is an important strategy for plant acclimatization and establishment [5]. These traits can be understood as morpho-physiophenological characteristics arising from the plants response to the environmental conditions [2,16,17]. Variation in leaf traits has been shown to be involved in countless plant processes, including plant growth, and can also impact the herbivory rates observed along the succession stages [2,7]. Furthermore, through morphological and physiological adjustments, herbivore-inflicted pressure can also be considered an important factor for the development and evolution of these traits, including physical and chemical defenses, and tolerance of future herbivore attacks [17,18].

Among the leaf traits that have been widely used in studies about ecological attributes and plant defense we find at least six distinctive examples: a) leaf area (LA) which represents the balance between carbon assimilation and water loss through transpiration; b) specific leaf area (SLA), a trait that reflects the leaf building cost – especially in carbon – and is directly related to plant function [18]; c) thickness and density, scleromorphic characteristics that have important relationships the leaf mechanical properties and are related to the leaf life span and plant–herbivore interactions [19–21]; d) nutrient composition, such as nitrogen and phosphorus; e) secondary compounds – essential to protect the leaf epidermis against high light irradiance and herbivore attack [22]. Thus, it is expected that the ecological strategies adopted by plant species under different circumstances can be inferred based on these leaf traits [23,24].

Studies linking the effects of herbivory-inflicted leaf damage with the different traits of the plants that manage to persist in different succession stages are still scarce. This is even more evident in the case of tropical dry forests. As herbivory can have a direct influence on forest ecosystems’ function and structure, it is crucial to learn how these interactions act in a landscape scale, as this information is useful for restoration and conservation strategies [25]. Thus, the current study aims to compare the expression of leaf defenses of two tree species established over different successional stages and to investigate how these defenses modify the plant–herbivore interactions based on herbivory rates. We hypothesize that the leaf damage inflicted by herbivore insects in a tropical dry forest would be higher in early-successional plants due to a higher investment in growth than in the development of leaf defenses.

Material and methods

Study area

This study was conducted at Fazenda Tamanduá (06°59ʹ13” to 07°00ʹ14” S e 37°18ʹ08” to 37°20ʹ38” W, 240 m above sea level), a farm in the countryside of Santa Terezinha, Paraíba, Brazil (Figure 1). With 3,073 ha, an average annual temperature of 32.8°C and an average annual precipitation of 600 mm. The rainfall episodes show an irregular pattern, which usually is condensed between February and May. During 2016 and 2017, the accumulated rainfall until April was, respectively, 590.3 mm and 531.6 mm (Figure 1). The Lithic Neosol (i.e. soil class of little developed and shallow, reaching up to 50 cm in depth, being normally associated with rocky outcrops) is the predominant type in this area [26].

Figure 1. Map of Brazil locating Paraiba state and the study area map (Fazenda Tamanduá). In details a google map of Fazenda Tamandua with the successional stages sites points. Historical series of annual rainfall (mm), accumulated precipitation from January to April (mm) and annual average air temperature (°C) over the last 10 years in the tropical dry forest (Santa Terezinha, Paraíba, Brazil)

The farm has one of the biggest Caatinga conservation and protection areas in the Brazilian Northeast, with almost 1,000 ha of land (one-third of its total area) comprising a Legal Reserve with 614 ha, and a Private Natural Reserve – RPPN Tamanduá (350 ha) recognized by IBAMA, the country’s Federal body in charge of environment and natural resources [27].

Nowadays, the land area is the property of Fazenda Tamanduá; therefore, it was originally composed by many small properties that cultivated perennial cotton plants and were abandoned during a widespread infestation by the cotton boll weevil in the 1980s [27,28]. The abandoned areas comprise a mosaic of different successional stages of tropical dry forests. The land-use history is well known, and there are no records of fertilizer usage in these places [5,6,28]. Our research group established five permanent plots since 2006 in each different stages of ecological succession (early, intermediate and late) when all the trees (DBH > 3 cm) were identified and got an aluminum labels to mark each tree with one code (number and letter). The early-successional stage has a history of anthropic disturbance such as clearcutting, cotton cultivation, and cattle pasture, undergoing natural regeneration for approximately the last 20 years. It also displays a high prevalence species of herbs during the rainy season alongside species of shrubs and trees with a relatively sparse distribution. The richness of trees in the early area is represented by seven species and three families, among them Aspidosperma pyrifolium Art. & Zucc, Cenostigma pyramidale Tull. and Bauhinia cheilantha (Bong.) Steud. The mean shoot biomass production for the early succession stage area was 29.9 Mg ha⁻¹ [15,28]. On the other hand, the late succession stage area (Figure 2(a)) is remarkably free of clearcutting and signals of anthropic action or big disturbance since the 1950s. This area is inside a Private Natural Reserve, and the vegetation is highly diverse when compared to the early successional stage, displaying many late-successional species. The richness is represented by 21 tree species distributed in 12 families, among them Anadenanthera colubrina (Vell.) Brenan, A. pyrifolium, C. pyramidale and Commiphora leptophloeos (Mart.) J.B. Gillett. The mean shoot biomass production for the late-successional stage area was 49.4 Mg ha⁻¹ [15,28].

Figure 2. General view of late-successional stage of Seasonally dry tropical forest (Caatinga) at rainy season (a); Aspidosperma pyrifolium (Apocynaceae) and a branch with leaves and green fruits in details (b) and Cenostigma pyramidale (Fabaceae) and the leaves in details (c)

Collection and sampling of the target tree species

A concomitant presence in both succession stages was the main criterion for the selection of model species and was based on information by Cabral and collaborators [15]. Two pioneer tree species were selected in this study: Aspidosperma pyrifolium (Apocynaceae) (Figure 2(b)), a native species with a massive presence in Brazilian tropical dry forests, such as Caatinga and Cerrado. These trees are deciduous and can grow up to 8 m, displaying straight trunks with smooth bark and alternate simple leaves [29,30] and Cenostigma pyramidale (Fabaceae) (Figure 2(c)), a semi-deciduous species endemic to the Caatinga. Individuals grow up from 4 to 6 m and display bipinnate leaves, light grey trunk, and high wood density [31]. The sampling criteria were defined as trees with a diameter at breast height (DBH) > 7 cm, and at least 10 individuals/species in each collection area.

In order to obtain an estimation of the accumulated effect (i.e. damage) caused by herbivores during the leaf lifespan, the sampling was performed at the end of the rainy season, in April 2016 and 2017. We collected and measured leaves from the two species located in the early- and late- succession stage areas. Each collection area has 2,000 m² and is fenced with barbed wire to avoid the presence of livestock.

Leaf herbivory

The herbivory percentage was obtained from a set of 20 leaves of each tree for both species. Sample collections were performed in two consecutive years, always at the end of the rainy season, guaranteeing the estimate of the damage accumulated during leaf’s lifespan. Leaves were scanned for damage measurement. We considered damage as the removal of the leaf blade caused by folivores (chewers and cutters leaves). We used ImageJ [32] to quantify the total leaf area for each leaf, as well as the areas removed by the folivores. Each image was converted to a binary file (black and white) from which the herbivory percentage was obtained using the following formula: [(area lost to herbivory/total leaf area) * 100]. The herbivory percentage was subsequently analyzed individually for each tree.

Leaf traits

Ten leaves from trees that satisfied the sampling criteria for each species in both successional stages were collected in two consecutive years and used for leaf trait measurements. We used the same trees for measurements every year. Leaf thickness was obtained using a digital calliper at leaf’s widest section while avoiding the midrib, as described by Pérez-Harguindeguy et al. [33]. Specific leaf area (SLA) was calculated using leaf area divided by leaf dry matter. The leaves were digitalized individually using a scanner, and each leaf area was obtained using ImageJ [32]. Subsequently, the same leaves were oven-dried individually at 70°C for 72 hours and the respective weights were obtained using a precision scale (HR-200, E, Tokyo, Japan) [34]. Leaf density was determined by the ratio between specific leaf mass and thickness [35].

The nitrogen (N), phosphorus (P), and potassium (K) content in leaves were obtained using 250 mg of oven-dry and ground leaf mass. This dry matter underwent acid digestion (H₂SO₄) in a digestion block at 350°C for extraction. Total N content was quantified from titration of the digested plant extract with HCl after the addition of boric acid and a coloured indicator [36]. Spectrophotometry was used for the analysis of P content (Spectrophotometer 600S, FEMTO, São Paulo, BR) following the methodology described by Murphy and Riley [37]. The N:P ratio was also calculated. The content of K was derived by flame photometry (DM-62, Digimed, São Paulo, BR) using a standard solution of 5 ppm of K [38].

We used 20 mg of dry leaf material for the determination of total phenolic compounds, boiling it in methanol 80% for extraction and further submitting it to a reaction with Folin-Ciocalteu 10% as described by Amorim et al. [39]. Results were derived by interpolating samples’ absorbance against a standard curve, which was obtained projecting tannic acid (0.1 mg ml⁻¹) in a spectrophotometer (Genesys 10S UVVIS, Thermo Scientific, Waltham, EUA) adjusted to a wavelength of 760 nm. The results from all measured traits were averaged individually for each sampled tree from both species.

Data analysis

The leaf trait variations in this study (specific leaf area, thickness, density, nutrient content, phenolic compounds, and herbivory rates) were analyzed using repeated measures MANOVA, with the succession stage set as an independent variable (Table 1). Averages were further compared using the Bonferroni test (5%). The analyses were carried out using R 2.15 [40]. A principal component analysis (PCA) was performed for each tree species to identify which leaf traits were more responsive to successional stage, using software Fitopac 2.1 [41].

Table 1. Multivariate analysis of variance (MANOVA) of functional leaf traits: specific leaf area (SLA), leaf density (DEN), leaf thickness (THICK), total phenolic compounds (PC), nutrient content (N, P, K, N:P), and leaf tissue lost to herbivory (HERB) in tropical dry forest (Santa Terezinha, Paraíba, Brazil)

Species	Aspidosperma pyrifolium				Cenostigma pyramidale
Response variables	F	df	P		F	df	P
SLA	73.034	18	0.000		11.02	18	0.004
DENS	49.73	18	0.000		52.77	18	0.000
THICK	124.764	18	0.000		162.16	18	0.000
N	36.292	18	0.000		286.64	18	0.000
P	110.077	18	0.000		29.54	18	0.000
K	321.59	18	0.000		28.63	18	0.000
N:P	339.232	18	0.000		165.22	18	0.000
PC	9.851	18	0.006		146.37	18	0.000
HERB	73.034	18	0.000		173.15	18	0.000

Significance of bold values: P>0.05

Results

The average percentage of leaf loss due to herbivory in A. pyrifolium was 32% higher in the early-successional area compared to the late stage (P < 0.05, Figure 3(a)). C. pyramidale trees also showed higher herbivory rates in the early stage (P < 0.05), with an increase of about 42% with respect to the values observed in the late stage (Figure 3(b)). The highest percentage of herbivory was exhibited by A. pyrifolium.

Figure 3. Leaf tissue lost to herbivory for Aspidosperma pyrifolium (a) and Cenostigma pyramidale (b) individuals from early and late succession stages of tropical dry forest over two consecutive years (Santa Terezinha, Paraíba, Brazil). Bars ± standard error followed by similar letters indicate a non-significant difference according to the Bonferroni test (5%), n = 10

Leaves from A. pyrifolium and C. pyramidale individuals in the early succession stage had lower levels of carbon-based defenses, as observed through sclerophylly results (i.e. contrary to SLA). SLA was significantly higher in the early stage for both species and in the two sampled years. This difference between the stages represented a minimum increase of 25% in the SLA in the early area (P < 0.05, Figure 4(a,b)). The highest value of density between species was observed in C. pyramidale in 2016 and 2017. However, A. pyrifolium and C. pyramidale trees from the late succession stage showed higher leaf density values (in general, this increase was over 80%) in the two sampled years (P < 0.05; Figure 4(c,d)). Leaf thickness showed the same trend, with values up to 42% higher for individuals from both species in the late succession stage (P < 0.05; Figure 4(e,f)). In turn, A. pyrifolium leaves had higher values of thickness. The species A. pyrifolium and C. pyramidale showed a similar behavior concerning the concentration of phenolic compounds, indicating a significant increase from early to late succession stages in the two sampling years (Figure 4(g,h)). Total phenolic compounds in A. pyrifolium and C. pyramidale had a difference of approximately 1.4 and 1.7 times between the areas, respectively.

Figure 4. Leaf traits for Aspidosperma pyrifolium and Cenostigma pyramidale individuals from early and late succession stages of tropical dry forest over two consecutive years (Santa Terezinha, Paraíba, Brazil). (a; b) specific leaf area; (c; d) leaf density; (e; f) leaf thickness; (g; h) phenolic compounds. Bars ± standard error followed by similar letters indicate a non-significant difference according to the Bonferroni test (5%), n = 10

Leaf N and P contents decreased from the early to the late succession stage for both species (P < 0.05, Table 2). This was observed in both years of sampling. Leaf K content was higher for A. pyrifolium individuals from the early succession stage (P < 0.05, Table 2). On the other hand, C. pyramidale trees did not show any differences in K content (P > 0.05, Table 2). We observed an increase in the N:P ratio for A. pyrifolium and C. pyramidale showed for both succession stages (P < 0.05, Table 2).

Table 2. Leaf nitrogen, phosphorus and potassium content, and N:P ratio for Aspidosperma pyrifolium and Cenostigma pyramidale individuals from early and late succession stages of tropical dry forest over two consecutive years (Santa Terezinha, Paraíba, Brazil). Averages ± standard error followed followed by similar letters indicate a non-significant difference according to the Bonferroni test (5%), n = 10

Sucessional Stage	Year	N (g.Kg¯^1)	P (g.Kg¯^1)	K (g.Kg¯^1)	N:P (ratio)
Aspidosperma pyrifolium
Early	2016	22.4 ± 0.5 c	3.0 ± 0.1 a	2.0 ± 0.1 a	7.3 ± 0.1 c
	2017	29.2 ± 0.1 a	3.0 ± 0.1 a	1.9 ± 0.2 a	9.8 ± 0.2 b
Late	2016	13.9 ± 0.3 d	2.7 ± 0.1 b	1.7 ± 0.1 b	5.1 ± 0.2 d
	2017	24.0 ± 0.2 b	2.2 ± 0.1 c	1.6 ± 0.1 b	10.9 ± 0.1 a
Cenostigma pyramidale
Early	2016	26.6 ± 0.4 b	3.9 ± 0.1 a	1.7 ± 0.1 a	6.8 ± 0.1 a
	2017	29.2 ± 0.2 a	3.8 ± 0.1 a	1.9 ± 0.1 b	7.6 ± 0.1 b
Late	2016	13.0 ± 0.2 c	3.5 ± 0.1 b	1.8 ± 0,1 ab	3.7 ± 0.1 c
	2017	25.8 ± 0.1 b	3.4 ± 0.1 b	1.9 ± 0.1 b	7.6 ± 0.1 bd

The PCA shows that year of collection and successional stage are the main drivers of the identified groups. In A. pyrifolium, the analysis showed a greater relationship between the successional stage to axis 1 (50.43%), while axis 2 was more related to the annual variation of attributes (32.91%). Leaf traits such as the total phenolic content, thickness, and density, with higher values in the late area, explained more strongly the differences presented between the areas in axis 1. In parallel, axis 2 was explained by the increase in the mobilization of nutrients, especially N and SLA (Figure 5(a)). For C. pyramidale, the successional stages also formed evident groups and more related to axis 1 of the PCA (63.68%), especially due to characteristics such as density, thickness and phenolic content. The variation in axis 2 was explained by 21.90% and also due to these interannual differences (Figure 5(b)).

Figure 5. Principal component analysis of the data collected from Aspidosperma pyrifolium (a) and Cenostigma pyramidale (b) in the tropical dry forest (Santa Terezinha, Paraíba, Brazil). Early stage (circle), late stage (square). The size of the symbols represents leaf tissue lost to herbivory. Specific leaf area (SLA), leaf density (DEN), leaf thickness (THICK), total phenolic compounds (PC) and nutrient content (N, P, K, N:P)

Discussion

The percentage of plant tissue lost to herbivory found for the species in this study is within the average values observed for tropical dry forests [42]. For instance, Fonseca et al. [2] found leaf herbivory rates between 2.5% and 4.4% in a tropical dry forest in Minas Gerais, Brazil. Aquino et al. [43] found leaf damage rates between 8.8% and 12% for Aspidosperma pyrifolium. Variations in herbivory percentages between species reflect the use of different strategies, such as physical and chemical defenses by plants, their characteristics, present environmental conditions, or even its tolerance to herbivory [18,44,45]. In our study, we observed a higher percentage of herbivory in A. pyrifolium. In general, the species showed a lower leaf density, which may have contributed to this difference. In addition, deciduous species have less investment in other defenses types such as chemical compounds that demand higher metabolic costs (i.e. tannins and alkaloids), generating herbivory percentages up to 2.8 times higher than those reported in perennial species for dry forests [46]. Our results corroborate this premise. Dourado et al. [47] studying herbivory in Caatinga also reported more than twice the leaf loss in deciduous species and with less leaf hardness, compared to perennial species. The herbivory increment between the analyzed years for the two species can also be justified by the increase in the leaf traits, such as SLA and nitrogen content and the decrease in phenolic compounds. All of these are associated with the plant’s palatability [18,33,43]

In general, the distinct responses shown by the measured leaf traits in the different stages (early and late) for both tree species corroborate predictions based on the resource availability hypothesis [4,7]. This hypothesis postulates that there is a higher investment in growth than in the production of defenses, as the costs of building leaves and replacing the tissue lost to herbivores are lower due to nutrient quantities and high availability of light in the early succession stages. Similar responses are also documented in other studies in tropical dry forests [2,48,49]. The area lost to herbivores in leaves for both species from the late stage showed a decrease of over 40% when compared to individuals from the early stage (Figure 3).

SLA increments in the early stages for both species are directly linked to plants photosynthetic capacity and efficiency in resource usage – like nutrients and light – lowering the overall cost of replacing tissues lost to herbivores [3,33,50]. These leaves are more vulnerable to herbivory, as they usually are tenderer and less resistant [48]. Increases in thickness and density in the late stage may have played a role in reducing herbivory’s incidence [51], but also elevate leaves’ production costs and often their lifespan time [22,44,52]. The leaf density is also considered a type of initial defense against many kinds of generalist herbivores, rendering harder leaves, with higher biomechanical resistance and reducing the palatability and/or tissue nutrition content [3,44,53]. In the early stage, an increase in nutrient content and a significant decrease in the content of phenolic compounds favor higher palatability and represent a higher metabolic activity when compared to leaves with opposing traits [3,10,53].

Individuals from the early stage adopt a strategy of low-cost investment (lower sclerophylly, higher nutrient content and less phenolic compounds), while the late-successional stage individuals allocate resources in long-lived tissue (higher sclerophylly, lower nutrient content) with higher levels of phenolic compounds, that are frequently related to higher protection against herbivores [10,44]. These results for both species suggest the same trend in leaf trait changes, possibly influenced by environmental differences between stages, for example, the light interception, the amount of soil nutrients, and water availability. Leaf traits can vary due to numerous environmental pressures, such as changes in SLA due to water conservation, but these adaptations can also favor the plant’s defenses against herbivores [44,54,55]. Furthermore, while leaf traits are prevalent over the inflicted herbivory percentage [56], higher levels of herbivory pressure could also be playing an essential role to elevate the herbivory rates observed in the early-successional stages [57]. In the same area of the current study, Vieira [58] found a higher abundance of arthropods in the early stage during wet and dry seasons when compared to the late-successional stage, and Silva et al. [59] found higher levels of richness and abundance of herbivore insects in early succession stage in tropical dry forests.

The tendency of producing simpler and lower cost tissue is corroborated by the higher concentrations of leaf N and P in trees from the early stage [4,7]. The increment in these nutrients also indicates a nitrogen mobilization for protein synthesis, in particular for proteins linked to photosynthesis as Rubisco, with the energy used to higher production of leaf biomass [6,60,61]. This perspective is correct if we analyze the increase in SLA together with the increase in N content during the years of sampling in this study. On the other hand, the lower content of N and P found in the plant species from the late succession stage might indicate a shift of energy use to metabolize phenolic compounds through a limitation in protein synthesis [61,62]. In fact, the increase observed in N content for species in the late stage in 2017 was followed by a decrease in the phenolic compounds, especially in C. pyramidale. The N:P ratio for both tree species also corroborates these data. These features reinforce the behavior of producing leaves that are costly and have more defense compounds, and they were really efficient, given the lower values of herbivory found in the late stage.

It is well known that the content of secondary metabolite, nutrients, and physical traits are central to how the herbivore perceives the quality of the plant [8,63,64]. For example, the negative relationship between phenolic compounds content and leaf herbivory has been pointed out by many authors, including in tropical dry forests, and indicates the important role these substances play on providing the plant with effective defensive mechanisms [7,18,53,65]. The higher content of N in the leaf seems to be a fundamental aspect for leaf consumption. This relationship has been observed in many ecosystems, including tropical dry forests [58,66–68]. The Principal Component Analysis showed the formation of groups in relation to successional stages; thus, the herbivory rate indicates to be higher in environments where the investment in defensive traits was lower and the nutritional quality of the leaves was higher [2,44].

Conclusions

Our results show that A. pyrifolium and C. pyramidale showed different strategies related to herbivore defense depending on successional stage. On the one hand, there is the production of low-cost leaves in the early stage that resulted in higher percentages of herbivory. On the other hand, there is allocation of resources in tissues with a higher life span, in addition to an increase in defense metabolites, which require more energy from the plant, but result in lower percentages of herbivory. In summary, we observed that plant defenses are strongly linked to the leaf economic spectrum. That is, chemical (nutrient and phenolic content) and structural (SLA, density, and thickness) properties in leaves determine how plants may adopt a growth-pattern that runs from quick to slow return on investments of nutrients and dry mass in leaves.

The understanding of how plants allocate resources according to the needs imposed by herbivory, and also by the abiotic environment, can contribute to determine which species are a priority for assisted restoration programs in degraded areas.

Acknowledgments

We gratefully acknowledge for Instituto Fazenda Tamanduá for logistical support and Mr Pierre Landolt for allowing us to develop research at Fazenda Tamanduá. Claudio Magalhães to manuscript English first review. Seu Jorge to field assistant. The authors thank FACEPE and Capes (financial code 001) for the research grant. The authors acknowledge the Brazilian Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq –PQ 309965/2016-0), for scientific productivity fellowship. We would also like to thank the two anonymous reviewers for commenting on manuscript.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was supported by the Fundação de Apoio a Pesquisa de Pernambuco (FACEPE); Conselho Nacional de Desenvolvimento Científico e Tecnológico (BR) [PQ 309965/2016-0]; Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).