Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 22 August 2025

Sec. Terrestrial Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1570294

This article is part of the Research TopicForest Microbiome: Dynamics and Interactions in the Anthropocene EraView all 21 articles

Soil microbial community responses to active and passive restoration of selectively logged Bornean tropical forest

  • 1UK Centre for Ecology and Hydrology, Lancaster, United Kingdom
  • 2Lancaster Environment Centre, Lancaster University, Lancaster, United Kingdom
  • 3UK Centre for Ecology and Hydrology, Wallingford, United Kingdom
  • 4School of Natural Sciences, Bangor University, Bangor, United Kingdom
  • 5Sabah Forestry Department, Forest Research Centre, Sandakan, Malaysia
  • 6Faculty of Tropical Forestry, Universiti Malaysia Sabah, Kota Kinabalu, Malaysia

Tropical rainforests support critical biogeochemical cycles regulated by complex plant-soil microbial interactions but are threatened by global change. Much of the uniquely biodiverse and carbon rich forest on Borneo has been lost through extensive conversion to monoculture plantation, and a significant proportion of the remaining forest has been heavily modified by selective logging. Ecological restoration of tropical forest aims to return forests to a near pristine state, but restoration initiatives are hindered by limited understanding of the underpinning plant-soil feedbacks, and impacts on soil microbial communities are unresolved. We characterized soil properties and soil bacterial and fungal communities using amplicon sequencing across adjacent old-growth and selectively logged lowland dipterocarp forest in Borneo undergoing either natural regeneration or restoration by enrichment planting. While many soil properties were similar across forest types, we found contrasting responses of different soil microbial groups to active and passive restoration. Bacterial and fungal community composition were generally distinct in old-growth forest and more similar in logged forest. Bacterial alpha diversity and rate of spatial turnover appeared to recover toward old-growth forest with active restoration, while fungal alpha diversity showed slower signs of recovery. The composition and rate of spatial turnover in mycorrhizal communities was most different between old-growth and actively restored forest, possibly resulting from mycorrhizal associations of tree species planted during restoration. Surprisingly, old-growth forest shared fewer microbial taxa with actively restored forest than with naturally regenerating forest, suggesting current restoration practices (removal of lianas and understorey vegetation) may be selecting for different microbial communities. Taken together, our findings show that certain attributes of key soil microbial groups remain distinct from old-growth forest almost two decades after logging disturbance, and some may diverge with active restoration. Changes in enrichment planting practices to promote rehabilitation of belowground communities may be required for successful biodiversity conservation and recovery of vital ecosystem functions.

1 Introduction

Old-growth forests are rapidly being replaced by human-modified secondary forest worldwide, with highest conversion rates in the tropics (Keenan et al., 2015), affecting the crucial role of tropical forests as global reservoirs of biodiversity and carbon (C) (Baccini et al., 2017; Myers et al., 2000; Pan et al., 2011; Powers and Jetz, 2019; Qie et al., 2017). The forests of Borneo are a hotspot of forest disturbance and loss, driven by timber extraction and conversion to oil palm plantation. A reduction of forest cover by more than 30% since the early 1970s means increasing pressure on the remaining forest to provide vital ecosystem services. However, almost half (16.8 Mha) of the remaining forest has been heavily modified through selective logging for the commercially valuable and canopy-dominant dipterocarp trees (Appanah and Turnbull, 1998; Gaveau et al., 2014; Gaveau et al., 2016). Their targeted removal and creation of landings and skid trails affects forest structure, plant and soil microbial communities, ecosystem C cycling and C sequestration (Asner et al., 2018; Both et al., 2019; Ellis et al., 2016; Marsh et al., 2025; Riutta et al., 2018; Robinson et al., 2020; 2024). This practice reflects regional and global trends, as selective logging is widespread across Southeast Asia (Stibig et al., 2014) and the primary driver of tropical forest modification, affecting more than half of tropical forest worldwide (Potapov et al., 2017; Potapov et al., 2008).

Forest disturbance can significantly affect soil physicochemical properties, often reducing soil C pools (Don et al., 2011; Wei et al., 2014), altering nutrient availability and overall fertility (Daljit Singh et al., 2013; Paul et al., 2010). Changes in soil properties can impact on soil microbial communities (Jesus et al., 2009; Tripathi et al., 2012) by determining availability of effective resources and creation of different ecological niches (Zhang et al., 2018). Soil microbial diversity may increase with a certain level of disturbance (Ferrenberg et al., 2013; Galand et al., 2016; Wilkinson, 1999; Zhang et al., 2011), while changes in community structure and function may directly relate to disturbance intensity (Berga et al., 2012). Soil fungal communities, particularly mycorrhizae, have been shown to be highly sensitive to logging disturbance, likely as dipterocarp trees are an ectomycorrhizal-associating species (Kerfahi et al., 2014; McGuire et al., 2015; Robinson et al., 2020; Robinson et al., 2024). Logging effects on soil bacterial communities are less clear; two studies in Malaysian Borneo found no impact of logging on either local bacterial alpha (Shannon) diversity or beta diversity (i.e., spatial heterogeneity in bacterial communities, evaluated using dispersion of community dissimilarities) (Lee-Cruz et al., 2013; Tripathi et al., 2016), while other studies revealed shifts in bacterial community composition between intact forest and selective logging gaps (Robinson et al., 2024) and effects of logging on change in bacterial community composition down the soil profile (Tin et al., 2018). This also highlights the importance of fine spatial resolution in assessment of soil bacterial communities which may vary at the meter- or even centimeter-scale (O’ Brien et al., 2016). Disturbance impacts on soil microbial communities may have consequences for biogeochemical cycling that is regulated by complex reciprocal feedbacks between plants and soil (Bever et al., 2010; Cortois et al., 2016; van der Heijden et al., 2008; van der Putten et al., 2013; Wardle et al., 2004). For example, a reduction in cycling rates of certain key soil nutrients has been associated with selective logging canopy gaps, with soil heterotrophic respiration negatively related to disturbance intensity (Robinson et al., 2024).

Although the preservation of primary forest is crucial for biodiversity conservation and the maintenance of ecosystem services (Gibson et al., 2011), natural and managed restoration of secondary tropical forest has great potential for the recovery of ecosystem functions including C storage (Chazdon, 2008; Melo et al., 2013; Pan et al., 2011; Philipson et al., 2020; Wright, 2010). Rehabilitation of disturbed tropical forest via planting programmes is now a widely used strategy to recover vegetation structure and diversity (Bonner et al., 2019; Celis and Jose, 2011; Shoo et al., 2016). Enrichment planting, the reintroduction of native tree species removed through human disturbance, has been adopted across Southeast Asia as a strategy to restore forest floristic composition toward that of old-growth forest (Axelsson et al., 2024; Banin et al., 2023; Perumal et al., 2016). This practice is often accompanied by liana and climber removal, liberation cutting and clearing of understorey vegetation along planting lines, to facilitate seedling establishment by reducing competition for resources. This approach is currently employed in large-scale restoration projects, for example over the 25,000 ha Ulu-Segama area in the Malaysian state of Sabah where forest degradation has been most extreme (Axelsson et al., 2024; Bartholomew et al., 2024; CIFOR-ICRAF, 2024; Face the Future, 2025; Gaveau et al., 2014; Hayward et al., 2021; Reynolds et al., 2011). However, success of tropical forest planting programmes is generally limited by lack of context-specific knowledge of the ecology of planted tree species, including plant-soil interactions (e.g., Rodrigues et al., 2009). Most enrichment planting studies in Southeast Asian tropical forest have focused on tree survival, growth and productivity, with less consideration of soil biodiversity and functioning (Banin et al., 2023; Perumal et al., 2016; Veryard et al., 2023). Following replanting of clear-felled forest, some studies observed increases in microbial biomass carbon (MBC) (Deng et al., 2010; Nurulita et al., 2016) and bacterial alpha diversity accompanied by bacterial community shifts that may indicate ecosystem recovery (Deng et al., 2010). The small number of studies undertaken in enrichment planted secondary Malaysian forest highlight increases in microbial biomass toward that of old-growth forest (Daisuke et al., 2013; Daljit Singh et al., 2013; Perumal et al., 2016). However, microbial indicators can be system-specific and depend on type, intensity and duration of disturbance (Banning et al., 2011). Understanding of the capacity to restore soil microbial diversity, physicochemical properties and biogeochemical cycling in hyper-diverse tropical forest is broadly limited (Bonner et al., 2019), reflecting a general lack of understanding of the patterns of microbial responses to ecosystem restoration (Banning et al., 2011; Strickland et al., 2017).

To address these knowledge gaps, we surveyed soil microbial community attributes, soil physicochemical properties and forest structural characteristics across old-growth and selectively logged forest undergoing either passive natural regeneration or active restoration by enrichment planting. We tested the following specific hypotheses:

1. Soil microbial communities (composition and diversity) will differ across forest types, with greater similarity between old-growth and actively restored forest.

2. Microbial biomass will increase with active restoration toward old-growth forest.

3. Soil physicochemical properties, microclimate and forest structural characteristics will differ across forest types, corresponding to expected differences in soil microbial community attributes.

2 Materials and methods

2.1 Study sites

This study was conducted in the state of Sabah, northern Malaysian Borneo. The climate is moist tropical (annual precipitation 2,600–2,700 mm, average daily temperature 27°C) and generally non-seasonal, but may undergo irregular inter-annual dry periods averaging ∼1.4 months of the year (Kumagai and Porporato, 2012; Walsh and Newbery, 1999). Sampling was conducted in March 2018 at nine sites distributed across old-growth and naturally regenerating or actively restored selectively logged lowland dipterocarp rainforest (3 sampling sites in each) (Figure 1). Logged forest sites were situated within the Innoprise Face Foundation Rainforest Rehabilitation Project (INFAPRO) area (4.99°N, 117.86°E). This large-scale restoration initiative, in partnership with the Yayasan Sabah Foundation, aims to restore 25,000 ha within the Ulu-Segama forest management unit (Face the Future, 2025). All logged forest sites were situated within the same logging coupe that was selectively logged once in 1989. The actively restored forest has undergone rehabilitation by enrichment planting since 2000, with mixtures of dipterocarp species, non-dipterocarp canopy-forming species and various fruit tree species planted at 3 m intervals along 10 m spaced parallel transects. Planting lines were maintained by regular liana cutting and removal of competing understorey vegetation for the initial 3 years following enrichment planting (Bartholomew et al., 2024; CIFOR-ICRAF, 2024; Hayward et al., 2021; Moura Costa, 1996). Old-growth forest was located in the adjacent Danum Valley Conservation Area (DVCA) (4.95°N, 117.79°E), a 438 km2 rainforest reserve that has undergone little or no anthropogenic disturbance having been legally protected from commercial timer operations since 1976 (Marsh and Greer, 1992). Recent extensive tree community surveys across the Ulu-Segama Forest Reserve (including the coupe sampled in this study) and adjacent DVCA have shown persistent shifts in tree species composition and reduced basal area 23–35 years after logging disturbance, regardless of restoration method (Hayward et al., 2021), with corresponding negative impacts on recruitment and diversity of tree seedlings (Bartholomew et al., 2024).

FIGURE 1
Four-panel figure showing the location and details of study sites in Sabah, Malaysia. Panel (a) highlights Southeast Asia with a focus on Borneo. Panel (b) shows a detailed map of Sabah with the Danum Valley Conservation Area marked. Panel (c) illustrates the locations of nine sampling sites using different shapes for old-growth forest (Danum Valley Conservation Area) and actively restored or naturally regenerating logged forest within the INFAPRO Project Area and includes the Segama River and access roads. Panel (d) illustrates the transect sampling design employed at each sampling site, with example distances between sampling points used for analysis of spatial turnover in soil microbial taxa.

Figure 1. Map of sampling sites in northern Malaysian Borneo (a) in the state of Sabah (b), located within old-growth forest (Danum Valley Conservation Area), and adjacent naturally regenerating and actively restored selectively-logged forest (INFAPRO forest rehabilitation project) (c). At each site, soil was collected following a geospatial sampling design (d), at increasing distances along three transects radiating out from one central point. Dotted lines illustrate example distances between sampling points of different transects (here between transects 1 and 3) used to calculate distance decay indices for evaluation of spatial turnover of soil microbial taxa, in addition to distances along each respective transect.

2.2 Sampling design for soil, microclimate and forest structural characteristics

Nine sampling sites were distributed across a contiguous area of forest that had undergone contrasting logging disturbance and restoration. Three sampling sites were chosen in each of three forest types of old-growth, naturally regenerating and actively restored logged forest, situated between a minimum of 500 m and maximum of 5.7 km apart. At each site, a geospatial transect design was used for soil sampling and measurement of microclimate and forest structural characteristics. This comprised three transects radiating out from one center point positioned at 120° to one another, with the first transect oriented North (Figure 1d). Soil and microclimate sampling points were located at the center point, then at six points along each of the three transects at increasing distances of 10 cm, 30 cm, 90 cm, 2.7 m, 8.4 m and 24.3 m relative to the center. This resulted in nineteen soil samples per site, totaling 171 across all sites, which were analysed individually. Soil cores were collected at each sampling point a using a 3 cm diameter gouge auger to a depth of approximately 10 cm for analysis of soil microbial community attributes and physicochemical analysis. The organic soil layer was collected by separating from underlying mineral soil and sealing in a Ziploc bag, and was then transported to a laboratory. Multiple cores were taken at each sampling point to ensure enough material for biological and physicochemical analyses due to the shallow organic layer depth. Samples were hand-homogenized and ∼10 g subsamples taken for analysis of soil microbial community attributes. These were frozen at −20°C on the day of collection and transported on ice to the United Kingdom for analysis of soil microbial communities: 5 g was transported to the United Kingdom Centre for Ecology and Hydrology, Wallingford for amplicon sequencing and 5 g to Lancaster University for phospholipid fatty acid (PLFA) analysis. The remaining soil was transported to the Sabah Forest Research Centre, Sepilok for physicochemical analysis.

Microclimate variables were measured on the day of soil sample collection. Soil temperature (approximate depth 0–10 cm) and understory air temperature at the soil surface (5–15 cm) were measured with using a thermistor (Salter, United Kingdom). Understory photosynthetically active radiation (PAR) was measured using a light meter (PP Systems, United States) with the sensor held just above the soil surface.

Forest structure was evaluated by recording and measuring circumferences of all stems with diameter at breast height (DBH) >5 cm within a 2.5 m buffer of all transects, for calculation of stem density and basal area. This DBH was chosen to capture finer-scale variation in tree abundance. Stem density and basal area were calculated for each sampling location (n = 9).

2.3 Soil physicochemical analysis

pH in water was measured on fresh soils using a pH meter with a combination glass-calomel electrode (a ratio of 1:2.5 soil to deionized water) after shaking overnight at 100 rev min–1 on an orbital shaker and standing for 30 min (Landon, 1984). The remaining soils were air-dried at 40°C to constant weight and passed through a 2 mm sieve for homogenization and removal of roots and stones. Subsamples for total C and N analysis were dried at 65°C for 48 h and milled to a fine powder with a pestle and mortar. Total soil C and N contents were determined by dry combustion at 900°C using an Elementar Vario Max CN analyser (Elementar Analysensysteme, Hanau, Germany). For soil total P, samples were digested using sulphuric acid-hydrogen peroxide (Allen, 1989). Inorganic P was extracted using a Bray No. 1 extractant (Bray and Kurtz, 1945). P contents of extracts and digests were determined using the molybdenum-blue method (Anderson and Ingram, 1993), read at 880 nm on a spectrophotometer (HITACHI-UV-VIS, Japan).

2.4 Soil phospholipid fatty acid (PLFA) analysis

A subset of soil samples were analysed for PLFAs to provide indicators of total microbial biomass and relative abundances of bacteria and fungi across forest types. Due to the intensive extraction requirements, 36 of the 171 samples were analysed (four samples per sampling location). These corresponded in each site to the transect center point and 30 cm, 2.7 m and 24.3 m sampling points along the first transect arm. PLFAs were extracted from 1.8 g freeze dried soil after removal of coarse roots and stones using a modified Bligh and Dyer extraction method (White et al., 1979). Extracts were analysed using an Agilent 6890 Gas Chromatograph with Flame Ionization Detector (GC-FID; Agilent Technologies, Unites States) using an RTx-1 capillary column (60 m × 0.32 mm ID, 0.25 μm film thickness). PLFA peaks were identified using retention times calibrated against known standards. As indicators of Gram-positive bacterial biomass, the branched-chained fatty acids C15:0i, C15:0a, C16:0i, 7Me-C17:0, C17:0i and C17:0a were used (Haack et al., 1994; Lechevalier and Lechevalier, 1988; O’Leary and Wilkinson, 1988; Rinnan and Bååth, 2009; Whitaker et al., 2014; Zelles, 1999). For Gram-negative bacteria, the monounsaturated fatty acids C16:1ω7c, C16:1ω5, C18:1ω7c and cyclopropane fatty acids cyC17:0 and cyC19:0 were used (Rinnan and Bååth, 2009; Whitaker et al., 2014; Zelles, 1999). For fungi, the fatty acids C18:2ω6,9 and C18:1ω9 were used (Bååth and Anderson, 2003; de Deyn et al., 2011). Total bacterial biomass was calculated as the sum of Gram-positive and Gram-negative PLFAs and the fatty acid C15:0 (de Deyn et al., 2011). Fungal to bacterial ratio (F:B) was calculated as the proportion of total bacterial relative to total fungal PLFAs. Total microbial PLFAs were determined as the sum of all identified PLFAs, including those above and the additional fatty acids C14:0, C16:1, C16:0, C17:1ω8, C17:0br, C18:0br, C18:1ω5, C18:0 and C19:1. PLFA contents were expressed as μg g–1 dry soil. Four samples (two samples from old-growth and two from naturally regenerating logged forest, all from different sites) were omitted from analysis due to extraction error, resulting in a total of 32 PLFA samples.

2.5 Molecular analysis of soil microbial communities and data pre-processing

Molecular analyses and bioinformatics were conducted using all soil samples collected (n = 171; 57 per forest type) following methods described by Robinson et al. (2024; 2020). DNA was extracted from 0.2 g soil using the PowerSoil® DNA Isolation Kit and protocol (MoBio Laboratories). Amplicon libraries were constructed according to a dual indexing strategy with each primer consisting of the appropriate Illumina adapter, 8-nt index sequence, a 10-nt pad sequence, a 2-nt linker and the amplicon specific primer (Kozich et al., 2013). Bacteria were targeted using V3-V4 16S rRNA amplicon primers CCTACGGGAGGCAGCAG and GCTATTGGAGCTGGAATTAC (Kozich et al., 2019). For fungi, the ITS2 region was amplified using primers fITS7 GTGARTCATCGAATCTTTG (Ihrmark et al., 2012) and ITS4 TCCTCCGCTTATTGATATGC (White et al., 1990). Although the capability of detecting AM fungi using ITS primers is debated (Hart et al., 2015), recent studies have shown that patterns in diversity and community composition can be adequately identified within sample types such as soil (Berruti et al., 2017; Lekberg et al., 2018). Amplicons were generated using a high-fidelity DNA polymerase (Q5 Taq, New England Biolabs). After an initial denaturation at 95°C for 2 min, PCR conditions were as follows: Denaturation at 95°C for 15 s; annealing at 55°C (bacteria) 52°C (fungi); annealing times were 30 s with extension at 72°C for 30 s; cycle numbers were 25 for bacteria and fungi; a final extension of 10 min at 72°C was included. Amplicon sizes were determined using an Agilent 2200 TapeStation system, samples were normalized using SequalPrep Normalization Plate Kit (Thermo Fisher Scientific) and pooled. The pooled library was quantified using a Qubit dsDNA HS kit (Thermo Fisher Scientific) prior to sequencing with an Illumina MiSeq using V3 600 cycle reagents at a concentration of 8 pM with a 5% PhiX Illumina control library. Sequences were processed in R using DADA2 to quality filter, merge, de-noise and assign taxonomies (Callahan et al., 2016). Forward sequence reads were used for 16S (trimmed to 250 bases), while forward and reverse were used for ITS (trimmed to 225 and 160 bases, respectively). Filtering settings were maximum number of Ns (maxN) = 0, maximum number of expected errors (maxEE) = 1. Sequences were dereplicated and the DADA2 core sequence variant inference algorithm applied. mergePairs and remove BimeraDenovo functions were used at default settings to merge ITS forward and reverse reads and remove chimeric sequences. The amplicon sequence variants (ASVs) were subject to taxonomic assignment using assignTaxonomy with default bootstrapping (50) and the training database UNITE version 7.2 (Abarenkov et al., 2010).

Fungal functional guild classifications were assigned to ASVs using the FUNGuild annotation tool (Nguyen et al., 2016). Only ASVs with unambiguous (non-multiple) classifications of “probable” or “highly-probable” confidence rankings were considered for analysis. These were used for calculating relative abundances of fungal guilds and sub-setting saprotrophic, mycorrhizal, ectomycorrhizal and pathogenic fungal datasets for assessment of diversity and community dissimilarity.

Sequencing data were pre-processed and alpha diversity indices (ASV richness, Shannon index) and fungal guild relative abundances calculated in R (R Core Team, 2025) using the phyloseq package (McMurdie and Holmes, 2013). PCRs failed for three bacterial samples (two from different naturally regenerating logged forest sites, one from actively restored logged forest) as indicated by abnormally low read counts, and these samples were excluded from all subsequent analyses. Only ASVs assigned to the kingdoms of Bacteria or Fungi were retained for downstream analysis, and all singleton ASVs were removed. Sub-setting by fungal guilds was conducted on the full unrarefied dataset to maximize the number of ASV reads available for analysis of functional groups. Sample sequencing depth was normalized for each group by rarefying to the minimum read counts per sample for bacterial (3,778 reads), and total fungal (3,868), saprotrophic (472), mycorrhizal (20), ectomycorrhizal (6) and pathogenic (69) fungal groups.

2.6 Statistical analyses

All statistical analyses were conducted in R (R Core Team, 2025), and significance of all tests was considered at the p ≤ 0.05 level. To test the differences in univariate soil microbial community attributes (alpha diversity metrics and fungal guild relative abundances) and soil, microclimate and forest structural characteristics between forest types, linear mixed effects regression models (LMMs) were constructed in the lme4 R package (Bates et al., 2015) with site ID included as a random intercept term to control for within-site pseudoreplication. Pairwise comparisons between forest types were conducted with the emmeans R package with Bonferroni correction (Lenth et al., 2019). Normality of model residuals were evaluated using Shapiro-Wilk tests and Q-Q plots, and variables were log-, square root- or exp- transformed where necessary to improve model fit.

Soil microbial community data were Hellinger-transformed prior to analysis (Legendre and Borcard, 2018) to control for the effect of rare taxa and merged at the site level (n = 9) to control for spatial pseudoreplication. Soil microbial community compositions across forest types were visualized with PCoA using Bray-Curtis dissimilarities. Differences in soil microbial community composition between forest types were tested with PERMANOVA in vegan (Oksanen et al., 2019) and homogeneity of multivariate dispersion between forest types was evaluated. All permutational tests were run with 9,999 permutations. Pairwise comparisons of soil microbial community dissimilarities between forest types could not be carried out due to the low number of true replicates, restricting the number of possible permutations for calculating significance level. UPGMA (unweighted pair-group method with mathematic average) hierarchical cluster analysis was performed using the hclust R function (R Core Team, 2025) to identify groups of more (dis)similar sites across forest types.

We employed a distance-decay approach to assess the response of microbial beta diversity to forest management, a biogeographical method widely used to study spatial turnover in macro-organism communities that is now being applied to soil microbes due to advances in molecular techniques (Green et al., 2004; Martiny et al., 2006; Morlon et al., 2008; Nekola and White, 1999). Distance-decay relationships can be derived using geographic distances from centimeters to kilometers simultaneously (e.g., Barreto et al., 2014). This allowed us to evaluate microbial community turnover with environmental heterogeneity across scales not previously studied in old-growth and disturbed tropical forest, and to address unresolved disturbance effects on soil microbial diversity that may be spatially dependent. Pairwise Bray-Curtis community dissimilarities for soil microbial groups and corresponding geographic distances were calculated between the 19 sampling points in each sampling location. This provided a total of 171 pairs per location, with geographic distances ranging from 10 cm to 42.09 m. Linear regression was used to obtain the coefficient of the distance decay relationship (Y), i.e., rate of spatial turnover in soil microbial taxa, between log-transformed Bray-Curtis dissimilarities and log-transformed geographic distances for each location and soil microbial group (Barreto et al., 2014; Nekola and White, 1999). ANOVA with Tukey HSD post hoc tests were used to test differences in Y-values between forest types (n = 9).

Numbers of shared and distinct soil microbial ASVs between forest types were visualized with Venn diagrams using the RAM R package (Chen et al., 2018). Indicator analysis was conducted to identify specific soil microbial taxa uniquely associated with different forest types using the mulitpatt function in the indicspecies R package (de Cáceres and Legendre, 2009). This provided an index of association (Indicator Value) between forest type and soil microbial ASVs and p-values denoting significant indicator taxa.

3 Results

3.1 Soil microbial diversity, community composition and biomass across old-growth, naturally regenerating and actively restored logged forest

In total, 37,379 bacterial ASVs (representing 44 identified phyla; 749 genera) and 21,298 fungal ASVs (12 phyla; 611 genera) were detected across all forest types (Figure 2 and Supplementary Figures 1, 2). In all soil microbial groups, a greater number of ASVs were only found in old-growth forest (14,000 bacterial, 7,695 fungal ASVs) compared to naturally regenerating (9,479 bacterial, 5,247 fungal) and actively restored logged forest (8,321 bacterial, 5,323 fungal), and generally more ASVs were shared between old-growth and naturally regenerating forest than between old-growth and actively restored forest (Supplementary Figure 3). Differences in alpha diversity indices or rate of spatial turnover of taxa (or both) were detected between forest types for all microbial groups studied (Figure 3). Mean bacterial Shannon alpha diversity was significantly higher in naturally regenerating forest relative to old-growth, while actively restored forest was similar to both (Figure 3g and Tables 1, 2). Total fungal richness and Shannon alpha diversity were significantly higher in old-growth relative to naturally regenerating and actively restored forest, while total fungal Shannon alpha diversity was higher in actively restored relative to naturally regenerating forest (Figures 3b, h). Saprotrophic fungal richness and Shannon alpha diversity were significantly higher in old-growth compared to both logged forest types (Figures 3c, i). Pathogenic fungal richness was higher in old-growth compared to naturally regenerating forest, while actively restored forest was similar to both (Figure 3f). Rates of spatial turnover of bacterial, total mycorrhizal fungal and ectomycorrhizal fungal taxa (distance decay index) differed by forest type. For bacteria, significantly lower Y-values (indicating slower spatial turnover) were found in naturally regenerating logged forest compared to old-growth and actively restored forest, which were similar (Figure 3m). Slower spatial turnover of both total mycorrhizal and ectomycorrhizal fungi was found in actively restored forest relative to old-growth forest. In naturally regenerating forest, spatial turnover of total mycorrhizal fungal taxa was similar to both other forest types, while spatial turnover of ectomycorrhizal fungal taxa was slower than old-growth and similar to actively restored forest (Figures 3p, q).

FIGURE 2
Two bar charts showing relative abundance (as percentages) of different phyla across three forest types of old-growth forest and actively restored or naturally regenerating logged forest. Chart (a) represents bacteria with phyla such as Proteobacteria and Acidobacteria. Chart (b) represents fungi with phyla including Ascomycota and Basidiomycota. Each color corresponds to a different phylum, detailed in the legends.

Figure 2. Relative abundances of bacterial (a) and fungal (b) phyla as percentages of total amplicon sequence variants (ASVs) across old-growth forest (OG), naturally regenerating logged forest (NR) and actively restored logged forest (AR) (n = 168 and n = 171 for bacterial and fungal phyla, respectively). Phyla with <1% relative abundance (rel. ab.) across all forest types are represented as one group.

FIGURE 3
Nine-panel figure showing variation of microbial diversity measures across forest types (old-growth forest and naturally regenerating or actively restored logged forest). Panels (a-f) represent richness in amplicon sequence variants for bacterial and fungal types. Panels (g-l) show Shannon alpha diversity indices for the same groups. Panels (m-r) depict rate of spatial turnover of microbial taxa for the same groups. Data points are marked by color: black (old-growth forest), orange (naturally regenerating logged forest), and blue (actively restored logged forest), with error bars and labelled significance levels.

Figure 3. Indices of soil microbial alpha diversity [richness as no. observed ASVs per 10 reads (a–f), and Shannon index (g–l)) and rate of spatial turnover of taxa (distance decay index; Y (m–r)] derived from amplicon sequence variants (ASVs) across old-growth forest (OG; black), naturally regenerating logged forest (NR; orange) and actively restored logged forest (AR; blue). Values are means ± 1 standard error. Lower case letters indicate statistically different or similar groups at the p < 0.05 level identified in post hoc tests after linear mixed model analysis. Sample numbers analysed were n = 168 for bacterial richness and Shannon diversity, n = 171 for fungal richness, Shannon diversity and fungal guild relative abundances, and n = 9 for both bacterial and fungal distance decay indices. See Table 2 for a summary of statistical test results.

TABLE 1
www.frontiersin.org

Table 1. Soil microbial community attributes in old-growth forest, naturally regenerating logged forest and actively restored logged forest (means ± 1 standard deviation). Superscript letters indicate statistically different or similar groups at the p < 0.05 level identified in post hoc tests after linear mixed model or Kruskal-Wallis analysis. Sample numbers analysed were n = 32 for Total PLFAs and Fungal: bacterial ratios, n = 168 for bacterial richness and Shannon diversity, n = 171 for fungal richness, Shannon diversity and fungal guild relative abundances, and n = 9 for both bacterial and fungal distance decay indices.

TABLE 2
www.frontiersin.org

Table 2. Linear model or Kruskal-Wallis* test statistics for significant differences in soil microbial alpha diversity metrics (richness and Shannon index), spatial turnover of taxa (distance decay index) and fungal guild relative abundances across old-growth forest (OG), naturally regenerating logged forest (NR) and actively restored logged forest (AR). Summaries are given for overall models and post hoc comparisons between forest types. p-values for pairwise tests for LMM and Kruskal-Wallis analyses were adjusted using the Tukey and Bonferroni methods, respectively. Significant p-values (> 0.05) are highlighted in bold. Sample numbers analysed were n = 168 for bacterial richness and Shannon diversity, n = 171 for fungal richness, Shannon diversity and fungal guild relative abundances, and n = 9 for both bacterial and fungal distance decay indices.

Bray-Curtis community dissimilarities were significantly affected by forest type for all soil microbial groups, with the exception of ectomycorrhizal fungi, which was marginally non-significant (Figure 4 and Table 3). Community dissimilarity dispersions were homogenous between all land-use types for all fungal groups (p > 0.05). UPGMA hierarchical clustering analysis (Supplementary Figure 4) identified communities in old-growth sites to be most dissimilar to logged forest sites across all microbial groups, with the exception of ectomycorrhizal fungi (bacteria: 86.52% dissimilarity; overall fungi: 89.84%; saprotrophic fungi: 87.46%; mycorrhizal fungi: 97.27%; pathogenic fungi: 78.91%). For ectomycorrhizal fungi, the greatest dissimilarity was found between OG-1 and all other sites (98.44%), followed by dissimilarity between OG-2 and OG-3 sites and all others (97.59% dissimilarity). Logged forest sites were generally more similar, with naturally regenerating and actively restored sites clustering together for all soil microbial groups.

FIGURE 4
Six plots show Principal Coordinates Analysis (PCoA) results for different microbial groups across old-growth forest, naturally regenerating logged forest and actively restored logged forests. Panels (a-f) represent bacterial, total fungal, saprotrophic fungal, mycorrhizal fungal, ectomycorrhizal fungal, and pathogenic fungal communities. The plots display axes labelled with percentages indicating the variance explained, and p-values indicating statistical significance of the effect of forest type on community dissimilarity.

Figure 4. Principal coordinates analysis (PCoA) ordinations of Bray-Curtis dissimilarities for soil bacterial (a) and fungal (b–f) communities studied across old-growth forest (black), naturally regenerating logged forest (orange) and actively restored logged forest (blue) using data merged at the site level (n = 9). p-values denote significance of differences in community composition across forest types in PERMANOVA tests. See Table 3 for a full summary of statistical test results.

TABLE 3
www.frontiersin.org

Table 3. PERMANOVA test statistics for the effect of forest type on soil microbial community Bray-Curtis dissimilarities using data merged at the site level (n = 9). Significant p-values (< 0.05) are highlighted in bold.

Fungal guild relative abundances for total pathogens, plant pathogens, animal pathogens and lichens significantly differed by forest type (Figure 5 and Tables 1, 2). Total pathogenic and lichenized fungal relative abundances were significantly higher in old-growth relative to actively restored forest (naturally regenerating forest was similar to both old-growth and actively restored forest types for both fungal guilds), while plant pathogenic fungal relative abundance was significantly higher in old-growth relative to both logged forest types. Animal pathogenic fungal relative abundance was higher in naturally regenerating forest relative to actively restored forest (old-growth was similar to both logged forest types). PLFA analysis indicated no differences in total microbial, fungal or bacterial biomass or F:B between forest types (p > 0.05 in overall and pairwise tests).

FIGURE 5
Bar chart comparing the relative abundance of fungal ASV reads across different forest types: old-growth forest (black), naturally regenerating logged forest (orange), and actively restored logged forest (blue). Categories include saprotrophic, mycorrhizal and pathogenic fungi. Statistical significance of the effect of forest type is marked with asterisks, and statistically different or similar forest types for fungal guilds are denoted by lower case letters. Error bars indicate variability within forest types (one standard error).

Figure 5. Relative abundances of fungal guilds (means ± 1 standard error) across old-growth forest (black), naturally regenerating logged forest (orange) and actively restored logged forest (blue) (n = 171). Horizontal bars and asterisks indicate significant overall effects of forest type on relative abundances of fungal guilds identified by linear mixed model or Kruskal-Wallis analysis, **p < 0.01, *p < 0.05. Lower case letters indicate statistically different or similar groups by forest type within fungal guilds at the p < 0.05 level identified in post hoc tests. See Table 2 for a summary of statistical test results. Means ± 1 standard error for fungal guilds with low relative abundances in old-growth forest (OG), naturally regenerating logged forest (NR) and actively restored logged forest (AR) are as follows: ericoid mycorrhizal fungi – OG: 0.09 ± 0.05, NR: 0.10 ± 0.03, AR: 0.14 ± 0.04; endophytic fungi - OG: 0.22 ± 0.07, NR: 0.47 ± 0.15, AR: 0.14 ± 0.08; epiphytic fungi - OG: 0.14 ± 0.11, NR: 0.05 ± 0.03, AR: 0.04 ± 0.02; lichenized fungi - OG: 1.45 ± 0.16, NR: 0.17 ± 0.04, AR: 0.05 ± 0.01.

3.2 Microbial indicators of forest types

Indicator analysis of bacterial taxa identified 566 significant indicator ASVs for old-growth forest, 99 for naturally regenerating forest and 139 for actively restored forest (Supplementary Data Sheet 2, see also for indicator ASV taxonomic classifications). At the phylum level, a larger proportion of indicator ASVs unique to old-growth appeared to belong to the Firmicutes in comparison to both logged forest types (Supplementary Figure 5a). Only two significant bacterial indicator ASVs were shared by old-growth and actively restored logged forest, both belonging to the phylum Firmicutes, class Bacilli (Alicyclobacillus sp. and Bacillus foraminis). 205 significant fungal indicator ASVs were identified for old-growth forest, 45 for naturally regenerating forest and 94 for actively restored forest. Notable differences in composition of fungal indicator taxa phyla was Mortierellomycota appearing to occur mostly only in old-growth forest, and Mucoromycota mostly only in actively restored forest (Supplementary Figure 5b). Four significant indicator ASVs were shared between old-growth and actively restored forest from the two phyla Ascomycota and Basidiomycota, and two classes Sordariomycetes and Tremellomycetes (Trichoderma deliquescens, Clonostachys rosea, Castanediella sp. and Saitozyma podzolica).

3.3 Soil physicochemical properties, microclimate and forest structural characteristics across old-growth, naturally regenerating and actively restored logged forest

Of all the soil physicochemical properties, microclimate and forest structural characteristics measured, only soil pH significantly differed between forest types (Table 4; overall model: R2 = 0.83, F = 50.56, p < 0.001). Post hoc tests identified significantly higher soil pH in old-growth forest relative to both logged forest types (p < 0.001), which did not significantly differ between themselves (p > 0.05). Mean understorey PAR was highest in naturally regenerating logged forest, although not significantly different due to high heterogeneity in this forest type.

TABLE 4
www.frontiersin.org

Table 4. Soil physicochemical properties, microclimate and vegetation characteristics in old-growth forest, naturally regenerating logged forest and actively restored logged forest (means ± 1 standard deviation). Superscript letters indicate statistically different or similar groups at the p < 0.05 level identified in post hoc tests after linear mixed model analysis. Sample numbers analysed were n = 171 for soil and microclimate variables and n = 9 for forest structural characteristics (i.e., one value per site).

4 Discussion

We evaluated responses in soil microbial community attributes, soil physicochemical properties, microclimate and forest structural characteristics to active and passive restoration of logged Bornean lowland dipterocarp rainforest relative to old-growth forest across different spatial scales. Bacterial community attributes differed by forest type, including community composition (Figure 4a), alpha diversity metrics and rate of spatial turnover of taxa (Figures 3a, g, m), broadly supporting our first hypothesis. Specifically, results showed higher bacterial alpha diversity and lower rates of spatial turnover of taxa in naturally regenerating forest, in agreement with the observation that forest disturbance increases local bacterial diversity while homogenizing communities over larger spatial scales (Rodrigues et al., 2013). Although tropical forest disturbance is often accompanied by declines in both alpha and beta diversity of aboveground organisms (Bierregaard et al., 2001; Sodhi et al., 2009), this is not necessarily the case for belowground microbial communities where the opposite trend may be seen in local diversity patterns (Petersen et al., 2019). As discussed by Rodrigues et al. (2013), this “decoupling” of alpha and beta diversity with disturbance may depend on the relative effects of disturbance on ecosystem productivity. For example, alpha diversity in aboveground communities has been shown to increase with disturbance when productivity rates are higher in resulting ecosystems (Smart et al., 2006). In microbial terms, old-growth rainforests may be characterized by relatively low belowground productivity compared to adjacent open ecosystems created through anthropogenic disturbance (Cenciani et al., 2009; Cerri et al., 2004). In the current study system, increased bacterial alpha diversity with SL may potentially result from changes in vegetation characteristics, e.g., changes in tree community composition, or increasing understorey vegetation through creation of canopy gaps (Denslow, 1995) affecting the quality and quantity of plant inputs to the soil, improving effective resource availability (Cenciani et al., 2009; Cerri et al., 2004) and creating different ecological niches for bacterial communities (Zhang et al., 2018). This observation is congruent with the intermediate disturbance hypothesis, which predicts increases in diversity with a certain level of disturbance (Ferrenberg et al., 2013; Galand et al., 2016; Wilkinson, 1999; Zhang et al., 2011). While the measured soil properties were largely similar between forest types, soil pH, an important determinant of bacterial communities in Southeast Asian tropical forests (Tripathi et al., 2012), was significantly greater in old-growth forest sites (Table 4). The observed differences in soil bacterial community attributes disagree with some previous work in the same region of Borneo that found bacterial communities and diversity to be broadly resilient to selective logging (Lee-Cruz et al., 2013; Tripathi et al., 2016). This is likely due to the coarse spatial sampling resolution used in these studies to evaluate alpha and beta diversity (composite samples comprising soil collected up to 200 m apart). As bacterial community structure can vary considerably over meter-and centimeter-scales (O’ Brien et al., 2016), sampling resolution of previous surveys may be inappropriate for evaluating bacterial diversity and biogeographical patterns in response to forest ecosystem disturbance. The findings of the present study highlight the need for landscape studies of soil microbial diversity to incorporate fine spatial scale approaches to identify impacts and implications for biogeochemical cycling. Although clear differences in bacterial alpha diversity and spatial turnover of taxa were detected in naturally regenerating forest, findings suggest that restoration of selectively forest by enrichment planting can recover these metrics to levels comparable to old-growth forest. Bacterial alpha diversity may have potential as an indicator of rehabilitation of selectively logged forest, with lower values representing ecosystem recovery. This is opposite to trends in bacterial alpha diversity used to monitor progress of forest rehabilitation after total clearance, which can increase with replanting (Nurulita et al., 2016). This emphasizes the importance of disturbance history in identifying appropriate context-specific recovery indicators. Indicator analysis revealed old-growth forest bacterial communities harbor a large number of unique taxa that are not found in nearby logged forest, also evident for soil fungi (Supplementary Figure 5a, b), suggesting some aspects of ecosystem complexity are lost through disturbance and not recovered after almost two decades of active restoration or natural regeneration. A large proportion of unique bacterial taxa in old-growth forest belonged to the Firmicutes, a phylum associated with high soil C availability and resilience to environmental (microclimatic) perturbations (Battistuzzi and Hedges, 2009; Rodrigues et al., 2013), that were largely absent in unique taxa of both logged forests. While we found no differences in total soil carbon between forest types, greater abundance of unique taxa in this phylum may result from changes in the quality of plant-derived carbon inputs to the soil, e.g., logging disturbance has been linked to increased recalcitrance of tree litter carbon and reduced functional breadth of microbial decomposer communities relative to old-growth forest in the same area (Elias et al., 2020).

Some fungal community attributes also differed by forest type. While results broadly agree with observed logging disturbance effects of fungal communities (Kerfahi et al., 2014; McGuire et al., 2015; Robinson et al., 2020; Robinson et al., 2024), our results did not corroborate our prediction that overall fungal community composition would be more similar between old-growth and actively restored forest. Total, saprotrophic, and mycorrhizal fungal communities significantly differed by forest type, with old-growth forest communities appearing distinct from naturally regenerating and actively restored forest which were generally more similar (Figures 4b–f). Differences observed in the composition and rate of spatial turnover of mycorrhizal communities across forest types is consistent with studies showing high sensitivity of mycorrhizae to logging disturbance, likely due to the targeted extraction of ectomycorrhizal-associating dipterocarp trees (Kerfahi et al., 2014; Robinson et al., 2020; Robinson et al., 2024). The composition of mycorrhizal and ectomycorrhizal fungal communities in actively restored forest appeared even more distant to old-growth forest than naturally regenerating forest (Figures 4d, e), possibly reflecting differences in the mycorrhizal associations of tree species selected for planting during active restoration and those removed during timber extraction. Similarly, the rate of spatial turnover of mycorrhizal taxa was most different between old-growth and actively restored forest (the latter approximately three times lower), with rates in naturally regenerating forest similar to both (Figure 3p). Mean overall mycorrhizal and ectomycorrhizal richness, Shannon alpha diversity and relative abundances were also highest in actively restored forest, but did not significantly differ due to high variability (and likely resulting from low number of true spatial replicates in each forest type). There were surprisingly few fungal taxa shared between old-growth and actively restored forest (420 ASVs; Supplementary Figure 3b) relative to those shared between the two logged forest types (1,334 ASVs), or even old-growth and naturally regenerating forest (659 ASVs) – a pattern also evident in bacteria (384 shared between old-growth and actively restored forest; 2,843 between logged forest types; 1,047 between old-growth and naturally regenerating forest; Supplementary Figure 3a). This suggests that some taxa may be lost due to current restoration practices that are otherwise present in both old-growth and naturally regenerating forest. This is reflected in the lower total fungal alpha diversity in actively restored forest relative to old-growth (Figures 3b, h), which is likely driven by lower saprotrophic fungal alpha diversity (Figures 3c, i) as this guild represented the largest proportion of total fungal reads (47.17%). These findings may be potentially related to long-term control of liana species which have their own soil microbial associations (McGuire et al., 2008; Schnitzer et al., 2005), and removal of understorey vegetation which is practiced in these study sites. Vegetation removal may alter plant litter inputs, which in turn can affecting resulting microbial decomposer communities (Shi et al., 2019). Further study of the effects of liana and understorey vegetation removal through controlled field experiments (i.e., enrichment planting with and without additional vegetation clearance) is required to unpick the underlying drivers of these observations, and possible role in the impediment of fungal community recovery toward characteristics of old-growth forest.

No differences were found in indicators of overall microbial biomass between forest types, refuting our second hypothesis and contrasting with previous studies observing clear reductions in MBC in degraded forest relative to old-growth (Deng et al., 2010; Nurulita et al., 2016), or higher MBC in restored versus unrestored forest (Daljit Singh et al., 2013). In the present study, it is possible either microbial biomass was unaffected by logging in these forests, or returned to comparable levels with old-growth in both logged forest types with natural or managed regeneration during time since disturbance. The small number of samples used for analysis (n = 32) may also have contributed to lack of differences found due to a large amount of within- forest type (and site) variation.

Of all the soil physicochemical properties, microclimate and forest structural characteristics measured, only soil pH was affected by forest type (third hypothesis) with more acidic soils found in both logged forest types relative to old-growth (Table 4). Interestingly, bacterial alpha diversity was found to be higher in unrestored forest soils relative to old-growth despite lower pH, contrasting with studies across multiple biomes (including Malaysia) which found bacterial alpha diversity to increase with soil neutrality over land-use gradients (Lauber et al., 2009; Tripathi et al., 2012). Bacterial alpha diversity may therefore be influenced more strongly by altered litter inputs while soil microbial successional processes are still underway, rather than soil pH which may govern attributes of climax bacterial communities (also see Robinson et al., 2024).

5 Conclusion

In conclusion, while many soil properties and soil microbial attributes were similar across forest types, our results demonstrate contrasting responses of different soil microbial groups to active and passive restoration of selectively logged forest. Bacterial and fungal community composition remained generally more similar between logged forest types and more distinct in old-growth forest following 18 years of natural regeneration or enrichment planting. Bacterial alpha diversity and spatial turnover of bacterial taxa may recover toward old-growth forest levels with active restoration, while fungal alpha diversity showed slower signs of recovery largely due to saprotrophic fungal alpha diversity remaining lower in both logged forest types relative to old-growth forest. The composition and rate of spatial turnover in mycorrhizal communities was most different between old-growth forest and actively restored forest, possibly resulting from discrepancies in respective mycorrhizal associations of tree species planted during restoration and those removed during timber harvest. Few fungal taxa were shared between old-growth and actively restored forest, indicating recovery of soil microbial communities may be impeded by current management practices, with implications for carbon cycling. Further study into the effects of liana and understorey vegetation removal through controlled experimentation is required to test underlying mechanisms. Taken together our findings emphasize the importance of evaluating belowground microbial communities during forest restoration, particularly at fine spatial (cm to m) scales, to evaluate and predict recovery of biodiversity and ecosystem functions in human modified tropical forest.

Data availability statement

The molecular datasets generated for this study are available through the National Center for Biotechnology Information Sequence Read Archive (NCBI-SRA), Project no.: PRJNA1218937 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1218937). The datasets of soil microbial community attributes (including ASV abundance tables, associated taxonomic and functional classifications, diversity indices, relative abundances of fungal guilds), soil properties, forest structural and microclimate parameters are available through the UKCEH Environmental Information Data Centre (EIDC) as part of the NERC Environmental Data Service (Robinson et al., 2025; https://doi.org/10.5285/08bfe302-d33f-490d-97be-27bb83a0f38d).

Author contributions

SR: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review and editing. DE: Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review and editing. TG: Data curation, Methodology, Writing – original draft, Writing – review and editing. NiM: Conceptualization, Supervision, Writing – original draft, Writing – review and editing. RG: Supervision, Writing – original draft, Writing – review and editing, Methodology. NoM: Conceptualization, Supervision, Writing – original draft, Writing – review and editing. YW: Conceptualization, Supervision, Writing – original draft, Writing – review and editing. NO: Conceptualization, Supervision, Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded through the UK Natural Environment Research Council (NERC) Human Modified Tropical Forests Programme. This publication is a contribution from the NERC Biodiversity and Land-use Impacts on Tropical Ecosystem Function (BALI) consortium [NE/K016253/1]. SR received additional funding from the NERC ENVISION: Developing next generation leaders in environmental science Doctoral Training Scheme [NE/L002604/1]. SR, DE and NM were also supported by NERC as part of the NC International Programme [NE/X006247/1] delivering National Capability.

Acknowledgments

We express our gratitude to the Sabah Biodiversity Council (SaBC), Yayasan Sabah (Sabah Foundation), Glen Reynolds and the South East Asia Rainforest Research Partnership (SEARRP) for support and granting permission to access field sites [SaBC Access License no. JKM/MBS.1000-2/2 JLD.6 (16)]. We thank Asni Bin Bulot at INFAPRO for his assistance in the field. We are extremely grateful to Laura Kruitbos and Unding Jami for project support and all research assistants in Sabah, without whom this study would not have been possible. We thank James Guy for assistance in soil PLFA analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1570294/full#supplementary-material

References

Abarenkov, K., Nilsson, R. H., Larsson, K.-H., Alexander, I. J., Eberhardt, U., Erland, S., et al. (2010). The UNITE database for molecular identification of fungi – recent updates and future perspectives. New Phytol. 186, 281–285. doi: 10.1111/j.1469-8137.2009.03160.x

PubMed Abstract | Crossref Full Text | Google Scholar

Allen, S. E. (1989). Chemical analysis of ecological materials, 2nd Edn. Oxford: Blackwell Scientific.

Google Scholar

Anderson, J. M. and Ingram, J. S. I. (1993). Tropical soil biology and fertility: a handbook of methods. 2nd Edn. Wallingford: CAB International.

Google Scholar

Appanah, S., and Turnbull, J. M. (1998). A Review of Dipterocarps: Taxonomy, Ecology and Silviculture. Available online at: http://www.cifor.org/publications/pdf_files/Books/Dipterocarps.pdf (accessed July 23, 2025).

Google Scholar

Asner, G. P., Brodrick, P. G., Philipson, C., Vaughn, N. R., Martin, R. E., Knapp, D. E., et al. (2018). Mapped aboveground carbon stocks to advance forest conservation and recovery in Malaysian Borneo. Biol. Conserv. 217, 289–310. doi: 10.1016/j.biocon.2017.10.020

Crossref Full Text | Google Scholar

Axelsson, E. P., Grady, K. C., Alloysius, D., Falck, J., Lussetti, D., Vairappan, C. S., et al. (2024). Lessons learned from 25 years of operational large-scale restoration: The Sow-A-Seed project, Sabah, Borneo. Ecol. Eng. 206:107282. doi: 10.1016/j.ecoleng.2024.107282

Crossref Full Text | Google Scholar

Bååth, E., and Anderson, T. H. (2003). Comparison of soil fungal/bacterial ratios in a pH gradient using physiological and PLFA-based techniques. Soil Biol. Biochem. 35, 955–963. doi: 10.1016/S0038-0717(03)00154-8

Crossref Full Text | Google Scholar

Baccini, A., Walker, W., Carvalho, L., Farina, M., Sulla-Menashe, D., and Houghton, R. A. (2017). Tropical forests are a net carbon source based on aboveground measurements of gain and loss. Science 358:230. doi: 10.1126/science.aam5962

PubMed Abstract | Crossref Full Text | Google Scholar

Banin, L. F., Raine, E. H., Rowland, L. M., Chazdon, R. L., Smith, S. W., Rahman, N. E. B., et al. (2023). The road to recovery: A synthesis of outcomes from ecosystem restoration in tropical and sub-tropical Asian forests. Philos. Trans. R. Soc. B Biol. Sci. 378:20210090. doi: 10.1098/rstb.2021.0090

PubMed Abstract | Crossref Full Text | Google Scholar

Banning, N. C., Gleeson, D. B., Grigg, A. H., Grant, C. D., Andersen, G. L., Brodie, E. L., et al. (2011). Soil microbial community successional patterns during forest ecosystem restoration. Appl. Environ. Microbiol. 77:6158. doi: 10.1128/AEM.00764-11

PubMed Abstract | Crossref Full Text | Google Scholar

Barreto, D. P., Conrad, R., Klose, M., Claus, P., and Enrich-Prast, A. (2014). Distance-decay and taxa-area relationships for bacteria, archaea and methanogenic archaea in a tropical lake sediment. PLoS One 9:e110128. doi: 10.1371/journal.pone.0110128

PubMed Abstract | Crossref Full Text | Google Scholar

Bartholomew, D. C., Hayward, R., Burslem, D. F. R. P., Bittencourt, P. R. L., Chapman, D., Bin Suis, M. A. F., et al. (2024). Bornean tropical forests recovering from logging at risk of regeneration failure. Glob. Change Biol. 30:e17209. doi: 10.1111/gcb.17209

PubMed Abstract | Crossref Full Text | Google Scholar

Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting linear mixed-effects models using lme4. J. Stat. Softw. 67, 1–48. doi: 10.18637/jss.v067.i01

Crossref Full Text | Google Scholar

Battistuzzi, F. U., and Hedges, S. B. (2009). A major clade of prokaryotes with ancient adaptations to life on land. Mol. Biol. Evol. 26, 335–343. doi: 10.1093/molbev/msn247

PubMed Abstract | Crossref Full Text | Google Scholar

Berga, M., Székely, A. J., and Langenheder, S. (2012). Effects of disturbance intensity and frequency on bacterial community composition and function. PLoS One 7:e36959. doi: 10.1371/journal.pone.0036959

PubMed Abstract | Crossref Full Text | Google Scholar

Berruti, A., Desirò, A., Visentin, S., Zecca, O., and Bonfante, P. (2017). ITS fungal barcoding primers versus 18S AMF-specific primers reveal similar AMF-based diversity patterns in roots and soils of three mountain vineyards. Environ. Microbiol. Rep. 9, 658–667. doi: 10.1111/1758-2229.12574

PubMed Abstract | Crossref Full Text | Google Scholar

Bever, J. D., Dickie, I. A., Facelli, E., Facelli, J. M., Klironomos, J., Moora, M., et al. (2010). Rooting theories of plant community ecology in microbial interactions. Trends Ecol. Evol. 25, 468–478. doi: 10.1016/j.tree.2010.05.004

PubMed Abstract | Crossref Full Text | Google Scholar

Bierregaard, R. O. Jr., Gascon, C., Lovejoy, T. E., and Mesquita, R. (2001). Lessons from Amazonia: The ecology and conservation of a fragmented forest. New Haven, CT: Yale University Press.

Google Scholar

Bonner, M. T. L., Allen, D. E., Brackin, R., Smith, T. E., Lewis, T., Shoo, L. P., et al. (2019). Tropical rainforest restoration plantations are slow to restore the soil biological and organic carbon characteristics of old growth rainforest. Microb. Ecol. 79, 432–442. doi: 10.1007/s00248-019-01414-7

PubMed Abstract | Crossref Full Text | Google Scholar

Both, S., Riutta, T., Paine, C. E. T., Elias, D. M. O., Cruz, R. S., Jain, A., et al. (2019). Logging and soil nutrients independently explain plant trait expression in tropical forests. New Phytol. 221, 1853–1865. doi: 10.1111/nph.15444

PubMed Abstract | Crossref Full Text | Google Scholar

Bray, R. H. and Kurtz, L. T. (1945). Determination of total, organic, and available forms of phosphorus in soils. Soil Sci. 59, 39–46. doi: 10.1097/00010694-194501000-00006

Crossref Full Text | Google Scholar

Callahan, B. J., Mcmurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869

PubMed Abstract | Crossref Full Text | Google Scholar

Celis, G., and Jose, S. (2011). Restoring abandoned pasture land with native tree species in Costa Rica: Effects of exotic grass competition and light. Forest Ecol. Manag. 261, 1598–1604. doi: 10.1016/j.foreco.2010.10.005

Crossref Full Text | Google Scholar

Cenciani, K., Lambais, M. R., Cerri, C. C., Azevedo, L. C. B. D., and Feigl, B. J. (2009). Bacteria diversity and microbial biomass in forest, pasture and fallow soils in the Southwestern Amazon basin. Rev. Bras. Ciênc. Solo 33, 907–916. doi: 10.1590/s0100-06832009000400015

Crossref Full Text | Google Scholar

Cerri, C. E. P., Paustian, K., Bernoux, M., Victoria, R. L., Melillo, J. M., and Cerri, C. C. (2004). Modeling changes in soil organic matter in Amazon forest to pasture conversion with the Century model. Glob. Change Biol. 10, 815–832. doi: 10.1111/j.1365-2486.2004.00759.x

Crossref Full Text | Google Scholar

Chazdon, R. L. (2008). Beyond deforestation: Restoring forests and ecosystem services on degraded lands. Science 320, 1458–1460. doi: 10.1126/science.1155365

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, W., Simpson, J., and Levesque, C. A. (2018). RAM: R for Amplicon-Sequencing-Based Microbial-Ecology. Available online at: https://CRAN.R-project.org/package=RAM (accessed April 27, 2022).

Google Scholar

CIFOR-ICRAF. (2024). INFAPRO Rehabilitation of Logged-over Dipterocarp Forest in Sabah, Malaysia. Available online at: https://www.reddprojectsdatabase.org/558-infapro-rehabilitation-of-logged-over-dipterocarp-forest-in-sabah-malaysia/ (accessed July 23, 2025).

Google Scholar

Cortois, R., Schröder-Georgi, T., Weigelt, A., van der Putten, W. H., and de Deyn, G. B. (2016). Plant–soil feedbacks: Role of plant functional group and plant traits. J. Ecol. 104, 1608–1617. doi: 10.1111/1365-2745.12643

Crossref Full Text | Google Scholar

Daisuke, H., Tanaka, K., Joseph Jawa, K., Ikuo, N., and Katsutoshi, S. (2013). Rehabilitation of degraded tropical rainforest using dipterocarp trees in Sarawak, Malaysia. Int. J. For. Res. 2013:683017. doi: 10.1155/2013/683017

Crossref Full Text | Google Scholar

Daljit Singh, K. S., Arifin, A., Radziah, O., Shamshuddin, J., Hazandy, A. H., Majid, N. M., et al. (2013). Status of soil microbial population, enzymatic activity and biomass of selected natural, secondary and rehabilitated forests. Am. J. Environ. Sci. 9, 301–309. doi: 10.3844/ajessp.2013.301.309

Crossref Full Text | Google Scholar

de Cáceres, M., and Legendre, P. (2009). Associations between species and groups of sites: Indices and statistical inference. Ecology 90, 3566–3574. doi: 10.1890/08-1823.1

PubMed Abstract | Crossref Full Text | Google Scholar

de Deyn, G. B., Shiel, R. S., Ostle, N. J., McNamara, N. P., Oakley, S., Young, I., et al. (2011). Additional carbon sequestration benefits of grassland diversity restoration. J. Appl. Ecol. 48, 600–608. doi: 10.1111/j.1365-2664.2010.01925.x

Crossref Full Text | Google Scholar

Deng, H., Zhang, B., Yin, R., Wang, H.-l, Mitchell, S., Griffiths, B., et al. (2010). Long-term effect of re-vegetation on the microbial community of a severely eroded soil in sub-tropical China. Plant Soil 328, 447–458. doi: 10.1007/s11104-009-0124-9

Crossref Full Text | Google Scholar

Denslow, J. S. (1995). Disturbance and diversity in tropical rain forests: The density effect. Ecol. Appl. 5, 962–968. doi: 10.2307/2269347

Crossref Full Text | Google Scholar

Don, A., Schumacher, J., and Freibauer, A. (2011). Impact of tropical land-use change on soil organic carbon stocks - A meta-analysis. Glob. Change Biol. 17, 1658–1670. doi: 10.1111/j.1365-2486.2010.02336.x

Crossref Full Text | Google Scholar

Elias, D. M. O., Robinson, S., Both, S., Goodall, T., Majalap-Lee, N., Ostle, N. J., et al. (2020). Soil microbial community and litter quality controls on decomposition across a tropical forest disturbance gradient. Front. For. Glob. Change 3:81. doi: 10.3389/ffgc.2020.00081

Crossref Full Text | Google Scholar

Ellis, P., Griscom, B., Walker, W., Gonçalves, F., and Cormier, T. (2016). Mapping selective logging impacts in Borneo with GPS and airborne lidar. For. Ecol. Manage. 365, 184–196. doi: 10.1016/j.foreco.2016.01.020

Crossref Full Text | Google Scholar

Face the Future. (2025). INFAPRO Restoration of Tropical Rainforest: Sabah, Malaysia. Available online at: https://facethefuture.com/projects/sabah-maleisie-bosherstel-en-bescherming (accessed July 23, 2025).

Google Scholar

Ferrenberg, S., O’Neill, S., Knelman, J. E., Todd, B., Duggan, S., Bradley, D., et al. (2013). Changes in assembly processes in soil bacterial communities following a wildfire disturbance. ISME J. 7, 1102–1111. doi: 10.1038/ismej.2013.11

PubMed Abstract | Crossref Full Text | Google Scholar

Galand, P. E., Lucas, S., Fagervold, S. K., Peru, E., Pruski, A. M., Vétion, G., et al. (2016). Disturbance increases microbial community diversity and production in marine sediments. Front. Microbiol. 7:1950. doi: 10.3389/fmicb.2016.01950

PubMed Abstract | Crossref Full Text | Google Scholar

Gaveau, D. L. A., Sheil, D., Husnayaen, Salim, M. A., Arjasakusuma, S., Ancrenaz, M., et al. (2016). Rapid conversions and avoided deforestation: Examining four decades of industrial plantation expansion in Borneo. Sci. Rep. 6:32017. doi: 10.1038/srep32017

PubMed Abstract | Crossref Full Text | Google Scholar

Gaveau, D. L. A., Sloan, S., Molidena, E., Yaen, H., Sheil, D., Abram, N. K., et al. (2014). Four decades of forest persistence, clearance and logging on Borneo. PLoS One 9:e101654. doi: 10.1371/journal.pone.0101654

PubMed Abstract | Crossref Full Text | Google Scholar

Gibson, L., Lee, T. M., Koh, L. P., Brook, B. W., Gardner, T. A., Barlow, J., et al. (2011). Primary forests are irreplaceable for sustaining tropical biodiversity. Nature 478:378. doi: 10.1038/nature10425

PubMed Abstract | Crossref Full Text | Google Scholar

Green, J. L., Holmes, A. J., Westoby, M., Oliver, I., Briscoe, D., Dangerfield, M., et al. (2004). Spatial scaling of microbial eukaryote diversity. Nature 432:747. doi: 10.1038/nature03034

PubMed Abstract | Crossref Full Text | Google Scholar

Haack, S. K., Garchow, H., Odelson, D. A., Forney, L. J., and Klug, M. J. (1994). Accuracy, reproducibility, and interpretation of fatty acid methyl ester profiles of model bacterial communities. Appl. Environ. Microbiol. 60:2483. doi: 10.1128/aem.60.7.2483-2493.1994

PubMed Abstract | Crossref Full Text | Google Scholar

Hart, M. M., Aleklett, K., Chagnon, P. L., Egan, C., Ghignone, S., Helgason, T., et al. (2015). Navigating the labyrinth: A guide to sequence-based, community ecology of arbuscular mycorrhizal fungi. New Phytol. 207, 235–247. doi: 10.1111/nph.13340

PubMed Abstract | Crossref Full Text | Google Scholar

Hayward, R. M., Banin, L. F., Burslem, D. F. R. P., Chapman, D. S., Philipson, C. D., Cutler, M. E. J., et al. (2021). Three decades of post-logging tree community recovery in naturally regenerating and actively restored dipterocarp forest in Borneo. For. Ecol. Manage. 488:119036. doi: 10.1016/j.foreco.2021.119036

Crossref Full Text | Google Scholar

Ihrmark, K., Bödeker, I. T. M., Cruz-Martinez, K., Friberg, H., Kubartova, A., Schenck, J., et al. (2012). New primers to amplify the fungal ITS 2 region – Evaluation by 454-sequencing of artificial and natural communities. FEMS Microbiol. Ecol. 82, 666–677. doi: 10.1111/j.1574-6941.2012.01437.x

PubMed Abstract | Crossref Full Text | Google Scholar

Jesus, E. D. C., Marsh, T. L., Tiedje, J. M., and Moreira, F. M. D. S. (2009). Changes in land use alter the structure of bacterial communities in Western Amazon soils. ISME J. 3, 1004–1011. doi: 10.1038/ismej.2009.47

PubMed Abstract | Crossref Full Text | Google Scholar

Keenan, R. J., Reams, G. A., Achard, F., De Freitas, J. V., Grainger, A. and Lindquist, E. (2015). Dynamics of global forest area: Results from the FAO Global Forest Resources Assessment 2015. For. Ecol. Manage. 352, 9–20. doi: 10.1016/j.foreco.2015.06.014

Crossref Full Text | Google Scholar

Kerfahi, D., Tripathi, B. M., Lee, J., Edwards, D. P., and Adams, J. M. (2014). The impact of selective-logging and forest clearance for oil palm on fungal communities in Borneo. PLoS One 9:e111525. doi: 10.1371/journal.pone.0111525

PubMed Abstract | Crossref Full Text | Google Scholar

Kozich, J., Schloss, P., Baxter, N., Jenior, M., Koumpouras, C., and Bishop, L. (2019). 16S rRNA Sequencing with the Illumina MiSeq: Library Generation, QC, & Sequencing - Version 6.0. Available online at: https://github.com/SchlossLab/MiSeq_WetLab_SOP/blob/master/MiSeq_WetLab_SOP.md (accessd Auguest 8, 2019).

Google Scholar

Kozich, J. J., Westcott, S. L., Baxter, N. T., Highlander, S. K., and Schloss, P. D. (2013). Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl. Environ. Microbiol. 79:5112. doi: 10.1128/AEM.01043-13

PubMed Abstract | Crossref Full Text | Google Scholar

Kumagai, T. O., and Porporato, A. (2012). Drought-induced mortality of a Bornean tropical rain forest amplified by climate change. J. Geophys. Res. 117:G02032. doi: 10.1029/2011JG001835

Crossref Full Text | Google Scholar

Landon, J. R. (1984). Booker tropical soil manual: A handbook for soil survey and agricultural land evaluation in the tropics and subtropics. London: Booker Agriculture International.

Google Scholar

Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009). Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75:5111. doi: 10.1128/AEM.00335-09

PubMed Abstract | Crossref Full Text | Google Scholar

Lechevalier, H., and Lechevalier, M. P. (1988). “Chemotaxonomic use of lipids - An overview,” in Microbial Lipids, eds C. Ratledge and S. G. Wilkinson (Cambridge, MA: Academic Press), 869–902.

Google Scholar

Lee-Cruz, L., Edwards, D. P., Tripathi, B. M., and Adams, J. M. (2013). Impact of logging and forest conversion to oil palm plantations on soil bacterial communities in Borneo. Appl. Environ. Microbiol. 79:7290. doi: 10.1128/AEM.02541-13

PubMed Abstract | Crossref Full Text | Google Scholar

Legendre, P., and Borcard, D. (2018). Box-Cox-chord transformations for community composition data prior to beta diversity analysis. Ecography 41, 1820–1824. doi: 10.1111/ecog.03498

Crossref Full Text | Google Scholar

Lekberg, Y., Vasar, M., Bullington, L. S., Sepp, S. K., Antunes, P. M., Bunn, R., et al. (2018). More bang for the buck? Can arbuscular mycorrhizal fungal communities be characterized adequately alongside other fungi using general fungal primers? New Phytol. 220, 971–976. doi: 10.1111/nph.15035

PubMed Abstract | Crossref Full Text | Google Scholar

Lenth, R. S., Henrik, S., Love, J., Buerkner, P., and Herve, M. (2019). Estimated Marginal Means, aka Least-Squares Means. Available online at: https://CRAN.R-project.org/package=emmeans (accessed July 11, 2025).

Google Scholar

Marsh, C. J., Turner, E. C., Wong Blonder, B., Bongalov, B., Both, S., Cruz, R. S., et al. (2025). Tropical forest clearance impacts biodiversity and function, whereas logging changes structure. Science 387, 171–175. doi: 10.1126/science.adf9856

PubMed Abstract | Crossref Full Text | Google Scholar

Marsh, C. W., and Greer, A. G. (1992). Forest land-use in Sabah, Malaysia: An introduction to Danum Valley. Philos. Trans. R. Soc. B Biol. Sci. 335, 331–339. doi: 10.1098/rstb.1992.0025

Crossref Full Text | Google Scholar

Martiny, J. B. H., Bohannan, B. J. M., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., et al. (2006). Microbial biogeography: Putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. doi: 10.1038/nrmicro1341

PubMed Abstract | Crossref Full Text | Google Scholar

McGuire, K., D’Angelo, H., Brearley, F., Gedallovich, S., Babar, N., Yang, N., et al. (2015). Responses of soil fungi to logging and oil palm agriculture in Southeast Asian tropical forests. Microb. Ecol. 69, 733–747. doi: 10.1007/s00248-014-0468-4

PubMed Abstract | Crossref Full Text | Google Scholar

McGuire, K., Henkel, T., Granzow de la Cerda, I., Villa, G., Edmund, F., and Andrew, C. (2008). Dual mycorrhizal colonization of forest-dominating tropical trees and the mycorrhizal status of non-dominant tree and liana species. Mycorrhiza 18, 217–222. doi: 10.1007/s00572-008-0170-9

PubMed Abstract | Crossref Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2013). phyloseq: An R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217

PubMed Abstract | Crossref Full Text | Google Scholar

Melo, F. P. L., Arroyo-Rodríguez, V., Fahrig, L., Martínez-Ramos, M., and Tabarelli, M. (2013). On the hope for biodiversity-friendly tropical landscapes. Trends Ecol. Evol. 28, 462–468. doi: 10.1016/j.tree.2013.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

Morlon, H., Chuyong, G., Condit, R., Hubbell, S., Kenfack, D., Thomas, D., et al. (2008). A general framework for the distance–decay of similarity in ecological communities. Ecol. Lett. 11, 904–917. doi: 10.1111/j.1461-0248.2008.01202.x

PubMed Abstract | Crossref Full Text | Google Scholar

Moura Costa, P. (1996). Tropical forestry practices for carbon sequestration: A review and case study from Southeast Asia. Ambio 25, 279–283.

Google Scholar

Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A., and Kent, J. (2000). Biodiversity hotspots for conservation priorities. Nature 403:853. doi: 10.1038/35002501

PubMed Abstract | Crossref Full Text | Google Scholar

Nekola, J. C., and White, P. S. (1999). The distance decay of similarity in biogeography and ecology. J. Biogeogr. 26, 867–878. doi: 10.1046/j.1365-2699.1999.00305.x

Crossref Full Text | Google Scholar

Nguyen, N. H., Song, Z., Bates, S. T., Branco, S., Tedersoo, L., Menke, J., et al. (2016). FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 20, 241–248. doi: 10.1016/j.funeco.2015.06.006

Crossref Full Text | Google Scholar

Nurulita, Y., Adetutu, E. M., Gunawan, H., Zul, D., and Ball, A. S. (2016). Restoration of tropical peat soils: The application of soil microbiology for monitoring the success of the restoration process. Agricult. Ecosyst. Environ. 216, 293–303. doi: 10.1016/j.agee.2015.09.031

Crossref Full Text | Google Scholar

O’ Brien, S. L., Gibbons, S. M., Owens, S. M., Hampton-Marcell, J., Johnston, E. R., Jastrow, J. D., et al. (2016). Spatial scale drives patterns in soil bacterial diversity. Environ. Microbiol. 18, 2039–2051. doi: 10.1111/1462-2920.13231

PubMed Abstract | Crossref Full Text | Google Scholar

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., O’Hara, R. B., Simpson, G. L., et al. (2019). Vegan: Community Ecology Package. Available online at: https://cran.r-project.org/package=vegan (accessed May 6, 2025).

Google Scholar

O’Leary, W. M., and Wilkinson, S. G. (1988). “Gram-positive bacteria,” in Microbial lipids, eds C. Ratledge and S. G. Wilkinson (Cambridge, MA: Academic Press), 117–201.

Google Scholar

Pan, Y., Birdsey, R. A., Fang, J., Houghton, R., Kauppi, P. E., Kurz, W. A., et al. (2011). A large and persistent carbon sink in the world’s forests. Science 333, 988–993. doi: 10.1126/science.1201609

PubMed Abstract | Crossref Full Text | Google Scholar

Paul, M., Catterall, C. P., Pollard, P. C., and Kanowski, J. (2010). Recovery of soil properties and functions in different rainforest restoration pathways. For. Ecol. Manage. 259, 2083–2092. doi: 10.1016/j.foreco.2010.02.019

Crossref Full Text | Google Scholar

Perumal, M., Wasli, M. E., and Lat, J. (2016). “Effects of enrichment planting on the soil physicochemical properties at reforestation sites planted with Shorea macrophylla de Vriese in Sampadi Forest Reserve, Sarawak, Malaysia,” in Proceedings of the symposium “Frontier in tropical forest research: Progress in joint projects between the Forest Department Sarawak and the Japan Research Consortium for Tropical Forests in Sarawak”, (Kuching).

Google Scholar

Petersen, I. A. B., Meyer, K. M., and Bohannan, B. J. M. (2019). Meta-analysis reveals consistent bacterial responses to land use change across the tropics. Front. Ecol. Evol. 7:391. doi: 10.3389/fevo.2019.00391

Crossref Full Text | Google Scholar

Philipson, C. D., Cutler, M. E. J., Brodrick, P. G., Asner, G. P., Boyd, D. S., Moura Costa, P., et al. (2020). Active restoration accelerates the carbon recovery of human-modified tropical forests. Science 369, 838–841. doi: 10.1126/science.aay4490

PubMed Abstract | Crossref Full Text | Google Scholar

Potapov, P., Hansen, M. C., Laestadius, L., Turubanova, S., Yaroshenko, A., Thies, C., et al. (2017). The last frontiers of wilderness: Tracking loss of intact forest landscapes from 2000 to 2013. Sci. Adv. 3:e1600821. doi: 10.1126/sciadv.1600821

PubMed Abstract | Crossref Full Text | Google Scholar

Potapov, P., Yaroshenko, A., Turubanova, S., Dubinin, M., Laestadius, L., Thies, C., et al. (2008). Mapping the world’s intact forest landscapes by remote sensing. Ecol. Soc. 13:51. doi: 10.5751/ES-02670-130251

PubMed Abstract | Crossref Full Text | Google Scholar

Powers, R. P., and Jetz, W. (2019). Global habitat loss and extinction risk of terrestrial vertebrates under future land-use-change scenarios. Nat. Clim. Change 9, 323–329. doi: 10.1038/s41558-019-0406-z

Crossref Full Text | Google Scholar

Qie, L., Lewis, S., Sullivan, M., Lopez-Gonzalez, G., Pickavance, G., Sunderland, T., et al. (2017). Long-term carbon sink in Borneo’s forests halted by drought and vulnerable to edge effects. Nat. Commun. 8, 1966–1966. doi: 10.1038/s41467-017-01997-0

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team. (2025). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Reynolds, G., Payne, J., Sinun, W., Mosigil, G., and Walsh, R. P. D. (2011). Changes in forest land use and management in Sabah, Malaysian Borneo, 1990–2010, with a focus on the Danum Valley region. Philos. Trans. R. Soc. B 366, 3168–3176. doi: 10.1098/rstb.2011.0154

PubMed Abstract | Crossref Full Text | Google Scholar

Rinnan, R., and Bååth, E. (2009). Differential utilization of carbon substrates by bacteria and fungi in tundra soil. Appl. Environ. Microbiol. 75:3611. doi: 10.1128/AEM.02865-08

PubMed Abstract | Crossref Full Text | Google Scholar

Riutta, T., Malhi, Y., Kho, L. K., Marthews, T. R., Huaraca Huasco, W., Khoo, M., et al. (2018). Logging disturbance shifts net primary productivity and its allocation in Bornean tropical forests. Glob. Change Biol. 24, 2913–2928. doi: 10.1111/gcb.14068

PubMed Abstract | Crossref Full Text | Google Scholar

Robinson, S. J. B., Elias, D., Johnson, D., Both, S., Riutta, T., Goodall, T., et al. (2020). Soil fungal community characteristics and mycelial production across a disturbance gradient in lowland dipterocarp rainforest in Borneo. Front. For. Glob. Change 3:64. doi: 10.3389/ffgc.2020.00064

Crossref Full Text | Google Scholar

Robinson, S. J. B., Elias, D. M. O., Goodall, T., McNamara, N., Griffiths, R., Majalap, N., et al. (2025). Soil microbial communities, soil, environmental and forest structural properties across old-growth forest and naturally regenerating or actively restored selectively logged tropical forest in Sabah, Borneo, March 2018 v2. Lancaster: NERC EDS Environmental Data Centre. doi: 10.5285/08bfe302-d33f-490d-97be-27bb83a0f38d

Crossref Full Text | Google Scholar

Robinson, S. J. B., Elias, D. M. O., Goodall, T., Nottingham, A. T., McNamara, N. P., Griffiths, R., et al. (2024). Selective logging impacts on soil microbial communities and functioning in Bornean tropical forest. Front. Microbiol. 15:1447999. doi: 10.3389/fmicb.2024.1447999

PubMed Abstract | Crossref Full Text | Google Scholar

Rodrigues, J. L. M., Pellizari, V. H., Mueller, R., Baek, K., Jesus, E. D. C., Paula, F. S., et al. (2013). Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities. Proc. Natl. Acad. Sci. U. S. A. 110, 988–993. doi: 10.1073/pnas.1220608110

PubMed Abstract | Crossref Full Text | Google Scholar

Rodrigues, R., Lima, R., Gandolfi, S., and Nave, A. (2009). On the restoration of high diversity forests: 30 years of experience in the Brazilian Atlantic Forest. Biol. Conserv. 142, 1242–1251. doi: 10.1016/j.biocon.2008.12.008

Crossref Full Text | Google Scholar

Schnitzer, S. A., Kuzee, M. E., and Bongers, F. (2005). Disentangling above- and below-ground competition between lianas and trees in a tropical forest. J. Ecol. 93, 1115–1125. doi: 10.1111/j.1365-2745.2005.01056.x

Crossref Full Text | Google Scholar

Shi, L., Dossa, G. G. O., Paudel, E., Zang, H., Xu, J., and Harrison, R. D. (2019). Changes in fungal communities across a forest disturbance gradient. Appl. Environ. Microbiol. 85:e00080-19. doi: 10.1128/AEM.00080-19

Crossref Full Text | Google Scholar

Shoo, L. P., Freebody, K., Kanowski, J., and Catterall, C. P. (2016). Slow recovery of tropical old-field rainforest regrowth and the value and limitations of active restoration. Conserv. Biol. 30, 121–132. doi: 10.1111/cobi.12606

PubMed Abstract | Crossref Full Text | Google Scholar

Smart, S. M., Thompson, K., Marrs, R. H., Le Duc, M. G., Maskell, L. C., and Firbank, L. G. (2006). Biotic homogenization and changes in species diversity across human-modified ecosystems. Proc. R. Soc. B Biol. Sci. 273, 2659–2665. doi: 10.1098/rspb.2006.3630

PubMed Abstract | Crossref Full Text | Google Scholar

Sodhi, N. S., Lee, T. M., Koh, L. P., and Brook, B. W. (2009). A meta-analysis of the impact of anthropogenic forest disturbance on Southeast Asia’s biotas. Biotropica 41, 103–109. doi: 10.1111/j.1744-7429.2008.00460.x

Crossref Full Text | Google Scholar

Stibig, H. J., Achard, F., Carboni, S., Raši, R., and Miettinen, J. (2014). Change in tropical forest cover of Southeast Asia from 1990 to 2010. Biogeosciences 11, 247–258. doi: 10.5194/bg-11-247-2014

Crossref Full Text | Google Scholar

Strickland, M. S., Callaham, M. A., Gardiner, E. S., Stanturf, J. A., Leff, J. W., Fierer, N., et al. (2017). Response of soil microbial community composition and function to a bottomland forest restoration intensity gradient. Appl. Soil Ecol. 119, 317–326. doi: 10.1016/j.apsoil.2017.07.008

Crossref Full Text | Google Scholar

Tin, H. S., Palaniveloo, K., Anilik, J., Vickneswaran, M., Tashiro, Y., Vairappan, C. S., et al. (2018). Impact of land-use change on vertical soil bacterial communities in Sabah. Microb. Ecol. 75, 459–467. doi: 10.1007/s00248-017-1043-6

PubMed Abstract | Crossref Full Text | Google Scholar

Tripathi, B., Kim, M., Singh, D., Lee-Cruz, L., Lai-Hoe, A., Ainuddin, A., et al. (2012). Tropical soil bacterial communities in Malaysia: pH dominates in the equatorial tropics too. Microb. Ecol. 64, 474–484. doi: 10.1007/s00248-012-0028-8

PubMed Abstract | Crossref Full Text | Google Scholar

Tripathi, B. M., Edwards, D. P., Mendes, L. W., Kim, M., Dong, K., Kim, H., et al. (2016). The impact of tropical forest logging and oil palm agriculture on the soil microbiome. Mol. Ecol. 25, 2244–2257. doi: 10.1111/mec.13620

PubMed Abstract | Crossref Full Text | Google Scholar

van der Heijden, M. G. A., Bardgett, R. D., and van Straalen, N. M. (2008). The unseen majority: Soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. doi: 10.1111/j.1461-0248.2007.01139.x

PubMed Abstract | Crossref Full Text | Google Scholar

van der Putten, W., Bardgett, R. D., Bever, J. D., Bezemer, T. M., Casper, B., Fukami, T., et al. (2013). Plant-soil feedbacks: The past, the present and future challenges. J. Ecol. 101, 265–276. doi: 10.1111/1365-2745.12054

Crossref Full Text | Google Scholar

Veryard, R., Wu, J., O’Brien, M. J., Anthony, R., Both, S., Burslem, D. F. R. P., et al. (2023). Positive effects of tree diversity on tropical forest restoration in a field-scale experiment. Sci. Adv. 9:eadf0938. doi: 10.1126/sciadv.adf0938

PubMed Abstract | Crossref Full Text | Google Scholar

Walsh, R. P. D., and Newbery, D. M. (1999). The ecoclimatology of Danum, Sabah, in the context of the world’s rainforest regions, with particular reference to dry periods and their impact. Philos. Trans. R. Soc. B Biol. Sci. 354, 1869–1883. doi: 10.1098/rstb.1999.0528

PubMed Abstract | Crossref Full Text | Google Scholar

Wardle, D., Bardgett, R. D., Klironomos, J., Setala, H., van der Putten, W., and Wall, D. (2004). Ecological linkages between aboveground and belowground biota. Science 304, 1629–1633. doi: 10.1126/science.1094875

PubMed Abstract | Crossref Full Text | Google Scholar

Wei, X., Shao, M., Gale, W., and Li, L. (2014). Global pattern of soil carbon losses due to the conversion of forests to agricultural land. Sci. Rep. 4:4062. doi: 10.1038/srep04062

PubMed Abstract | Crossref Full Text | Google Scholar

Whitaker, J., Ostle, N., McNamara, N. P., Nottingham, A. T., Stott, A. W., Bardgett, R. D., et al. (2014). Microbial carbon mineralization in tropical lowland and montane forest soils of Peru. Front. Microbiol. 5:720. doi: 10.3389/fmicb.2014.00720

PubMed Abstract | Crossref Full Text | Google Scholar

White, D., Davis, W., Nickels, J., King, J., and Bobbie, R. (1979). Determination of the sedimentary microbial biomass by extractible lipid phosphate. Oecologia 40, 51–62. doi: 10.1007/BF00388810

PubMed Abstract | Crossref Full Text | Google Scholar

White, T. J., Bruns, T., Lee, S., and Taylor, J. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics,” in PCR Protocols: A Guide to Methods and Applications, eds M. A. Innis, D. H. Gelfland, J. J. Sninsky, and T. J. White (San Diego, CA: Academic Press), 315–322.

Google Scholar

Wilkinson, D. M. (1999). The disturbing history of intermediate disturbance. Oikos 84, 145–147. doi: 10.2307/3546874

Crossref Full Text | Google Scholar

Wright, S. J. (2010). The future of tropical forests. Ann. N. Y. Acad. Sci. 11951, 1–27. doi: 10.1111/j.1749-6632.2010.05455.x

PubMed Abstract | Crossref Full Text | Google Scholar

Zelles, L. (1999). Fatty acid patterns of phospholipids and lipopolysaccharides in the characterisation of microbial communities in soil: A review. Biology and Fertility of Soils 29, 111–129. doi: 10.1007/s003740050533

Crossref Full Text | Google Scholar

Zhang, X., Johnston, E. R., Barberán, A., Ren, Y., Wang, Z., and Han, X. (2018). Effect of intermediate disturbance on soil microbial functional diversity depends on the amount of effective resources. Environ. Microbiol. 20, 3862–3875. doi: 10.1111/1462-2920.14407

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, X., Liu, W. E. I., Bai, Y., Zhang, G., and Han, X. (2011). Nitrogen deposition mediates the effects and importance of chance in changing biodiversity. Mol. Ecol. 20, 429–438. doi: 10.1111/j.1365-294X.2010.04933.x

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: 16S, ITS, ectomycorrhizal fungi, forest rehabilitation, liana cutting, climber cutting, liberation thinning, distance decay

Citation: Robinson SJB, Elias DMO, Goodall T, McNamara NP, Griffiths R, Majalap N, Sau Wai Y and Ostle NJ (2025) Soil microbial community responses to active and passive restoration of selectively logged Bornean tropical forest. Front. Microbiol. 16:1570294. doi: 10.3389/fmicb.2025.1570294

Received: 03 February 2025; Accepted: 16 July 2025;
Published: 22 August 2025.

Edited by:

Antonio Castellano-Hinojosa, Universidad de Granada, Spain

Reviewed by:

Rosa María Pineda-Mendoza, Instituto Politécnico Nacional (IPN), Mexico
Juan P. Frene, New Mexico State University, United States

Copyright © 2025 Robinson, Elias, Goodall, McNamara, Griffiths, Majalap, Sau Wai and Ostle. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Samuel J. B. Robinson, c2Ftcm9iMzlAY2VoLmFjLnVr

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.