The Rhizosphere Selects for Particular Groups of Acidobacteria and Verrucomicrobia

There is a lack in our current understanding on the putative interactions of species of the phyla of Acidobacteria and Verrucomicrobia with plants. Moreover, progress in this area is seriously hampered by the recalcitrance of members of these phyla to grow as pure cultures. The purpose of this study was to investigate whether particular members of Acidobacteria and Verrucomicrobia are avid colonizers of the rhizosphere. Based on previous work, rhizosphere competence was demonstrated for the Verrucomicrobia subdivision 1 groups of Luteolibacter and Candidatus genus Rhizospheria and it was hypothesized that the rhizosphere is a common habitat for Acidobacteria subdivision 8 (class Holophagae). We assessed the population densities of Bacteria, Verrucomicrobia subdivision 1 groups Luteolibacter and Candidatus genus Rhizospheria and Acidobacteria subdivisions 1, 3, 4, 6 and Holophagae in bulk soil and in the rhizospheres of grass, potato and leek in the same field at different points in time using real-time quantitative PCR. Primers of all seven verrucomicrobial, acidobacterial and holophagal PCR systems were based on 16S rRNA gene sequences of cultivable representatives of the different groups. Luteolibacter, Candidatus genus Rhizospheria, subdivision 6 acidobacteria and Holophaga showed preferences for one or more rhizospheres. In particular, the Holophaga 16S rRNA gene number were more abundant in the leek rhizosphere than in bulk soil and the rhizospheres of grass and potato. Attraction to, and colonization of, leek roots by Holophagae strain CHC25 was further shown in an experimental microcosm set-up. In the light of this remarkable capacity, we propose to coin strain CHC25 Candidatus Porrumbacterium oxyphilus (class Holophagae, Phylum Acidobacteria), the first cultured representative with rhizosphere competence.


Introduction
The phyla Acidobacteria and Verrucomicrobia are among the most dominant bacterial groups present in most soils [1][2][3][4][5][6]. The phylum Acidobacteria consists of at least 26 monophyletic groups, so called subdivisions, whereas the Verrucomicrobia have at least seven subdivisions [5,7]. Both phyla are intriguing prokaryotes given their presumed roles in soil ecosystems and also because the vast majority of species belonging to these phyla remains uncultured to date [4,8]. In particular, very little is known about the ecological roles of members of the Acidobacteria and Verrucomicrobia in plant-soil ecosystems.
Most of the available data on the ecology of Acidobacteria and Verrucomicrobia in plant-soil ecosystems comes from studies in which cultivation-independent (metagenomic) approaches were applied. For instance, members of the Verrucomicrobia have been shown to be present in varying plant-soil ecosystems [9][10][11][12][13][14]. Also, representatives of the phylum Acidobacteria were found in these systems, although they tended to be more associated with bulk, than with rhizosphere soils [13,15]. In soil, pH seemed to play an important role as a determinant of acidobacterial assemblages [5,[16][17][18][19][20]. In addition, mineral composition [21,22], temperature [18] and nutrient availability [15,[23][24][25] were important. The relative abundances of Acidobacteria in clone libraries from pasture soil were found to be between 7 -14% [1] and of Verrucomicrobia around 1.9 % [26]. In Brazilian Atlantic forest soil, such relative numbers were higher, i.e. 29 -54% for Acidobacteria and 0.6 -14% for Verrucomicrobia [27]. Using barcoded pyrosequencing on soil-extracted DNA, Verrucomicrobia operational taxonomic units (OTUs) comprised 35% of all bacterial OTUs retrieved from 112 undisturbed mineral soils sampled across different continents (North and South America, Europe and Antarctica) [28]. In that study, Spartobacteria and subdivision 3 were found to be the most abundant subdivisions among Verrucomicrobia. In 88 different soils, Acidobacteria contributed to 35% of total bacterial OTUs, whereas for Verrucomicrobia it was 0.9% [5]. In that study, Acidobacteria subdivision 1 accounted for 7.4% of all bacterial and for 17.6% of all acidobacterial OTUs. Considering all data, there appears to be wide variation in the contribution of both phyla and their separate subdivisions to the bacterial communities present in soil environments.
A major task is the cultivation of some of the diversity of the Acidobacteria and Verrucomicrobia detected by direct molecular approaches in soils. Thus, based on cultured representatives from soil, subdivisions 1, 2, 3, 4 and 6 of Acidobacteria have been found in pasture soil [4,16,[29][30][31], and subdivision 8 (class of Holophagae) in the leek rhizosphere [32]. Verrucomicrobia subdivision 2 (Spartobacteria) and subdivision 3 strains were isolated from pasture soils [4,33], whereas Verrucomicrobia subdivision 1 strains were found in the rhizospheres of potato and leek [12,32,34]. There appears to be a distinction in verrucomicrobial and acidobacterial compositions between bulk and rhizosphere soils, as Verrucomicrobia subdivision 1 and Holophagae isolates were found in close proximity to plant roots. Rhizosphere competence was demonstrated for two distinct Verrucomicrobia subdivision 1 representatives, provisionally denominated Luteolibacter and Candidatus genus Rhizospheria [12]. Rhizosphere competence has never been demonstrated for representatives of the phylum Acidobacteria. It was, therefore, a challenge to explore the distribution of Holophagae, in relation with other groups of Acidobacteria belonging to subdivisions 1, 3, 4, 6 and of Verrucomicrobia subdivision 1 Luteolibacter and Candidatus genus Rhizospheria, in the rhizosphere of different plants. For that purpose, we used real time quantitative PCR (qPCR), accurately calibrated by making use of culturable representatives of these groups [12,32].
In this study, we addressed the question whether the cultured Holophagae strains are rhizosphere-competent, in analogy to the study performed on two different Verrucomicrobia subdivision 1 groups [12]. We define rhizosphere competence as the ability of bacteria to move towards plant roots and to grow on root-released nutrients [35]. The study includes data from the field, focusing on different Acidobacteria subdivisions and Verrucomicrobia subdivision 1, and from a rhizosphere microcosm study. The latter was set up to further explore the possible competence in the leek rhizosphere of a soil-introduced Holophagae strain.

Field site, soil and plant sampling, sample processing and analysis procedures
The field site was located at the experimental farm 'De Droevendaal' (51°59'32"N, 5°40'12''E), Wageningen, The Netherlands. The soil was a loamy sand containing 2% organic matter, with a water holding capacity of 25% and a pH (KCl) of 4.8. Before onset of the experiment, the field (21 by 25 m) was covered with a permanent grass ley (commercial mix, containing Lolium perenne as the main plant species) and maintained under agricultural management practices. Then the field was divided into 16 plots of 4 by 5 m in size with a distance of 1 m between subplots and the margins of the field. Four treatments, i.e. fallow, grass, potato and leek, were in fourfold applied over the field according to a randomized scheme. Therefore, grass was removed from 12 plots, whereas it was maintained on four plots (grass). Two fallow plots were immediately planted with potato (Solanum tuberosum L. cultivar Agria) or leek (Allium porrum cultivar Kenton, Nunhems Seeds BV, The Netherlands) and the other four plots were kept fallow (non-rooted bulk soil).
Seed potatoes and leek nursery plants were planted in May 2009. Organic agricultural management practices were continued for all plots, which specifically comprehends no use of pesticides or chemical fertilizers and weekly removal of weed plants by hand. Samples from each subplot were taken in June, July and September. Samples from the potato and leek plots were taken as individual plants, whereas those from the grass and fallow soil (one per plot) were taken with a soil bore (diameter size of 7 cm) to a maximal depth of 15 cm in fallow soil. All samples were directly processed in the laboratory, where soil adhering to the grass, leek and potato roots after manual shaking of the plants was considered as rhizosphere soils. For (non-rooted) bulk soil, samples from the 5-10 cm horizons of the fallow plots were singled out. Soil pH was measured in all bulk and rhizosphere soils in 0.01 M CaCl 2 (1: 10 w/v ratio) according to the procedure described in [36].

Leek Rhizospere colonization by soil-indigenous Holophagae cells and strain CHC25 in a plant-soil microcosm
The behavior of Holophagae species and their representative strain, CHC25 [32], was studied near leek roots in non-sterile and sterilized Vredepeel soil with or without leek plants, using the same microcosm set up (Kuchenbuch-style) as previously described for Luteolibacter and Candidatus genus Rhizospheria strains in Nunes da Rocha et al. [12]. In short, non-sterile and non-inoculated soil (set up A), or sterilized soil with approximately 10 5 strain CHC25 cells per g dry soil (set up B), or with a 1 cm non-inoculated and sterilized soil layer placed between strain CHC25-inoculated soil and the membrane separating leek roots from soil (set up C), or the same as set up C, but then without leek plants (set up D). After 35 days, rings (in triplicate) were destructively sampled and soils at 0-2 mm and 10-12 mm from the nylon membrane with roots were singled out and homogenized. One-gram subsamples were drawn for later DNA extraction and Holophagae-specific real-time qPCR analysis [32].

DNA extraction from soils and real-time quantitative PCR analyses
DNA from all bulk and rhizosphere soils (Vredepeel and Droevendaal soils) was extracted using the PowerSoil Isolation Kit (MO BIO Laboratories, Inc., CA, USA) following the instructions provided by the manufacturer. Quantitative PCR primers Eub338 [37] and Eub518 [38] were used for quantification of bacteria, representing 'total bacteria' within the domain of Bacteria (Table 1). Primer combinations VS1Af/ VS1Ar, VS1Bf /VS1Br and Acg8f/ Acg8r were, respectively, used for quantification of Verrucomicrobia subdivision 1 groups of Luteolibacter and Candidatus genus Rhizospheria, and of Holophaga (representing the class Holophagae) [12,32]. Four new primer systems, used for quantitative detection of subdivisions 1, 3, 4 and 6 acidobacteria (Table 1), were designed based on almost entire (> 1300 bp) 16S rRNA gene sequences of cultured strains IGEO12 (subdivision 1, accession number GU187028), IGEO15 (subdivision 3, GU187034), IGEO17 (subdivision 4, GU187032) and IGEO01 (subdivision 6, GU187036) (all strains are described in George et al. [29]), according to the procedure described in Nunes et al. [32]. In short, primers were validated in three steps. The first step comprehended in silico validation of primers. Therefore, alignments were made for each acidobacterial subdivision using 16S rRNA gene sequences of these strains and those of related bacterial groups retrieved from the SILVA database, release 102 [39]. Primers, specific for each subdivision, were designed based on conserved sequences and checked for absence on possible occurrences of mispriming events using Primer-BLAST software (http://www.ncbi.nlm.nih.gov/tools/ primer-blast/). The second step comprehended validation by PCR on DNA extracts from pure culture strains. Therefore, specificity of designed primers, per subdivision, was checked with DNA from corresponding (target) and non-corresponding (non-target) strains by standard PCR amplification ( Table 1). As non-target strains, all non-corresponding Acidobacteria and Verrucomicrobia subdivision 1 strains were chosen, supplemented with Agrobacterium tumefaciens UBAPF2 (Alphaproteobacteria), Burkholderia cepacia LMG 1222T (Betaproteobacteria), Escherichia coli E1 (Gammaproteobacteria), Streptomyces griseus IPO 857 (Actinobacteria), Flavobacterium columnar 2003/035 (Bacteroidetes) and Bacillus subtilis Bs4 (Firmicutes). All these strains were derived from the strain collection of Plant Research International (Wageningen, The Netherlands). The third step comprehended specificity checks on amplicon sequences derived by standard PCR with these primers from Droevendaal soil DNA extracts. Therefore, soil extracted DNA was PCR amplified and individual amplicons were cloned into the pGEM-T easy vector (Promega, WI, USA) for sequencing. A total of 192 sequences (48 per primer system) from randomly selected clones were aligned using MEGA 4 software [40] and individually compared by BlastN-assisted database searches (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Finally, standard curves for each of the four primer systems were made by qPCR, based on ranges of between 10 and 10 9 cells per corresponding strain, whereas for the bacterial primer system, cells of Pseudomonas fluorescens Pf5 were used. Calibration curves were made in triplicate by plotting measured threshold cycle (Ct) values against 10 log cell number for each qPCR system. Line slopes and intercepts were calculated by linear regression analysis (Genstat 15th edition, Hemel Hempstaed, UK) and the amplification efficiency (Ae) of the different primer systems was calculated using the formula Ae = 10 (-1/slope) . Theoretical dynamic ranges for all qPCR systems were determined according to Nunes da Rocha et al. [32].
Quantitative PCR systems (Table 1) were applied for molecular quantification of different Acidobacteria and Verrucomicrobia subdivisions in Droevendaal (with all eight primer systems) and Vredepeel (with only the one of Holophagae) soils. Therefore, DNA extracts were 10-fold diluted to approximately 5 ng per 25 μL reaction mixture, prior to running under the conditions previously described in Nunes da Rocha et al. [32]. A total of three qPCRs per primer system were run for each sample and obtained Ct values were averaged prior to conversion to log cell equivalent numbers using the appropriate regression equation for each primer system.

Statistical comparisons and multivariate analyses
Statistical comparisons, based on 10 log-transformed bacteria, subdivisions 1, 3, 4, 6 acidobacteria and Holophaga, Luteolibacter and Candidatus genus Rhizospheria cell equivalent (Ceq) numbers (expressed per g dry soil) were made between (1) different bulk soils over time, (2) between rhizospheres and bulk soils for calculation of Δ Ceq rs, b values, and (3) between fractions of total bacteria (individual population size as fraction of total bacteria within the same sample) in rhizopheres and bulk soil over time. Comparisons between rhizosphere and bulk soils sampled over time were also made on the basis of pH values. All comparisons were based on four replicate samples per treatment (rhizosphere type or bulk soil sampled over three time points).
The effects of grass, potato and leek roots in soil on the eight different populations, expressed as Δ Ceq rs, b , were calculated for each population by subtraction of the log Ceq number (per g dry soil) in bulk soil from each of the corresponding rhizosphere soils. Values were presented as 'positive' when Ceq numbers were significantly higher in rhizosphere than in bulk soil, 'negative' when significantly lower, and 'zero' when statistically indistinguishable.
In experimentation with the plant-soil microcosms, comparisons in log-transformed Holophaga cell equivalent numbers (per g dry soil) were made between: (1) 0-2 and 10-12 mm layers of set up A, (2) 0-2 and 10-12 mm layers of set up B, (3) 0-2 and 10-12 mm layers of set up C, (4) 0-2 mm layers of set ups C and D. Averages per soil layer were based on triplicate values for each of the four microcosm set ups. Significance of differences were calculated with two-way ANOVA (Genstat 15th edition). Least significant differences were calculated from standard errors of difference. All differences were considered to be significant at levels of P < 0.05.
Multivariate analysis (CANOCO for Windows version 4.5, Biometris, Plant Research International, The Netherlands) was performed on all soil samples using sample type (rooted versus non-rooted soils), period of sampling, plant species (all nominal variables) and pH (numerical variable) as the 'environmental' variables and log Ceq numbers for each group (per g dry soil), as 'species' variables. Gradient lengths were calculated by detrended correspondence analysis (DCA) in a first step, and correlations between 'environmental' and 'species' variables in a second step by redundancy analysis (RDA). Monte Carlo permutation test (499 permutations) was included to calculate significance of effects on species variables.

Specificity of Acidobacteria subdivisions 1, 3, 4 and 6 quantitative PCR primer systems
Four, of eight, qPCR primer systems (Table 1) were newly designed for the purpose of this study. From Primer-Blast analysis it was predicted that all primers targeting subdivisions 1, 3, 4 and 6 of Acidobacteria would specifically amplify 16S rRNA gene sequences of the targeted subdivisions. Standard PCR amplifications using the four primer systems on genomic DNA extracts from respective target strains invariably resulted in single amplicons of the expected sizes (Table 1) in the absence of any visible primer diming or other products resulting from primer mismatching (data not shown). Standard PCR amplifications with these four primer systems on genomic DNA extracts from non-target strains from different bacterial phyla (Firmicutes, Proteobacteria, Bacteriodetes) as well as from non-corresponding Acidobacteria and/or Verrucomicrobia subdivision 1 strains resulted in absence of any bands under the applied amplification conditions. Sequence comparisons of 192 amplicons, made with the four primer systems in standard PCRs with Droevendaal soil DNA extract as template, consistently revealed similarities of 96% and over with database sequences belonging to the expected subdivisions, with the exception of the primer set designed for detection of subdivision 6 species that revealed matches with the corresponding subdivision at 46 occasions, whereas at the other two occasions, sequences showed closest matches with subdivision 10 species. In total 22 distinguishable sequence groups (five of subdivision 1, seven of subdivision 3, four of subdivision 4 and six of subdivision 6), containing one to up to 18 identical sequences per group, were deposited in the EMBL Nucleotide Sequence Database and available under accession numbers FN994868 to FN994889.
Quantitative PCR on a density range of 10 -10 9 cells per target strain for each qPCR system resulted in linear regression equations with R 2 values of 0.9841 and higher. Calculated amplification efficiency values (in %) ranged

Plant growth and soil pH in the experimental field plot
Potato and leek plants grew normally in the field plots during the experimental period from May -September 2009 in the absence of any visible harm caused by pests, pathogens or abiotic stressors. The average pH values over the different samples revealed significant effects of plant released protons in the rhizosphere (related to plant growth), but no effect of time. The average pH values were significantly (n=4, P≤ 0.05) lower in bulk soil (4.77 ± 0.09) and in the rhizospheres of potato (4.78 ± 0.33) and leek (4.65 ± 0.19) than in that of grass (5.07 ± 0.16).

Dynamics of Bacteria, Acidobacteria subdivision 1, 3, 4, 6 and Holophagae, Luteolibacter and Candidatus genus Rhizopheria groups in bulk soil
The qPCR analyses (data expressed in log cell equivalents, Ceq, per g dry soil), revealed that total bacterial numbers significantly (n=4, P≤ 0.05) declined over time from 9.70 (May) to 8.61 (September) (Figure 1). Luteolibacter and Candidatus genus Rhizospheria numbers also declined over time, respectively, from 5.67 and 5.92 in May to 3.16 and 5.02 in September. The dynamics of the acidobacterial subdivisions 1, 3, 4, 6 and Holophaga was diverse. Subdivisions 3, 4 and 6 declined over time, respectively, from 8.37, 8.62 and 9.55 in May to 7.80, 6.74 and 7.78 in September. This in contrast to subdivision 1 acidobacteria, of which the numbers roughly remained the same (between 8.32 in May to 8.57 in September) and to Holophaga, whose numbers after an initial decrease, significantly increased, from 5.17 in May to 5.42 in September. Removal of the grass layer from the fallow plot thus led to a unique increase in estimated Holophagae cell numbers.

Effect of different Rhizospheres on Bacteria, Acidobacteria subdivision 1, 3, 4, 6, Holophagae, Luteolibacter and Candidatus genus Rhizospheria estimated cell numbers
Positive Δ Ceq rs, b values (log cell equivalent numbers from bulk soil subtracted from those from corresponding rhizospheres) were found for bacteria across all three plant species and sampling periods, indicating that plant roots stimulated bacterial growth in soil ( Figure 2). For Luteolibacter, positive Δ Ceq rs, b values were also found at all occasions, i.e. in the rhizospheres of all three plant species in June, July and September. For Candidatus genus Rhizospheria, Δ Ceq rs, b values were positive in all rhizospheres taken in July and September and in the leek rhizosphere in June. However, these were negative in the June samples from the grass and potato rhizospheres. Remarkably, Candidatus genus Rhizospheria was specifically enhanced in the rhizosphere of leek as compared to bulk soil and the other two rhizospheres.
The different acidobacterial subdivisions did not always prefer rhizosphere over bulk soils. In concrete terms, the Δ Ceq rs, b values of subdivision 1 acidobacteria were negative in all grass and potato rhizospheres, zero in the leek rhizosphere in June and July, indicating no effect of plant roots on subdivision 1 acidobacteria, and slightly positive in the leek rhizosphere in September (about two-fold higher in the leek rhizosphere than in bulk soil). Representatives of this subdivision thus grossly remained unaffected in the leek rhizosphere where they were  stimulated in their growth later in the season. For subdivision 3 acidobacteria, Δ Ceq rs, b values were negative in all rhizospheres in June and in the grass and potato rhizospheres in July, were zero in the leek rhizosphere in July and in the grass and potato rhizospheres in September. Again these were positive in the leek rhizosphere in September (about four-fold higher than in bulk soil). For subdivision 4 acidobacteria, the Δ Ceq rs, b values were negative in the leek rhizosphere in June, positive in the grass and leek rhizospheres in September (respectively two and eight-fold higher than in bulk soil) and close to zero in all other samples. Representatives of this subdivision thus remained grossly unaffected in the potato rhizosphere. For subdivision 6 acidobacteria, the Δ Ceq rs, b values were close to zero in June and July and positive in all three rhizospheres in September (between five and 25-fold higher than in bulk soil). Moreover, members of this subdivision had a stronger preference for grass and potato rhizospheres than for the one of leek. For Holophaga, the Δ Ceq rs, b values were positive in all leek rhizospheres across time (between four and 16-fold higher than in bulk soil), whereas they were negative or zero in the ones of grass and potato. The behavior of Holophagae in the three rhizospheres was thus different from that of all other subdivisions of the phylum Acidobacteria, in the sense that this group showed a strong preference for the leek rhizosphere throughout the experimental time period.

Contribution of Bacteria, Acidobacteria subdivision 1, 3, 4, 6, Holophagae, Luteolibacter and Candidatus genus Rhizospheria to total bacterial community in bulk and rhizosphere soils
Luteolibacter, Candidatus genus Rhizospheria and acidobacteria/ Holophaga numbers, expressed as percentage of total bacteria, were between 2.00 . 10 -5 and . 91.6 over all groups (Figure 3). The relative abundances of subdivisions 1, 3, 4 and 6 acidobacteria in all three rhizospheres were equal to, or significantly lower, than those in corresponding bulk soils. For Luteolibacter, Candidatus genus Rhizospheria and Holophaga, the relative abundances in the rhizospheres of grass and potato were also equal to, or significantly lower than those in corresponding bulk soil, the exception being the grass rhizosphere in July, where the fraction of Candidatus genus Rhizospheria was significantly higher than in bulk soil. However, in all leek rhizospheres, the relative abundances of these three groups were always significantly higher than in bulk soil, with one exception (i.e. Holophaga in September, whose relative abundance was equal to the one in bulk soil). This indicates that Luteolibacter, Candidatus genus Rhizospheria and Holophaga are competitive towards other bacteria in the rhizosphere of growing leek plants where their abundances, relative to the total bacterial community, is higher than in bulk soil.

Factors affecting different bacterial populations in field soil
The effects of sample type, time and pH as environmental variables on cell estimates of all studied groups, as species variables, were calculated by multivariate analysis (RDA). A total of 93.3% of all variation was explained by the first two RDA axes (Figure 4). The rhizosphere of leek versus that of grass was discriminatory for most of the different studied populations. Subdivisions 1 and 3 acidobacteria and Holophaga correlated with the leek rhizosphere, whereas bacteria and subdivision 6 acidobacteria correlated with the grass rhizosphere.
The factors soil pH and 'grass rhizosphere' correlated with each other, indicating that either one or both are discriminative for bacteria and subdivision 6 acidobacteria. Subdivision 4 acidobacteria, Luteolibacter and Candidatus genus Rhizospheria did not show strong correlations with grass or leek rhizospheres.

Selection of Holophagae and of strain CHC25 in experimental leek-soil microcosms
A Kuchenbuch-style experimental plant-soil microcosms, the same used to assess the rhizosphere competence of Luteolibacter and Candidatus genus Rhizospheria strains [12], was applied to assess the competence of indigenous Holophagae (non-sterile soil) and of Holophagae strain CHC25 (following introduction into sterilized soil).
In non-sterilized non-inoculated soil, the average Holophaga cell number estimate (expressed as log Ceq per g dry soil) in the zone between 0 and 2 mm from the membrane that separated leek roots from the soil was 4.90 (range between 4.87 and 4.95). In a zone beneath, between 10 and 12 mm, the average log cell estimate was significantly (n=3, P≤ 0.05) lower, i.e. 4.24 (4.19 -4.29) ( Figure 5). This indicates that Holophagae naturally present in the soil increase in number when proximate to leek roots, confirming the observations made in the field.
Upon introduction into sterilized soil without leek plants, the strain CHC25 cell numbers persisted between estimated average log values of 4.32 (after 1 d) and 5.14 (after 21 d) for the duration of the experiment (35 days). In sterilized soil with added strain CHC25 cells planted with leek, the average cell number estimate in the 0-2 mm zone was 5.90 (5.86 -5.93), i.e. significantly higher than in the same layer of the system without growing leek plants (5.24, range between 5.03 -5.27). This number was also significantly higher than in the 10-12 mm layers of both systems (5.14, range between 4.91 -5.36). In the experiment in which a one-centimeter layer of sterilized non-inoculated soil was placed between the strain CHC25- inoculated soil and the membrane with leek roots, the estimated average log value in the 0 -2 mm layer after 35 days was 5.24 (5.08 -5.45) per g. In the same system without leek roots, this value was at background level, i.e. 3.52 (3.43 -3.65) per g. This background estimate was statistically indistinguishable from the one measured in sterilized noninoculated soil without plants (3.57, range between 3.39 -3.75). These data strongly indicate that the strain CHC25 cells in soil are selected by leek roots at a distance of at least 10 mm in soil. Likely, cells migrated to direct influence of leek roots (rhizosphere/rhizoplane), where they occupied available sites and multiplied. We posit, therefore, that particular member of class Holophagae, as represented by strain CHC25, should be considered as rhizosphere-competent.

Discussion
A field experiment was designed to explore the responses of five subdivisions of Acidobacteria and of two distinct groups within Verrucomicrobia subdivision 1 to the roots of different plant species. The selected subdivisions were found to be erratically present in one or more of the rhizospheres studied. Strikingly, we obtained compelling evidence for the contention that members of the Holophagae are competent in the leek rhizosphere. Leek rhizosphere competence has been shown before for Verrucomicrobia subdivision 1, exemplified by Candidatus genus Rhizospheria [12]. Hence, particular Acidobacteria can be common in rhizospheres, which is consistent with earlier reports on the presence of members of this phylum in the rhizospheres of Lolium perenne and Trifolium repens [41], Lodgepool pine [9], different grasses (Stipa hymenoides and Hilaria jamesii) [42], taxus [43], Thlaspi goesingense [44], chestnut [6] and oilseed rape (Brassica napus) [45]. This also implies that the roles of Acidobacteria in the rhizosphere can be complex. Further work will need to address the precise roles of different Acidobacteria that are found to be competent in the rhizospheres of particular plant species.
Our data are relevant for the current understanding on the interactions of predominant soil bacteria with the roots of different plants, as hardly anything is known about the association of the often numerically dominant members of the Acidobacteria and Verrucomicrobia with these. Representatives of both groups are often difficult to culture and hence most ecological studies in plant-soil environments have been performed with molecular tools that target entire phyla, thereby ignoring the behavior of specific subgroups within such phyla. An important message from this study, and the ones of Nunes da Rocha et al. [12] and Navarrete et al [22], is that, given their widely divergent ecological behavior, more attention needs to be paid to the behavior of the individual groups within the Acidobacteria and Verrucomicrobia, e.g. via isolation and re-introduction strategies. The strains that are isolated can be used for studying interactions with plants under selected experimental conditions. Validation of our subdivision and subgroup-specific quantification systems by making use of cultivable representatives of the different groups, allowed us to proximate actual cell number in the studied soil compartments over time. In the assumption that maximally one to two 16S ribosomal gene copies will be present in the genomes of different Verrucomicrobia [46] and Acidobacteria groups [47,48], cell equivalent numbers may proximate actual cell numbers if the genome numbers per cell remain constant for the different groups in the different soil compartments. So far, it is unknown to which extent the genome copy number per cell of the typical rhizosphere-responsive groups, such as Candidatus genus Rhizospheria and Holophagae, increases in the neighborhood of plant roots. Eventual increases in genome copy numbers in these groups may lead to an over-estimation of cell numbers near leek roots. Other confounding factors like presence of plant-derived (chloroplast) DNA in rhizosphere extracts can be excluded to influence bacterial quantities in different rhizospheres. Namely, no plant-specific amplicon sequences were found upon PCR amplification of rhizosphere soil DNA with the same bacterial primers as was used in our study [49], and only a small fraction of amplicons of nonbacterial origin were found after bacterial PCR amplification and high throughput sequencing from rhizosphere soil processed according the same procedure as applied in our study [50].
Remarkably, we found evidence supporting the fact that Holophagae as group, or a particular subset thereof, specifically responded to leek roots by an increase in 16S rRNA gene copy number and not to the ones of potato and grass growing in the same field. Moreover, removal of the grass layer covering the field resulted in an increase in the Holophagae 16S rRNA gene copy number later during the season in bulk soil, which allows the hypothesis that grass roots can be suppressive towards Holophagae. Grass and leek both are monocotyledonous plant species and hence the preference of Holophagae for cannot be explained along the monocot/dicot dichotomy. The lack of a stimulatory effect of the dicot potato in the field indicated that local conditions established by the roots were not propitious to holophagal cell growth. This in spite of the fact that pH in the potato rhizosphere was indistinguishable from that in the leek rhizosphere.
Our observations thus shed new light on the lifestyles of particular soil Acidobacteria. Acidobacteria commonly are considered to encompass mainly oligotrophic or K strategist forms [15,24,25,51]. The likely presence of low-specificity / high-affinity substrate uptake systems, as evidenced from analyses of the genomes of Acidobacteria subdivision1 and subdivision 3 strains [48], may indicate that these strains indeed exhibit oligotrophy as a major lifestyle in soil. This stands in sharp contrast to the here-defined Holophagae lifestyle, which was clearly responsive to leek roots, either by increased cell division and/ or by increase in genome quantity per cell, showing typical r-strategist behavior. As is the case for many other lineages within the bacterial domain [24], a clear niche differentiation exist among species of the phylum Acidobacteria.
The rhizosphere of leek thus appears to represent a specific niche for the Holophagae species that were studied. Two strains, CHC25 and ORAC (> 1300 bp stretches of the rRNA gene sequences were deposited in the EMBL Nucleotide Sequence Database, respectively, under accession numbers FN554392 and FN689719) , were able to grow on simple organic acids common in root exudates like oxalic acid, malic acid, succinic acid and citric acid [52]. Both strains closely resembled each other and substantially differed in taxonomy and physiology from Holophagae strains Geothrix fermentans H-5 [53] and Holophaga foetida TMBS4 [54]. Whereas the latter two strains are obligatory anaerobic, our strains were aerobic [52]. In their physiologies, strains CHC25 and ORAC resembled the aerobic Holophagae strain Acanthopleuribacter pedis FYK2218 T [55]. Cells of this strain are also motile as was the case for our strains CHC25 and ORAC [52]. However, the taxonomical distance between our strains and A. pedis strain FYK22218 T was larger than with G. fermentans H-5/ H. foetida TMBS4 [32]. Thus, strains CHC25 and ORAC potentially represent a new species, clearly distinct from any previously described Holophagae species. We tentatively propose the name for strains CHC25 and ORAC as Candidatus Porrumbacterium oxyphilus, a leek bacterium that prefers oxygen; the ecologically-defined traits that clearly distinguish strains CHC25 and ORAC from other strains within the class Holophagae.
In conclusion, we found specialized groups within several subdivisions of Verrucomicrobia and Acidobacteria that are rhizosphere competent. Their lifestyles with the plant may suggest that these bacteria either interact with the plants themselves or with communities associated with plants. This novel insight extends our current understanding of bacteria that associate with plants and may be a basis for further exploration of their putative roles in other habitats.