Bacterial structure and dynamics in mango (Mangifera indica) orchards after long term organic and conventional treatments under subtropical ecosystem

Bacterial strains isolation and identification

Fifty six bacterial cultures were isolated from both management systems (G1 and G2) of mango orchards (rhizosphere) at CISH, Lucknow, India. Isolation of microorganisms using spread plate methods revealed that the Nutrient agar medium had the highest number of colony appearances compared to the Rose Bengal Agar medium. Microbial enumeration showed organic system enriched with higher bacterial and fungal population than conventional system (Fig. 1). From organic system, thirty seven bacteria were isolated out of which, twenty-three isolates were (G⁺), and fourteen were (G⁻). While, in the conventional system, nineteen bacteria were isolated, out of which fifteen were (G⁺) and four were (G⁻) isolates.

Plant growth promotion properties

For plant growth promotory properties out of fifty-six bacterial isolates total, ten bacterial cultures (2, 3, 4, 8, 15, 23 and 31) from the organic system showed positive results for phosphate solubilization. In contrast, three bacterial cultures (I1, I8 and I9) from the inorganic system (conventional system) showed positive phosphate solubilization in Pikovaskya’s agar medium. For siderophore production, bacterial cultures (2, 3, 4, 8, 12 and 26) from the organic system showed positive results, while four bacterial cultures (I1, I6, I8 and I9) inorganic system showed positive results. Bacterial cultures (2, 3, 4 and 8) from the organic system showed positive results for K-solubilization, while five bacterial cultures (I1, I2, I7, I8 and I9) from the inorganic system showed positive K-solubilization. A total of ten isolates (7 from organic and 3 from the inorganic system) possessed Zn-solubilizing activity. The test isolated from the organic system showed better Zn (ZnO), Zn₃ (PO₄)_2, and (ZnCO₃) solubilization as compared to test culture isolated from the inorganic system (Supplementary S1.8).

Acetylene reduction assay (ARA)

Results from acetylene reduction assay showed in aerophilic condition, bacterial isolates 1, 3, 4 (from organic treated soil) and I1, I8 and I9 (conventional system) showed 134.8, 37.70, 36.73, 13.15, 16.70 and 12.87 ppm of ethylene tube⁻¹ h⁻¹, respectively. In case of microaerophilic condition, bacterial isolates 4, 9, I9 showed 24.17, 19.14, and 12.71 ppm ethylene, respectively. Results indicate possible use of these bacterial isolates as a bioinoculant agent for horticultural crops, especially mango and other subtropical climate fruit crops.

Soil enzymatic study

The soil enzymatic activity in the organic system (G1) showed better dehydrogenase activity than the conventional system (G2). For both methods, alkaline phosphatase almost showed similar activity (at pH 11), while in the case of acid phosphatase showed better activity in the inorganic system (G2) as compared to the organic system (G1) at pH level 6.5 (Fig. 2). The dehydrogenase enzyme oxidizes the organic matter, and it belongs to the oxidoreductase type of enzyme. In the process of respiration of soil microorganisms, the dehydrogenase enzyme facilitates the transfer of protons and electrons from the substrate to the acceptor. It was significant to observe that the dehydrogenase activity was higher in organic treated soils (0.784 µg TPF g⁻¹ h⁻¹) than in conventional system (0.053 µg TPF g⁻¹ h⁻¹).

Alpha biodiversity with samples and rarefaction curves

In this segment, by measuring Shannon, Chao1, and observed species metrics, we analyze the microbial diversity within the samples. The chao1 metric measures the richness of the ecosystem, while the Shannon metric is the formula for calculating reported OTU abundances and accounts for both prosperity and equality. The rarefaction curve is provided in Fig. 3 for each metric. Using QIIME software, the metric measurement was done. The impact of both treatments on the microbial complexity and abundance in the sample was also revealed using the Shannon diversity Index (depicting richness and evenness) and Chao 1 representing only richness. Shannon’s diversity index of the bacterial community in the treatment (G1 and G2) was 8.06 and 8.12. The Simpson index in ecology is used to quantify biological diversity in a region, which was also nearly similar in both the treatments. Chao 1 richness estimator showed an increase in species richness. Rarefaction analysis conducted to confirm species richness revealed a difference in the number of reads and OTUs between the samples. The Rare fraction curve had a similar pattern for both samples and showed an impact on the bacterial population in the experiment (Fig. 3a–c).

Bacterial diversity analysis at phyla level

Taxonomic study of the 16S rRNA gene amplicon reads yielded seven classifiable bacterial phyla. Six phyla, namely Acidobacteria, Actinobacteria, Bacteroides, Proteobacteria, Firmicutes, and Chloroflexi were dominant in both the systems. The Organically treated soil (G1) sample harbored a higher percentage of Bacteroidetes (14.55%), Actinobacteria (7.45%), and Proteobacteria (10.82%) as compared to conventional treatment (G2) 8.98%, 5.71%, and 6.64%, respectively. However, phylum Acidobacteria(13.6%), Firmicutes(4.84%), and Chloroflexi (2.56) were higher abundance in conventional treatment as compared to the organic treatment, which showed the same phyla with lesser quantity, i.e., 5.63%, 0.91%, and 0.79% respectively (Fig. 4a).

Distribution of bacterial community at order level

The bacterial orders in both systems were diversified. The most abundant orders in organic and conventional systems were Chitinophagales (Organic-11.32%, Conventional-43%), Elev-16S-573 (Organic-3.09%, Conventional-8.69%), Pedosphaerales (Organic-1.56%, Conventional-3.55%), Opitutales (Organic-2.46%, Conventional-0.27%), Chthoniobacterales (Organic-1.35%, Conventional-2.84%), Bacillales (Organic-0.91%, Conventional-4.84%) and Solibacterales (Organic-1.39%, Conventional-2.26%) (Fig. 4b).

Bacterial community distribution at family level

Bacterial family members were identified and enriched including Pedosphaeraceae (O-1.56%, C-3.55%), Opitutaceae (O-2.46%, C-0.27%), Chthoniobacteraceae (O-1.03%, C-2.68%), Steroidobacteraceae (O-2.05%, C-0.73%), Bacillaceae (O-0.77%, C-4.55%), Chitinophagaceae (O-10.99%, C-5.06%), and Xanthomonadaceae (O-1.39%, C-0.06%) and other families (Fig. 5a).

Bacterial community distribution at the genus level

Comparative abundance of unidentified genus in organic system were uncultured soil bacterium, Glycomyces, Chitinophaga, Lysobacter, Udaeobacter, Bacillus (not detected, 1.85%, 4.77%, 1.19%,1.03% and 0.75% respectively) whereas same genus-group were observed in conventional system with different percentage i.e., 0.11%, not detected, 0.56%, 0.04%, 2.67%, 4.54% respectively (Fig. 5b).

Bacterial communities at species level

Because most of the species were unidentified and uncultured bacterium based on relative abundance, they could not be assigned a species name in either sample. Few species are identified in both systems, like Sphingomonas sp. (O-1.57%, C-1.05%), Bacillus drentensis (O-0.25%, C-2.65%), and Chitinophaga sp. (O-4.64%, C-0.11%) (Fig. 6).

Heat map and PCA analysis

Under long-term exposure of organic and conventional treatments, a microbial shift was observed in the rhizosphere microbiome of mango orchards. Based on percent abundance, nine different microbial genera Acidobacteria, Actinobacteria, Bacteroidetes and Proteobacteria formed Cluster I. While, Firmicutes, Chloroflexi and Opitutales were abundances in cluster II. Cluster III includes Chitinobacterales, Bacillales, Chitinophagarales and Otherales genera. Whereas cluster IV (Elev7-16S-573, Otherales, Solibacterales and Pedobacteriaceae), cluster V (Opitutaceae, Chitnobacteraceae, Bacillaceae, Chitinophagaceae and Otherales), cluster VI (Xanthomonadaceae, Uncultured soil bacterium, Candidatus-Udaeobacter, Lysobacter and Bacillus), cluster VII (Chitinophaga, Glycomyces and Other), cluster VIII (Uncultured bacterium and Others) and cluster IX (Bacillus drentensis and Others) (Fig. 7). The cluster I observed with the highest abundance was closely related to clusters II and III. Cluster IV to IX created large groups and is distantly related to cluster I to III of the microbial groups in organic and conventional systems (Fig. 7). In the organic system (G1), microbial groups like Proteobacteria, Actinobacteria, Bacteroidetes, and Opitutaceae were largely dominated and provided benefits to the mango rhizosphere in terms of nutrient availability, plant growth promotion, and protection against biotic and abiotic stress. Phylum Proteobacteria and Actinobacteria are closely linked with the rhizosphere and identified as potential PGPR. Acidobacteria and firmicutes, on the other hand, were dominated primarily by conventional systems and serve as a bio-indicator of anthropogenic stress caused by excessive chemical fertilizer application. Undefined Acidobacteria is oligotrophic in nature and considered as an indicator of low organic carbon and acidic environment. To desire higher productivity, the indiscriminate use of chemical fertilizers or pesticides in conventional systems leads to low nutrient availability, microbial shift, less PGPR, and developing the environment for Acidobacteria, Firmicutes and Chloroflexia group of microorganisms. Principal component analysis (PCA) was performed for both systems (organic-component 1; conventional-component 2). The total variables of principal component analysis were the percentage of different parameters such as alkaline phosphatase, acid phosphatase, DHA, Acetylene reduction assay (ARA1, ARA2, ARA3), and CFU mL⁻¹ (bacteria and fungi). The results of PCA yielded two components that explained 100% of the total variance in the data and had an Eigen value of 6.1 for component 1. In contrast, 1.8 for component 2 and together they described 100% of the total variance in the data (Fig. 8). In the organic system, the loading factor with score plot indicates that component-1 is positively associated with DHA, ARA1, ARA2, alkaline phosphatase, acid phosphatase while negatively correlated with CFU ARA3 activity. Component-1 explains the 76.42% variance of the experimental data, while component-2 explains 23.58%. The second component (PC2) represents the positive association with DHA, ARA1, ARA2, ARA3 activity, and CFU while negatively correlated with alkaline phosphatase and acid phosphatase. In the conventional system, the loading factor with score plot indicates that component-1 is positively associated with single variable acid phosphatise while negatively correlated with DHA, ARA1, ARA2, ARA3, CFU, and alkaline phosphatase activity. The second component (PC2) of the conventional system showed positive association with DHA, ARA1, ARA2, ARA3 activity, and CFU, while the negative association with alkaline phosphatase and acid phosphatase.

Source: Ecology - nature.com