## 3 Results
### Overall Relationships
Based on the results presented on the heatmap (Figure 5A), it was observed that the control and TM were similar in terms of characteristics and contained the lowest overall physiological feature loads. The next group consisted of VC and ZM-d, which were characterized by the higher activities of their traits. The highest potential activity of physiological features was demonstrated by the rhizobiome treated with ZM-w, which was most distinct from the other variants. The control showed activities of sulphur respiration and moderate activity of magnesium, iron, nitrogen, nitrate, aerobic ammonia oxidation, and nitrate reduction. The TM variant showed activities of thiosulphate, sulphite, and arsenate respiration, as well as activity typical of gut microbiota and mammalian and invertebrate pathogens. In the ZM-d variant, activities were observed in traits such as gas-producing bacteria, arsenic detoxification, fumarate respiration, chitinolysis, xylanolysis, aerobic ammonia oxidation, and nitrate reduction. Additionally, ureolytic activity, moderate nitrogen binding, phytopathogenicity, a set of organic compound degradation traits, a set of traits associated with sulfur metabolism in the dark, denitrification capability, and a set of traits associated with respiration and nitrogen compound reduction and chemotherapotropy were observed. The VC variant was partly similar to ZM-d, but did not show such clear metabolic activity in the dark, instead being characterized by thiosulphate and arsenate respiration activity, as well as invertebrate pathogenicity. In contrast to ZM-d, it showed higher activity of parasitic bacteria, dark sulphite and iron oxidation, and high cellulolytic and aromatic compound degradation activities. The ZM-w variant showed high activities in 42 out of 60 analyzed physiological traits. Traits with low potential load included thiosulphate, arsenate, manganese, iron, and sulphur-containing compound respiration, as well as invertible parasites, nitrate, and nitrite ammonification. Based on the results of rigorous network analysis (Figure 5B), it was observed that individual OTUs of bacteria were directly correlated with a small number of metabolic features. Among the few traits correlated with a specific OTU were anoxygenic photoautotrophy S and Fe oxidation, methanomorph, and intercellular parasites correlated with Rhizobiaceae and Vicinamibacteraceae. Most OTUs were directly associated with typical gut bacteria and correlated with sulphate and sulphite respiration. Two exceptions to the relationship observed in other rhizobiome OTUs were Rhiziales and _Bacillus_, whose changes in abundance were not related to rhizobiome abundance. Representatives of the Rhizoides order were a completely different variable, unrelated to any observed variable. The _Bacillus_ was correlated with the detoxification, oxidation, and respiration of arsenic, the respiration of thiosulphate and sulphur compounds, aerobic ammonification, as well as nitrite and nitrate ammonification. Figure 5: Heatmap (**A**) and network analysis (**B**) (co-occurrence analysis), representing relationships between bacterial microbiomes and their metabolic features (r > 0.8).