Low- to Moderate-Level Forest Disturbance Effects on Plant Functional Traits and Associated Soil Microbial Diversity in Western Himalaya

This study quantifies the effect of low- to moderate-intensity anthropogenic activities and natural activities affecting less than 40% of the canopy cover on major plant functional traits and associated soil microbial diversity in western Himalayan temperate forests. Plots of 0.1 ha were placed in the temperate forests along disturbance gradient and were categorized based on lopping intensity, deadwood counts, grazing, and litter removal. Plots were classified into three classes, i.e., low disturbance intensity (LDI) and moderate disturbance intensities (MDI1 and MDI2) based on lopping activities and canopy cover. The study was conducted on functional traits related to growth and survival strategies of a species in a complex forest ecosystem. Furthermore, DNA was extracted and metagenome of soil samples was performed using Illumina MiSeq platform from three disturbance classes to study the effect of disturbance and plant traits on microbial diversity. Tree basal area was found to be most significantly affected by disturbance intensity. Total density was found to be greater for the LDI site. Specific leaf area (SLA) and crown cover (CC) were the most affected traits in the moderately disturbed sites, whereas maximum plant height (HT) and seed mass (SM) were least affected by disturbance. Soil microbial diversity was found to be negatively associated with disturbance index. Microbial biomass carbon (MBC) and microbial biomass nitrogen (MBN) were found to be significantly higher in LDI sites. Disturbance level elicited variation in microbial community composition. Proteobacteria was the most abundant phylum and Phenylobacterium, DA101, and Candidatus solibacter were mainly abundant at the genus level. A decreasing level of disturbance due to the absence of human residences in the LDI site led to the dominance of Phenylobacterium (27%) which reduced to 18% in MDI1 site. Shannon alpha bacterial diversity and plant species diversity (H′) were found to be greatest for MDI2 site. In forests with varying levels of management, treefall gaps due to low levels of logging intensity might have a similar effect to those of reduced lopping intensities and deadwood count (due to natural disturbances) in the study sites. The study concludes that moderate disturbance is important for promoting species diversity and species richness (SR), but species having conservative ecological strategies would be more prone to continued disturbance intensity. Therefore, low- to moderate-level disturbance in such forests can be used as a model for natural treefall gaps, and moderate-level disturbance intensity plays a powerful role in buffering ecosystem processes.


INTRODUCTION
Natural and anthropogenic forms of disturbances are known to alter forest structure and functions, but variation in disturbance intensity and frequency causes variable outcomes (Atkins et al., 2020). The anthropogenic forms of disturbances have abruptly increased with increase in population and urbanization and increased demand for land and forest resources causing adverse effects on climate changes. Anthropogenic disturbances such as overgrazing, logging, and litter removal directly or indirectly affect the biodiversity and ecosystem processes (Heywood and Watson, 1995). Most of the human populations reside in the range of 1,000-2,000 masl in the Indian Himalayan region (Singh, 1998). The forest ecosystem in these regions is severely impacted by human activities such as deforestation, lopping for fuel and fodder, litter removal, and frequent forest fire. While there have been numerous studies on the impact of long-term chronic disturbances on forest ecosystem in the Indian Himalayan region (Thakur et al., 2020;Kaushal et al., 2021), there has been limited evaluation of moderate to low severity disturbances (e.g., lopping that is known to alter canopy cover and litter removal that causes denudation of the ground) in the region. Low to moderate levels of disturbances could lead to a wide range of compositional and structural variation throughout the region (Fahey et al., 2015), therefore have different consequences for ecosystem functioning and services (Nave et al., 2011).
Functional traits could help in better understanding the variation of species responses to external factors (McGill et al., 2006). Specific leaf area (SLA) is known to define resource allocation strategies for acquisition and storage, thus reflecting important trade-offs between construction cost and survival under varied conditions (Poorter et al., 2008Iida et al., 2014). Also, leaves with low SLA exhibit high thickness and toughness that suggests increased resource allocation into defense than into growth (Wright and Cannon, 2001). This leads to greater leaf life span (Sterck et al., 2006), thus affecting overall leaf economic spectrum. Such species also tend to produce small seeds having less generation time (Salguero-Gómez et al., 2016). Plant functional traits are also known to be an important predictor to assess the interaction between aboveground and belowground ecosystem processes (Grigulis et al., 2013;Cantarel et al., 2015;Moreau et al., 2015). Plant diversity along with functional traits such as plant height (HT), basal area, and SLA is known to affect soil bacterial communities (Pei et al., 2016), which directly or indirectly affects belowground ecosystem processes (Cao et al., 2017) specifically in the mountainous soils (Singh and Gupta, 2018). Carbon and nitrogen intake by variation in plant litter and fine root composition are known to further influence the soil microbial biomass. The community composition altered by disturbances can further affect soil abiotic properties by inducing changes in microbial community characteristics (Bagchi et al., 2014;Delgado-Baquerizo et al., 2018) and vice versa (Demenois et al., 2018).
Disturbance can cause variation in plant functional trait composition but not community composition or vice versa (Mayfield et al., 2010). Studies in the Indian Himalayan region have focused on the effect of disturbance on species composition (Kumar and Ram, 2005;Rawal et al., 2012), but variations in functional traits have limited information. Also, despite such insights and awareness, studies linking varied levels of disturbances and plant community composition and soil microbial communities are limited in the Indian Himalayan regions (Rawat et al., 2020). The studies majorly focus on species composition and diversity; hence, the role of plant traits and community composition on soil microbial community composition remains unclear. Thus, inclusion of several plant traits and plant diversity might play an important role in governing soil bacterial diversity and composition.
This study explains the role of low-to moderate-intensity disturbance activities on plant functional traits and associated soil microbial diversity in the temperate forests of western Himalaya. Driven by the aforementioned facts and circumstances, a systematic study in reserved temperate forests of western Himalaya was conducted considering the following broad hypotheses: (i) variation in disturbance intensity would significantly impact community characteristics, plant traits, and associated microbial diversity in the studied region; and (ii) species richness (SR) and diversity will be greater following intermediate disturbance hypothesis. The present study is based upon the following objectives: (i) to observe the variation in different community characteristics across varying intensity of disturbances, (ii) to study the patterns of change in the plant functional traits and soil microbial diversity with disturbance intensity, and (iii) to investigate consistent patterns in association of soil microbial diversity and plant factors with variation in disturbance intensity.

Study Site
The study site was located in Sawana Reserved Forest which lies in Rajgarh range of Sirmour District, Himachal Pradesh (Figure 1). Its geographical range lies within 30 • 86 N-31 • 01 N and 77 • 38 E and 77 • 49 E. It covers an area of 1,395.38 ha with an altitudinal variation of 1,537-3,164 masl. Climate varies from moderate and cold to very cold at higher elevations. Temperature range from 10 to 30 • C and rainfall is about 125 cm annually (District Survey Report, 2016). The forest is classified under group 12-Himalayan moist-temperate forest as per the forest classification of Champion and Seth (1968). These forests occur in pure as well as mixed forms. The soil type is clayey loam type. Granite, gneiss, micaceous, schist, and quartzite are the main rock types. Moderate to steep slope occurs throughout the region. Lopping and grazing remain to be the most important primary disturbance factors that affect the plant species diversity in the study region.

Plot Selection
The reserved forest was classified into three different categories of disturbance based on varied anthropogenic intensity and human settlements in the region, i.e., low disturbance intensity site (LDI) and moderately disturbed intensity sites (MD1 and MD2) ( Table 1). Eighteen plots of 0.1 ha (31 m × 32 m) were established along the forest to sample live trees having diameter at breast height (DBH) ≥ 10 cm. Six plots were established in each disturbed site to account for better representation within each category. Plant traits were sampled from each site to study the effect of low-to moderate-level disturbances on plant growth and associated soil microbial diversity in the region. Plant functional traits were selected given that they have an impact when grown in disturbed sites. Ocular estimation based on field survey suggested deadwood counts in the region to be associated with natural disturbances, whereas anthropogenic disturbances were mainly due to human disturbances and livestock population. However, livestock populations are known to severely impact canopy species in the area. This was beyond the scope of this study because the major focus was only for tree layer (mature individuals). Lopping, defined as cutting off branches from the main stem of the tree, could be seen as the most common anthropogenic activity in the region. Lopped branches are generally used for fuel and fodder purposes. Lopping intensity was calculated in each plot by counting the lopped tree individuals against total tree individuals. The sites were also classified based on canopy cover. Canopy cover in the forest was measured using a densiometer. The categorization of canopy cover was done following the Forest Survey of India (2015). Grazing intensities were estimated based on the number of livestock dungs found in each plot during the survey, and litter removal was calculated as the percentage of bare land to the total ground cover in each plot. The information regarding different  Disturbances to microbial biomass may be attributable to differences in soil disruption and organic carbon removal from forests among disturbance types (Holden and Treseder, 2013) Microbial biomass nitrogen (MBN) µg g −1 270 disturbance activities was collected following Rawal et al. (2012). Deadwood counts (as a source of natural disturbance) were categorized as veteran tree, standing dead tree, wind-thrown tree, fallen deadwood tree, and stump in each plot against the total number of individuals (Humphrey et al., 2002).

Floristic Survey and Trait Measurements
Tree species in the study sites were identified with the help of the glossary of forest flora provided by the Rajgarh Forest Department, Sirmour. CBH (circumference at breast height) was measured during the field survey which was later converted into DBH for determination of tree basal area. The data were quantitatively analyzed to calculate stem density, frequency, and abundance following Curtis and McIntosh (1950) for the three disturbance categories. Similarly, the Shannon-Wiener diversity index (H ) (Shannon and Weaver, 1963), species evenness (Pielou, 1966), and SR (Margalef, 1958) were also calculated for each site. All species that contributed to 90% of plot basal area (in decreasing order of species basal area) were identified for sampling. In each plot, 10 individuals per species were sampled to study all the FTs chosen for the study. Sampling was done from April to June 2019 on 287 mature tree individuals. Selection of traits has been done following Westoby (1998) and Grime (1977) who suggest summarizing plant strategy through three major axes: resource exploitation, competition ability, and response to disturbance. Leaf traits represent resource use, and plant height and crown cover (CC) for competitive ability and seed mass embody the response to perturbations (Westoby, 1998). Table 2 represents the rationale behind the selection of plant traits and associated variables. Plant height was measured using a clinometer. CC for each selected tree individual was calculated using the line intercept method (USDA, 1996). Samples were collected aiming for leaves fully expanded with little or no damage by herbivores or pathogens. Fresh weight was taken by a portable weighing machine after drying through gently blotting with tissue paper to remove any surface water. After sampling, the leaves were air dried and transported in sealed plastic bags to the laboratory. SLA and SM were calculated following Cornelissen et al. (2003).

Soil Sampling
Five randomly located soil core samples at depths of 0-10, 10-20, and 20-30 cm were collected from each plot using a soil corer of 5 cm diameter. Soil microbial carbon and soil microbial nitrogen were estimated following the chloroform fumigation-extraction method (Brookes et al., 1984;Vance et al., 1987). Soil samples from each disturbed class were collected to make composite samples (from varying depths) that were sent for DNA extraction and purification in the laboratory.

DNA Extraction and Purification
Metagenomics sequencing was performed on the collected soil sample (0.25 g) to evaluate the relationship between microbial diversity and plant functional traits. DNA extraction was done using commercially available kits from Qiagen (Hilden, Germany), Zymo Research (Irvine, CA, United States), and Thermo Fisher (Waltham, MA, United States) followed by PCR amplification of the V3-V4 region. Forty nanograms of extracted DNA were used for amplification along with 10 pM of each primer. Twenty-five cycles each of initial denaturation were done at 95 • C for 15 s, annealing at 60 • C for 15 s, elongation at 72 • C for 2 min, final extension at 72 • C for 10 min, and hold at 4 • C. The amplicons from each sample were purified with Ampure beads to remove unused primers and eight additional cycles of PCR were performed using Illumina barcoded adapters to prepare the sequencing libraries. Libraries were purified using Ampure beads and quantitated using Qubit dsDNA High Sensitivity assay kit. Sequencing was performed using Illumina MiSeq with 2 × 300 PE v3 sequencing kit. Raw data QC was done using FASTQC and MULTIQC, followed by trimming of adapters and low-quality reads by TRIMGALORE. The trimmed reads are further taken for processing which includes merging of paired end reads, chimera removal, and operational taxonomic unit (OTU) abundance calculation and estimation correction-this is achieved by QIIME/MOTHUR/KRAKEN/BRACKEN workflows (Supplementary Table 1). This workflow enables highly accurate investigations at the genus level. Each read was classified based on percent coverage and identity. Extracted DNA from the samples was subjected to NanoDrop and GEL Check before being taken for PCR amplification: the NanoDrop readings of 260/280 at an approximate value of 1.8-2 were used to determine the quality of DNA. The amplified 16S PCR product was purified and subjected to GEL Check and NanoDrop QC. The Nano Drop readings of 260/280 at an approximate value of 1.8-2 were used to determine the quality of DNA.

Statistical Analysis
Data normality was checked with Shapiro-Wilk's test and homogeneity of variance of data with Levene's test. Kernel density plots were drawn for each plant trait and soil variables to understand the variability within each disturbed site. Oneway ANOVA followed by Tukey post hoc test was conducted to evaluate the effect of disturbance on plant functional traits, community composition, and associated soil microbial diversity, and the differences were reported significant at p < 0.001. Principal component analysis (PCA) was performed including the plant traits for 287 individuals in the study sites and soil properties from each plot. The importance value index of plant species was calculated following Curtis (1959), and its association with different plant traits and soil parameters was determined using canonical correspondence analysis (CCA). CCA was also used to explore changes in traits and diversity with forest disturbances. Pearson correlations were performed to examine the relationships between soil microbial communities, PFTs, and plant community composition. All statistical analyses were conducted using R version 3.5.0 (R Core Team, 2018) using the packages ggplot, vegan, ggbiplot, imputeTs, and corrplot respectively.

Changes in Plant Community Composition and Functional Traits With Varying Disturbance Intensity
After phytosociological survey, 22 tree species having 1,220 individuals were found throughout the study site. Total basal area and tree density were found to be maximum for the LDI site (71.45 m 2 ha −1 and 834.16 individuals ha −1 ) and least for the MDI1 site (22.86 m 2 ha −1 and 546.12 individuals ha −1 ), respectively. Total basal area of mature individuals was found to be significantly varying among the study sites (Table 3). However, the differences were not significant for tree density. Plant SR varied from 0.88 to 2.09 being maximum for the MDI2 site and least for the LDI site (Supplementary Table 2). Species diversity (H ) did not vary much throughout the study sites. However, values for H were also greater for moderately disturbed sites. Values for plant traits (CC and SLA) varied significantly among the study sites, whereas there was no significant difference found for HT and SM (Table 3). HT as a function of growth and competitive vigor was found to be maximum for the MDI1 site and least for the LDI site. CC was found to be varying throughout the study site due to the effect of lopping. SM showed greater variation within the study sites reflecting fundamental divergence in the phylogenetic distance between deciduous and evergreen species (Figure 2  were found. Analysis of soil bacterial community composition revealed sequences belonging to 19 phyla, 44 classes, 70 orders, 127 families, and 221 genera in soil samples across all the stands. At the phylum level, the soil bacterial communities were found to be dominated by Proteobacteria, Acidobacteria, Actinobacteria, Firmicutes, Planctomycetes, Verrucomicrobia, Bacteroidetes, and Cyanobacteria together accounting for more than 80% of the bacterial community composition (Figure 3). Furthermore, Nitrospirae, Gemmatimonadetes, Spirochetes, Chlamydiae, Nitrospirae, Chlorofelxi, Planctomycetes, and Fusobacteria were also present in the soil samples with varying proportions across the study sites. Proteobacteria was the most abundant phylum and Phenylobacterium, DA101, and Candidatus solibacter were mainly abundant at the genus level. Decreasing the level of disturbance with increasing elevation due to the absence of human residences in higher elevation led to the dominance of Phenylobacterium (27%) that reduced to 18% in moderately disturbed sites. The alpha diversity indices (Shannon and Simpson) revealed similar patterns to plant species diversity and richness, i.e., the values were found to be maximum for the MDI2 site and lowest for the LDI site (Supplementary Table 2). Phylotypes within Actinobacteria and Acidobacteria were found to be greater for moderately disturbed sites, whereas phylotypes within Proteobacteria, Planctomycetes, Verrucomicrobia, and Bacteroidetes were greater for LDI sites. Soil microbial biomass carbon and nitrogen were found to increase around twofold with values being significantly greater for the LDI site and least for the MDI1 site (Figure 2). One-way ANOVA test revealed that soil microbial biomass carbon and nitrogen were the most affected variables with increase in disturbance intensity ( Table 3).

Relationship Between Plant Community Characteristics and Soil Microbial Diversity
Both plant species diversity and bacterial species diversity followed the same pattern being maximum for the MDI2 site and least for the LDI site. Pearson correlation coefficients calculated for plant traits, community characteristics, and soil microbial diversity along the disturbance gradient revealed that total basal area (TBA) and tree density (DEN) were found to have a positive correlation among the study sites (Figure 4). They were also found to have a strong positive correlation with soil MBC and MBN and bacterial phyla Proteobacteria, Planctomycetes, and Verrucomicrobia and negatively correlated with phyla Actinobacteria, Firmicutes, and Acidobacteria. Phyla Proteobacteria, Planctomycetes, and Verrucomicrobia were found to have a strong positive correlation with soil MBC and MBN, whereas phyla Actinobacteria, Firmicutes, and Acidobacteria were found to have a negative correlation with soil MBC and MBN. Similarly, CC was found to have negative correlation with bacterial phyla Actinobacteria, Firmicutes, and Acidobacteria and positively correlated with soil MBC and MBN and bacterial phyla Proteobacteria, Planctomycetes, Bacteroidetes, and Verrucomicrobia. Principal component analysis suggested that the first and second principal components explained 36.93 and 30.74% of data variability, respectively, while the third component accounted for 15.42% of data variability (Figure 5). The first PCA axis was associated with SLA and SM and opposed by CC, whereas the second axis was associated with MBC and MBN. HT was found to be negatively correlated with SM. CCA explained the relationship of plant community characteristics with soil variables in all the disturbed sites ( Figure 6). MBN, CC, and SLA were found to be the major contributing factors in the LDI and MDI2 sites, respectively. Eigenvalues of CCA axis 1 and CCA axis 2 were 0.80 and 0.46, respectively. Variances explained by CCA axis 1 and CCA axis 2 were was 42 and 24%, respectively. Correlation between CCA axis 1 and CC was found to be positive with CC (0.66) and MBN (0.92) and negative with SLA (−0.21), SM (−0.12), and HT (−0.08).

DISCUSSION
In the present study, the effect of low-to moderate-level disturbances on plant community characteristics and functional traits along with soil microbial diversity in temperate reserved forests in the western Himalayan region was explored. Tree basal area was found to be significantly varying throughout the study sites (p < 0.01). CC (p < 0.01) and SLA (p < 0.001) along with MBC (p < 0.001) and MBN (p < 0.001) were found to be significantly varying among the study sites. However, there was no effect on HT and SM based on disturbance intensity. The findings indicate that low to moderate level of disturbance does not affect the community composition, but with increasing and continued disturbance intensity, species with conservative growth strategies promoting resistance of species to environmental stresses will be replaced with species having higher trait values for acquisitive traits such as SLA, indicating higher competitive ability and faster resource acquisition. Here, the ecological significance of the present findings and their proper implications for sustainable forest management are discussed.

Changes in Plant Community Composition and Functional Traits With Varying Disturbance Intensity
Neither species diversity (H ) nor any of the plant functional traits were significantly affected by low to moderate levels of disturbance intensity. This holds for the hypothesis suggesting that SR and diversity will not be affected because of limited disturbance activities throughout the study sites. This lack of effect was probably due to regular management implications that do not allow the local people to barely denude the tree individuals which affects their growth and survival. However, SR was found to be increased by more than twofold in MDI2 sites that could be attributed to differences in management practices (Toledo et al., 2012) as well as recruitment of new species in moderately disturbed sites due to scattered canopy (Kumar, 2000). Studies have proven that mild disturbances provide greater opportunity for species turnover, establishment, and colonization, thus leading to higher SR (Whittaker and Likens, 1975;Connell, 1978;Mishra et al., 2004). Additionally, other factors might contribute to the differences among studies, such as topographical heterogeneity, environmental conditions, and forest types in the study site. Local disturbance causes tree gaps resulting in scattered canopy that creates patches of greater light availability and availability of optimum resources and nutrients for lower subcanopy and regenerating species (Raghubanshi and Tripathi, 2009). The findings tend to approach closer to the intermediate disturbance hypothesis theory which suggests an increase in species diversity with a low to moderate amount of disturbance (Connell, 1978). The hypothesis that the increasing level of lopping intensity impacts forest community characteristics holds for tree basal area that tends to decrease with increase in disturbance intensities but remains insignificant for tree density and other community characteristics. Similarly, low to moderate levels of disturbances caused no effect on species composition and diversity in Bornean rainforest (Verburg and van Eijk-Bos, 2003), tropical rainforest (Carreño-Rocabado et al., 2012), and semiarid forest (Mayfield et al., 2010). Varying levels of disturbance intensity had a significant impact on SLA and CC, whereas there was no effect on HT and SM. With an increase in disturbance intensity, SLA was found to be maximum at MDI2 sites suggesting shifting to acquisitive traits rather than conservative traits (Kitajima and Poorter, 2010;Onoda et al., 2011), whereas for HT and SM, low to moderate disturbance promoted species with conservative traits rather than acquisitive (Loehle and Namkoong, 1987). Similar results were found in a study conducted by Hall et al. (2020) in European vineyards indicating higher photosynthetic capacity and growth rates in bare soil vineyards. This suggests a faster turnover of plants as a response to disturbance (Kazakou et al., 2016). LDI sites were found to have the lowest value for SLA, and such individuals tend to have longer leaf life spans causing lower growth rates (Louault et al., 2005). SM was not found to be significantly affected by varying disturbance intensity, but within site, variation was mainly due to phylogenetic distances between the species. Conifers in the study site are being characterized by a higher number of seeds with lower seed mass as compared with broadleaved trees that are known to have increased seed mass. CC is considered to be another trait that is highly affected by disturbance activities, especially lopping. Variation in CC further affects light intensity variation in the study site as well loss of structural components, i.e., branches and photosynthetic organs, i.e., leaves (Sterck et al., 2001). CCs are also known to be affected by topographic and edaphic conditions (Lang, 2010). HT did not vary much with disturbance but was found to be lowest for the LDI site being less light demanding hence tends to be shorter, whereas for species at the MDI1 site being higher light demanding, successional species tends to be taller (Falster and Westoby, 2005). This suggests that additional factors such as topography (slope aspect and elevation), forest types, and environmental conditions (Castro-Luna et al., 2011) might be important in governing plant traits and processes in the region. Therefore, the inclusion of the combined effect of disturbances, topographical heterogeneity, forest types, and environmental conditions might play an important role in the regulation of plant species diversity and functional traits in the study site.

Effect of Varied Disturbance Intensity on Soil Microbial Diversity
The study reported that soil bacterial communities were mostly composed of Proteobacteria, Acidobacteria, Actinobacteria, Planctomycetes, Verrucomicrobia, and Bacteroidetes across the study sites. Various previous studies confirm the present findings conducted in temperate forest ecosystems (Shen et al., 2013), subtropical forest ecosystems (Delgado-Baquerizo et al., 2018), and temperate grassland ecosystems . It was found that plant species diversity and bacterial species diversity showed a similar pattern of increment with increment in disturbance intensity. Thus, soil bacterial diversity was found to be positively correlated with plant species diversity. This might be due to an increase in aboveground plant communities and belowground microbial communities. Plant communities are also known to modify the resource availability and microclimatic conditions for belowground soil level which affects microbial growth responses (Zumsteg et al., 2012). The study was also in confirmation with other studies conducted in temperate grassland ecosystems (Chai et al., 2016) and tropical forest ecosystems (Barberán et al., 2015). Proteobacteria remained to be the most abundant bacterial phylum through the study site. These are represented by fast-growing copiotrophs that are acclimatized to increased carbon content and nutrient availability (Cong et al., 2015). This was in support of this study that showed increased density and basal area at LDI sites leading to higher carbon content and nutrient accumulation than other sites. This supports the suggested hypothesis that soil bacterial diversity would be largely determined by soil attributes supporting a strong correlation between soil microbial biomass and bacterial phyla. Soil microbial biomass carbon and nitrogen remained to be the only variables that varied significantly among study sites. The values were found to be maximum for the LDI site followed by the MDI2 and MDI1 sites. Values for soil MBC and soil MBN showed a similar pattern of increment with bacterial phyla such as Proteobacteria, Verrucomicrobia, and Planctomycetes, suggesting a strong link between soil microbial biomass and community composition. Thus, the present study confirms that the role of moderate disturbances is in one way or another helpful in the regulation of aboveground and belowground ecosystem processes for sustainable growth and diversity in the region. The increased basal area and tree density at LDI sites might be the most important reason for increased MBC and MBN. This could be attributed to an increase in litter inputs which governs soil organic matter and nutrients (Tiwari et al., 2019). The values for MBC and MBN decreased with increasing depth throughout the study site.

Relationship Between Plant Community Characteristics and Soil Microbial Diversity
Many studies have been reported confirming the important link between aboveground community characteristics and belowground microbial composition at regional and landscape levels (de Vries et al., 2012;Legay et al., 2014;Sterkenburg et al., 2015). Plant community characteristics and soil microbial biomass were important predictors of soil bacterial community composition and diversity in the sampling site. The Pearson correlation coefficient diagram ( Figure 5) suggested a strong positive correlation between plant community composition and soil MBC and MBN. This suggests an increased amount of litterfall due to increased basal area at LDI sites and optimum decomposition rates might lead to increased soil MBC and MBN. This was found to be consistent with other studies that suggested that soil MBC and MBN are highly influenced by carbon content in pine forests (Arunachalam et al., 1996). TBA and density were also found to have a negative correlation with bacterial phyla Actinobacteria and Acidobacteria. Actinobacteria maximally dominated in dryland regions thriving under scarcenutrient conditions and relatively drier environments (Maestre et al., 2015). Hence, MDI1 and MDI2 sites being lower in tree basal area and density have higher counts of Actinobacteria. Similarly, CC was found to have a negative correlation with bacterial phylum Actinobacteria suggesting reduced canopy cover leads to increased solar radiation intensity and, thus, suitable conditions for phylum growth. The effect of plant traits in governing belowground microbial communities was not found in this study. This was in confirmation with other studies that found no link between plant traits and microbial communities (Barberán et al., 2015;Leff et al., 2018). A strong positive link between SMN and bacterial phyla suggests that soil nitrogen plays an important role in governing microbial composition since nitrogen is considered to be the most important limiting nutrient in the forest ecosystem (Chai et al., 2016). The first PCA axis was associated with SLA and SM and opposed by CC, whereas the second axis was associated with MBC and MBN. This suggests that species with higher SLA tend to have lower CC due to the effect of disturbance. PCA revealed that the three sites having variation in disturbance intensity (low and moderate) differed in trait values; however, a partial overlap between the MDI1 and MDI2 sites was still observed. This shows that the trait values for the LDI site are different from the MDI1 and MDI2 sites FIGURE 6 | CCA plot diagram to analyze the major gradients among the combination of explanatory variables. The length of vectors in the ordination diagram represents the strength of the correlation, while the angle between two vectors shows the degree of correlation among variables with each axis.
( Figure 6). The analyses confirmed that the plant traits (CC and SLA) along with MBN explained around one-third of the variation in disturbance intensities among species in all the sites. CCA explained that MBN, CC, and SLA were found to be major contributing factors in the LDI and MDI2 sites, respectively. This suggests overdependence of species present in the MDI2 site. For instance, Quercus leucotrichophora, one of the most dominant species in the region and utilized for local livelihood, is known to be heavily lopped for fuel and fodder (Kumar and Ram, 2005). Since the species remains to be restricted in certain environmental conditions, the same species could not be found in another environmental condition to account for the effect of disturbance within species.

Implications and Conclusions
Varying levels of disturbance intensity are known to affect ecosystem functioning and productivity in forest ecosystems. Positive effects of disturbance could be accounted as increased resource availability for recruitment of new species; thus, increased diversity and SR in MDI2 sites and negative effects could be considered through variation in plant functional traits throughout the sites. The study concludes that moderate disturbance is important for promoting species diversity and SR by recruiting new species to the community leading to mixing of pioneer and climax-community species in the region. Plant traits related to competitive ability (height) and resource exploitation (SLA) were suggested to be good predictors of disturbance effect in the study. Our study suggests that species with more conservative ecological strategies would be more prone to continued disturbance intensity. The Himalayan forests being highly fragile should be protected following periodic forest management with initiation at the site level, to maintain sustainable ecosystem services and biodiversity in the region.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI (accession: PRJNA691764, PRJNA692097, and PRJNA692098).

AUTHOR CONTRIBUTIONS
The author confirms being the sole contributor of this work and has approved it for publication.

ACKNOWLEDGMENTS
The author is grateful to Late Shri Ambadutt Sharma (Forest Guard) for his continuous support during the sampling period. The author is also thankful to Abhishek K. Verma (Forest Ecology and Climate Change Division, FRI, Dehradun) for his help in statistical analysis.