Skip to main content
  • Research article
  • Open access
  • Published:

The composition of braconid wasp communities in three forest fragments in a tropical lowland forest of Panama

Abstract

Background

In the last 171 years, the forests along the eastern bank of the Panama Canal have been pressured by anthropic activities. Studies of the influence of habitat fragmentation on braconid wasp communities in Central America is scarce, showing the existing information gap on these communities required to implement strategic plans for ecosystem sustainability and conservation. This study investigated how fragmentation affects braconid wasp communities in three areas in Panama City: Metropolitan Natural Park, Albrook and Corozal. Two permanent Malaise Traps were installed in the center of each fragment and were reviewed weekly from May 2019 to March 2020. Alpha and beta diversity indices and the similarity index were used to demonstrate the composition of braconid wasp communities in three forest fragments.

Results

A similarity of 94% was estimated for the subfamily composition and 74% was estimated for the morphospecies composition of wasp community in the fragments studied. Wasp subfamily and morphospecies assemblages were more similar between fragments of Albrook and Metropolitan Natural Park. Richness and abundance of braconid wasps observed were statistically different between the fragments studied.

Conclusion

Richness, abundance, and composition of braconid wasps differ among habitat fragments with high similarity between subfamilies and morphospecies. Therefore, the fragments studied can be used as stepping stones to maintain remaining populations of braconid wasp communities. Monitoring is recommended to assess the effect of fragmentation on the remaining forests.

Background

Habitat fragmentation represents one of the most serious threats to global biodiversity [1,2,3,4,5,6]. Global assessments have shown that habitat fragmentation is provoking a decrease in population size and an increased risk of disappearance of many flora and fauna [7, 8]. More than 50% of tropical forests have been degraded globally and with the use of satellite images, it is revealed that almost 43% of the terrestrial surface has been converted from its natural state for anthropogenic purposes [9, 10]. From 2000 to 2010, there was a net loss of forest cover of 7 million hectares (ha) per year in the tropical countries of the world [11]. In Panama, the forest cover has decreased from 5,245,000 ha in 1947 to 2,481,658 ha in 2019, which represents 47% forest cover loss [12]. If the fragmentation process continues at an exponential rate, the world’s tropical forests could disappear completely [13].

Habitat fragmentation often occurs due to some disturbance mechanism (agriculture, deforestation, urbanization, fires, etc.) as a result of topographic differences [14]. As the world’s human population continually increases, urban areas are growing rapidly and threaten the habitat of native species of flora and fauna. Studies investigating the influence of habitat fragmentation on biodiversity and the risk of species extinction are of the utmost importance and have been a main focus of biodiversity conservation research [15,16,17,18,19,20]. For instance, the abundance of dung beetle decreased with increased urban land cover [21]. Both common and rare species of social wasps are threatened by forest fragmentation in Central Amazon [22].

In tropical forests, evidence indicates that one of the taxa that respond faster to environmental changes are insects [23, 24]. Insects play important roles in almost all trophic levels; therefore, it is important to understand the response of these organisms to fragmentation. Within this taxon, parasitoids as a group may be used to assess the effects of habitat fragmentation because they play an ecological role in regulating populations of other insects due to prey denso-dependence [13]. González and Ruíz [25] proposed the use of braconid wasps (parasitoids) as indicators of biological diversity in deciduous forests and in evaluating and monitoring the effects of anthropogenic activities on ecosystems. Braconid wasps are regulatory agents of various groups of herbivorous insects that indicate the presence or absence of other species through the food chain [26]. Most braconid wasps are endo and exo parasitoids which feed on the larval stages of Coleoptera, Diptera and Lepidoptera [27]. This makes braconid wasps good biological indicators of habitat disturbances [25].

Smith and Mayfield [28] studied the taxonomic and functional diversity of bees visiting flowers of three tree species in small and large tropical forest fragments in tropical Australian landscapes. Species and functional diversity were found to differ significantly between small and large fragments. There was less taxonomic diversity of bees visiting flowers in small fragments. Additionally, native eusocial stingless bees were not common on small remains despite the presence of floral resources similar to those sampled on large remains [28].

In a similar study, Ruiz-Guerra et al. [29] studied the abundance, species richness, similarity, and prevalence of braconid parasitoid wasps for four types of land use (secondary forest, plantations, live fences, and pastures) and preserved tropical humid forest remnants in southern Mexico. Species richness and abundance were found to be higher in preserved and secondary forests than in other land use types.

To establish more direct links with ecosystem processes, there is a need to investigate patterns of functional trait and taxonomic diversity with biological indicators. Fragmentation experiments are useful tools used to provide clear evidence of the strong and typically degrading impacts of biodiversity loss [30]. Studies on the influence of habitat fragmentation on braconid wasp communities in lowland forests of Central America are scarce, portraying the existing information gap on these communities required to implement strategic plans for ecosystem sustainability and conservation. It is expected that landscapes with an intermediate degree of fragmentation will cause separation of braconid wasp communities which is reflected by a low similarity between communities. This research study contributes to the construction of baseline data by evaluating the vulnerability of braconid wasp communities leading to strategic plans for the sustainability and conservation of ecosystems. We sought to determine how habitat fragmentation may affect braconid wasp communities in fragmented lowland forest locations in Panama using similarity and fragmentation indices.

Results

A total of 1697 individual wasps belonging to 77 morphospecies and 16 subfamilies were recorded. Of the 1697 individuals, approximately 39% were collected in the PNM fragment, 36% in the fragment of COR and 25% in the fragment of ALB (Table 1). Among the three fragments, Rogadinae was found to be the most abundant subfamily with 456 individual wasps, followed by Alysiinae with 391 individual wasps, Adeliinae with 254 individual wasps, Doryctinae with 168 individual wasps, Aphidiinae with 154 individual wasps and Microgastrinae with 108 individual wasps. In the PNM fragment, 664 individual wasps and 14 subfamilies were observed (Table 1); with the most abundant subfamily being Adeliinae with 119 individual wasps, followed by Doryctinae with 114 individual wasps, Rogadinae with 106 individual wasps, Alysiinae with 89 individual wasps and Microgastrinae with 85 individual wasps (Fig. 1). In the COR fragment, 603 individual wasps and 13 subfamilies were observed (Table 1); with the most abundant subfamily being Alysiinae with 268 individual wasps, followed by Rogadinae with 160 individual wasps, Aphidiinae with 53 individual wasps and Adeliinae with 48 individual wasps (Fig. 1). In the ALB fragment, 430 individual wasps and 16 subfamilies were observed (Table 1); with the most abundant subfamily being Rogadinae with 190 individual wasps, followed by Adeliinae with 87 individual wasps, Alysiinae with 34 individual wasps and Aphidiinae with 27 individual wasps (Fig. 1).

Table 1 Numbers of subfamilies, morphospecies and individuals of braconid wasps in the three fragments studied during the years 2019–2020
Fig. 1
figure 1

Variation in individual abundance of braconid wasp in the three habitat fragments studied in lowland tropical forests of Panama. Each panel shows data for a habitat fragment and is grouped by subfamily. Stacked bar plots show the dominant subfamilies within each habitat fragment

Among the three fragments, the most abundant morphospecies belonged to the subfamilies Rogadinae (M 112 y M1), Adeliinae (M 117), Microgastrinae (M 118), Aphidiinae (M 130) and Alysiinae (M 152), respectively (Table 2). The aforementioned morphospecies account for approximately 65% of the total individual wasp separated as morphospecies. In the PNM fragment, the most abundant morphospecies belonged to the subfamilies Rogadinae (M 117) and Microgastrinae (M 118) with 80 individual wasps for each, followed by Alysiinae (M 152) with 75 individual wasps, Doryctinae (M 128) with 73 individual wasps and Aphidiinae (M 130) with 73 individual wasps (Table 2). In the COR fragment, the most abundant morphospecies belonged to the subfamilies Rogadinae (M 112, M 1 and M 29) with 387 individual wasps and Aphidiinae (M 113) with 32 individual wasps (Table 2). In the ALB fragment, the most abundant morphospecies belonged to the subfamilies Rogadinae (M 1, M117, M5 and M2) with 215 individual wasps and Adelinae (M 116) with 21 individual wasps (Table 2).

Table 2 Abundance and relative abundance of morphospecies in the three fragments studied during the years 2019–2020

The Margalef index value indicated a higher species richness in the PNM, followed by ALB and COR, respectively (Table 3). The range of D´ values were close to 1 in the three fragments studied, which represents a lower degree of dominance/predominance of part of one or two species, which is also interpreted as a high diversity (Table 3). The values found with this index showed that the fragments with the highest diversity are ALB and PNM and that of COR showed a lower diversity. The H´ index indicated that the fragments with the greatest diversity were ALB and PNM and the lowest value was found in the COR fragment. According to the range of J´values all the fragments were proportional to the diversity and those that presented a more equitable distribution/pairing were ALB and PNM (Table 3). The fragment that presented less equality was COR.

Table 3 Indices of α diversity of braconid wasps in three habitat fragments in Panama

According to the Whittaker index, the highest turnover value in the composition of morphospecies recorded was between the COR and PNM fragments, followed by COR and ALB, and finally the ALB and PNM fragments (Table 4). The highest turnover value in the composition of subfamilies recorded was between the COR and PNM fragments, followed by COR and ALB, and finally the ALB and PNM fragments (Table 5).

Table 4 Whittaker’s β diversity index of morphospecies of braconid wasps in three habitat fragments in Panama
Table 5 Whittaker’s β diversity index of subfamilies of braconid wasps in three habitat fragments in Panama

The Bray–Curtis similarity analysis showed a morphospecies similarity and variation of approximately 0.5706 among the three fragments studied (Fig. 2). In the dendrogram it can be seen that there are two groupings where the morphospecies composition is more similar between the COR and ALB fragments with 0.51 (Fig. 2). The Bray–Curtis similarity analysis showed a similarity and variation of subfamilies of approximately 0.6496 among the three fragments studied (Fig. 3). In the dendrogram it can be seen that there are two groupings where the subfamily composition is more similar between the PNM and ALB fragments with 0.63 (Fig. 3).

Fig. 2
figure 2

Similarity dendrogram (Bray–Curtis) of the Braconidae morphospecies found in the three fragments studied

Fig. 3
figure 3

Similarity dendrogram (Bray–Curtis) of the Braconidae subfamilies found in the three fragments studied

The results obtained by means of the Kruskal–Wallis test determined that richness (H = 8.22, gl = 2, p = 0.015) and abundance (H = 12.95, gl = 2, p < 0.05) of braconid subfamilies observed were statistically different between the fragments. The results obtained by means of the Kruskal–Wallis test determined that richness (H = 6.053, gl = 2, p < 0.05) and abundance (H = 12.95, gl = 2, p < 0.05) of braconid morphospecies observed were statistically different between the fragments.

In the correspondence analysis (CA), 100% cumulative percentage was explained for both morphospecies and subfamilies in the first two axes. A total inertia of 0.286 was calculated for braconid subfamilies and a total inertia of 0.469 for braconid morphospecies. Wasp assemblages (both morphospecies and subfamilies) were more similar between the ALB and PNM fragment. The data obtained demonstrated distinct clustering in the three fragments, as depicted in Fig. 4 for subfamilies. The data also demonstrated that the subfamily Euphorinae was only present in the fragment of ALB and Meteorinae was present in the remnants of ALB and COR.

Fig. 4
figure 4

Correspondence analysis (CA) to compare the similarity of braconid wasp communities between habitat remnants of Corozal (COR), Albrook (ALB) and Metropolitan Natural Park (PNM). The two replicates of each habitat remnant are enveloped to make similarities among habitat remnants more apparent

Of the 16 subfamilies, 12 were shared between all fragments studied, two between the ALB and PNM fragments, and one between the ALB and COR fragments (Fig. 5). One subfamily was registered in the fragment of ALB which was not shared between the fragments studied. Of the 77 morphospecies, 29 were shared between all fragments, 11 between the ALB and PNM fragments, three between the COR and PNM fragments and three between the ALB and COR fragments (Fig. 6). Nine registered morphospecies were only found in the ALB fragment, nine different morphospecies were only found in the COR fragment and 13 different morphospecies were found in the PNM fragment, all of which were not shared between the fragments.

Fig. 5
figure 5

Venn diagram of the number of shared subfamilies of braconid wasps in the three fragments studied in lowland tropical forests of Panama: Corozal (COR), Albrook (ALB) and Metropolitan Natural Park (PNM). Of the 16 subfamilies, 12 were shared between the three habitat fragments studied

Fig. 6
figure 6

Venn diagram of the number of shared morphospecies of braconid wasps in the three fragments studied in lowland tropical forests of Panama: Corozal (COR), Albrook (ALB) and Metropolitan Natural Park (PNM). Of the 77 morphospecies, 29 were shared between the three habitat fragments studied

According to the Diserud-Odegaard Index, a similarity of 0.9418 (94%) was estimated for subfamily composition in the fragments studied. Likewise, a similarity of 0.7401 (74%) was estimated for morphospecies composition in the fragments studied. These results depicted similar subfamily and morphospecies composition within and between the fragments studied.

Discussion

Braconid wasp communities in fragmented lowland tropical forests of Panama are still very similar, both at the subfamily and morphospecies level. The fragments studied are considered to be areas large enough to maintain the biodiversity of wasp communities. The results of this study are consistent with those obtained by Valdés [31] and Ruiz-Guerra et al. [29]. Ruiz-Guerra et al. [29] demonstrated that Braconidae communities are very similar in the remnants of conserved tropical forests and secondary forests of Mexico. Valdés [31] depicted that butterfly communities are very similar in lowland tropical forests of Panama. In the study by Valdés [31], three of the four habitat fragments compared correspond to those used in the present study. The results of Valdés [31], indicated that when these three sites were compared, a similarity of 97% was obtained for the butterfly community. The proposed explanation stems from the species’ dispersal ability; since these taxa have the ability to fly, these organisms are able to move freely between fragments and as a result maintain high similarities between communities.

Richness and abundance of braconid subfamilies were statistically different between the fragments studied. The higher richness of braconid wasps found in the fragment of ALB can be explained by the heterogeneous vegetation, composed of open grassland, stubble, border vegetation and secondary forests in this fragment. The higher abundance and heterogeneity in the fragment of PNM can be explained by its proximity to a source; such as the remainder of the PNM. This is seen with the island-continent model, where local populations may be unequal in size and longevity, with one large fragment from which dispersers migrate to other fragments [32, 33]. The PNM habitat fragment is a part of the PNM; however, it has been isolated due to the development of highways. The habitat fragment however, is still considered a part of the PNM; where the remainder of the park is considered as the continent/source and the habitat fragment as the island/sink. Both fragments offer a wide range of microhabitats for organisms, allowing the survival of more individuals and increasing the availability and diversity of hosts [22]..The low richness in the fragment of COR demonstrates the introduction of biodiversity loss due to the process of fragmentation. This highlights the importance of maintaining continuous forests close to other remnants and the need to conserve fragments which provide various habitats for maintenance of species diversity.

The theory provided by the finding of braconid wasp communities is that of metapopulations [34], which assumes that species are distributed over a heterogenous space and not all territories are habitable for each species; as seen with the subfamilies of Euphorinae, Braconinae, Helconinae and Meterorinae. The fragments studied were divided into remnants/patches at a given moment primarily due to the effects of urbanization and industrialization, which separated the populations of braconid wasps forming their metapopulations. The model of the theory that supports this finding is that of patched populations [32, 34], where there are similar patches without clear distinction between sources and sinks. Two suggestions for which there is no distinction may be that the fragments have the same probability of being colonized by braconid wasps due to their dispersal ability and all fragments studied were of similar size. The fragments studied are considered to be areas large enough to keep the metapopulations of wasps in equilibrium. However, it would be necessary to prove how the dispersal dynamics of braconid wasps is carried out among fragments to prove that the fragments are indeed patched populations and determine the minimum critical size of the ecosystem to preserve the diversity and composition of species.

The results of this study also indicate that the communities of wasps (Hymenoptera: Braconidae) in the lowland forest fragments of Panama are mainly constituted by the subfamilies of Rogadinae, Alysiinae, Adeliinae, Doryctinae, Aphidiinae and Microgastrinae. The presence of these subfamilies is noteworthy since it implies interspecific relationships by parasitism between the species and other arthropods and implies that the fragments are diverse. The Microgastrinae, Adeliinae and Rogadinae subfamilies indicate interspecific relationships with Lepidoptera as they are endo parasitoids of larvae of the order Lepidoptera [35]. The Alysiinae subfamily indicates interspecific relationships with Diptera as they are endo parasitoids of the larvae or eggs of the order Diptera [35]. The Aphidiinae subfamily indicates interspecific relationships with the Stenorrhyncha order, particularly the Aphididae family [35]. The Doryctinae subfamily indicates both interspecific relationships and diversification in the fragments since they are ecto parasitoids of the larvae of the orders Coleoptera, Lepidoptera, Hymenoptera, Embioptera and Phytophagia [35]. The Doryctinae subfamily has been reported to dominate in a dry tropical forest in Mexico [36] and in a semi-deciduous forest in Venezuela [37].

There are two major limitations in this study that could be addressed in future research. First, although the appropriate trap was used, it is recommended to complement the research with another technique. An interesting technique for this community may be light traps [38, 39]. Species that are nocturnal are always attracted to light. Secondly, although the research demonstrated the composition of braconid wasp communities in three forest fragments, the research did not demonstrate the effect of habitat fragmentation on species communities. This can be demonstrated by comparing fragmented landscapes to a continuous forest or a gradient of fragment size.

Taking this into account, effective planning for the conservation and preservation of the fragments studied require learning and adaptation. A substantial set of theoretical and practical guides have been developed to maintain biodiversity and ecosystem function [40], as well as operational models that use these guides for systematic conservation ([33]; Fig. 7). With the aid of the operational model for conservation planning by Knight et al. [41], it is suggested that the fragments be considered as stepping stones. Stepping stones can improve the persistence of metapopulations allowing the flow of individuals between fragments, ensuring the exchange of stochastic local extinction and recolonization [42]. The preservation of these fairly conserved fragments, with heterogenous vegetation, can favor the presence of organisms by offering various microhabitats to ensure the viability of metapopulations. Baum et al. [43] showed that a matrix can determine if, and to what extent, corridors and stepping stones, increase the connectivity of a landscape for the survival of species using Prokelisia crocea and Spartina pectinata as indicators. Spatial configuration is particularly important for regional dynamics and must be taken into account in management plans. Small and medium fragments play a fundamental role; in such a way that it is necessary to identify and preserve these fragments.

Fig. 7
figure 7

Operational model for conservation planning, incorporating assessment, and management phases (Grantham et al. 2010; Knight et al. 2006)

Conclusions

The results of this study indicate that the richness, abundance, and composition of braconid wasps differ among habitat fragments in lowland tropical forests of Panama. Forest fragments have relatively high similarity of braconid wasp subfamilies and morphospecies. This finding was interpreted as an indication that a species’ dispersal ability plays a major role in its survival.

Annual inventories to assess the real change in braconid wasp communities or other organisms could provide critical monitoring of the effects of habitat fragmentation. These results can help to increase the understanding of the influence of habitat fragmentation on braconid wasp communities in Panama to develop successful conservation strategies. This would make it possible to address the question of the effect of fragmentation on the braconid wasp species.

Methods

Study sites

The three selected study sites in Panama City were Metropolitan Natural Park (PNM), public land in the town of Albrook (ALB) and public land in the town of Corozal (COR) (Fig. 8). Each fragment was bordered and labeled; as depicted in the map constructed with Google Earth. The PNM fragment is situated 8°59′41.55′′N and 79°32′35.22′′ W with an area of approximately 18.12 ha and a perimeter of 1.756 km. Vegetation is characterized as a mixture of tropical humid forest and lowland tropical dry forest, with few areas of stubble and grasslands, and a well-defined stratum. The ALB fragment is situated 8°58′37.49′′N and 79°33′43.82′′W with an area of approximately 34.79 ha and a perimeter of 5.003 km. Vegetation is characterized as heterogenous, composed of open grasslands, stubble and secondary forests. The COR fragment is situated 8°59′19.34′′N and 79°34′11.83′′W with an area of approximately 56.31 ha and a perimeter of 3.028 km. Vegetation is characterized as herbaceous with late secondary forest and some open areas. In this study, the habitat fragmentation definition defined by Saunders et al. [44] was used, where they described a fragment as any patch of native vegetation around which most of the original vegetation has been removed. For this reason, well-preserved fragments that are separated by a matrix of urban areas were considered as fragments in Panama City. The selected fragments are close in proximity to each other, which guarantees that the results obtained are effects of the fragmentation process of an original matrix and not the natural result caused by the distances between them. The distance between the fragments of ALB and COR is 1.129 km. The distance between the fragments of COR and PNM is 2.565 km. The distance between the fragments of ALB and PNM is 2.329 km.

Fig. 8
figure 8

Map of the Pacific Basin of the Panama Canal where the habitat fragments studied in lowland tropical forests of Panama are located: Corozal (COR), Albrook (ALB) and Metropolitan Natural Park (PNM). Google Earth was used in order to construct the map

The study sites form part of a biological corridor that runs along the eastern bank of the Panama Canal [45]. In accordance with the UNESCO classification, the three study sites are characterized as lowland tropical semi-deciduous forests. The three sites present annual average temperatures of 26.4 °C, with an annual average precipitation between 1501 to 1800 mm and altitudes between 20 to 150 masl. In the last 171 years, the forests along the eastern bank of the Panama Canal have been pressured by anthropic activities [46]. One of the first anthropic impact along the banks of the Panama Canal was the construction of Panama’s railway in 1850. The second anthropic impact along the banks of the Panama Canal was the construction of the French Canal in 1881. Nevertheless, the greatest anthropic impact started in the 1900’s when the forests along the eastern bank of the Panama Canal were intervened by military bases [46]. During this time, the remaining fragments maintained advanced secondary forest vegetation. In the last 50 years, these sites have been anthropically pressured due to the construction of neighborhoods along the banks of the Panama Canal [46]. The COR fragment has been used for Panama government security training activities. The ALB fragment has minimal human traffic; however, it contains a water reservoir at the peak of the fragment. The PNM fragment is along a hiking trail in a protected area but with visitor traffic [47].

Sampling protocol

Two permanent lightweight Malaise Trap, Townes Style separated by approximately 0.2 km were installed in the center of each fragment and wasp samples were collected weekly from May 2019 to March 2020. The traps were made with organza fabric, with dimensions of 5.8 ft tall by 5.4 ft long, and contained a polyethylene collector bottle with 95% alcohol at its highest end [48]. This is a trap designed to collect fast-flying insects whose behavior is to fly upwards when it touches a surface. The Malaise Trap, Townes Style has been widely used in braconid wasps diversity studies in Central America and the world [49,50,51,52,53]. A reasonable flat, log-less area of approximately 2 m2 was chosen for the placement location of each trap in the forest. The collection bottle was placed facing the magnetic north and half-filled with 95% ethanol. The traps were placed in open areas which served as insect corridors. To exclude the edge effect, a safety distance of approximately 0.2 km was defined from the border of each fragment towards the center. Samples were collected every 7 days; resulting in a sampling effort of 44 weeks per trap for the period between 2019 and 2020.

Samples were taken to the Master’s degree laboratory at the University of Panama, where braconid specimens were separated from the rest of the sample and placed in vials of one dran with 95% ethanol. All braconid individuals were mounted on entomological pins, number 2 or 3. Braconid individuals were first sorted as morphospecies (M) using the method of Oliver and Beattie [54]. Oliver and Beattie [54] demonstrated that comparisons can be made using morphospecies assemblages as long as each contain a unique identification. Samples were then identified to subfamily level using the Sharkey and Campos taxonomic key from the book Aguilar et al. [35] and the manual by Sharkey et al. [26].

Statistical analyses

In order to measure the diversity of subfamilies and morphospecies over a special scale, alpha (α) and beta (β) diversity were calculated. α diversity is the species richness of a particular community that is consider to be homogeneous. Currently there are many indices to measure α diversity. In this study, to determine the α of each fragment, the Simpson dominance index (D′), Shannon-Weiner Index (H′), Equity index or Pielou equity (J′) and Margalef diversity index were used. D′ takes into account the most important species without considering the rest of the species.

$${D}^{^{\prime}}=1-\sum ({{p}_{i})}^{2},$$

where pi is the proportional abundance of species i (the number of individuals of species i divided by the total number of individuals in the sample).

H' combines information on species richness and equity in what is called diversity or heterogeneity [55, 56]. The average degree of uncertainty is measured in predicting the species to which a given randomly chosen individual within a biotic community belongs.

$$H^{\prime}= -\sum_{i=1}{p}_{i} ln {p}_{i},$$

where ln is the natural logarithm and pi is the proportional abundance of species i (the number of individuals of species i divided by the total number of individuals in the sample).

J′ is the maximum possible diversity for a given number of species that occurs if all species are present in equal numbers. Its value ranges from 0 to 1, so that 1 corresponds to situations where all species are equally abundant.

$${J}^{^{\prime}}=\frac{H{^{\prime}}}{{H}^{^{\prime}}max},$$

where H' max is the lnS, H’ is the Shannon-Weiner index.

Margalef diversity index estimates the biodiversity of a community based on the numerical distribution of the individuals of the different species based on the number of individuals in the analyzed sample.

$${D}_{Mg}=\frac{S-1}{\mathrm{ln}N},$$

where S is the number of species and N is the total number of individuals.

The β diversity is the diversity of species between communities. It is the degree of species replacement or biotic change through environmental gradients [57, 58]. The indices used to determine β diversity are approaches based on pairwise similarities. To determine the β of each fragment, the Whittaker and Bray–Curtis Indices were used. The Bray Curtis Index is a modified version of the Sorensen Index [59]

$${S}_{S}=\frac{2a}{2a+b+c},$$

where a is the number of species present in both samples, b is the number of species found in community A and c is the number of species found in community B. Whittaker index estimates the degree of species replacement or biotic change through environmental gradients [60]. The α and β diversity indices were estimated using the EstimateS and the “vegan” package from RStudio.

For analysis of the data, all samples were pooled monthly, resulting in one sample per site. A Kruskal–Wallis test was performed using the program STATISTICA [61], to determine differences between the number of individuals registered in each study site. Correspondence analysis (CA) calculated by the program XL-Stat [62] was used to characterize the braconid community based on the number of individuals and site. A comparison of the composition of the braconid community was calculated using formula for the Multiple Similarity Index by Diserud and Odegaard [63], as given by the Eq. (2):

$${\text{C}}_{\text{S}}^{\text{T}}\text{=}\frac{\text{T}}{\text{T-1}}\left(\text{1} - \frac{{\text{S}}_{\text{T}}}{\sum_{\text{i}}{{\text{a}}}_{\text{i}}}\right),$$

where ai is the number of species in site Ai, T is the number of sites and ST = \(\sum_{i}{a}_{i}-\sum_{i<j}{a}_{ij}+ \sum_{i<j<k}{a}_{ijk}\dots ,\) aij is the number of species shared by sites Ai and Aj; and aijk is the number of species shared by sites Ai, Aj and Ak, etc. The multiple similarity index takes into consideration the information of species shared by two or more sites and avoids the problem of covariance between pairwise similarities [63], reducing the probability of Type II errors. The multiple similarity index avoids the loss of information concerning the number of species shared among three or more sites and the lack of independence between pairwise similarities due to the repetition of each site in several pairs as seen with β diversity indices [63]. The calculation of this index was done manually.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ALB:

Albrook

CA:

Corresponding analysis

COR:

Corozal

FRAG:

Fragmentation Index

M:

Morphospecies

M_#:

Morphospecies number

PNM:

Metropolitan Natural Park

References

  1. Fahrig L. Effects of habitat fragmentation on biodiversity. Annu Rev Ecol Evol Syst. 2003;34:487–515.

    Article  Google Scholar 

  2. Fontúrbel FE. Evaluation of the loss of the Chaqueño Dry Forest coverage in Torotoro’s Municipality and in Torotoro National Park (Potosí, Bolivia), through Remote Sensing. Ecol Apl. 2007;6(12):59–66.

    Article  Google Scholar 

  3. Rands MRW, Adams WM, Bennun L, Butchart SHM, Clements A, Coomes D, et al. Biodiversity conservation: challenges beyond 2010. Science (80). 2010;329:1298–303.

    Article  CAS  Google Scholar 

  4. Zurita MEC. Análisis de la fragmentación en el Parque Nacional Cerro Azul Meámbar (PANACAM). 2014.

  5. Macedo-Reis LE, Quesada M, de Siqueira NF. Forest cover drives insect guild diversity at different landscape scales in tropical dry forests. For Ecol Manage. 2019;443:36–42.

    Article  Google Scholar 

  6. Rogan JE, Lacher TE. Impacts of habitat loss and fragmentation on terrestrial biodiversity. In: Reference module in earth systems and environmental sciences. Amsterdam: Elsevier Inc; 2018. p. 1–18.

    Google Scholar 

  7. Mas J-F, Sandoval JC. Analysis of landscape fragmentation in the Petenes protected area, Campeche, Mexico. Investig Geográficas. 2000;6:42–59.

    Google Scholar 

  8. Tittensor DP, Walpole M, Hill SLL, Boyce DG, Britten GL, Burgess ND, et al. A mid-term analysis of progress toward international biodiversity targets. Science (80-). 2014;346(6206):241–4.

    Article  CAS  Google Scholar 

  9. Chazdon RL, Harvey CA, Komar O, Griffith DM, Ferguson BG, Martínez-Ramos M, et al. Beyond reserves: a research agenda for conserving biodiversity in human-modified tropical landscapes. Biotropica. 2009;41(2):142–53.

    Article  Google Scholar 

  10. Barnosky AD, Hadly EA, Bascompte J, Berlow EL, Brown JH, Fortelius M, et al. Approaching a state shift in Earth’s biosphere. Nature. 2012;486(7401):52–8.

    Article  CAS  PubMed  Google Scholar 

  11. FAO. El estado de los bosques del mundo. 2016.

  12. SINIA. Diagnóstico de cobertura de Bosques y Otras Tierras Boscosas en la República de Panamá: año 2019 [Internet]. 2020 [cited 2020 Jul 30]. Available from: https://www.sinia.gob.pa/index.php/cobertura-boscosa/ano-2019/a-nivel-nacional.

  13. Haddad NM, Brudvig LA, Clobert J, Davies KF, Gonzalez A, Holt RD, et al. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Appl Ecol. 2015.

  14. Franklin AB, Noon BR, George TL. What is habitat fragmentation? Stud Avian Biol. 2002;25:20–9.

    Google Scholar 

  15. Laurance WF. Comparative responses of five arboreal marsupials to tropical forest fragmentation. J Mammal. 1990;71(4):641–53.

    Article  Google Scholar 

  16. Conner RN, Rudolph DC. Forest habitat loss, fragmentation, and red-cockaded woodpecker populations. NCASI Tech Bull. 1991;103(3):446–57.

    Google Scholar 

  17. Newmark WD. Tropical forest fragmentation and the local extinction of understory birds in the Eastern Usambara Mountains, Tanzania. Conserv Biol. 1991;5:67–78.

    Article  Google Scholar 

  18. Didham RK, Hammond PM, Lawton JH, Eggleton P, Stork NE. Beetle species responses to tropical forest fragmentation. Ecol Monogr. 1998;68(3):295–323.

    Article  Google Scholar 

  19. Bergman KO, Dániel-Ferreira J, Milberg P, Öckinger E, Westerberg L. Butterflies in Swedish grasslands benefit from forest and respond to landscape composition at different spatial scales. Landsc Ecol. 2018;33:2189–204.

    Article  Google Scholar 

  20. Andrade AC, Franzini LD, Mesquita DO. Assessing the effect of urbanization on tropical forest dwelling teiid lizards. Ecol Indic. 2019;99:225–9.

    Article  Google Scholar 

  21. Salomão RP, Alvarado F, Baena-Díaz F, Favila ME, Iannuzzi L, Liberal CN, et al. Urbanization effects on dung beetle assemblages in a tropical city. Ecol Indic. 2019;103:665–75.

    Article  Google Scholar 

  22. Graça MB, Somavilla A. Effects of forest fragmentation on community patterns of social wasps (Hymenoptera: Vespidae) in Central Amazon. Austral Entomol. 2018;58(3):657–65.

    Article  Google Scholar 

  23. Basset Y, Novotny V, Miller SE, Kitching RL. Arthropods of tropical forests: spatio-temporal dynamics and resource use in the canopy. In: Basset Y, Novotny V, Miller SE, Kitching RL, editors. 1st edn. Cambridge University Press, Cambridge, UK; 2003.

  24. Basset Y, Lamarre GPA. Toward a world that values insects. Science (80-). 2019;364(6447):1230–1.

    Article  CAS  Google Scholar 

  25. González HD, Ruíz DB. Los braconidos (Hymenoptera: Braconidae) como grupo parametro de biodiversidad en las selvas deciduas del tropico: una discusion acerca de su posible uso. Acta Zool Mex. 2000;79:43–56.

    Article  Google Scholar 

  26. Sharkey M, Dangerfield P, Fernández-Triana J, Kittel R, Kula R, Marsh P, et al. Manual of the new world genera of the family braconidae (Hymenoptera). In: Wharton R, Marsh P, Sharkey M, editors. Annals of the Entomological Society of America. 2nd edn. Washington, DC: The International Society of Hymenopterists; 1997.

  27. Whitfield JB. Estimating the age of the polydnavirus/braconid wasp symbiosis. Proc Natl Acad Sci. 2002;99(11):7508–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Smith TJ, Mayfield MM. The effect of habitat fragmentation on the bee visitor assemblages of three Australian tropical rainforest tree species. Ecol Evol. 2018;8:8204–16.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Ruiz-Guerra B, López-Acosta JC, Zalsivar-Riverón A, Velázquez-Rosas N. Braconidae (Hymenoptera: Ichneumonoidea) abundance and richness in four types of land use and preserved rain forest in southern Mexico. Rev Mex Biodivers. 2015;86:164–71. https://doi.org/10.7550/rmb.43865.

    Article  Google Scholar 

  30. Wright SJ, Muller-Landau HC. The future of tropical forest species. Biotropica. 2006;38(3):287–301.

    Article  Google Scholar 

  31. Valdés L. Estructura de las comunidades de mariposas diurnas (lepidoptera: rhopalocera) en fragmentos de vegetación en la ciudad de panamá. Universidad de Panamá; 2018.

  32. Leibold MA, Holyoak M, Mouquet N, Amarasekare P, Chase JM, Hoopes MF, et al. The metacommunity concept: a framework for multi-scale community ecology. Ecol Lett. 2004;7(7):601–13.

    Article  Google Scholar 

  33. Chase J, Amarasekare P, Cottenic K, Gonzalez A, Holt R, Holyoal M, et al. Competing theories for competitive metacommunities. In: Holyoak M, Leibold MA, Holt RD, editors. Metacommunity ecology: emerging views of community structure and dynamics above the local scale. University of Chicago Press; 2005. p. 335–54.

  34. Hanski I, Gilpin ME. Metapopulation dynamics: from concepts and observations to predictive models. In: Metapopulation biology. Academic Press; 1997. p. xv–xvi.

  35. Aguilar A, Brothers DJ, Amarante S, Buffington M, Arias-Penna D, Campos D, et al. Introducción a los Hymenoptera de la Región Neotropical. In: Fernández F, Sharkey MJ, editors. vol. 53. 2006; p. 93–114.

  36. Coronado-Blanco JM, Zaldívar-Riverón A. Biodiversidad de Braconidae (Hymenoptera: Ichneumonoidea) en México. Rev Mex Biodivers [Internet]. 2014;85:372–8.

    Article  Google Scholar 

  37. Briceño R, Torres D, Zaldívar-Riveron A. Primer reporte de la familia braconidae (hymenoptera: ichneumonoidea) en los parques cerro saroche y terepaima, Estado Lara. Venezuela Bioagro. 2009;21(3):223–6.

    Google Scholar 

  38. Basset Y et al. IBISCA-Panama, a large-scale study of arthropod beta-diversity and vertical stratification in a lowland rainforest: rationale, study sites and field protocols. Bull Inst R Des Sci Nat Belgique. 2007;77.

  39. Basset Y, Cizek L, Cuénoud P, Didham RK, Guilhaumon F, Missa O, et al. Arthropod diversity in a tropical forest. Science (80-). 2012;338:1481–4.

    Article  CAS  Google Scholar 

  40. Margules CR, Sarkar S. Systematic conservation planning. Cambridge, UK: Cambridge University Press; 2007.

    Google Scholar 

  41. Knight AT, Cowling RM, Campbell BM. An operational model for implementing conservation action. Conserv Biol. 2006;20(2):408–19.

    Article  PubMed  Google Scholar 

  42. Sawchik J, Dufrêne M, Lebrun P, Schtickzelle N, Baguette M. Metapopulation dynamics of the bog fritillary butterfly: modelling the effect of habitat fragmentation. Acta Oecolog. 2002;23(5):287–96.

    Article  Google Scholar 

  43. Baum KA, Haynes KJ, Dillemuth FP, Cronin JT. The matrix enhances the effectiveness of corridors. Ecology. 2004;85(10):2671–6.

    Article  Google Scholar 

  44. Saunders DA, Hobbs RJ, Margules CR. Biological consequences of ecosystem fragmentation: a review. Conserv Biol. 1991;5(1):18–32.

    Article  Google Scholar 

  45. Miambiente. Áreas protegidas y Vida Silvestre. Parque Natural Metropolitano [Internet]. 2017. Available from: http://www.miambiente.gob.pa/index.php/sistema-nacional-de-areas-protegidas-sinap.

  46. Pinilla ELG. Old North American military bunkers located in the corners of the Panama Canal for tourist purposes. Contact. 2021;1:60–77.

    Google Scholar 

  47. Correa M, Stapf M, DeSedas A, Hernández F, Carranza R. Árboles y Arbustos del Parque Natural Metropolitano, Panamá. 1a ed. Morales FJ, editor. Santo Domingo de Heredia: Instituto Nacionalde Biodiversidad, INBio; 2010.

  48. Townes H. A light-weight malaise trap. Ent News. 1972;83:239–47.

    Google Scholar 

  49. Marchiori CH, Penteado-Dias AM. Braconidae (Hymenoptera: Ichneumonoidea) collected in a native forest area in Itumbiara, Goiás, Brazil. Neotrop Entomol. 2002;31:647–9.

    Article  Google Scholar 

  50. Querino R, Couceiro S, Queiroz L, Penteado-Dias A. The spatial distribution of Hymenoptera parasitoids in a forest reserve in Central Amazonia, Manaus, AM, Brazil. Braz J Biol. 2011;71(4):865–71.

    Article  Google Scholar 

  51. An T, Belokobylskij SA, Kim H, Kim K, Seo H, Ku D. New faunistic data on the family Braconidae (Hymenoptera) from Korea. J Asia-Pacific Biodivers. 2014;7:489–92.

    Article  Google Scholar 

  52. Reséndiz-Flores A, Nunes JF, García-Paris M, Zaldivar-Riverón A. Six new species of the parasitoid wasp genus Notiospathius (Hymenoptera: Braconidae: Doryctinae) from Mexico. Rev Mex Biodivers. 2014;85:391–401.

    Article  Google Scholar 

  53. Rodriguez-Mota AJ, Ruiíz-Cancino E, Ivanovich-Khalaim A, Coronado-Blanco JM, Treviño-Carreón J. Diversidad de Ichneumonidae (Hymenoptera) en un bosque de Pinus spp y. Juniperus flaccida en Jaumave, Tamaulipas, México. Rev Mex Biodivers. 2015;86:972–80.

    Article  Google Scholar 

  54. Oliver I, Beattie AJ. Invertebrate morphospecies as surrogates for species: a case study. Conserv Biol. 1996;10(1):99–109.

    Article  Google Scholar 

  55. Magurran AE. Measuring biological diversity. Oxford: Blackwell Publishing; 2004.

    Google Scholar 

  56. Moreno C, Zuria I, García-zenteno M, Sánchez-rojas G, Castellanos I, Martínez-morales M, et al. Trends in the measurement of alppha diversity in the last two decades. Interciencia. 2006;31(1).

  57. Peet R. The measurement of species diversity. Annu Rev Ecol Syst. 1974;5:285–307.

    Article  Google Scholar 

  58. Moreno CE, Patrón JC, Regalado NM, Falcón APM, Martinez IJO, Díaz CLR, et al. Measuring species diversity in the tropics: a review of methodological approaches and framework for future studies. Biotropica. 2018;50(6):929–41.

    Article  Google Scholar 

  59. Chao A, Chazdon RL, Shen TJ. A new statistical approach for assessing similarity of species composition with incidence and abundance data. Ecol Lett. 2005;8:148–59.

    Article  Google Scholar 

  60. Koleff P. Conceptos y medidas de la diversidad. In: Sobre diversidad biológica: el significado de las diversidades alfa, beta y gamma. 2005.

  61. Statsoft. STATISTICA [Internet]. 2016. Available from: http://www.statsoft.com/Products/STATISTICA.

  62. Addinsoft. XLStat statistical and data analysis solution [Internet]. New York, USA; 2021. Available from: https://www.xlstat.com.

  63. Diserud OH, Odegaard F. A multiple-site similarity measure. Biol Lett. 2007;3:20–2.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

To the University of Panama, Secretaría Nacional de Ciencia Tecnología e Innovación (SENACYT) and Sistema Nacional de Investigación (SNI) for providing funds and resources to carry out research activities. To MIAmbiente (Ministry of Environment) and Metropolitan Natural Park for issuing collection permits. To Marcelo Mack and JeanCarlos Abrego for their help with sample collection.

Funding

This research was funded by Secretaría Nacional de Ciencia, Tecnología e Innovación (SENACYT), Sistema Nacional de Investigación (SNI) and the University of Panama (UP). The funding bodies played no role in the design of the study collection, analysis and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

LAR and EM contributed to the study conception and design. Material preparation, data collection and analysis were performed by LAR. The first draft of the manuscript was written by LAR. EM and LAR commented on previous versions of the manuscript and have substantively revised it. EM and LAR read and approved the final manuscript.

Corresponding author

Correspondence to Enrique Medianero.

Ethics declarations

Ethics approval and consent to participate

A scientific collecting permit was obtained from the Ministry of Environment to collect specimens. A letter of authorization was also obtained from the manager of Metropolitan Natural Park to collect specimens inside the park.

Consent to publication

Not applicable.

Competing interests

Not applicable.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rodríguez, L.A., Medianero, E. The composition of braconid wasp communities in three forest fragments in a tropical lowland forest of Panama. BMC Ecol Evo 22, 98 (2022). https://doi.org/10.1186/s12862-022-02051-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-022-02051-4

Keywords