PULDB: the expanded database of Polysaccharide Utilization Loci

Carbohydrate Hydrolytic Potential and Redundancy of an Anaerobic Digestion Microbiome Exposed to Acidosis, as Uncovered by Metagenomics

The enzymatic hydrolysis of lignocellulosic biomass is mainly driven by the action of carbohydrate-active enzymes. By characterizing the gene profiles at the different stages of the anaerobic digestion experiment, we showed that the microbiome retains its hydrolytic functional redundancy even during severe acidosis, despite significant changes in taxonomic composition. By analyzing reconstructed bacterial genomes, we demonstrate that Bacteroidetes hydrolytic gene diversity likely favors the abundance of this phylum in some anaerobic digestion systems. Further, we observe genetic redundancy within the Bacteroidetes group, which accounts for the preserved hydrolytic potential during acidosis. This work also uncovers new polysaccharide utilization loci involved in the deconstruction of various biomasses and proposes the model of acetylated glucomannan degradation by Bacteroidetes. Acetylated glucomannan-enriched biomass is a common substrate for many industries, including pulp and paper production. Using naturally evolved cocktails of enzymes for biomass pretreatment could be an interesting alternative to the commonly used chemical pretreatments.

Increased hydrolysis of easily digestible biomass may lead to acidosis of anaerobic reactors and decreased methane production. Previously, it was shown that the structure of microbial communities changed during acidosis; however, once the conditions are back to optimal, biogas (initially CO₂) production quickly restarts. This suggests the retention of the community functional redundancy during the process failure. In this study, with the use of metagenomics and downstream bioinformatics analyses, we characterize the carbohydrate hydrolytic potential of the microbial community, with a special focus on acidosis. To that purpose, carbohydrate-active enzymes were identified, and to further link the community hydrolytic potential with key microbes, bacterial genomes were reconstructed. In addition, we characterized biochemically the specificity and activity of selected enzymes, thus verifying the accuracy of the in silico predictions. The results confirm the retention of the community hydrolytic potential during acidosis and indicate Bacteroidetes to be largely involved in biomass degradation. Bacteroidetes showed higher diversity and genomic content of carbohydrate hydrolytic enzymes that might favor the dominance of this phylum over other bacteria in some anaerobic reactors. The combination of bioinformatic analyses and activity tests enabled us to propose a model of acetylated glucomannan degradation by Bacteroidetes.

IMPORTANCE The enzymatic hydrolysis of lignocellulosic biomass is mainly driven by the action of carbohydrate-active enzymes. By characterizing the gene profiles at the different stages of the anaerobic digestion experiment, we showed that the microbiome retains its hydrolytic functional redundancy even during severe acidosis, despite significant changes in taxonomic composition. By analyzing reconstructed bacterial genomes, we demonstrate that Bacteroidetes hydrolytic gene diversity likely favors the abundance of this phylum in some anaerobic digestion systems. Further, we observe genetic redundancy within the Bacteroidetes group, which accounts for the preserved hydrolytic potential during acidosis. This work also uncovers new polysaccharide utilization loci involved in the deconstruction of various biomasses and proposes the model of acetylated glucomannan degradation by Bacteroidetes. Acetylated glucomannan-enriched biomass is a common substrate for many industries, including pulp and paper production. Using naturally evolved cocktails of enzymes for biomass pretreatment could be an interesting alternative to the commonly used chemical pretreatments.

In marine Bacteroidetes the bulk of glycan degradation during algae blooms is mediated by few clades using a restricted set of genes

We investigated Bacteroidetes during spring algae blooms in the southern North Sea in 2010–2012 using a time series of 38 deeply sequenced metagenomes. Initial partitioning yielded 6455 bins, from which we extracted 3101 metagenome-assembled genomes (MAGs) including 1286 Bacteroidetes MAGs covering ~120 mostly uncultivated species. We identified 13 dominant, recurrent Bacteroidetes clades carrying a restricted set of conserved polysaccharide utilization loci (PULs) that likely mediate the bulk of bacteroidetal algal polysaccharide degradation. The majority of PULs were predicted to target the diatom storage polysaccharide laminarin, alpha-glucans, alpha-mannose-rich substrates, and sulfated xylans. Metaproteomics at 14 selected points in time revealed expression of SusC-like proteins from PULs targeting all of these substrates. Analyses of abundant key players and their PUL repertoires over time furthermore suggested that fewer and simpler polysaccharides dominated early bloom stages, and that more complex polysaccharides became available as blooms progressed.

Investigating Host Microbiota Relationships Through Functional Metagenomics

The human Intestinal mucus is formed by glycoproteins, the O- and N-linked glycans which constitute a crucial source of carbon for commensal gut bacteria, especially when deprived of dietary glycans of plant origin. In recent years, a dozen carbohydrate-active enzymes from cultivated mucin degraders have been characterized. But yet, considering the fact that uncultured species predominate in the human gut microbiota, these biochemical data are far from exhaustive. In this study, we used functional metagenomics to identify new metabolic pathways in uncultured bacteria involved in harvesting mucin glycans. First, we performed a high-throughput screening of a fosmid metagenomic library constructed from the ileum mucosa microbiota using chromogenic substrates. The screening resulted in the isolation of 124 clones producing activities crucial in the degradation of human O- and N-glycans, namely sialidases, β-D-N-acetyl-glucosaminidase, β-D-N-acetyl-galactosaminidase, and/or β-D-mannosidase. Thirteen of these clones were selected based on their diversified functional profiles and were further analyzed on a secondary screening. This step consisted of lectin binding assays to demonstrate the ability of the clones to degrade human intestinal mucus. In total, the structural modification of several mucin motifs, sialylated mucin ones in particular, was evidenced for nine clones. Sequencing their metagenomic loci highlighted complex catabolic pathways involving the complementary functions of glycan sensing, transport, hydrolysis, deacetylation, and deamination, which were sometimes associated with amino acid metabolism machinery. These loci are assigned to several Bacteroides and Feacalibacterium species highly prevalent and abundant in the gut microbiome and explain the metabolic flexibility of gut bacteria feeding both on dietary and human glycans.

Structural and functional analyses of glycoside hydrolase 138 enzymes targeting chain A galacturonic acid in the complex pectin rhamnogalacturonan II

The metabolism of carbohydrate polymers drives microbial diversity in the human gut microbiome. The selection pressures in this environment have spurred the evolution of a complex reservoir of microbial genes encoding carbohydrate-active enzymes (CAZymes). Previously, we have shown that the human gut bacterium Bacteroides thetaiotaomicron (Bt) can depolymerize the most structurally complex glycan, the plant pectin rhamnogalacturonan II (RGII), commonly found in the human diet. Previous investigation of the RGII-degrading apparatus in Bt identified BT0997 as a new CAZyme family, classified as glycoside hydrolase 138 (GH138). The mechanism of substrate recognition by GH138, however, remains unclear. Here, using synthetic substrates and biochemical assays, we show that BT0997 targets the d-galacturonic acid-α-1,2-l-rhamnose linkage in chain A of RGII and that it absolutely requires the presence of a second d-galacturonic acid side chain (linked β-1,3 to l-rhamnose) for activity. NMR analysis revealed that BT0997 operates through a double displacement retaining mechanism. We also report the crystal structure of a BT0997 homolog, BPA0997 from Bacteroides paurosaccharolyticus, in complex with ligands at 1.6 Å resolution. The structure disclosed that the enzyme comprises four domains, including a catalytic TIM (α/β)8 barrel. Characterization of several BT0997 variants identified Glu-294 and Glu-361 as the catalytic acid/base and nucleophile, respectively, and we observed a chloride ion close to the active site. The three-dimensional structure and bioinformatic analysis revealed that two arginines, Arg-332 and Arg-521, are key specificity determinants of BT0997 in targeting d-galacturonic acid residues. In summary, our study reports the first structural and mechanistic analyses of GH138 enzymes.

Bacteroidetes use thousands of enzyme combinations to break down glycans

Unlike proteins, glycan chains are not directly encoded by DNA, but by the specificity of the enzymes that assemble them. Theoretical calculations have proposed an astronomical number of possible isomers (> 10¹² hexasaccharides) but the actual diversity of glycan structures in nature is not known. Bacteria of the Bacteroidetes phylum are considered primary degraders of polysaccharides and they are found in all ecosystems investigated. In Bacteroidetes genomes, carbohydrate-degrading enzymes (CAZymes) are arranged in gene clusters termed polysaccharide utilization loci (PULs). The depolymerization of a given complex glycan by Bacteroidetes PULs requires bespoke enzymes; conversely, the enzyme composition in PULs can provide information on the structure of the targeted glycans. Here we group the 13,537 PULs encoded by 964 Bacteroidetes genomes according to their CAZyme composition. We find that collectively Bacteroidetes have elaborated a few thousand enzyme combinations for glycan breakdown, suggesting a global estimate of diversity of glycan structures much smaller than the theoretical one.

Bacteroidetes genomes contain polysaccharide utilization loci (PULs), each of which encodes enzymes for the breakdown of one particular glycan. By analyzing the enzyme composition of 13,537 PULs, the authors suggest that the natural glycan diversity is orders of magnitude lower than previously proposed.

The agar-specific hydrolase ZgAgaC from the marine bacterium Zobellia galactanivorans defines a new GH16 protein subfamily

Agars are sulfated galactans from red macroalgae and are composed of a d-galactose (G unit) and l-galactose (L unit) alternatively linked by α-1,3 and β-1,4 glycosidic bonds. These polysaccharides display high complexity, with numerous modifications of their backbone (e.g. presence of a 3,6-anhydro-bridge (LA unit) and sulfations and methylation). Currently, bacterial polysaccharidases that hydrolyze agars (β-agarases and β-porphyranases) have been characterized on simple agarose and more rarely on porphyran, a polymer containing both agarobiose (G-LA) and porphyranobiose (GL6S) motifs. How bacteria can degrade complex agars remains therefore an open question. Here, we studied an enzyme from the marine bacterium Zobellia galactanivorans (ZgAgaC) that is distantly related to the glycoside hydrolase 16 (GH16) family β-agarases and β-porphyranases. Using a large red algae collection, we demonstrate that ZgAgaC hydrolyzes not only agarose but also complex agars from Ceramiales species. Using tandem MS analysis, we elucidated the structure of a purified hexasaccharide product, L6S-G-LA2Me-G(2Pentose)-LA2S-G, released by the activity of ZgAgaC on agar extracted from Osmundea pinnatifida By resolving the crystal structure of ZgAgaC at high resolution (1.3 Å) and comparison with the structures of ZgAgaB and ZgPorA in complex with their respective substrates, we determined that ZgAgaC recognizes agarose via a mechanism different from that of classical β-agarases. Moreover, we identified conserved residues involved in the binding of complex oligoagars and demonstrate a probable influence of the acidic polysaccharide's pH microenvironment on hydrolase activity. Finally, a phylogenetic analysis supported the notion that ZgAgaC homologs define a new GH16 subfamily distinct from β-porphyranases and classical β-agarases.

Unraveling the subtleties of β-(1→3)-glucan phosphorylase specificity in the GH94, GH149, and GH161 glycoside hydrolase families

Glycoside phosphorylases (GPs) catalyze the phosphorolysis of glycans into the corresponding sugar 1-phosphates and shortened glycan chains. Given the diversity of natural β-(1→3)-glucans and their wide range of biotechnological applications, the identification of enzymatic tools that can act on β-(1→3)-glucooligosaccharides is an attractive area of research. GP activities acting on β-(1→3)-glucooligosaccharides have been described in bacteria, the photosynthetic excavate Euglena gracilis, and the heterokont Ochromonas spp. Previously, we characterized β-(1→3)-glucan GPs from bacteria and E. gracilis, leading to their classification in glycoside hydrolase family GH149. Here, we characterized GPs from Gram-positive bacteria and heterokont algae acting on β-(1→3)-glucooligosaccharides. We identified a phosphorylase sequence from Ochromonas spp. (OcP1) together with its orthologs from other species, leading us to propose the establishment of a new GH family, designated GH161. To establish the activity of GH161 members, we recombinantly expressed a bacterial GH161 gene sequence (PapP) from the Gram-positive bacterium Paenibacillus polymyxa ATCC 842 in Escherichia coli. We found that PapP acts on β-(1→3)-glucooligosaccharide acceptors with a degree of polymerization (DP) ≥ 2. This activity was distinct from that of characterized GH149 β-(1→3)-glucan phosphorylases, which operate on acceptors with DP ≥ 1. We also found that bacterial GH161 genes co-localize with genes encoding β-glucosidases and ATP-binding cassette transporters, highlighting a probable involvement of GH161 enzymes in carbohydrate degradation. Importantly, in some species, GH161 and GH94 genes were present in tandem, providing evidence that GPs from different CAZy families may work sequentially to degrade oligosaccharides.

Interspecies cross-feeding orchestrates carbon degradation in the rumen ecosystem

Because of their agricultural value, there is a great body of research dedicated to understanding the microorganisms responsible for rumen carbon degradation. However, we lack a holistic view of the microbial food web responsible for carbon processing in this ecosystem. Here, we sampled rumen-fistulated moose, allowing access to rumen microbial communities actively degrading woody plant biomass in real time. We resolved 1,193 viral contigs and 77 unique, near-complete microbial metagenome-assembled genomes, many of which lacked previous metabolic insights. Plant-derived metabolites were measured with NMR and carbohydrate microarrays to quantify the carbon nutrient landscape. Network analyses directly linked measured metabolites to expressed proteins from these unique metagenome-assembled genomes, revealing a genome-resolved three-tiered carbohydrate-fuelled trophic system. This provided a glimpse into microbial specialization into functional guilds defined by specific metabolites. To validate our proteomic inferences, the catalytic activity of a polysaccharide utilization locus from a highly connected metabolic hub genome was confirmed using heterologous gene expression. Viral detected proteins and linkages to microbial hosts demonstrated that phage are active controllers of rumen ecosystem function. Our findings elucidate the microbial and viral members, as well as their metabolic interdependencies, that support in situ carbon degradation in the rumen ecosystem.

A combination of proteomics, metagenome-assembled genomes and heterologous gene expression experiments reveals a trophic system for carbon utilization in the moose rumen microbiome and provides insights into phage dynamics in this ecosystem.

Bacteria in the ageing gut: did the taming of fire promote a long human lifespan?

Unique among animals as they evolved towards Homo sapiens, hominins progressively cooked their food on a routine basis. Cooked products are characterized by singular chemical compounds, derived from the pervasive Maillard reaction. This same reaction is omnipresent in normal metabolism involving carbonyls and amines, and its products accumulate with age. The gut microbiota acts as a first line of defence against the toxicity of cooked Maillard compounds, that also selectively shape the microbial flora, letting specific metabolites to reach the blood stream. Positive selection of metabolic functions allowed the body of hominins who tamed fire to use and dispose of these age-related compounds. I propose here that, as a hopeful accidental consequence, this resulted in extending human lifespan far beyond that of our great ape cousins. The limited data exploring the role of taming fire on the human genetic setup and on its microbiota is discussed in relation with ageing.

dbCAN2: a meta server for automated carbohydrate-active enzyme annotation

Complex carbohydrates of plants are the main food sources of animals and microbes, and serve as promising renewable feedstock for biofuel and biomaterial production. Carbohydrate active enzymes (CAZymes) are the most important enzymes for complex carbohydrate metabolism. With an increasing number of plant and plant-associated microbial genomes and metagenomes being sequenced, there is an urgent need of automatic tools for genomic data mining of CAZymes. We developed the dbCAN web server in 2012 to provide a public service for automated CAZyme annotation for newly sequenced genomes. Here, dbCAN2 (http://cys.bios.niu.edu/dbCAN2) is presented as an updated meta server, which integrates three state-of-the-art tools for CAZome (all CAZymes of a genome) annotation: (i) HMMER search against the dbCAN HMM (hidden Markov model) database; (ii) DIAMOND search against the CAZy pre-annotated CAZyme sequence database and (iii) Hotpep search against the conserved CAZyme short peptide database. Combining the three outputs and removing CAZymes found by only one tool can significantly improve the CAZome annotation accuracy. In addition, dbCAN2 now also accepts nucleotide sequence submission, and offers the service to predict physically linked CAZyme gene clusters (CGCs), which will be a very useful online tool for identifying putative polysaccharide utilization loci (PULs) in microbial genomes or metagenomes.

Ancient acquisition of "alginate utilization loci" by human gut microbiota

In bacteria from the phylum Bacteroidetes, the genes coding for enzymes involved in polysaccharide degradation are often colocalized and coregulated in so-called “polysaccharide utilization loci” (PULs). PULs dedicated to the degradation of marine polysaccharides (e.g. laminaran, ulvan, alginate and porphyran) have been characterized in marine bacteria. Interestingly, the gut microbiome of Japanese individuals acquired, by lateral transfer from marine bacteria, the genes involved in the breakdown of porphyran, the cell wall polysaccharide of the red seaweed used in maki. Sequence similarity analyses predict that the human gut microbiome also encodes enzymes for the degradation of alginate, the main cell wall polysaccharide of brown algae. We undertook the functional characterization of diverse polysaccharide lyases from family PL17, frequently found in marine bacteria as well as those of human gut bacteria. We demonstrate here that this family is polyspecific. Our phylogenetic analysis of family PL17 reveals that all alginate lyases, which have all the same specificity and mode of action, cluster together in a very distinct subfamily. The alginate lyases found in human gut bacteria group together in a single clade which is rooted deeply in the PL17 tree. These enzymes were found in PULs containing PL6 enzymes, which also clustered together in the phylogenetic tree of PL6. Together, biochemical and bioinformatics analyses suggest that acquisition of this system appears ancient and, because only traces of two successful transfers were detected upon inspection of PL6 and PL17 families, the pace of acquisition of marine polysaccharide degradation system is probably very slow.

Cultivation and sequencing of rumen microbiome members from the Hungate1000 Collection

Rumen microbiome biology gets a boost with the release of 410 high-quality reference genomes from the Hungate1000 project.

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B₁₂, ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen.

Identification of Euglena gracilis β-1,3-glucan phosphorylase and establishment of a new glycoside hydrolase (GH) family GH149

Glycoside phosphorylases (EC 2.4.x.x) carry out the reversible phosphorolysis of glucan polymers, producing the corresponding sugar 1-phosphate and a shortened glycan chain. β-1,3-Glucan phosphorylase activities have been reported in the photosynthetic euglenozoan Euglena gracilis, but the cognate protein sequences have not been identified to date. Continuing our efforts to understand the glycobiology of E. gracilis, we identified a candidate phosphorylase sequence, designated EgP1, by proteomic analysis of an enriched cellular protein lysate. We expressed recombinant EgP1 in Escherichia coli and characterized it in vitro as a β-1,3-glucan phosphorylase. BLASTP identified several hundred EgP1 orthologs, most of which were from Gram-negative bacteria and had 37-91% sequence identity to EgP1. We heterologously expressed a bacterial metagenomic sequence, Pro_7066 in E. coli and confirmed it as a β-1,3-glucan phosphorylase, albeit with kinetics parameters distinct from those of EgP1. EgP1, Pro_7066, and their orthologs are classified as a new glycoside hydrolase (GH) family, designated GH149. Comparisons between GH94, EgP1, and Pro_7066 sequences revealed conservation of key amino acids required for the phosphorylase activity, suggesting a phosphorylase mechanism that is conserved between GH94 and GH149. We found bacterial GH149 genes in gene clusters containing sugar transporter and several other GH family genes, suggesting that bacterial GH149 proteins have roles in the degradation of complex carbohydrates. The Bacteroidetes GH149 genes located to previously identified polysaccharide utilization loci, implicated in the degradation of complex carbohydrates. In summary, we have identified a eukaryotic and a bacterial β-1,3-glucan phosphorylase and uncovered a new family of phosphorylases that we name GH149.

The 2018 Nucleic Acids Research database issue and the online molecular biology database collection

The 2018 Nucleic Acids Research Database Issue contains 181 papers spanning molecular biology. Among them, 82 are new and 84 are updates describing resources that appeared in the Issue previously. The remaining 15 cover databases most recently published elsewhere. Databases in the area of nucleic acids include 3DIV for visualisation of data on genome 3D structure and RNArchitecture, a hierarchical classification of RNA families. Protein databases include the established SMART, ELM and MEROPS while GPCRdb and the newcomer STCRDab cover families of biomedical interest. In the area of metabolism, HMDB and Reactome both report new features while PULDB appears in NAR for the first time. This issue also contains reports on genomics resources including Ensembl, the UCSC Genome Browser and ENCODE. Update papers from the IUPHAR/BPS Guide to Pharmacology and DrugBank are highlights of the drug and drug target section while a number of proteomics databases including proteomicsDB are also covered. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been updated, reviewing 138 entries, adding 88 new resources and eliminating 47 discontinued URLs, bringing the current total to 1737 databases. It is available at http://www.oxfordjournals.org/nar/database/c/.