3.1. Structural shifts in bacterial microbiota
Both Shannon and Simpson indexes indicated significant differences in bacterial alpha diversity between all groups, with I roots showing higher values compared to both seeds and NI roots (Fig. 1A and B). Furthermore, PCoA based on Bray–Curtis and Jaccard dissimilarity followed by PERMANOVA test highlighted a significant difference between the composition of the communities (p = 0.01, Supplemetaty Table S1a, Fig. 1C). Distribution of taxa based on their patterns of appearance across treatments (Fig. 1D) revealed only one core genus shared between seeds, NI roots and I roots, which showed an overall relative abundance of 75% across treatments. Four and seven genera were exclusive in NI and I roots, respectively. In line with these observations, relative abundance profiles at the genus level revealed clear shifts in bacterial microbiota between seeds, NI roots and I roots (Fig. 1E). Seeds were dominated by the genus Pantoea which accounted for 70% of the bacterial population followed, in decreasing order, by the genus Anaerobacillus, Staphylococcus and Cutibacterium. In the transition from seeds to seedlings, Pantoea abundance increased up to 95% in NI roots, with the remaining 5% including, in decreasing order, members from the genus Curtobacterium, Paenibacillus, Luteibacter, Massilia and Staphylococcus. Compared to NI roots, an increase in the relative abundance of Massilia was observed in I ones, accounting for almost 10% of the bacterial microbiota. Furthermore, I roots showed a 40% reduction in the relative abundance of Pantoea compared to NI ones, concomitant with the emergence of 5 genera from the Pseudomonadota phylum, including, Brevundimonas, Brucella, Luteibacter, and Advenella. Importantly, the inoculant H. seropedicae accounted for nearly 13% of the bacterial microbiota, suggesting the ability of this inoculant to penetrate the roots and multiply within wheat root tissues. Finally, two additional genera from the Actinomycetota phylum, Brevibacterium sp. and Microbacterium sp., were also exclusively detected in I roots.
3.2. Structural shifts in fungal microbiota
Regarding the fungal microbiota, I roots showed significantly higher alpha diversity indexes compared with both seeds and NI roots, which did not differ significantly (Fig. 2A and B). Furthermore, PCoA based on Bray–Curtis and Jaccard dissimilarity showed a significant difference in the composition of the fungal microbiota (p = 0.01, Table S1, Supplementary Table S1b, Fig. 2C). Three genera were shared between seeds, NI roots and I roots, which showed an overall relative abundance of 50% across treatments. In the transition from seeds to seedlings, NI and I roots hosted three common fungal genera, while exclusive fungal genera were found only in I roots (Fig. 2D). Relative abundance profiles at the genus level (Fig. 2E) showed that seeds were dominated by the genus Alternaria which accounted for 82% of the seed microbiota, followed, in decreasing order of abundance, by Fusarium sp., Epicoccum sp. and Meyerozyma sp. Minor taxa (relative abundance < 1%) were also identified in the seed, including Stemphylium, Wallemia, Pleurotus, and Malassezia. In comparison with seeds, both NI and I roots showed a decrease in the population of Alternaria of 27 and 46%, respectively. In NI roots, this decrease in Alternaria populations was accompanied by the emergence of the unicellular yeast Rhodotorula, which accounted for almost 50% of the fungal microbiota, followed members from the genera Meyerozyma, Malassezia and Wallemia. Compared to NI roots, Rhodotorula population decreased drastically and accounted for less than 5% of the population in I roots. In turn, I roots exclusively hosted the yeast from the genus Filobasidium, which accounted for over 40% of the total fungal community.
3.3. Impact of H. seropedicae on biomarker taxa and microbiota functionality
The relative abundance profiles generated with the number of reads were used to detect genera with potential biological relevance related to H. seropedicae inoculation. The LDA scores resulting from LEfSe analysis highlighted that 12 out of 28 genera significantly explained the greatest differences between communities (LDA > 2.0 and p-value < 0.05; Fig. 3A). Most of these genera belonged to the Pseudomonadota and Ascomycota phyla, followed by Actinomycetota. While most of the biomarker taxa included those exclusively found either in NI or in I roots, LDA scores also indicated that the relative abundance of Rhodotorula was significantly higher in NI roots, while Massilia was significantly higher in I roots.
Furthermore, the functional profiles of each microbiota were predicted by comparing the sequencing data with FAPROTAX database. A total of 8 functional sub-categories within four major categories, “Energy source”, “N-cycle”, “Trophic mode” and “Guild”, were identified and linked to the microbial communities across seeds, NI roots and I roots (Fig. 3B). Compared to seeds and NI roots, I roots harbored higher relative abundance of bacterial OTUs involved in N-cycle functions, specifically in NO3 respiration/reduction and N fixation, and in Energy source functions, specifically in aerobic chemoheterotrophy. Regarding the fungal microbiota, higher relative abundance of potential pathotrophs was found in the seeds, while their abundance in NI and I roots was generally lower and more variable between replicates. The Guild category was exclusively affected in seeds and I roots, where soil and litter saprotrophs were exclusively linked with the seed or with I roots, respectively, while the abundance of fungal parasites was similar in both environments.