Soil bacterial communities play crucial roles in ecosystem functions and biogeochemical cycles of fundamental elements and are sensitive to environmental changes. However, the response of soil bacterial communities to chronosequence in tropical ecosystems is still poorly understood. This study characterized the structures and co-occurrence patterns of soil bacterial communities in rhizosphere and bulk soils along a chronosequence of teak plantations and adjacent native grassland as control. Stand ages significantly shifted the structure of soil bacterial communities but had no significant impact on bacterial community diversity. Bacterial community diversity in bulk soils was significantly higher than that in rhizosphere soils. The number of nodes and edges in the bacterial co-occurrence network first increased and then decreased with the chronosequence. The number of strongly positive correlations per network was much higher than negative correlations. Available potassium, total potassium, and available phosphorus were significant factors influencing the structure of the bacterial community in bulk soils. In contrast, urease, total potassium, pH, and total phosphorus were significant factors affecting the structure of the bacterial community in the rhizosphere soils. These results indicate that available nutrients in the soil are the main drivers regulating soil bacterial community variation along a teak plantation chronosequence.

Arbuscular mycorrhizal (AM) fungi play a crucial role in promoting plant growth, enhancing plant stress resistance, and sustaining a healthy ecosystem. However, little is known about the mycorrhizal status of teak plantations. Here, we evaluated how the AM fungal communities of rhizosphere soils and roots respond to different stand ages of teak: 22, 35, 45, and 55-year-old from the adjacent native grassland (CK). A high-throughput sequencing method was used to compare the differences in soil and root AM fungal community structures. In combination with soil parameters, mechanisms driving the AM fungal community were revealed by redundancy analysis and the Mantel test. Additionally, spore density and colonization rates were analyzed. With increasing stand age, the AM fungal colonization rates and spore density increased linearly. Catalase activity and ammonium nitrogen content also increased, and soil organic carbon, total phosphorous, acid phosphatase activity, available potassium, and available phosphorus first increased and then decreased. Stand age significantly changed the structure of the AM fungal community but had no significant impact on the diversity of the AM fungal community. However, the diversity of the AM fungal community in soils was statistically higher than that in the roots. In total, nine and seven AM fungal genera were detected in the soil and root samples, respectively. The majority of sequences in soils and roots belonged to Glomus. Age-induced changes in soil properties could largely explain the alterations in the structure of the AM fungal community along a chronosequence, which included total potassium, carbon-nitrogen ratio, ammonium nitrogen, catalase, and acid phosphatase levels in soils and catalase, acid phosphatase, pH, and total potassium levels in roots. Soil nutrient availability and enzyme activity were the main driving factors regulating the shift in the AM fungal community structure along a chronosequence of the teak plantations.