Statistical analysis of plant phenotypic data
The relative changes in shoot biomass, root biomass in the combinations of the three plants and three rhizobacterial strains was analyzed with R Studio (Version 3.5.2). First, the normality and homogeneity of variance of the data was assessed and when the two assumptions were not met, the data was subjected to Box-Cox transformation using the package MASS. Differences were tested by two-way analysis of variance (ANOVA). A Tukey-HSD test was used to separate group mean values when the ANOVA was significant at p  < 0.05. The ANOVA table is shown in Supplementary Material, Table S1 . Differences between rhizobacterial treatments and non-treated control on raw phenotype parameters were compared by Student’s t -Test.