# Glycan complexity dictates microbial resource allocation in the large intestine

Artur Rogowski, Jonathon A. Briggs, Jennifer C. Mortimer, Theodora Tryfona, Nicolas Terrapon, Elisabeth C. Lowe, Arnaud Basle, Carl Morland, Alison M. Day, Hongjun Zheng, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Artur Rogowski, Jonathon A. Briggs, Jennifer C. Mortimer, Theodora Tryfona, Nicolas Terrapon, et al.. Glycan complexity dictates microbial resource allocation in the large intestine. Nature Communications, 2015, 6, 10.1038/ncomms8481. hal-01439017

HAL Id: hal-01439017

https://hal.science/hal-01439017

Submitted on 7 Jan 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Glycan complexity dictates microbial resource allocation in the large intestine

Artur Rogowski<fnref n="1" />,*, Jonathon A. Briggs<fnref n="1" />,*, Jennifer C. Mortimer², Theodora Tryfona², Nicolas Terrapon³, Elisabeth C. Lowe<fnref n="1" />, Arnaud Baslé<fnref n="1" />, Carl Morland<fnref n="1" />, Alison M. Day<fnref n="1" />, Hongjun Zheng<fnref n="1" />, Theresa E. Rogers⁴, Paul Thompson<fnref n="1" />, Alastair R. Hawkins<fnref n="1" />, Madhav P. Yadav⁵, Bernard Henrissat³,⁶, Eric C. Martens⁴, Paul Dupree², Harry J. Gilbert<fnref n="1" /> & David N. Bolam<fnref n="1" />

The structure of the human gut microbiota is controlled primarily through the degradation of complex dietary carbohydrates, but the extent to which carbohydrate breakdown products are shared between members of the microbiota is unclear. We show here, using xylan as a model, that sharing the breakdown products of complex carbohydrates by key members of the microbiota, such as *Bacteroides ovatus*, is dependent on the complexity of the target glycan. Characterization of the extensive xylan degrading apparatus expressed by *B. ovatus* reveals that the breakdown of the polysaccharide by the human gut microbiota is significantly more complex than previous models suggested, which were based on the deconstruction of xylans containing limited monosaccharide side chains. Our report presents a highly complex and dynamic xylan degrading apparatus that is fine-tuned to recognize the different forms of the polysaccharide presented to the human gut microbiota.

<footnote n="1">¹Institute for Cell and Molecular Biosciences, The Medical School, Newcastle University, Newcastle upon Tyne NE2 4HH, UK. ²Department of Biochemistry, University of Cambridge, Tennis Court Road, Cambridge CB2 1QW, UK. ³Architecture et Fonction des Macromolécules Biologiques (AFMB), UMR 7257 CNRS, Université Aix-Marseille, 163 Avenue de Luminy, 13288 Marseille, France. ⁴Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, Michigan 48109, USA. ⁵United States Department of Agriculture- Agricultural Research Service, Eastern Regional Research Center, 600 East Mermaid Lane, Wyndmoor, Pennsylvania 19038, USA. ⁶Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia.</footnote>

<footnote>* These authors contributed equally to this work. Correspondence and requests for materials should be addressed to H.J.G. (email: harry.gilbert@ncl.ac.uk) or to D.N.B. (email: david.bolam@ncl.ac.uk).</footnote>
Genomic and metagenomic data show that the catabolism of complex dietary carbohydrates (CDCs), also known as dietary fibre, play a central role in defining the structure of the human gut microbiota (HGM)¹,². The Bacteroidetes of the HGM, which include the genus *Bacteroides*, are viewed as glycan-degrading generalists capable of utilizing an extensive array of both dietary (plant and animal) and host complex carbohydrates³,⁴. The other main phyla of the gut, Firmicutes and Actinobacteria (particularly *Bifidobacterium* spp.), also have the capacity to degrade glycans, but appear more specialized with a preference for plant storage polysaccharides (starch and fructan) and oligosaccharides, such as those derived from arabinoxylans (AXs)⁵,⁶. Current strategies deployed to manipulate the composition of the HGM to maximize human health are centred on the use of complex carbohydrates (prebiotics)⁷ to encourage the expansion of beneficial (probiotic) microorganisms⁸. The effects of these strategies are ultimately dependent on the glycan food web in the HGM⁹,¹⁰. However, the extent to which members of the HGM act in synergy to degrade and utilize CDCs is currently unclear. Furthermore, while genomic and metagenomic sequencing has identified an extensive array of putative glycan degrading enzymes encoded by the HGM¹,¹¹,¹², the lack of functional signposts, obtained through enzyme characterization, greatly restricts our capacity to identify the activities displayed by these candidate glycanases and thus predict the glycan preferences of specific members of the HGM.

Xylan, a major component of many plant cell walls, is one of the most variable plant structural polysaccharides, and is a common component of cereal-derived human foods¹³⁻¹⁵. The

main classes of xylan are the glucuronoxylans (GXs), arabinoxylans (AXs) and glucuronoarabinoxylans (GAXs). The designation of these three polysaccharides reflects the nature of the monosaccharides that decorate the conserved β1,4-xylose (Xyl) backbone (Fig. 1). Many GXs and AXs are relatively simple. In contrast, the number and complexity of the side chains on cereal GAXs, exemplified by corn bran xylan (CX) (Fig. 1) explains why these glycans are recalcitrant to breakdown by known xylan-degrading enzyme systems¹³.

*Bacteroides ovatus*, a prominent member of the HGM, was previously shown to utilize an extensive range of plant polysaccharides, including xylans⁴. Here we explored the mechanisms by which this bacterium degrades the different forms of the polysaccharide presented to the HGM. Our data show that *B. ovatus* has evolved a highly dynamic xylan degrading system that can respond to the different forms of the xylose polymer. The bacterium has recruited several enzymes from glycoside hydrolase (GH) families¹⁶ not previously associated with xylan degradation. We also show that, although polysaccharide breakdown products (PBPs) are released during degradation of both simple and complex xylans by *B. ovatus*, the ability of other members of the HGM to use these PBPs is dictated by the complexity of the target xylan.

## Results

**B. ovatus** contains two discrete loci activated by xylans. The genomes of Bacteroidetes contain numerous polysaccharide utilization loci (PULs) that are optimized to orchestrate the degradation of specific polysaccharides⁴,¹⁷⁻¹⁹. The glycans

<figure><img src="image_4.png" /></figure>

**Figure 1** | Schematic of the structures of the main classes of xylan. The monosaccharides and linkages in the main classes of xylan are shown and are represented in their Consortium for Functional Glycomics format⁶⁹. The xylans used in this study were from birchwood (birch glucuronoxylan; BGX), wheat flour (wheat arabinoxylan; WAX) and corn bran (corn glucuronoarabinoxylan; CX).
targeted by the PULs activate transcription of their cognate loci. *Bacteroides ovatus* strain ATCC 8483 contains two PULs, spanning locus tags *bacova_03417-50* (large xylan PUL or PUL-XylL) and *bacova_04385-94* (small xylan PUL or PUL-XylS), that are activated when the organism is grown on wheat arabinoxylan⁴ (WAX) (Fig. 2a,b). In contrast, linear xylooligosaccharides (preferentially xylotetraose) induce transcription of PUL-XylS alone⁴. Here we show that growth on CX leads to significant activation of PUL-XylL and only very low activation of PUL-XylS, while glucuronoxylan (from birch; BGX) activates only PUL-XylS (Fig. 2c). These data suggest that the degradative system encoded by PUL-XylS targets simple linear or glucurono-substituted xylans, while the products of PUL-XylL drive the metabolism of more complex decorated forms of the hemicellulose. Bacteroidetes spp. are known to express enzymes capable of breaking down simple forms of xylan²⁰,²¹; however, the mechanisms by which these microorganisms utilize complex cereal GAXs (for example, CX) are not known. Both PULs encode glycanases belonging to GH families¹⁶ that are known to contribute to the degradation of simple xylans (GH3, GH10, GH30, GH43, GH67 and GH115)²² (Fig. 2a,b). Significantly, PUL-XylL also directs the synthesis of enzymes belonging to families (GH31, GH95, GH97 and GH98)⁴ not previously implicated in the deconstruction of the hemicellulose (Fig. 2a).

**Surface acting xylan-degrading enzymes.** Bioinformatic analysis, in conjunction with cell localization studies and whole-cell assays, was used to identify the enzymes encoded by PUL-XylS and PUL-XylL that are displayed on the bacterial surface. The presence of a

predicted Type II lipoprotein signal sequence or uncleaved transmembrane (TM) anchor in seven GH family proteins indicated that they were membrane attached. The cell localization of these seven proteins was evaluated experimentally by (i) fluorescence microscopy of whole cells treated with labelled antibodies directed against the appropriate proteins and (ii) western blotting after removal of surface proteins with Proteinase K. The data indicated that BACOVA_03419 (GH3), BACOVA_03421 (GH43), BACOVA_03431 (inactive GH10, see later results), BACOVA_03432 (GH30), BACOVA_03433 (GH98) and BACOVA_04390 (GH10) were presented on the surface of *B. ovatus* (Fig. 3). The only lipoprotein facing the periplasm was the GH43 enzyme BACOVA_03425. These conclusions were further supported by whole-cell assays of the *B. ovatus* mutant lacking PUL-XylS (ΔPUL-XylS) incubated with BGX. The product profile of ΔPUL-XylS was extremely similar to that of recombinant BACOVA_03432 GH30 against the polysaccharide (Supplementary Fig. 1a). Thus, BGX degradation is mediated predominantly by BACOVA_03432 in the ΔPUL-XylS mutant, confirming the extracellular location of the GH30 enzyme. Whole-cell assays also showed that the PUL-XylL encoded GH43 arabinofuranosidase, BACOVA_03421 (see results below), was extracellular, as arabinose was released from WAX by the ΔPUL-XylS mutant (Supplementary Fig. 1b). It is surprising that BACOVA_03419, a GH3 β-xylosidase, is located on the surface as its obvious role is to contribute to the degradation of the xylooligosaccharides generated in the periplasm (Supplementary Fig. 2a). BACOVA_03419, however, displays only low activity against linear β1,4-xylooligosaccharides, but, unlike the cytoplasmic xylosidase, BACOVA_04386 (see results below), will

<figure><img src="image_5.png" /></figure>

**Figure 2 | Schematic of the *B. ovatus* xylan PULs and differential expression during growth on different xylans. (a,b) schematic of PUL-XylL and PUL-XylS, respectively. Each gene is drawn to scale as a rectangle with its orientation indicated by the arrow head. The numbers below each gene is its locus tag (*bacova_XXXXX*). Genes encoding known or predicted functionalities are colour-coded and, where appropriate, are also annotated according to their CAZy family number: glycoside hydrolase (GH; green), carbohydrate esterase (CE; purple), carbohydrate-binding module (CBM). Surface located proteins are marked with an asterisk. SGBP = surface glycan binding proteins are coloured orange, or, if also in a CAZy family (BACOVA_03431; inactive GH10), are coloured half orange, half green. UNK = unknown (purple), but distant similarity to CE6 carbohydrate esterases. HTCS = hybrid two component system (light or dark blue). MFS = transporter of the major facilitator superfamily (pink). Grey = unknown function (note, there is a structure of BACOVA_03430, PDB accession code 3N91). SusC-like (yellow) and SusD-like (light tan) proteins are a defining feature of PULs and are responsible for import of complex glycans across the outer membrane¹⁸. SusC-like proteins are TonB-dependent transporters, while SusD-like proteins are surface lipoproteins that likely function to deliver the target glycan to their partner SusC. (c) Cells were grown on minimal media with polysaccharide as the sole carbon source, and levels of different *susC* transcripts (locus tags shown; used as a proxy for expression of the whole PUL) from each PUL were quantified by qRT-PCR. The y-axis shows the Log fold-change relative to a minimal media-glucose reference; x-axis labels indicate the xylans used. WAX = wheat arabinoxylan, CX = corn xylan, BGX = birchwood glucuronoxylan. CX is a highly complex xylan compared to WAX and BGX. All data were analysed by one-way ANOVA followed by Tukey's multiple comparison test (*** = *P* ≤ 0.001).
<figure><img src="image_6.png" /><figcaption>**Figure 3 | Cellular location of Xylan PUL enzymes containing a predicted N-terminal lipidation site or uncleaved transmembrane anchor.**</figcaption></figure>

(a) Fluorescent microscopic images of *B. ovatus* cells cultured on WAX and incubated with polyclonal antibodies raised against recombinant BACOVA_03419 (GH3), BACOVA_03421 and BACOVA_03425 (GH43s), BACOVA_03431 (SGBP), BACOVA_03432 (GH30), BACOVA_03433 (GH98) and BACOVA_04390 (GH10). (b) Western blots of the *B. ovatus* cells shown in panel a, either untreated with Proteinase K (-) or incubated with 2 mg ml⁻¹ Proteinase K for 16 h (+). The blots were probed with antibodies against the *B. ovatus* proteins indicated. The full blots are shown in Supplementary Fig. 8.

cleave linkages in which the leaving group xylose is decorated (Supplementary Table 1 and Supplementary Fig. 2b,c). While the roles of BACOVA_03419 and BACOVA_03421 in xylan degradation are unclear, we speculate that the enzymes contribute to the surface 'trimming' of complex xylans.

Three of the surface enzymes (BACOVA_04390 GH10, BACOVA_03432 GH30 and BACOVA_03433 GH98 (Fig. 3a,b) exhibited endo activity and therefore are likely involved in cleaving the xylan backbone before importing into the periplasm. Studies on recombinant BACOVA_04390 GH10 showed that the endo-xylanase displayed a catalytic efficiency against xylan ~30-fold lower than typical GH10 enzymes²³,²⁴ (Supplementary Table 2). The limit products against sparsely decorated xylan contains oligosaccharides with a degree of polymerization (d.p.) >4 (Supplementary Fig. 3a). Typical GH10 xylanases²³,²⁴, exemplified by the periplasmic enzyme BACOVA_04387 generate xylose and xylobiose as the final products (Supplementary Fig. 3b). These data suggest that BACOVA_04390 contains an extended substrate binding cleft that requires a high degree of occupancy for catalysis to occur. Substrate depletion data using xylooligosaccharides confirmed that the enzyme had very poor activity against substrates with a d.p. ≤5 compared with typical GH10 enzymes (Supplementary Table 2). BACOVA_04390 has an unusual primary sequence in which the canonical GH10 (β/α)₈-fold barrel is interrupted by two copies of a non-catalytic carbohydrate-binding module (CBM) that is distantly related to CBM family 4, suggesting that the two CBMs contribute to an extended substrate binding cleft²¹,²⁵. Indeed a recent study has shown that the internal CBMs from a related *B. intestinalis* GH10 bind xylans and long xylooligosaccharides²¹. This view is consistent with the observation that endo-acting enzymes, which specifically target large polymers, contain additional subsite(s) likely compensating for relatively weak binding at subsites proximal to the active site²⁶,²⁷.

The GH30 enzyme BACOVA_03432 produced a range of large oligosaccharides from BGX (Supplementary Fig. 1a and Supplementary Table 2), but was not active on WAX or CX. This specificity is typical of GH30 glucuronoxylanases and is driven by recognition of the (Me)GlcA side chain at the -2 subsite of these enzymes²⁸. While the lack of activity against CX is somewhat surprising as the enzyme is part of PUL-XylL, it may be that less heavily decorated GAXs are the target substrates for this enzyme.

The surface GH98 enzyme BACOVA_03433 also generates very large oligosaccharides but exclusively from CX (Fig. 4a,b). As BACOVA_03433 is a member of GH98, a family that hitherto contained only endo-β-D-galactosidases²⁹, it was formally possible that the enzyme targets rare β-D-Gal residues in the backbone of CX. Converting the reducing end aldehyde of the products generated by BACOVA_03433 to the alditol followed by acid hydrolysis revealed the presence of xylitol but not galactitol (Fig. 4c). These data show that BACOVA_03433 is an endo-xylanase that hydrolyses xylosidic bonds in a context found only in very complex xylans. It is possible that the enzyme targets an unusual linkage in the xylan backbone; however, this is unlikely as the enzyme was not active against β1,4-β1,3-mixed linked xylans and CX is reported to contain a backbone comprising exclusively β1,4-Xyl units¹³. To explore the specificity determinants of the GH98 enzyme further, CX was pre-treated with cocktails of side-chain cleaving hydrolyses (arabinofuranosidases, xylosidases, glucuronidases and galactosidases) encoded by PUL-XylL, and the debranched forms of the xylan were then incubated with the GH98 xylanase. Only GH consortia lacking BACOVA_03417 generated a variant of CX that was hydrolysed by the GH98 enzyme (Fig. 4d). As described below BACOVA_03417 is an α-L-arabinofuranosidase that targets the O3 linked Araf on backbone Xyl residues decorated at both O2 and O3 with α-L-Araf. It would appear, therefore, that the double Araf decoration is a component of the specificity determinant for BACOVA_03433. Wheat arabinoxylan, however, also contains this double arabinose decoration and is not a substrate for the GH98 xylanase (Fig. 4a). Thus, the specificity of BACOVA_03433 is more complex than simply Xyl units linked O2 and O3 to arabinose. Interestingly, the CX backbone contains low levels of doubly substituted residues with Xylp as the O3 linked side chain (see Fig. 1 and results below). This structure is also a substrate for BACOVA_03417 (Supplementary Fig. 4a), suggesting Xylp at O3 (and Araf at O2) may be the critical specificity determinant of the GH98 enzyme. These Xylp decorations are absent in WAX (Supplementary Fig. 4a-c), providing a possible explanation for the lack of activity of the GH98 enzyme against this substrate.

The large products generated by all three surface xylanases indicate that these enzymes recognize infrequent specificity determinants in the xylan, which, as discussed below, likely
<figure><img src="image_7.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_10.png" /></figure>

**Figure 4** | The activity of the GH98 enzyme BACOVA_03433. (a) Different xylans at 1% (w/v) were incubated with 1 µM of the GH98 enzyme in PBS for 16 h at 37 °C. The reaction products were subjected to thin layer chromatography and compared with xylooligosaccharide standards with a d.p. ranging from 1 to 5 (X1-X5). (b) HPAEC-PAD trace showing the high molecular weight oligosaccharide products of CX digestion by the GH98 endo-xylanase, separated using a sodium acetate gradient of 0-500 mM. The reaction conditions used are the same as in panel a. Control reactions of substrate and enzyme only are shown. (c) CX was incubated with 5 µM of the GH98 enzyme for 24 h in 50 mM sodium phosphate buffer, pH 7.0, at 37 °C. The reaction products, and CX that was not incubated with the GH98 enzyme, were then treated with sodium borohydride, to convert the reducing end monosaccharide unit into its alditol, followed by acid hydrolysis to convert the oligosaccharides (or in the case of control CX alone) into their monosaccharide constituents. The hydrolysed products were subjected to HPAEC and compared to the migration of standard monosaccharides as well as the alditol of xylose (xylitol) and galactose (galactitol). (d) The GH98 enzyme was incubated with CX (lane CX vs GH98), or CX that had been pre-treated with all the side chain cleaving enzymes encoded by PUL-XylL (lane CX vs MIX), or CX that had been pre-treated with all the side chain cleaving enzymes encoded by PUL-XylL, except the GH43 enzyme BACOVA_03417 that removes O3 linked arabinose from double substituted xylose (lane CX vs MIX no double abfase).

contributes to the capacity of *B. ovatus* to efficiently utilize the oligosaccharides generated by these GHs. It is particularly intriguing that the specificities of the three endo-acting enzymes are adapted to recognize the variable decorations appended to the backbone of different xylans, likely reflecting the multiple forms of the polysaccharide presented to the human large intestine.

**Surface xylan-binding proteins.** PUL-XylS and PUL-XylL encode one (BACOVA_04392) and two (BACOVA_03427 and BACOVA_03429) SusD-like proteins, respectively. SusD-like proteins appear to be involved in binding and delivering oligosaccharides to their cognate SusC-like transporter for import into the periplasm⁴,¹⁸,¹⁹,³⁰. Isothermal titration calorimetry showed that BACOVA_04392 bound to the relatively simple xylans WAX

and BGX with moderate affinity ($K_d$ ~ 0.3–0.4 mM), but did not interact with highly decorated CX (Supplementary Table 3). The protein bound xylohexaose, but not shorter xyloooligosaccharides, consistent with the long products generated by the surface GH10 xylanase encoded by PUL-XylS, BACOVA_04390 (Supplementary Fig. 3a). The SusD-like protein, derived from PUL-XylL, BACOVA_03427, bound tightly to BGX and CX ($K_d$ ~ 7–12 µM), but not WAX or linear xyloooligosaccharides, although the precise target ligand for this protein is unclear (Supplementary Table 3). The other PUL-XylL SusD-like protein, BACOVA_03429, could not be expressed in a recombinant form.

In addition to the SusD-like proteins, PULs often encode surface-located glycan binding proteins (SGBPs), which are functionally analogous to SusE and F from the prototypic starch utilization system¹⁷⁻¹⁹,³¹. The capacity of three non-catalytic
lipoproteins encoded by the xylan PULs (candidate SGBPs) to bind xylans and undecorated xylooligosaccharides was thus also explored. The data show that BACOVA_04391, from PUL-XylS, bound xylans and xylooligosaccharides with affinity increasing with chain length (Supplementary Table 3). The SGBP displayed similar specificity to the PUL-XylS SusD-like protein, albeit with higher affinity for all ligands. The crystal structure of BACOVA_04391 (PDB accession code 3ORJ) revealed a protein that contains three β-sandwich domains that are aligned in a linear orientation with respect to each other (Supplementary Fig. 5). The loops that connect the two β-sheets of the C-terminal domain form a wide cleft in which two aromatic residues, Trp241 and Tyr404, are orientated parallel to each other and are in an ideal position to make hydrophobic interactions with the two faces of a xylopyranose ring. The topology of this site centred on the parallel aromatic residues is very similar to non-catalytic CBMs that also target xylan²⁵. The relative positioning of the domains in the lipoprotein BACOVA_04391 positions the ligand binding site away from the cell surface, likely to maximize efficiency of initial glycan capture.

PUL-XylL encodes two potential SGBPs, BACOVA_03430 and BACOVA_03431. BACOVA_03430 did not bind to any of the xylans or xylooligosaccharides evaluated and thus is unlikely to be a SGBP. BACOVA_03431, a surface located lipoprotein (Fig. 2), has sequence similarity to GH10 enzymes but was inactive against xylans, reflecting the loss of the catalytic nucleophile and disruption of the -2 subsite³². This inactive GH10 was shown to be a xylan binding protein recognizing both BGX and CX. BACOVA_03431 displayed similar specificity (but lower affinity) to BACOVA_03427, one of the SusD homologues encoded by PUL-XylL. Although no other naturally occurring 'inactive' GH10 members have been described, the function of BACOVA_03431 has some resonance with GH18 chitin binding lectins that lack the catalytic residues present in chitinases within this family³³. Collectively, the data indicate that BACOVA_03431 plays a role in xylan acquisition at the cell surface.

It would appear, therefore, that the binding specificities of the SGBPs and SusD-like proteins are entirely consistent with the xylans targeted by the enzymes encoded by the two PULs. Thus, the PUL-XylS SusD-like protein and SGBP target GXs and long chain undecorated xylooligosaccharides, whereas the equivalent PUL-XylL proteins recognize the more complex structures found in cereal xylans such as CX.

**Enzymatic debranching of xylans in the periplasm of B. ovatus.**
The xylan-derived xylooligosaccharides entering the periplasm are degraded by a suite of primarily exo-acting enzymes encoded by PUL-XylS and PUL-XylL, which are described below.

BACOVA_03417 GH43 generated arabinose from AXs, CX and sugar beet arabinan, before and after treatment with an arabinofuranosidase that removes only single Araf decorations (CjAbf51A³⁴; Abfase-S), (Supplementary Fig. 4a,b and Supplementary Table 4), while the enzyme displayed no activity against WAX pre-treated with Abfase-D (*HiAXHd3*³⁵), a fungal enzyme that targets the O3 arabinose of Xyl units decorated at both O2 and O3 with Araf (Supplementary Fig. 4c). BACOVA_03417 shows ~30% sequence identity with *HiAXHd3*, with all the residues that contribute to specificity in the active site and the O2 Araf specificity pocket conserved in the two enzymes. It is likely, therefore, that the *B. ovatus* enzyme targets the O3 linkage in the double substitutions, using the O2 Araf as a critical specificity determinant as occurs in *HiAXHd3*³⁵.

Incubation of BACOVA_03421 and BACOVA_03425 GH43s with WAX or CX released arabinose (Supplementary Fig. 4d,e and Supplementary Table 4), while neither enzyme was active

against WAX pre-treated with the mono-specific GH51 arabino-
furanosidase CjAbf51A³⁴ (Supplementary Fig. 4f). These data
show that both BACOVA_03421 and BACOVA_03425 are
arabinofuranosidases that target single substitutions (either O2
or O3 linked) on the xylan backbone.

It should be noted that the three GH43 arabinofuranosidases (BACOVA_03417, BACOVA_03424 and BACOVA_03425) also released a small amount of xylose from CX but not from WAX (Supplementary Fig. 4), reflecting β1,3-Xyl side chains in CX (as single decorations and as a component of an O2 Araf and O3 Xylp double decoration, Fig. 1). These data are consistent with the observation that the active sites of GH43 enzymes bind Xylp and Araf³⁵.

PUL-XylL encoded two other GH43 enzymes, BACOVA_03424 and BACOVA_03436, that contained a C-terminal family 6 CBM¹⁶,³⁶ and displayed ~ 46% sequence identity with each other. Both enzymes displayed limited activity against CX only, generating very low concentrations of a short oligosaccharide. The functional significance of these two enzymes is unclear.

BACOVA_03449 is a GH115 α-glucuronidase from PUL-XylL that was previously shown to cleave internally substituted glucuronic and methyl-glucuronic acid from GXs, GAXs and glucurono-xylooligosaccharides³⁷. These data are consistent with the predicted periplasmic location of BACOVA_03449. This supports the view that the target substrates for the GH115 α-glucuronidase are internally substituted oligosaccharide products generated by the surface-located GH30 or GH98 endo-xylanases, which are also from PUL-XylL. BACOVA_03485 from PUL-XylS, is a GH67 α-glucuronidase and, like other members of this family, is active only on terminally substituted xylo-oligosaccharides (Supplementary Table 5)³⁸. This is consistent with the products generated by the surface GH10 endo-xylanase, BACOVA_04390, which is also encoded by PUL-XylS.

**Novel xylan debranching enzymes.** As discussed above the two xylan PULs encode debranching enzymes located in GH families typically associated with xylan degradation (GH43, GH67 and GH115) that cleave α-L-arabinofuranosyl and α-D-glucuronyl linkages. In addition, PUL-XylL encodes three periplasmic proteins that belong to GH31, GH97 and GH95, families that have not previously been implicated in the deconstruction of the polysaccharide. Assays showed that the GH31 enzyme, BACOVA_03422, released xylose from CX but displayed no activity against a range of other xylans (Supplementary Fig. 6). As GH31 enzymes target α-D-glycosides³⁹, BACOVA_03422 is likely an α-xylosidase, consistent with the presence of α-D-Xyl in CX¹³. BACOVA_03422 is ~3,000-fold more active against CX than a xyloglucan oligosaccharide (Supplementary Table 6), demonstrating that the enzyme is a xylan-specific α-xylosidase. In contrast, a previously characterized GH31 α-xylosidase, CjXyl31A⁴⁰, exhibited minimal activity against CX, strongly preferring xyloglucan oligosaccharide (Supplementary Table 6). These data show that the xylan specificity displayed by the BACOVA_03422 enzyme is not a general property of GH31 enzymes. The GH97 enzyme, BACOVA_03423, was active against 4-nitrophenyl-α-D-galactopyranoside, consistent with the presence of α-D-Gal in some highly complex xylans¹³. The GH95 enzyme BACOVA_03438, although located in an 'α-L-fucosidase specific family'⁴¹, released exclusively L-Gal from CX, displaying only trace activity against oligosaccharides containing terminal α-L-Fuc residues (Fig. 5a-c and Supplementary Table 7). Conversely, a previously characterized GH95 α-1,2-L-fucosidase (AfcA)⁴² exhibited minimal α-L-galactosidase activity against CX (Supplementary Table 7), showing that GH95 enzymes can be
<figure><img src="image_11.png" /><figcaption>**Figure 5** | Activity and structure of BACOVA_03438 GH95 α-L-galactosidase and comparison to the *Bifidobacterium bifidum* GH95 α-L-fucosidase. (a,b) TLC and HPAEC profiles, respectively, of the reaction product (L-galactose) generated by incubating CX at 1 mg ml⁻¹ with 1 µM BACOVA_03438 for 16 h in 20 mM sodium phosphate buffer, pH 7.0, at 37 °C. In (a) the lanes are as follows: (–) CX no enzyme; (+) CX incubated with BACOVA_03438. The asterisk marks the position of the single reaction product released from CX by the GH95. (c) The product of the GH95 is resistant to oxidation by D-galactose dehydrogenase (GDH). (d) Schematic of BACOVA_03438, colour-ramped from blue (N-terminus) to red (C-terminus), with L-Gal bound in the active site (carbohydrate carbons shown as green sticks). (e) Solvent exposed surface representation of *B. ovatus* GH95. The active site pocket that houses β-L-Gal product is boxed. (f) Side chains of the active site residues of BACOVA_03438 (magenta carbons) and the GH95 α-L-fucosidase (blue carbons; PDB: 2EAE) from *Bifidobacterium bifidum*, bound to L-Gal (green carbons) and L-Fuc (orange carbons), respectively. Predicted hydrogen bonds between the amino-acid side chains and carbohydrates are shown as yellow dotted lines, except for the H-bond between the Oγ of BACOVA_03438 Thr370 and O6 of L-Gal, which is shown as a red dotted line. (g) Overlay of the active site residues displayed in panel (f) with the side-chains of the Thr/His polymorphism that may play a role in specificity for L-Gal over L-Fuc labelled. A stereo image of a portion of the electron density map is shown in Supplementary Fig. 9.</figcaption></figure>

either L-Gal or L-Fuc specific. Although these data provide an initial example of a glycosidase specific for L-galactose, it should be noted that the GH117 family contains enzymes that target the 3,6-anhydro form of the sugar⁴³.

**Crystal structure of the GH95 α-L-galactosidase.** The crystal structure of the α-L-galactosidase, BACOVA_03438, displays a

central (α/α)₆-helical barrel domain flanked by β-sandwich domains at the N- and C-termini, respectively. In the structure of the protein-product complex, β-L-Gal (GH95s are inverting enzymes), which adopts a ¹C₄ geometry typical of L-mono-saccharides, is buried in the central cavity of the helical barrel domain, (Fig. 5d–g). Structural analyses in combination with mutational experiments revealed that the highly conserved Glu501 acts as the general acid catalyst in a single displacement
(inverting) mechanism (Fig. 5f and Supplementary Table 7). No carboxylic acid residue was found at the appropriate position for a general base catalyst, similar to the *Bifidobacterium bifidum* GH95 α-L-fucosidase, AfcA. It was proposed that in AfcA two asparagine residues activate a water molecule that is suitably located to perform a nucleophilic attack on the C1 atom of L-fucose in the active site⁴², which is in sharp contrast with typical GHs where Asp or Glu act as the catalytic base⁴⁴. The proposed catalytic asparagines in AfcA (Asn421 and Asn423) are structurally conserved in the GH95 α-L-galactosidase (Asn372 and Asn374), suggesting a similar role in BACOVA_03438 (Fig. 5f,g). Almost all of the residues that interact with the pyranose ring of L-Gal in BACOVA_03438 are conserved in the AfcA fucosidase, although, intriguingly, Thr370 in the L-galactosidase is replaced with His419 in the fucosidase (Fig. 5f,g). In AfcA His419 contributes to the hydrophobic sheath in which the C6 methyl group of the Fuc is housed, while in BACOVA_03438 the γO of Thr370 makes a direct polar contact with the O6 of L-Gal, the only part of the monosaccharide that differs from L-Fuc. The importance of the Thr370 in BACOVA_03438 is demonstrated by the complete loss of activity of the T370H mutant against CX (Supplementary Table 7). While the Thr/His polymorphism may be a specificity determinant for galactosidase and fucosidase activity within GH95, the paucity of enzymes characterized in this family prevents a definitive conclusion regarding its functional significance. Indeed, the observation that the T370H mutation does not increase fucosidase activity suggests that L-Gal to L-Fuc specificity is driven by more than just the –1 site (Supplementary Table 7). Significantly, the leaving group subsite (+1 subsite) is not conserved in the two enzymes (inferred in the case of BACOVA_03438), and it is possible that substrate binding in this region contributes to enzyme specificity. A key role for the +1 subsite in the two enzymes is consistent with their lack of activity against 4-nitrophenyl-α-L-fucopyranoside (AfcA) and 4-nitrophenyl-α-L-galactopyranoside (BACOVA_03438).

**Hydrolysis of undecorated xylooligosaccharides.** Breakdown of debranched xylooligosaccharides, derived through the action of

the exo-acting enzymes described above, is mediated by enzymes encoded by PUL-XylS. A typical GH10 endo-acting xylanase (BACOVA_04387) rapidly generates oligosaccharides in the periplasm with a d.p. of primarily 1 to 4 (Supplementary Fig. 3b and Supplementary Table 2). These small xylooligosaccharides are substrates for the cytoplasmic GH43 β-xylosidase, BACOVA_04386 (Supplementary Table 1 and Supplementary Fig. 2c). The spatial separation of the β-xylosidase from the enzymes that generate its substrate may reflect a selection advantage as less energy would be required to transport a small number of xylooligosaccharides across the inner membrane (likely via the PUL-encoded MFS class transporter, BACOVA_04388, Fig. 2b), compared with a larger number of xylose units if the β-xylosidase was periplasmic. It is also possible that the different cellular locations of the xylosidase and xylanase maximizes the signal concentration for the PUL-XylS HTCS sensor-regulator (BACOVA_04394), which recognizes short xylooligosaccharides⁴.

**The function of the two xylan PULs in B. ovatus.** The role of PUL-XylS and PUL-XylL in xylan degradation is illustrated by the growth profiles of mutant strains. Deletion of PUL-XylS (ΔPUL-XylS) did not influence the growth of *B. ovatus* on complex xylans (GAXs) from corn, rice and sorghum, but prevented growth on simpler xylans (BGX or WAX) (Fig. 6). By contrast, inactivation of PUL-XylL (ΔPUL-XylL) prevented growth on the three complex GAXs but did not significantly affect the utilization of simpler BGX or WAX (Fig. 6). The somewhat higher final density observed for the wild-type bacterium on WAX vs ΔPUL-XylL is likely due to the inability of this mutant to remove the arabinose decorations from this substrate. The growth profile of the mutants are consistent with the biochemical data presented above showing that PUL-XylL encodes the only endo-acting enzyme (GH98) able to cleave the backbone of CX (and likely other complex GAXs), a prerequisite for the utilization of these polysaccharides (see below). Similarly, the inability of the PUL-XylS mutant to utilize simple xylans such as GXs and AXs, where the xylan backbone is the major component, is consistent with the locus encoding the key endo-acting surface GH10 (BACOVA_04390) that hydrolyses the backbone of sparsely

<figure><img src="image_12.png" /></figure>

**Figure 6 | Growth profile of wild type and mutants of B. ovatus on xylans.** WT and PUL deletion mutants were cultured in minimal media containing different xylans at 0.5% (w/v) as the sole carbon source. Birchwood glucuronoxylan (BGX; (a)) and wheat arabinoxylan (WAX; (b)) are relatively simple, sparsely decorated structures, while corn (c), sorghum (d) and rice (e) xylans are more complex heavily decorated glucuronoarabinoxylans (GAXs). The high starting OD for BGX is due to the background turbidity of the polysaccharide. A total of 6 replicate cultures were monitored at 20 min intervals for each substrate and used to generate the average growth curves shown.
decorated xylans. It is somewhat surprising, however, that ΔPUL-XylS is able to efficiently utilize GAX even though the mutant is unable to fully depolymerize the xylan backbone. This likely reflects the extensive side chains present in CX and other complex GAXs, where the xylose backbone constitutes only a relatively minor component of the total monosaccharide units in these polysaccharides.

To summarize, using the expression, localization and bio-
chemical data presented above Fig. 7 shows the proposed models
for the acquisition and breakdown by *B. ovatus* of the different
classes of xylan presented to the colonic microbiota.

The conservation of PUL-XylS and PUL-XylL in Bacteroidetes.

To identify PULs similar to PUL-XylL and -XylS in other genomes, we examined not only the gene family composition but also gene synteny in 260 Bacteroidetes species (Supplementary Table 8). Comparative genomics revealed a high degree of conservation of PUL-XylL in *B. ovatus* and *B. xylanisolvens* strains. Indeed, PULs with a high level of synteny with PUL-XylL was observed in 22 of the 23 strains examined. In these PULs, 17 genes show strict conservation with those present in PUL-XylL, which are arranged into three syntenic blocks (Fig. 8a). These 22 PUL-XylL homologues can be subdivided into two groups that

a

<figure><img src="image_13.png" /></figure>

**Figure 7** | Models for the degradation of different forms of xylan by *B. ovatus*. In the upper part of each panel the monosaccharides and linkages in the main classes of xylan are shown and are represented in their Consortium for Functional Glycomics format⁶⁹. The lower parts of each panel show the model for degradation of complex GAXs like CX (a) AXs like WAX (b) and GXs like BGX (c) The red asterisk next to BACOVA_03425 GH43 indicates that the single specific arabinofuranosidase acts after BACOVA_03417 (α1,3 double-specific arabinofuranosidase) to remove the remaining α1,2 linked L-arabinose. The black arrows indicate examples of the linkages cleaved by the enzymes, shown in green ovals and identified by their glycoside hydrolase family and locus tag. In panel a, the dotted arrows from the GH98 enzyme identify possibly linkages hydrolysed by the xylanase. Although not shown on the CX model, the PUL-XylS apparatus will be present, albeit at low levels (see Fig. 2c). A schematic of the crystal structures, colour-ramped from blue (N-terminus) to red (C-terminus), of two of the xylan degrading enzymes, determined in this study (BACOVA_03438 GH95; PDB accession code 4UFC) and in a recent report³⁷ (BACOVA_03449 GH115; PDB accession code 4C91) are shown. The red arrow heads indicate glycan transport between cellular locations, while the yellow arrow heads show the ligand (or potential ligand in case of the question mark) that binds and activates the hybrid two component system (HTCS) encoded by PUL-XylS (BACOVA_04394, light blue HTCS; activating ligand xylotetraose⁴) or PUL-XylL (BACOVA_03437, dark blue HTCS; activating ligand likely an arabino-xylooligosaccharide). MFS is a member of the Major Facilitator Superfamily of transport proteins, and SGBP signifies the two surface xylan binding proteins encoded by the xylan PULs, BACOVA_03431 (inactive GH10, PUL-XylL) and BACOVA_04391 (PUL-XylS). A schematic of the crystal structure of BACOVA_04391 (PDB accession code 3ORJ), colour-ramped from blue (N-terminus) to red (C-terminus), is shown. Note only one of the two SusC/D pairs from PUL-XylL is shown for clarity.

b

<figure><img src="image_14.png" /></figure>

C

<figure><img src="image_15.png" /></figure>
a

<figure><img src="image_16.png" /></figure>

b

<figure><img src="image_17.png" /></figure>

**Figure 8** | Conservation of the *B. ovatus* ATCC 8483 xylan PULs in other *Bacteroides* species. CAZymes (green, CBM modules were omitted for clarity), susC-like (yellow), susD-like (orange) and the HTCS sensor-regulator (blue) genes are colour-coded. Genes encoding proteins of unknown function or non-CAZymes (for example, surface glycan binding protein and MFS carbohydrate importer in PUL-XylS) are shown in grey. (a) PUL-XylL. Vertical dotted lines indicate split gene models, while missing gene models are shown in red. Microsyntenic segments are depicted by grey connectors, and include adjacent conserved non-PUL members (shown with dotted lines). Eighteen genes are conserved across all strains, and 17 of these form three highly conserved blocks indicated on top. The major difference between the two groups of the GH98 for GH5_21 swap is indicated. Inset shows a TLC of the activity of the recombinant *B. xylanisolvens* XBA1 GH5_21 enzyme (BXY_29320; boxed in red) against different xylans (1 µM final enzyme vs 0.5% xylan in 20 mM sodium phosphate buffer, pH 7.0, 16 h at 37 °C). The *B. xylanisolvens* GH5 enzyme displays endo-xylanase activity against CX, producing a smear of large oligosaccharide products similar to the product profile observed for the *B. ovatus* GH98 enzyme, BACOVA_03438. The GH5 xylanase also hydrolyses WAX, unlike the GH98 enzyme, but displays no activity against BGX. X1-X5 indicates xylooligosaccharide standards. (b) PUL-XylS. The lower part of the panel shows the strains that contain only the second half of PUL-XylS that encodes the intracellular (periplasmic and cytoplasmic) enzymes only.

differ by the presence of a tandem array of three genes. The three genes in group one PUL-XylLs, which encode GH30_8 (subfamily 8)⁴⁵, GH98 and GH115 enzymes (including *B. ovatus* ATCC 8483), are replaced in group two PUL-XylLs by a single gene encoding a GH5 subfamily 21 (GH5_21)⁴⁶ glycanase (Fig. 8a). Significantly, the GH5_21 enzyme from *B. xylanisolvens* XBA1 displays activity on CX (Fig. 8a inset), similar to the *B. ovatus* GH98 enzyme BACOVA_03433. This

suggests that both forms of PUL-XylL retain the capacity to hydrolyse the backbone of complex GAXs without the need for extensive debranching. The high synteny described above for PUL-XylL in the *B. ovatus* and *B. xylanisolvens* strains was not found in any of the other 260 Bacteroidetes genomes analysed. In particular, the distribution of the GH98-encoding gene is extremely narrow, being present in only five of the 260 genomes analysed.
With respect to PUL-XylS there is perfect conservation of the synteny of the genes encoding the GH67 α-glucuronidase, GH43 β-xylosidase and periplasmic GH10 xylanase in all 23 *B. ovatus* and *B. xylanisolvens* strains (Fig. 8b). In 11 of these strains, however, the genes encoding the surface GH10 xylanase, SusD-like, SusC-like and HTCS proteins were absent. The loss of these critical components in these strains is likely to render the PULs non-functional. No other candidate PUL with complete conservation of all PUL-XylS components could be found in the other Bacteroidetes that we examined. The closest candidate PULs were found in *B. intestinalis*, *B. gallinarum*, *B. oleiciplenus* and *Proteiniphilum acetatigenes*. Interestingly, some of these PULs are adjacent to gene families that are found in PUL-XylL that encode GH5\_21, GH10, GH31, GH43, GH95 or GH115 enzymes (the degree of synteny is uncertain due to limited scaffold assemblies), suggesting a possible common evolutionary origin for the *B. ovatus* XylL- and XylS- PULs. Overall, of the 23 *B. ovatus* and *B. xylanisolvens* strains, only five have both XylL- and XylS-PULs, illustrating the dynamic acquisition (and possibly loss and/or reshuffling) of PULs by Bacteroidetes.

**Loss of xylan degradation products from the cell surface.**
Similar to other characterized *Bacteroides* polysaccharide degrading systems, the xylan degrading apparatus of *B. ovatus* appears to be optimized to maximize intracellular breakdown, which may indicate that the bacterium adopts a 'selfish' strategy when deconstructing the hemicellulose¹⁷⁻¹⁹. In contrast with this selfish hypothesis is the observation that on simple xylans, such as BGX and WAX, *B. ovatus* was able to support the growth of *Bifidobacterium adolescentis* strain ATCC 15703, which utilizes simple linear and arabino-xylooligosaccharides but not xylans⁴⁷ (Fig. 9a,b and Supplementary Fig. 7). The *Bacteroides*, however, was unable to promote the growth of the *Bifidobacterium* on

highly complex xylans such as CX (Fig. 9c). It could be argued that this differential capacity of *B. ovatus* to support the growth of *B. adolescentis* reflects the extent to which oligosaccharides from complex (CX) and simple (BGX and WAX) xylans are released into the culture medium and are thus available to the *Bifidobacterium*. It should be emphasized, however, that as we show here, CX breakdown requires a suite of enzymes not associated with canonical xylan degrading systems. Indeed only a limited number of closely related species of *Bacteroides* appear to contain these additional enzymes (see above and Fig. 8a). Thus, the inability of *B. ovatus* to support growth of *B. adolescentis* on CX may simply reflect the fact that the *Bifidobacterium* lacks the necessary apparatus to utilize this complex xylan. To explore this possibility a BACOVA_03433 GH98 mutant was constructed in which the key catalytic residues of the GH98 enzyme (Glu361 and Asp467) had been replaced with alanine. The mutant (ΔGH98) lacked an active surface GH98 xylanase and was unable to cleave the backbone of CX, preventing growth on the complex GAX (Fig. 9e). The *B. ovatus* variant, however, retained the other components of the xylan degrading apparatus, and was thus able to grow well on CX that had been pre-digested with recombinant GH98 (Fig. 9e). Co-culture of ΔGH98 with wild-type *B. ovatus* on intact CX supported the growth of both strains, demonstrating that PBPs are indeed released into the culture medium by wild-type *B. ovatus* during growth on CX (Fig. 9f). Thus, the inability of *B. adolescentis* to grow on CX in the presence of *B. ovatus* does not reflect the release of an inadequate supply of xylan-derived PBPs, but indicates that the *Bifidobacterium* lacks the apparatus required to utilize these highly complex oligosaccharides. This is consistent with the finding that *B. adolescentis* is unable to grow on CX that had been pre-digested with either the GH98 alone or all the surface located enzymes encoded by PUL-XylL (Supplementary Fig. 7).

<figure><img src="image_18.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_20.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_23.png" /></figure>

**Figure 9 | The ability of other members of the microbiota to use PBPs released by *B. ovatus* during growth on xylan is determined by the complexity of the polysaccharide.** *Bacteroides ovatus* (*Bo*) and *Bifidobacterium adolescentis* were co-cultured on BGX (a) WAX (b) and CX (c) and the number of CFUs of each determined at different points on the growth curve. (d) Shows the CFU ml⁻¹ of *B. ovatus* alone at different phases of growth. Note, when the *Bifidobacterium* alone was grown on digested xylans (see Supplementary Fig. 8), the CFU ml⁻¹ at late exponential phase (OD₆₀₀ ~ 1.2) was ~8.0 × 10⁸, indicating that ~25% of the total xylan is used by *B. adolescentis*. (e) Growth of *B. ovatus* wild type and the ΔGH98 mutant (BACOVA_03433) on CX (0.5% w/v). Pre-digested indicates the CX was digested to completion with the GH98 xylanase prior to addition to the media. (f) Tagged strains of wild-type *B. ovatus* and the ΔGH98 mutant were co-cultured on CX (see panel e for growth curve). Samples were taken at different time points and qPCR with primers unique to each strain was used to determine the ratio of each in the culture. Each data point is the average and s.d. of triplicate growths.
Currently the degree to which glycan degradation is mediated by synergistic interactions between different members of the gut microbiota is unclear. It has been proposed that there are keystone organisms belonging to the Firmicute phylum that make recalcitrant glycans, such as resistant starch, available to other members of the ecosystem⁴⁸. Similarly, a recent study showed that the products released during degradation of highly accessible glycans (branched starch and fructans) by *Bacteroides* spp. were shared with recipients that alone are unable to utilize these polysaccharides⁴⁹. In contrast, *B. thetaiotaomicron* utilizes complex yeast mannan in a selfish manner as breakdown occurs almost entirely intracellularly, with no evidence of sharing of PBPs with other *Bacteroides* spp.²⁷. This study extends our understanding of resource allocation in the microbiota by showing that selfish glycan utilization is not dependent solely on an ability to retain PBPs within the cell, but can be determined by the complexity of the enzyme systems required to access the extracellular oligosaccharides generated.

## Discussion

This report provides a model for the microbial utilization of xylan, a key polysaccharide present in the human diet. The data demonstrate that current understanding of xylan degradation, based primarily on the breakdown of relatively simple xylans by soil microbial ecosystems, is inadequate when considering the complex GAXs found in the human diet²². It is evident that *B. ovatus* has acquired an extensive range of enzymes, several from GH families that had not previously been implicated in xylan breakdown, to assemble an efficient catalytic apparatus capable of deconstructing highly complex forms of the polysaccharide. It is significant that the genetic basis for xylan degradation is organized into two distinct loci regulated by different signals. This provides *B. ovatus* with a flexible strategy to degrade the hemicellulose. When presented with a simple linear xylan or GX a minimum enzyme cocktail is assembled (that is, PUL-XylS), while PUL-XylL is activated only in response to complex GAXs (Fig. 2c). It is likely, therefore, that the distribution of these two PULs within human microbiomes provides insight into the nature of the diet consumed by the respective host.

Using xylans that contain variable decorations but a conserved core backbone, we have shown that it is the complexity of the target polysaccharide that determines whether other members of the HGM can utilize the PBPs released by *B. ovatus* during growth on different xylans. Thus, the relatively simple oligosaccharides that are released during growth on WAX and BGX can be used by *B. adolescentis*, while the more complex molecules released from CX cannot be accessed by the *Bifidobacterium*. Indeed, our analysis of Bacteroidetes genomes suggests that very few organisms, outwith *B. ovatus* and *B. xylanisolvens*, have a sufficiently complex xylan degrading system to utilize the PBPs released from CX by these two *Bacteroides* species. It is possible that these complex glycans are only used in a nutritional crisis when competition for resources is intense within the microbiota. Recently, it has been shown that *Bacteroides thetaiotaomicron* displays a degradative hierarchy in which some simple glycans are prioritized above more complex polysaccharides⁵⁰. Understanding the mechanisms by which different glycans are degraded by the microbiota, specifically whether they are public-good products or are used exclusively by the degrading organism, is critical to developing nutraceutical strategies to manipulate the structure of the microbiota to maximize human health. For example, complex GAXs will favour the growth of *B. ovatus*, which produces propionate, a molecule that reduces lipogenesis in the host⁵¹. In contrast, the degradation of simple xylans by *B. ovatus* will also promote the growth of

*Bifidobacteria*, which can lead to the production of butyrate, a molecule that is known to maintain the health of the intestinal epithelium⁵². Thus, it may be possible to develop bespoke xylan-based diets that are designed to have different beneficial effects on human health.

## Methods

**Isolation of corn xylan.** Oven-dried corn fibre sample was kindly provided by ADM Research (North America), which was further oven-dried and ground to a 20-mesh particle size using a Wiley mill⁵³. The ground corn fibre sample was extracted with hexane to remove oil and treated with thermostable α-amylase (a gift from Novo Nordisk Bioindustrials, Danbury, CT) at 90-95 °C to hydrolyse the starch present in the fibre⁵⁴,⁵⁵. CX was isolated from de-oiled and de-starched corn fibre following the alkaline hydrogen peroxide technology of Yadav *et al.*⁵⁶ with the following modifications. De-oiled and de-starched corn fibre (50 g) was mechanically stirred into water (1.0 l) and NaOH (8 ml from 50% solution), Ca(OH)₂ (3.8 g) and 42 ml of 30% H₂O₂ were carefully added in an open beaker in a fume hood. The mixture was boiled with mechanical stirring for 1 h. During the reaction, its pH was kept at 11.5 by addition of 50% NaOH as needed. After cooling the hot reaction mixture by stirring at room temp for an additional half an hour, it was centrifuged at 6,000 g for 20 min and the supernatant was separated from the residue by decantation. The pH of the alkaline H₂O₂ extract was then adjusted to 4.0-4.5 by adding conc. HCl to precipitate Hemicellulose A (acid-insoluble arabinoxylan, 'Hemi. A'), which was collected by centrifugation at 10,000 g for 30 min and discarded. Two volumes of ethanol (2.0 l) were gradually added to the supernatant (1.0 l), obtained after Hemi. A removal, with stirring to precipitate the major arabinoxylan fraction (CX). The CX was allowed to settle out as a white flocculent precipitate at the bottom of the beaker for 10-15 min. The clear alcohol/water mixture above the precipitate was removed by decantation. The white flocculent precipitate was collected by filtration and dried in a vacuum oven at 50 °C overnight.

The other xylans used in this study were low viscosity wheat arabinoxylan (WAX, Megazyme International) and birchwood glucuronoxylan (BGX; Sigma Aldrich).

**Cloning, expression and purification of recombinant proteins.** The open-reading frames encoding mature forms of the proteins used in the study were amplified from *Bacteroides ovatus* ATCC 8483 genomic DNA by PCR using appropriate primers (Supplementary Table 9), which introduce BamHI/HindIII and NdeI/XhoI sites to the flanks of the target genes. The amplified DNA was cloned into *Escherichia coli* expression vectors (see Supplementary Table 9) such that the encoded recombinant proteins contain either an N- or C-terminal His₆-tag. Site-directed mutants were made using the Quikchange kit (Agilent), using primers shown in Supplementary Table 10. The recombinant proteins were produced in *E. coli* BL21 (DE3) or Tuner (Novagen) cells cultured in LB broth containing ampicillin (100 µg ml⁻¹) or kanamycin (50 µg ml⁻¹) at 37 °C. Cells were grown to mid-exponential phase (OD₆₀₀ ~ 0.6), at which point isopropyl β-D-thiogalactopyranoside was added to a final concentration of 1 mM (BL21) or 0.2 mM (Tuner), and the cultures were incubated for a further 16 h at 16 °C. The cells were harvested by centrifugation, sonicated and His₆-tagged recombinant protein were purified from cell-free extracts by immobilized metal ion affinity chromatography using Talon resin (Clontech) and standard methodology as described previously³⁷. For crystallographic studies BACOVA_03438 was further purified by size exclusion chromatography using a Superdex 75 column (GE Healthcare Life Sciences). All proteins were purified to electrophoretic homogeneity as judged by SDS-PAGE and their concentrations determined from their calculated molar extinction coefficient at 280 nm using a NanoDrop 2000c (Thermo Scientific).

**Growth of anaerobic bacteria.** *Bacteroides ovatus* ATCC 8483 was routinely grown in tryptone-yeast extract-glucose (TYG) medium or on brain-heart infusion (BHI; Beckton Dickinson) agar plus porcine hematin (1.2 µg ml⁻¹; Sigma-Aldrich)⁴,⁵⁷. *Bifidobacterium adolescentis* ATCC 15703 was grown on clostridial nutrient medium (CNM; Sigma-Aldrich) plus hematin. Antibiotics were added as appropriate: erythromycin (25 µg ml⁻¹) and gentamicin (200 µg ml⁻¹). *Bacteroides* minimal medium (MM) was formulated as described previously³³ as was *Bifidobacterium* MM⁵⁸. Cells were grown in an anaerobic cabinet (Whitley A35 Workstation; Don Whitley, UK). Cultures (5 ml) were routinely grown in glass test tubes and the OD₆₀₀ monitored using a Biochrom WPA cell density meter. Growth of *B. ovatus* wild type and PUL-knockout mutants was measured on various xylan preparations by mixing the autoclave-sterilised polysaccharides (0.5% final concentration) with minimal medium and monitoring growth continuously as previously described⁴.

**Quantitative RT-PCR (qPCR).** Comparison of the levels of *susC* transcript expression from each of the xylan PULs was performed by qPCR. Previous studies have shown the *susC* genes are a good proxy for expression of their cognate PUL⁴. *B. ovatus* was cultured in 5 ml of minimal media containing 0.5% (w/v) carbon
source, as described above. Triplicate bacterial cultures were harvested at mid-log phase (OD~600~ ~0.8) and placed in RNAprotect (Qiagen), then stored at -80 °C overnight, before purification with RNEasy kit (Qiagen). RNA purity was assessed spectrophotometrically, and 1 µg of RNA was used immediately for reverse transcription (QuantiTect Reverse Transcription kit, Qiagen). Quantitative RT-PCR was performed in a 96-well plate on a LightCycler 480 System (Roche) with FastStart Essential DNA Green Master (Roche) using the primers shown in Supplementary Table 11. Reactions were carried out in 10 µl, consisting of 5 µl SYBR Green mix, 20 ng of cDNA, and 1 µM (susC genes) or 0.125 µM (16 S rRNA) primer mix. Reaction conditions were 95 °C 600 s, followed by 45 cycles of 95 °C for 10 s, 55 °C for 10 s, 72 °C for 10 s. Cq values were calculated using LightCycler 480 SW 1.5. Data were normalized to 16 S rRNA transcript levels, and change in expression level calculated as fold-change compared with minimal media, glucose cultures.

**Cell localization studies.** Potential surface-located enzymes (containing either a predicted lipoprotein signal sequence or uncleaved TM anchor) encoded by the *B. ovatus* xylan PULs were initially identified using the LipoP 1.0 (http://www.cbs.dtu.dk/services/LipoP/) and SigP 4.1 (http://www.cbs.dtu.dk/services/SignalP/) servers. Note, the database N-termin term methionine for BACOVA_03425 GH34 is mis-annotated; the correct N-terminus actually starts 24aa downstream at the sequence MKN.

*Proteinase K treatment.* Cultures of *B. ovatus* (50 ml) were grown in minimal media on WAX (0.5% w/v) as a sole carbon source, to mid-exponential growth phase (OD~600~ ~0.8). Cells were harvested by centrifugation and washed twice in 10 ml phosphate buffered saline pH 7.2 (PBS), before being resuspended in 2.5 ml of the buffer. The cells were split into four 0.5 ml aliquots. To three of the aliquots 2 mg ml⁻¹ Proteinase K was added and incubated at 37 °C for up to 16 h, the fourth sample was left as an untreated control for 16 h. Following incubation with the proteinase the samples were centrifuged at 5,000 g for 10 min and the supernatant discarded. The cell pellets were resuspended in 1 ml PBS and remaining ProteinaseK precipitated by the addition of 200 µl trichloroacetic acid and incubation on ice for 30 min. The cells were pelleted by centrifugation and washed four times in 1 ml ice cold acetone. The cell pellets were resuspended in 250 µl Laemml buffer and subjected to SDS-PAGE. Proteins were transferred to Whatman Protram BA 85 nitrocellulose membrane using a wet transfer system (BioRad Mini Trans-Blot). Proteins of interest were detected using anti-sera raised in rats (Eurogentec) against the corresponding recombinant form of the protein, diluted 1/1,000 or 1/2,000 in 1% milk powder in TBS buffer. The secondary antibody used was a chicken anti-rat conjugated to horseradish peroxidase (Santa Cruz) diluted 1/5,000 in 1% milk powder in TBS buffer. Antibodies were detected by chemiluminescence using Biorad Clarity Western ECL Substrate.

*Immunofluorescence microscopy.* *B. ovatus* suspensions of mid-exponential phase cells (0.5 ml in PBS, produced as for Proteinase K treatment above) were fixed with an equal volume of 2 × formalin (9% formaldehyde in PBS), and rocked for 90 min at 25 °C in Eppendorf tubes. The cells were then pelleted by centrifugation for 3 min at 7,000 g and washed twice with 1 ml of PBS. The bacterial cell pellet was resuspended in 1 ml of blocking solution (2% goat serum, 0.02% NaN₃ in PBS) and incubated at 4 °C for 16 h. After incubation cells were centrifuged again at 7,000 g and the supernatant discarded. For labelling, the bacteria were incubated with 0.5 ml of primary rat IgG (1/500 dilution of IgG in 1% milk powder in PBS, blocking solution) for 2 h at 25 °C. The cells were then pelleted, washed in 1 ml of PBS and resuspended in 0.4 ml goat anti-rat IgG Alexa-Fluor 488 (Sigma-Aldrich), diluted 1/500 in blocking solution and incubated 1 h at 25 °C in the dark. The cells were again pelleted, washed with PBS, resuspended in 1 ml of PBS containing ProLong Gold antifade reagent (Life Technologies). Aliquots (50 µl) of labelled bacterial cells were mounted onto glass slides and secured with coverslips. Fluorescence was visualized using a Leica SP2 UV microscope (Leica Microsystems, Heidelberg, GmbH) with ×63 NA1.32 lens. Alexa-Fluor 488-labelled bacteria were visualized under a UV view and compared with bright-field phase-contrast of the same image.

**Whole-cell assays.** Whole-cell assays under aerobic conditions were used to determine surface enzyme activity uncoupled from glycan import. *B. ovatus* was grown in 5 ml minimal media containing WAX (0.5% w/v; activates both xylan PULs⁴), or glucose as the sole carbon source to mid exponential phase (OD~600~ ~0.8). Cells were harvested by centrifugation at 4,000 g for 10 min at room temperature and washed in 5 ml PBS, pH 7.2, before being resuspended in 500 µl PBS. Cells (50 µl) were assayed against BGX, WAX and CX at 37 °C for 16 h. Assays were analysed by thin layer chromatography, 5 µl of each sample was spotted onto silica plates and resolved in butanol/acetic acid/water buffer (2:1:1). The plates were dried and carbohydrates visualized by orcinol/sulphuric acid heated to 70 °C for 10 min.

**Glycan-binding studies.** The binding of proteins to their glycan ligands was quantified by isothermal titration calorimetry (ITC), as described previously⁵⁹. Titrations were carried out in 50 mM Na-HEPES buffer, pH 7.5 at 25 °C. The reaction cell contained protein at 50–100 µM, while the syringe contained either the oligosaccharide at 1–10 mM or the polysaccharide at 3–10 mg ml⁻¹. Integrated heats were fit to a single-site model using Microcal Origin v7.0 to derive *n*, *K*<sub>a</sub>, and *ΔH* values.

*Glycoside hydrolase assays.* *Spectrophotometric assays.* The activity of the GH31 α-xylosidase and GH3 and GH34 β-xylosidases against xylans and xylooligosaccharides were quantified by the continuous monitoring of xylose release using a D-xylose detection kit (Megazyme International). The reaction was carried out at 37 °C in 50 mM sodium phosphate buffer, pH 7.0 containing 2 mM MgCl₂, 1 mM ATP and 1 mM NAD⁺, excess concentrations of linked enzymes (hexokinase, xylose mutarotase and β-xylose dehydrogenase) and 1 mg ml⁻¹ BSA. The amount of NADH formed in this reaction is stoichiometric with the amount of D-xylose released and NADH levels were monitored at 340 nm using an extinction coefficient of 6,220 M⁻¹ cm⁻¹. The xyloglucan heptasaccharide and linear xylooligosaccharide substrates used in this assay were purchased from Megazyme International. The glucurono-xylooligosaccharides, which were also used as substrates, were generated as described previously³⁷. In brief, BGX was digested with CjXyn10A or NpXyn11A to generate UXX or XUXX, respectively, as the main charged oligosaccharides. These were purified from neutral xylooligosaccharides using Dowex chromatography and freeze dried. For kinetic assays of the GH34 α-L-arabinofuranosidases and GH67 α-D-glucuronidase, the L-arabinose/n-D-galactose assay and D-glucuronic acid assay kits (Megazyme International) were used, respectively³⁷. The rye arabinoxylan used to assay the arabinofuranosidases was from Megazyme International.

The activity of the GH10 and GH30 endo-xylanases against BGX was determined in 50 mM sodium phosphate buffer, pH 7.0 at 37 °C containing an appropriate concentration of the polysaccharide and 1 mg ml⁻¹ BSA. Reactions were incubated at 37 °C and at regular time intervals 500 µl aliquots were removed and the amount of reducing sugar was quantified using the dinitrosalicylic acid reagent (Sigma Aldrich) at 575 nm and a standard curve of xylose in the reaction conditions used. HPAEC-PAD was used to determine the profile of the xylooligosaccharides generated by the endo-xylanases from these reactions and to measure the relative activity of the GH10 enzymes against specific xylooligosaccharides.

To measure the activity of the GH95 enzymes, the L-fucose and L-galactose products released from 2'-fucosyllactose (Carbosynth) and CX, respectively, were monitored by HPAEC-PAD. The reactions were carried out in 20 mM sodium phosphate buffer, pH 7.0 at 37 °C. In brief substrates were incubated with an appropriate concentration of each enzyme (10 nM to 5 µM). Aliquots were removed at regular intervals for up to 1 h and, after boiling for 10 min to inactivate the enzyme, the amount of the product was quantified by comparison to known standards. The activity of wild type and mutant BoGH95 (10 µM final enzyme in 20 mM sodium phosphate buffer, pH 7.0 at 37 °C) against 2-chloro-4-nitrophenyl α-L-fucose (Carbosynth) was determined by continuously monitoring the release of 2-chloro-4-nitrophenol product at 400 nm.

The molar concentration of polysaccharide available to exo acting enzymes was determined by incubating the glycans at 1 mg ml⁻¹ with a large excess of enzyme for 1 h to ensure the reaction goes to completion. The molar concentration of product at this point was then determined using the appropriate monosaccharide detection kit (detailed above). This equates to the molar concentration of total available substrate, thus enabling mg ml⁻¹ of substrate to be converted to molarity.

Where *k*<sub>cat</sub>/<i>K</i><sub>m</sub> alone is reported this was determined by monitoring the initial rate at a range of substrate concentrations significantly below the *K*<sub>m</sub> such that a linear plot of *V*₀ vs [S] was obtained. Thus, the initial rate of substrate hydrolysis gives a direct readout of *k*<sub>cat</sub>/<i>K</i><sub>m</sub> using the equation *k*<sub>cat</sub>/<i>K</i><sub>m</sub> = *V*₀/[S]<sup>*n</sup>[E].

*High-performance anion-exchange chromatography.* Products generated by glycoside hydrolases used in the study were analysed using an analytical CARBOPAC™ PA-100 anion-exchange column (Dionex) equipped with a guard column. The system had a loop size of 100 µl and was run at a flow rate of 1 ml min⁻¹ at a pressure of ~2,300 lbf in⁻². Carbohydrates were detected by pulsed amperometric detection (PAD) with electrode settings of E₁ = +0.05, E₂ = +0.6 and E₃ = -0.6. The elution conditions were 0–10 min, 100 mM NaOH; 10–25 min, 100 mM NaOH with a 0–75 mM sodium acetate gradient; and 25–35 min, 100 mM NaOH containing 500 mM sodium acetate. For analysis of the very large oligosaccharides produced by the GH98 xylanase a 0-100% sodium acetate gradient from 0 to 90 mins was used. The standards used to identify the chromatographic peaks were arabinose, xylose, xylobiose, xylotriose, xylotetraose, xylopentaose and xylohexaose.

*Analysis of the GH98 (BACOVA_03433) products reduced with sodium borohydride.* To analyse the specificity of BACOVA_03433, the recombinant protein (5 µM final) was incubated with 20 mg of CX resuspended in 0.5 ml 50 mM sodium phosphate buffer, pH 7.0 at 37 °C for 24 h. The same amount of CX was incubated separately without the enzyme as a control. Both samples were then incubated overnight with NaBH₄ (10 mg per sample) in the presence of NH₄OH, pH 10.8 at 25 °C. After incubation, samples were freeze-dried, resuspended in 500 µl of water and loaded separately on Bio-gel P2 size-exclusion columns (Bio-Rad laboratories Inc.). Water was applied to the column (0.5 ml min⁻¹) using a Pharmacia CKB peristaltic pump to elute carbohydrate products. A 4 µl aliquot of each of the 120 fractions (2 ml) was spotted on silica gel TLC plates, which were developed in butanol:acetic acid:water 2:1:1 and carbohydrate products detected by spraying with 0.5% (w/v) orcinol in 10% sulphuric acid and heating to 70 °C for 10 min. Those containing oligosaccharide were pooled together and dried using a speed-vac. Each sample (enzyme and control) was then resuspended in 20 µl of water and hydrolysed with 0.3 M HCl for 3 h at 100 °C. The products of acid
hydrolysis were neutralized with 1 M NaOH and analysed by HPAEC using standard methodology (see above). The identity of the reduced sugar product was determined by comparison to known standards (galactitol and xylitol).

**X-ray crystallography.** BACOVA_03438 GH95, at 10 mg ml⁻¹, was crystallized from 160 mM calcium acetate, 80 mM sodium cacodylate pH 6.5, 14.4% w/v polyethylene glycol 8,000 and 20% v/v glycerol using the vapour diffusion sitting drop method. Crystals of BACOVA_03438 in complex with L-galactose were obtained by preparing protein:monosaccharide mixtures at 10 mg ml⁻¹; 300 mM and mixing 1:1 or 2:1 with the reservoir solution grown at 20 °C. Crystals were mounted in mother liquor and cryocooled in liquid N₂. Data for BACOVA_03438 were collected at Diamond Light Source, UK, beamline I02, to 2.81 Å resolution. Data were processed and integrated with XDS and scaled using Aimless⁶⁰,⁶¹. All other computing used the CCP4 suite of programs⁶². The space group was determined to be P2₁ 2₁ 2₁ with 2 molecules in the asymmetric unit (giving a Matthews coefficient of 2.61 Å³ Da⁻¹ and a solvent content of 53%). The structure of BACOVA_03438 was solved by molecular replacement in MOLREP using pdb file with accession code 2RDY as the search model. The model underwent cycles of model building in COOT⁶³ and refinement in REFMAC⁶⁴. Five percent of the observations were used to monitor the progress of refinement. The model was submitted to the PDB_REDO server to improve final model quality⁶⁵. The model was validated using Molprobity⁶⁶. The crystal and data collection statistics are reported in Supplementary Table 12.

**Genetic manipulation of *Bacteroides ovatus***. Knockout strains lacking either PUL-XylL, PUL-XylS or both PULs were created using counterselectable allelic exchange as described previously³⁰. The inactive mutant GH98 (BACOVA_03433) in pET21 lacking both catalytic residues was made using the Quikchange site directed mutagenesis kit (Agilent) following the manufacturer's protocol (Supplementary Table 10). The mutated gene was then cloned into pExchange and integrated into the genome of *B. ovatus* using the same counterselectable allelic exchange methodology as for the PUL knockouts³⁰. Signature-tagged strains of wild type and the ΔGH98 mutant, for determining the ratios of each during co-culture by qPCR, were made by integration of unique sequences into one of two *att* sites in the *Bacteroides ovatus* genome through homologous recombination of pNBU2 conjugated into the *Bacteroides* from S17-1 λpir *E. coli* as described previously³. Following positive selection on BHI-agar contain gentamycin (200 µg ml⁻¹) and tetracycline (1 µg ml⁻¹), genomic DNA was extracted (QIAamp DNA Mini Kit; Qiagen) and PCR with primers flanking each possible *att* insertion site was used to select for *Bacteroides* strains that had successfully integrated the tag sequence. Primers used for all genetic manipulations are shown in Supplementary Table 13.

**Sharing assays.** *Bifidobacterium adolescentis* and *Bacteroides ovatus* cells were grown anaerobically to mid/late-log phase in CNM (Sigma-Aldrich) and TYG rich media, respectively. Cells were washed in PBS and used, in approximately equal proportions, to inoculate Minimal Medium (based on medium devised by Van der Meulen et al.⁵⁸) containing 0.5% (w/v) BGX, WAX or CX as the sole carbon source. During co-culture of *Bacteroides* and *Bifidobacterium* aliquots were taken over 18 h. Serial dilutions of the samples were plated onto CNM agar and incubated anaerobically for 48 h at 37 °C to obtain CFU of each sample. Single colonies were picked and replica plated onto CNM and CNM with 200 µg ml⁻¹ gentamycin to differentiate between *Bacteroides* and *Bifidobacterium* as *Bacteroides* spp. are resistant to this antibiotic. The percentage of each bacterium present was used, in conjunction with total CFUs was used to calculate the CFU of each.

To compare ratios of the ΔGH98 when grown in co-culture with wild-type *B. ovatus* on CX, 5 ml minimal medium containing 0.5% CX was inoculated with tagged ΔGH98 and wild type overnights grown on TYG in a 50:50 ratio (estimated by OD₅₀₀nm). Samples (0.5 ml) were taken during growth, 50 µl of which was serial diluted and plated onto to BHI agar, while genomic DNA was isolated from the remaining cells. gDNA samples were subject to quantitative PCR (Roche LightCycler 480 System; see above) using tag-specific primers (Supplementary Table 13) to give relative quantities of the mutant and wild-type cells in the culture. The ratio obtained by qPCR with total CFU counts was used to deduce CFU of mutant and wild type at the different phases of growth.

To analyse the growth of *B. adolescentis* on digested xylans, *Bifidobacterium* minimal medium containing with 0.5% WAX, BGX or CX that had been digested to completion with either recombinant BACOVA_04390 GH10 for WAX and BGX, or for CX with BACOVA_03433 GH98 alone, or all surface-located enzymes from PUL-XylL (BACOVA_03419 GH3, 03421 GH43, 03432 GH30 and 03433 GH98), was inoculated with a *B. adolescentis* overnight that had been grown on rich medium (CNM). Aliquots were taken at different phases of growth, serial diluted and plated onto CNM agar to determine CFU at each growth phase.

**Comparative genomics analysis.** PULs similar to PUL-XylS and PUL-XylL were searched for in 283 Bacteroidetes genomes (complete list provided in Supplementary Table 8). The identification of similar PULs was based on PUL alignments. Gene composition and order of Bacteroidetes PULs were first computed using the PUL predictor described in PULDB⁶⁷. Then, in a manner similar to amino-acid sequence alignments, the predicted PULs were aligned to PUL-XylS

and PUL-XylL according to their modularity as proposed in RADS/RAMPAGE method⁶⁸. Modules taken into account include CAZy families, sensor-regulators and susCD-like genes. Finally, PUL boundaries and limit cases were refined by BLASTP-based analysis.

References

1. Arumugam, M. et al. Enterotypes of the human gut microbiome. *Nature* **473**, 174–180 (2011).

2. Backhed, F., Ley, R. E., Sonnenburg, J. L., Peterson, D. A. & Gordon, J. I. Host-bacterial mutualism in the human intestine. *Science* **307**, 1915–1920 (2005).

3. Martens, E. C., Chiang, H. C. & Gordon, J. I. Mucosal glycan foraging enhances fitness and transmission of a saccharolytic human gut bacterial symbiont. *Cell Host Microbe* **4**, 447–457 (2008).

4. Martens, E. C. et al. Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts. *PLoS Biol.* **9**, e1001221 (2011).

5. Flint, H. J., Scott, K. P., Duncan, S. H., Louis, P. & Forano, E. Microbial degradation of complex carbohydrates in the gut. *Gut Microbes* **3**, 289–306 (2012).

6. De Vuyst, L., Moens, F., Selak, M., Riviere, A. & Leroy, F. Summer Meeting 2013: growth and physiology of bifidobacteria. *J. Appl. Microbiol.* **116**, 477–491 (2014).

7. Roberfroid, M. et al. Prebiotic effects: metabolic and health benefits. *Br. J. Nutr.* **104**(J): S1–63 (2010).

8. Chapman, C. M., Gibson, G. R. & Rowland, I. Health benefits of probiotics: are mixtures more effective than single strains? *Eur. J. Nutr.* **50**, 1–17 (2011).

9. Gibson, M. K., Pesesky, M. W. & Dantas, G. The Yin and Yang of Bacterial Resilience in the Human Gut Microbiota. *J. Mol. Biol.* **426**, 3866–3876 (2014).

10. Martens, E. C., Kelly, A. G., Tauzin, A. S. & Brumer, H. The devil lies in the details: how variations in polysaccharide fine-structure impact the physiology and evolution of gut microbes. *J. Mol. Biol.* **426**, 3851–3865 (2014).

11. El Kaoutari, A., Armougom, F., Gordon, J. I., Raoult, D. & Henrissat, B. The abundance and variety of carbohydrate-active enzymes in the human gut microbiota. *Nat. Rev. Microbiol.* **11**, 497–504 (2013).

12. Flint, H. J., Bayer, E. A., Rincon, M. T., Lamed, R. & White, B. A. Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. *Nat. Rev. Microbiol.* **6**, 121–131 (2008).

13. Agger, J., Vikso-Nielsen, A. & Meyer, A. S. Enzymatic xylose release from pretreated corn bran arabinoxylan: differential effects of deacetylation and deferuloylation on insoluble and soluble substrate fractions. *J. Agric. Food Chem.* **58**, 6141–6148 (2010).

14. Allerdings, E., Ralph, J., Steinhart, H. & Bunzel, M. Isolation and structural identification of complex feruloylated heteroxylan side-chains from maize bran. *Phytochemistry* **67**, 1276–1286 (2006).

15. Brett, C. T. & Waldren, K. *Physiology and Biochemistry of Plant Cell Walls. Topics in Plant Functional Biology* vol. 1Chapman and Hall, 1996).

16. Lombard, V., Golaconda Ramulu, H., Drula, E., Coutinho, P. M. & Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Res.* **42**, D490–D495 (2014).

17. Larsbrink, J. et al. A discrete genetic locus confers xyloglucan metabolism in select human gut Bacteroidetes. *Nature* **506**, 498–502 (2014).

18. Martens, E. C., Koropatkin, N. M., Smith, T. J. & Gordon, J. I. Complex glycan catabolism by the human gut microbiota: The bacteroidetes Sus-like paradigm. *J. Biol. Chem.* **284**, 24673–24677 (2009).

19. Sonnenburg, E. D. et al. Specificity of polysaccharide use in intestinal bacteroides species determines diet-induced microbiota alterations. *Cell* **141**, 1241–1252 (2010).

20. Dodd, D., Moon, Y. H., Swaminathan, K., Mackie, R. I. & Cann, I. K. O. Transcriptomic analyses of xylan degradation by prevotella bryantii and insights into energy acquisition by xylanolytic bacteroidetes. *J. Biol. Chem.* **285**, 30261–30273 (2010).

21. Zhang, M. et al. Xylan utilization in human gut commensal bacteria is orchestrated by unique modular organization of polysaccharide-degrading enzymes. *Proc. Natl Acad. Sci. USA* **111**, E3708–E3717 (2014).

22. Gilbert, H. J. The biochemistry and structural biology of plant cell wall deconstruction. *Plant Physiol.* **153**, 444–455 (2010).

23. Charnock, S. J. et al. Key residues in subsite F play a critical role in the activity of Pseudomonas fluorescens subspecies cellulosa xylanase A against xylooligosaccharides but not against highly polymeric substrates such as xylan. *J. Biol. Chem.* **272**, 2942–2951 (1997).

24. Pell, G. et al. The mechanisms by which family 10 glycoside hydrolases bind decorated substrates. *J. Biol. Chem.* **279**, 9597–9605 (2004).

25. Charnock, S. J. et al. The X6 "thermostabilizing" domains of xylanases are carbohydrate-binding modules: structure and biochemistry of the *Clostridium thermocellum* X6b domain. *Biochemistry* **39**, 5013–5021 (2000).

26. Tailford, L. E. et al. Understanding how diverse beta-mannanases recognize heterogeneous substrates. *Biochemistry* **48**, 7009–7018 (2009).

27. Cuskin, F. et al. Human gut Bacteroidetes can utilize yeast mannan through a selfish mechanism. *Nature* **517**, 165–169 (2015).
28. St John, F. J., Hurlbert, J. C., Rice, J. D., Preston, J. F. & Pozharski, E. Ligand bound structures of a glycosyl hydrolase family 30 glucuronoxylan xylanohydrolase. *J. Mol. Biol.* **407**, 92–109 (2011).

29. Higgins, M. A. et al. Differential recognition and hydrolysis of host carbohydrate antigens by Streptococcus pneumoniae family 98 glycoside hydrolases. *J. Biol. Chem.* **284**, 26161–26173 (2009).

30. Koropatkin, N. M., Martens, E. C., Gordon, J. I. & Smith, T. J. Starch catabolism by a prominent human gut symbiont is directed by the recognition of amylose helices. *Structure* **16**, 1105–1115 (2008).

31. Cameron, E. A. et al. Multidomain carbohydrate-binding proteins involved in Bacteroides thetaiatoamicron starch metabolism. *J. Biol. Chem.* **287**, 34614–34625 (2012).

32. Pell, G. et al. Structural and biochemical analysis of Cellvibrio japonicus xylanase 10C: how variation in substrate-binding cleft influences the catalytic profile of family GH-10 xylanases. *J. Biol. Chem.* **279**, 11777–11788 (2004).

33. Patil, D. N. et al. Structural investigation of a novel N-acetyl glucosamine binding chi-lectin which reveals evolutionary relationship with class III chitinases. *PLoS ONE* **8**, e63779 (2013).

34. Beylot, M. H., McKie, V. A., Voragen, A. G., Doeswijk-Voragen, C. H. & Gilbert, H. J. The Pseudomonas cellulosa glycoside hydrolase family 51 arabinofuranosidase exhibits wide substrate specificity. *Biochem. J.* **358**, 607–614 (2001).

35. McKee, L. S. et al. Introducing endo-xylanase activity into an exo-acting arabinofuranosidase that targets side chains. *Proc. Natl Acad. Sci. USA* **109**, 6537–6542 (2012).

36. Boraston, A. B., Bolam, D. N., Gilbert, H. J. & Davies, G. J. Carbohydrate-binding modules: fine-tuning polysaccharide recognition. *Biochem. J.* **382**, 769–781 (2004).

37. Rogowski, A. et al. Evidence that GH115 alpha-glucuronidase activity, which is required to degrade plant biomass, is dependent on conformational flexibility. *J. Biol. Chem.* **289**, 53–64 (2014).

38. Nurizzo, D., Nagy, T., Gilbert, H. J. & Davies, G. J. The structural basis for catalysis and specificity of the Pseudomonas cellulosa alpha-glucuronidase, GlcA67A. *Structure* **10**, 547–556 (2002).

39. Frandsen, T. P. & Svensson, B. Plant alpha-glucosidases of the glycoside hydrolase family 31. Molecular properties, substrate specificity, reaction mechanism, and comparison with family members of different origin. *Plant Mol. Biol.* **37**, 1–13 (1998).

40. Larsbrink, J. et al. Structural and enzymatic characterization of a glycoside hydrolase family 31 alpha-xylosidase from Cellvibrio japonicus involved in xyloglucan saccharification. *Biochem. J.* **436**, 567–580 (2011).

41. Katayama, T. et al. Molecular cloning and characterization of Bifidobacterium bifidum 1,2-alpha-L-fucosidase (AfcA), a novel inverting glycosidase (glycoside hydrolase family 95). *J. Bacteriol.* **186**, 4885–4893 (2004).

42. Nagae, M. et al. Structural basis of the catalytic reaction mechanism of novel 1,2-alpha-L-fucosidase from *Bifidobacterium bifidum*. *J. Biol. Chem.* **282**, 18497–18509 (2007).

43. Pluvinage, B., Hehemann, J. H. & Boraston, A. B. Substrate recognition and hydrolysis by a family 50 exo-beta-agarase, Aga50D, from the marine bacterium Saccharophagus degradans. *J. Biol. Chem.* **288**, 28078–28088 (2013).

44. Davies, G. & Henrissat, B. Structures and mechanisms of glycosyl hydrolases. *Structure* **3**, 853–859 (1995).

45. St John, F. J., Gonzalez, J. M. & Pozharski, E. Consolidation of glycosyl hydrolase family 30: a dual domain 4/7 hydrolase family consisting of two structurally distinct groups. *FEBS Lett.* **584**, 4435–4441 (2010).

46. Aspeborg, H., Coutinho, P. M., Wang, Y., Brumer, 3rd H. & Henrissat, B. Evolution, substrate specificity and subfamily classification of glycoside hydrolase family 5 (GH5). *BMC Evol. Biol.* **12**, 186 (2012).

47. Pastell, H., Westermann, P., Meyer, A. S., Tuomainen, P. & Tenkanen, M. *In vitro* fermentation of arabinoxylan-derived carbohydrates by bifidobacteria and mixed fecal microbiota. *J. Agric. Food Chem.* **57**, 8598–8606 (2009).

48. Ze, X., Le Mougen, F., Duncan, S. H., Louis, P. & Flint, H. J. Some are more equal than others: the role of "keystone" species in the degradation of recalcitrant substrates. *Gut Microbes* **4**, 236–240 (2013).

49. Rakoff-Nahoum, S., Coyne, M. J. & Comstock, L. E. An ecological network of polysaccharide utilization among human intestinal symbionts. *Curr. Biol.* **24**, 40–49 (2014).

50. Rogers, T. E. et al. Dynamic responses of Bacteroides thetaiatoamicron during growth on glycan mixtures. *Mol. Microbiol.* **88**, 876–890 (2013).

51. Hosseini, E., Grootaert, C. & Verstraete, W. Van de Wiele T. Propionate as a health-promoting microbial metabolite in the human gut. *Nutr. Rev.* **69**, 245–258 (2011).

52. Wachterhauser, A. & Stein, J. Rationale for the luminal provision of butyrate in intestinal diseases. *Eur. J. Nutr.* **39**, 164–171 (2000).

53. Yadav, M. P., Moreau, R. A., Hotchkiss, A. T. & Hicks, K. B. A new corn fiber gum polysaccharide isolation process that preserves functional components. *Carbohydr. Polym.* **87**, 1169–1175 (2012).

54. Yadav, M. P., Johnston, D. B., Hotchkiss, A. T. & Hicks, K. B. Corn fiber gum: A potential gum arabic replacer for beverage flavor emulsification. *Food Hydrocarbon* **21**, 1022–1030 (2007).

55. Moreau, R. A., Powell, M. J. & Hicks, K. B. Extraction and quantitative analysis of oil from commercial corn fiber. *J. Agric. Food Chem.* **44**, 2149–2154 (1996).

56. Yadav, M. P., Cooke, P., Johnston, D. B. & Hicks, K. B. Importance of protein-rich components in emulsifying properties of corn fiber gum. *Cereal Chem.* **87**, 89–94 (2010).

57. Eley, A., Greenwood, D. & O'Grady, F. Comparative growth of Bacteroides species in various anaerobic culture media. *J. Med. Microbiol.* **19**, 195–201 (1985).

58. Van der Meulen, R., Makras, L., Verbrugghe, K., Adriany, T. & De Vuyst, L. In vitro kinetic analysis of oligofructose consumption by Bacteroides and Bifidobacterium spp. indicates different degradation mechanisms. *Appl. Environ. Microbiol.* **72**, 1006–1012 (2006).

59. Szabo, L. et al. Structure of a family 15 carbohydrate-binding module in complex with xylopentaose. Evidence that xylan binds in an approximate 3-fold helical conformation. *J. Biol. Chem.* **276**, 49061–49065 (2001).

60. Kabsch, W. XDS. *Acta. Crystallogr. D. Biol. Crystallogr.* **66**, 125–132 (2010).

61. Evans, P. R. & Murshudov, G. N. How good are my data and what is the resolution? *Acta. Crystallogr. D. Biol. Crystallogr.* **69**, 1204–1214 (2013).

62. Winn, M. D. et al. Overview of the CCP4 suite and current developments. *Acta. Crystallogr. D. Biol. Crystallogr.* **67**, 235–242 (2011).

63. Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. *Acta. Crystallogr. D. Biol. Crystallogr.* **66**, 486–501 (2010).

64. Vagin, A. A. et al. REFMAC5 dictionary: organization of prior chemical knowledge and guidelines for its use. *Acta. Crystallogr. D. Biol. Crystallogr.* **60**, 2184–2195 (2004).

65. Joosten, R. P., Joosten, K., Cohen, S. X., Vriend, G. & Perrakis, A. Automatic rebuilding and optimization of crystallographic structures in the Protein Data Bank. *Bioinformatics* **27**, 3392–3398 (2011).

66. Chen, V. B. et al. MolProbity: all-atom structure validation for macromolecular crystallography. *Acta. Crystallogr. D. Biol. Crystallogr.* **66**, 12–21 (2010).

67. Terrapon, N., Lombard, V., Gilbert, H. J. & Henrissat, B. Automatic prediction of polysaccharide utilization loci in Bacteroidetes species. *Bioinformatics* **31**, 647–655 (2014).

68. Terrapon, N., Weiner, J., Grath, S., Moore, A. D. & Bornberg-Bauer, E. Rapid similarity search of proteins using alignments of domain arrangements. *Bioinformatics* **30**, 274–281 (2014).

69. Raman, R. et al. Advancing glycomics: implementation strategies at the consortium for functional glycomics. *Glycobiology* **16**, 82R–90R (2006).

## Acknowledgements

We thank Dr Bruce Hamaker (Purdue University, IN, USA) for the kind gift of rice and sorghum xylans. This work was supported in part by grants to D.N.B. (BBSRC BB/G016186/1) and H.J.G. (Wellcome Trust WT097907AIA).

## Author contributions

D.N.B. and H.J.G. conceived the study. A.R., J.A.B, N.T., E.C.L, A.B., A.M.D., H.Z., T.E.R. E.C.M, P.D., B.H., H.J.G. and D.N.B. designed the experiments. A.R., J.A.B, N.T., E.C.L, A.B., C.M., A.M.D., H.Z. and T.E.R. performed the experiments. A.R., J.A.B J.C.M., T.T., N.T., E.C.L, A.B., A.M.D., H.Z., T.E.R., E.C.M, P.D., B.H., H.J.G. and D.N.B. analysed data. P.T, A.R.H. and M.P.Y provided new reagents. B.H., N.T., H.J.G. and D.N.B. wrote the manuscript.

## Additional information

**Accession codes:** Coordinates and structure factors for BACOVA_03438 have been deposited in the RCSB Protein Data Bank under accession code 4UFC.

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission information** is available online at http://npg.nature.com/reprintsandpermissions/

**How to cite this article:** Rogowski, A. *et al.* Glycan complexity dictates microbial resource allocation in the large intestine. *Nat. Commun.* **6**:7481 (2015).

<figure><img src="image_24.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
# Corrigendum: Glycan complexity dictates microbial resource allocation in the large intestine

Artur Rogowski, Jonathon A. Briggs, Jennifer C. Mortimer, Theodora Tryfona, Nicolas Terrapon, Elisabeth C. Lowe, Arnaud Baslé, Carl Morland, Alison M. Day, Hongjun Zheng, Theresa E. Rogers, Paul Thompson, Alastair R. Hawkins, Madhav P. Yadav, Bernard Henrissat, Eric C. Martens, Paul Dupree, Harry J. Gilbert & David N. Bolam

*Nature Communications* 6:7481 doi: 10.1038/ncomms8481 (2015); Published 26 Jun 2015; Updated 5 Feb 2016

The financial support for the work described in this Article was not fully acknowledged. The Acknowledgements should have included the following:
This work was supported in part by a grant to H.J.G., B.H. and N.T. (European Union's Seventh Framework Program (FP/2007/2013)/European Research Council (ERC) Grant Agreement 322820).

<figure><img src="image_25.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/