http://orcid.org/0000-0001-8903-8053
Figures
Abstract
The Acidobacteria is one of the major bacterial phyla in soils and peatlands. The currently explored diversity within this phylum is assigned to 15 class-level units, five of which contain described members. The ecologically relevant traits of acidobacteria from different classes remain poorly understood. Here, we compared the patterns of acidobacterial diversity in sandy soils of tundra, along a gradient of increasing vegetation–unfixed aeolian sand, semi-fixed surfaces with mosses and lichens, and mature soil under fully developed plant cover. The Acidobacteria-affiliated 16S rRNA gene sequences retrieved from these soils comprised 11 to 33% of total bacterial reads and belonged mostly to members of the classes Acidobacteriia and Blastocatellia, which displayed opposite habitat preferences. The relative abundance of the Blastocatellia was maximal in unfixed sands and declined in soils of vegetated plots, showing positive correlation with soil pH and negative correlation with carbon and nitrogen availability. An opposite tendency was characteristic for the Acidobacteriia. Most Blastocatellia-affiliated reads belonged to as-yet-undescribed members of the family Arenimicrobiaceae, which appears to be characteristic for dry, depleted in organic matter soil habitats. The pool of Acidobacteriia-affiliated sequences, apart from Acidobacteriaceae- and Bryobacteraceae-related reads, had a large proportion of sequences from as-yet-undescribed families, which seem to specialize in degrading plant-derived organic matter. This analysis reveals sandy soils of tundra as a source of novel acidobacterial diversity and provides an insight into the ecological preferences of different taxonomic groups within this phylum.
Citation: Ivanova AA, Zhelezova AD, Chernov TI, Dedysh SN (2020) Linking ecology and systematics of acidobacteria: Distinct habitat preferences of the Acidobacteriia and Blastocatellia in tundra soils. PLoS ONE 15(3): e0230157. https://doi.org/10.1371/journal.pone.0230157
Editor: Ying Ma, Universidade de Coimbra, PORTUGAL
Received: August 20, 2019; Accepted: February 22, 2020; Published: March 17, 2020
Copyright: © 2020 Ivanova et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The 16S rRNA gene sequence dataset used in this study is deposited with GenBank under the Bioproject accession number PRJNA497067 and is publicly available.
Funding: This study was supported by the Russian Science Foundation (RSF -http://www.rscf.ru/en/ Field studies and 16S rRNA gene sequencing were performed by A.Z. and T.C. and supported by the project No 17-16-01057. Bioinformatic and taxonomic analyses were performed by A.I. and S.D. and were supported by the project No 16-14-10210. The funder had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
The Acidobacteria is one of the most abundant and highly diverse bacterial phyla in soils and peatlands [1–6]. The proportion of Acidobacteria-affiliated 16S rRNA gene reads in sequence pools retrieved from various soil habitats ranges between 5 and 50% of the total bacterial community [3,7–10]. Our knowledge of the roles of acidobacteria in soils includes decomposition of various biopolymers and participation in the global cycling of carbon, iron and hydrogen, but this list of functional capabilities remains far from being complete and is attributed to several sub-groups of this phylum only.
The currently explored diversity within the Acidobacteria is commonly addressed as corresponding to 26 major 16S rRNA gene sequence clades or subdivisions (SD) [11]. Recently, these 26 subdivisions were assigned to 15 class-level units, five of which contain described members [12]. These include three earlier established classes Acidobacteriia, Blastocatellia and Holophagae [13–15] as well as two recently proposed classes Vicinamibacteria and Thermoanaerobaculia [12]. The phylogenetic range of the class Acidobacteriia accommodates 16S rRNA gene sequences from several SDs including 1, 2, 3, 5, 11, 12, 13, 14, 15, and 24. The classes Blastocatellia, Vicinamibacteria and Thermoanaerobaculia correspond to one subdivision each, i.e. SDs 4, 6, and 23, while the class Holophagae includes SDs 8 and 22. In soil, acidobacterial SDs 1, 2, 3, 4 and 6 are the most abundant ones [8]. This is true for a wide range of ecosystem types, including boreal and tropical forests, grasslands and pastures, as well as arid landscapes [3–5, 8–10]. The abundances of subdivisions 1, 2 and 3 show negative correlation with soil pH, while the opposite tendency is characteristic for subdivisions 4 and 6 [8]. Peatlands are also among the preferred habitats of Acidobacteria. Acidic Sphagnum-dominated boreal peatlands are colonized mainly by members of SDs 1 and 3 [16–18]. Subarctic peatlands, in addition to SDs 1 and 3, may contain a relatively high proportion of SD2 Acidobacteria [19,20].
Despite their wide distribution in various soils, acidobacteria remain strongly underrepresented in culture collections due to difficulties in their cultivation and laboratory maintenance [6]. Much of the currently described diversity (ca. 15 genera) belong to the class Acidobacteriia, which includes the orders Acidobacteriales and Bryobacterales (SDs 1 and 3), and accommodates acidophilic or acidotolerant, mesophilic and psychrotolerant, chemoheterotrophic bacteria that utilize various sugars and polysaccharides and possess a number of hydrolytic capabilities [12]. The class Blastocatellia contains 7 genera of aerobic, mesophilic or thermophilic, chemo-heterotrophic bacteria that specialize on degradation of complex proteinaceous compounds and one genus of microaerophilic thermophilic anoxygenic photoheterotrophs, Chloracidobacterium [21]. Other classes of Acidobacteria contain only a limited number of characterized representatives and, therefore, their ecologically relevant traits remain poorly understood.
In the course of our recent study of prokaryotic community succession following sand fixation and soil formation in the tundra zone, Acidobacteria were identified as one of the most abundant bacterial groups, which was present along a whole gradient of increasing vegetation, from unfixed aeolian sands to mature soils under fully developed plant cover [22]. Since the phenotype of bacteria that colonize dry and nutrient-poor sands is clearly different from that of bacteria in organic-rich soils under plant cover, the acidobacterial populations in these sites with contrasting characteristics should have been quite distinct with regard to their lifestyles and environmental adaptations. The analysis presented here was undertaken in order to verify this hypothesis and to examine diversity patters of the Acidobacteria along a gradient of increasing vegetation in tundra.
Materials and methods
Sampling sites
The 16S rRNA gene sequence dataset retrieved by Zhelezova et al. [22] and deposited with GenBank under the Bioproject accession number PRJNA497067 was used for the analysis. Two chronosequences of soil formation on aeolian sands with similar initial stages and different mature vegetation (typical tundra and wooded tundra) corresponded to the two sampling sites located on the shores of the Pechora River (Northwestern Russia, Nenetsia region). The two sites differed in the plant cover that developed on mature soil, i.e. typical tundra vegetation with subshrubs (Site I) and wooded tundra with rare trees and subshrubs (Site II). These sampling sites are shown in Fig 1 and the characteristics of soil samples used for the analysis are given in Table 1. Detailed description of the sampling sites are reported elsewhere [22]. No specific permits were required for sampling at these tundra sites. These locations are not privately-owned or protected in any way and are also not parts of national parks or reserves. Our sampling did not involve endangered or protected species.
(a) flat sand hills near Nelmin Nos (site I) and (b) small sand dunes near Naryan-Mar (site II). White flags indicate different types of plots along a gradient of increasing vegetation: US—unfixed sand, SF—semi-fixed surface, MS—mature soil. (c) Example of vegetation gradient from sand with gravel to cover of lichens and mosses (left to right) in site I near Nelmin Nos.
For both sites, sampling was performed on three types of surfaces: 1- unfixed aeolian sand (US), 2- semi-fixed surface with mosses and lichens (SF), and 3—mature soil under developed plant cover (MS) (Fig 1, Table 1). Sampling plots of different types were located on a transect with 3–5 m between each plot. For every surface type on each site, five soil samples (~50 g of soil free of plants, mosses and lichens) were taken from depths of 1–5 cm. For molecular analyses, samples were stored at -70 ˚C. The total organic carbon (TOC) and total nitrogen (TN) contents were determined for the pooled sample of the five replicates taken for the molecular analysis from each plot using a Vario MACRO Cube CN-analyzer (Elementar Analysensysteme GmbH, Germany) [22]. pH was measured in soil suspension in distilled water in the ratio of 1:2.5.
Analysis of SSU rRNA gene sequences
The sequences retrieved by Zhelezova et al. [22] were represented by the fragments of V3–V4 regions of bacterial and archaeal 16S rRNA genes. This pool of sequences was reanalyzed with QIIME 2 v.2018.8 (https://qiime2.org) [23]. DADA2 plugin was used for sequence quality control, merging of paired-end reads and chimera filtering [24]. Since the number of reads obtained from two replicates of MS plot in site I was too low (below 5.000 sequences), the comparative analysis of all three plots in site I was performed based on three replicates.
Operational Taxonomic Units (OTUs) were clustered applying VSEARCH plugin [25] with open-reference function using db Silva132 [26,27] with 97% identity. Taxonomy assignment was performed using BLAST against db Silva132 with 97% identity. 16S rRNA sequences affiliated with the Acidobacteria were extracted and used for further detailed analyses. The alpha-diversity indices were calculated using the core-metrics-phylogenetic method implemented in QIIME 2 v.2018.8. Phylogenetic analysis was carried out using the ARB program package [28] (version 6.0.3) and the sequence alignment based on dbSilva132. Short sequences were added to the existing tree using the quick ARB parsimony insertion tool. The network diagram was constructed using Gephi [29].
Statistical analyses
Statistical evaluations were made with GraphPad Prism (v. 7.0) applying multiple t-tests. Corrections for multiple comparisons were made using Holm-Sidak method. The significance level alpha was set at 0.05. Pearson correlation test was performed to check correlations between soil chemical properties and sequence abundances of different acidobacterial groups.
Results
Characteristics of sampled substrates along a gradient of increasing vegetation
The sampling plots corresponding to the three successional stages of soil formation in tundra—unfixed aeolian sand, semi-fixed surfaces with mosses and lichens, and mature soil under fully developed plant cover—differed from each other with regard to various physico-chemical parameters (Table 1). Most pronounced differences were observed with regard to moisture and total organic carbon contents. Unfixed and semi-fixed sands were extremely dry, while soils under plant cover contained more water. The content of organic carbon was nearly non-measurable in unfixed sands and increased gradually with formation of vegetation cover.
Community composition of the Acidobacteria
A total of 1,019,619 partial 16S rRNA gene sequences (mean amplicon length 260 bp) were retrieved from the examined soil samples (Table 2). Of these, 463,507 sequences were retained after merging of paired-end reads, quality filtering, removing chimeras and singletons. The pool of Acidobacteria-affiliated reads included 74037 sequences, which accounted for 9–31% of all bacterial reads in different samples. Overall, the acidobacterial alpha-diversity in unfixed aeolian sands was slightly lower than that in soils of vegetated plots (see mean values of Shannon indexes in Table 2) but this difference was not statistically significant.
As revealed by principle coordinate analyses, the acidobacterial communities corresponding to the three successional stages of soil formation in tundra were also clearly distinct (Fig 2).
PCoA plot is based on the weighted UniFrac distance of the sequencing dataset.
The Acidobacteria-affiliated sequences retrieved from the two tundra sites belonged to members of the classes Acidobacteriia (53–100% of all acidobacterial reads in different samples), Blastocatellia (0–34%), and Holophagae (0–13%) (Fig 3). The proportion of sequences that could not be assigned to these classes was in the range of 0–6% in unfixed sands and 0–1% in mature soils under plant cover. The pool of sequences from the Acidobacteriia was composed of reads from representatives of the orders Acidobacteriales (SD1) and Bryobacterales (SD3) as well as SD2 acidobacteria, which comprise an as-yet-undescribed order of this class. Notably, a large proportion of Acidobacteriales-affiliated sequences (6–37% of all acidobacterial reads in different samples) could not be assigned to the only currently described family of this order, i.e. Acidobacteriaceae. This group of sequences is further addressed as ‘uncultured Acidobacteriales’ in accordance with the taxonomic classification used in db Silva132 and indicated by black hatching in Fig 3.
A–site I, B–site II. The taxonomic analysis was performed according to dbSilva 132. The values of the relative abundance of different acidobacterial groups in individual samples are given in S1 Table. Significant differences in relative abundances of particular acidobacterial groups between unfixed sand and different stages of soil formation were revealed for the Acidobacteriaceae (P-value < 0.01), Bryobacterales (P-value < 0.001) and Blastocatellia (P-value < 0.001) (S2 Table).
Taxonomy-based analysis revealed a clear shift in the Acidobacteria community composition along a gradient of increasing vegetation in both study sites (Fig 3). The most characteristic feature of acidobacterial communities in unfixed sands was a prominent presence of the Blastocatellia (SD4), which comprised 31±2% and 14±5% of all Acidobacteria-affiliated reads obtained from sites I and II, respectively (S1 Table). The relative abundance of these bacteria declined dramatically in vegetated soils of both sites, down to 7±4% of reads in mature soils from the site I and a complete disappearance in mature soils from site II. The same trend was observed for members of the Holophagae, which accounted for 12±1% and 7±3% of all acidobacterial reads in unfixed sands from sites I and II, respectively, but were not detected in mature soils (Fig 3). The opposite trend was detected for members of the Acidobacteriaceae and SD2, which were either absent (Acidobacteriaceae in site I) or present in a low abundance in unfixed sands (13±5% for Acidobacteriaceae in site I and 4–7% for SD2 in both sites) but became major community members (with relative abundances of approximately 20–30%) in vegetated mature soils (S1 Table). The Bryobacterales-affiliated acidobacteria (SD3) displayed highest relative abundances in semi-fixed soils, accounting for 36±9% and 50±4% of all acidobacterial reads in sites I and II, respectively. No clear trend could be observed for the uncultured group within the order Acidobacteriales (Fig 3).
Most abundant OTUs of Acidobacteria and their distribution patterns
Using 97% sequence identity, a total of 232 acidobacterial OTUs were identified in all studied samples. Of these, 157 OTUs were detected in the site I and 159 OTUs were identified in the site II (Table 2); 80 OTUs were shared between sites I and II (S1 Fig). The numbers of OTUs shared between the two sites within each of the succession stages were 23, 24 and 42 for US, SF and MS plots, respectively (S1 Fig).
Our further, detailed analysis was focused on five major taxonomic groups within the Acidobacteria. Two of these groups, Blastocatellia (SD4) and Acidobacteriaceae (SD1), were chosen because they displayed two clearly opposite trends along a gradient of increasing vegetation in tundra soils. These groups of Acidobacteria contain a number of cultured and characterized representatives, thus offering a possibility of analyzing specific reasons behind their environmental distribution. Two other groups, SD2 and uncultured Acidobacteriales, do not contain described representatives and were chosen, therefore, in order to find out their eco-niche preferences. One additional group that was included in the analysis, Holophagae (SD8), contains only few cultured representatives; the physiology and functional capabilities of these acidobacteria remain poorly understood. The OTUs of these bacteria comprising ≥1% of all Acidobacteria-affiliated reads in at least one of the examined sites are listed in Table 3 and are displayed in S2 Fig.
Although the same taxonomic representation of Acidobacteria at the order and family levels was observed for both sites, diversity of OTUs identified in the sites I and II was clearly different (S2 Fig). For example, Blastocatellia-affiliated OTUs 1–3, which were among the most abundant OTUs in site I, were absent from site II. Vice versa, the most abundant Blastocatellia-affiliated OTU in site II, OTU 7, was not detected in site I. Other groups of Acidobacteria in the two sites were also represented by different OTUs, which may be due to some differences between the sites I and II, such as the plant cover composition and the moisture level (Table 1).
The most abundant OTUs of the Blastocatellia (OTUs 1–3 and 7–9) belonged to the as-yet-uncultivated clade RB41 within the family Pyrinomonadaceae and occurred mainly in unfixed sands in both tundra sites. Only three OTUs were classified at the genus level as representing the genera Blastocatella (OTUs 11 and 12) and Stenotrophobacter (OTU 5).
By contrast, nearly all of the Acidobacteriaceae-affiliated OTUs could be classified at the genus level and represented the genera Granulicella (OTUs 24–28, 30, 33, 34), Acidipila (OTUs 23, 32) and Bryocella (OTU 29). The most abundant OTU from the Acidobacteriaceae (No 31), however, could not be assigned to any of the described genera and displayed 96.0–96.5% similarity with 16S rRNA gene sequences from ‘Acidisarcina’/Acidipila group.
The most abundant SD2-affiliated OTUs (No 13–15) displayed high similarity with environmental 16S rRNA gene sequences retrieved from soils or sediments of various northern locations (GenBank numbers FJ004757, KM200371, FJ625349). These were detected mainly in mature soils, although OTU 13 was also highly abundant in semi-fixed sand from site I.
Finally, OTUs comprising the group of uncultured Acidobacteriales were present in all examined samples but varied with regard to their identity and relative abundance. Thus, OTUs 38 and 48 were abundant in unfixed sands but were not detected in semi-fixed sands and soils under plant cover. By contrast, OTU 39 was highly representative in semi-fixed sands but was also present in other samples. Environmental 16S rRNA gene sequences representing this group of acidobacteria were retrieved mainly from various grasslands and forest soils but also from microbial mats of lava tube walls, orthoquartzite caves, uranium and potassium mines (Table 3).
Major trends of changes in the Acidobacteria community composition in sandy soils of tundra along the gradient of increasing vegetation are seen in the network diagram illustrating the most abundant OTU distribution between the three examined types of plots (Fig 4). An apparent shift from the community composed of the Blastocatellia and uncultured Acidobacterales in aeolian sands to the community dominated by various subgroups of the Acidobacteriia in mature soils is observed in this diagram.
The size of the OTU nodes is weighted according to the relative abundance of the particular OTU. The diagram was constructed using gephi [29].
Specific phylogenetic sub-groups of Acidobacteria revealed in tundra soils
A large group of OTUs affiliated with the class Blastocatellia, which were identified mostly in unfixed sands (OTUs 1–4, 6–10), were classified as belonging to the as-yet-uncultivated clade RB41 within the family Pyrinomonadaceae (Table 3). This group was named after the environmental clone sequence RB41 (GenBank accession No Z95722), which was retrieved from the Roggenstein field site soil near Munich (Germany) in the seminal work of Ludwig et al. [1]. Detailed phylogenetic analysis of these sequences revealed that they form a distinct phylogenetic group, which clusters separately from the other described families of the order Blastocatellales, i.e. Blastocatellaceae, Arenimicrobiaceae and Pyrinomonadaceae (Fig 5). These sequences displayed 91–92% similarity to 16S rRNA gene sequences from members of the genera Brevitalea and Arenimicrobium, while their sequence similarity to other described members of the Blastocatellia was in the range of 86–89%.
Various families within the Blastocatellia are indicated as follows: 1—Blastocatellaceae, 2 - ‘Chloracidobacteriaceae’, 3—Arenimicrobiaceae, 4—Pyrinomonadaceae. Bootstrap values are derived from 1000 pseudoreplicates. An outgroup was composed of three 16S rRNA gene sequence from members of the Holophagae, Holophaga foetida (X77215) and two related environmental sequences (FQ658676 and FQ659446). The scale bar indicates 10% estimated sequence divergence.
Another large group of numerically abundant OTUs identified in tundra soils affiliated with the order Acidobacteriales but did not belong to the only currently described family Acidobacteriaceae (Fig 6). In db Silva 132, this group is addressed as ‘uncultured Acidobacteriales’, which is not fully correct since it comprises a number of isolates named “Ellin’ (Ellin7137, 7522, 6547, 6528, 6527, 5106, 323 and others). These isolates were obtained by Janssen and co-workers [30–33]; none of them, however, was characterized, so that no information about these bacteria is currently available. 16S rRNA gene sequence similarity of these bacteria to described members of the family Acidobacteriaceae is below 91%.
Various families and groups within the Acidobacteriia are indicated as follows: 1—Acidobacteriaceae, 2 –‘Koribacteraceae’, 3 –uncultured Acidobacteriales, 4 –SD2. Bootstrap values are derived from 1000 pseudoreplicates. An outgroup was composed of three 16S rRNA gene sequence from members of the Holophagae, Holophaga foetida (X77215) and two related environmental sequences (FQ658676 and FQ659446). The scale bar indicates 10% estimated sequence divergence.
Correlations between chemical soil properties and the abundances of different acidobacterial groups
To reveal the influence of several chemical soil properties (pH, C and N contents) on the relative abundances of different acidobacterial groups, Pearson correlation test was performed. The positive correlation was revealed between pH and the relative abundances of the Holophagae (SD8) (R = 0.82, p = 0.05) and Blastocatellia (SD4) (R = 0.80, p = 0.06). The highest relative abundance of these groups was detected in the acidobacterial communities of unfixed sands in both sites. The negative correlation was revealed between pH and the relative abundances of the Acidobacteriaceae (SD1) (R = -0.83, p = 0.04) and SD2 (R = -0.73, p = 0.1), which tended to dominate in mature soils. These correlations could be explained by the gradual decrease of pH during the process of soil formation. No significant correlations, however, were revealed between C and N contents and the abundance of particular acidobacterial groups (Table 4).
Discussion
All previously available knowledge of Acidobacteria diversity in the tundra zone refers to the ecosystems with developed plant communities, i.e. tundra heaths with a mosaic vegetation of dwarf shrubs and alpine grasses [19,34], wetlands with mixed cover of lichens and mosses [20] or forested tundra with lichen cover [35]. All of these studies report SDs 1, 2 and 3 (belonging to the class Acidobacteriia) as the major groups of tundra-inhabiting Acidobacteria. The results obtained in our study for vegetated plots (Fig 3) are in agreement with the previous reports. The search for information on Acidobacteria diversity in aeolian sand dunes and the corresponding chronosequence of soil formation on sands in the tundra zone, however, yielded no results. As revealed in our study, sandy soils of tundra are characterized by a prominent presence of the Blastocatellia (SD4), which may comprise up to one third of the acidobacterial community. This finding was somewhat unexpected because none of the currently described members of the Blastocatellia were obtained from cold environments. Thermophilic representatives of this class, Pyrinomonas methylaliphatogenes (family Pyrinomonadaceae) and Chloracidobacterium thermophilum (family ‘Chloracidobacteriaceae’), were isolated from a hot spring and a thermal soil, respectively [36,37]. All other described representatives, i.e. members of the genera Blastocatellia, Aridibacter, Tellurimicrobium, Stenotrophobacter (family Blastocatellaceae) and the genera Arenimicrobium, Brevitalea (family Arenimicrobiaceae), were obtained from Namibian semiarid savanna soils [14,38,39]. These bacteria were characterized as possessing wide growth temperature ranges (from 8–11°C to 45–52°C) and, yet, their growth optima were recorded at 28–45°C. The 16S rRNA gene sequences retrieved from aeolian sands of tundra affiliated only with the families Arenimicrobiaceae and Blastocatellaceae. Most of these sequences displayed low similarity (below 91–92%) to those from described members of these families, suggesting that the phenotypes of these bacteria are also different from those of acidobacteria from Namibian savanna soils. One particular group of sequences characteristic for sands of tundra, which were classified by using dbSilva 132 as ‘uncultured Pyrinomonadaceae RB41’ (Table 3), formed a common clade with 16S rRNA gene sequences from the Arenimicrobium and Brevitalea (Fig 5), suggesting that these sequences should be addressed as affiliating with the family Arenimicrobiaceae. Final conclusions about the taxonomic position of these bacteria, however, should await their isolation and characterization. Future efforts in culturing these bacteria could possibly benefit from the evidence for their ability to survive drought, nutrient limitation and low temperatures (as suggested by Wüst et al. [40] and our study).
The trend observed for members of the family Acidobacteriaceae (Table 4) agrees well with our current knowledge of their preference for low pH and high availability of plant-derived organic matter [4,6,8,10]. Notably, a large group of OTUs in semi-fixed sands or mature soils were represented by Granulicella species (OTUs 24–28, 30, 33, 34). These bacteria were isolated from tundra soils or northern boreal wetlands [40–42] and are commonly associated with mosses and lichens. Good tolerance of low temperatures, wide repertoires of carbohydrate-active enzymes encoded in their genomes and pronounced hydrolytic capabilities explain wide distribution of Granulicella species in tundra habitats [42–45]. The OTU31, which was shared between SF and MS plots in site I and was highly abundant in MS plots of site II, is represented by ‘Acidisarcina’- like bacteria, which were also isolated from lichen-dominated tundra soils and displayed chitinolytic and xylanolytic capabilities [35].
A notable group of OTUs in unfixed sands was represented by Holophagae-affiliated 16S rRNA gene sequences (OTUs 61–68). This finding was somewhat unexpected because these acidobacteria are commonly not abundant in soils [9]. At present, the class Holophagae contains only four described representatives [12]. Three of these acidobacteria are strict anaerobes and none of them was obtained from soils. 16S rRNA gene sequences retrieved in our study from upper, aerobic layers of sandy tundra soils were only distantly related (83–86% sequence similarity) to those of described Holophagae members. Most likely, these sequences belong to aerobic, psychrotolerant bacteria with as-yet-unknown functional potential.
An insight into community changes of the Acidobacteriia may also provide valuable information regarding some as-yet-uncultivated sub-groups within this class. One example is SD2 Acidobacteria which, by now, does not contain characterized representatives. Similar to members of the Acidobacteriaceae, the relative abundance of SD2 Acidobacteria was negatively correlated with pH and positively correlated with organic carbon content (Table 4). However, their correlation with nitrogen availability was more strongly pronounced than that of Acidobacteriaceae, thus suggesting the need for increasing nitrogen content in cultivation media. Another numerically abundant group of Acidobacteriia-affiliated 16S rRNA gene sequences, which were retrieved from all examined samples and addressed as ‘uncultured Acidobacteriales’, was classified as belonging to an as-yet-undescribed family (Fig 6). Phylogenetic diversity within this group is relatively wide and no clear correlation was observed with any of the examined parameters (Table 4). As noticed above, several representatives of this group (isolates named Ellin7137, 7522, 6547, 6528, 6527, 5106, 323) were earlier obtained and reported as being represented by extremely slow-growing bacteria that formed mini-colonies (with diameters of 25–200 μm) after long (> 12 weeks) periods of incubation [33]. It is hardly surprising that none of these isolates were characterized and taxonomically described. The occurrence of these slow-growing and presumably oligotrophic bacteria in aeolian sands of tundra appears to be logical based on their phenotype. Cultivation and characterization of these microorganisms represents a challenge for further taxonomic studies on Acidobacteria.
Conclusions
In summary, our analysis revealed two clearly distinct profiles of ‘ecological fitness’ of the Acidobacteriia and Blastocatellia. The latter, in particular members of the family Arenimicrobiaceae, appear to be characteristic for dry, depleted in organic matter sandy soils. An opposite habitat preference was demonstrated by the Acidobacteriia, which seem to specialize in degrading plant-derived organic matter and, therefore, become a major group of acidobacteria in soils under fully developed plant cover. This linkage between the taxonomic affiliation and the potential functional capabilities could help interpreting the results of molecular diversity surveys and navigate our efforts in obtaining representatives of as-yet-uncultivated groups of soil Acidobacteria.
Supporting information
S1 Table. Relative abundance of different acidobacterial groups.
https://doi.org/10.1371/journal.pone.0230157.s001
(PDF)
S2 Table. Statistical values calculated by multiple t-test comparisons for acidobacterial groups from Sites 1 and 2.
https://doi.org/10.1371/journal.pone.0230157.s002
(PDF)
S1 Fig. Venn diagrams showing the number of shared and unique OTUs in Sites 1 and 2, and in three types of successional stages examined in this study.
https://doi.org/10.1371/journal.pone.0230157.s003
(PDF)
S2 Fig. Heat map showing the relative abundances of well-represented acidobacterial OTUs.
https://doi.org/10.1371/journal.pone.0230157.s004
(PDF)
References
- 1. Ludwig W, Bauer SH, Bauer M, Held I, Kirchhof G, Schulze R, et al. Detection and in situ identification of representatives of a widely distributed new bacterial phylum. FEMS Microbiol Lett. 1997;153: 181–190. pmid:9252585
- 2. Hugenholtz P, Pitulle C, Hershberger KL, Pace NR. Novel division level bacterial diversity in a Yellowstone hot spring. J Bacteriol. 1998;180: 366–76. Available: http://www.ncbi.nlm.nih.gov/pubmed/9440526 pmid:9440526
- 3. Janssen PH. Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Applied and Environmental Microbiology. 2006. pp. 1719–1728. pmid:16517615
- 4. Kielak AM, Barreto CC, Kowalchuk GA, van Veen JA, Kuramae EE. The ecology of Acidobacteria: moving beyond genes and genomes. Frontiers in Microbiology. 2016. p. 744. pmid:27303369
- 5. Eichorst SA, Trojan D, Roux S, Herbold C, Rattei T, Woebken D. Genomic insights into the Acidobacteria reveal strategies for their success in terrestrial environments. Environ Microbiol. 2018;20: 1041–1063. pmid:29327410
- 6.
Dedysh SN, Sinninghe Damsté JS. Acidobacteria. eLS. Chichester, UK: John Wiley & Sons, Ltd; 2018. pp. 1–10. https://doi.org/10.1002/9780470015902.a0027685
- 7. Lee SH, Ka JO, Cho JC. Members of the phylum Acidobacteria are dominant and metabolically active in rhizosphere soil. FEMS Microbiol Lett. 2008;285: 263–269. pmid:18557943
- 8. Jones RT, Robeson MS, Lauber CL, Hamady M, Knight R, Fierer N. A comprehensive survey of soil acidobacterial diversity using pyrosequencing and clone library analyses. ISME J. 2009;3: 442–453. pmid:19129864
- 9. Lauber CL, Hamady M, Knight R, Fierer N. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol. 2009;75: 5111–5120. pmid:19502440
- 10. Foesel BU, Nägele V, Naether A, Wüst PK, Weinert J, Bonkowski M, et al. Determinants of Acidobacteria activity inferred from the relative abundances of 16S rRNA transcripts in German grassland and forest soils. Environ Microbiol. 2014;16: 658–675. pmid:23802854
- 11. Barns SM, Cain EC, Sommerville L, Kuske CR. Acidobacteria phylum sequences in uranium-contaminated subsurface sediments greatly expand the known diversity within the phylum. Appl Environ Microbiol. 2007;73: 3113–6. pmid:17337544
- 12. Dedysh SN, Yilmaz P. Refining the taxonomic structure of the phylum Acidobacteria. Int J Syst Evol Microbiol. 2018; 3796–3806. pmid:30325293
- 13.
Thrash J, Coates J. Class I. Acidobacteriia Cavalier-Smith 2002, 12VP. In: Krieg NR, Staley JT, Brown DR, Hedlund BPBPNW, editor. Bergey’s Manual of Systematic Bacteriology, 2nd edn, vol 4. Springer, New York; 2011. p. 727.
- 14. Pascual J, Wüst PK, Geppert A, Foesel BU, Huber KJ, Overmann J. Novel isolates double the number of chemotrophic species and allow the first description of higher taxa in Acidobacteria subdivision 4. Syst Appl Microbiol. 2015;38: 534–544. pmid:26460220
- 15. Thrash J, Coates J. Class II. Holophagae classis nov. Bergey’s Manual of Systematic Bacteriology, 2nd edn, vol 4. 2011. p. 731.
- 16. Dedysh SN. Cultivating uncultured bacteria from northern wetlands: knowledge gained and remaining gaps. Front Microbiol. 2011;2: 184. pmid:21954394
- 17. Pankratov TA, Ivanova AO, Dedysh SN, Liesack W. Bacterial populations and environmental factors controlling cellulose degradation in an acidic Sphagnum peat. Environ Microbiol. 2011;13: 1800–1814. pmid:21564458
- 18. Serkebaeva YM, Kim Y, Liesack W, Dedysh SN. Pyrosequencing-based assessment of the bacteria diversity in surface and subsurface peat layers of a northern wetland, with focus on poorly studied phyla and candidate divisions. PLoS One. 2013;8: e63994. pmid:23700443
- 19. Männistö MK, Kurhela E, Tiirola M, Häggblom MM. Acidobacteria dominate the active bacterial communities of Arctic tundra with widely divergent winter-time snow accumulation and soil temperatures. FEMS Microbiol Ecol. 2013;84: 47–59. pmid:23106413
- 20. Danilova O V., Belova SE, Gagarinova I V., Dedysh SN. Microbial community composition and methanotroph diversity of a subarctic wetland in Russia. Microbiology. 2016;85: 583–591.
- 21. Huber KJ, Pascual J, Foesel BU, Overmann J. Blastocatellia. Bergey’s Manual of Systematics of Archaea and Bacteria. 2017.
- 22. Zhelezova A, Chernov T, Tkhakakhova A, Xenofontova N, Semenov M, Kutovaya O. Prokaryotic community shifts during soil formation on sands in the tundra zone. Esposito A, editor. PLoS One. Public Library of Science; 2019;14: e0206777. pmid:30939175
- 23. Caporaso J, Kuczynski J, Stombaugh J. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. Nature Publishing Group; 2010;7: 335–336.
- 24. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13: 581–583. pmid:27214047
- 25. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4: e2584. pmid:27781170
- 26. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41: D590–D596. pmid:23193283
- 27. Glöckner FO, Yilmaz P, Quast C, Gerken J, Beccati A, Ciuprina A, et al. 25 years of serving the community with ribosomal RNA gene reference databases and tools. J Biotechnol. 2017;261: 169–176. pmid:28648396
- 28. Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar A, et al. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32: 1363–1371. pmid:14985472
- 29. Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Third International AAAI Conference on Weblogs and Social Media. 2009. pp. 361–362.
- 30. Sait M, Hugenholtz P, Janssen PH. Cultivation of globally distributed soil bacteria from phylogenetic lineages previously only detected in cultivation-independent surveys. Environ Microbiol. 2002;4: 654–666. pmid:12460273
- 31. Joseph SJ, Hugenholtz P, Sangwan P, Osborne CA, Janssen PH. Laboratory cultivation of widespread and previously uncultured soil bacteria. Appl Environ Microbiol. 2003;69: 7210–7215. pmid:14660368
- 32. Davis KER, Joseph SJ, Janssen PH. Effects of growth medium, inoculum size, and incubation time on culturability and isolation of soil bacteria. Appl Environ Microbiol. 2005;71: 826–834. pmid:15691937
- 33. Davis KER, Sangwan P, Janssen PH. Acidobacteria, Rubrobacteridae and Chloroflexi are abundant among very slow-growing and mini-colony-forming soil bacteria. Environ Microbiol. 2011;13: 798–805. pmid:21108723
- 34. Männistö MK, Tiirola M, Häggblom MM. Bacterial communities in Arctic fjelds of Finnish Lapland are stable but highly pH-dependent. FEMS Microbiology Ecology. 2007. pp. 452–465. pmid:17328122
- 35. Belova SE, Ravin N V., Pankratov TA, Rakitin AL, Ivanova AA, Beletsky A V., et al. Hydrolytic capabilities as a key to environmental success: chitinolytic and cellulolytic Acidobacteria from acidic sub-arctic soils and boreal peatlands. Front Microbiol. 2018;9: 1–14.
- 36. Crowe MA, Power JF, Morgan XC, Dunfield PF, Lagutin K, Rijpstra WIC, et al. Pyrinomonas methylaliphatogenes gen. nov., sp. nov., a novel group 4 thermophilic member of the phylum Acidobacteria from geothermal soils. Int J Syst Evol Microbiol. 2014;64: 220–227. pmid:24048862
- 37. Tank M, Bryant DA. Chloracidobacterium thermophilum gen. nov., sp. nov.: an anoxygenic microaerophilic chlorophotoheterotrophic acidobacterium. Int J Syst Evol Microbiol. 2015;65: 1426–1430. pmid:25667398
- 38. Foesel BU, Rohde M, Overmann J. Blastocatella fastidiosa gen. nov., sp. nov., isolated from semiarid savanna soil—the first described species of Acidobacteria Subdivision 4. Syst Appl Microbiol. 2013;36: 82–89. pmid:23266188
- 39. Wüst PK, Foesel BU, Geppert A, Huber KJ, Luckner M, Wanner G, et al. Brevitalea aridisoli, B. deliciosa and Arenimicrobium luteum, three novel species of Acidobacteria subdivision 4 (class Blastocatellia) isolated from savanna soil and description of the novel family Pyrinomonadaceae. Int J Syst Evol Microbiol. 2016;66: 3355–3366. pmid:27255677
- 40. Pankratov TA, Dedysh SN. Granulicella paludicola gen. nov., sp. nov., Granulicella pectinivorans sp. nov., Granulicella aggregans sp. nov. and Granulicella rosea sp. nov., acidophilic, polymer-degrading acidobacteria from Sphagnum peat bogs. Int J Syst Evol Microbiol. 2010;60: 2951–9. pmid:20118293
- 41. Männistö MK, Rawat S, Starovoytov V, Häggblom MM. Granulicella arctica sp. nov., granulicella mallensis sp. nov., granulicella tundricola sp. nov. and granulicella sapmiensis sp. nov., novel acidobacteria from tundra soil. Int J Syst Evol Microbiol. 2012;62: 2097–2106. pmid:22058325
- 42. Oshkin IY, Kulichevskaya IS, Rijpstra WIC, Sinninghe Damsté JS, Rakitin AL, Ravin N V., et al. Granulicella sibirica sp. nov., a psychrotolerant acidobacterium isolated from an organic soil layer in forested tundra, West Siberia. Int J Syst Evol Microbiol. Microbiology Society; 2019;69: 1195–1201. pmid:30775959
- 43. Rawat SR, Männistö MK, Bromberg Y, Häggblom MM. Comparative genomic and physiological analysis provides insights into the role of Acidobacteria in organic carbon utilization in Arctic tundra soils. FEMS Microbiol Ecol. 2012;82: 341–355. pmid:22486608
- 44. Rawat SR, Männistö MK, Starovoytov V, Goodwin L, Nolan M, Hauser LJ, et al. Complete genome sequence of Granulicella mallensis type strain MP5ACTX8T, an acidobacterium from tundra soil. Stand Genomic Sci. 2013;9: 71–82. pmid:24501646
- 45. Rawat SR, Männistö MK, Starovoytov V, Goodwin L, Nolan M, Hauser L, et al. Complete genome sequence of Granulicella tundricola type strain MP5ACTX9T, an Acidobacteria from tundra soil. Stand Genomic Sci. 2014;9: 449–461. pmid:25197431