<figure><img src="image_1.png" /></figure>

HAL
open science

# Large-scale sampling of beetle communities in Laos shows that conversion of natural forests into plantations leads to a decline in family richness and abundance

Bounsanong Chouangthavy, Yoan Fourcade

► To cite this version:

Bounsanong Chouangthavy, Yoan Fourcade. Large-scale sampling of beetle communities in Laos shows that conversion of natural forests into plantations leads to a decline in family richness and abundance. Ecology and Evolution, 2023, 13 (7), pp.e10258. 10.1002/ece3.10258. hal-04199799

HAL Id: hal-04199799

https://hal.u-pec.fr/hal-04199799

Submitted on 8 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Large-scale sampling of beetle communities in Laos shows that conversion of natural forests into plantations leads to a decline in family richness and abundance

Bounsanong Chouangthavy<fnref n="1" /> ID |

Yan Fourcade<fnref n="2" /> ID

### Correspondence

Bounsanong Chouangthavy, Entomology Laboratory, Faculty of Agriculture, National University of Laos, P.O. Box 7322, Vientiane, Laos.

Email: bsnchouangthavy@gmail.com

### Abstract

Rapid economic development can pose a threat to the biodiversity of tropical countries. In Laos, this is manifested by the conversion of natural forests into plantations, even though this area is one of the biodiversity hotspots of Southeast Asia. Beetle communities can be good indicators of the impact of anthropogenic pressure on natural ecosystems. In this study, we analyzed for the first time a large-scale inventory of Coleoptera to assess the ecological and anthropogenic drivers of beetle communities in Laos. We examined beetle communities (described at the family level) across the country, located in distinct habitat types, in order to understand the impact of the conversion of natural forest into plantations. We found that beetle abundance had declined in plantations compared to natural forests. At the same time, we observed fewer beetle families in plantations overall, but at the scale of sampling sites there was no difference in local richness compared to natural forests, suggesting a homogenization of beetle communities in anthropogenic habitats. Although results are certainly sensitive to our coarse classification of beetle specimens into families, the negative impact of the conversion of natural tropical forests into agriculture area can still be clearly demonstrated. Our findings highlight that it is possible to make use of unstructured large-scale inventories to explore how beetle communities responds to landscape changes induced by human activities. We suggest that sampling beetle communities can be used as an ecological indicator to monitor anthropogenic impacts on tropical ecosystems.

### KEYWORDS

agriculture, beetles, biodiversity, deforestation, landscape contexts, rapid change

### TAXONOMY CLASSIFICATION

Biodiversity ecology

This is an open access article under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2023 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.

<footnote n="1">¹Entomology Laboratory, Faculty of Agriculture, National University of Laos, Vientiane, Laos</footnote>

<footnote n="2">²Univ. Paris Est Creteil, Sorbonne Université, Univ Paris Cité, CNRS, IRD, INRAE, Institut d'écologie et des sciences de l'environnement, IEES, Créteil, France</footnote>
# 1 | INTRODUCTION

Southeast Asian tropical forests are home to a rich biodiversity, and are essential habitats and food supply for humans and animals (Estoque et al., 2019). Southeast Asia's tropical lowland rainforests are among the most diverse ecoregions on earth, where a high proportion of endemic species and a high rate of habitat degradation coexist (Myers et al., 2000). However, huge amounts of forests are being destroyed for agricultural cash crops and industrial tree plantations, leading to land use change and intensification that are primarily driven by an increase in the worldwide demand for agricultural products (Kusuma et al., 2018). Insects, for which a global decline has been recently documented (Sánchez-Bayo & Wyckhuys, 2019), are among the taxa that are known to be affected by deforestation in a tropical context (Correa-Carmona et al., 2022). Generally, deforestation, agricultural intensification, and climate change, including more frequent extreme weather events, have been suggested as being the major drivers of the global insect decline (Eggleton, 2020; Wagner, 2020). However, this assessment mostly comes from population trends estimated in the Global North (Sánchez-Bayo & Wyckhuys, 2019), while the current state of insect diversity in tropical contexts remains poorly known.

Laos covers a large part of the Indochinese limestone belt and is one of the most biodiversity-rich countries of Southeast Asia (Kumar et al., 2016). Its dominant land cover is tropical dense forests, of which approximately 80% are located in mountainous areas with steep to moderate slopes. It includes several National Protect Areas (NPAs) that are considered as biodiversity hotspots, such as the Hin Nam No NPA, which has officially been submitted to become the first natural World Heritage Site of Laos. Many studies report a highly diverse fauna in the region including amphibians, reptiles, birds, bats, and over 100 species of large mammals, new species being frequently discovered (Ceballos & Ehrlich, 2006; MAF & STEA, 2003; MoNRE, 2016; Myers et al., 2000). However, as of yet, no single survey has attempted to describe the whole diversity of Coleoptera, or even insects in general, in Laos. Until fairly recently, the insect fauna of Laos remained one of the most poorly known in Southeast Asia (Sekerka & Geiser, 2016) and existing knowledge mostly comes from specimens collected by foreign visitors before the 1920s. Recently, though, we observed an increase in the number of entomological expeditions, permitted by the country becoming more accessible to foreigners (Chouangthavy et al., 2020).

The extremely rapid economic growth that Laos is experiencing comes at the expense of biodiversity, which is facing a growing number of significant challenges associated with land use changes (World Bank national accounts data, 2017). For example, during the 1990s and 2000s, the land area dedicated to rubber plantations has increased exponentially from 115,732 ha to reach an evaluated surface of 450,000 ha in 2015 (Smith et al., 2020). Such conversion of natural tropical forests into rubber plantations occurs in several tropical regions of the world where it is recognized to negatively impact biodiversity and ecosystems (Warren-Thomas et al., 2015). Moreover, the economic growth of the region is likely to continue

or even accelerate in the near future, as the railway that connects Kunming, China to Bangkok, Thailand, passing through much of Laos, is completed (Chen & Haynes, 2017; Ng et al., 2020). Indeed, infrastructure development will increase the general appeal of the region and encourage foreign investment, contributing to direct and indirect threats to local ecosystems (Borda-de Água et al., 2017; Torres et al., 2016). Knowledge of the influence of human impact, through an effect on landscape structure, on insect diversity in southern Asia, and in Laos in particular, remains insufficient (but see e.g., Chouangthavy et al., 2020).

Beetles (Coleoptera) are the most diverse taxonomic order on Earth. Because they exhibit rich abundance, biomass, and diversity, beetles are often used as indicator species to estimate anthropogenic impact on ecosystems, including in tropical forests (Ghannem et al., 2018; Parikh et al., 2021; Zödl & Wittmann, 2003). For example, dung beetles (Scarabaeidae: Sacarabaeinae) play an important role in the functioning of tropical forest ecosystems while being also sensitive to human disturbance and environmental changes, making them ideal focal species for investigating conservation issues (Slade et al., 2011; Spector, 2006). In southern Asia and in Laos in particular, despite the fact that beetle diversity is high (Moodley et al., 2022) and human activity is growing, the approach of employing beetle community composition, richness, and abundance as surrogates for estimating the impact of agricultural intensification and anthropogenic disturbance has been rarely carried out, especially in Laos. The few studies addressing the question of beetle community richness in relation to human factors were restricted to specific local contexts (Chouangthavy et al., 2020), and were insufficient to estimate more broadly the actual impact of human pressures on beetle biodiversity in Laos. There is thus a need for large-scale assessments of beetle diversity conducted in natural versus anthropogenic landscapes, in order to estimate the effect of Laos' land use change on its rich biodiversity.

Besides human activity, there is evidence that beetle diversity is also partly structured by climate at large spatial scales (Andrew & Hughes, 2004; Hortal et al., 2011). Even at a more regional scale, beetle assemblages appear to be structured across elevation gradients following the corresponding climatic variation (Dolson et al., 2021; Gebert et al., 2020). This implies that (i) beetle species may be affected by climate change in the recent past and in the future (Harris et al., 2019), and (ii) any attempt to characterize the impact of anthropogenic factors on beetle diversity at a large geographical scale must also account for climatic gradients that may influence the richness, abundance, and composition of beetle communities. Therefore, it is likely that the beetle fauna of Laos, a country that covers a latitudinal gradient of ca. 900km and hosts four different climate zones (Am, Aw, Cwa, Cwd; Essenwanger, 2001), is somewhat influenced by variation in temperature and precipitation across the country.

In order to understand how the rapid land use change that follows the economic development of the region affects its biodiversity, we investigated beetle community composition, diversity, and abundance across a large spatial scale in Laos (focusing on northern
and southern areas, and excluding three provinces in the center), focusing on two contrasting landscape contexts. Specifically, we made use of an unprecedentedly large inventory of beetles carried out in various locations in the country to compare beetle assemblages (characterized at the family level) in natural forests and in plantations, accounting also for climatic gradients that may be an additional driving force of the composition and diversity of beetle communities at macrogeographical scales. We hypothesize that the disturbance caused by anthropogenization results in differences in beetle community composition between natural forests compared to plantations, and that family richness and abundance are higher in natural forests compared to plantations. This work provides an assessment of the effect of landscape context and anthropization on beetle diversity that is still rarely carried out at such a large scale in this region of the world.

## 2 | MATERIALS AND METHODS

### 2.1 | Landscape context

The study was conducted in 12 locations that cover the northern and southern parts of Laos, leaving only three provinces in central Laos unsampled (Figure 1), in two contrasting landscape contexts (Table 1). Five sampling locations correspond to relatively undis turbed natural forests, which have long been recognized for their outstanding biodiversity. There, the traditional human activity consists mainly of logging, food searching, and hunting. However, natural forests have recently been facing land use intensification pressures, due to extremely rapid economic growth leading to the

development of multiple aspects of human activities, including the expansion of agricultural lands. Part of the natural forest in the study areas has been influenced by the construction of a railroad going from the north to the center of Laos, which will be part of a larger railway linking China to Thailand through Laos.

In contrast, seven sampling locations were located in different types of agricultural landscapes, hereafter referred to as plantations. This type of landscape structure had different agricultural and deforestation histories, and is mostly dedicated to rubber and eucalyptus plantations that have increased exponentially in the country. Part of the present study was undertaken on a large rub ber plantation in northern Laos (site 4), which covered 33,000ha in 2016 and was occupied by various ethnic groups (Kusakabe & Chanthoumphone, 2021); more rubber plantation areas now extend to the center and to the south. Most of the deforestation happened several years before sampling, resulting in a single unfragmented rubber plantation. In addition, sampling also occurred in a eucalyptus plantation planted with a monocrop of eucalyptus trees in rows, the natural forest around plantations being used for extensive livestock and rice paddy field (site 1, 2, and 3). Furthermore, different other types of agricultural activities such as coffee, cabbage, strawberry farms, and grassland mixed with rice paddy fields were covered in the present study across the country (Figure 1).

### 2.2 | Beetle sampling and taxonomic assignment

We explored beetle communities across the country in the two distinct habitat types, applying different trapping methods (pitfall traps, window traps, light trapping, and hand collection, in order to

<figure><img src="image_3.png" /><figcaption>FIGURE 1 Map showing the distribution of sampling areas in 2018-19-20-21 across Laos, with land cover (from the Copernicus Global Land Service project; data from 2019) shown as background. Note that some sites were sampled using different methods and are thus figured with different shapes.</figcaption></figure>
capture different aspects of community composition) throughout four consecutive years (2018-19-20-21, depending on sampling site). Pitfall traps consisted of plastic cups (300cm³) with a diameter of 8cm at the top, buried in the ground so that the top rim was flush with the soil surface. The window traps were made of 1.5-L transparent plastic drinking water bottles with one window and were suspended upside-down. The traps were hung on wooden poles at approximately 150cm above ground level. Both trap types were filled with about 50mL of 70% alcohol and had a cover for rain protection during the collections carried out in the rainy season. Traps were left 5days in the field, then all captured insects were brought back to the laboratory, where beetles were separated from other species and debris and sorted into families under a microscope. Light trapping consisted of a 125v bulb and white clothing of 4.5m² hung up between trees. All beetle specimens were directly stored in 95% alcohol after sorting. Furthermore, hand collection was done by the first author during visits to the field. Because of taxonomic uncertainties and difficulties to identify species precisely, all beetle specimens were assigned to the family level. Note that like several other studies have reported the use of higher insect taxa (i.e., family level) to investigate the impact of habitats and environmental changes (Báldi, 2003; González et al., 2015; Parikh et al., 2021). We then described diversity and community composition across landscape contexts based on the identity and abundance of beetle families in each combination of sampling site, trapping method, and sampling date.

Some sites were sampled each month for several years using different trapping methods and many traps, whereas in some sampling sites collection occurred only occasionally during a few days in a single year, using one or a couple of complementary trapping methods, resulting in a small number of traps. For example, four sampling sites were equipped with pitfall (2 × 30) and window (2 × 30) traps in 2020 (see Table 1). Then, in these sites, sampling occurred each month of the year for 5 days, resulting in a particularly large sampling effort (60 traps × 5 days × 12 months). In contrast, a single sampling site was equipped with 15 window traps for 9 days in 2018. Six sites were sampled by hand collection during a few days in 2019 and 2021. Finally, light trapping was used for 3–5 days in four sampling sites in 2019 and 2021 (Table 1). All beetle specimens were housed at the Entomology Laboratory, Faculty of Agriculture, National University of Laos for future identification.

## 2.3 | Climate data

Climate variables used in this study were based on mean annual temperature and precipitation obtained from the Worldclim web portal version 2.1 (Fick & Hijmans, 2017). The variables of Worldclim are raster surfaces derived from the interpolation of weather stations' data collected across the period 1970–2000. We imported a shapefile of Laos boundaries from the GADM database (version 3.6, https://gadm.org/), the coordinates of sampling sites and the rasters of annual precipitation and mean annual temperature from Worldclim into QGIS v.3.22.7. Then, we extracted temperature and

precipitation at each sampling locations as a proxy for the climatic variation that exists across sampling sites distributed all over Laos.

## 2.4 | Statistical analyses

We first tested the effect of four independent response variables (landscape context [i.e., natural forest vs. plantations], trapping method, temperature, and precipitation) on beetle community composition using a permutational multivariate analysis of variance (PERMANOVA) as implemented in the "vegan" R package (Oksanen et al., 2022), using Jaccard distance as dissimilarity index. We then used non-metric multidimensional scaling (NMDS) to represent the dissimilarity of beetle communities between natural forest and plantation areas. These analyzes were conducted using species incidence by sampling site and trapping method to reduce the bias due to varying sampling effort.

Subsequently, using now abundance-based data, we compared family richness between natural forest and plantation areas, by plotting rarefaction and extrapolation curves with the number of collected individuals as a measure of sampling intensity ("iNEXT" package in R v.4.2.1; Hsieh et al., 2016). Additionally, the same procedure was followed to compare family richness obtained by each of the four trapping methods used in our sampling. We also used the asymptotic estimators provided by iNEXT as a measure of the total family richness (including unobserved families) and diversity (Shannon and Simpson indices), in each sampling site, and each trapping method and collection date per site when appropriate. Additionally, we reported sampling completeness by site.

The drivers of the diversity of beetle assemblages were explored with a linear mixed model. The model included the log-transformed asymptotic estimates of richness obtained in each sampling site, trapping method, and collection date, as a dependent variable. Landscape context, trapping methods, temperature, and precipitation were included as independent fixed variables, while random intercepts included the identity of sampling site nested within the province, and the sampling date. The same approach was used for Shannon and Simpson indices as estimates of the diversity of beetle communities, again using asymptotic estimators from iNEXT. For comparison, we also refitted the same models for an estimator of species richness/Shannon diversity/Simpson diversity that was rarefied/extrapolated for a common sample size equal to the median of observed sample sizes across all sites (n = 204).

We used negative binomial generalized linear mixed-effects model to test the effect of the same variables (landscape context, trapping methods, temperature, and precipitation) on beetle abundances recorded in each sampling site at each date and for each trapping method, using again sampling site identify nested within province and sampling date as random intercepts. Negative binomial models were used because a Poisson generalized linear model that we fitted first showed evidence of overdispersion (dispersion parameter = 0.012, p-value = .032), and the negative binomial model had a lower AIC than the Poisson model (ΔAIC = 2208). Here, we
<table><caption>TABLE 1 Characteristics of sampling sites, including sampling method and intensity (number of traps and sampling days), environmental variables (mean annual climate and elevation), location, and landscape context (natural forest vs. plantation).</caption><thead><tr><th>Sites</th><th>Trapping methods</th><th>Number of traps</th><th>Number of sampling days</th><th>Mean annual temperature (°C)</th><th>Total annual precipitation (mm/m²)</th><th>Latitude</th><th>Longitude</th><th>Elevation (m.a.s.l.)</th><th>Landscape context</th><th>Number of sampled individuals</th><th>Number of sampled families</th><th>Sampling completeness</th></tr></thead><tbody><tr><td>Site 1</td><td>Window</td><td>30</td><td>60</td><td>25.85</td><td>1995</td><td>18.12</td><td>102.79</td><td>173</td><td>Plantation</td><td>4461</td><td>30</td><td>0.999</td></tr><tr><td>Site 1</td><td>Pitfall</td><td>30</td><td>60</td><td>25.85</td><td>1995</td><td>18.12</td><td>102.79</td><td>173</td><td>Plantation</td><td>1820</td><td>21</td><td>0.999</td></tr><tr><td>Site 2</td><td>Window</td><td>30</td><td>60</td><td>25.54</td><td>2344</td><td>18.43</td><td>102.55</td><td>190</td><td>Plantation</td><td>3926</td><td>39</td><td>0.999</td></tr><tr><td>Site 3</td><td>Pitfall</td><td>30</td><td>60</td><td>25.54</td><td>2344</td><td>18.43</td><td>102.55</td><td>205</td><td>Natural forest</td><td>6062</td><td>53</td><td>1.000</td></tr><tr><td>Site 4</td><td>Light</td><td>1</td><td>3</td><td>23.15</td><td>1831</td><td>21.19</td><td>101.14</td><td>655</td><td>Plantation</td><td>1023</td><td>19</td><td>0.999</td></tr><tr><td>Site 4</td><td>Hand</td><td>1</td><td>3</td><td>23.15</td><td>1831</td><td>21.19</td><td>101.14</td><td>655</td><td>Plantation</td><td>41</td><td>4</td><td>1.000</td></tr><tr><td>Site 5</td><td>Light</td><td>2</td><td>5</td><td>22.28</td><td>1362</td><td>19.99</td><td>104.62</td><td>368</td><td>Natural forest</td><td>97</td><td>11</td><td>0.969</td></tr><tr><td>Site 6</td><td>Light</td><td>1</td><td>3</td><td>24.52</td><td>1454</td><td>19.75</td><td>102.18</td><td>354</td><td>Natural forest</td><td>142</td><td>3</td><td>1.000</td></tr><tr><td>Site 6</td><td>Hand</td><td>1</td><td>3</td><td>24.52</td><td>1454</td><td>19.75</td><td>102.18</td><td>354</td><td>Natural forest</td><td>45</td><td>7</td><td>0.912</td></tr><tr><td>Site 7</td><td>Window</td><td>15</td><td>9</td><td>26.05</td><td>1347</td><td>17.69</td><td>101.33</td><td>289</td><td>Plantation</td><td>487</td><td>26</td><td>0.988</td></tr><tr><td>Site 8</td><td>Light</td><td>1</td><td>2</td><td>24.32</td><td>2552</td><td>18.4</td><td>103.07</td><td>787</td><td>Natural forest</td><td>784</td><td>39</td><td>0.991</td></tr><tr><td>Site 9</td><td>Hand</td><td>1</td><td>3</td><td>20.49</td><td>2551</td><td>15.24</td><td>106.3</td><td>1232</td><td>Plantation</td><td>23</td><td>2</td><td>1.000</td></tr><tr><td>Site 10</td><td>Hand</td><td>1</td><td>3</td><td>20.26</td><td>2559</td><td>15.14</td><td>106.28</td><td>1313</td><td>Plantation</td><td>93</td><td>2</td><td>1.000</td></tr><tr><td>Site 11</td><td>Hand</td><td>1</td><td>3</td><td>20.35</td><td>2548</td><td>15.18</td><td>106.25</td><td>1290</td><td>Natural forest</td><td>168</td><td>4</td><td>1.000</td></tr><tr><td>Site 12</td><td>Hand</td><td>1</td><td>3</td><td>20.35</td><td>2548</td><td>15.19</td><td>106.27</td><td>1305</td><td>Plantation</td><td>48</td><td>9</td><td>0.939</td></tr></tbody></table>
had to account for the very variable sampling effort that produced the observed variation in beetle abundances; therefore, models also included as an offset the number of traps that was used in each sampling site. The same approach was applied first for the total beetle abundance, then for the most abundant beetle families separately: Carabidae, Scarabaeidae, Nitidulidae, Curculionidae, and Chrysomelidae. All statistical analyses were performed in the R platform (version 4.2.1, R Development Core Team, 2022), mixed models being fitted using the “lme4” package (Bates et al., 2015). Overdispersion and model residuals were checked using the “DHARMa” R package (Hartig, 2022). The significance of each term was estimated using Wald Chi-square tests.

3 | RESULTS

A total of 19,220 beetle individuals were recorded across all sampling sites during the study period, which represents 63 beetle families, including 58 recorded in natural forests and 48 in plantations. All sites had a sampling completeness >0.91 (Table 1). The rarefaction and extrapolation curves, adjusting for the number of specimens collected in natural forests (7339) and plantation areas (11,881), clearly confirmed that beetle community composition in natural for- ests had a higher overall cumulative diversity of families compared to plantations (Figure 2). As expected, sampling conducted using dif- ferent trapping methods did not provide the same number of beetle families, the lowest being achieved by hand collection (Figure S1).

Beetle community composition appeared to be dependent on the trapping method ($F_{3,56}=2.571$, *p*-value < .001), landscape con- text ($F_{1,56}=0.821$, *p*-value < .001), and precipitation ($F_1=0.426$, *p*-value = .026), but not on temperature ($F_{1,56}=0.369$, *p*-value = .0824), according to the PERMANOVA analysis. We noted that the tempera- ture variable becomes significant when this term is included before trapping methods in the PERMANOVA, which may reflect an effect

of temperature on community composition. Similarly, visual exam-
ination of NMDS confirmed that the beetle families collected using
different trapping methods, or collected either in natural forests
or in plantation areas, differed but were also largely overlapping
(Figure 3).

The richness of beetle communities (in terms of the estimated number of families per site, per date and per trapping method) was found to be dependent on the trapping method (Chi-square=8.947, df=3, p-value=.030; Figure S2). Specifically, the highest richness was sampled in pitfall traps (marginal mean family richness=31.8; CI 95%=[1.473, 687.6]), while window traps provided the lowest beetle family richness (marginal mean family richness=10.0; CI 95%=[0.404, 248.9]). It was not, however, influenced by the other three variables of interest (landscape context (Chi-square=0.098, df=1, p-value=.754), temperature (Chi-square=0.293, df=1, p-value=.588), precipitation (Chi-square=0.218, df=1, p-value=.640)). In addition, there was also no significant effect of any variable on beetle community diversity as estimated by Shannon index, but it appeared that Simpson diversity was slightly higher in natural forests compared to plantations (Chi-square=4.091, df=1, p-value=.043; marginal mean=3.95; CI 95%=[1.890, 8.230] vs. 2.48; CI 95%=[1.160, 5.310], respectively; Table S1 and Figure S2). There was no significant effect of any of the variables tested when richness and diversity indices were estimated for a common sample size.

There was a strongly significant effect of the trapping method (Chi-square=40.319, df=3, *p*-value < .001) on the variation of beetle abundance across sites, with light traps being predicted to provide higher abundance (Figure 4). Moreover, the abundance of beetle communities was significantly higher (Chi-square=3.892, df=1, *p*-value = .049) in natural forests (estimated marginal mean abundance=60.0; CI 95%=[31.3, 114.7]) compared to plantations (estimated marginal mean abundance=29.8; CI 95%=[14.2, 62.7]). It was not, however, associated with temperature (Chi-square = 0.047, df=1, *p*-value = 0.829) or precipitation (Chi-square = 0.860, df=1, *p*-value = .354) (Figure 4). When analyzes were computed for five families separately, we similarly identified a significant effect of trap- ping methods. We were not, however, able to distinguish the effect of landscape context on beetle abundance, except for Scarabaeidae (Table S2) where the predicted abundance were actually higher in plantations (estimated marginal mean abundance = 6.735; CI 95%=[1.959, 23.155]) than in natural forests (estimated marginal mean abundance = 1.007; CI 95%=[0.555, 5.410]).

<figure><img src="image_4.png" /><figcaption>FIGURE 2 Rarefaction and extrapolation curves of beetle family richness for natural forests and plantation areas.</figcaption></figure>

4 | DISCUSSION

Despite the fact that tropical areas host the majority of insect diver-
sity, they are generally vastly undersampled compared to temper-
ate ecosystems (Hellmann & Sanders, 2007). This knowledge gap is
particularly undesirable in Southeastern Asia, which is experiencing
a rapid land use change that may imperil its rich insect biodiver-
sity. In this study, we report for the first time a large sampling of
<figure><img src="image_5.png" /><figcaption>FIGURE 3 Non-metric multidimensional scaling (NMDS) representing the pairwise dissimilarity of beetle communities (computed as the Bray-Curtis distance). Colored dots show the position of each combination of sampling site, date of collection and trapping method in the NMDS ordination, grouped into natural forest and plantation areas.</figcaption></figure>

the whole order Coleoptera carried out across a large part of Laos, which enabled the investigation of the drivers of beetle diversity at a large spatial and taxonomical scale. This is a step toward a better understanding of insect diversity in Asian tropical forests, and of the threats they may face.

The present study first confirmed our expectation that beetle diversity is high in the country. From the ca. 20,000 collected specimens, we recorded 63 beetle families. Moreover, more than 50 beetle specimens remained impossible to assign to any known family and were thus excluded from this study. Some of them may belong to additional beetle families, suggesting that the diversity of beetle families in the region may be even larger. It is however hard to reliably estimate the actual number of families living across the various land cover types of the country without a meticulous approach involving standardized sampling protocols and molecular taxonomy (García-Robledo et al., 2020). The data presented here reveals slightly higher number of families than that of earlier recorded data in surrounding countries such as Thailand, Vietnam, and Hong Kong (Rattanawannee et al., 2013; Thinh et al., 2004;

Zhao et al., 2022). Overall, 177 Coleoptera families were recorded globally, meaning that we recorded 36% of all known beetle families (Moodley et al., 2022). However, due to practical constraints, the present study described beetle communities at the family level only. Indeed, beetle taxonomy is notoriously difficult in the absence of detailed identification keys or molecular tools (see e.g., Jin et al., 2020; Sabatelli et al., 2021). The actual diversity at the species level, that is species richness, is thus much more important, and may include endemic or undescribed species.

According to our results, beetle abundance was reduced in plantations compared to natural forests, which is in line with global patterns which show that insects are sensitive to habitat disturbance from human activities such as agricultural expansions or settlements (Hansen et al., 2012; New et al., 2021; Sánchez-Bayo & Wyckhuys, 2019). Natural forest ecosystems play an important role in species diversity worldwide (Gibson et al., 2011), whereas intensification of agriculture is identified as a major cause of insect diversity decline and extinction (Sánchez-Bayo & Wyckhuys, 2019), but also soil carbon loss (Guo et al., 2022).
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

FIGURE 4 Partial dependence plots showing the predicted beetle abundance for different levels of the four response variables tested (only trapping method and landscape context are significant at the $α=.05$ level), modeled using a negative binomial mixed-effect model.

A strong impact of land use intensification has been reported for beetles and other insects in tropical forests in Asia, as well as in Africa and America (Phillips et al., 2017). Several observational and experimental studies have revealed that the conversion of natural forests into plantations is harmful to species that cannot adapt to their new environmental conditions (Uribe et al., 2021; Warren-Thomas et al., 2015); our results show that Laos is no exception. Therefore, a large number of insect communities may be, currently or in the near future, at risk from land use intensification in Laos, even though the region is still mainly covered by mountains and forests.

A more in-depth understanding of the response of beetle diversity to current and future economic development in Laos is needed to implement practical conservation actions. In this regard, this study clearly demonstrates that the rapid and continuous land use changes in the country is experiencing may threaten beetle communities; not only their abundance declined in plantations, but across the whole survey we also sampled fewer families in plantations compared to natural forests (48 vs. 58 observed families, respectively). This finding is in line with other studies showing the impact human-modified landscapes can have on beetle biodiversity; for example, dung beetle communities are well recognized as a good indicator to estimate the influence of anthropogenic habitats in tropical forests (Gardner et al., 2008; Halffter & Arellano, 2002). In the present study, we investigated beetle communities across a long latitudinal gradient, in which the northern to center parts are facing a modernization of the road network in addition to the conversion of forests into plantations. In this regard, a recent report by Danyo et al. (2018) pointed out the potential risks of forest and biodiversity loss resulting from road improvement in Laos. We believe that our findings will

therefore be useful and important in order to predict the conservation issues arising from land use changes in the region.

The observed difference in richness and abundance between plantations and natural forests were associated with differences in terms of family composition, although it was difficult to detect families that were strongly associated with a type of landscape. Similarly, although fewer families were sampled in plantations at the scale of the whole country, anthropization seems to have relatively little impact on community diversity when analyzed at the scale of each sampling site. Indeed, we did not detect difference in family richness, but only on Simpson diversity index (which confirms the effect on abundance). Logically, this must be caused by a higher homogeneity of communities located in plantations, that is, the same set of families are found in all plantations, while forest communities are more diverse from each other (despite a similar local diversity). Biotic homogenization, in which a few common species takeover specialist species, is frequently observed in human-modified landscapes (McKinney & Lockwood, 1999), and has also been observed in beetles (Ramirez-Ponce et al., 2019). However, results may vary at a lower taxonomic level; it is likely that species diversity is actually reduced along human activities, in accordance with previous studies (e.g., Jung & Lee, 2016; Vanbergen et al., 2005). Here, the beetle collection contained a huge number of specimens, and the taxonomy of many groups is challenging due to their complex diagnostic morphological characters and their small body size. An improved dataset that would distinguish individuals at the species-level may reveal a slightly picture, including perhaps an effect of agricultural development on species diversity at the local scale.

Due to the geographical scale of the study, we found an effect of climate, specifically precipitation, on community composition. It was
not, however, reflected by differences in terms of family richness or abundance along the climate gradient. However, again, it could be because the present study has investigated families only and not species. Generally, assessments of climatic niches are considerably more precise when carried out at lower taxonomic levels (Bayliss et al., 2022; Gonzalez et al., 2011), meaning that each individual species may have vastly different responses to climatic variables that would be masked when merged into whole families.

It was clear that both composition and abundance differed depending on the trapping method. Here, we showed that, when all other variables are taken into account, pitfall traps were able to capture a much larger number of beetle families compared to all other methods. However, while light traps were way behind in terms of diversity, they provided a larger sample (in terms of abundance). In this case, a good sampling strategy should probably use several complementary approaches to sample the whole diversity of beetles. For instance, while light traps are probably suitable for sampling flying insects, pitfall traps are adapted for ground-dwelling species.

## 5 | CONCLUSIONS

The present study provides the first approach that attempts at investigating the effect of various independent variables on beetle community composition and richness along contrasted landscape types across Laos. A potential limitation of the study is that sampling did not follow a standardized protocol across the country. On the contrary, we compiled here a dataset of beetle specimens collected using different sampling efforts and methods. We employed various approaches to account for these unequal sampling strategies: country scale and local family richness were estimated from accumulation curves, trapping methods were always included as covariates, and sampling effort was incorporated as an offset in statistical models. The main outcome of our results was that the conversion of natural forests to plantations appeared to be harmful for beetle communities, since less families were found in plantations compared to natural forests overall, and their abundance was reduced locally.

We may gain additional knowledge by comparing beetle communities in pairs of sites that were sampled in similar conditions. Here, site 9 (rubber plantation) and site 11 (natural forest) were located in the same province of Champasak, and were sampled in the same year and using the same method of hand collection. Results of the sampling confirmed our general conclusions, as we found a reduced number of families (2 vs. 4) and abundance (23 vs. 168) in the plantation area compared to the nearby natural forest. More studies are needed to better understand this pattern. This may be achieved in the future by implementing long-term monitoring of beetles across Laos following a simple protocol that can be used by many volunteers, and by incorporating a better taxonomic resolution (i.e., species) in analyzes. The central part of the country, which we did not sample in this study, would also benefit from a proper inventory of beetles. In this regard, our study shed light on the potential usefulness and difficulties of sampling beetle diversity in southern Asia

and in Laos in particular. It is widely accepted that Laos is a hotspot for the biodiversity of beetles, insects, and other organisms. By providing a first large-scale view of beetle family diversity of Laos, we aim to provide a basis for future studies investigating the impact of an extremely rapid economic change, associated with land use change, on insect diversity in Laos.

## AUTHOR CONTRIBUTIONS

**Bounsanong Chouangthavy:** Conceptualization (equal); data curation (equal); formal analysis (equal); methodology (equal); visualization (equal); writing – original draft (equal); writing – review and editing (equal). **Yaan Fourcade:** Conceptualization (equal); data curation; formal analysis (equal); methodology; software (equal); visualization (equal); writing – review and editing (equal).

## ACKNOWLEDGMENTS

The authors thank the students of the Plant Protection Unit at NUoL for their help with field collections and for working on the specimens in the laboratory. We are thankful to Dr. Christian Hartmann for arranging the first author's hosting at IEES Paris, Université Paris-Est Créteil (UPEC), Institut de Recherche pour le Développement (IRD), and French National Center for Scientific Research/Sorbonne Université from May to August 2022. We also thank Dr. Phoutthasone Sibounnavong, Vice Head Department of Plant Science, Faculty of Agriculture, NUoL for providing the necessary facilities. Thanks are extended to Dr. Naoise Nunan, Dr. Thomas Lerch, Dr. Virginie Roy, and Dr. Lise Dupont for their kind care during the first author's stay in France, to Mrs. Laurence Goury for providing a training course, to Mr. Thierry Pilorge for his French language course, and to Mr. Johnny Louangphan for his help. This research was partly funded by an Erasmus+ grant funded through UPEC and IRD, which supported the first author's visit to France.

## CONFLICT OF INTEREST STATEMENT

The authors do not have any conflicts of interest.

## DATA AVAILABILITY STATEMENT

The beetle sampling data is available in the Figshare repository: https://doi.org/10.6084/m9.figshare.23428457.v1.

## ORCID

Bounsanong Chouangthavy https://orcid.org/0000-0003-1715-6434

Yaan Fourcade https://orcid.org/0000-0003-3820-946X

## REFERENCES

Andrew, N. R., & Hughes, L. (2004). Species diversity and structure of phytophagous beetle assemblages along a latitudinal gradient: Predicting the potential impacts of climate change. *Ecological Entomology*, *29*, 527-542. https://doi.org/10.1111/j.0307-6946.2004.00639.x

Báldi, A. (2003). Using higher taxa as surrogates of species richness: A study based on 3700 Coleoptera, Diptera, and Acari species in central-Hungarian reserves. *Basic and Applied Ecology*, *4*, 589-593. https://doi.org/10.1078/1439-1791-00193
Bates, D., Machler, M., Boler, B. M., & Walker, S. C. (2015). Fitting linear mixed-effects models using lme4. *Journal of Statistical Software*, **67**, 1–48. https://doi.org/10.18637/jss.v067.i01

Bayliss, S. L. J., Papeş, M., Schweitzer, J. A., & Bailey, J. K. (2022). Aggregate population-level models informed by genetics predict more suitable habitat than traditional species-level model across the range of a widespread riparian tree. *PLoS One*, 17, e0274892. https://doi.org/10.1371/journal.pone.0274892

Borda-de Água, L., Barrientos, R., Beja, P., & Pereira, H. M. (2017). *Railway ecology* (p. 320). Springer Nature.

Ceballos, G., & Ehrlich, P. R. (2006). Global mammal distributions, biodiversity hotspots, and conservation. *Proceedings of the National Academy of Sciences of the United States of America*, **103**, 19374–19379. https://doi.org/10.1073/pnas.0609334103

Chen, Z., & Haynes, K. E. (2017). Impact of high-speed rail on regional economic disparity in China. *Journal of Transport Geography*, **65**, 80–91. https://doi.org/10.1016/j.jtrangeo.2017.08.003

Chouangthavy, B., Bouttavong, K., Louangphan, J., Phewphanh, P., Sibounnavong, P., Souksavat, S., & Babbreier, D. (2020). Beetle biodiversity in forest habitats in Laos depends on the level of human exploitation. *Journal of Insect Conservation*, **24**, 833–840. https://doi.org/10.1007/s10841-020-00255-x

Correa-Carmona, Y., Rougerie, R., Arnal, P., Ballesteros-Mejia, L., Beck, J., Dolédec, S., & Decaëns, T. (2022). Functional and taxonomic responses of tropical moth communities to deforestation. *Insect Conservation and Diversity*, **15**, 236–247. https://doi.org/10.1111/icad.12549

Danyo, S., Dasgupta, S., & Wheeler, D. (2018). Potential forest loss and biodiversity risks from road improvement in Lao PDR. World Bank Policy Research Working Paper, (8569).

Dolson, S. J., Loewen, E., Jones, K., Jacobs, S. R., Solis, A., Hallwachs, W., Brunke, A. J., Janzen, D. H., & Smith, M. A. (2021). Diversity and phylogenetic community structure across elevation during climate change in a family of Hyperdiverse neotropical beetles (Staphylinidae). *Ecography*, **44**, 740–752. https://doi.org/10.1111/ecog.05427

Eggleton, P. (2020). The state of the world's insects. *Annual Review of Environment and Resources*, **45**, 7–20. https://doi.org/10.1146/annurev-environ-012420-050035

Essenwanger, O. M. (2001). Classification of climates, in world survey of climatology 1C (p. 102). Elsevier Science.

Estoque, R. C., Ooba, M., Avitabile, V., Hijioka, Y., DasGupta, R., Togawa, T., & Murayama, Y. (2019). The future of Southeast Asia's forests. *Nature Communications*, **10**, 1829. https://doi.org/10.1038/s41467-019-09646-4

Fick, S. E., & Hijmans, R. J. (2017). WorldClim 2: New 1 km spatial resolution climate surfaces for global land areas. *International Journal of Climatology*, **37**, 4302–4315. https://doi.org/10.1002/joc.5086

García-Robledo, C., Kuprewicz, E. K., Baer, C. S., Clifton, E., Hernández, G. G., & Wagner, D. L. (2020). The Erwin equation of biodiversity: From little steps to quantum leaps in the discovery of tropical insect diversity. *Biotropica*, **52**, 590–597. https://doi.org/10.1111/btp.12811

Gardner, T. A., Hernández, M. I., Barlow, J., & Peres, C. A. (2008). Understanding the biodiversity consequences of habitat change: The value of secondary and plantation forests for neotropical dung beetles. *Journal of Applied Ecology*, **45**, 883–893. https://doi.org/10.1111/j.1365-2664.2008.01454.x

Gebert, F., Steffan-Dewenter, I., Moretto, P., & Peters, M. K. (2020). Climate rather than dung resources predict dung beetle abundance and diversity along elevational and land use gradients on Mt. Kilimanjaro. *Journal of Biogeography*, **47**, 371–381. https://doi.org/10.1111/jbi.13710

Ghannem, S., Touaylia, S., & Boumaiza, M. (2018). Beetles (Insecta: Coleoptera) as bioindicators of the assessment of environmental pollution. *Human and Ecological Risk Assessment: An International Journal*, **24**, 456–464. https://doi.org/10.1080/10807039.2017.1385387

Gibson, L., Lee, T. M., Koh, L. P., Brook, B. W., Gardner, T. A., Barlow, J., Peres, C. A., Bradshaw, C. J. A., Laurance, W. F., Lovejoy, T. E., & Sodhi, N. S. (2011). Primary forests are irreplaceable for sustaining tropical biodiversity. *Nature*, **478**(7369), 378–381. https://doi.org/10.1038/nature10425

González, E., Salvo, A., & Valladares, G. (2015). Arthropods on plants in a fragmented neotropical dry forest: A functional analysis of area loss and edge effects. *Insect Sci.*, **22**, 129–138. https://doi.org/10.1111/1744-7917.12107

Gonzalez, S. C., Soto-Centeno, J. A., & Reed, D. L. (2011). Population distribution models: Species distributions are better modeled using biologically relevant data partitions. BMC Ecology, **11**, 20. https://doi.org/10.1186/1472-6785-11-20

Guo, J., Feng, H., McNie, P., Wang, W., Peng, C., Feng, L., & Yu, Y. (2022). The effect of the conversion from natural broadleaved forests into Chinese fir (*Cunninghamia lanceolata* (Lamb.) Hook.) plantations on soil microbial communities and nitrogen functional genes. *Forests*, **13**, 158. https://doi.org/10.3390/f13020158

Halffter, G., & Arellano, L. (2002). Response of dung beetle diversity to human-induced changes in a tropical landscape. *Biotropica*, **34**, 144–154. https://doi.org/10.1111/j.1744-7429.2002.tb00250.x

Hansen, A. J., DeFries, R. S., & Turner, W. (2012). Land use change and biodiversity. In G. Gutman, A. C. Janetos, C. O. Justice, E. F. Moran, J. F. Mustard, R. R. Rindfuss, D. Skole, B. L. Turner II, & M. A. Cochrane (Eds.), *Land change science* (pp. 277–299). Springer.

Harris, J. E., Rodenhouse, N. L., & Holmes, R. T. (2019). Decline in beetle abundance and diversity in an intact temperate forest linked to climate warming. *Biological Conservation*, **240**, 108219. https://doi.org/10.1016/j.biocon.2019.108219

Hartig, F. (2022). DHARMA: Residual diagnostics for hierarchical (Multi-Level/Mixed) regression models. R package version 0.4.6. https://CRAN.R-project.org/package=DHARMA

Hellmann, J. J., & Sanders, N. J. (2007). Chapter 2: The extent and future of global insect diversity. In R. E. Hester & R. M. Harrison (Eds.), *Biodiversity under threat* (pp. 33–55). RSC Publishing. https://doi.org/10.1039/9781847557650-00033

Hortal, J., Diniz-Filho, J. A. F., Bini, L. M., Rodríguez, M. Á., Baselga, A., Nogués-Bravo, D., Rangel, T. F., Hawkins, B. A., & Lobo, J. M. (2011). Ice age climate, evolutionary constraints and diversity patterns of European dung beetles. *Ecology Letters*, **14**, 741–748. https://doi.org/10.1111/j.1461-0248.2011.01634.x

Hsieh, T. C., Ma, K. H., & Chao, A. (2016). iNEXT: An R package for rarefaction and extrapolation of species diversity (Hill numbers). *Methods in Ecology and Evolution*, **7**, 1451–1456. https://doi.org/10.1111/2041-210X.12613

Jin, M., Zwick, A., Ślipiński, A., de Keyzer, R., & Pang, H. (2020). Museomics reveals extensive cryptic diversity of Australian prionine longhorn beetles with implications for their classification and conservation. *Systematic Entomology*, **45**, 745–770. https://doi.org/10.1111/syen.12424

Jung, J. K., & Lee, J. H. (2016). Forest-farm edge effects on communities of ground beetles (Coleoptera: Carabidae) under different landscape structures. *Ecological Research*, **31**, 799–810. https://doi.org/10.1007/s11284-016-1388-1

Kumar, P., Gale, S. W., Bouamanivong, S., & Fischer, G. A. (2016). Identifying orchid hotspots for biodiversity conservation in Laos: The limestone karst vegetation of Vang Vieng District, Vientiane Province. *Journal of Threatened Taxa*, **8**, 9397–9417. https://doi.org/10.11609/jott.2826.8.12.9397-9417

Kusakabe, K., & Chanthoumphone, C. (2021). Transition from subsistence agriculture to rubber plantations in northern Laos: Analysis of household livelihood strategies by ethnicity and gender. SAGE Open, **11**, 1–13. https://doi.org/10.1177/21582440211011463
Kusuma, Y. W. C., Rembold, K., Tjitrosoedirdjo, S. S., & Kreft, H. (2018). Tropical rainforest conversion and land use intensification reduce understorey plant phylogenetic diversity. *Journal of Applied Ecology*, 55, 2216-2226. https://doi.org/10.1111/1365-2664.13201

MAF, & STEA. (2003). *Biodiversity country report*. Ministry of Agriculture and Forestry (MAF).

McKinney, M. L., & Lockwood, J. L. (1999). Biotic homogenization: A few winners replacing many losers in the next mass extinction. *Trends in ecology & evolution*, 14, 450-453. https://doi.org/10.1016/S0169-5347(99)01679-1

Ministry of Natural Resources and the Environment (MoNRE). (2016). *National biodiversity strategy and action plan for Lao PDR 2016-2025*. MoNRE.

Moodley, S., Procheş, Ş., Perera, S. J., Lubbe, E., Ramdhani, S., & Leschen, R. A. (2022). Analysis of the diversity and distributional patterns of coleopteran families on a global scale. *Zootaxa*, 5138, 575-583. https://doi.org/10.11646/zootaxa.5138.5.5

Myers, N., Mittermeier, R. A., Mittermeier, C. G., Da Fonseca, G. A., & Kent, J. (2000). Biodiversity hotspots for conservation priorities. *Nature*, 403, 853-858. https://doi.org/10.1038/35002501

New, T. R., Sands, D. P., & Taylor, G. S. (2021). Roles of roadside vegetation in insect conservation in Australia. *Austral Entomology*, 60, 128-137. https://doi.org/10.1111/aen.12511

Ng, L. S., Campos-Arceiz, A., Sloan, S., Hughes, A. C., Tiang, D. C. F., Li, B. V., & Lechner, A. M. (2020). The scale of biodiversity impacts of the belt and road initiative in Southeast Asia. *Biological Conservation*, 248, 108691. https://doi.org/10.1016/j.biocon.2020.108691

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., Simpson, G. L., Solymos, P., Stevens, M. H. H., & Wagner, H. (2022). *Vegan: Community ecology package*. https://CRAN.R-proje ct.org/package=vegan

Parikh, G., Rawtani, D., & Khatri, N. (2021). Insects as an indicator for environmental pollution. *Environmental Claims Journal*, 33, 161-181. https://doi.org/10.1080/10406026.2020.1780698

Phillips, H. R., Newbold, T., & Purvis, A. (2017). Land-use effects on local biodiversity in tropical forests vary between continents. *Biodiversity and Conservation*, 26, 2251-2270. https://doi. org/10.1007/s10531-017-1356-2

R Core Team. (2022). R: A language and environment for statistical computing. Version 4.2.1. R Foundation for Statistical Computing. https:// www.R-project.org/

Ramírez-Ponce, A., Calderón-Patrón, J. M., Vásquez, H. M. G., & Moreno, C. E. (2019). Biotic heterogeneity among scarab beetle communities in an anthropized landscape in the Central Valleys of Oaxaca, Mexico. *Journal of Insect Conservation*, 23, 765-776. https://doi. org/10.1007/s10841-019-00169-3

Rattanawannee, A., Duangpukdee, O., & Poolprasert, P. (2013). Insect diversity during different stages of Asiatic elephant dung deterioration in eastern Thailand. *Agriculture and Natural Resources*, 47, 387-397.

Sabatelli, S., Ruspantini, P., Cardoli, P., & Audisio, P. (2021). Underestimated diversity: Cryptic species and phylogenetic relationships in the subgenus *Cobalius* (Coleoptera: Hydraenidae) from marine rockpools. *Molecular Phylogenetics and Evolution*, 163, 107243. https://doi.org/10.1016/j.ympev.2021.107243

Sánchez-Bayo, F., & Wyckhuys, K. A. G. (2019). Worldwide decline of the entomofauna: A review of its drivers. *Biological Conservation*, 232, 8-27. https://doi.org/10.1016/j.biocon.2019.01.020

Sekerka, L., & Geiser, M. (2016). Sagrinae of Laos (Coleoptera: Chrysomelidae). *Entomologica Basiliensia et Collectionis Frey*, 35, 443-453.

Slade, E. M., Mann, D. J., & Lewis, O. T. (2011). Biodiversity and ecosystem function of tropical forest dung beetles under contrasting logging regimes. *Biological Conservation*, 144, 166-174. https://doi.org/10.1016/j.biocon.2010.08.011

Smith, H., Lu, J., To, P. X., Mienmany, S., & Soukphaxay, K. (2020). Rubber plantation value chains in Laos: Opportunities and constraints in policy, legality and wood processing, report produced for ACIAR project FST/2016/151 - Advancing enhanced wood manufacturing industries in Laos and Australia and Forest Trends.

Spector, S. (2006). Scarabaeine dung beetles (Coleoptera: Scarabaeidae: Scarabaeinae): An invertebrate focal taxon for biodiversity research and conservation. *The Coleopterists Bulletin*, 60, 71-83. https://doi.org/10.1649/0010-065X(2006)60[71:SDBCSS]2.0.CO;2

The World Bank. (2017). *Lao People's Democratic Republic: First programmatic green growth development policy operation*. The World Bank.

Thinh, T. H., Tru, H. V., Du, T. T., & Thai, P. H. (2004). The insect diversity at some national parks and nature reserves of Vietnam. *Academia Journal of Biology*, 24, 1-12. https://doi.org/10.15625/0866-7160/v24n4.6905

Torres, A., Jaeger, J. A., & Alonso, J. C. (2016). Assessing large-scale wildlife responses to human infrastructure development. *Proceedings of the National Academy of Sciences of the United States of America*, 113, 8472-8477. https://doi.org/10.1073/pnas.1522488113

Uribe, S. V., Garcia, N., & Estades, C. F. (2021). Effect of land use history on biodiversity of pine plantations. *Frontiers in Ecology and Evolution*, 9, 430. https://doi.org/10.3389/fevo.2021.609627

Vanbergen, A. J., Woodcock, B. A., Watt, A. D., & Niemelä, J. (2005). Effect of land-use heterogeneity on carabid communities at the landscape scale. *Ecography*, 28, 3-16. https://doi.org/10.1111/j.0906-7590.2005.03991.x

Wagner, D. L. (2020). Insect declines in the Anthropocene. *Annual Review of Entomology*, 65, 457-480. https://doi.org/10.1146/annur ev-ento-011019-025151

Warren-Thomas, E., Dolman, P. M., & Edwards, D. P. (2015). Increasing demand for natural rubber necessitates a robust sustainability initiative to mitigate impacts on tropical biodiversity. *Conservation Letters*, 8, 230-241. https://doi.org/10.1111/conl.12170

Zhao, S., Tong, Y., Teng, B., Chen, X., Yang, X., Li, J., & Bai, M. (2022). A species diversity dataset of beetles by three passive acquisition methods in Tei Tong Tsai (Hong Kong). *Scientific Data*, 9, 1-7. https://doi.org/10.1038/s41597-022-01310-9

Zödl, B., & Wittmann, K. J. (2003). Effects of sampling, preparation and defecation on metal concentrations in selected invertebrates at urban sites. *Chemosphere*, 52, 1095-1103. https://doi.org/10.1016/S0045-6535(03)00442-9

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Chouangthavy, B., & Fourcade, Y. (2023). Large-scale sampling of beetle communities in Laos shows that conversion of natural forests into plantations leads to a decline in family richness and abundance. *Ecology and Evolution*, 13, e10258. https://doi.org/10.1002/ ece3.10258
Supporting Information

<table><caption>TABLE S1 Effect of trapping method, landscape context (natural forest vs. plantation), temperature and precipitation on Shannon or Simpson diversity indices of beetle communities.</caption><thead><tr><th>Variables</th><th>LR Chi²</th><th>df</th><th>p-value</th></tr></thead><tbody><tr><td><em>Shannon</em></td><td></td><td></td><td></td></tr><tr><td>Trapping method</td><td>1.2619</td><td>3</td><td>0.7382</td></tr><tr><td>Landscape context</td><td>0.3129</td><td>1</td><td>0.5759</td></tr><tr><td>Temperature</td><td>0.1302</td><td>1</td><td>0.7182</td></tr><tr><td>Precipitation</td><td>0.0594</td><td>1</td><td>0.8075</td></tr><tr><td><em>Simpson</em></td><td></td><td></td><td></td></tr><tr><td>Trapping method</td><td>1.5638</td><td>3</td><td>0.6676</td></tr><tr><td>Landscape context</td><td>4.0912</td><td>1</td><td>0.0431</td></tr><tr><td>Temperature</td><td>0.0014</td><td>1</td><td>0.9699</td></tr><tr><td>Precipitation</td><td>0.3368</td><td>1</td><td>0.5617</td></tr></tbody></table>
<table><caption>**TABLE S2** Effect of trapping method, landscape context (natural forest vs. plantation), temperature and precipitation on the abundance of five beetle families, which included Carabidae, Scarabaeidae, Nitidulidae, Curculionidae, and Chrysomelidae.</caption><thead><tr><th>Family</th><th>Variables</th><th>LR Chi²</th><th>df</th><th>p-value</th></tr></thead><tbody><tr><td rowspan="4">Carabidae</td><td>Trapping methods</td><td>42.3270</td><td>3</td><td>&lt;0.0001</td></tr><tr><td>Landscape context</td><td>1.2381</td><td>1</td><td>0.2658</td></tr><tr><td>Temperature</td><td>2.4874</td><td>1</td><td>0.1148</td></tr><tr><td>Precipitation</td><td>0.0001</td><td>1</td><td>0.9907</td></tr><tr><td rowspan="4">Scarabaeidae</td><td>Trapping methods</td><td>30.0213</td><td>3</td><td>&lt;0.0001</td></tr><tr><td>Landscape context</td><td>5.7986</td><td>1</td><td>0.0160</td></tr><tr><td>Temperature</td><td>0.3972</td><td>1</td><td>0.5286</td></tr><tr><td>Precipitation</td><td>0.2440</td><td>1</td><td>0.6213</td></tr><tr><td rowspan="4">Nitidulidae</td><td>Trapping methods</td><td>47.4962</td><td>3</td><td>&lt;0.0001</td></tr><tr><td>Landscape context</td><td>1.2096</td><td>1</td><td>0.2714</td></tr><tr><td>Temperature</td><td>0.1611</td><td>1</td><td>0.6882</td></tr><tr><td>Precipitation</td><td>0.0775</td><td>1</td><td>0.7806</td></tr><tr><td rowspan="4">Curculionidae</td><td>Trapping methods</td><td>27.7715</td><td>3</td><td>&lt;0.0001</td></tr><tr><td>Landscape context</td><td>0.7580</td><td>1</td><td>0.3839</td></tr><tr><td>Temperature</td><td>3.5495</td><td>1</td><td>0.0596</td></tr><tr><td>Precipitation</td><td>4.0963</td><td>1</td><td>0.0430</td></tr><tr><td rowspan="4">Chrysomelidae</td><td>Trapping methods</td><td>43.8536</td><td>3</td><td>&lt;0.0001</td></tr><tr><td>Landscape context</td><td>0.8676</td><td>1</td><td>0.3516</td></tr><tr><td>Temperature</td><td>3.0045</td><td>1</td><td>0.0830</td></tr><tr><td>Precipitation</td><td>0.0961</td><td>1</td><td>0.7566</td></tr></tbody></table>
<figure><img src="image_10.png" /></figure>

**FIGURE S1:** Rarefaction and extrapolation curves of beetle family richness for each of the four trapping methods used in our sampling.
<figure><img src="image_11.png" /></figure>

**FIGURE S2:** Partial dependence plots showing the predicted beetle richness (A), Simpson diversity (B) and Shannon diversity (C) for different levels of the four response variables tested. Only the effects of trapping method on richness and of landscape context on Simpson diversity are significant at the α=0.05 level.