Building a Cell House from Cellulose: The Case of the Soil Acidobacterium Acidisarcina polymorpha SBC82T

Acidisarcina polymorpha SBC82T is a recently described representative of the phylum Acidobacteriota from lichen-covered tundra soil. Cells of this bacterium occur within unusual saccular chambers, with the chamber envelope formed by tightly packed fibrils. These extracellular structures were most pronounced in old cultures of strain SBC82T and were organized in cluster-like aggregates. The latter were efficiently destroyed by incubating cell suspensions with cellulase, thus suggesting that they were composed of cellulose. The diffraction pattern obtained for 45-day-old cultures of strain SBC82T by using small angle X-ray scattering was similar to those reported earlier for mature wood samples. The genome analysis revealed the presence of a cellulose biosynthesis locus bcs. Cellulose synthase key subunits A and B were encoded by the bcsAB gene whose close homologs are found in genomes of many members of the order Acidobacteriales. More distant homologs of the acidobacterial bcsAB occurred in representatives of the Proteobacteria. A unique feature of bcs locus in strain SBC82T was the non-orthologous displacement of the bcsZ gene, which encodes the GH8 family glycosidase with a GH5 family gene. Presumably, these cellulose-made extracellular structures produced by A. polymorpha have a protective function and ensure the survival of this acidobacterium in habitats with harsh environmental conditions.


Introduction
The Acidobacteriota (former name Acidobacteria) is one of the phylogenetically diverse and cosmopolitan phyla of the domain Bacteria [1][2][3][4][5]. Known diversity within this phylum includes 26 major 16S rRNA gene sequence clades or subdivisions (SDs) [2], which were assigned to 15 class-level groups, five of which contain taxonomically described representatives [6]. Members of the phylum Acidobacteriota colonize a large variety of terrestrial habitats with highly diverse environmental conditions, such as pristine and agricultural soils, wetlands of different trophic status, caves, unfixed aeolian sands, freshwater and subsurface sediments, hot springs, and many other ecosystems [2,5,[7][8][9][10][11][12][13]. Soils of all climatic zones, from sandy soils of arctic tundra to organic-rich soils of tropical forests, are one of the major habitats of these bacteria. As evidenced by the results of cultivation-independent microbial diversity surveys, Acidobacteriota-affiliated sequences comprise between 5 and 50% of 16S rRNA gene reads commonly retrieved from various soils [7,8,10,[13][14][15][16][17]. Notably, soils colonized by the acidobacteria differ dramatically regarding their physicochemical encoding the inner-membrane catalytic subunit A of cellulose synthase (EC 2.4.1.12), which contains the cytosolic domain of the family 2 of glycosyltransferases (GT2). It is generally accepted that this gene is the first one in the bcsABCD operon [35]. The next bcsB gene encodes the periplasmic B subunit with a carbohydrate-binding domain. It is indispensable for cellulose polymerization, thus making BcsB a co-catalytic subunit. A number of Proteobacteria synthesize the two-domain BcsAB protein. The tetratricopeptide repeat-rich BcsC protein is located in the outer membrane and acts as a porin. The BcsD protein is facultative and contributes to the regular packing of glucan chains in species secreting crystalline cellulose. The same locus hosts the bcsZ gene encoding periplasmic cellulase, which belongs to the GH8 family of glycoside hydrolases. Additional Bcs components, such as BcsH, BcsO, BcsP, BcsQ, BcsR, or BcsS can also be encoded. The neighboring bcsEFG operon is also involved in cellulose synthesis. The structure and composition of the bcs locus vary greatly in different Proteobacteria [35][36][37][38][39]. Cellulose synthesis loci have also been found in Acidobacteriota [23], Actinobacteria [40], Cyanobacteria [36], and Firmicutes [41]. Clostridioides difficile cellulose synthesis locus has ccsZ gene encoding the GH5_25 subfamily endoglucanase (GenPept, AJP12293.1).
In this study, we followed up the formation of the unusual extracellular saccular chambers produced by Acidisarcina polymorpha SBC82 T to elucidate their nature and the potential functional role. We also performed a detailed analysis of the cellulose synthase locus in the genome of strain SBC82 T and compared it to that in other described acidobacteria.

Strain and Culture Conditions
Acidisarcina polymorpha SBC82 T (=KCTC 82304 T = VKM B-3225 T ) was used in this study. Strain SBC82 T was maintained on the solid medium MA containing (per liter distilled water): 0.5 g fructose, 0.1 g yeast extract, 0.2 g KNO 3 , 0.04 g MgSO 4 ·7H 2 O, 0.04 g KH 2 PO 4 , 0.02 g CaCl 2 ·2H 2 O, 9 g phytagel, 0.03 g alginic acid, pH 4.6-5.0. For the electron microscopy, X-ray microanalysis, and the analysis using small angle X-ray scattering, strain SBC82 T was grown in the liquid medium MA with shaking at 180 rpm under room temperature.

Examination of Cell Morphology
Morphological observations of 10-, 20-and 40-day-old cells of strain SBC82 T were made with a Zeiss Axioplan 2 microscope and Axiovision 4.2 software (Zeiss, Germany).
For preparation of ultrathin sections, 14-day-old and 30-day-old cells of strain SBC82 T were collected from plates and pre-fixed with 1.5% (w/v) glutaraldehyde in 0.05 M cacodylate buffer (pH 6.5) for 1 h at 4 • C and then fixed with 1% (w/v) OsO 4 in the same buffer for 4 h at 20 • C. After dehydration in an ethanol series, the samples were embedded into Epon 812 epoxy resin. Sections were obtained using the Reichert-Jung Ultracut ultramicrotome (Austria). Sections were mounted on the support grids, contrasted for 30 min in a 3% solution of uranyl acetate in 70% alcohol, additionally contrasted by lead citrate at 20 • C for 4-5 min, and examined in the JEM-1400 transmission electron microscope (JEOL, Japan) at accelerating voltage of 80 kV at the UNIQEM Collection Core Facility, Research Center of Biotechnology of the Russian Academy of Science.

X-ray Microanalysis
X-ray microanalysis of unstained ultrathin cell sections was performed using a JEM-1400 microscope (JEOL, Akishima, Japan) equipped with an X-ray microanalyzer INCA Energy TEM 350 EDS system (Oxford Instruments) at the UNIQEM Collection Core Facility, Research Center of Biotechnology of the Russian Academy of Science.

Small Angle X-ray Scattering (SAXS) and Data Processing
In the study, 14-and 45-day-old cultures of strain SBC82 T were used for small-angle X-ray scattering detection. Similar culture aliquots were centrifuged at 10,000× g for 5 min. A thin layer of cellular biomass was placed in the special mesh sample holder, which was installed opposite the X-ray beam.
SAXS measurements were performed at synchrotron stations ID23-1 and ID23-2 at the European Synchrotron Radiation Facility-ESRF (Grenoble, France), with the following parameters: Wavelength-2 Å, beam spot-4 microns, flux-5 × 10 10 photons/second, exposition time-5 s per shot. Scattered radiation was recorded by Pilatus 6M and Pilatus 3X 2M detectors, placed at a 982.8 mm distance behind the sample holder. All measurements were held at 100 K temperature. The scheme of the SAXS experiment is depicted in Figure 1A.
Energy TEM 350 EDS system (Oxford Instruments) at the UNIQEM Collection Core Facility, Research Center of Biotechnology of the Russian Academy of Science.

Small Angle X-ray Scattering (SAXS) and Data Processing
In the study, 14-and 45-day-old cultures of strain SBC82 T were used for small-angle X-ray scattering detection. Similar culture aliquots were centrifuged at 10,000× g for 5 min. A thin layer of cellular biomass was placed in the special mesh sample holder, which was installed opposite the X-ray beam.
SAXS measurements were performed at synchrotron stations ID23-1 and ID23-2 at the European Synchrotron Radiation Facility-ESRF (Grenoble, France), with the following parameters: Wavelength-2 Å, beam spot-4 microns, flux-5 × 10 10 photons/second, exposition time-5 s per shot. Scattered radiation was recorded by Pilatus 6M and Pilatus 3X 2M detectors, placed at a 982.8 mm distance behind the sample holder. All measurements were held at 100 K temperature. The scheme of the SAXS experiment is depicted in Figure 1A. For a detailed analysis of the angle distribution of individual fibrils respectively to each other ( ), values of radial averaged intensity profile ( ) were calculated by using Equation (1). The relationship between ( ) and ( ) is shown in (2) and Figure   1B:

Enzymatic Treatment of Cells of Strain SBC82 T with Cellulase
Cellulase (other names: endo-1,4-β-glucanase, 1,4-(1,3:1,4)-β-D-glucan-4-glucanohydrolase) (Serva Feinbiochemica, Heidelberg, Germany) was used to treat the cellular aggregates of strain SBC82 T . For this, 6 mg cellulase (1.55 U/mg) was dissolved in 300 μL of 20 mM acetic buffer (pH 4.5). The resulting cellulase solution (0.03 U/μL) was used to treat the cells collected from 45-day-old culture of strain SBC82 T at 40 °C for 1-2 h. The results of this treatment were verified by using phase-contrast microscopy. For a detailed analysis of the angle distribution of individual fibrils respectively to each other f (µ), values of radial averaged intensity profile H(χ) were calculated by using Equation (1). The relationship between H(χ) and f (µ) is shown in (2) and Figure 1B:

Phylogenetic Analysis
The full-length amino acid sequence of the BcsAB protein from Acidisarcina polymorpha SBC82 T (ACPOL_3662, AXC12945.1) was used to search for close homologs using the blastp algorithm on the NCBI website (http://www.ncbi.nlm.nih.gov/ (accessed on 9 April 2022). Proteins of acidobacteria Candidatus Sulfotelmatomonas gaucii (GenPept, SPE19250.1) and Granulicella aggregans M8UP14 (MBB5059869.1) were excluded from further analysis, as they are incomplete and, presumably, correspond to pseudogenes. For the construction of a multiple alignment and further phylogenetic analysis, 74 nearest (almost) full-length homologs of ACPOL_3662 were selected. Multiple sequence alignment was prepared manually using the program BioEdit (https://bioedit.software.informer.com/7.2 (accessed on 9 April 2022) on the basis of blastp pairwise alignments.
After removing the most variable regions, the multiple sequence alignment was used to implement phylogenetic inference programs, using the maximum parsimony method. Programs SEQBOOT, PROTPARS, and CONSENSE from the PHYLIP package (http:// evolution.gs.washington.edu/phylip.html (accessed on 9 April 2022) were successively used to derive confidence limits estimated by 1000 bootstrap replicates. The program TreeView Win32 (https://colab.research.google.com/drive/11lKiuyH6X2o2XifQKxiH5 JjLAskO_3dX (accessed on 9 April 2022) was used for drawing the trees, while the most divergent group among the analyzed proteins (which corresponded to proteobacteria) was chosen as the outer one. The list of experimentally characterized representatives of the GT2 glycosyltransferase family was taken from the CAZy database [34].

Cell Morphology
Examination of cells of strain SBC82 T at different growth stages by phase-contrast microscopy revealed highly variable cell morphology that varied from single polymorphic cells to large aggregates composed of sarcina-like cell clusters ( Figure 2). At the early growth stages, the cells were odd in shape, elongated, and slightly convoluted ( pha SBC82 T (ACPOL_3662, AXC12945.1) was used to search for close homologs using the blastp algorithm on the NCBI website (http://www.ncbi.nlm.nih.gov/ (accessed on 9 April 2022). Proteins of acidobacteria Candidatus Sulfotelmatomonas gaucii (GenPept, SPE19250.1) and Granulicella aggregans M8UP14 (MBB5059869.1) were excluded from further analysis, as they are incomplete and, presumably, correspond to pseudogenes. For the construction of a multiple alignment and further phylogenetic analysis, 74 nearest (almost) full-length homologs of ACPOL_3662 were selected. Multiple sequence alignment was prepared manually using the program BioEdit (https://bioedit.software.informer.com/7.2 (accessed on 9 April 2022) on the basis of blastp pairwise alignments.
After removing the most variable regions, the multiple sequence alignment was used to implement phylogenetic inference programs, using the maximum parsimony method. Programs SEQBOOT, PROTPARS, and CONSENSE from the PHYLIP package (http://evolution.gs.washington.edu/phylip.html (accessed on 9 April 2022) were successively used to derive confidence limits estimated by 1000 bootstrap replicates. The program TreeView Win32 (https://colab.research.google.com/drive/11lKiuyH6X2o2Xi-fQKxiH5JjLAskO_3dX (accessed on 9 April 2022) was used for drawing the trees, while the most divergent group among the analyzed proteins (which corresponded to proteobacteria) was chosen as the outer one. The list of experimentally characterized representatives of the GT2 glycosyltransferase family was taken from the CAZy database [34].

Cell Morphology
Examination of cells of strain SBC82 T at different growth stages by phase-contrast microscopy revealed highly variable cell morphology that varied from single polymorphic cells to large aggregates composed of sarcina-like cell clusters ( Figure 2). At the early growth stages, the cells were odd in shape, elongated, and slightly convoluted (

Cell Ultrastructure
Electron microscopy of ultrathin cell sections of strain SBC82 T revealed a cell wall structure typical of Gram-negative bacteria. The cell wall consisted of a cytoplasmic membrane, a thin electron-dense peptidoglycan layer, and an outer membrane ( Figure 3). Many small heterogeneous electron-dense inclusions that looked similar to polyphosphate granules were observed in the cytoplasm.

Cell Ultrastructure
Electron microscopy of ultrathin cell sections of strain SBC82 T revealed a cell wall structure typical of Gram-negative bacteria. The cell wall consisted of a cytoplasmic membrane, a thin electron-dense peptidoglycan layer, and an outer membrane ( Figure 3). Many small heterogeneous electron-dense inclusions that looked similar to polyphosphate granules were observed in the cytoplasm.  The unique ultrastructural feature of this acidobacterium, however, was the presence of extracellular structures appearing as saccular chambers, with the cells occurring and dividing inside these structures ( Figure 3). The chamber envelope was composed of tightly packed fibrils layered in parallel to each other (see insert in Figure 3b). These characteristic structures were observed in cells from both young ( Figure 3a) and old cultures of strain SBC82 T (Figure 3b), but the chamber envelope was significantly thicker in cells from old cultures. The cells occurring inside the chambers were separated from the chamber envelope by the extensive electron-transparent matrix. The latter contained a fibrillar matter of low density as well as numerous vesicles, which were released from the outer cell membrane in the chamber matrix and transferred to the inner surface of the chamber envelope. Notably, the chamber division by means of the fibrillar septum formation was also observed (Figure 3b).
X-ray microanalysis of the ultrathin cell sections revealed the presence of silicon (Si) in the chamber wall (Figure 3c).

X-ray Diffraction Patterns from Cells of Strain SBC82 T
SAXS measurements have an advantage over electron microscopy due to the ability to analyze samples in the native hydrated state, which minimizes the risk of artifacts and allows quick data acquisition. This method has been used to detect the orientation distribution of cellulose microfibrils in natural objects [42]. In this work, the SAXS method was applied to analyze bacterial cell conglomerates in a 45-day-old culture of strain SBC82 T (Figure 4(Ia)). A 14-day-old culture consisting mainly of single cells was used as a control (Figure 4(IIa)). The diffraction pattern of the detector-recorded two-dimensional distribution of scattering intensity from a 45-day-old culture of strain SBC82 T was represented by symmetrical radial streaks (Figure 4(Ib)). This indicated the presence of a rod-like scattering object resembling ordered fibrils in the sample. In contrast, the scattering intensity recorded for a 14-day-old culture had circular symmetry, indicating the absence of ordered fibrils which, however, did not exclude the presence of a small amount of stochastically oriented fibrils (Figure 4(IIb)).
Further processing of the obtained diffraction data was performed to elucidate the structure of the chamber envelope produced by strain SBC82 T . Based on the two-dimensional distributions of scattering intensity of 45-and 14-day-old cell samples, one-dimensional curves of averaged intensity depending on the distance from the detector center I(q r ) were plotted (Figure 4(Ic,IIc)). For the 45-day-old culture, an averaging was performed only over the detector region where symmetrical streaks were recorded. The obtained curve was approximated (in the region of q values corresponding to intraparticle interference) by a power law decay~q −2.03 and contained no modulations and peaks (Figure 4(Ic)). In the SAXS experiment, this decay of the scattering intensity indicates the presence of lamellar (layered) structures in the solution [43]. In our case, this indicated the layered structure of the acidobacterial chamber envelope, with each layer consisting of densely packed fibrils. Averaging of the scattering intensity for the 14-day-old culture was performed over the entire area (360-degree angle), since the two-dimensional scattering pattern possessed a circular symmetry. The one-dimensional curve of the averaged intensity dependence on the distance from the detector center clearly differed from that of 45-day-old culture ( Figure 4(IIc)). The curve could not be approximated by a power decay function but had a very broad peak centered at the q = 0.85 nm −1 value. This suggested the presence of some weakly periodic structure in the examined material with a characteristic period of 7.5 nm.   To characterize the arrangement of fibrous structures in a 45-day-old culture, the obtained diffraction data were processed in a different way via averaging of two-dimensional distributions of scattering intensity by each detector pixel corresponding to the same angle χ (the angle is indicated in Figure 1) (H(χ)). The curve displayed in Figure 4(Id) shows the numerical expression of the sum of scattering intensity recorded by each detector pixel located along the angle χ. Determined H(χ) values give the possibility to calculate the positional arrangements of individual microfibrils in the chamber envelope of strain SBC82 T .
The H(χ) plot (Figure 4(Id)) shows the intensity peak corresponding to the position of the symmetric streak in the two-dimensional diffraction pattern from a 45-day-old culture. Notably, the central part of the streak was much brighter than its surroundings. This suggests that the microfibrils in the chamber envelope deviate from each other by a small angle. The H(χ) values were used to calculate the distribution of the angle µ between two individual microfibrils (see Formula (2) in the Section 2.4.). The corresponding probability distribution curve f (µ) is shown in Figure 4(Ie). According to these calculations, the microfibrils from the chamber envelope were arranged nearly parallel to each other, because the root-mean-square deviation of µ was only 3.3 • .
Analysis of the averaged intensity H(χ) curve for a 14-day-old culture, which had no peaks (Figure 4(IId)), confirmed the conclusion about the absence of fibrous structures or their stochastic arrangement relative to each other in young cells.
Thus, the obtained SAXS data are in good agreement with the observations made by using electron microscopy, which showed the occurrence of densely packed fibrous structures in the cell chambers of old cultures of strain SBC82 T (Figure 3). Although these fibrous structures were also present in minor amounts in cells from young cultures, they were not yet structurally organized and densely packed.

Enzymatic Treatment of Cell Conglomerates with Cellulase
The treatment of 45-day-old culture of strain SBC82 T composed of tightly packed cell conglomerates with cellulase resulted in an apparent degradation of cell clusters and related extracellular structures so that the obtained suspension was composed of single cells and small incoherent cell aggregates ( Figure 5). This result provided additional evidence for the presence of cellulose in the extracellular saccular structures of strain SBC82 T .

Structure of the Cellulose Synthesis Locus
According to Table S1 from the supplement to Belova et al. (2018) [21], this locus contains five genes in Acidisarcina polymorpha SBC82 T : ACPOL_3660-ACPOL_3664. The ACPOL_3660 encodes GH5 family cellulase, ACPOL_3661 and ACPOL_3662 encode BcsQ and BcsAB proteins, respectively, and ACPOL_3663 and ACPOL_3664 encode two fragments of the BcsC protein ( Figure 6).

Structure of the Cellulose Synthesis Locus
According to Table S1 from the supplement to Belova et al. (2018) [21], this locus contains five genes in Acidisarcina polymorpha SBC82 T : ACPOL_3660-ACPOL_3664. The ACPOL_3660 encodes GH5 family cellulase, ACPOL_3661 and ACPOL_3662 encode BcsQ and BcsAB proteins, respectively, and ACPOL_3663 and ACPOL_3664 encode two fragments of the BcsC protein ( Figure 6).

Figure 5.
Phase-contrast images of cells aggregates of Acidisarcina polymorpha SBC82 T after cellulase treatment (b) compared to intact aggregates (a). Clear zones that appeared due to the degradation of cell conglomerates are indicated by black arrows. Bar, 5 μm.

Structure of the Cellulose Synthesis Locus
According to Table S1 from the supplement to Belova et al. (2018) [21], this locus contains five genes in Acidisarcina polymorpha SBC82 T : ACPOL_3660-ACPOL_3664. The ACPOL_3660 encodes GH5 family cellulase, ACPOL_3661 and ACPOL_3662 encode BcsQ and BcsAB proteins, respectively, and ACPOL_3663 and ACPOL_3664 encode two fragments of the BcsC protein ( Figure 6).  Possible functionality of the BcsC protein is unclear: the corresponding genome region is considered a pseudogene by Belova et al. [21] and GenBank (CP030840.1). The structure of the bcs locus of A. polymorpha is similar to that described previously for several other representatives of acidobacteria [23], except that instead of the GH5 family glycoside hydrolase gene, other acidobacteria posses the BscZ endoglucanase (GH8 family) gene located between the bscAB and bscC genes ( Figure 6).
Our screening of the NCBI database using the ACPOL_3662 protein as a query revealed that its close homologs are present in many members of the order Acidobacteriales, with some of them having two paralogs. Analysis of the genome context showed that the bcsAB gene is usually located in acidobacteria as part of the cellulose synthesis locus (Figure 6), but if there are two paralogs, then this is true for only one of them. The figure shows the structure of this locus for three such organisms, each having two bcsAB paralogs. The second paralogous gene is located in Acidobacteria bacterium AB60 (DYQ86_07050), Granulicella rosea DSM 18704 (SAMN05421770_10970), and Terriglobus roseus DSM 18391 (Terro_0419) in a variable genome context. In the case of T. roseus, the neighboring gene (Terro_0418) turned out to be the gene of yet another glycosyltransferase (GT4 family). Our analysis confirmed the data of Rawat et al. [23] about the high conservation of the structure of the cellulose synthesis locus among acidobacteria. The statement of Rawat et al. (2012) about the absence of the bcsC gene in Acidobacterium capsulatum ATCC 51196 appears to be incorrect since it was annotated as a pseudogene (ACP_0077) in the corresponding genome sequence (GenBank, CP001472.1).
Among the host organisms of more distant homologs of the ACPOL_3662 protein, it was possible to detect representatives of a number of proteobacteria that possess the two-domain BcsAB protein. In the case of Escherichia coli GE3 (GenBank, CP012376.1), the structure of the cellulose synthesis locus characteristic of acidobacteria was observed (AKN41_3846-AKN41_3849), while in Shewanella avicenniae FJAT-51800 (CP071503.1) there was a slightly different order and set of genes ( Figure 6). In both cases, the presence of the gene encoding BcsZ endoglucanase from the GH8 family was unchanged. This indicates that it was the original component of this locus in acidobacteria but was lost in the lineage leading to A. polymorpha SBC82 T . Its ACPOL_3660 protein turned out to belong to the GH5_1 subfamily, which is not characteristic of most acidobacteria: Among genome-wide sequences, its fairly close homolog (56% amino acid sequence identity) was found in Edaphobacter aggregans EB153 (EDE15_2663), but its gene located in the genome far from the cellulose synthesis locus (EDE15_1043-EDE15_1046).

Phylogenetic Analysis of Cellulose Synthase BcsAB
Multiple alignment of bacterial cellulose synthase amino acid sequences showed that they consist of two conservative regions (BcsA and BcsB), which are connected by a linker of variable size and composition. These two regions were used to reconstruct BcsAB phylogenetic tree (Figure 7). All analyzed acidobacterial proteins compose a very stable monophyletic group (95.7% bootstrap support), whose most divergent representative is ACPOL_3662 from A. polymorpha SBC82 T . Paralogous BcsAB proteins encoded in eight acidobacteria genomes belong to two distinct clusters (99.8% and 99.2%) on the tree. One of them includes proteins encoded by genes located in cellulose synthase loci, but the other one contains only their outcast paralogs.

Discussion
The first evidence for exopolysaccharide biosynthesis by Acidobacteriota was obtained by Eichorst and co-authors [44], who isolated several representatives of these bacteria from soils and described them as representing a novel genus and species, Terriglobus roseus. These soil isolates produced an extracellular matrix of yet unknown chemical composition, which caused cells to stick together and form visible clumps in liquid culture. As was suggested in that study, the formation of such extracellular material in soils may serve as a form of protection from predation or as a web to trap water or nutrients and may be involved in the formation of soil aggregates [44]. Later, extensive production of an amorphous extracellular polysaccharide-like substance was described for representatives of another acidobacterial genus, Granulicella, which was isolated from boreal peatlands and Arctic tundra soils [16,45]. Similar observations regarding the exopolysaccharide biosynthesis were reported for many other members of the order Acidobacteriales, which corresponds to the formerly recognized subdivision 1 (SD 1) of the phylum Acidobacteriota [6]. In addition to Granulicella and Terriglobus, the list of exopolysaccharide-producing members of the Acidobacteriales includes the genera Acidobacterium, Acidipila, Acidisoma, Bryocella, Edaphobacter, and Terracidiphilus. The only study that characterized exopolysaccharides produced by acidobacteria was performed for two strains of the genus Granulicella [46]. However, only a fraction of water-soluble polysaccharides was examined in that work.
The first genome analysis of acidobacteria was performed for Acidobacterium capsulatum DSM 11244 T , Candidatus Koribacter versatilis Ellin345, and Candidatus Solibacter usitatus Ellin6076 [22]. Of these, only Acidobacterium capsulatum DSM 11244 T belongs to the order Acidobacteriales. The A. capsulatum genome contained an operon with all the genes identified as necessary for cellulose synthesis and, at that time, was one of the very few non-proteobacterial genomes reported to contain this operon. The latter appeared to be incomplete in two other examined acidobacteria since the cellulose synthesis gene was not identified. Later, the genome analysis of three Arctic tundra soil acidobacteria of the genera Granulicella and Terriglobus also revealed their ability for biosynthesis of diverse structural polysaccharides, including cellulose [23]. The ability to synthesize structural carbohydrates involved in cell envelope biogenesis was suggested to be of special importance for maintaining cell integrity in tundra soil habitats that are characterized by low temperatures and multiple freeze-thaw cycles [23]. As we see from Figure 7, the ability to produce cellulose is encoded in the genomes of many representatives of the order Acidobacteriales, including both strains from tundra habitats and other soils. The presence of cellulose synthesis genes suggests potential traits for desiccation resistance, biofilm formation, and/or contribution to soil structure [22].
The experimental evidence for cellulose biosynthesis by acidobacteria, however, was lacking till now. In our study, the SAXS method using the synchrotron radiation source ESRF in Grenoble was applied for the first time to analyze the structure of the capsular envelope produced by the soil acidobacterium. Previously, this method was used to study the structure of wood from various plant species [42], cellulose-containing cell walls of Arabidopsis thaliana and the algae Chara corallina, and cells of Dictyostelium discoideum, which contain actin protein fibrils [47]. Using the SAXS method allows obtaining averaged data for an extremely small object determined by the size of the X-ray beam and, in the presence of a fibrous structure, to determine the orientation of microfibrils. Application of this method to 45-day-old cellular conglomerates produced by strain SBC82 T revealed the occurrence of a layered structure consisting of densely packed fibrils, presumably of cellulose nature. Calculations have shown that the packing of fibrils in the cell envelope is very dense, even denser than that in the wood of two softwood species (Norwegian spruce and Scots pine) and two hardwood species (pedunculated oak and copper beech) [42]. In contrast, the results obtained for 14-day-old culture demonstrated a lack of an orderly stacked fibril structure which, however, does not rule out the presence of fibrils in their stochastic orientation. A similar situation was described for algal cells [48] in which dense packing of cellulose fibrils in the cell wall was evident in six-day-old individuals and was not evident in young one-day-old cells. Although the streak pattern obtained in our study was clearly different from that obtained for actin-containing cells of Dictyostelium discoideum [47], these results served only as indirect evidence for cellulose presence in chamber envelopes of strain SBC82 T . The treatment of 45-day-old culture of strain SBC82 T with cellulase and the observed degradation of cell clusters provided additional direct evidence for the presence of cellulose in the extracellular saccular structures of this soil acidobacterium.
Notably, the catalytic subunit BcsAB of cellulose synthase in Acidisarcina polymorpha SBC82 T occupies the most divergent position in the corresponding phylogenetic tree (Figure 7). The chamber-like structures described in this study have not been observed in other members of the Acidobacteriales as well as in other cultivated bacteria. The only study that described somewhat similar extracellular structures was focused on radiation-resistant pseudomonads [49]. In that case, the chamber-like envelopes played a protective role by integrating cells and providing vegetative growth in conditions of continuous intensive UV irradiation. The material used to construct the chamber-like structures in these pseudomonads, however, has not been analyzed.
The analysis of BcsAB proteins included in Figure 7B confirms that all of them belong exclusively to members of the order Acidobacteriales. It is likely that acidobacteria of this order obtained the fused version of bcsAB gene from a common ancestor that received it by lateral transfer from some proteobacteria. The conservative structure of the acidobacterial bcs locus clearly indicates that it was wholly transferred from there. The exact origin of the bcsAB gene remains unclear: most likely, the fusion of the bcsA and bcsB genes in proteobacteria has happened several times independently. Even within the same species, several cases have been found of both the presence of individual bcsA and bcsB genes, and the fused bcsAB gene: Enterobacter cloacae (GenPept, CZX88434.1 and CZW39489.1), Escherichia coli (CQR82949.1 and ASO85436.1), and Salmonella enterica (EBH9976197.1 and SUG18009.1).
In the lineage leading to Acidisarcina polymorpha SBC82 T , a non-orthologous displacement of the glycoside hydrolase gene has occurred. As a result, its genome encodes the ACPOL_3660 protein belonging to the GH5_1 subfamily. It should be noted that Clostridioides difficile (GenPept, AJP12293.1) and Clostridium chromiireducens (OPJ65959.1) have GH5_25 encoding genes in their cellulose synthesis loci (53% amino acid sequence identity). Therefore, we can suppose that GH8/GH5 displacement in the bacterial cellulose synthesis system can provide some adaptive advantage to the cells.
Several acidobacterial species have two paralogs of the fused bcsAB gene. Only one of them is located in the cellulose synthesis locus. The second paralog exists in a variable genome context and corresponds to a distinct cluster on the phylogenetic tree indicating its independent evolutionary history. This group of paralogs has appeared in acidobacteria either by ancient gene duplication or by a lateral gene transfer.
In summary, we confirmed the occurrence of genome-encoded capability for cellulose biosynthesis in Acidisarcina polymorpha SBC82 T and obtained evidence for cellulose presence in the fibrillar envelope produced by this soil acidobacterium. Apparently, these cellulosemade extracellular structures produced by strain SBC82 T have a protective function and ensure the survival of this acidobacterium in habitats with harsh environmental conditions. A more detailed analysis of the material composing the cellular envelope that, most likely, includes not only cellulose but also some other polysaccharides, was not feasible due to the extremely slow and poor growth of strain SBC82 T , which made it impossible to collect sufficient cell material for the conventional chemical analysis. The latter could potentially be applied for the analysis of extracellular material produced by fast-growing acidobacteria of the genera Terriglobus and Granulicella, which possess genome-encoded capability for cellulose biosynthesis.