- Open Access
Glycoside hydrolase gene transcription by Alicyclobacillus acidocaldarius during growth on wheat arabinoxylan and monosaccharides: a proposed xylan hydrolysis mechanism
Biotechnology for Biofuelsvolume 11, Article number: 110 (2018)
Metabolism of carbon bound in wheat arabinoxylan (WAX) polysaccharides by bacteria requires a number of glycoside hydrolases active toward different bonds between sugars and other molecules. Alicyclobacillus acidocaldarius is a Gram-positive thermoacidophilic bacterium capable of growth on a variety of mono-, di-, oligo-, and polysaccharides. Nineteen proposed glycoside hydrolases have been annotated in the A. acidocaldarius Type Strain ATCC27009/DSM 446 genome. Experiments were performed to understand the effect of monosaccharides on gene expression during growth on the polysaccharide, WAX.
Molecular analysis using high-density oligonucleotide microarrays was performed on A. acidocaldarius strain ATCC27009 when growing on WAX. When a culture growing exponentially at the expense of arabinoxylan saccharides was challenged with glucose or xylose, most glycoside hydrolases were downregulated. Interestingly, regulation was more intense when xylose was added to the culture than when glucose was added, showing a clear departure from classical carbon catabolite repression demonstrated by many Gram-positive bacteria. In silico analyses of the regulated glycoside hydrolases, along with the results from the microarray analyses, yielded a potential mechanism for arabinoxylan metabolism by A. acidocaldarius. Glycoside hydrolases expressed by this strain may have broad substrate specificity, and initial hydrolysis is catalyzed by an extracellular xylanase, while subsequent steps are likely performed inside the growing cell.
Glycoside hydrolases, for the most part, appear to be found in clusters, throughout the A. acidocaldarius genome. Not all of the glycoside hydrolase genes found at loci within these clusters were regulated during the experiment, indicating that a specific subset of the 19 glycoside hydrolase genes found in A. acidocaldarius were used during metabolism of WAX. While specific functions of the glycoside hydrolases were not tested as part of the research discussed, many of the glycoside hydrolases found in the A. acidocaldarius Type Strain appear to have a broader substrate range than that represented by the glycoside hydrolase family in which the enzymes were categorized.
Alicyclobacillus acidocaldarius is a Gram-positive, spore-forming, thermophilic acidophile that grows optimally in strictly aerobic conditions at 60 °C and at pH levels between 3 and 4. This bacterium was isolated from a hot spring in Nymph Creek area of Yellowstone National Park . Originally classified as Bacillus acidocaldarius, it was reclassified as A. acidocaldarius based on the prevalence of ω-cyclic fatty acids in the cell wall and an abbreviated helix 6 of the 16S rRNA . A. acidocaldarius has subsequently been isolated from diverse habitats, including water and soil from geothermal sites, submarine hot springs, and orchard soils, and also as a contaminant in heat-processed foods (e.g., fruit juices) [3,4,5,6,7,8]. A. acidocaldarius has demonstrated the ability to gain cellular carbon and energy from a wide variety of 5- and 6-carbon sugars, including l-arabinose, ribose, d-xylose, d-galactose, d-fructose, d-mannose, rhamnose, mannitol, and tagatose; the disaccharides d-turanose, melibiose, cellobiose, lactose, maltose, sucrose, and trehalose, as well as the more complex polysaccharides: cellulose, hemicellulose (xylan), starch, and glycogen .
Lignocellulosic biomass is currently being studied as a feedstock for producing fuels and chemicals [10, 11]. Chemically, lignocellulose is a heterogeneous, three-dimensional matrix made of carbohydrate, primarily hemicellulose (20–30%) and cellulose (40–50%), and lignin. These components are intertwined, providing a structure that is highly resistant to microbial degradation. The hemicellulose fraction accounts for nearly one-third of the renewable carbon found on Earth [12, 13]. Hemicellulose is composed predominantly of xylan, which consists of a conserved backbone of 1,4-linked β-d-xylose residues, decorated with side groups composed of arabinose, glucuronic acid, 4-O-methyl-glucuronic acid, and sometimes galactose .
A thorough understanding of polysaccharide, monosaccharide, and mixed saccharide metabolism by A. acidocaldarius will help determine the utility of this bacterium for lignocellulose depolymerization. In part, this can be accomplished by understanding the effect of transcriptional regulation [i.e., carbon catabolite repression (CCR) as one example] on expression of genes encoding glycoside hydrolases, which represents the first step in metabolism of lignocellulose. Genome analysis of A. acidocaldarius has revealed the presence of all components of Gram-positive CCR, as well as other regulators of genes encoding enzymes for monosaccharide and polysaccharide metabolism. The purpose of this study was to monitor gene transcription during growth on the polysaccharide, WAX, and to monitor the regulatory effects of monosaccharides on transcription of A. acidocaldarius genes related to xylan use (i.e., glycoside hydrolases). Analyses were performed using batch chemostat studies and global transcriptome analysis using high-density oligonucleotide microarray studies. This research represents the first transcriptome analysis of A. acidocaldarius growing on plant polysaccharides, such as WAX.
Glycoside hydrolase inventory of A. acidocaldarius
Alicyclobacillus acidocaldarius represents a source of thermostable glycoside hydrolases for application as industrial catalysts for the hydrolysis of the cellulose and hemicellulose components of lignocellulose. To date, 19 glycoside hydrolase genes have been identified in the A. acidocaldarius genome , and β-galactosidase, α-amylase, cellulase, neopullulanase, exo-pectinase, mannanase, β-glycosidase, and endoglucanase enzymes have been expressed and characterized [6, 9, 16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33]. Enzymes required for efficient depolymerization of lignocellulose, and glycoside hydrolase activities found encoded in the genome of A. acidocaldarius or have been characterized by other research laboratories using recombinant A. acidocaldarius enzymes, are shown (Table 1). Analysis of the A. acidocaldarius genome shows 19 gene loci that are annotated as glycoside hydrolases (EC 3.2.1.-), supporting the observed broad substrate range of A. acidocaldarius related to oligo- and polysaccharides for this bacterium (Table 2). Only three of the loci (Aaci_0048, Aaci_2457, and Aaci_2874) encode proteins that contain a signal peptide, which would indicate extracellular activity. The remaining enzymes are likely located and active in the periplasm or cytoplasm of A. acidocaldarius.
Extent of regulation of each of these genes during the experiment is presented, followed by an in-depth discussion of the bioinformatic analysis to assign function of the enzyme product from these genes. Microarray data discussed in this publication have been deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) and are accessible through GEO Series Accession Number GSE89078.
Regulation of proposed arabinoxylan hydrolysis genes
Gene transcription by A. acidocaldarius grown to mid-exponential phase on WAX was compared to transcription following the addition of glucose or xylose. Following initial exposure of the A. acidocaldarius growing on WAX to the monomer sugars, transcription of 13 of the 19 genes that encode putative glycoside hydrolases were downregulated, regardless of whether glucose or xylose was added (Fig. 1). In general, downregulation of glycoside hydrolase gene transcription was greater when xylose was added to the culture than when glucose was added. This is puzzling because in most Gram-positive bacteria, the opposite has been shown [34, 35]. Two glycoside hydrolase genes were upregulated when glucose was added (Aaci_0912 and Aaci_2891), while these same genes were downregulated when xylose was added. Aaci_0912, which is a proposed amylo-α-1,6-glycosidase, was upregulated slightly over onefold (or a 2 times increase comparing final to initial gene transcription) when glucose was added. Since this locus was downregulated when xylose was added, the gene was being expressed at low levels when A. acidocaldarius was growing on WAX. Similarly, Aaci_2891, which encodes a β-galactosidase, was positively regulated when glucose was added, but downregulated when xylose was added. Again, regulation of Aaci_2891 was approximately onefold, regardless of sugar or direction of regulation. Two other glycoside hydrolase genes (Aaci_2869 and Aaci_2887) behaved in the opposite manner of the two genes just discussed and were upregulated when xylose was added but downregulated with glucose. One of these genes, Aaci_2887, which is proposed to encode a putative α-xylosidase, was upregulated nearly fourfold with xylose, but was downregulated nearly twofold when glucose was added. The extent of downregulation (1.4-fold) of Aaci_2869 when glucose was added was greater than the magnitude of upregulation on xylose (onefold). Transcription of the remaining 15 glycoside hydrolase genes was downregulated following addition of either xylose or glucose. Gene locus Aaci_0048, which encodes a 959-amino acid protein annotated as a α-l-arabinofuranosidase-like protein, was downregulated 5.6-fold when xylose was added and 2.8-fold when glucose was added. Gene transcription from locus Aaci_0060 from A. acidocaldarius was downregulated 2.5-fold when glucose was added and 1.7-fold when xylose was added. The putative protein product is annotated as an α-1,2-glucuronidase, which is most likely involved in hydrolysis of 4-O-methyl-d-glucuronic acid or glucuronic acid, both of which are found as modifications of the xylan backbone of WAX. Transcription of gene loci Aaci_0786, Aaci_0789, and Aaci_0797 was downregulated during this experiment, and regulation was greater when xylose was supplied as the inducing sugar. These genes encode a putative α-galactosidase, a hypothetical protein with a glycoside hydrolase catalytic domain, and a second α-galactosidase, respectively. When xylose was added, transcription of gene locus Aaci_1218, a putative β-galactosidase, was downregulated nearly 14-fold, but less than twofold when glucose was added. This gene locus was the most highly downregulated gene proposed to encode a glycoside hydrolase. A second gene (Aaci_1895) that has been annotated as a β-galactosidase/β-glucosidase in glycoside hydrolase family 1 was downregulated 7.9-fold when xylose was added and 2.3-fold when glucose was added. A gene encoding a family 9 glycoside hydrolase (Aaci_2475) was downregulated equally (~ 2.8-fold) when the exponentially growing A. acidocaldarius was exposed to xylose or glucose. The transcript for a 782-amino acid protein, encoded by genome locus Aaci_2630, was downregulated nearly twofold when glucose was added and greater than fourfold when xylose was added. This gene encodes a glycoside hydrolase family 3 domain protein with β-glucosidase or β-xylosidase activity. While transcription of genes encoding enzymes that would be active toward glucose polymers would not be expected during growth on a xylan substrate, gene loci Aaci_2874, which encodes a protein annotated as an α-amylase, was downregulated nearly threefold when xylose was added to the exponentially growing culture, but not with glucose. A gene annotated to encode an α-arabinofuranosidase (Aaci_2894) was downregulated 3.4-fold when xylose was added to the logarithmically growing A. acidocaldarius culture, but less than twofold when glucose was added. Gene Aaci_2887, annotated as a glycoside hydrolase family 31 enzyme, was upregulated 3.8-fold when xylose was added, but downregulated 1.7-fold when glucose was added.
Regulation of other accessory genes associated with arabinoxylan hydrolysis
When considering arabinoxylan compounds, the xylose backbone and sugar moieties attached to this backbone can also have acetyl and ferulic acid groups chemically bound in the arabinoxylan structure. For optimal liberation of the sugars in arabinoxylan, a number of carbohydrate esterases (CE) or deacetylase enzymes are also required for removal of these groups. Putative CE enzymes attributed to A. acidocaldarius strain DSM 446 as annotated in the CAZy database are listed in Table 3. Typically, CE activity is attributed to families 1–7 and 16, acting on fragments generated by the activity of endo-β-1,4-xylanases . Proteins found in CE family 14, such as LmbE family proteins, have broad substrate specificity and can act as CE [37, 38]. Deacetylation activity by this group of enzymes is primarily associated with cell wall rearrangement and sporulation in bacteria, but the broad substrate specificity may allow these enzymes to act on acetyl moieties found in arabinoxylan or even chitin . Six of the eight esterase-encoding genes were regulated during the experiments (Fig. 2).
Two (Aaci_1372 and Aaci_1443) of the four genes encoding CE family 4 enzymes were upregulated during the experiments. Aaci_1372 was upregulated to a greater extent when xylose was added to the A. acidocaldarius strain DSM 446 culture growing on WAX, while the converse was true for Aaci_1443, where greater upregulation was seen when glucose was added. The protein encoded by Aaci_1372 also contains a signal peptide indicating extracellular localization of this enzyme. The protein encoded by Aaci_1372 is homologous to putative xylanase and chitin deacetylase enzymes found in Desulfosporosinus species and to polysaccharide deacetylase enzymes in other Alicyclobacilli (Fig. 3). The putative CE family 4 esterase encoded by Aaci_1443 is homologous to polysaccharide deacetylases in other Alicyclobacilli, as well as xylanase/chitin deacetylase enzymes found in a variety of Anoxybacillus species. While substrate specificity of the proteins encoded by these two loci has not been determined, upregulation during the experiment and homology to other CE family 4 proteins may indicate that these enzymes are involved in deacetylation of arabinoxylan. One possibility is that these genes were expressed to remove the acetyl groups from small oligosaccharides, allowing further metabolism of the arabinoxylan fractions.
All three genes encoding LmbE family proteins found in the A. acidocaldarius strain DSM 446 genome were regulated during the experiments (see Table 3). Aaci_0881 was upregulated to similar levels when either glucose or xylose was added, while Aaci_0954 was only upregulated when xylose was added. In contrast, Aaci_2070 was downregulated during the experiments, 8.3-fold when xylose was added and threefold when glucose was added. Aaci_0881 and Aaci_0954 show homology to annotated deacetylases in other Alicyclobacilli and a number of Geobacillus species (Fig. 4). Aaci_2070 shows homology to deacetylases found in many hyperthermophiles. While this group of enzymes has been annotated as deacetylases, their specific function in arabinoxylan depolymerization, if any, is not known.
Lignocellulose has three components—a hexose-containing component (cellulose), a primarily pentose-containing component (hemicellulose), and an aromatic-containing component (lignin)—and accounts for a vast amount of carbon in the Earth’s surface environment [14, 40]. Due to the massive amount of this material in the environment, many microorganisms have developed the ability to hydrolyze the various components into oligomer and monomer units that are further metabolized for cellular processes, or fermented to produce alcohols and organic acids . An important first step in metabolism of lignocellulose is the production of glycoside hydrolases that hydrolyze chemical linkages between the sugars, establishing the secondary structure of the polymeric compounds (i.e., cellulose, hemicellulose, and lignin). The primary objective of the current study was to understand gene transcription during growth on xylan (i.e., hemicellulose) by the Gram-positive thermoacidophile, A. acidocaldarius. To accomplish this objective, a physiological, molecular, and in silico approach was carried out to identify the specific glycoside hydrolases active when A. acidocaldarius is grown on hemicellulose provided in the form of WAX.
Molecular analysis was carried out using high-density oligonucleotide microarray studies, comparing gene expression of A. acidocaldarius during logarithmic growth on WAX, to gene expression after a monosaccharide, glucose, or xylose, was spiked into the growth medium. If genes were downregulated when the inducing sugar was added, the assumption was made that the protein product of the gene was being expressed during growth on WAX. Conversely, if genes were upregulated when the monosaccharide was added, the assumption was that the gene product was induced by the sugar. In silico analysis of the glycoside hydrolases expressed by A. acidocaldarius was used to determine how xylan metabolism in A. acidocaldarius is related to hemicellulose metabolism in other Bacteria and Archaea. Finally, the whole information was combined to propose a model for xylan metabolism by A. acidocaldarius (Fig. 5).
Steps in arabinoxylan metabolism hypothesized in this model include an initial extracellular depolymerization, followed by numerous steps inside the cell. The first step in the process is extracellular hydrolysis of β-1,4-bonds between internal xylopyranose residues in the arabinoxylan xylose backbone by the α-arabinofuranosides-like enzyme encoded by Aaci_0048. Studies by Eckert and Schneider have shown xylanase activity by this enzyme . Since the remaining enzymes that were downregulated during the experiment do not contain signal sequences, these enzymes are most likely active in the cytoplasm. Inside the cell, smaller β-d-xylopyranose oligosaccharides, including β-d-xylobiose, would be depolymerized by an endo-acting xylosidase (Aaci_0789 and Aaci_2328) and Aaci_2630 and Aaci_2887 for xylobiose molecules. The upregulation of Aaci_2887 indicates a shift in metabolism to hydrolyze xylobiose residues already inside the cell when the second sugar was added. Removal of α-l-arabinofuranose side groups would then be catalyzed by the α-l-arabinofuranosidase encoded by Aaci_2894. Approximately six enzymes that are annotated for putative activity toward 4-O-methyl glucuronic acid or galactose side groups were regulated during the experiments. While galactose is not shown in the schematic (Fig. 5), galactose is often found in arabinoxylan molecules . Finally, a number of proposed esterase enzymes are thought to be active on acetyl- and feruloyl- (not shown in model) groups on the xylan backbone. Polysaccharide deacetylase (Aaci_1372 and Aaci_1443) and LmbE-like proteins (Aaci_0881, Aaci_0952, and Aaci_2070) were regulated during the experiment. Phylogenetic analysis of enzymes encoded by these genes indicated homology to other xylanase type enzymes (Fig. 2) and polysaccharide deacetylase (Fig. 3), which would be involved in removing acetyl side groups from the arabinoxylan molecule.
Bioinformatic analysis comparing the regulated A. acidocaldarius glycoside hydrolase genes to other bacteria allowed for putative verification of function for these genes in A. acidocaldarius. Comparison of gene locus Aaci_0048, which encodes a 959-amino acid protein annotated as a α-l-arabinofuranosidase-like protein, shows conserved domains from the amino acid sequence are homologous to family 44 glycoside hydrolases, which indicate possible activity toward cellulose and xylan (i.e., endoglucanase activity). The N-terminal domain of the protein also contains a sequence with close homology to carbohydrate binding domains of numerous glycoside hydrolases. Prior to publication of the A. acidocaldarius genome, Eckert and Schneider  expressed and characterized a 959-amino acid thermophilic endoglucanase (CelB; Accession AJ551527.1) from A. acidocaldarius with 100% identity to the gene product from Aaci_0048. The recombinant CelB enzyme demonstrated endo-acting activity toward carboxy methyl cellulose, steam-exploded cellulose and oat spelt xylan, demonstrating a broad substrate range for the enzyme. Sequencing and characterization of the enzyme classified it as an extracellular endoglucanase, with high sequence similarity to GH family 51 of arabinofuranosidases. While this enzyme has been annotated as a α-l-arabinofuranosidase-like protein, no activity was demonstrated using p-nitrophenyl-α-l-arabinofuranoside. These results have been verified in unrelated studies at the Idaho National Laboratory for the native enzyme from A. acidocaldarius and a recombinant enzyme expressed from Pichia pastoris (data not shown). In total, these data indicate that the gene product from Aaci_0048 is a glycoside hydrolase protein that catalyzes hydrolysis of β-1,4-linked xylose subunits in the WAX xylan backbone. It should be noted that the glycoside hydrolase encoded by Aaci_0048 shares 43% homology with a β-1,4-glucanase expressed from Alicyclobacillus sp. A4 .
The Homology Toolkit of the Integrated Microbial Genomes (IMG) Tool at JGI was used to compare the α-l-arabinofuranosidase-like protein (Aaci_0048) from A. acidocaldarius to other genes within IMG. As would be expected, the protein showed the highest percent identity to similarly annotated genes in other species of Alicyclobacillus. The most closely related protein sequence from non-Alicyclobacillus was a cellulose binding family II protein from Ktedonobacter racemifer SOSP1-21. A distance tree generated from sequence information shows that, based on percent identity, Aaci_0048 had the highest homology to α-arabinofuranosidases or cellulose binding family II proteins from numerous species of bacteria from the Order Actinomycetales (Fig. 6). Differences in activity between this enzyme and other α-arabinofuranosidases in the database may be due to the additional 200–400 amino acids found in Aaci_0048 and the other similar proteins found in other strains of A. acidocaldarius. For comparison, the β-1,4-glucanase expressed from Alicyclobacillus sp. A4 still has fewer amino acids than the product of Aaci_0048, which contains an additional ~ 250 amino acids at the C-terminus of the protein sequence. There is very little homology between Aaci_0048 and another annotated α-arabinofuranosidase encoded by Aaci_2894 (discussed below). Evolutionarily, while the exact mechanism of catalysis is not known, an α-arabinofuranosidase could have been transferred from an Actinomycete, followed by fusion with a protein that supplied the endoglucanase and endoxylanase activity.
The putative protein product from genome locus Aaci_0060 is annotated as an α-1,2-glucuronidase, which is most likely involved in hydrolysis of 4-O-methyl-d-glucuronic acid or glucuronic acid, both of which are found as modifications of the xylan backbone of WAX. Multidomain analysis of the protein product of Aaci_0060 supports the annotation as an α-1,2-glucuronidase, and more detailed analysis shows the presence of C-terminal and middle domains (probable catalytic region) found in glycoside hydrolase family 67 proteins . Aaci_0060 and Aaci_0048 may act in concert, since their loci are in the same gene neighborhood of the A. acidocaldarius genome; however, the relationship of transcription of these two genes has not been thoroughly studied. The proximity of these two genes appears to be conserved in a variety of Alicyclobacillus species, as demonstrated by comparing with the genomes of A. acidocaldarius subsp. acidocaldarius (Tc-4-1) and A. hesperidum URH17-3-68 (Fig. 7).
Phylogenetic analysis of the protein encoded by Aaci_0060 indicates relatedness to proteins from a number of other cellulolytic bacteria found in the phylum Firmicutes (Fig. 8). This protein from A. acidocaldarius shows greater than 50% homology to proteins that have been annotated as α-glucuronidases in both aerobic and anaerobic bacteria in the genera Bacillus, Clostridium, Paenibacillus, Geobacillus, and Thermoanaerobacterium.
Regulation of gene loci Aaci_0786, Aaci_0789, and Aaci_0797 was also stimulated during this experiment. The first of these genes, Aaci_0786, has been annotated as a glycoside hydrolase clan GH-D protein, specifically an α-galactosidase. A conserved domain comparison with other glycoside hydrolases in GenBank indicates specific hits to melibiase, with nonspecific hits to glycoside hydrolase family 31 (GH31) enzymes. Enzymes in the GH31 superfamily are present in all three domains of life, and activities include α-glucosidase, α-xylosidase, 6-α-glucosyltransferase, 3-α-isomaltosyltransferase, and α-1,4-glucan lyase. While this enzyme has been annotated as an α-galactosidase, and could be involved in removal of galactose from the xylan backbone of WAX, expression of this enzyme during growth on WAX may also indicate α-xylosidase activity.
Protein sequence analysis and classification (IPR013781) of Aaci_0789 indicates the presence of a catalytic TIM beta/alpha barrel common to many different families of glycoside hydrolases. This protein is greater than 50% identical with hypothetical proteins containing glycoside hydrolase catalytic domains found in other species of Alicyclobacillus. The proposed protein product for Aaci_0789 shows homology to an endo-β-1,4-mannanase that was characterized in Alicyclobacillus species A4 and A. acidocaldarius strains Tc-4-1 and Tc-12-31 [26, 45] (Fig. 9). The glycoside hydrolase encoded by this gene locus also demonstrates up to 50% homology to endo-β-1,4-xylanases from a number of other Firmicutes, as well as hypothetical proteins from Firmicutes known to degrade xylan. In silico analysis of this gene along with expression results indicate a possible broad substrate range for this enzyme that would include mannose- and xylose-containing polysaccharides.
The final gene locus in this group, Aaci_0797, has been annotated as a glycoside hydrolase family 4, α-galactosidase. The protein encoded by Aaci_0797 is 99% identical to a β-galactosidase cloned from A. acidocaldarius . This recombinant enzyme hydrolyzed o-nitrophenyl-β-d-galactopyranoside and o-nitrophenyl-β-d-fucopyranoside, but did not catalyze hydrolysis of o-nitrophenyl-β-d-glucopyranoside, o-nitrophenyl-β-d-xylopyranoside, or o-nitrophenyl-β-d-arabinopyranoside. These results indicate that Aaci_0797 is likely a β-galactosidase active toward galactose moieties on side chains of the WAX. Aaci_0797 also has greater than 70% homology to α-galactosidase from a number of other Gram-positive bacteria, including numerous species of Paenibacillus.
Arrangement of the genes encoding the glycoside hydrolase clan GH-D (Aaci_0786) and the hypothetical protein containing a glycoside hydrolase catalytic domain (Aaci_0789) are nearly identical in A. acidocaldarius subsp. acidocaldarius (DSM 446), A. acidocaldarius subsp. acidocaldarius (Tc-4-1), and A. hesperidum URH17-3-68 (Fig. 10). However, this gene arrangement is not conserved in A. herbarius DSM 12609. DSM 446 appears to have a cassette containing the glycoside hydrolase family 4 α-galactosidase (Aaci_0797), genes encoding ABC-transporter components (Aaci_0794-0796), and a LacI family transcriptional regulator (Aaci_0793) that is not present in the Tc-4-1 genome. Two β-galactosidase genes can be found in the Tc-4-1 genome, but translated protein sequences show low homology to Aaci_0797, indicating that Tc-4-1 may lack α-galactosidase activity.
Annotation within the JGI-IMG/M categorizes Aaci_1218 with a product name of a glycoside hydrolase family 2 TIM barrel. Analysis of the protein sequence showed no signal sequence, indicating activity inside the cell. Gene Ontology (GO:0004565), clusters of orthologous groups (COG3250), and KEGG Orthology (KO:01190) indicate that this protein may be a β-galactosidase/β-glucuronidase. This 1041 amino acid protein shows the highest homology to glycoside hydrolases found in other Bacilli. This enzyme most likely acts on either galactose or 4-O-methyl-d-glucuronic acid side chains found in WAX. These results indicate that A. acidocaldarius was expressing genes related to side chain hydrolysis of the WAX, which were regulated upon addition of xylose. Aaci_1218 appears to be in an apparent operon with genes for an ABC-transport system and a LacI transcriptional regulator, all of which were highly downregulated.
The 453 amino acid A. acidocaldarius protein encoded by Aaci_1895 has greater than 60% homology to β-galactosidase/β-glucosidase enzymes found in other Firmicutes, namely numerous species of Thermoanaerobacter (Fig. 11). Di Lauro et al.  characterized a recombinant β-glycosidase from A. acidocaldarius that showed wide substrate specificity and was able to hydrolyze β-d-gluco-, -galacto-, and fucosides. The protein encoded by Aaci_1895 is 98% identical to the characterized β-glycosidase, and thus is likely the same protein. Like many other glycoside hydrolase encoding genes in the A. acidocaldarius genome, Aaci_1895 is arranged in an operon with an ABC-type transport system; however, this specific ABC-type transporter has been annotated as an oligopeptide/dipeptide transporter. Utilization of oligopeptide/dipeptide ABC-type transporters for carbohydrate transport is not uncommon in other thermophilic lignocellulose-degrading bacteria, including Thermotoga maritima and Caldicellulosiruptor saccharolyticus [46, 47]. This gene arrangement, including the oligopeptide/dipeptide ABC-type transporter genes, is also common in a number of Thermoanaerobacter species. The primary difference between A. acidocaldarius and the Thermoanaerobacter species is the presence of an additional gene, which encodes a family 2 glycoside hydrolase, situated between the β-galactosidase/β-glucosidase and the transporter genes (Fig. 12). The family 2 glycoside hydrolase gene encoded in the Thermoanaerobacter genomes may have a similar function to the family 2 glycoside hydrolase encoded by Aaci_1218, which was also highly downregulated in A. acidocaldarius, as discussed above. Comparing the family 2 glycoside hydrolase from Thermoanaerobacter mathranii to the A. acidocaldarius subsp. acidocaldarius DSM 446 genome yielded the closest homology to Aaci_1218 (E value = 9e−10); other proteins had E values greater than 1.
Domain analysis of the 537 amino acid protein encoded by Aaci_2475 indicates that the enzyme may have endoglucanase activity, but since it appears to be active on WAX, the enzyme may have an expanded substrate range, and act on substrates other than cellulose or cellobiose. Endoglucanases with broad substrate specificity, including activity toward xylans, have been demonstrated in other thermophiles, including Pyrococcus furiosus and T. maritima [48, 49]. The glycoside hydrolase encoded by Aaci_2475 is next to an ABC-type oligopeptide transport system, which was also highly downregulated when xylose was added, as was the case with Aaci_1895 and neighboring genes. Analysis of the gene neighborhood where this family 9 glycoside hydrolase is found indicates that this apparent operon structure is common in numerous Alicyclobacillus species, including A. acidocaldarius subsp. acidocaldarius (DSM 446), A. acidocaldarius subsp. acidocaldarius (Tc-4-1), and A. hesperidum URH17-3-68 (Fig. 13). Xylan-hydrolyzing activity by an endoglucanase is not unprecedented; Hall et al.  characterized an endoglucanase from Clostridium thermocellum with xylanase activity. More recently, a bifunctional xylanase–glucanase from Paenibacillus sp. Strain E18 was characterized and showed activity toward xylan and glucan . Likewise, bifunctional xylanase/endoglucanase enzymes have been generated from the metagenome libraries from a yak and bovine rumen microbial community and showed activity toward a variety of xylans and glucans [52, 53].
Aaci_2630 encodes a glycoside hydrolase family 3 domain protein with β-glucosidase or β-xylosidase activity. Regulation during growth on WAX further supports the possibility of β-xylosidase activity. Homology analysis of Aaci_2630 to other β-xylosidases genes deposited in the National Center for Biotechnology Information (NCBI) database shows greater than 50% similarity to family 3 glycoside hydrolases from other thermophilic, lignocellulose-degrading bacteria, including Thermoanaerobacter, Caldicellulosiruptor and Thermotoga species (Fig. 14). This proposed β-xylosidase also has ~ 60% identity with three enzymes expressed from recombinant DNA showing specificity toward xylo-oligosaccharides for which activity assays have been performed [54,55,56]. Xyl3A from Caldanaerobius polysaccharolyticus demonstrated both β-glucosidase and β-xylosidase activity, but activity toward xylooligosaccharides was the highest . Enzymes that showed homology to Aaci_2630 from two species of Thermoanaerobacter have also been characterized as β-xylosidases. A bifunctional xylosidase–arabinosidase (XarB) from T. ethanolicus JW200 demonstrated activity toward a variety of substrates, but the highest activity/affinity was toward xylo- and arabinopyranoside molecules . Similarly, a xylodextrin-hydrolyzing xylo-β-glucosidase from T. brockii was shown to possess both β-glucosidase and β-xylosidase activity . Homology of Aaci_2630 with these other proteins that have demonstrated β-xylosidase activity and the downregulation in these experiments support classification of the encoded enzyme as a β-xylosidase.
A section of the A. acidocaldarius genome from gene locus Aaci_2868 through Aaci_2894 appears to be a ‘hot spot’ for glycoside hydrolase activity because 6 of the 19 glycoside hydrolase encoding genes in A. acidocaldarius are found in this region (Fig. 15). Annotation in this region indicates genes encoding both cellulose- and hemicellulose-hydrolyzing enzymes. Half of the glycoside hydrolase genes were regulated during the experiment, maintaining the trend of higher regulation when xylose was added to the exponentially growing culture of A. acidocaldarius.
The translated protein from Aaci_2874 shows over 99% homology to an enzyme cloned from A. acidocaldarius that had α-amylase and pullulanase catalytic activity . A thermoacidophilic α-amylase expressed from Alicyclobacillus sp. A4 which showed 65% identity with the protein encoded by Aaci_2874 also showed high activity toward starch compounds, but the substrate range of the enzyme was not determined . Aaci_2874 is located in what may be an operon containing two other glycoside hydrolases—Aaci_2868 and Aaci_2869—neither of which was regulated to the extent of Aaci_2874 when xylose was added. Since there have been no reports of bifunctional α-amylases, the relationship of this enzyme to WAX depolymerization is not known, but regulation of this gene when A. acidocaldarius was grown on WAX indicates that this enzyme has adapted for hydrolysis of xylan polysaccharides.
During metabolism of WAX, α-arabinofuranosidase catalyzes the hydrolysis of terminal α-arabinofuranosidic linkages to xylose residues in the xylan backbone. Gene locus Aaci_2894 encodes a putative α-arabinofuranosidase that is associated with an operon involved in arabinose metabolism, for which all genes were downregulated, regardless of whether xylose or glucose was used as the inducing sugar. The protein product from this gene has greater than 60% homology to α-arabinofuranosidases found in a number of Gram-positive bacteria, including Geobacillus, Paenibacillus, and Thermoanaerobacterium. As with other genes previously discussed, this or similar operons appear to be conserved among Alicyclobacillus species.
Domain analysis for the protein encoded by Aaci_2887 is consistent with α-xylosidase/α-glucosidase activity, which consists of cleaving the terminal carbohydrate moiety from a variety of substrates. Since the gene was upregulated upon the addition of xylose, this 779-amino acid glycoside hydrolase may be involved in xylose metabolism, or even cleaving xylose that is linked to other sugars present in WAX. Comparison to other genes in the NCBI database indicates that Aaci_2887 is 50% identical to α-xylosidase/α-glucosidase enzymes in other Gram-positive lignocellulose utilizing bacteria, including numerous Paenibacillus and Clostridium species (Fig. 16). The protein sequence from Aaci_2887 is 99% identical to an enzyme expressed from recombinant DNA, thermostable α-glucosidase (ABI81478.1) from A. acidocaldarius that was deposited into GenBank, but no information on enzyme activity or specificity has been published. Other glycoside hydrolase family 31 enzymes that display α-xylosidase activity have been found in the Crenarchaeota Sulfolobus solfataricus (i.e., XylS) and the Gram-positive Cellovibrio japonicus [58, 59].
In general, A. acidocaldarius appears to have enzymes necessary for the efficient depolymerization of WAX. Likewise, these enzymes were expressed when A. acidocaldarius strain ATCC27009 was grown on WAX. A number of the enzymes also appear to have a broad substrate range, allowing for hydrolysis of various bonds present in WAX even though the annotated activity of the enzyme would not predict this hydrolysis.
Glycoside hydrolases, for the most part, appear to be found in clusters, throughout the A. acidocaldarius genome. Not all of the glycoside hydrolase genes found at loci within these clusters were regulated during the experiment, indicating that a specific subset of the 19 glycoside hydrolase genes found in A. acidocaldarius were used during metabolism of WAX. While specific functions of the glycoside hydrolases was not tested as part of the research discussed, many of the glycoside hydrolases found in the A. acidocaldarius Type Strain appear to have a broader substrate range than represented by the glycoside hydrolase family in which the enzymes were categorized.
Materials and methods
Alicyclobacillus acidocaldarius ATCC 27009 was purchased from the American Type Culture Collection (ATCC) and used for all experiments. Chemostat inoculum was prepared by inoculating 1 mL of the frozen stock into 25 mL of Modified 402 Medium, which contained the following (g/L): (NH4)2SO4 (1.3), Fe(III) EDTA (0.047), CaCl2·2H2O (0.07), MgSO4·7H2O (0.25), KH2PO4 (3.0), and xylose (4.0). In addition, 1 mL of a mineral (solutions A and B) and vitamin stock (solution C) were added. Solution A contained (g/L): MgCl2 (25), CaCl2·2H2O (6.6), H3BO3 (0.58), FeCl3·6H2O (5), Co(NO3)2·6H2O (0.05), NiCl2·6H2O (0.02). Solution B contained (g/L): MnSO4·H2O (2.0), ZnSO4·7H2O (0.5), CuSO4·5H2O (0.15), Na2MoO4·2H2O (0.025). Solution C contained (g/L): pyridoxine hydrochloride (0.08), folic acid (0.012), thiamine hydrochloride (0.13), riboflavin (0.042), nicotinamide (0.084), p-aminobenzoate (0.088), biotin (0.01), cyanocobalamin (0.0004), d-pantothenic acid, calcium salt (0.086), myo-inositol (0.021), choline bromide (0.053), orotic acid, sodium salt (0.021), and spermidine (0.1). Base medium was autoclaved (121 °C, 20 psi) for 30 min prior to use; KH2PO4 solution was adjusted to pH 4.0, autoclaved separately, and added once the base medium had cooled. The 25 mL A. acidocaldarius culture was grown overnight and then used to inoculate 250 mL of Modified 402 Medium containing 4 g/L xylose. This overnight culture grown at pH 4 at a temperature of 60 °C was then used to inoculate the chemostat.
Experiments were performed in a BioFlo 3000 chemostat system (New Brunswick Scientific, Enfield, CT). Medium was added to the reactor, and the entire reactor was autoclaved at 121 °C, and a pressure of 20 psi for 1 h. WAX (Megazyme, Inc.) was prepared by first wetting 2 g with 95% ethanol. Once in a slurry, 190 mL of sterile distilled water (pH 4.0) was added, and the solution was heated for 30 min to evaporate residual ethanol. Pre-dissolved WAX (2 g/190 mL medium) was then added to the 2 L of Modified 402 Medium in the chemostat, giving a final WAX concentration of 1 g/L. The chemostat was run in batch mode and the temperature was held at 60 °C and the pH of the growth medium was automatically controlled to 4.0 using the addition of 1 N NaOH. To ensure that the cultures were not oxygen limited, dissolved oxygen was controlled to 10%.
Once the A. acidocaldarius culture reached an OD600 of 0.5 while growing on WAX, a sample was taken for RNA extraction, glucose or xylose was added at a concentration of 2 g/L, and then a second sample was immediately (within 5 min) taken for RNA extraction. Three biological replicates for each condition were performed.
Isolation of total RNA
Samples were taken from the chemostat, and RNA was stabilized using RNA Protect Bacteria Reagent (Qiagen, Valencia, CA), followed by flash freezing of the cell pellet in liquid nitrogen and storage at − 80 °C. Total RNA was extracted from the A. acidocaldarius cells using an RNeasy Midi Kit (Qiagen, Valencia, CA) with slight modification of the manufacturer’s protocol. A. acidocaldarius cells were thawed, and lysis was accomplished by adding 200 µL of Tris–EDTA buffer containing 15 mg/mL lysozyme and 0.1 mg/mL proteinase K. Samples were vortexed for 10 s and then incubated at room temperature for 15 min with shaking. Manufacturer’s instructions supplied with the kit were followed up to elution of the RNA. To increase RNA yield, the flow through was re-applied to the spin filter and centrifuged. Residual DNA in the samples was removed by treatment with Ambion TURBO DNA-free kit (Life Technologies, Grand Island, NY). RNA was purified to remove compounds that might interfere with cDNA synthesis and concentrated using ethanol precipitation. To inhibit RNA degradation during storage, 1 µL of Ambion Superase-In RNase Inhibitor (Life Technologies, Grand Island, NY) was added. RNA concentration and purity were determined using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE). RNA integrity was determined using an RNA Nano Chip Kit run on an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA).
Synthesis of cDNA
Double-stranded cDNA was synthesized from total RNA using the Invitrogen Superscript Double-Stranded cDNA Synthesis Kit (Life Technologies, Grand Island, NY) according to manufacturer’s instructions, with the following modifications. Following cDNA synthesis, residual RNA was degraded by adding RNase A, followed by proteins removal using phenol:chloroform:isoamyl alcohol, and separation using Phase Lock Tubes (5 Prime, Inc., Gaithersburg, MD). cDNA in the aqueous phase was then precipitated and concentrated by ethanol precipitation. The pellet was dried and then resuspended in 20 µL of DNase/RNase-Free water and allowed to solubilize overnight. cDNA concentration in each reaction was determined using a NanoDrop ND-1000 spectrophotometer (Thermo Scientific, Wilmington, DE). A DNA 7500 Chip Kit run on an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA) was used to verify that most of the cDNA was ≥ 400 bp.
Microarray experiments and data analysis
Microarrays were designed and synthesized by NimbleGen using their 4 × 72 K Custom Gene Expression Array format from the complete genome sequence information for A. acidocaldarius ATCC 27009 (DSM 446). Seven probes, each 60 nt long, were designed for each of the 3554 identified open reading frames in the genome, based on a genome sequence generated by Integrated Genomics, Inc for the Idaho National Laboratory. Each probe was synthesized on the microarray in triplicate. Control probes were also included to ensure that there was no intra-quadrant contamination during the hybridization process.
One color cDNA labeling using Cy3, hybridization to the A. acidocaldarius microarrays, array imaging, and initial analysis of the array data were performed by NimbleGen. Data were normalized using Nimblescan software, which normalizes probe response using quantile normalization and gene calls generation using Robust Multichip Averaging (RMA) [60, 61]. Log2-transformed RMA data files were imported into ArrayStar 4 software (DNASTAR, Inc., Madison, WI), and the mean expression levels of three replicate arrays for each condition were considered. For comparison of gene expression, statistical significance was determined with a Bonferroni corrected moderate t test, and only genes where expression had greater than 95% confidence (p ≤ 0.05) were considered significant.
Comparative analysis of the A. acidocaldarius genome was performed using the Integrated Microbial Genomes feature within the Joint Genome Institute . Homology determinations for A. acidocaldarius proteins were accomplished using the Basic Local Alignment Search Tool (BLAST) for protein sequences using the ‘blastp’ algorithm [63, 64]. The nonredundant protein sequence database was used for searches, and uncultured and environmental sample sequences were excluded from the search. Phylogenetic trees comparing protein sequences were generated using the Molecular Evolutionary Genetics Analysis (MEGA6) software program . Protein sequences were aligned using the MUSCLE program within MEGA6 [66, 67]. Phylogenetic reconstruction was accomplished using the maximum likelihood statistical method, and distances between sequences were determined using 1000 bootstrap replicates.
National Center for Biotechnology Information
carbon catabolite repression
ATP binding cassette
American Type Culture Collection
Deutsche Sammlung von Mikroorganismen
Kyoto Encyclopedia of Genes and Genomes
Joint Genome Institute
Integrated Microbial Genomes
cluster of orthologous genes
Basic Local Alignment Search Tool
Robust Multichip Averaging
Molecular Evolutionary Genetics Alignment
Multiple sequence comparison by log-expectation
Darland G, Brock TD. Bacillus acidocaldarius sp. nov., an acidophilic thermophilic spore-forming bacterium. J Gen Microbiol. 1971;67(1):9–15.
Wisotzkey J, Jurtshuk P Jr, Fox G, Deinhard G, Poralla K. Comparative sequence analyses on the 16S rRNA (rDNA) of Bacillus acidocaldarius, Bacillus acidoterrestris, and Bacillus cycloheptanicus and proposal for creation of a new genus, Alicyclobacillus gen. nov. Int J System Bacteriol. 1992;42:263–9.
Nicolaus B, Improta R, Manca MC, Lama L, Esposito E, Gambacorta A. Alicyclobacilli from an unexplored geothermal soil in Antarctica: mount Rittmann. Polar Biol. 1998;19(2):133–41.
Goto K, Mochida K, Asahara M, Suzuki M, Yokota A. Application of the hypervariable region of the 16S rDNA sequence as an index for the rapid identification of species in the genus Alicyclobacillus. J Gen Appl Microbiol. 2002;48(5):243–50.
Gouws PA, Gie L, Pretorius A, Dhansay N. Isolation and identification of Alicyclobacillus acidocaldarius by 16S rDNA from mango juice and concentrate. Int J Food Sci Technol. 2005;40(7):789–92.
Di Lauro B, Rossi M, Moracci M. Characterization of a β-glycosidase from the thermoacidophilic bacterium Alicyclobacillus acidocaldarius. Extremophiles. 2006;10(4):301–10.
Groenewald WH, Gouws PA, Witthuhn RC. Isolation and identification of species of Alicyclobacillus from orchard soil in the Western Cape, South Africa. Extremophiles. 2008;12(1):159–63.
Groenewald WH, Gouws PA, Witthuhn RC. Isolation, identification and typification of Alicyclobacillus acidoterrestris and Alicyclobacillus acidocaldarius strains from orchard soil and the fruit processing environment in South Africa. Food Microbiol. 2009;26(1):71–6.
Eckert K, Schneider E. A thermoacidophilic endoglucanase (CelB) from Alicyclobacillus acidocaldarius displays high sequence similarity to arabinofuranosidases belonging to family 51 of glycoside hydrolases. Eur J Biochem. 2003;270(17):3593–602.
Chundawat SP, Beckham GT, Himmel ME, Dale BE. Deconstruction of lignocellulosic biomass to fuels and chemicals. Annu Rev Chem Biomol Eng. 2011;2:121–45.
Sorek N, Yeats TH, Szemenyei H, Youngs H, Somerville CR. The implications of lignocellulosic biomass chemical composition for the production of advanced biofuels. Bioscience. 2014;64:192–201.
Olcay H, Subrahmanyam AV, Xing R, Lajoie J, Dumesic JA, Huber GW. Production of renewable petroleum refinery diesel and jet fuel feedstocks from hemicellulose sugar streams. Energy Environ Sci. 2013;6(1):205–16.
Jae J, Tompsett GA, Lin Y-C, Carlson TR, Shen J, Zhang T, Yang B, Wyman CE, Conner WC, Huber GW. Depolymerization of lignocellulosic biomass to fuel precursors: maximizing carbon efficiency by combining hydrolysis with pyrolysis. Energy Environ Sci. 2010;3(3):358–65.
de Souza WR. Microbial degradation of lignocellulosic biomass. In: Silva AK, Da Silva SS, editors. Sustainable degradation of lignocellulosic biomass - techniques, applications and commercialization. InTech; 2013. p. 207–47.
Mavromatis K, Sikorski J, Lapidus A, Rio TGD, Copeland A, Tice H, Cheng J-F, Lucas S, Chen F, Nolan M, et al. Complete genome sequence of Alicyclobacillus acidocaldarius type strain (104-IAT). Stand Genomic Sci. 2010;2(1):9–18.
La Cara F, Scarffi MR, D’Auria S, Massa R, d’Ambrosio G, Franceschetti G, Rossi M, De Rosa M. Different effects of microwave energy and conventional heat on the activity of a thermophilic ß-galactosidase from Bacillus acidocaldarius. Bioelectromagnetics. 1999;20(3):172–6.
Gul-Guven R, Guven K, Poli A, Nicolaus B. Purification and some properties of a β-galactosidase from the thermoacidophilic Alicyclobacillus acidocaldarius subsp. Rittmannii isolated from Antarctica. Enzyme Microb Technol. 2007;40(6):1570–7.
Di Lauro B, Strazzulli A, Perugino G, La Cara F, Bedini E, Corsaro MM, Rossi M, Moracci M. Isolation and characterization of a new family 42 β-galactosidase from the thermoacidophilic bacterium Alicyclobacillus acidocaldarius: identification of the active site residues. Biochim Biophys Acta (BBA) Proteins Proteomics. 2008;1784(2):292–301.
Yuan T, Yang P, Wang Y, Meng K, Luo H, Zhang W, Wu N, Fan Y, Yao B. Heterologous expression of a gene encoding a thermostable β-galactosidase from Alicyclobacillus acidocaldarius. Biotechnol Lett. 2008;30(2):343–8.
Koivula Y, Hemila H, Pakkanen R, Sibakov M, Palva I. Cloning and sequencing of a gene encoding acidophilic amylase from Bacillus acidocaldarius. J Gen Microbiol. 1993;139(10):2399–407.
Schwermann B, Pfau K, Liliensiek B, Schleyer M, Fischer T, Bakker EP. Purification, properties and structural aspects of a thermoacidophilic α-amylase from Alicyclobacillus acidocaldarius ATCC 27009. Eur J Biochem. 1994;226(3):981–91.
Matzke J, Schwermann B, Bakker EP. Acidostable and acidophilic proteins: the example of the α-amylase from Alicyclobacillus acidocaldarius. Comp Biochem Physiol Part A Physiol. 1997;118(3):475–9.
Morlon-Guyot J, Ordonez RG, Gasparian S, Guyot JP. Pre-harvesting treatments to recover ia a soluble form the cell-bound α-amylase of Alicyclobacillus acidocaldarius grown in liquid culture media containing soluble and granular starch. J Food Sci Technol. 1998;35(2):117–21.
Satheesh Kumar G, Chandra M, Mallaiah K, Sreenivasulu P, Choi Y-L. Purification and characterization of highly thermostable α-amylase from thermophilic Alicyclobacillus acidocaldarius. Biotechnol Bioprocess Eng. 2010;15(3):435–40.
Morana A, Esposito A, Maurelli L, Ruggiero G, Ionata E, Rossi M, Cara FL. A novel thermoacidophilic cellulase from Alicyclobacillus acidocaldarius. Protein Pept Lett. 2008;15:1017–21.
Zhang Y, Ju J, Peng H, Gao F, Zhou C, Zeng Y, Xue Y, Li Y, Henrissat B, Gao GF, et al. Biochemical and structural characterization of the intracellular mannanase AaManA of Alicyclobacillus acidocaldarius reveals a novel glycoside hydrolase family belonging to clan GH-A. J Biol Chem. 2008;283(46):31551–8.
Ordoñez R, Morlon-Guyot J, Gasparian S, Guyot J. Occurrence of a thermoacidophilic cell-bound exo-pectinase in Alicyclobacillus acidocaldarius. Folia Microbiol. 1998;43(6):657–60.
Matzke J, Herrmann A, Schneider E, Bakker EP. Gene cloning, nucleotide sequence and biochemical properties of a cytoplasmic cyclomaltodextrinase (neopullulanase) from Alicyclobacillus acidocaldarius, reclassification of a group of enzymes. FEMS Microbiol Lett. 2000;183(1):55–61.
Eckert K, Zielinski F, Leggio LL, Schneider E. Gene cloning, sequencing, and characterization of a family 9 endoglucanase (CelA) with an unusual pattern of activity from the thermoacidophile Alicyclobacillus acidocaldarius ATCC27009. Appl Microbiol Biotechnol. 2002;60(4):428–36.
Eckert K, Vigouroux A, Lo Leggio L, Moréra S. Crystal structures of A. acidocaldarius endoglucanase Cel9A in complex with cello-oligosaccharides: strong − 1 and − 2 subsites mimic cellobiohydrolase activity. J Mol Biol. 2009;394(1):61–70.
Pereira JH, Sapra R, Volponi JV, Kozina CL, Simmons B, Adams PD. Structure of endoglucanase Cel9A from the thermoacidophilic Alicyclobacillus acidocaldarius. Acta Crystallogr Sect D Biol Crystallogr. 2009;D65(8):744–50.
Thompson DN, Apel WA, Thompson VS, Reed DW, Lacey JA, Henriksen ED. Thermophilic and thermoacidophilic biopolymer-degrading genes and enzymes from alicyclobacillus acidocaldarius and related organisms, methods. Google Patents; 2013.
Thompson DN, Thompson VS, Schaller KD, Apel WA, Reed DW, Lacey JA. Thermal and acid tolerant beta xylosidases, arabinofuranosidases, genes encoding, related organisms, and methods. Google Patents; 2013.
Deutscher J. The mechanisms of carbon catabolite repression in bacteria. Curr Opin Microbiol. 2008;11(2):87–93.
Gorke B, Stulke J. Carbon catabolite repression in bacteria: many ways to make the most out of nutrients. Nat Rev Microbiol. 2008;6(8):613–24.
Biely P. Microbial carbohydrate esterases deacetylating plant polysaccharides. Biotechnol Adv. 2012;30(6):1575–88.
Deli A, Koutsioulis D, Fadouloglou VE, Spiliotopoulou P, Balomenou S, Arnaouteli S, Tzanodaskalaki M, Mavromatis K, Kokkinidis M, Bouriotis V. LmbE proteins from Bacillus cereus are de-N-acetylases with broad substrate specificity and are highly similar to proteins in Bacillus anthracis. FEBS J. 2010;277(13):2740–53.
Viars S, Valentine J, Hernick M. Structure and function of the LmbE-like superfamily. Biomolecules. 2014;4(2):527–45.
Mine S, Niiyama M, Hashimoto W, Ikegami T, Koma D, Ohmoto T, Fukuda Y, Inoue T, Abe Y, Ueda T, et al. Expression from engineered Escherichia coli chromosome and crystallographic study of archaeal N,N′-diacetylchitobiose deacetylase. FEBS J. 2014;281(11):2584–96.
Beringer T, Lucht W. Bioenergy and biospheric carbon. In: Lal R, Lorenz K, Hüttl RF, Schneider BU, von Braun J, editors. Recarbonization of the biosphere. Dordrecht: Springer; 2012. p. 481–92.
Bhalla A, Bansal N, Kumar S, Bischoff KM, Sani RK. Improved lignocellulose conversion to biofuels with thermophilic bacteria and thermostable enzymes. Bioresour Technol. 2013;128:751–9.
Mandalari G, Faulds CB, Sancho AI, Saija A, Bisignano G, LoCurto R, Waldron KW. Fractionation and characterisation of arabinoxylans from brewers’ spent grain and wheat bran. J Cereal Sci. 2005;42(2):205–12.
Bai Y, Wang J, Zhang Z, Shi P, Luo H, Huang H, Feng Y, Yao B. Extremely acidic β-1,4-glucanase, CelA4, from thermoacidophilic Alicyclobacillus sp. A4 with high protease resistance and potential as a pig feed additive. J Agric Food Chem. 2010;58(3):1970–5.
Nurizzo D, Nagy T, Gilbert HJ, Davies GJ. The structural basis for catalysis and specificity of the Pseudomonas cellulosa α-glucuronidase, GlcA67A. Structure. 2002;10(4):547–56.
Zhang Y, Gao F, Xue Y, Zeng Y, Peng H, Qi J, Ma Y. Crystallization and preliminary X-ray study of native and selenomethionyl [beta]-1,4-mannanase AaManA from Alicyclobacillus acidocaldarius Tc-12-31. Acta Crystallogr Sect F Struct Biol Crystall Commun. 2008;64(3):209–12.
Nanavati DM, Thirangoon K, Noll KM. Several archaeal homologs of putative oligopeptide-binding proteins encoded by Thermotoga maritima bind sugars. Appl Environ Microbiol. 2006;72(2):1336–45.
VanFossen AL, Verhaart MRA, Kengen SMW, Kelly RM. Carbohydrate utilization patterns for the extremely thermophilic bacterium Caldicellulosiruptor saccharolyticus reveal broad growth substrate preferences. Appl Environ Microbiol. 2009;75(24):7718–24.
Bauer MW, Driskill LE, Callen W, Snead MA, Mathur EJ, Kelly RM. An endoglucanase, EglA, from the hyperthermophilic archaeon Pyrococcus furiosus hydrolyzes β-1,4 bonds in mixed-linkage (1 → 3), (1 → 4)-β-d-glucans and cellulose. J Bacteriol. 1999;181(1):284–90.
Liebl W, Ruile P, Bronnenmeier K, Riedel K, Lottspeich F, Greif I. Analysis of a Thermotoga maritima DNA fragment encoding two similar thermostable cellulases, CelA and CelB, and characterization of the recombinant enzymes. Microbiology. 1996;142(9):2533–42.
Hall J, Hazlewood GP, Barker PJ, Gilbert HJ. Conserved reiterated domains in Clostridium thermocellum endoglucanases are not essential for catalytic activity. Gene. 1988;69(1):29–38.
Shi P, Tian J, Yuan T, Liu X, Huang H, Bai Y, Yang P, Chen X, Wu N, Yao B. Paenibacillus sp. strain E18 bifunctional xylanase–glucanase with a single catalytic domain. Appl Environ Microbiol. 2010;76(11):3620–4.
Chang L, Ding M, Bao L, Chen Y, Zhou J, Lu H. Characterization of a bifunctional xylanase/endoglucanase from yak rumen microorganisms. Appl Microbiol Biotechnol. 2011;90(6):1933–42.
Rashamuse KJ, Visser DF, Hennessy F, Kemp J, Roux-van der Merwe MP, Badenhorst J, Ronneburg T, Francis-Pope R, Brady D. Characterisation of two bifunctional cellulase-xylanase enzymes isolated from a bovine rumen metagenome library. Curr Microbiol. 2013;66(2):145–51.
Breves R, Bronnenmeier K, Wild N, Lottspeich F, Staudenbauer WL, Hofemeister J. Genes encoding two different beta-glucosidases of Thermoanaerobacter brockii are clustered in a common operon. Appl Environ Microbiol. 1997;63(10):3902–10.
Han Y, Agarwal V, Dodd D, Kim J, Bae B, Mackie RI, Nair SK, Cann IKO. Biochemical and structural insights into xylan utilization by the thermophilic bacterium Caldanaerobius polysaccharolyticus. J Biol Chem. 2012;287(42):34946–60.
Mai V, Wiegel J, Lorenz WW. Cloning, sequencing, and characterization of the bifunctional xylosidase–arabinosidase from the anaerobic thermophile Thermoanaerobacter ethanolicus. Gene. 2000;247(1–2):137–43.
Bai Y, Huang H, Meng K, Shi P, Yang P, Luo H, Luo C, Feng Y, Zhang W, Yao B. Identification of an acidic α-amylase from Alicyclobacillus sp. A4 and assessment of its application in the starch industry. Food Chem. 2012;131(4):1473–8.
Larsbrink J, Izumi A, Ibatullin FM, Nakhai A, Gilbert HJ, Davies GJ, Brumer H. Structural and enzymatic characterization of a glycoside hydrolase family 31 α-xylosidase from Cellvibrio japonicus involved in xyloglucan saccharification. Biochem J. 2011;436(3):567–80.
Trincone A, Cobucci-Ponzano B, Di Lauro B, Rossi M, Mitsuishi Y, Moracci M. Enzymatic synthesis and hydrolysis of xylogluco-oligosaccharides using the first archaeal α-xylosidase from Sulfolobus solfataricus. Extremophiles. 2001;5(4):277–82.
Bolstad BM, Irizarry RA, Åstrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19(2):185–93.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003;4(2):249–64.
Markowitz VM, Chen IMA, Chu K, Szeto E, Palaniappan K, Grechkin Y, Ratner A, Jacob B, Pati A, Huntemann M, et al. IMG/M: the integrated metagenome data management and comparative analysis system. Nucleic Acids Res. 2012;40(D1):D123–9.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Gish W, States DJ. Identification of protein coding regions by database similarity search. Nat Genet. 1993;3(3):266–72.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004;5(1):113.
BL designed and performed the study, analyzed the data, and drafted the manuscript. WA, PS, and LD assisted in interpretation of data, provided direction for the studies, and revised the manuscript. All authors read and approved the final manuscript.
This work was supported by the Idaho National Laboratory’s Laboratory Directed Research and Development program under the Department of Energy Idaho Operations Office Contract DE-AC07-05ID14517.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and analyzed during the current study are available in the Gene Expression Omnibus (GEO) repository within NCBI.
Consent for publication
Ethics approval and consent to participate
Funding for the research was provided through the Idaho National Laboratory’s Laboratory Directed Research and Development Program.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.