BACKGROUND
the plant cell wall consists of four major structural polymer components, namely cellulose, hemicellulose, pectin and lignin
<cit> .
in contrast to cellulose, which is a homopolymer of β--linked glucose monomers, hemicelluloses incorporate a wide range of heterogeneous polymers containing a backbone of xylose, arabinose, galactose or, mannose monomers
<cit> .
the most common hemicelluloses are the xylans, which have a backbone of β--linked xylose monomers, and are frequently decorated with side chains of arabinose, glucuronic acid or its 4-o-methyl derivative, and acetyl side groups
<cit> .
the structure of the hemicellulose, its backbone monomers and degree of branching and type of side chains varies extensively between different plant species and even tissues within the same plant.
for example, xyloglucans dominate the primary cell walls of dicotyledonous plants, while glucuronoarabinoxylans are prevalent in both the primary and secondary cell walls of commelinid monocots .
while the primary cell walls of conifers such as pine are dominated by xyloglucans, the hemicellulose fraction of their secondary cell walls is comprised mainly of galactoglucomannans
<cit> .
similarly, pectin represents a family of heterogeneous polymers exhibiting complex structures which are found in the plant primary cell wall
<cit> .
in nature, plant biomass represents the most abundant carbohydrate source for microorganisms.
as such, many bacteria, fungi and yeast can degrade hemicellulose polymers through the production of a battery of hydrolytic enzymes, including xylanases, xylosidases, arabinofuranosidases and glucuronidases
<cit> .
in recent decades there has been increasing interest in the use of microorganisms and their enzymes for the hydrolysis of hemicelluloses, as a result of the various biotechnological and industrial applications in which the hydrolytic enzymes, the hemicelluloytic process and its products can be exploited.
for example, xylanases can be applied in the bleaching of pulp and paper, increasing the digestibility of animal feeds, clarification of fruit juices and the bioconversion of plant bionass
.
the latter can be exploited, in combination with the fermentative capacities of microorganisms in the conversion of lignocellulosic biomass to high yields of the renewable energy source bioethanol and other value-added products
<cit> .
the oligosaccharides generated from hemicellulose hydrolysis can also be applied as prebiotics that inhibit harmful microorganisms and stimulate the growth of bifidobacteria and lactobacilli in the human intestine
<cit> .
there has been increasing interest in the use of thermophilic microorganisms for the hemicellulytic hydrolysis of plant biomass, as they display a number of advantages over mesophilic microorganisms
<cit> .
the rate of biomass conversion is enhanced at higher growth temperatures.
furthermore, as high temperatures are required for the pre-treatment and delignification of plant biomass in many industrial processes, the ability of thermophilic microorganisms to tolerate higher temperatures, as well as producing thermostable hydrolytic enzymes, can potentially reduce process costs
<cit> .
obligately thermophilic species within the genus bacillus have been accommodated in the genus geobacillus, which currently encompasses nineteen species
<cit> .
members of this genus are rod-shaped, aerobic endo-spore formers and have optimum growth temperatures of between  <dig> and 75°c
<cit> .
geobacillus spp.
serve as an attractive resource for production of a wide range of thermostable enzymes for a variety of biotechnological applications
<cit> .
geobacillus strains have been shown to be highly effective in the degradation of hemicellulose polymers, and considerable attention has focused on the thermostable extra- and intra-cellular xylanases and β-xylosidases
.
the molecular determinants underlying the capacity of this genus to utilize hemicellulose have, however, been largely elucidated in only one strain, geobacillus stearothermophilus t- <dig>  and have been shown to reside on a single genetic locus
.
on the basis of characterization of some of the gene clusters within this locus, a model for the hemicellulolytic capacity of this strain has been developed
.
the model hemicellulose substrate, acetylglucuronoarabinoxylan , is cleaved by the extracellular gh <dig> xylanase xyna <dig> into large xylooligomers, including xylooligosaccharides and aldotetraouronic acid, which are subsequently transport across the cellular membrane by specific oligosaccharide transport systems
<cit> .
the glucuronic acid is cleaved from the aldotetrauronic acid -xylotriose or meglcaxyl3) polymer by the action of a gh <dig> α-glucuronidase and is metabolized via the uronate metabolic pathway, also encoded in the hemicellulose utilization locus
<cit> .
the xylooligosaccharides and xylotriose derived from meglcaxyl <dig> polymers are further cleaved into d-xylose monomers through the concerted action of a second intracellular gh <dig> xylanase, xyna <dig>  and three β-xylosidases xynb <dig>  xynb <dig> and xynb <dig>  which are then metabolized via the pentose phosphate and glycolysis pathways
<cit> .
pathways for the hydrolysis and utilization of arabinofuranosyl substituents and for the main chain deacetylation have also been described
<cit> .
a gene cluster within the hemicellulose utilization locus of g. stearothermophilus t- <dig> has also been shown to play a role in the degradation of the pectin-associated l-arabinan polymer
<cit> , further highlighting the capacity of geobacillus spp.
to degrade and utilize polymers in plant biomass.
to date, the genomes of twenty-four geobacillus spp.
have been sequenced.
here, using the g. stearothermophilus t- <dig> hemicellulose utilization system  as a model, we have identified and characterized loci coding for the degradation and utilization of hemicellulose in most of the sequenced strains.
our analyses reveal the extensive variability in the geobacillus hemicellulose utilization loci, which could be exploited to generate recombinant strains with improved hemicellulolytic capacities.
RESULTS
a hemicellulose utilization system is a common feature among most geobacillusspp
the capacity of g. stearothermophilus t- <dig> to degrade and utilize hemicellulose is dependent on a single complex gene locus, which has recently been sequenced .
this locus is ~76 kb in size, has a g + c content of  <dig> % and encodes  <dig> proteins.
the hus locus of t- <dig> can be roughly divided into thirteen distinct gene clusters .
twelve of the gene clusters have been characterized functionally and/or in silico and have been shown to play roles in the concerted degradation and utilization of the hemicellulose polymer glucuronoarabinoxylan in this strain
<cit> .
of these, seven can be considered as central to hemicellulose degradation, with the proteins encoded in these clusters driving the extracellular degradation of arabinoglucuronoxylan , transport of the cleaved oligosaccharide products into the cell  and further intracellular degradation of these intermediates into metabolizable monosaccharides .
four additional gene clusters  are not essential for the degradation of the hemicellulose polymer, but rather encode pathways for the metabolism of the end-product pentose sugars arabinose and xylose and for uronic acids, as well as a putative alternative pathway for pentose sugar metabolism
<cit> .
an additional cluster which is not essential for hemicellulose degradation, encodes a transporter of l-arabinose monosaccharides .
a gene cluster with a role in the degradation and utilization of the pectin-associated polymer l-arabinan has also been described in g. stearothermophilus t-6
<cit>  and is positioned in the hemicellulose utilization cluster of this strain .figure  <dig>
schematic diagram of the
g. stearothermophilus
t- <dig> hemicellulose utilization locus.
each arrow represents a gene in the locus, with genes encoding predicted transposons colored in black, while open reading frames interrupted by transposons are colored in red.
genes encoding glycosyl hydrolases are colored in blue.
the g. stearothermophilus hus locus was subdivided into thirteen gene clusters on the basis of their predicted function.
the genomes of  <dig> geobacillus strains, thirteen complete and eleven draft genomes, have been sequenced
.
the twelve gene clusters of t- <dig> with a reported role in hemicellulose utilization were used as units for the identification of orthologous clusters in the sequenced geobacillus genomes, by localized blastn analysis.
between three and twelve partial or complete orthologous clusters were found in twenty of the sequenced strains, while no orthologous clusters could be found in the genomes of the remaining four strains .
a further three strains, namely g. thermoglucosidasius ccb_us3_uf <dig> and b <dig>  as well as g. kaustophilus hta <dig> only carry three orthologous gene clusters, namely those for xylose metabolism, arabinose transport and l-arabinose metabolism.
the lack of orthologs of all t- <dig> clusters required for hemicellulose degradation suggests that these strains do not have this capacity, and rather make use of l-arabinose and d-xylose monomers that may be present in the environment.
the remaining seventeen sequenced geobacillus strains contain between eight and twelve gene clusters sharing orthology with the thirteen clusters that make up the hemicellulose utilization locus of g. stearothermophilus t- <dig> , including between five and seven of the clusters with a role in hemicellulose degradation.
the genomic locations of the orthologous clusters indicate that, as is the case of t- <dig>  they are arranged into a single chromosomal hemicellulose utilization locus.
thus, including g. stearothermophilus t- <dig>  eighteen geobacillus strains  can be considered to contain a hemicellulose utilization system .
in addition to t- <dig>  the genomes of three further sequenced geobacillus strains, namely y412mc <dig>  y412mc <dig> and wsucf <dig>  were observed to encode an l-arabinan utilization system  incorporated within the hemicellulose utilization locus.table  <dig>
genome characteristics of twenty-four sequenced
geobacillus
strains
g. thermoleovorans
g. thermoleovorans
g. kaustophilus
g. kaustophilus
g. thermodenitrificans
g. thermodenitrificans
g. caldoxylolyticus
g. thermoglucosidasius
g. thermoglucosidasius
g. thermoglucosidasius
ncbi acc # refers to the ncbi accession numbers of the genomes for the specified organisms.
arefers to the integrated microbial genomes  object identifier for the genome of g. thermoglucosidasius m10exg .
the number  of orthologous gene clusters from the
g
.
stearothermophilus
t- <dig> hus loci in sequenced
geobacillus
strains
g. thermodenitrificans dsm465tcomparative diagram of the
geobacillus
hemicellulose loci.
the aligned echd-npd islands of the twenty-four sequenced hus+ and hus-
geobacillus strains and that of g. stearothermophilus t- <dig> are shown.
the flanking echd and npd genes are colored in yellow, genes coding for transposons in black, and genes in which the reading frames are transposon-disrupted are in red.
those genes conserved among >70% of the hus+ strains are colored in green, those conserved among >50% and <70% in light green, while those conserved among <50% of the hus+ strains are colored in white.
a phylogenetic tree including all sequenced geobacillus spp.
and twelve geobacillus spp.
type strains was constructed using the recn gene.
phylogeny using the recn gene has been demonstrated to have great resolving power at the species and subs-species level for the genus geobacillus
<cit> .
the topology  showed that the hemicellulose utilization system is encoded across all major branches of the geobacillus genus phylogeny, although clusters of hus- strains can be observed in the g. thermoleovorans-kaustophilus and g. thermoglucosidasius-thermantarcticus-toebii branches.figure  <dig>
phylogeny of the sequenced
geobacillus
strains showing the distribution of hemicellulose loci.
a neighbour-joining phylogeny was constructed on the basis of the recn gene for the twenty-four sequenced geobacillus strains as well as twelve geobacillus spp.
type strains.
the presence or absence of hus loci in the sequenced strains is indicated by green and red dots, respectively.
the blue star next to g. stearothermophilus atcc 12980t indicates the predicted phylogenetic position of g. stearothermophilus t- <dig> for which a recn sequence is not available.
the recn gene sequence of bacillus subtilis  <dig> was used as outgroup.
bootstrap values  are shown.
the geobacillushemicellulose system forms part of a genomic island
genomic islands were identified in the genomes of the thirteen complete geobacillus genomes using the islandviewer server
<cit> .
between thirteen and thirty genomic islands were predicted for each strain.
one predicted island, flanked by genes encoding an enoyl-coa hydratase  and a nitropropane dioxygenases , is present on the chromosome of twelve out of thirteen complete sequenced geobacillus strains.
blastn analysis with these genes identified similar islands in all eleven draft genomes , suggesting that an echd-npd island is present in twenty-three of the twenty-four sequenced strains.
while orthologs of the echd and npd genes were detected in the genome sequence of geobacillus sp.
wch <dig>  these are separated by ~540 kb of sequence and were not identified as an island.
the echd-npd islands range in size from  <dig>  to  <dig>  kilobases and carry between eleven and sixty-eight protein coding sequences .
the g + c content of the echd-npd island in all geobacillus strains is on average  <dig>  ±  <dig> %  lower than the genomic g + c content, indicative of horizontal acquisition of this island.table  <dig>
characteristics of the
echd
-
npd
islands of twenty-five
geobacillus
strains
g. thermoleovorans
g. thermoleovorans
g. kaustophilus
g. kaustophilus
geobacillus stearothermophilus
g. thermodenitrificans
g. thermodenitrificans
g. caldoxylolyticus
g. thermoglucosidasius
g. thermoglucosidasius
g. thermoglucosidasius
blastn analysis with the echd and npd genes from the complete genomes against the g. stearothermophilus t- <dig> hus locus revealed that a complete copy of echd and a partial npd gene were present at the 5′ and 3′ ends of this locus, respectively.
a search up- and downstream of the hemicellulose degradation and utilization gene clusters in sixteen of the sequenced hus+ strains showed that, similarly, all of the hemicellulose utilization loci in these strains were maintained within the echd-npd genomic islands, the three geobacillus strains which contained orthologs of only three of the thirteen clusters in the t- <dig> hemicellulose utilization locus also carry these clusters within the echd-npd island, indicating that the echd-npd genomic islands serve as a centralized repository for the hemicellulose degradation and utilization genes within geobacillus spp.
one exception is the hus+ strain g. thermoglucosidasius c56ys <dig>  where the orthologous clusters for hemicellulose degradation and utilization are localized downstream of the npd-echd island and are flanked by transposon genes .
the average difference in g + c content between the echd-npd island and the genome is  <dig>  ±  <dig> % for hus- strains, while the mean difference for hus+ is somewhat larger .
an exception is again observed for g. thermoglucosidasius c56ys <dig>  where the g + c content of the echd-npd island is only  <dig> % below that of the genome, while that of the island plus the hus locus is  <dig> % above that the genome.
of note is the great variability in genomic g + c contents observed for the different sequenced geobacillus strains.
for example, the mean g + c content of the g. caldoxylolyticus cic <dig> genome is  <dig> %, and that of geobacillus sp.
a <dig> is  <dig> %, while the mean g + c contents of their echd-npd islands are  <dig> % and  <dig> %, respectively.
we suggest that these differences, given that the islands are highly variable, might be due to the non-conserved content of the island.
the mean g + c contents for  <dig> conserved genes and those of the non-conserved genes in the echd-npd islands of these two strains were determined.
this showed that g + c content differences existed between the islands of the a <dig> and cic <dig> strains for both conserved and non-conserved genes, but that there was negligible variation in the g + c contents of both the conserved and non-conserved gene complement in the two strains .
the concept of g + c amelioration
<cit>  implies that while the g + c contents of genomic islands that have recently been horizontally acquired reflect those of the donor strain, the genes in these islands are subject to the same mutational processes as all other genes in the recipient strain and will become more similar to that of the recipient genomes over time.
we suggest that the observed differences of genomic and island g + c contents among geobacillus strains is due to the fact that the echd-npd islands may have derived from a common ancestor with a much lower genomic g + c content, and that these islands show g + c contents which are more, or less, ameliorated, depending on the time since acquisition.
comparison of the hemicellulose utilization loci reveals extensive variation among geobacillusstrains
the protein coding sequences  on the echd-npd islands of all twenty-five geobacillus strains were predicted using fgenesb
<cit> .
the translated cds sets from both hus- and hus+ strains were compared by localized blastp analyses and the proteins were annotated on the basis of orthologous proteins identified by comparison against the ncbi non-redundant protein database .
no proteins were conserved among the echd-npd islands of all twenty-five analyzed strains, while only seven proteins are conserved among the hus+ strains.
as between  <dig> and  <dig> proteins are encoded on the echd-npd islands of hus+ strains, the low number of conserved proteins  highlights the extensive variability among the geobacillus hus loci.
cdss conserved among the eighteen hus+ strains include two coding for a two component system , three coding for a xylotriose transporter  and two genes  required for xylose metabolism.
the translated proteins products share an average amino acid identity of  <dig> % between pair-wise compared strains, although lower orthology is observed for the xyndcefg-xylab sets of g. caldoxylolyticus cic <dig> and geobacillus sp.
c56-t <dig>  which share  <dig> % and  <dig> % average amino acid identity with each of the other sixteen hus+ strains, respectively, and  <dig> % amino acid identity between them.
a neighbour-joining phylogeny was constructed on the basis of the concatenated xyndcefg-xylab amino acid sequences of the eighteen hus+geobacillus strains, as well as sixteen phylogenetic relatives for which whole genome sequences are available and in which orthologous proteins are present.
the resultant tree  shows congruence in the clustering of most hus+geobacillus strains with the associated recn tree, with the exception of geobacillus sp.
c56-t <dig> and g. caldoxylolyticus cic <dig>  suggesting that the hemicellulose utilization loci of sixteen strains have been derived from a common ancestor and have subsequently diverged in parallel with their genomes, while those of the other two strains are derived from a separate, but possibly common, ancestor.figure  <dig>
phylogenetic comparison of the hus
+
geobacillus
strains on the basis of recn and xyndcefg-xylab.
a neighbour-joining tree was constructed on the basis of the recn amino acid sequences of the hus+
geobacillus strains as well as sixteen representatives of related genera in which orthologs of xyndcefg-xylab are present.
similarly a neighbour-joining tree was constructed based on alignment of the concatenated xyndcefg-xylab amino acid sequences.
desmospora sp.
<dig> was used as outgroup and bootstrap values are shown .
the branches in red and in green, represent the phylogenetic positions of geobacillus sp.
c56-t <dig> and g. caldoxylolyticus cic <dig>  respectively, in both trees.
variability in the extracellular hydrolytic enzyme complement encoded on the geobacillus hus loci
the subcellular localizations of all the hus-encoded cdss were predicted by comparison of their translated protein products against the psortb  <dig>  server
<cit>  and the detection of signal peptides by submission to the signalp  <dig>  server
<cit> .
signalp analysis predicted seventeen distinct hdi-encoded proteins with typical gram-positive signal peptide/cleavage sites required for the sec-dependent transport of proteins across the cellular membrane.
psortb analysis showed that the majority of these  are integrated into the cellular membrane and we propose that these form part of transport systems for the uptake of oligosaccharides.
on the basis of psort and signalp results, two hus-encoded proteins are predicted to be secreted.
the protein sequences were compared against the carbohydrate active enzymes  database using the dbcan server’s blast tool
<cit> .
this showed that these two proteins represent glycoside hydrolases belonging to the glycoside hydrolase  families gh <dig> and gh <dig>  respectively.
orthologs of a thermostable extracellular gh <dig> family xylanase xyna <dig> that degrades xylan to its component xylooligosaccharides have been identified in a number of geobacillus spp.
.
an orthologous xyna <dig> gene is found in the hemicellulose utilization locus of 15/ <dig> hus+ strains.
this gene is ~ <dig>  kb in size and is transcribed on the opposite strand to all other genes in the hus loci.
truncated genes  showing sequence homology to the 3′ end of the xyna <dig> gene are also present in the hus loci of geobacillus sp.
a <dig> and wsucf <dig>  while no orthologs is present in g. caldoxylolyticus cic <dig> hdi, suggesting that the latter three strains lack the ability to degrade xylan extracellularly into its component xylooligosaccharides.
the translated protein products of the complete xyna <dig> genes share an average amino acid identity of  <dig> % among the hdi+ and  <dig> % to xyna <dig> of g. stearothermophilus t- <dig> .
an extracellular gh <dig> endo α- <dig> -arabinanase  has been shown to hydrolyze pectin-associated l-arabinan polymers into shorter arabino-saccharides and arabinose in g. stearothermophilus t-6
<cit>  and forms part of the aus locus in this strain.
orthologs of abna are also found in the other three aus+ strains and their encoded proteins share  <dig> % average amino acid identity to the t- <dig> arabinanase.
additional predicted orthologous glycoside hydrolases are encoded in the hus loci of geobacillus sp.
jf <dig>  mas <dig>  and g. thermoglucosidasius c <dig> ys <dig>  which share  <dig> % average amino acid identity among them.
the ~ <dig>  kb gene for this enzyme is localized in a non-conserved region of the hus locus between the genes coding for a predicted acetyl-esterase  and a β-xylosidase  .
a truncated copy of this gene is also present within this region in geobacillus sp.
camr <dig>  comparison of the translated protein products of the three complete gene copies against the cazy database using the dbcan blast tool
<cit>  showed that they code for a gh <dig> glycoside hydrolase, showing greatest orthology to a glycoside hydrolase family protein in verminephrobacter eiseniae ef01- <dig> .
the gh <dig> family includes enzymes active on a wide range of carbohydrate substrates, including chitosanases, β-mannosidases, xyloglucan-specific endo- β- <dig> -glucanases and cellulases
<cit> .
the target carbohydrate for a glycoside hydrolases can sometimes be further elucidated by the identification of carbohydrate binding modules , which promote the association of the enzyme with a particular substrate, within the protein sequence of the enzyme
<cit> .
comparison against the cazy database using the dbcan blast tool identified a module showing weak orthology to a cbm family  <dig> module of a predicted endoglucanase in ignavibacterium album jcm <dig>  at the n-terminal end  of the gh <dig> protein orthologs of all four geobacillus strains.
the cbm <dig> modules have been demonstrated to have binding specificities towards a number of substrates, including both branched and debranched xylan, β- <dig> -glucan , and agarose
<cit> .
it is thus difficult to distinguish the target carbohydrate substrate for the geobacillus gh <dig> enzyme.
similarly, it is difficult to determine the subcellular localizations of these enzymes.
signalp analysis shows that no typical gram-positive signal peptides are present.
the psortb subcellular localization scl-blast + tool, however, showed that the geobacillus proteins share orthology with the cellodextrinase cela of ruminococcus flavefaciens fd- <dig> , for which an extracellular localization has been determined experimentally
<cit> .
the limited number of secreted hydrolytic enzymes has been proposed to provide a selective advantage to geobacillus spp.
and related hemicellulolytic organisms, as the relatively large oligosaccharides which result from extracellular hemicellulose hydrolysis will require further cleavage by intracellular enzymes and are thus not readily available to competing, non-hemicellulolytic, microorganisms
<cit> .
the assimilation of xylo- and arabinoxylo-oligosaccharides may furthermore serve in the conservation of energy.
variability in the membrane transport systems encoded on the geobacillus hus locus
following the extracellular degradation of xylan into its component oligosaccharides, these intermediates are typically transported into the cell by means of specific abc sugar transporters
<cit> .
carbohydrate uptake transporters  have been classified into two distinct families, with cut <dig> signifying those transporters that facilitate the uptake of di- and oligosaccharides as well as glycerol phosphate and polyols, while members of the cut <dig> family transport only monosaccharides
<cit> .
eight gene clusters encoding distinct cut <dig> transporters could be identified among the hus+ strains, with between two and four of these present in each of the strains, while a single cut <dig> family transporter is encoded in most of the hus+ as well as three hus- strains.
three cut <dig> transporters with a role in transport of hemicellulose/pectic polymer degradation products have been described and functionally characterized in g. stearothermophilus t- <dig>  these are required for the uptake of arabinosaccharides , xylooligosaccharides  and aldotetraouronic acid -xylotriose – meglcaxyl3)
.
orthologs of the xynefg xylooligosaccharide transporter, which preferentially transports xylotriose, are found in all eighteen hus+ strains, while the aguefg aldotetraouronic acid transporter is present in sixteen hus+ strains.
the reading frame of the ague gene coding for an aldotetraouronic acid extracellular sugar-binding protein in t- <dig> is disrupted by a transposon, suggesting that this transporter is non-functional in this strain.
in contrast to the prevalent aguefg and xynefg transport systems, the abnefj arabinosaccharide transporter, which forms part of the l-arabinan utilization system, is restricted in distribution to three strains, namely t- <dig>  y412mc <dig> and y412mc <dig>  one cut <dig> transporter has also been identified in the hus locus of g. stearothermophilus t- <dig>  which is predicted to transport arabinose monosaccharide
<cit> .
the polycistronic operon  coding for this transporter is also disrupted in t- <dig> by a transposon insertion.
orthologs are present in ten out of the  <dig> hdi+ strains, and also in three hdi- strains, suggesting that this transporter is not a characteristic of hemicelluloses degrading strains, but may form a part of an arabinose metabolic pathway present in most geobacillus spp.
the five remaining distinct cut <dig> transporters encoded in the hus loci of the different sequenced geobacillus strains share only limited orthology to characterized transporters.
for further reference, the gene clusters encoding these transporters were named cutabc, cutdef, cutghi, cutjkl, cutmno in order to distinguish between them and to avoid duplicity of gene names .
a locus in bacillus subtilis encodes a cut <dig> transporter for the uptake of α- <dig> -arabinooligosaccharides
<cit> .
orthologous loci  are present in 10/ <dig> hus+ strains, all of which are missing abnefj orthologs, with the translated protein products sharing  <dig> % average amino acid identity with the aranpq proteins of b. subtilis  <dig>  in geobacillus sp.
camr <dig>  only a partial cutb and complete cutc gene are present, while an ortholog of cuta is missing.
the b. subtilis aranpq transport system is required for the uptake of linear α- <dig> -arabinosaccharide homopolymers consisting of three or four arabinosyl units residues, but not the branched arabinooligosaccharides decorated with α- <dig> - and α- <dig> -linked arabinosyl side chains.
these branched arabinooligosaccharides are transported into the cell by an as yet unknown mechanism
<cit> .
the nature of the substrate transported by the geobacillus cutabc system is yet to be experimentally determined.
a distinct cut <dig> transporter  is encoded within the aus locus of geobacillus sp.
wsucf <dig> in the location occupied by abnefg in the other aus+ strains.
the translated products for these genes share only  <dig> % average amino acid identity to g. stearothermophilus t- <dig> abnefj, but share greater orthology  with an abc transporter of an unknown sugar substrate in bacillus halodurans c- <dig> .
the location of the cut- <dig> encoding genes between the arabinanase-encoding abna and abnb genes in wsucf <dig>  however, suggests that as in the case of t- <dig>  they encode an arabinosaccharide transporter, although whether this cut <dig> system facilitates the uptake of linear or branched saccharides is unknown.
the cutdef genes encoding a cut <dig> transporter are located between the axe <dig> and xynb <dig> genes in the loci of 5/ <dig> hus+ strains.
their translated products share  <dig> % average amino acid identity with a characterized aldotetrouronate transporter in paenibacillus jdr-2
<cit> , but share only  <dig> % average amino acid identity with the aguefg aldotetrouronate transporter of g. stearothermophilus t- <dig>  four of the strains with a cutabc cluster also encode aguefg orthologs, suggesting they possess two alternative pathways for the uptake of aldotetraouronic acid.
on the other hand, c56-t <dig> lacks aguefg orthologs, but does appear to carry a functional cutdef system.
the translated protein products of the predicted cutmno cut <dig> transporter  in the hus loci of in geobacillus sp.
mas <dig> and jf <dig> share  <dig> % amino acid identity between the two strains but show limited orthology  to the xynefg transporter in g. stearothermophilus t- <dig>  the localization of the cutmno genes adjacent to the gene encoding the predicted gh <dig> hydrolytic enzyme in these two strains, suggests they likely code for a transporter for the uptake of the degraded saccharide products of the latter enzyme.
the other two strains which encode a gh <dig> ortholog, namely geobacillus sp.
camr <dig> and g. thermoglucosidasius c56ys <dig> do not contain an orthologous transporter.
three genes encoding an additional orthologous cut <dig> transporter  are localized between arardba and abfa-arajklmn in geobacillus sp.
hh <dig> and between aram and aran in g. caldoxylolyticus cic <dig>  respectively .
the encoded proteins share  <dig> % average amino acid identity between these two strains and  <dig> % average amino acid identity with a predicted transporter of an unknown sugar substrate in paenibacillus mucilaginosis  <dig> .
variability in the intracellular enzyme complement of hus+ geobacillus spp
typically, following the intracellular transport of the component oligosaccharides they are further degraded into monomeric sugars through the concerted action of a battery of intracellular glycosyl hydrolases
<cit> .
the intracellular enzyme complement for the cleavage of xylooligosacharides, methylglucuronate side chains and l-arabinan arabinosaccharides have been characterized in g. stearothermophilus t-6
.
the aguefg-transported substrate meglcaxyl <dig> is cleaved into d-glucuronic acid and xylotriose by the gh <dig> family α-glucuronidase agua
<cit> .
orthologs of the t- <dig> agua gene are located in the loci of all hus+ strains, with the exception of geobacillus sp.
jf <dig>  as this strain is also missing orthologs of aguefg it is likely that this strain cannot utilize this substrate.
the geobacillus sp.
c56-t <dig> hus locus does encode an agua ortholog, but is also missing orthologs of the aguefg genes, suggesting this strain can also not utilize aldotetraouronic acid as a metabolic intermediate, although it cannot be excluded that an alternative transporter for this substrate may exist in this strain.
the xylotriose derived from cleavage of meglcuax <dig>  as well as the xylooligosaccharides transported into the cell via the xynefg transporter, are degraded to d-xylose through the concerted action of an intracellular xylanase and β-xylosidases
<cit> .
the xyna <dig> gene in g. stearothermophilus t- <dig> encodes an intracellular gh <dig> family xylanase that degrades xylooligosaccharides, including xylotriose, into xylose and xylobiose subunits
<cit> .
orthologs are encoded in all hus+ strains with the exception of geobacillus sp.
jf <dig> and g. caldoxylolyticus cic <dig>  the structure, biochemistry and mechanism of action of three distinct β-xylosidases encoded in the hus locus of g. stearothermophilus t- <dig>  xynb <dig> , xynb <dig>  and xynb <dig> , have been described
.
the natural substrate of xynb <dig> and xynb <dig> is xylobiose, while that of xynb <dig> is xylotriose.
the xynb <dig> gene forms part of the xylooligosaccharide degradation cluster  of the loci of all eighteen hus+geobacillus strains, and as such probably degrades xylobiose released through the intracellular hydrolysis of xynefg-transported xylooligosaccharides by xyna <dig>  the xynb <dig> gene forms part of the aldotetraouronic acid utilization cluster  and the ability of xynb <dig> to degrade xylotriose, the cleavage product of meglcuax <dig>  may preclude the requirement of an additional intracellular xylanase for the cleavage of this substrate.
a copy of xynb <dig> is present in 17/ <dig> hus+ strains.
the presence of a xynb <dig> ortholog in g. caldoxylolyticus cic <dig> suggests its protein product may complement the xylotriose cleaving action of the missing xyna <dig> enzyme.
in geobacillus sp.
jf <dig>  both xynb <dig> and xyna <dig> orthologs are missing, implying that this strain is incapable of utilizing the xylotriose substrate.
complete copies of xynb <dig> are found in 14/ <dig> hus+geobacillus strains.
in mas <dig>  y412mc <dig>  y412mc61and c56t <dig> the xynb <dig> reading frame is disrupted by a transposon insertion, but given the similar function of the two β-xylosidases xynb <dig> and xynb <dig> in the hydrolysis of xylobiose, the complete xynb <dig> copy in these strains may complement the disrupted xynb <dig> product.
three intracellular enzymes in the aus locus of g. stearothermophilus t- <dig> play a role in the concerted release of component arabinose monosaccharides from the arabinan-derived oligosaccharides
<cit> .
the gh <dig> arabinofuranosidase abf <dig> of b. subtilis  <dig>  has been shown to remove o- <dig> and o-3-linked arabinofuranosyl side chains of branched arabinan polymers
<cit> .
the translated protein products of the abfb gene found in all four aus+geobacillus strains share  <dig> % average amino acid identity with abf <dig>  suggesting a similar function for abfb proteins.
while the majority of arabinose in plant cell wall polymers exists in the furanose  conformation, a small fraction of arabinose in the side-chain terminal ends of arabinan is in the six-membered ring pyranose form
<cit> .
orthologs of a gh <dig> family arabinopyranosidase  which can cleave the terminal arabinopyranose residues
<cit>  are encoded in the aus loci of all four strains, but the presence of a transposon integrated in the t- <dig> abp gene suggests this capacity is disrupted in this strain.
finally, the resultant debranched linear α- <dig> -arabinofuranose chain is cleaved into its component arabinose monomers by the intracellular gh <dig> arabinanase abnb
<cit> , orthologs of which are encoded in all four aus+ strains.
orthologs of a gh <dig> family arabinofuranosidase  are encoded in seventeen hus+ strains.
the absence of additional encoded enzymes for the utilization of l-arabinan polymers in most strains suggests that this enzyme plays a role in the liberation of o- <dig> and o-3-linked arabinofuranosyl residues from arabinoxylan
<cit> .
the presence of both abfa and abfb orthologs in the aus+ strains thus suggests that they may have a capacity for debranching both arabinan and arabinoxylan polymers.
the hus locus of c56-t <dig> lacks orthologs of both abfa and abfb, but instead carries a gene encoding a distinct gh <dig> family glycoside hydrolase.
the encoded protein shares only  <dig> % and  <dig> % amino acid identity with the abfa and abfb proteins of g. stearothermophilus t- <dig>  respectively.
however, it shares orthology with the abf43b protein of paenibacillus sp.
e <dig> , which was characterized as a gh <dig> family α-arabinofuranosidase with a role in degradation of arabinoxylans
<cit> .
in nature, hemicellulose polymers are frequently substituted with acetyl side chains.
for example, the main hemicellulose of hardwoods, 4-o-methyl-d-glucuronoxylan, contains acetyl side chains at positions two or three of 70% of the xylose backbone subunits
<cit> .
these acetyl side chains contribute towards the recalcitrance of hemicelluloses by restricting access of hydrolytic enzymes to the backbone sugars.
as a ‘counteractive’ measure, many microorganisms produce acetylesterase enzymes, which can cleave the acetyl side chains
<cit> .
two genes encoding predicted acetylesterases have been identified within the g. stearothermophilus t- <dig> hemicellulose utilization locus, axe <dig> and axe2
<cit> , with orthologous genes present in 16/ <dig> and 17/ <dig> hus+ strains, respectively.
the axe <dig> gene forms part of xylooligosaccharide utilization cluster  and the encoded axe <dig> protein products share  <dig> % average amino acid identity among the hus+ strains.
comparison against the cazy database with the dbcan blast tool shows that they belong to carbohydrate esterase family ce4
<cit> , which incorporates acetylxylan esterases and chitin deacetylases, and show greatest orthology with the predicted polysaccharide deacetylase of paenibacillus sp.
jdr- <dig> .
the incorporation of the axe <dig> gene in the xylooligosaccharide utilization cluster  suggests that this esterase plays a role in the deacetylation of xylan or its internalized derivative oligosaccharides.
the crystal structure of the intracellular axe <dig> protein of g. stearothermophilus t- <dig> has been determined, and this functionally characterized acetylxylan esterase has been shown to belong to a novel ce family
<cit> .
orthologs of the axe <dig> proteins of sixteen sequenced hus+ strains share  <dig> % average amino acid identity with the t- <dig> axe <dig> protein, and negligible sequence similarity to the axe <dig> protein in this strain.
the axe <dig> gene forms part of the second intracellular xylooligosaccharide degradation cluster .
while the axe <dig> protein has been shown to be active on short acetylated xylo-oligmers
<cit> , the natural target xylooligosaccharide is unknown.
the native substrate of axe <dig> is also unknown, but it can be hypothesized that these two esterases act in concert to completely deacetylate the xylooligomers resulting from xylan degradation.
a further putative intracellular hydrolytic enzyme is encoded in the hus locus of jf <dig>  while a transposon-disrupted copy is also found in mas <dig>  the translated jf <dig> protein product lacks a recognizable gram-positive signal peptide, suggesting an intracellular localization, and shares 73% amino acid identity with a gh <dig> family glycosyl hydrolase of thermoanaerobacterium thermosaccharolyticum dsm <dig> .
furthermore, the amino acid sequence contains a β-glucosidase conserved domain .
the predicted jf <dig> and mas <dig> β-glucosidase-encoding  gene is located downstream of gh <dig> and cutdef genes, suggesting the encoded enzyme is likely involved in the intracellular hydrolysis of the substrate derived through the extracellular hydrolysis by the gh <dig> enzyme and intracellular transport by the potentially associated cutabc transporter.
orthologs of the jf <dig> bgl gene and the cutdef transporter are, however, absent in geobacillus sp.
camr <dig> and g. thermoglucosidasius c56ys <dig>  which contain a partial and complete copy of the gh5-encoding gene, respectively.
members of the gh <dig> family include glucosidases and xylosidases, and as such the target substrate of this gene cluster and its encoded proteins, whether it is a glucose- or xylose-based polymer, remains undetermined.
variability in the metabolic genes in the hus loci of geobacillus spp
the monosaccharides derived from intracellular oligosaccharide cleavage are typically catabolised into intermediate substrates that enter the pentose phosphate and embden meyerhof parnas  glycolysis pathways
<cit> .
the catabolic systems for d-xylose , uronic acids  and l-arabinose  monosaccharides have been well-characterized in geobacillus stearothermophilus and its phylogenetic relative b. subtilis
<cit> .
orthologs of the xylab and kdgka-uxac-uxuab loci are present in all eighteen hus+ strains, while arardba orthologs are present in seventeen of the strains.
the latter locus is absent from g. thermoglucosidasius c56y <dig>  which is also missing orthologs of all additional genes encoding enzymes required for the transport and catabolism of arabinosaccharides, suggesting that the hemicellulose degradation system functions on arabinofuranose-free hemicellulose substrates.
in y412mc <dig> and y412mc <dig>  a transposase is inserted into the reading frame of uxac, which encodes a uronate isomerase that catalyzes the first step in the conversion of uronic acids to glycerol-3-phosphate and pyruvate, which can enter the emp pathway
<cit> .
this suggests that these two strains are unable to catabolise this substrate.
an additional gene cluster, arajklmn, has been identified in the g. stearothermophilus t- <dig> hus locus.
complete arajklmn loci are present in eleven of the eighteen hus+ strains, while partial or transposon-disrupted loci are present in 6/ <dig> hus+ strains.
both arak and aran orthologs are absent from strains ng- <dig>  g11mc <dig> and dsm <dig>  while an aran ortholog is missing from the c56-t <dig> hus locus.
in c56-t <dig> and camr <dig>  the aran reading frame is disrupted by a transposon.
a second copy of the aran gene, arano is located downstream of the cutabc locus of ghh <dig>  with its translated protein product sharing  <dig> % amino acid identity with the aran protein in this strain.
similarly, a second copy  is located upstream of arapst in the g. caldoxylolyticus cic <dig> hus locus, but its translated product shares lower sequence similarity  with the cic <dig> aran protein.
the genes in this cluster have been shown to encode orthologs of a nadp sugar dehydrogenase , aldose-1-epimerase , sugar phosphatase , nadp-dependent glycerol-1-dehydrogenase  and hypothetical protein  and have been postulated to constitute an alternative pathway for the catabolism of pentose sugars
<cit> .
additional genes in the geobacillus echd-npdislands may play a role in hemicellulose utilization and alcohol fermentation
additional non-conserved genes can be observed at the 5′ end of the echd-npd islands of several strains, flanking the hemicellulose utilization locus .
their translated protein products were annotated by blastp comparison against the ncbi non-redundant protein database to identify orthologs for which functions have been predicted or experimentally determined .
a gene cluster  at the 5′ end of the echd-npd island of five hus+geobacillus strains, namely y412mc <dig>  y412mc <dig>  mas <dig>  ghh <dig> and c56ys <dig>  encodes a predicted oligopeptide transporter.
orthologs could furthermore be identified in the three hus- strains, geobacillus sp.
y <dig> mc <dig>  g. thermoglucosidasius tno- <dig>  and m10exg.
the echd-npd island 5′ end regions of g. caldoxylolyticus cic <dig>  geobacillus sp.
g11mc <dig>  and g. thermodenitrificans ng80- <dig> and dsm465t harbor large inserts,  <dig>  to 28 kb in size.
the genes in these inserts encode a system for the catabolism of the polyol myo-inositol , as well as two predicted abc sugar transporters.
strains g11mc <dig>  ng80- <dig> and dsm465t also encode a  <dig>  amino acid protein which, by comparison against the cazy database, was identified as a putative gh <dig> family glycoside hydrolase.
the gh <dig> protein products show extensive orthology to the functionally characterized α-mannosidase  of bacillus sp.
gl <dig> .
one of these genes codes for an ortholog of the aldehyde dehydrogenase aldh-t of g. stearothermophilus sic <dig> .
this enzyme has been shown to play a role in the conversion of acetyl coa to aldehydes which are further converted into alcohols by alcohol dehydrogenases
<cit> .
furthermore, another gene in this cluster encodes an orthologs of the aldo/keto-reductase ytbe of bacillus sp.
ecu <dig> , which has been shown to catalyze the nadh-dependent reduction of aliphatic and aromatic aldehydes into their corresponding alcohols
<cit> .
the co-localization of genes with a putative role in alcohol fermentation and hemicellulose utilization further suggests a role for the echd-npd island in the centralization, and potentially maintenance and regulation, of the metabolic and fermentative pathways of geobacillus spp.
CONCLUSIONS
our comparative genomic analyses have revealed the presence of a centralized hemicellulose utilization locus in most sequenced geobacillus strains.
however, this locus, which is localized on a common genomic island, displays extensive genetic variability.
genes or pathways central to the utilization of the hemicellulose backbone and substituents are absent from the hus loci of many strains, while other genes, or operons, are frequently disrupted by transposon insertion.
from this we may deduce that different geobacillus strains have a highly variable capacity to degrade hemicellulose polymers.
it is possible that, in their natural environments, geobacillus strains form part of consortia, whereby strains deficient in genes or pathways required for hemicellulose degradation and utilization may be complemented by other members of the consortia.
however, given the limited number of extracellular hydrolytic enzymes produced by geobacillus, which are generally conserved across their hus loci, and the specific trans-membrane transport systems they utilize for the internalization of relatively large oligosaccharides, we argue that such complementation is unlikely.
comparisons of the hus loci of the sequenced strains to that of the well-characterized g. stearothermophilus t- <dig> revealed the presence of a number of additional genes, some of which have orthologs in related genera and code for both extra- and intracellular hydrolytic enzymes and transmembrane transporters.
these findings further support the concept of variable hemicellulolytic capacities among geobacillus strains, and imply that they have evolved to degrade distinct hemicellulose substrates, such as are found in different plant species and tissues.
for example, the lack of genes encoding an α-glucuronidase , its associated transport system  and uronate metabolic genes in geobacillus sp.
jf <dig> implies that this strain utilizes a xylan which is not substituted with glucuronic acid, while a xylan polymer devoid of arabinofuranosyl substituents is the likely target of g. thermoglucosidasius c56ys <dig>  this novel substrate-guided ‘host-parasite’ relationship has not been demonstrated experimentally, but is potentially testable.
a recent review highlighted the extensive geographical and environmental diversity of members of the genus geobacillus and their important ecological roles in composting of plant biomass under thermophilic conditions
<cit> .
their ability to degrade distinct hemicellulose polymers may thus reflect their environmental exposure to particular hemicelluloses and potentially, their role in carbon cycling in places where hot conditions persist.
the extensive genetic variability among the geobacillus hemicellulose utilization systems highlighted in this study can potentially serve as a backbone for the development of recombinant geobacillus strain with improved hemicellulolytic capacities or with the capacity to degrade structurally diverse hemicellulose substrates.
such engineered strains, with the added advantage of the thermostability of many of the geobacillus extra- and intracellular hydrolytic enzymes, could then be exploited in various biotechnological applications, such as the production of bioethanol and other value-added products.

 the plant cell wall consists of four major structural polymer components, namely cellulose, hemicellulose, pectin and lignin. in contrast to cellulose, which is a homopolymer of β--linked glucose monomers, hemicelluloses incorporate a wide range of heterogeneous polymers containing a backbone of xylose, arabinose, galactose or, mannose monomers. the most common hemicelluloses are the xylans, which have a backbone of β--linked xylose monomers, and are frequently decorated with side chains of arabinose, glucuronic acid or its 4-o-methyl derivative, and acetyl side groups. the structure of the hemicellulose, its backbone monomers and degree of branching and type of side chains varies extensively between different plant species and even tissues within the same plant. for example, xyloglucans dominate the primary cell walls of dicotyledonous plants, while glucuronoarabinoxylans are prevalent in both the primary and secondary cell walls of commelinid monocots. while the primary cell walls of conifers such as pine are dominated by xyloglucans, the hemicellulose fraction of their secondary cell walls is comprised mainly of galactoglucomannans. similarly, pectin represents a family of heterogeneous polymers exhibiting complex structures which are found in the plant primary cell wall. in nature, plant biomass represents the most abundant carbohydrate source for microorganisms. as such, many bacteria, fungi and yeast can degrade hemicellulose polymers through the production of a battery of hydrolytic enzymes, including xylanases, xylosidases, arabinofuranosidases and glucuronidases. in recent decades there has been increasing interest in the use of microorganisms and their enzymes for the hydrolysis of hemicelluloses, as a result of the various biotechnological and industrial applications in which the hydrolytic enzymes, the hemicelluloytic process and its products can be exploited. for example, xylanases can be applied in the bleaching of pulp and paper, increasing the digestibility of animal feeds, clarification of fruit juices and the bioconversion of plant bionass. the latter can be exploited, in combination with the fermentative capacities of microorganisms in the conversion of lignocellulosic biomass to high yields of the renewable energy source bioethanol and other value-added products. the oligosaccharides generated from hemicellulose hydrolysis can also be applied as prebiotics that inhibit harmful microorganisms and stimulate the growth of bifidobacteria and lactobacilli in the human intestine. there has been increasing interest in the use of thermophilic microorganisms for the hemicellulytic hydrolysis of plant biomass, as they display a number of advantages over mesophilic microorganisms. the rate of biomass conversion is enhanced at higher growth temperatures. furthermore, as high temperatures are required for the pre-treatment and delignification of plant biomass in many industrial processes, the ability of thermophilic microorganisms to tolerate higher temperatures, as well as producing thermostable hydrolytic enzymes, can potentially reduce process costs. obligately thermophilic species within the genus bacillus have been accommodated in the genus geobacillus, which currently encompasses nineteen species. members of this genus are rod-shaped, aerobic endo-spore formers and have optimum growth temperatures of between and 75°c. geobacillus spp. serve as an attractive resource for production of a wide range of thermostable enzymes for a variety of biotechnological applications. geobacillus strains have been shown to be highly effective in the degradation of hemicellulose polymers, and considerable attention has focused on the thermostable extra- and intra-cellular xylanases and β-xylosidases. the molecular determinants underlying the capacity of this genus to utilize hemicellulose have, however, been largely elucidated in only one strain, geobacillus stearothermophilus t- and have been shown to reside on a single genetic locus. on the basis of characterization of some of the gene clusters within this locus, a model for the hemicellulolytic capacity of this strain has been developed. the model hemicellulose substrate, acetylglucuronoarabinoxylan, is cleaved by the extracellular gh xylanase xyna into large xylooligomers, including xylooligosaccharides and aldotetraouronic acid, which are subsequently transport across the cellular membrane by specific oligosaccharide transport systems. the glucuronic acid is cleaved from the aldotetrauronic acid -xylotriose or meglcaxyl3) polymer by the action of a gh α-glucuronidase and is metabolized via the uronate metabolic pathway, also encoded in the hemicellulose utilization locus. the xylooligosaccharides and xylotriose derived from meglcaxyl polymers are further cleaved into d-xylose monomers through the concerted action of a second intracellular gh xylanase, xyna and three β-xylosidases xynb xynb and xynb which are then metabolized via the pentose phosphate and glycolysis pathways. pathways for the hydrolysis and utilization of arabinofuranosyl substituents and for the main chain deacetylation have also been described. a gene cluster within the hemicellulose utilization locus of g. stearothermophilus t- has also been shown to play a role in the degradation of the pectin-associated l-arabinan polymer, further highlighting the capacity of geobacillus spp. to degrade and utilize polymers in plant biomass. to date, the genomes of twenty-four geobacillus spp. have been sequenced. here, using the g. stearothermophilus t- hemicellulose utilization system as a model, we have identified and characterized loci coding for the degradation and utilization of hemicellulose in most of the sequenced strains. our analyses reveal the extensive variability in the geobacillus hemicellulose utilization loci, which could be exploited to generate recombinant strains with improved hemicellulolytic capacities.  a hemicellulose utilization system is a common feature among most geobacillusspp the capacity of g. stearothermophilus t- to degrade and utilize hemicellulose is dependent on a single complex gene locus, which has recently been sequenced. this locus is ~76 kb in size, has a g + c content of % and encodes proteins. the hus locus of t- can be roughly divided into thirteen distinct gene clusters. twelve of the gene clusters have been characterized functionally and/or in silico and have been shown to play roles in the concerted degradation and utilization of the hemicellulose polymer glucuronoarabinoxylan in this strain. of these, seven can be considered as central to hemicellulose degradation, with the proteins encoded in these clusters driving the extracellular degradation of arabinoglucuronoxylan, transport of the cleaved oligosaccharide products into the cell and further intracellular degradation of these intermediates into metabolizable monosaccharides. four additional gene clusters are not essential for the degradation of the hemicellulose polymer, but rather encode pathways for the metabolism of the end-product pentose sugars arabinose and xylose and for uronic acids, as well as a putative alternative pathway for pentose sugar metabolism. an additional cluster which is not essential for hemicellulose degradation, encodes a transporter of l-arabinose monosaccharides. a gene cluster with a role in the degradation and utilization of the pectin-associated polymer l-arabinan has also been described in g. stearothermophilus t-6 and is positioned in the hemicellulose utilization cluster of this strain.figure schematic diagram of the g. stearothermophilus t- hemicellulose utilization locus. each arrow represents a gene in the locus, with genes encoding predicted transposons colored in black, while open reading frames interrupted by transposons are colored in red. genes encoding glycosyl hydrolases are colored in blue. the g. stearothermophilus hus locus was subdivided into thirteen gene clusters on the basis of their predicted function. the genomes of geobacillus strains, thirteen complete and eleven draft genomes, have been sequenced. the twelve gene clusters of t- with a reported role in hemicellulose utilization were used as units for the identification of orthologous clusters in the sequenced geobacillus genomes, by localized blastn analysis. between three and twelve partial or complete orthologous clusters were found in twenty of the sequenced strains, while no orthologous clusters could be found in the genomes of the remaining four strains. a further three strains, namely g. thermoglucosidasius ccb_us3_uf and b as well as g. kaustophilus hta only carry three orthologous gene clusters, namely those for xylose metabolism, arabinose transport and l-arabinose metabolism. the lack of orthologs of all t- clusters required for hemicellulose degradation suggests that these strains do not have this capacity, and rather make use of l-arabinose and d-xylose monomers that may be present in the environment. the remaining seventeen sequenced geobacillus strains contain between eight and twelve gene clusters sharing orthology with the thirteen clusters that make up the hemicellulose utilization locus of g. stearothermophilus t-, including between five and seven of the clusters with a role in hemicellulose degradation. the genomic locations of the orthologous clusters indicate that, as is the case of t- they are arranged into a single chromosomal hemicellulose utilization locus. thus, including g. stearothermophilus t- eighteen geobacillus strains can be considered to contain a hemicellulose utilization system. in addition to t- the genomes of three further sequenced geobacillus strains, namely y412mc y412mc and wsucf were observed to encode an l-arabinan utilization system incorporated within the hemicellulose utilization locus.table genome characteristics of twenty-four sequenced geobacillus strains g. thermoleovorans g. thermoleovorans g. kaustophilus g. kaustophilus g. thermodenitrificans g. thermodenitrificans g. caldoxylolyticus g. thermoglucosidasius g. thermoglucosidasius g. thermoglucosidasius ncbi acc # refers to the ncbi accession numbers of the genomes for the specified organisms. arefers to the integrated microbial genomes object identifier for the genome of g. thermoglucosidasius m10exg. the number of orthologous gene clusters from the g. stearothermophilus t- hus loci in sequenced geobacillus strains g. thermodenitrificans dsm465tcomparative diagram of the geobacillus hemicellulose loci. the aligned echd-npd islands of the twenty-four sequenced hus+ and hus- geobacillus strains and that of g. stearothermophilus t- are shown. the flanking echd and npd genes are colored in yellow, genes coding for transposons in black, and genes in which the reading frames are transposon-disrupted are in red. those genes conserved among >70% of the hus+ strains are colored in green, those conserved among >50% and <70% in light green, while those conserved among <50% of the hus+ strains are colored in white. a phylogenetic tree including all sequenced geobacillus spp. and twelve geobacillus spp. type strains was constructed using the recn gene. phylogeny using the recn gene has been demonstrated to have great resolving power at the species and subs-species level for the genus geobacillus. the topology showed that the hemicellulose utilization system is encoded across all major branches of the geobacillus genus phylogeny, although clusters of hus- strains can be observed in the g. thermoleovorans-kaustophilus and g. thermoglucosidasius-thermantarcticus-toebii branches.figure phylogeny of the sequenced geobacillus strains showing the distribution of hemicellulose loci. a neighbour-joining phylogeny was constructed on the basis of the recn gene for the twenty-four sequenced geobacillus strains as well as twelve geobacillus spp. type strains. the presence or absence of hus loci in the sequenced strains is indicated by green and red dots, respectively. the blue star next to g. stearothermophilus atcc 12980t indicates the predicted phylogenetic position of g. stearothermophilus t- for which a recn sequence is not available. the recn gene sequence of bacillus subtilis was used as outgroup. bootstrap values are shown. the geobacillushemicellulose system forms part of a genomic island genomic islands were identified in the genomes of the thirteen complete geobacillus genomes using the islandviewer server. between thirteen and thirty genomic islands were predicted for each strain. one predicted island, flanked by genes encoding an enoyl-coa hydratase and a nitropropane dioxygenases, is present on the chromosome of twelve out of thirteen complete sequenced geobacillus strains. blastn analysis with these genes identified similar islands in all eleven draft genomes, suggesting that an echd-npd island is present in twenty-three of the twenty-four sequenced strains. while orthologs of the echd and npd genes were detected in the genome sequence of geobacillus sp. wch these are separated by ~540 kb of sequence and were not identified as an island. the echd-npd islands range in size from to kilobases and carry between eleven and sixty-eight protein coding sequences. the g + c content of the echd-npd island in all geobacillus strains is on average  ±  % lower than the genomic g + c content, indicative of horizontal acquisition of this island.table characteristics of the echd - npd islands of twenty-five geobacillus strains g. thermoleovorans g. thermoleovorans g. kaustophilus g. kaustophilus geobacillus stearothermophilus g. thermodenitrificans g. thermodenitrificans g. caldoxylolyticus g. thermoglucosidasius g. thermoglucosidasius g. thermoglucosidasius blastn analysis with the echd and npd genes from the complete genomes against the g. stearothermophilus t- hus locus revealed that a complete copy of echd and a partial npd gene were present at the 5′ and 3′ ends of this locus, respectively. a search up- and downstream of the hemicellulose degradation and utilization gene clusters in sixteen of the sequenced hus+ strains showed that, similarly, all of the hemicellulose utilization loci in these strains were maintained within the echd-npd genomic islands, the three geobacillus strains which contained orthologs of only three of the thirteen clusters in the t- hemicellulose utilization locus also carry these clusters within the echd-npd island, indicating that the echd-npd genomic islands serve as a centralized repository for the hemicellulose degradation and utilization genes within geobacillus spp. one exception is the hus+ strain g. thermoglucosidasius c56ys where the orthologous clusters for hemicellulose degradation and utilization are localized downstream of the npd-echd island and are flanked by transposon genes. the average difference in g + c content between the echd-npd island and the genome is  ±  % for hus- strains, while the mean difference for hus+ is somewhat larger. an exception is again observed for g. thermoglucosidasius c56ys where the g + c content of the echd-npd island is only % below that of the genome, while that of the island plus the hus locus is % above that the genome. of note is the great variability in genomic g + c contents observed for the different sequenced geobacillus strains. for example, the mean g + c content of the g. caldoxylolyticus cic genome is %, and that of geobacillus sp. a is %, while the mean g + c contents of their echd-npd islands are % and %, respectively. we suggest that these differences, given that the islands are highly variable, might be due to the non-conserved content of the island. the mean g + c contents for conserved genes and those of the non-conserved genes in the echd-npd islands of these two strains were determined. this showed that g + c content differences existed between the islands of the a and cic strains for both conserved and non-conserved genes, but that there was negligible variation in the g + c contents of both the conserved and non-conserved gene complement in the two strains. the concept of g + c amelioration implies that while the g + c contents of genomic islands that have recently been horizontally acquired reflect those of the donor strain, the genes in these islands are subject to the same mutational processes as all other genes in the recipient strain and will become more similar to that of the recipient genomes over time. we suggest that the observed differences of genomic and island g + c contents among geobacillus strains is due to the fact that the echd-npd islands may have derived from a common ancestor with a much lower genomic g + c content, and that these islands show g + c contents which are more, or less, ameliorated, depending on the time since acquisition. comparison of the hemicellulose utilization loci reveals extensive variation among geobacillusstrains the protein coding sequences on the echd-npd islands of all twenty-five geobacillus strains were predicted using fgenesb. the translated cds sets from both hus- and hus+ strains were compared by localized blastp analyses and the proteins were annotated on the basis of orthologous proteins identified by comparison against the ncbi non-redundant protein database. no proteins were conserved among the echd-npd islands of all twenty-five analyzed strains, while only seven proteins are conserved among the hus+ strains. as between and proteins are encoded on the echd-npd islands of hus+ strains, the low number of conserved proteins highlights the extensive variability among the geobacillus hus loci. cdss conserved among the eighteen hus+ strains include two coding for a two component system, three coding for a xylotriose transporter and two genes required for xylose metabolism. the translated proteins products share an average amino acid identity of % between pair-wise compared strains, although lower orthology is observed for the xyndcefg-xylab sets of g. caldoxylolyticus cic and geobacillus sp. c56-t which share % and % average amino acid identity with each of the other sixteen hus+ strains, respectively, and % amino acid identity between them. a neighbour-joining phylogeny was constructed on the basis of the concatenated xyndcefg-xylab amino acid sequences of the eighteen hus+geobacillus strains, as well as sixteen phylogenetic relatives for which whole genome sequences are available and in which orthologous proteins are present. the resultant tree shows congruence in the clustering of most hus+geobacillus strains with the associated recn tree, with the exception of geobacillus sp. c56-t and g. caldoxylolyticus cic suggesting that the hemicellulose utilization loci of sixteen strains have been derived from a common ancestor and have subsequently diverged in parallel with their genomes, while those of the other two strains are derived from a separate, but possibly common, ancestor.figure phylogenetic comparison of the hus + geobacillus strains on the basis of recn and xyndcefg-xylab. a neighbour-joining tree was constructed on the basis of the recn amino acid sequences of the hus+ geobacillus strains as well as sixteen representatives of related genera in which orthologs of xyndcefg-xylab are present. similarly a neighbour-joining tree was constructed based on alignment of the concatenated xyndcefg-xylab amino acid sequences. desmospora sp. was used as outgroup and bootstrap values are shown. the branches in red and in green, represent the phylogenetic positions of geobacillus sp. c56-t and g. caldoxylolyticus cic respectively, in both trees. variability in the extracellular hydrolytic enzyme complement encoded on the geobacillus hus loci the subcellular localizations of all the hus-encoded cdss were predicted by comparison of their translated protein products against the psortb server and the detection of signal peptides by submission to the signalp server. signalp analysis predicted seventeen distinct hdi-encoded proteins with typical gram-positive signal peptide/cleavage sites required for the sec-dependent transport of proteins across the cellular membrane. psortb analysis showed that the majority of these are integrated into the cellular membrane and we propose that these form part of transport systems for the uptake of oligosaccharides. on the basis of psort and signalp results, two hus-encoded proteins are predicted to be secreted. the protein sequences were compared against the carbohydrate active enzymes database using the dbcan server’s blast tool. this showed that these two proteins represent glycoside hydrolases belonging to the glycoside hydrolase families gh and gh respectively. orthologs of a thermostable extracellular gh family xylanase xyna that degrades xylan to its component xylooligosaccharides have been identified in a number of geobacillus spp.. an orthologous xyna gene is found in the hemicellulose utilization locus of 15/ hus+ strains. this gene is ~  kb in size and is transcribed on the opposite strand to all other genes in the hus loci. truncated genes showing sequence homology to the 3′ end of the xyna gene are also present in the hus loci of geobacillus sp. a and wsucf while no orthologs is present in g. caldoxylolyticus cic hdi, suggesting that the latter three strains lack the ability to degrade xylan extracellularly into its component xylooligosaccharides. the translated protein products of the complete xyna genes share an average amino acid identity of % among the hdi+ and % to xyna of g. stearothermophilus t-. an extracellular gh endo α- -arabinanase has been shown to hydrolyze pectin-associated l-arabinan polymers into shorter arabino-saccharides and arabinose in g. stearothermophilus t-6 and forms part of the aus locus in this strain. orthologs of abna are also found in the other three aus+ strains and their encoded proteins share % average amino acid identity to the t- arabinanase. additional predicted orthologous glycoside hydrolases are encoded in the hus loci of geobacillus sp. jf mas and g. thermoglucosidasius c ys which share % average amino acid identity among them. the ~  kb gene for this enzyme is localized in a non-conserved region of the hus locus between the genes coding for a predicted acetyl-esterase and a β-xylosidase. a truncated copy of this gene is also present within this region in geobacillus sp. camr comparison of the translated protein products of the three complete gene copies against the cazy database using the dbcan blast tool showed that they code for a gh glycoside hydrolase, showing greatest orthology to a glycoside hydrolase family protein in verminephrobacter eiseniae ef01-. the gh family includes enzymes active on a wide range of carbohydrate substrates, including chitosanases, β-mannosidases, xyloglucan-specific endo- β- -glucanases and cellulases. the target carbohydrate for a glycoside hydrolases can sometimes be further elucidated by the identification of carbohydrate binding modules, which promote the association of the enzyme with a particular substrate, within the protein sequence of the enzyme. comparison against the cazy database using the dbcan blast tool identified a module showing weak orthology to a cbm family module of a predicted endoglucanase in ignavibacterium album jcm at the n-terminal end of the gh protein orthologs of all four geobacillus strains. the cbm modules have been demonstrated to have binding specificities towards a number of substrates, including both branched and debranched xylan, β- -glucan, and agarose. it is thus difficult to distinguish the target carbohydrate substrate for the geobacillus gh enzyme. similarly, it is difficult to determine the subcellular localizations of these enzymes. signalp analysis shows that no typical gram-positive signal peptides are present. the psortb subcellular localization scl-blast + tool, however, showed that the geobacillus proteins share orthology with the cellodextrinase cela of ruminococcus flavefaciens fd-, for which an extracellular localization has been determined experimentally. the limited number of secreted hydrolytic enzymes has been proposed to provide a selective advantage to geobacillus spp. and related hemicellulolytic organisms, as the relatively large oligosaccharides which result from extracellular hemicellulose hydrolysis will require further cleavage by intracellular enzymes and are thus not readily available to competing, non-hemicellulolytic, microorganisms. the assimilation of xylo- and arabinoxylo-oligosaccharides may furthermore serve in the conservation of energy. variability in the membrane transport systems encoded on the geobacillus hus locus following the extracellular degradation of xylan into its component oligosaccharides, these intermediates are typically transported into the cell by means of specific abc sugar transporters. carbohydrate uptake transporters have been classified into two distinct families, with cut signifying those transporters that facilitate the uptake of di- and oligosaccharides as well as glycerol phosphate and polyols, while members of the cut family transport only monosaccharides. eight gene clusters encoding distinct cut transporters could be identified among the hus+ strains, with between two and four of these present in each of the strains, while a single cut family transporter is encoded in most of the hus+ as well as three hus- strains. three cut transporters with a role in transport of hemicellulose/pectic polymer degradation products have been described and functionally characterized in g. stearothermophilus t- these are required for the uptake of arabinosaccharides, xylooligosaccharides and aldotetraouronic acid -xylotriose – meglcaxyl3). orthologs of the xynefg xylooligosaccharide transporter, which preferentially transports xylotriose, are found in all eighteen hus+ strains, while the aguefg aldotetraouronic acid transporter is present in sixteen hus+ strains. the reading frame of the ague gene coding for an aldotetraouronic acid extracellular sugar-binding protein in t- is disrupted by a transposon, suggesting that this transporter is non-functional in this strain. in contrast to the prevalent aguefg and xynefg transport systems, the abnefj arabinosaccharide transporter, which forms part of the l-arabinan utilization system, is restricted in distribution to three strains, namely t- y412mc and y412mc one cut transporter has also been identified in the hus locus of g. stearothermophilus t- which is predicted to transport arabinose monosaccharide. the polycistronic operon coding for this transporter is also disrupted in t- by a transposon insertion. orthologs are present in ten out of the hdi+ strains, and also in three hdi- strains, suggesting that this transporter is not a characteristic of hemicelluloses degrading strains, but may form a part of an arabinose metabolic pathway present in most geobacillus spp. the five remaining distinct cut transporters encoded in the hus loci of the different sequenced geobacillus strains share only limited orthology to characterized transporters. for further reference, the gene clusters encoding these transporters were named cutabc, cutdef, cutghi, cutjkl, cutmno in order to distinguish between them and to avoid duplicity of gene names. a locus in bacillus subtilis encodes a cut transporter for the uptake of α- -arabinooligosaccharides. orthologous loci are present in 10/ hus+ strains, all of which are missing abnefj orthologs, with the translated protein products sharing % average amino acid identity with the aranpq proteins of b. subtilis in geobacillus sp. camr only a partial cutb and complete cutc gene are present, while an ortholog of cuta is missing. the b. subtilis aranpq transport system is required for the uptake of linear α- -arabinosaccharide homopolymers consisting of three or four arabinosyl units residues, but not the branched arabinooligosaccharides decorated with α- - and α- -linked arabinosyl side chains. these branched arabinooligosaccharides are transported into the cell by an as yet unknown mechanism. the nature of the substrate transported by the geobacillus cutabc system is yet to be experimentally determined. a distinct cut transporter is encoded within the aus locus of geobacillus sp. wsucf in the location occupied by abnefg in the other aus+ strains. the translated products for these genes share only % average amino acid identity to g. stearothermophilus t- abnefj, but share greater orthology with an abc transporter of an unknown sugar substrate in bacillus halodurans c-. the location of the cut- encoding genes between the arabinanase-encoding abna and abnb genes in wsucf however, suggests that as in the case of t- they encode an arabinosaccharide transporter, although whether this cut system facilitates the uptake of linear or branched saccharides is unknown. the cutdef genes encoding a cut transporter are located between the axe and xynb genes in the loci of 5/ hus+ strains. their translated products share % average amino acid identity with a characterized aldotetrouronate transporter in paenibacillus jdr-2, but share only % average amino acid identity with the aguefg aldotetrouronate transporter of g. stearothermophilus t- four of the strains with a cutabc cluster also encode aguefg orthologs, suggesting they possess two alternative pathways for the uptake of aldotetraouronic acid. on the other hand, c56-t lacks aguefg orthologs, but does appear to carry a functional cutdef system. the translated protein products of the predicted cutmno cut transporter in the hus loci of in geobacillus sp. mas and jf share % amino acid identity between the two strains but show limited orthology to the xynefg transporter in g. stearothermophilus t- the localization of the cutmno genes adjacent to the gene encoding the predicted gh hydrolytic enzyme in these two strains, suggests they likely code for a transporter for the uptake of the degraded saccharide products of the latter enzyme. the other two strains which encode a gh ortholog, namely geobacillus sp. camr and g. thermoglucosidasius c56ys do not contain an orthologous transporter. three genes encoding an additional orthologous cut transporter are localized between arardba and abfa-arajklmn in geobacillus sp. hh and between aram and aran in g. caldoxylolyticus cic respectively. the encoded proteins share % average amino acid identity between these two strains and % average amino acid identity with a predicted transporter of an unknown sugar substrate in paenibacillus mucilaginosis. variability in the intracellular enzyme complement of hus+ geobacillus spp typically, following the intracellular transport of the component oligosaccharides they are further degraded into monomeric sugars through the concerted action of a battery of intracellular glycosyl hydrolases. the intracellular enzyme complement for the cleavage of xylooligosacharides, methylglucuronate side chains and l-arabinan arabinosaccharides have been characterized in g. stearothermophilus t-6. the aguefg-transported substrate meglcaxyl is cleaved into d-glucuronic acid and xylotriose by the gh family α-glucuronidase agua. orthologs of the t- agua gene are located in the loci of all hus+ strains, with the exception of geobacillus sp. jf as this strain is also missing orthologs of aguefg it is likely that this strain cannot utilize this substrate. the geobacillus sp. c56-t hus locus does encode an agua ortholog, but is also missing orthologs of the aguefg genes, suggesting this strain can also not utilize aldotetraouronic acid as a metabolic intermediate, although it cannot be excluded that an alternative transporter for this substrate may exist in this strain. the xylotriose derived from cleavage of meglcuax as well as the xylooligosaccharides transported into the cell via the xynefg transporter, are degraded to d-xylose through the concerted action of an intracellular xylanase and β-xylosidases. the xyna gene in g. stearothermophilus t- encodes an intracellular gh family xylanase that degrades xylooligosaccharides, including xylotriose, into xylose and xylobiose subunits. orthologs are encoded in all hus+ strains with the exception of geobacillus sp. jf and g. caldoxylolyticus cic the structure, biochemistry and mechanism of action of three distinct β-xylosidases encoded in the hus locus of g. stearothermophilus t- xynb, xynb and xynb, have been described. the natural substrate of xynb and xynb is xylobiose, while that of xynb is xylotriose. the xynb gene forms part of the xylooligosaccharide degradation cluster of the loci of all eighteen hus+geobacillus strains, and as such probably degrades xylobiose released through the intracellular hydrolysis of xynefg-transported xylooligosaccharides by xyna the xynb gene forms part of the aldotetraouronic acid utilization cluster and the ability of xynb to degrade xylotriose, the cleavage product of meglcuax may preclude the requirement of an additional intracellular xylanase for the cleavage of this substrate. a copy of xynb is present in 17/ hus+ strains. the presence of a xynb ortholog in g. caldoxylolyticus cic suggests its protein product may complement the xylotriose cleaving action of the missing xyna enzyme. in geobacillus sp. jf both xynb and xyna orthologs are missing, implying that this strain is incapable of utilizing the xylotriose substrate. complete copies of xynb are found in 14/ hus+geobacillus strains. in mas y412mc y412mc61and c56t the xynb reading frame is disrupted by a transposon insertion, but given the similar function of the two β-xylosidases xynb and xynb in the hydrolysis of xylobiose, the complete xynb copy in these strains may complement the disrupted xynb product. three intracellular enzymes in the aus locus of g. stearothermophilus t- play a role in the concerted release of component arabinose monosaccharides from the arabinan-derived oligosaccharides. the gh arabinofuranosidase abf of b. subtilis has been shown to remove o- and o-3-linked arabinofuranosyl side chains of branched arabinan polymers. the translated protein products of the abfb gene found in all four aus+geobacillus strains share % average amino acid identity with abf suggesting a similar function for abfb proteins. while the majority of arabinose in plant cell wall polymers exists in the furanose conformation, a small fraction of arabinose in the side-chain terminal ends of arabinan is in the six-membered ring pyranose form. orthologs of a gh family arabinopyranosidase which can cleave the terminal arabinopyranose residues are encoded in the aus loci of all four strains, but the presence of a transposon integrated in the t- abp gene suggests this capacity is disrupted in this strain. finally, the resultant debranched linear α- -arabinofuranose chain is cleaved into its component arabinose monomers by the intracellular gh arabinanase abnb, orthologs of which are encoded in all four aus+ strains. orthologs of a gh family arabinofuranosidase are encoded in seventeen hus+ strains. the absence of additional encoded enzymes for the utilization of l-arabinan polymers in most strains suggests that this enzyme plays a role in the liberation of o- and o-3-linked arabinofuranosyl residues from arabinoxylan. the presence of both abfa and abfb orthologs in the aus+ strains thus suggests that they may have a capacity for debranching both arabinan and arabinoxylan polymers. the hus locus of c56-t lacks orthologs of both abfa and abfb, but instead carries a gene encoding a distinct gh family glycoside hydrolase. the encoded protein shares only % and % amino acid identity with the abfa and abfb proteins of g. stearothermophilus t- respectively. however, it shares orthology with the abf43b protein of paenibacillus sp. e, which was characterized as a gh family α-arabinofuranosidase with a role in degradation of arabinoxylans. in nature, hemicellulose polymers are frequently substituted with acetyl side chains. for example, the main hemicellulose of hardwoods, 4-o-methyl-d-glucuronoxylan, contains acetyl side chains at positions two or three of 70% of the xylose backbone subunits. these acetyl side chains contribute towards the recalcitrance of hemicelluloses by restricting access of hydrolytic enzymes to the backbone sugars. as a ‘counteractive’ measure, many microorganisms produce acetylesterase enzymes, which can cleave the acetyl side chains. two genes encoding predicted acetylesterases have been identified within the g. stearothermophilus t- hemicellulose utilization locus, axe and axe2, with orthologous genes present in 16/ and 17/ hus+ strains, respectively. the axe gene forms part of xylooligosaccharide utilization cluster and the encoded axe protein products share % average amino acid identity among the hus+ strains. comparison against the cazy database with the dbcan blast tool shows that they belong to carbohydrate esterase family ce4, which incorporates acetylxylan esterases and chitin deacetylases, and show greatest orthology with the predicted polysaccharide deacetylase of paenibacillus sp. jdr-. the incorporation of the axe gene in the xylooligosaccharide utilization cluster suggests that this esterase plays a role in the deacetylation of xylan or its internalized derivative oligosaccharides. the crystal structure of the intracellular axe protein of g. stearothermophilus t- has been determined, and this functionally characterized acetylxylan esterase has been shown to belong to a novel ce family. orthologs of the axe proteins of sixteen sequenced hus+ strains share % average amino acid identity with the t- axe protein, and negligible sequence similarity to the axe protein in this strain. the axe gene forms part of the second intracellular xylooligosaccharide degradation cluster. while the axe protein has been shown to be active on short acetylated xylo-oligmers, the natural target xylooligosaccharide is unknown. the native substrate of axe is also unknown, but it can be hypothesized that these two esterases act in concert to completely deacetylate the xylooligomers resulting from xylan degradation. a further putative intracellular hydrolytic enzyme is encoded in the hus locus of jf while a transposon-disrupted copy is also found in mas the translated jf protein product lacks a recognizable gram-positive signal peptide, suggesting an intracellular localization, and shares 73% amino acid identity with a gh family glycosyl hydrolase of thermoanaerobacterium thermosaccharolyticum dsm. furthermore, the amino acid sequence contains a β-glucosidase conserved domain. the predicted jf and mas β-glucosidase-encoding gene is located downstream of gh and cutdef genes, suggesting the encoded enzyme is likely involved in the intracellular hydrolysis of the substrate derived through the extracellular hydrolysis by the gh enzyme and intracellular transport by the potentially associated cutabc transporter. orthologs of the jf bgl gene and the cutdef transporter are, however, absent in geobacillus sp. camr and g. thermoglucosidasius c56ys which contain a partial and complete copy of the gh5-encoding gene, respectively. members of the gh family include glucosidases and xylosidases, and as such the target substrate of this gene cluster and its encoded proteins, whether it is a glucose- or xylose-based polymer, remains undetermined. variability in the metabolic genes in the hus loci of geobacillus spp the monosaccharides derived from intracellular oligosaccharide cleavage are typically catabolised into intermediate substrates that enter the pentose phosphate and embden meyerhof parnas glycolysis pathways. the catabolic systems for d-xylose, uronic acids and l-arabinose monosaccharides have been well-characterized in geobacillus stearothermophilus and its phylogenetic relative b. subtilis. orthologs of the xylab and kdgka-uxac-uxuab loci are present in all eighteen hus+ strains, while arardba orthologs are present in seventeen of the strains. the latter locus is absent from g. thermoglucosidasius c56y which is also missing orthologs of all additional genes encoding enzymes required for the transport and catabolism of arabinosaccharides, suggesting that the hemicellulose degradation system functions on arabinofuranose-free hemicellulose substrates. in y412mc and y412mc a transposase is inserted into the reading frame of uxac, which encodes a uronate isomerase that catalyzes the first step in the conversion of uronic acids to glycerol-3-phosphate and pyruvate, which can enter the emp pathway. this suggests that these two strains are unable to catabolise this substrate. an additional gene cluster, arajklmn, has been identified in the g. stearothermophilus t- hus locus. complete arajklmn loci are present in eleven of the eighteen hus+ strains, while partial or transposon-disrupted loci are present in 6/ hus+ strains. both arak and aran orthologs are absent from strains ng- g11mc and dsm while an aran ortholog is missing from the c56-t hus locus. in c56-t and camr the aran reading frame is disrupted by a transposon. a second copy of the aran gene, arano is located downstream of the cutabc locus of ghh with its translated protein product sharing % amino acid identity with the aran protein in this strain. similarly, a second copy is located upstream of arapst in the g. caldoxylolyticus cic hus locus, but its translated product shares lower sequence similarity with the cic aran protein. the genes in this cluster have been shown to encode orthologs of a nadp sugar dehydrogenase, aldose-1-epimerase, sugar phosphatase, nadp-dependent glycerol-1-dehydrogenase and hypothetical protein and have been postulated to constitute an alternative pathway for the catabolism of pentose sugars. additional genes in the geobacillus echd-npdislands may play a role in hemicellulose utilization and alcohol fermentation additional non-conserved genes can be observed at the 5′ end of the echd-npd islands of several strains, flanking the hemicellulose utilization locus. their translated protein products were annotated by blastp comparison against the ncbi non-redundant protein database to identify orthologs for which functions have been predicted or experimentally determined. a gene cluster at the 5′ end of the echd-npd island of five hus+geobacillus strains, namely y412mc y412mc mas ghh and c56ys encodes a predicted oligopeptide transporter. orthologs could furthermore be identified in the three hus- strains, geobacillus sp. y mc g. thermoglucosidasius tno- and m10exg. the echd-npd island 5′ end regions of g. caldoxylolyticus cic geobacillus sp. g11mc and g. thermodenitrificans ng80- and dsm465t harbor large inserts, to 28 kb in size. the genes in these inserts encode a system for the catabolism of the polyol myo-inositol, as well as two predicted abc sugar transporters. strains g11mc ng80- and dsm465t also encode a amino acid protein which, by comparison against the cazy database, was identified as a putative gh family glycoside hydrolase. the gh protein products show extensive orthology to the functionally characterized α-mannosidase of bacillus sp. gl. one of these genes codes for an ortholog of the aldehyde dehydrogenase aldh-t of g. stearothermophilus sic. this enzyme has been shown to play a role in the conversion of acetyl coa to aldehydes which are further converted into alcohols by alcohol dehydrogenases. furthermore, another gene in this cluster encodes an orthologs of the aldo/keto-reductase ytbe of bacillus sp. ecu, which has been shown to catalyze the nadh-dependent reduction of aliphatic and aromatic aldehydes into their corresponding alcohols. the co-localization of genes with a putative role in alcohol fermentation and hemicellulose utilization further suggests a role for the echd-npd island in the centralization, and potentially maintenance and regulation, of the metabolic and fermentative pathways of geobacillus spp.  our comparative genomic analyses have revealed the presence of a centralized hemicellulose utilization locus in most sequenced geobacillus strains. however, this locus, which is localized on a common genomic island, displays extensive genetic variability. genes or pathways central to the utilization of the hemicellulose backbone and substituents are absent from the hus loci of many strains, while other genes, or operons, are frequently disrupted by transposon insertion. from this we may deduce that different geobacillus strains have a highly variable capacity to degrade hemicellulose polymers. it is possible that, in their natural environments, geobacillus strains form part of consortia, whereby strains deficient in genes or pathways required for hemicellulose degradation and utilization may be complemented by other members of the consortia. however, given the limited number of extracellular hydrolytic enzymes produced by geobacillus, which are generally conserved across their hus loci, and the specific trans-membrane transport systems they utilize for the internalization of relatively large oligosaccharides, we argue that such complementation is unlikely. comparisons of the hus loci of the sequenced strains to that of the well-characterized g. stearothermophilus t- revealed the presence of a number of additional genes, some of which have orthologs in related genera and code for both extra- and intracellular hydrolytic enzymes and transmembrane transporters. these findings further support the concept of variable hemicellulolytic capacities among geobacillus strains, and imply that they have evolved to degrade distinct hemicellulose substrates, such as are found in different plant species and tissues. for example, the lack of genes encoding an α-glucuronidase, its associated transport system and uronate metabolic genes in geobacillus sp. jf implies that this strain utilizes a xylan which is not substituted with glucuronic acid, while a xylan polymer devoid of arabinofuranosyl substituents is the likely target of g. thermoglucosidasius c56ys this novel substrate-guided ‘host-parasite’ relationship has not been demonstrated experimentally, but is potentially testable. a recent review highlighted the extensive geographical and environmental diversity of members of the genus geobacillus and their important ecological roles in composting of plant biomass under thermophilic conditions. their ability to degrade distinct hemicellulose polymers may thus reflect their environmental exposure to particular hemicelluloses and potentially, their role in carbon cycling in places where hot conditions persist. the extensive genetic variability among the geobacillus hemicellulose utilization systems highlighted in this study can potentially serve as a backbone for the development of recombinant geobacillus strain with improved hemicellulolytic capacities or with the capacity to degrade structurally diverse hemicellulose substrates. such engineered strains, with the added advantage of the thermostability of many of the geobacillus extra- and intracellular hydrolytic enzymes, could then be exploited in various biotechnological applications, such as the production of bioethanol and other value-added products. 
the genetic variability among the geobacillus hus loci implies that they have variable capacities to degrade hemicellulose polymers, or that they may degrade distinct polymers, as are found in different plant species and tissues. comparative analyses of these loci revealed extensive variability among the geobacillus hemicellulose utilization systems, with only seven out of 41– proteins encoded on these loci conserved among the hus+ strains. these hus loci are localized on a common genomic island. members of the thermophilic genus geobacillus can grow at high temperatures and produce a battery of thermostable hemicellulose hydrolytic enzymes, making them ideal candidates for the bioconversion of biomass to value-added products. this translates into extensive differences in the hydrolytic enzymes, transport systems and metabolic pathways employed by geobacillus spp. to date the molecular determinants for hemicellulose degradation and utilization have only been identified and partially characterized in one strain, namely geobacillus stearothermophilus t- where they are clustered in a single genetic locus. to degrade and utilize hemicellulose polymers. the data from this study can serve as a basis for the genetic engineering of a geobacillus strain with an improved capacity to degrade and utilize hemicellulose.
 members of the thermophilic genus geobacillus can grow at high temperatures and produce a battery of thermostable hemicellulose hydrolytic enzymes, making them ideal candidates for the bioconversion of biomass to value-added products. to date the molecular determinants for hemicellulose degradation and utilization have only been identified and partially characterized in one strain, namely geobacillus stearothermophilus t- where they are clustered in a single genetic locus.  using the g. stearothermophilus t- hemicellulose utilization locus as genetic marker, orthologous hemicellulose utilization loci were identified in the complete and partial genomes of 17/ geobacillus strains. these hus loci are localized on a common genomic island. comparative analyses of these loci revealed extensive variability among the geobacillus hemicellulose utilization systems, with only seven out of 41– proteins encoded on these loci conserved among the hus+ strains. this translates into extensive differences in the hydrolytic enzymes, transport systems and metabolic pathways employed by geobacillus spp. to degrade and utilize hemicellulose polymers.  the genetic variability among the geobacillus hus loci implies that they have variable capacities to degrade hemicellulose polymers, or that they may degrade distinct polymers, as are found in different plant species and tissues. the data from this study can serve as a basis for the genetic engineering of a geobacillus strain with an improved capacity to degrade and utilize hemicellulose. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords geobacillushemicellulosearabinoxylanarabinanxylanaseacetylesterasearabinofuranoseissue-copyright-statement© the author 2014 