Metatranscriptomics sheds light on the links between the functional traits of fungal guilds and ecological processes in forest soil ecosystems

Soil fungi belonging to different functional guilds, such as saprotrophs, pathogens, and mycorrhizal symbionts, play key roles in forest ecosystems. To date, no study has compared the actual gene expression of these guilds in different forest soils. We used metatranscriptomics to study the competition for organic resources by these fungal groups in boreal, temperate, and Mediterranean forest soils. Using a dedicated mRNA annotation pipeline combined with the JGI MycoCosm database, we compared the transcripts of these three fungal guilds, targeting enzymes involved in C- and N mobilization from plant and microbial cell walls. Genes encoding enzymes involved in the degradation of plant cell walls were expressed at a higher level in saprotrophic fungi than in ectomycorrhizal and pathogenic fungi. However, ectomycorrhizal and saprotrophic fungi showed similarly high expression levels of genes encoding enzymes involved in fungal cell wall degradation. Transcripts for N-related transporters were more highly expressed in ectomycorrhizal fungi than in other groups. We showed that ectomycorrhizal and saprotrophic fungi compete for N in soil organic matter, suggesting that their interactions could decelerate C cycling. Metatranscriptomics provides a unique tool to test controversial ecological hypotheses and to better understand the underlying ecological processes involved in soil functioning and carbon stabilization.

Bioactive glycans in a microbiome-directed food for children with malnutrition

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Stable functional structure despite high taxonomic variability across fungal communities in soils of old-growth montane forests

BACKGROUND

Major advances over the past decade in molecular ecology are providing access to soil fungal diversity in forest ecosystems worldwide, but the diverse functions and metabolic capabilities of this microbial community remain largely elusive. We conducted a field survey in montane old-growth broadleaved and conifer forests, to investigate the relationship between soil fungal diversity and functional genetic traits. To assess the extent to which variation in community composition was associated with dominant tree species (oak, spruce, and fir) and environmental variations in the old-growth forests in the Jade Dragon Snow Mountain in Yunnan Province, we applied rDNA metabarcoding. We also assessed fungal gene expression in soil using mRNA sequencing and specifically assessed the expression of genes related to organic matter decomposition and nutrient acquisition in ectomycorrhizal and saprotrophic fungi.

RESULTS

Our taxonomic profiling revealed striking shifts in the composition of the saprotrophic and ectomycorrhizal guilds among the oak-, fir-, and spruce-dominated forests. The core fungal microbiome comprised only ~ 20% of the total OTUs across all soil samples, although the overlap between conifer-associated communities was substantial. In contrast, seasonality and soil layer explained only a small proportion of the variation in community structure. However, despite their highly variable taxonomic composition, fungal guilds exhibited remarkably similar functional traits for growth-related and core metabolic pathways across forest associations, suggesting ecological redundancy. However, we found that the expression profiles of genes related to polysaccharide and protein degradation and nutrient transport notably varied between and within the fungal guilds, suggesting niche adaptation.

CONCLUSIONS

Overall, our metatranscriptomic analyses revealed the functional potential of soil fungal communities in montane old-growth forests, including a suite of specialized genes and taxa involved in organic matter decomposition. By linking genes to ecological traits, this study provides insights into fungal adaptation strategies to biotic and environmental factors, and sheds light on the importance of understanding functional gene expression patterns in predicting ecosystem functioning. Video Abstract.

Bioactive glycans in a microbiome-directed food for malnourished children

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Designing effective microbiome-directed therapeutic foods to repair these perturbations requires knowledge about how food components interact with the microbiome to alter its expressed functions. Here we use biospecimens from a randomized, controlled trial of a microbiome-directed complementary food prototype (MDCF-2) that produced superior rates of weight gain compared to a conventional ready-to-use supplementary food (RUSF) in 12-18-month-old Bangladeshi children with moderate acute malnutrition (MAM)4. We reconstructed 1000 bacterial genomes (metagenome-assembled genomes, MAGs) present in their fecal microbiomes, identified 75 whose abundances were positively associated with weight gain (change in weight-for-length Z score, WLZ), characterized gene expression changes in these MAGs as a function of treatment type and WLZ response, and used mass spectrometry to quantify carbohydrate structures in MDCF-2 and feces. The results reveal treatment-induced changes in expression of carbohydrate metabolic pathways in WLZ-associated MAGs. Comparing participants consuming MDCF-2 versus RUSF, and MDCF-2-treated children in the upper versus lower quartiles of WLZ responses revealed that two Prevotella copri MAGs positively associated with WLZ were principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilization of its component glycans. Moreover, the predicted specificities of carbohydrate active enzymes expressed by polysaccharide utilization loci (PULs) in these two MAGs correlate with the (i) in vitro growth of Bangladeshi P. copri strains, possessing differing degrees of PUL and overall genomic content similarity to these MAGs, cultured in defined medium containing different purified glycans representative of those in MDCF-2, and (ii) levels of carbohydrate structures identified in feces from clinical trial participants. In the accompanying paper5, we use a gnotobiotic mouse model colonized with age- and WLZ-associated bacterial taxa cultured from this study population, and fed diets resembling those consumed by study participants, to directly test the relationship between P. copri, MDCF-2 glycan metabolism, host ponderal growth responses, and intestinal gene expression and metabolism. The ability to identify bioactive glycan structures in MDCFs that are metabolized by growth-associated bacterial taxa will help guide recommendations about use of this MDCF for children with acute malnutrition representing different geographic locales and ages, as well as enable development of bioequivalent, or more efficacious, formulations composed of culturally acceptable and affordable ingredients.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

Glycoside hydrolase subfamily GH5_57 features a highly redesigned catalytic interface to process complex hetero-β-mannans

Glycoside hydrolase family 5 (GH5) harbors diverse substrate specificities and modes of action, exhibiting notable molecular adaptations to cope with the stereochemical complexity imposed by glycosides and carbohydrates such as cellulose, xyloglucan, mixed-linkage β-glucan, laminarin, (hetero)xylan, (hetero)mannan, galactan, chitosan, N-glycan, rutin and hesperidin. GH5 has been divided into subfamilies, many with higher functional specificity, several of which have not been characterized to date and some that have yet to be discovered with the exploration of sequence/taxonomic diversity. In this work, the current GH5 subfamily inventory is expanded with the discovery of the GH5_57 subfamily by describing an endo-β-mannanase (CapGH5_57) from an uncultured Bacteroidales bacterium recovered from the capybara gut microbiota. Biochemical characterization showed that CapGH5_57 is active on glucomannan, releasing oligosaccharides with a degree of polymerization from 2 to 6, indicating it to be an endo-β-mannanase. The crystal structure, which was solved using single-wavelength anomalous diffraction, revealed a massively redesigned catalytic interface compared with GH5 mannanases. The typical aromatic platforms and the characteristic α-helix-containing β6–α6 loop in the positive-subsite region of GH5_7 mannanases are absent in CapGH5_57, generating a large and open catalytic interface that might favor the binding of branched substrates. Supporting this, CapGH5_57 contains a tryptophan residue adjacent and perpendicular to the cleavage site, indicative of an anchoring site for a substrate with a substitution at the −1 glycosyl moiety. Taken together, these results suggest that despite presenting endo activity on glucomannan, CapGH5_57 may have a new type of substituted heteromannan as its natural substrate. This work demonstrates the still great potential for discoveries regarding the mechanistic and functional diversity of this large and polyspecific GH family by unveiling a novel catalytic interface sculpted to recognize complex heteromannans, which led to the establishment of the GH5_57 subfamily.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Mechanistic insights into consumption of the food additive xanthan gum by the human gut microbiota

Processed foods often include food additives such as xanthan gum, a complex polysaccharide with unique rheological properties, that has established widespread use as a stabilizer and thickening agent. Xanthan gum's chemical structure is distinct from those of host and dietary polysaccharides that are more commonly expected to transit the gastrointestinal tract, and little is known about its direct interaction with the gut microbiota, which plays a central role in digestion of other dietary fibre polysaccharides. Here we show that the ability to digest xanthan gum is common in human gut microbiomes from industrialized countries and appears contingent on a single uncultured bacterium in the family Ruminococcaceae. Our data reveal that this primary degrader cleaves the xanthan gum backbone before processing the released oligosaccharides using additional enzymes. Some individuals harbour Bacteroides intestinalis that is incapable of consuming polymeric xanthan gum but grows on oligosaccharide products generated by the Ruminococcaceae. Feeding xanthan gum to germfree mice colonized with a human microbiota containing the uncultured Ruminococcaceae supports the idea that the additive xanthan gum can drive expansion of the primary degrader Ruminococcaceae, along with exogenously introduced B. intestinalis. Our work demonstrates the existence of a potential xanthan gum food chain involving at least two members of different phyla of gut bacteria and provides an initial framework for understanding how widespread consumption of a recently introduced food additive influences human microbiomes.

Gut microbiome of the largest living rodent harbors unprecedented enzymatic systems to degrade plant polysaccharides

The largest living rodent, capybara, can efficiently depolymerize and utilize lignocellulosic biomass through microbial symbiotic mechanisms yet elusive. Herein, we elucidate the microbial community composition, enzymatic systems and metabolic pathways involved in the conversion of dietary fibers into short-chain fatty acids, a main energy source for the host. In this microbiota, the unconventional enzymatic machinery from Fibrobacteres seems to drive cellulose degradation, whereas a diverse set of carbohydrate-active enzymes from Bacteroidetes, organized in polysaccharide utilization loci, are accounted to tackle complex hemicelluloses typically found in gramineous and aquatic plants. Exploring the genetic potential of this community, we discover a glycoside hydrolase family of β-galactosidases (named as GH173), and a carbohydrate-binding module family (named as CBM89) involved in xylan binding that establishes an unprecedented three-dimensional fold among associated modules to carbohydrate-active enzymes. Together, these results demonstrate how the capybara gut microbiota orchestrates the depolymerization and utilization of plant fibers, representing an untapped reservoir of enzymatic mechanisms to overcome the lignocellulose recalcitrance, a central challenge toward a sustainable and bio-based economy.

Here, Cabral et al., perform a multi-omics analysis of the gut microbiome of capybara, the largest living rodent, unveiling enzymatic mechanisms for the breakdown of lignocellulosic biomass, and report two undescribed families of carbohydrate-active enzymes.

The carbohydrate-active enzyme database: functions and literature

Thirty years have elapsed since the emergence of the classification of carbohydrate-active enzymes in sequence-based families that became the CAZy database over 20 years ago, freely available for browsing and download at www.cazy.org. In the era of large scale sequencing and high-throughput Biology, it is important to examine the position of this specialist database that is deeply rooted in human curation. The three primary tasks of the CAZy curators are (i) to maintain and update the family classification of this class of enzymes, (ii) to classify sequences newly released by GenBank and the Protein Data Bank and (iii) to capture and present functional information for each family. The CAZy website is updated once a month. Here we briefly summarize the increase in novel families and the annotations conducted during the last 8 years. We present several important changes that facilitate taxonomic navigation, and allow to download the entirety of the annotations. Most importantly we highlight the considerable amount of work that accompanies the analysis and report of biochemical data from the literature.

Bioinformatic Analysis of Lytic Polysaccharide Monooxygenases Reveals the Pan-Families Occurrence of Intrinsically Disordered C-Terminal Extensions

Lytic polysaccharide monooxygenases (LPMOs) are monocopper enzymes secreted by many organisms and viruses. LPMOs catalyze the oxidative cleavage of different types of polysaccharides and are today divided into eight families (AA9-11, AA13-17) within the Auxiliary Activity enzyme class of the CAZy database. LPMOs minimal architecture encompasses a catalytic domain, to which can be appended a carbohydrate-binding module. Intriguingly, we observed that some LPMO sequences also display a C-terminal extension of varying length not associated with any known function or fold. Here, we analyzed 27,060 sequences from different LPMO families and show that 60% have a C-terminal extension predicted to be intrinsically disordered. Our analysis shows that these disordered C-terminal regions (dCTRs) are widespread in all LPMO families (except AA13) and differ in terms of sequence length and amino-acid composition. Noteworthily, these dCTRs have so far only been observed in LPMOs. LPMO-dCTRs share a common polyampholytic nature and an enrichment in serine and threonine residues, suggesting that they undergo post-translational modifications. Interestingly, dCTRs from AA11 and AA15 are enriched in redox-sensitive, conditionally disordered regions. The widespread occurrence of dCTRs in LPMOs from evolutionarily very divergent organisms, hints at a possible functional role and opens new prospects in the field of LPMOs.

Evaluating microbiome-directed fibre snacks in gnotobiotic mice and humans

Changing food preferences brought about by westernization that have deleterious health effects1,2-combined with myriad forces that are contributing to increased food insecurity-are catalysing efforts to identify more nutritious and affordable foods3. Consumption of dietary fibre can help to prevent cardiovascular disease, type 2 diabetes and obesity4-6. A substantial number of reports have explored the effects of dietary fibre on the gut microbial community7-9. However, the microbiome is complex, dynamic and exhibits considerable intra- and interpersonal variation in its composition and functions. The large number of potential interactions between the components of the microbiome makes it challenging to define the mechanisms by which food ingredients affect community properties. Here we address the question of how foods containing different fibre preparations can be designed to alter functions associated with specific components of the microbiome. Because a marked increase in snack consumption is associated with westernization, we formulated snack prototypes using plant fibres from different sustainable sources that targeted distinct features of the gut microbiomes of individuals with obesity when transplanted into gnotobiotic mice. We used these snacks to supplement controlled diets that were consumed by adult individuals with obesity or who were overweight. Fibre-specific changes in their microbiomes were linked to changes in their plasma proteomes indicative of an altered physiological state.

Strain-level functional variation in the human gut microbiota based on bacterial binding to artificial food particles

Greater understanding of the spatial relationships between members of the human gut microbiota and available nutrients is needed to gain deeper insights about community dynamics and expressed functions. Therefore, we generated a panel of artificial food particles with each type composed of microscopic paramagnetic beads coated with a fluorescent barcode and one of 60 different dietary or host glycan preparations. Analysis of 160 Bacteroides and Parabacteroides strains disclosed diverse strain-specific and glycan-specific binding phenotypes. We identified carbohydrate structures that correlated with binding by specific bacterial strains in vitro and noted strain-specific differences in the catabolism of glycans that mediate adhesion. Mixed in vitro cultures revealed that these adhesion phenotypes are maintained in more complex communities. Additionally, orally administering glycan beads to gnotobiotic mice confirmed specificity in glycan binding. This approach should facilitate analyses of how strains occupying the same physical niche interact, and it should advance the development of synbiotics, more nutritious foods, and microbiota-based diagnostics.

Characterization of three bacterial glycoside hydrolase family 9 endoglucanases with different modular architectures isolated from a compost metagenome

BACKGROUND

Environmental bacteria express a wide diversity of glycoside hydrolases (GH). Screening and characterization of GH from metagenomic sources provides an insight into biomass degradation strategies of non-cultivated prokaryotes.

METHODS

In the present report, we screened a compost metagenome for lignocellulolytic activities and identified six genes encoding enzymes belonging to family GH9 (GH9a-f). Three of these enzymes (GH9b, GH9d and GH9e) were successfully expressed and characterized.

RESULTS

A phylogenetic analysis of the catalytic domain of pro- and eukaryotic GH9 enzymes suggested the existence of two major subgroups. Bacterial GH9s displayed a wide variety of modular architectures and those harboring an N-terminal Ig-like domain, such as GH9b and GH9d, segregated from the remainder. We purified and characterized GH9 endoglucanases from both subgroups and examined their stabilities, substrate specificities and product profiles. GH9e exhibited an original hydrolysis pattern, liberating an elevated proportion of oligosaccharides longer than cellobiose. All of the enzymes exhibited processive behavior and a synergistic action on crystalline cellulose. Synergy was also evidenced between GH9d and a GH48 enzyme identified from the same metagenome.

CONCLUSIONS

The characterized GH9 enzymes displayed different modular architectures and distinct substrate and product profiles. The presence of a cellulose binding domain was shown to be necessary for binding and digestion of insoluble cellulosic substrates, but not for processivity.

GENERAL SIGNIFICANCE

The identification of six GH9 enzymes from a compost metagenome and the functional variety of three characterized members highlight the importance of this enzyme family in bacterial biomass deconstruction.

α-Galactosidase and Sucrose-Kinase Relationships in a Bi-functional AgaSK Enzyme Produced by the Human Gut Symbiont Ruminococcus gnavus E1

Plant α-galactosides belonging to the raffinose family oligosaccharides (RFOs) and considered as prebiotics, are commonly degraded by α-galactosidases produced by the human gut microbiome. In this environment, the Ruminococcus gnavus E1 symbiont-well-known for various benefit-is able to produce an original RgAgaSK bifunctional enzyme. This enzyme contains an hydrolytic α-galactosidase domain linked to an ATP dependent extra-domain, specifically involved in the α-galactoside hydrolysis and the phosphorylation of the glucose, respectively. However, the multi-modular relationships between both catalytic domains remained hitherto unexplored and has been, consequently, herein investigated. Biochemical characterization of heterologously expressed enzymes either in full-form or in separated domains revealed similar kinetic parameters. These results were supported by molecular modeling studies performed on the whole enzyme in complex with different RFOs. Further enzymatic analysis associated with kinetic degradation of various substrates followed by high pressure anionic exchange chromatography revealed that catalytic efficiency decreased as the number of D-galactosyl moieties branched onto the oligosaccharide increased, suggesting a preference of RgAgaSK for RFO's short chains. A wide prevalence and abundance study on a human metagenomic library showed a high prevalence of the RgAgaSK encoding gene whatever the health status of the individuals. Finally, phylogeny and synteny studies suggested a limited spread by horizontal transfer of the clusters' containing RgAgaSK to only few species of Firmicutes, highlighting the importance of these undispersed tandem activities in the human gut microbiome.

Investigating host-microbiome interactions by droplet based microfluidics

BACKGROUND

Despite the importance of the mucosal interface between microbiota and the host in gut homeostasis, little is known about the mechanisms of bacterial gut colonization, involving foraging for glycans produced by epithelial cells. The slow pace of progress toward understanding the underlying molecular mechanisms is largely due to the lack of efficient discovery tools, especially those targeting the uncultured fraction of the microbiota.

RESULTS

Here, we introduce an ultra-high-throughput metagenomic approach based on droplet microfluidics, to screen fosmid libraries. Thousands of bacterial genomes can be covered in 1 h of work, with less than ten micrograms of substrate. Applied to the screening of the mucosal microbiota for β-N-acetylgalactosaminidase activity, this approach allowed the identification of pathways involved in the degradation of human gangliosides and milk oligosaccharides, the structural homologs of intestinal mucin glycans. These pathways, whose prevalence is associated with inflammatory bowel diseases, could be the result of horizontal gene transfers with Bacteroides species. Such pathways represent novel targets to study the microbiota-host interactions in the context of inflammatory bowel diseases, in which the integrity of the mucosal barrier is impaired.

CONCLUSION

By compartmentalizing experiments inside microfluidic droplets, this method speeds up and miniaturizes by several orders of magnitude the screening process compared to conventional approaches, to capture entire metabolic pathways from metagenomic libraries. The method is compatible with all types of (meta)genomic libraries, and employs a commercially available flow cytometer instead of a custom-made sorting system to detect intracellular or extracellular enzyme activities. This versatile and generic workflow will accelerate experimental exploration campaigns in functional metagenomics and holobiomics studies, to further decipher host-microbiota relationships. Video Abstract.

Analysis of the diversity of the glycoside hydrolase family 130 in mammal gut microbiomes reveals a novel mannoside-phosphorylase function

Mannoside phosphorylases are involved in the intracellular metabolization of mannooligosaccharides, and are also useful enzymes for the in vitro synthesis of oligosaccharides. They are found in glycoside hydrolase family GH130. Here we report on an analysis of 6308 GH130 sequences, including 4714 from the human, bovine, porcine and murine microbiomes. Using sequence similarity networks, we divided the diversity of sequences into 15 mostly isofunctional meta-nodes; of these, 9 contained no experimentally characterized member. By examining the multiple sequence alignments in each meta-node, we predicted the determinants of the phosphorolytic mechanism and linkage specificity. We thus hypothesized that eight uncharacterized meta-nodes would be phosphorylases. These sequences are characterized by the absence of signal peptides and of the catalytic base. Those sequences with the conserved E/K, E/R and Y/R pairs of residues involved in substrate binding would target β-1,2-, β-1,3- and β-1,4-linked mannosyl residues, respectively. These predictions were tested by characterizing members of three of the uncharacterized meta-nodes from gut bacteria. We discovered the first known β-1,4-mannosyl-glucuronic acid phosphorylase, which targets a motif of the Shigella lipopolysaccharide O-antigen. This work uncovers a reliable strategy for the discovery of novel mannoside-phosphorylases, reveals possible interactions between gut bacteria, and identifies a biotechnological tool for the synthesis of antigenic oligosaccharides.

Proteome specialization of anaerobic fungi during ruminal degradation of recalcitrant plant fiber

The rumen harbors a complex microbial mixture of archaea, bacteria, protozoa, and fungi that efficiently breakdown plant biomass and its complex dietary carbohydrates into soluble sugars that can be fermented and subsequently converted into metabolites and nutrients utilized by the host animal. While rumen bacterial populations have been well documented, only a fraction of the rumen eukarya are taxonomically and functionally characterized, despite the recognition that they contribute to the cellulolytic phenotype of the rumen microbiota. To investigate how anaerobic fungi actively engage in digestion of recalcitrant fiber that is resistant to degradation, we resolved genome-centric metaproteome and metatranscriptome datasets generated from switchgrass samples incubated for 48 h in nylon bags within the rumen of cannulated dairy cows. Across a gene catalog covering anaerobic rumen bacteria, fungi and viruses, a significant portion of the detected proteins originated from fungal populations. Intriguingly, the carbohydrate-active enzyme (CAZyme) profile suggested a domain-specific functional specialization, with bacterial populations primarily engaged in the degradation of hemicelluloses, whereas fungi were inferred to target recalcitrant cellulose structures via the detection of a number of endo- and exo-acting enzymes belonging to the glycoside hydrolase (GH) family 5, 6, 8, and 48. Notably, members of the GH48 family were amongst the highest abundant CAZymes and detected representatives from this family also included dockerin domains that are associated with fungal cellulosomes. A eukaryote-selected metatranscriptome further reinforced the contribution of uncultured fungi in the ruminal degradation of recalcitrant fibers. These findings elucidate the intricate networks of in situ recalcitrant fiber deconstruction, and importantly, suggest that the anaerobic rumen fungi contribute a specific set of CAZymes that complement the enzyme repertoire provided by the specialized plant cell wall degrading rumen bacteria.

A catalog of microbial genes from the bovine rumen unveils a specialized and diverse biomass-degrading environment

Background

The rumen microbiota provides essential services to its host and, through its role in ruminant production, contributes to human nutrition and food security. A thorough knowledge of the genetic potential of rumen microbes will provide opportunities for improving the sustainability of ruminant production systems. The availability of gene reference catalogs from gut microbiomes has advanced the understanding of the role of the microbiota in health and disease in humans and other mammals. In this work, we established a catalog of reference prokaryote genes from the bovine rumen.

Results

Using deep metagenome sequencing we identified 13,825,880 non-redundant prokaryote genes from the bovine rumen. Compared to human, pig, and mouse gut metagenome catalogs, the rumen is larger and richer in functions and microbial species associated with the degradation of plant cell wall material and production of methane. Genes encoding enzymes catalyzing the breakdown of plant polysaccharides showed a particularly high richness that is otherwise impossible to infer from available genomes or shallow metagenomics sequencing. The catalog expands the dataset of carbohydrate-degrading enzymes described in the rumen. Using an independent dataset from a group of 77 cattle fed 4 common dietary regimes, we found that only <0.1% of genes were shared by all animals, which contrast with a large overlap for functions, i.e., 63% for KEGG functions. Different diets induced differences in the relative abundance rather than the presence or absence of genes, which explains the great adaptability of cattle to rapidly adjust to dietary changes.

Conclusions

These data bring new insights into functions, carbohydrate-degrading enzymes, and microbes of the rumen to complement the available information on microbial genomes. The catalog is a significant biological resource enabling deeper understanding of phenotypes and biological processes and will be expanded as new data are made available.

Discovery of hyperstable carbohydrate-active enzymes through metagenomics of extreme environments

The enzymes from hyperthermophilic microorganisms populating volcanic sites represent interesting cases of protein adaptation and biotransformations under conditions where conventional enzymes quickly denature. The difficulties in cultivating extremophiles severely limit access to this class of biocatalysts. To circumvent this problem, we embarked on the exploration of the biodiversity of the solfatara Pisciarelli, Agnano (Naples, Italy), to discover hyperthermophilic carbohydrate-active enzymes (CAZymes) and to characterize the entire set of such enzymes in this environment (CAZome). Here, we report the results of the metagenomic analysis of two mud/water pools that greatly differ in both temperature and pH (T = 85 °C and pH 5.5; T = 92 °C and pH 1.5, for Pool1 and Pool2, respectively). DNA deep sequencing and following in silico analysis led to 14 934 and 17 652 complete ORFs in Pool1 and Pool2, respectively. They exclusively belonged to archaeal cells and viruses with great genera variance within the phylum Crenarchaeota, which reflected the difference in temperature and pH of the two Pools. Surprisingly, 30% and 62% of all of the reads obtained from Pool1 and 2, respectively, had no match in nucleotide databanks. Genes associated with carbohydrate metabolism were 15% and 16% of the total in the two Pools, with 278 and 308 putative CAZymes in Pool1 and 2, corresponding to ~ 2.0% of all ORFs. Biochemical characterization of two CAZymes of a previously unknown archaeon revealed a novel subfamily GH5_19 β-mannanase/β-1,3-glucanase whose hemicellulose specificity correlates with the vegetation surrounding the sampling site, and a novel NAD⁺ -dependent GH109 with a previously unreported β-N-acetylglucosaminide/β-glucoside specificity. DATABASES: The sequencing reads are available in the NCBI Sequence Read Archive (SRA) database under the accession numbers SRR7545549 (Pool1) and SRR7545550 (Pool2). The sequences of GH5_Pool2 and GH109_Pool2 are available in GenBank database under the accession numbers MK869723 and MK86972, respectively. The environmental data relative to Pool1 and Pool2 (NCBI BioProject PRJNA481947) are available in the Biosamples database under the accession numbers SAMN09692669 (Pool1) and SAMN09692670 (Pool2).

A subfamily roadmap of the evolutionarily diverse glycoside hydrolase family 16 (GH16)

Glycoside hydrolase family (GH) 16 comprises a large and taxonomically diverse family of glycosidases and transglycosidases that adopt a common β-jelly-roll fold and are active on a range of terrestrial and marine polysaccharides. Presently, broadly insightful sequence–function correlations in GH16 are hindered by a lack of a systematic subfamily structure. To fill this gap, we have used a highly scalable protein sequence similarity network analysis to delineate nearly 23,000 GH16 sequences into 23 robust subfamilies, which are strongly supported by hidden Markov model and maximum likelihood molecular phylogenetic analyses. Subsequent evaluation of over 40 experimental three-dimensional structures has highlighted key tertiary structural differences, predominantly manifested in active-site loops, that dictate substrate specificity across the GH16 evolutionary landscape. As for other large GH families (i.e. GH5, GH13, and GH43), this new subfamily classification provides a roadmap for functional glycogenomics that will guide future bioinformatics and experimental structure–function analyses. The GH16 subfamily classification is publicly available in the CAZy database. The sequence similarity network workflow used here, SSNpipe, is freely available from GitHub.

PULDB: the expanded database of Polysaccharide Utilization Loci

The Polysaccharide Utilization Loci (PUL) database was launched in 2015 to present PUL predictions in ∼70 Bacteroidetes species isolated from the human gastrointestinal tract, as well as PULs derived from the experimental data reported in the literature. In 2018 PULDB offers access to 820 genomes, sampled from various environments and covering a much wider taxonomical range. A Krona dynamic chart was set up to facilitate browsing through taxonomy. Literature surveys now allows the presentation of the most recent (i) PUL repertoires deduced from RNAseq large-scale experiments, (ii) PULs that have been subjected to in-depth biochemical analysis and (iii) new Carbohydrate-Active enzyme (CAZyme) families that contributed to the refinement of PUL predictions. To improve PUL visualization and genome browsing, the previous annotation of genes encoding CAZymes, regulators, integrases and SusCD has now been expanded to include functionally relevant protein families whose genes are significantly found in the vicinity of PULs: sulfatases, proteases, ROK repressors, epimerases and ATP-Binding Cassette and Major Facilitator Superfamily transporters. To cope with cases where susCD may be absent due to incomplete assemblies/split PULs, we present ‘CAZyme cluster’ predictions. Finally, a PUL alignment tool, operating on the tagged families instead of amino-acid sequences, was integrated to retrieve PULs similar to a query of interest. The updated PULDB website is accessible at www.cazy.org/PULDB_new/

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.

A metagenomics investigation of carbohydrate-active enzymes along the goat and camel intestinal tract

Studies of the digestive microbiota of ruminant animals most often focus on the bacterial diversity in the rumen or the feces of the animals, but little is known about the diversity and functions of their distal intestine. Here, the bacterial microbiota of the distal intestinal tract of two goats and two camels was investigated by metagenomics techniques. The bacterial taxonomic diversity and carbohydrate-active enzyme profile were estimated for samples taken from the small intestine, the large intestine, and the rectum of each animal. The bacterial diversity and abundance in the small intestine were lower than in the rectal and large intestinal samples. Analysis of the carbohydrate-active enzyme profiles at each site revealed a comparatively low abundance of enzymes targeting xylan and cellulose in all animals examined, similar to what has been reported earlier for sheep and therefore suggesting that plant cell wall digestion probably takes place elsewhere, such as in the rumen.

Development and characterization of stable anaerobic thermophilic methanogenic microbiomes fermenting switchgrass at decreasing residence times

BACKGROUND

Anaerobic fermentation of lignocellulose occurs in both natural and managed environments, and is an essential part of the carbon cycle as well as a promising route to sustainable production of fuels and chemicals. Lignocellulose solubilization by mixed microbiomes is important in these contexts.

RESULTS

Here, we report the development of stable switchgrass-fermenting enrichment cultures maintained at various residence times and moderately high (55 °C) temperatures. Anaerobic microbiomes derived from a digester inoculum were incubated at 55 °C and fed semi-continuously with medium containing 30 g/L mid-season harvested switchgrass to achieve residence times (RT) of 20, 10, 5, and 3.3 days. Stable, time-invariant cellulolytic methanogenic cultures with minimal accumulation of organic acids were achieved for all RTs. Fractional carbohydrate solubilization was 0.711, 0.654, 0.581 and 0.538 at RT = 20, 10, 5 and 3.3 days, respectively, and glucan solubilization was proportional to xylan solubilization at all RTs. The rate of solubilization was described well by the equation r = k(C - C₀f_r), where C represents the concentration of unutilized carbohydrate, C₀ is the concentration of carbohydrate (cellulose and hemicellulose) entering the bioreactor and f_r is the extrapolated fraction of entering carbohydrate that is recalcitrant at infinite residence time. The 3.3 day RT is among the shortest RT reported for stable thermophilic, methanogenic digestion of a lignocellulosic feedstock. 16S rDNA phylotyping and metagenomic analyses were conducted to characterize the effect of RT on community dynamics and to infer functional roles in the switchgrass to biogas conversion to the various microbial taxa. Firmicutes were the dominant phylum, increasing in relative abundance from 54 to 96% as RT decreased. A Clostridium clariflavum strain with genetic markers for xylose metabolism was the most abundant lignocellulose-solubilizing bacterium. A Thermotogae (Defluviitoga tunisiensis) was the most abundant bacterium in switchgrass digesters at RT = 20 days but decreased in abundance at lower RTs as did multiple Chloroflexi. Synergistetes and Euryarchaeota were present at roughly constant levels over the range of RTs examined.

CONCLUSIONS

A system was developed in which stable methanogenic steady-states were readily obtained with a particulate biomass feedstock, mid-season switchgrass, at laboratory (1 L) scale. Characterization of the extent and rate of carbohydrate solubilization in combination with 16S rDNA and metagenomic sequencing provides a multi-dimensional view of performance, species composition, glycoside hydrolases, and metabolic function with varying residence time. These results provide a point of reference and guidance for future studies and organism development efforts involving defined cultures.

A model species for agricultural pest genomics: the genome of the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

The Colorado potato beetle is one of the most challenging agricultural pests to manage. It has shown a spectacular ability to adapt to a variety of solanaceaeous plants and variable climates during its global invasion, and, notably, to rapidly evolve insecticide resistance. To examine evidence of rapid evolutionary change, and to understand the genetic basis of herbivory and insecticide resistance, we tested for structural and functional genomic changes relative to other arthropod species using genome sequencing, transcriptomics, and community annotation. Two factors that might facilitate rapid evolutionary change include transposable elements, which comprise at least 17% of the genome and are rapidly evolving compared to other Coleoptera, and high levels of nucleotide diversity in rapidly growing pest populations. Adaptations to plant feeding are evident in gene expansions and differential expression of digestive enzymes in gut tissues, as well as expansions of gustatory receptors for bitter tasting. Surprisingly, the suite of genes involved in insecticide resistance is similar to other beetles. Finally, duplications in the RNAi pathway might explain why Leptinotarsa decemlineata has high sensitivity to dsRNA. The L. decemlineata genome provides opportunities to investigate a broad range of phenotypes and to develop sustainable methods to control this widely successful pest.

Comparative Metagenomics of Cellulose- and Poplar Hydrolysate-Degrading Microcosms from Gut Microflora of the Canadian Beaver (Castor canadensis) and North American Moose (Alces americanus) after Long-Term Enrichment

To identify carbohydrate-active enzymes (CAZymes) that might be particularly relevant for wood fiber processing, we performed a comparative metagenomic analysis of digestive systems from Canadian beaver (Castor canadensis) and North American moose (Alces americanus) following 3 years of enrichment on either microcrystalline cellulose or poplar hydrolysate. In total, 9,386 genes encoding CAZymes and carbohydrate-binding modules (CBMs) were identified, with up to half predicted to originate from Firmicutes, Bacteroidetes, Chloroflexi, and Proteobacteria phyla, and up to 17% from unknown phyla. Both PCA and hierarchical cluster analysis distinguished the annotated glycoside hydrolase (GH) distributions identified herein, from those previously reported for grass-feeding mammals and herbivorous foragers. The CAZyme profile of moose rumen enrichments also differed from a recently reported moose rumen metagenome, most notably by the absence of GH13-appended dockerins. Consistent with substrate-driven convergence, CAZyme profiles from both poplar hydrolysate-fed cultures differed from cellulose-fed cultures, most notably by increased numbers of unique sequences belonging to families GH3, GH5, GH43, GH53, and CE1. Moreover, pairwise comparisons of moose rumen enrichments further revealed higher counts of GH127 and CE15 families in cultures fed with poplar hydrolysate. To expand our scope to lesser known carbohydrate-active proteins, we identified and compared multi-domain proteins comprising both a CBM and domain of unknown function (DUF) as well as proteins with unknown function within the 416 predicted polysaccharide utilization loci (PULs). Interestingly, DUF362, identified in iron-sulfur proteins, was consistently appended to CBM9; on the other hand, proteins with unknown function from PULs shared little identity unless from identical PULs. Overall, this study sheds new light on the lignocellulose degrading capabilities of microbes originating from digestive systems of mammals known for fiber-rich diets, and highlights the value of enrichment to select new CAZymes from metagenome sequences for future biochemical characterization.

Feed in summer, rest in winter: microbial carbon utilization in forest topsoil

BACKGROUND

Evergreen coniferous forests contain high stocks of organic matter. Significant carbon transformations occur in litter and soil of these ecosystems, making them important for the global carbon cycle. Due to seasonal allocation of photosynthates to roots, carbon availability changes seasonally in the topsoil. The aim of this paper was to describe the seasonal differences in C source utilization and the involvement of various members of soil microbiome in this process.

RESULTS

Here, we show that microorganisms in topsoil encode a diverse set of carbohydrate-active enzymes, including glycoside hydrolases and auxiliary enzymes. While the transcription of genes encoding enzymes degrading reserve compounds, such as starch or trehalose, was high in soil in winter, summer was characterized by high transcription of ligninolytic and cellulolytic enzymes produced mainly by fungi. Fungi strongly dominated the transcription in litter and an equal contribution of bacteria and fungi was found in soil. The turnover of fungal biomass appeared to be faster in summer than in winter, due to high activity of enzymes targeting its degradation, indicating fast growth in both litter and soil. In each enzyme family, hundreds to thousands of genes were typically transcribed simultaneously.

CONCLUSIONS

Seasonal differences in the transcription of glycoside hydrolases and auxiliary enzyme genes are more pronounced in soil than in litter. Our results suggest that mainly fungi are involved in decomposition of recalcitrant biopolymers in summer, while bacteria replace them in this role in winter. Transcripts of genes encoding enzymes targeting plant biomass biopolymers, reserve compounds and fungal cell walls were especially abundant in the coniferous forest topsoil.

Ninety-nine de novo assembled genomes from the moose (Alces alces) rumen microbiome provide new insights into microbial plant biomass degradation

The moose (Alces alces) is a ruminant that harvests energy from fiber-rich lignocellulose material through carbohydrate-active enzymes (CAZymes) produced by its rumen microbes. We applied shotgun metagenomics to rumen contents from six moose to obtain insights into this microbiome. Following binning, 99 metagenome-assembled genomes (MAGs) belonging to 11 prokaryotic phyla were reconstructed and characterized based on phylogeny and CAZyme profile. The taxonomy of these MAGs reflected the overall composition of the metagenome, with dominance of the phyla Bacteroidetes and Firmicutes. Unlike in other ruminants, Spirochaetes constituted a significant proportion of the community and our analyses indicate that the corresponding strains are primarily pectin digesters. Pectin-degrading genes were also common in MAGs of Ruminococcus, Fibrobacteres and Bacteroidetes and were overall overrepresented in the moose microbiome compared with other ruminants. Phylogenomic analyses revealed several clades within the Bacteriodetes without previously characterized genomes. Several of these MAGs encoded a large numbers of dockerins, a module usually associated with cellulosomes. The Bacteroidetes dockerins were often linked to CAZymes and sometimes encoded inside polysaccharide utilization loci, which has never been reported before. The almost 100 CAZyme-annotated genomes reconstructed in this study provide an in-depth view of an efficient lignocellulose-degrading microbiome and prospects for developing enzyme technology for biorefineries.

How members of the human gut microbiota overcome the sulfation problem posed by glycosaminoglycans

The human microbiota, which plays an important role in health and disease, uses complex carbohydrates as a major source of nutrients. Utilization hierarchy indicates that the host glycosaminoglycans heparin (Hep) and heparan sulfate (HS) are high-priority carbohydrates for Bacteroides thetaiotaomicron, a prominent member of the human microbiota. The sulfation patterns of these glycosaminoglycans are highly variable, which presents a significant enzymatic challenge to the polysaccharide lyases and sulfatases that mediate degradation. It is possible that the bacterium recruits lyases with highly plastic specificities and expresses a repertoire of enzymes that target substructures of the glycosaminoglycans with variable sulfation or that the glycans are desulfated before cleavage by the lyases. To distinguish between these mechanisms, the components of the B. thetaiotaomicron Hep/HS degrading apparatus were analyzed. The data showed that the bacterium expressed a single-surface endo-acting lyase that cleaved HS, reflecting its higher molecular weight compared with Hep. Both Hep and HS oligosaccharides imported into the periplasm were degraded by a repertoire of lyases, with each enzyme displaying specificity for substructures within these glycosaminoglycans that display a different degree of sulfation. Furthermore, the crystal structures of a key surface glycan binding protein, which is able to bind both Hep and HS, and periplasmic sulfatases reveal the major specificity determinants for these proteins. The locus described here is highly conserved within the human gut Bacteroides, indicating that the model developed is of generic relevance to this important microbial community.

A Metagenomics Investigation of Carbohydrate-Active Enzymes along the Gastrointestinal Tract of Saudi Sheep

The digestive microbiota of humans and of a wide range of animals has recently become amenable to in-depth studies due to the emergence of DNA-based metagenomic techniques that do not require cultivation of gut microbes. These techniques are now commonly used to explore the feces of humans and animals under the assumption that such samples are faithful proxies for the intestinal microbiota. Sheep (Ovis aries) are ruminant animals particularly adapted to life in arid regions and in particular Najdi, Noaimi (Awassi), and Harrei (Harri) breeds that are raised in Saudi Arabia for milk and/or meat production. Here we report a metagenomics investigation of the distal digestive tract of one animal from each breed that (i) examines the microbiota at three intestinal subsites (small intestine, mid-colon, and rectum), (ii) performs an in-depth analysis of the carbohydrate-active enzymes genes encoded by the microbiota at the three subsites, and (iii) compares the microbiota and carbohydrate-active enzyme profile at the three subsites across the different breeds. For all animals we found that the small intestine is characterized by a lower taxonomic diversity than that of the large intestine and of the rectal samples. Mirroring this observation, we also find that the spectrum of encoded carbohydrate-active enzymes of the mid-colon and rectal sites is much richer than that of the small intestine. However, the number of encoded cellulases and xylanases in the various intestinal subsites was found to be surprisingly low, indicating that the bulk of the fiber digestion is performed upstream in the rumen, and that the carbon source for the intestinal flora is probably constituted of the rumen fungi and bacteria that pass in the intestines. In consequence we argue that ruminant feces, which are often analyzed for the search of microbial genes involved in plant cell wall degradation, are probably a poor proxy for the lignocellulolytic potential of the host.

Discovery of genes coding for carbohydrate-active enzyme by metagenomic analysis of lignocellulosic biomasses

In this study, a high-throughput sequencing approach was applied to discover novel biocatalysts for lignocellulose hydrolysis from three dedicated energy crops, Arundo donax, Eucalyptus camaldulensis and Populus nigra, after natural biodegradation. The microbiomes of the three lignocellulosic biomasses were dominated by bacterial species (approximately 90%) with the highest representation by the Streptomyces genus both in the total microbial community composition and in the microbial diversity related to GH families of predicted ORFs. Moreover, the functional clustering of the predicted ORFs showed a prevalence of poorly characterized genes, suggesting these lignocellulosic biomasses are potential sources of as yet unknown genes. 1.2%, 0.6% and 3.4% of the total ORFs detected in A. donax, E. camaldulensis and P. nigra, respectively, were putative Carbohydrate-Active Enzymes (CAZymes). Interestingly, the glycoside hydrolases abundance in P. nigra (1.8%) was higher than that detected in the other biomasses investigated in this study. Moreover, a high percentage of (hemi)cellulases with different activities and accessory enzymes (mannanases, polygalacturonases and feruloyl esterases) was detected, confirming that the three analyzed samples were a reservoir of diversified biocatalysts required for an effective lignocellulose saccharification.

Broad phylogeny and functionality of cellulosomal components in the bovine rumen microbiome

The cellulosome is an extracellular multi‐enzyme complex that is considered one of the most efficient plant cell wall‐degrading strategies devised by nature. Its unique modular architecture, achieved by high affinity and specific interaction between protein modules (cohesins and dockerins) enables formation of various enzyme combinations. Extensive research has been dedicated to the mechanistic nature of the cellulosome complex. Nevertheless, little is known regarding its distribution and abundance among microbes in natural plant fibre‐rich environments. Here, we explored these questions in bovine rumen microbial communities, specialized in efficient degradation of lignocellulosic plant material. We bioinformatically screened for cellulosomal modules in this complex environment using a previously published ultra‐deep fibre‐adherent rumen metagenome. Intriguingly, a large portion of the functions of the dockerin‐containing proteins were related to alternative biological processes, and not necessarily to the classic fibre degradation function. Our analysis was experimentally validated by characterizing specific interactions between selected cohesins and dockerins and revealed that cellulosome is a more generalized strategy used by diverse bacteria, some of which were not previously associated with cellulosome production. Remarkably, our results provide additional proof of similarity among rumen microbial communities worldwide. This study suggests a broader and widespread role for the cellulosomal machinery in nature.

Increased Gut Redox and Depletion of Anaerobic and Methanogenic Prokaryotes in Severe Acute Malnutrition

Severe acute malnutrition (SAM) is associated with inadequate diet, low levels of plasma antioxidants and gut microbiota alterations. The link between gut redox and microbial alterations, however, remains unexplored. By sequencing the gut microbiomes of 79 children of varying nutritional status from three centers in Senegal and Niger, we found a dramatic depletion of obligate anaerobes in malnutrition. This was confirmed in an individual patient data meta-analysis including 107 cases and 77 controls from 5 different African and Asian countries. Specifically, several species of the Bacteroidaceae, Eubacteriaceae, Lachnospiraceae and Ruminococceae families were consistently depleted while Enterococcus faecalis, Escherichia coli and Staphylococcus aureus were consistently enriched. Further analyses on our samples revealed increased fecal redox potential, decreased total bacterial number and dramatic Methanobrevibacter smithii depletion. Indeed, M. smithii was detected in more than half of the controls but in none of the cases. No causality was demonstrated but, based on our results, we propose a unifying theory linking microbiota specificity, lacking anaerobes and archaea, to low antioxidant nutrients, and lower food conversion.

Structural and biochemical characterization of the β-N-acetylglucosaminidase from Thermotoga maritima: toward rationalization of mechanistic knowledge in the GH73 family

Members of the GH73 glycosidase family cleave the β-1,4-glycosidic bond between the N-acetylglucosaminyl (GlcNAc) and N-acetylmuramyl (MurNAc) moieties in bacterial peptidoglycan. A catalytic mechanism has been proposed for members FlgJ, Auto, AcmA and Atl(WM) and the structural analysis of FlgJ and Auto revealed a conserved α/β fold reminiscent of the distantly related GH23 lysozyme. Comparison of the active site residues reveals variability in the nature of the catalytic general base suggesting two distinct catalytic mechanisms: an inverting mechanism involving two distant glutamate residues and a substrate-assisted mechanism involving anchimeric assistance by the C2-acetamido group of the GlcNAc moiety. Herein, we present the biochemical characterization and crystal structure of TM0633 from the hyperthermophilic bacterium Thermotoga maritima. TM0633 adopts the α/β fold of the family and displays β-N-acetylglucosaminidase activity on intact peptidoglycan sacculi. Site-directed mutagenesis identifies Glu34, Glu65 and Tyr118 as important residues for catalysis. A thorough bioinformatic analysis of the GH73 sequences identified five phylogenetic clusters. TM0633, FlgJ and Auto belong to a group of three clusters that conserve two carboxylate residues involved in a classical inverting acid-base mechanism. Members of the other two clusters lack a conserved catalytic general base supporting a substrate-assisted mechanism. Molecular modeling of representative members from each cluster suggests that variability in length of the β-hairpin region above the active site confers ligand-binding specificity and modulates the catalytic mechanisms within the GH73 family.

Extensive sampling of basidiomycete genomes demonstrates inadequacy of the white-rot/brown-rot paradigm for wood decay fungi

Basidiomycota (basidiomycetes) make up 32% of the described fungi and include most wood-decaying species, as well as pathogens and mutualistic symbionts. Wood-decaying basidiomycetes have typically been classified as either white rot or brown rot, based on the ability (in white rot only) to degrade lignin along with cellulose and hemicellulose. Prior genomic comparisons suggested that the two decay modes can be distinguished based on the presence or absence of ligninolytic class II peroxidases (PODs), as well as the abundance of enzymes acting directly on crystalline cellulose (reduced in brown rot). To assess the generality of the white-rot/brown-rot classification paradigm, we compared the genomes of 33 basidiomycetes, including four newly sequenced wood decayers, and performed phylogenetically informed principal-components analysis (PCA) of a broad range of gene families encoding plant biomass-degrading enzymes. The newly sequenced Botryobasidium botryosum and Jaapia argillacea genomes lack PODs but possess diverse enzymes acting on crystalline cellulose, and they group close to the model white-rot species Phanerochaete chrysosporium in the PCA. Furthermore, laboratory assays showed that both B. botryosum and J. argillacea can degrade all polymeric components of woody plant cell walls, a characteristic of white rot. We also found expansions in reducing polyketide synthase genes specific to the brown-rot fungi. Our results suggest a continuum rather than a dichotomy between the white-rot and brown-rot modes of wood decay. A more nuanced categorization of rot types is needed, based on an improved understanding of the genomics and biochemistry of wood decay.

The genome of the white-rot fungus Pycnoporus cinnabarinus: a basidiomycete model with a versatile arsenal for lignocellulosic biomass breakdown

BACKGROUND

Saprophytic filamentous fungi are ubiquitous micro-organisms that play an essential role in photosynthetic carbon recycling. The wood-decayer Pycnoporus cinnabarinus is a model fungus for the study of plant cell wall decomposition and is used for a number of applications in green and white biotechnology.

RESULTS

The 33.6 megabase genome of P. cinnabarinus was sequenced and assembled, and the 10,442 predicted genes were functionally annotated using a phylogenomic procedure. In-depth analyses were carried out for the numerous enzyme families involved in lignocellulosic biomass breakdown, for protein secretion and glycosylation pathways, and for mating type. The P. cinnabarinus genome sequence revealed a consistent repertoire of genes shared with wood-decaying basidiomycetes. P. cinnabarinus is thus fully equipped with the classical families involved in cellulose and hemicellulose degradation, whereas its pectinolytic repertoire appears relatively limited. In addition, P. cinnabarinus possesses a complete versatile enzymatic arsenal for lignin breakdown. We identified several genes encoding members of the three ligninolytic peroxidase types, namely lignin peroxidase, manganese peroxidase and versatile peroxidase. Comparative genome analyses were performed in fungi displaying different nutritional strategies (white-rot and brown-rot modes of decay). P. cinnabarinus presents a typical distribution of all the specific families found in the white-rot life style. Growth profiling of P. cinnabarinus was performed on 35 carbon sources including simple and complex substrates to study substrate utilization and preferences. P. cinnabarinus grew faster on crude plant substrates than on pure, mono- or polysaccharide substrates. Finally, proteomic analyses were conducted from liquid and solid-state fermentation to analyze the composition of the secretomes corresponding to growth on different substrates. The distribution of lignocellulolytic enzymes in the secretomes was strongly dependent on growth conditions, especially for lytic polysaccharide mono-oxygenases.

CONCLUSIONS

With its available genome sequence, P. cinnabarinus is now an outstanding model system for the study of the enzyme machinery involved in the degradation or transformation of lignocellulosic biomass.

The carbohydrate-active enzymes database (CAZy) in 2013

The Carbohydrate-Active Enzymes database (CAZy; http://www.cazy.org) provides online and continuously updated access to a sequence-based family classification linking the sequence to the specificity and 3D structure of the enzymes that assemble, modify and breakdown oligo- and polysaccharides. Functional and 3D structural information is added and curated on a regular basis based on the available literature. In addition to the use of the database by enzymologists seeking curated information on CAZymes, the dissemination of a stable nomenclature for these enzymes is probably a major contribution of CAZy. The past few years have seen the expansion of the CAZy classification scheme to new families, the development of subfamilies in several families and the power of CAZy for the analysis of genomes and metagenomes. This article outlines the changes that have occurred in CAZy during the past 5 years and presents our novel effort to display the resolution and the carbohydrate ligands in crystallographic complexes of CAZymes.

Gut microbiota from twins discordant for obesity modulate metabolism in mice

The role of specific gut microbes in shaping body composition remains unclear. We transplanted fecal microbiota from adult female twin pairs discordant for obesity into germ-free mice fed low-fat mouse chow, as well as diets representing different levels of saturated fat and fruit and vegetable consumption typical of the U.S. diet. Increased total body and fat mass, as well as obesity-associated metabolic phenotypes, were transmissible with uncultured fecal communities and with their corresponding fecal bacterial culture collections. Cohousing mice harboring an obese twin's microbiota (Ob) with mice containing the lean co-twin's microbiota (Ln) prevented the development of increased body mass and obesity-associated metabolic phenotypes in Ob cage mates. Rescue correlated with invasion of specific members of Bacteroidetes from the Ln microbiota into Ob microbiota and was diet-dependent. These findings reveal transmissible, rapid, and modifiable effects of diet-by-microbiota interactions.

Role of glycoside phosphorylases in mannose foraging by human gut bacteria

To metabolize both dietary fiber constituent carbohydrates and host glycans lining the intestinal epithelium, gut bacteria produce a wide range of carbohydrate-active enzymes, of which glycoside hydrolases are the main components. In this study, we describe the ability of phosphorylases to participate in the breakdown of human N-glycans, from an analysis of the substrate specificity of UhgbMP, a mannoside phosphorylase of the GH130 protein family discovered by functional metagenomics. UhgbMP is found to phosphorolyze β-D-Manp-1,4-β-D-GlcpNAc-1,4-D-GlcpNAc and is also a highly efficient enzyme to catalyze the synthesis of this precious N-glycan core oligosaccharide by reverse phosphorolysis. Analysis of sequence conservation within family GH130, mapped on a three-dimensional model of UhgbMP and supported by site-directed mutagenesis results, revealed two GH130 subfamilies and allowed the identification of key residues responsible for catalysis and substrate specificity. The analysis of the genomic context of 65 known GH130 sequences belonging to human gut bacteria indicates that the enzymes of the GH130_1 subfamily would be involved in mannan catabolism, whereas the enzymes belonging to the GH130_2 subfamily would rather work in synergy with glycoside hydrolases of the GH92 and GH18 families in the breakdown of N-glycans. The use of GH130 inhibitors as therapeutic agents or functional foods could thus be considered as an innovative strategy to inhibit N-glycan degradation, with the ultimate goal of protecting, or restoring, the epithelial barrier.

Dissecting partner recognition by an intrinsically disordered protein using descriptive random mutagenesis

In view of getting insights into the molecular determinants of the binding efficiency of intrinsically disordered proteins (IDPs), we used random mutagenesis. As a proof of concept, we chose the interaction between the intrinsically disordered C-terminal domain of the measles virus nucleoprotein (NTAIL) and the X domain (XD) of the viral phosphoprotein and assessed how amino acid substitutions introduced at random within NTAIL affect partner recognition. In contrast with directed evolution approaches, we did not apply any selection and used the gene library approach not for production purposes but for achieving a better understanding of the NTAIL/XD interaction. For that reason, and to differentiate our approach from similar approaches that make use of systematic (i.e., targeted) mutagenesis, we propose to call it "descriptive random mutagenesis" (DRM). NTAIL variants generated by error-prone PCR were picked at random in the absence of selection pressure and were characterized in terms of sequence and binding abilities toward XD. DRM not only identified determinants of NTAIL/XD interaction that were in good agreement with previous work but also provided new insights. In particular, we discovered that the primary interaction site is poorly evolvable in terms of binding abilities toward XD. We also identified a critical NTAIL residue whose role in stabilizing the NTAIL/XD complex had previously escaped detection, and we identified NTAIL regulatory sites that dampen the interaction while being located outside the primary interaction site. Results show that DRM is a valuable approach to study binding abilities of IDPs.

Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes

BACKGROUND

Since its inception, the carbohydrate-active enzymes database (CAZy; http://www.cazy.org) has described the families of enzymes that cleave or build complex carbohydrates, namely the glycoside hydrolases (GH), the polysaccharide lyases (PL), the carbohydrate esterases (CE), the glycosyltransferases (GT) and their appended non-catalytic carbohydrate-binding modules (CBM). The recent discovery that members of families CBM33 and family GH61 are in fact lytic polysaccharide monooxygenases (LPMO), demands a reclassification of these families into a suitable category.

RESULTS

Because lignin is invariably found together with polysaccharides in the plant cell wall and because lignin fragments are likely to act in concert with (LPMO), we have decided to join the families of lignin degradation enzymes to the LPMO families and launch a new CAZy class that we name "Auxiliary Activities" in order to accommodate a range of enzyme mechanisms and substrates related to lignocellulose conversion. Comparative analyses of these auxiliary activities in 41 fungal genomes reveal a pertinent division of several fungal groups and subgroups combining their phylogenetic origin and their nutritional mode (white vs. brown rot).

CONCLUSIONS

The new class introduced in the CAZy database extends the traditional CAZy families, and provides a better coverage of the full extent of the lignocellulose breakdown machinery.

The Paleozoic origin of enzymatic lignin decomposition reconstructed from 31 fungal genomes

Wood is a major pool of organic carbon that is highly resistant to decay, owing largely to the presence of lignin. The only organisms capable of substantial lignin decay are white rot fungi in the Agaricomycetes, which also contains non-lignin-degrading brown rot and ectomycorrhizal species. Comparative analyses of 31 fungal genomes (12 generated for this study) suggest that lignin-degrading peroxidases expanded in the lineage leading to the ancestor of the Agaricomycetes, which is reconstructed as a white rot species, and then contracted in parallel lineages leading to brown rot and mycorrhizal species. Molecular clock analyses suggest that the origin of lignin degradation might have coincided with the sharp decrease in the rate of organic carbon burial around the end of the Carboniferous period.

Complex carbohydrate utilization by the healthy human microbiome

The various ecological habitats in the human body provide microbes a wide array of nutrient sources and survival challenges. Advances in technology such as DNA sequencing have allowed a deeper perspective into the molecular function of the human microbiota than has been achievable in the past. Here we aimed to examine the enzymes that cleave complex carbohydrates (CAZymes) in the human microbiome in order to determine (i) whether the CAZyme profiles of bacterial genomes are more similar within body sites or bacterial families and (ii) the sugar degradation and utilization capabilities of microbial communities inhabiting various human habitats. Upon examination of 493 bacterial references genomes from 12 human habitats, we found that sugar degradation capabilities of taxa are more similar to others in the same bacterial family than to those inhabiting the same habitat. Yet, the analysis of 520 metagenomic samples from five major body sites show that even when the community composition varies the CAZyme profiles are very similar within a body site, suggesting that the observed functional profile and microbial habitation have adapted to the local carbohydrate composition. When broad sugar utilization was compared within the five major body sites, the gastrointestinal track contained the highest potential for total sugar degradation, while dextran and peptidoglycan degradation were highest in oral and vaginal sites respectively. Our analysis suggests that the carbohydrate composition of each body site has a profound influence and probably constitutes one of the major driving forces that shapes the community composition and therefore the CAZyme profile of the local microbial communities, which in turn reflects the microbiome fitness to a body site.

Insight into trade-off between wood decay and parasitism from the genome of a fungal forest pathogen

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.