Discovery and characterization of family 39 glycoside hydrolases from rumen anaerobic fungi with polyspecific activity on rare arabinosyl substrates

Mucosal Immunization with Spore-Based Vaccines against Mannheimia haemolytica Enhances Antigen-Specific Immunity

BACKGROUND

Mannheimia haemolytica is a bovine respiratory pathogen commonly associated with bacterial bronchopneumonia. Current vaccine strategies have shown variable efficacy in feedlot cattle, and therefore novel vaccines are needed. Bacillus subtilis spores have been investigated as a mucosal vaccine platform, due to their ability to bind and present antigens to the mucosa and act as an adjuvant. The aim of this study was to develop two spore-based mucosal vaccines targeting M. haemolytica and evaluate their immunogenicity in mice.

METHODS

Two antigen constructs composed of cholera toxin B subunit, M. haemolytica leukotoxin, and either the M. haemolytica outer membrane protein PlpE (MhCP1) or GS60 (MhCP2) were synthesized, purified and then bound to spores as vaccines. In two separate mice trials, the spore-bound vaccines (Spore-MhCP1 and Spore-MhCP2) were administered to mice through intranasal and intragastric routes, while free antigens were administered intranasally and intramuscularly. Unbound spores were also evaluated intranasally. Antigen-specific serum IgG and mucosal IgA from bronchoalveolar lavage, feces, and saliva were measured after vaccination. Mice sera from all treatment groups were assessed for their bactericidal activity against M. haemolytica.

RESULTS

In both mice experiments, intramuscular immunization induced the strongest serum IgG antibody response. However, the intranasal administration of Spore-MhCP1 and Spore-MhCP2 elicited the greatest secretory IgA-specific response against leukotoxin, PlpE, and GS60 in bronchoalveolar lavage, saliva, and feces (p < 0.05). Compared to the intranasal administration of free antigen, spore-bound antigen groups showed greater bactericidal activity against M. haemolytica (p < 0.05).

CONCLUSIONS

Since intranasally delivered Spore-MhCP1 and Spore-MhCP2 elicited both systemic and mucosal immune responses in mice, these vaccines may have potential to mitigate lung infection in cattle by restricting M. haemolytica colonization and proliferation in the respiratory tract. The efficacy of these mucosal spore-based vaccines merits further assessment against M. haemolytica in cattle.

Bifidobacterial GH146 β-L-arabinofuranosidase for the removal of β1,3-L-arabinofuranosides on plant glycans

L-Arabinofuranosides with β-linkages are present in several plant molecules, such as arabinogalactan proteins (AGPs), extensin, arabinan, and rhamnogalacturonan-II. We previously characterized a β-L-arabinofuranosidase from Bifidobacterium longum subsp. longum JCM 1217, Bll1HypBA1, which was found to belong to the glycoside hydrolase (GH) family 127. This strain encodes two GH127 genes and two GH146 genes. In the present study, we characterized a GH146 β-L-arabinofuranosidase, Bll3HypBA1 (BLLJ_1848), which was found to constitute a gene cluster with AGP-degrading enzymes. This recombinant enzyme degraded AGPs and arabinan, which contain Araf-β1,3-Araf structures. In addition, the recombinant enzyme hydrolyzed oligosaccharides containing Araf-β1,3-Araf structures but not those containing Araf-β1,2-Araf and Araf-β1,5-Araf structures. The crystal structures of Bll3HypBA1 were determined at resolutions up to 1.7 Å. The monomeric structure of Bll3HypBA1 comprised a catalytic (α/α)6 barrel and two β-sandwich domains. A hairpin structure with two β-strands was observed in Bll3HypBA1, to extend from a β-sandwich domain and partially cover the active site. The active site contains a Zn2+ ion coordinated by Cys3-Glu and exhibits structural conservation of the GH127 cysteine glycosidase Bll1HypBA1. This is the first study to report on a β1,3-specific β-L-arabinofuranosidase. KEY POINTS: • β1,3-l-Arabinofuranose residues are present in arabinogalactan proteins and arabinans as a terminal sugar. • β-l-Arabinofuranosidases are widely present in intestinal bacteria. • Bll3HypBA1 is the first enzyme characterized as a β1,3-linkage-specific β-l-arabinofuranosidase.

Biochemical unravelling of the endoxylanase activity in a bifunctional GH39 enzyme cloned and expressed from thermophilic Geobacillus sp. WSUCF1

The glycoside hydrolase family 39 (GH39) proteins are renowned for their extremophilic and multifunctional enzymatic properties, yet the molecular mechanisms underpinning these unique characteristics continue to be an active subject of research. In this study, we introduce WsuXyn, a GH39 protein with a molecular weight of 58 kDa, originating from the thermophilic Geobacillus sp. WSUCF1. Previously reported for its exceptional thermostable β-xylosidase activity, WsuXyn has recently demonstrated a significant endoxylanase activity (3752 U·mg-1) against beechwood xylan, indicating towards its bifunctional nature. Physicochemical characterization revealed that WsuXyn exhibits optimal endoxylanase activity at 70 °C and pH 7.0. Thermal stability assessments revealed that the enzyme is resilient to elevated temperatures, with a half-life of 168 h. Key kinetic parameters highlight the exceptional catalytic efficiency and strong affinity of the protein for xylan substrate. Moreover, WsuXyn-mediated hydrolysis of beechwood xylan has achieved 77 % xylan conversion, with xylose as the primary product. Structural analysis, amalgamated with docking simulations, has revealed strong binding forces between xylotetraose and the protein, with key amino acid residues, including Glu278, Tyr230, Glu160, Gly202, Cys201, Glu324, and Tyr283, playing pivotal roles in these interactions. Therefore, WsuXyn holds a strong promise for biodegradation and value-added product generation through lignocellulosic biomass conversion.

Development of a spore-based mucosal vaccine against the bovine respiratory pathogen Mannheimia haemolytica

Bovine respiratory disease (BRD) is a significant health issue in the North American feedlot industry, causing substantial financial losses due to morbidity and mortality. A lack of effective vaccines against BRD pathogens has resulted in antibiotics primarily being used for BRD prevention. The aim of this study was to develop a mucosal vaccine against the BRD pathogen, Mannheimia haemolytica, using Bacillus subtilis spores as an adjuvant. A chimeric protein (MhCP) containing a tandem repeat of neutralizing epitopes from M. haemolytica leukotoxin A (NLKT) and outer membrane protein PlpE was expressed to produce antigen for adsorption to B. subtilis spores. Adsorption was optimized by comparing varying amounts of antigen and spores, as well as different buffer pH and reaction temperatures. Using the optimal adsorption parameters, spore-bound antigen (Spore-MhCP) was prepared and administered to mice via two mucosal routes (intranasal and intragastric), while intramuscular administration of free MhCP and unvaccinated mice were used as positive and negative control treatments, respectively. Intramuscular administration of MhCP elicited the strongest serum IgG response. However, intranasal immunization of Spore-MhCP generated the best secretory IgA-specific response against both PlpE and NLKT in all samples evaluated (bronchoalveolar lavage, saliva, and feces). Since proliferation of M. haemolytica in the respiratory tract is a prerequisite to lung infection, this spore-based vaccine may offer protection in cattle by limiting colonization and subsequent infection, and Spore-MhCP warrants further evaluation in cattle as a mucosal vaccine against M. haemolytica.

Assimilation of arabinogalactan side chains with novel 3-O-β-L-arabinopyranosyl-α-L-arabinofuranosidase in Bifidobacterium pseudocatenulatum

Aim: Dietary plant fibers affect gut microbiota composition; however, the underlying microbial degradation pathways are not fully understood. We previously discovered 3-O-α-D-galactosyl-α-L-arabinofuranosidase (GAfase), a glycoside hydrolase family 39 enzyme involved in the assimilation of side chains of arabinogalactan protein (AGP), from Bifidobacterium longum subsp. longum (B. longum) JCM7052. Although GAfase homologs are not highly prevalent in the Bifidobacterium genus, several Bifidobacterium strains possess the homologs. To explore the differences in substrate specificity among the homologs, a homolog of B. longum GAfase in Bifidobacterium pseudocatenulatum MCC10289 (MCC10289_0425) was characterized. Methods: Gum arabic, larch, wheat AGP, and sugar beet arabinan were used to determine the substrate specificity of the MCC10289_0425 protein. An amino acid replacement was introduced into GAfase to identify a critical residue that governs the differentiation of substrate specificity. The growth of several Bifidobacterium strains on β-L-arabinopyranosyl disaccharide and larch AGP was examined. Results: MCC10289_0425 was identified to be an unprecedented 3-O-β-L-arabinopyranosyl-α-L-arabinofuranosidase (AAfase) with low GAfase activity. A single amino acid replacement (Asn119 to Tyr) at the catalytic site converted GAfase into AAfase. AAfase releases sugar source from AGP, thereby allowing B. pseudocatenulatum growth. Conclusion: Bifidobacteria have evolved several homologous enzymes with overlapping but distinct substrate specificities depending on the species. They have acquired different fitness abilities to respond to diverse plant polysaccharide structures.

Characterization of the Dynamic Changes of Ruminal Microbiota Colonizing Citrus Pomace Waste during Rumen Incubation for Volatile Fatty Acid Production

Rumen microorganisms are promising for efficient bioconversion of lignocellulosic wastes to biofuels and industrially relevant products. Investigating the dynamic changes of the rumen microbial community colonizing citrus pomace (CtP) will advance our understanding of the utilization of citrus processing waste by rumen fluid. Citrus pomace in nylon bags was incubated in the rumen of three ruminally cannulated Holstein cows for 1, 2, 4, 8, 12, 24, and 48 h. Results showed that total volatile fatty acids concentrations and proportions of valerate and isovalerate were increased over time during the first 12 h. Three major cellulose enzymes attached to CtP rose initially and then decreased during the 48-h incubation. Primary colonization happened during the initial hours of CtP incubation, and microbes compete to attach CtP for degrading easily digestible components and/or utilizing the waste. The 16S rRNA gene sequencing data revealed the diversity and structure of microbiota adhered to CtP were distinctly different at each time point. The increased abundance of Fibrobacterota, Rikenellaceae_RC9_gut_group, and Butyrivibrio may explain the elevated volatile fatty acids concentrations. This study highlighted key metabolically active microbial taxa colonizing citrus pomace in a 48-h in situ rumen incubation, which could have implications for promoting the biotechnological process of CtP. IMPORTANCE As a natural fermentation system, the rumen ecosystem of ruminants can efficiently degrade plant cellulose, indicating that the rumen microbiome offers an opportunity for anaerobic digestion to utilize biomass wastes containing cellulose. Knowledge of the response of the in situ microbial community to citrus pomace during anaerobic fermentation will help improve the current understanding of citrus biomass waste utilization. Our results demonstrated that a highly diverse rumen bacterial community colonized citrus pomace rapidly and continuously changed during a 48-h incubation period. These findings may provide a deep understanding of constructing, manipulating, and enriching rumen microorganisms to improve the anaerobic fermentation efficiency of citrus pomace.

Novel 3-O-α-d-Galactosyl-α-l-Arabinofuranosidase for the Assimilation of Gum Arabic Arabinogalactan Protein in Bifidobacterium longum subsp. longum

Gum arabic arabinogalactan (AG) protein (AGP) is a unique dietary fiber that is degraded and assimilated by only specific strains of Bifidobacterium longum subsp. longum Here, we identified a novel 3-O-α-d-galactosyl-α-l-arabinofuranosidase (GAfase) from B. longum JCM7052 and classified it into glycoside hydrolase family 39 (GH39). GAfase released α-d-Galp-(1→3)-l-Ara and β-l-Arap-(1→3)-l-Ara from gum arabic AGP and β-l-Arap-(1→3)-l-Ara from larch AGP, and the α-d-Galp-(1→3)-l-Ara release activity was found to be 594-fold higher than that of β-l-Arap-(1→3)-l-Ara. The GAfase gene was part of a gene cluster that included genes encoding a GH36 α-galactosidase candidate and ABC transporters for the assimilation of the released α-d-Galp-(1→3)-l-Ara in B. longum Notably, when α-d-Galp-(1→3)-l-Ara was removed from gum arabic AGP, it was assimilated by both B. longum JCM7052 and the nonassimilative B. longum JCM1217, suggesting that the removal of α-d-Galp-(1→3)-l-Ara from gum arabic AGP by GAfase permitted the cooperative action with type II AG degradative enzymes in B. longum The present study provides new insight into the mechanism of gum arabic AGP degradation in B. longumIMPORTANCE Bifidobacteria harbor numerous carbohydrate-active enzymes that degrade several dietary fibers in the gastrointestinal tract. B. longum JCM7052 is known to exhibit the ability to assimilate gum arabic AGP, but the key enzyme involved in the degradation of gum arabic AGP remains unidentified. Here, we cloned and characterized a GH39 3-O-α-d-galactosyl-α-l-arabinofuranosidase (GAfase) from B. longum JCM7052. The enzyme was responsible for the release of α-d-Galp-(1→3)-l-Ara and β-l-Arap-(1→3)-l-Ara from gum arabic AGP. The presence of a gene cluster including the GAfase gene is specifically observed in gum arabic AGP assimilative strains. However, GAfase carrier strains may affect GAfase noncarrier strains that express other type II AG degradative enzymes. These findings provide insights into the bifidogenic effect of gum arabic AGP.

Combined whole cell wall analysis and streamlined in silico carbohydrate-active enzyme discovery to improve biocatalytic conversion of agricultural crop residues

The production of biofuels as an efficient source of renewable energy has received considerable attention due to increasing energy demands and regulatory incentives to reduce greenhouse gas emissions. Second-generation biofuel feedstocks, including agricultural crop residues generated on-farm during annual harvests, are abundant, inexpensive, and sustainable. Unlike first-generation feedstocks, which are enriched in easily fermentable carbohydrates, crop residue cell walls are highly resistant to saccharification, fermentation, and valorization. Crop residues contain recalcitrant polysaccharides, including cellulose, hemicelluloses, pectins, and lignin and lignin-carbohydrate complexes. In addition, their cell walls can vary in linkage structure and monosaccharide composition between plant sources. Characterization of total cell wall structure, including high-resolution analyses of saccharide composition, linkage, and complex structures using chromatography-based methods, nuclear magnetic resonance, -omics, and antibody glycome profiling, provides critical insight into the fine chemistry of feedstock cell walls. Furthermore, improving both the catalytic potential of microbial communities that populate biodigester reactors and the efficiency of pre-treatments used in bioethanol production may improve bioconversion rates and yields. Toward this end, knowledge and characterization of carbohydrate-active enzymes (CAZymes) involved in dynamic biomass deconstruction is pivotal. Here we overview the use of common "-omics"-based methods for the study of lignocellulose-metabolizing communities and microorganisms, as well as methods for annotation and discovery of CAZymes, and accurate prediction of CAZyme function. Emerging approaches for analysis of large datasets, including metagenome-assembled genomes, are also discussed. Using complementary glycomic and meta-omic methods to characterize agricultural residues and the microbial communities that digest them provides promising streams of research to maximize value and energy extraction from crop waste streams.

Crystal structure of Dictyoglomus thermophilum β-d-xylosidase DtXyl unravels the structural determinants for efficient notoginsenoside R1 hydrolysis

Dictyoglomus thermophilum β-d-xylosidase DtXyl is attractive as a potential thermostable biocatalyst able to produce biologically active ginsenosides intermediates from β-(1,2)-D-xylosylated compounds, including Notoginsenoside-R1. DtXyl was expressed as an active N-terminal His-tagged protein, and its crystal structure was solved in presence or absence of d-xylose product. Modelling of notoginsenoside R1 in DtXyl active site led to the identification of several hydrophobic residues interacting in close contact to the substrate hydrophobic core. Unlike other residues involved in substrate binding, these residues are not conserved among GH39 xylosidase family, and their physico-chemical properties can be correlated to the efficient binding and subsequent hydrolysis of Notoginsenoside R1.

Exploring the Molecular Basis for Substrate Affinity and Structural Stability in Bacterial GH39 β-Xylosidases

The glycoside hydrolase family 39 (GH39) is a functionally expanding family with limited understanding about the molecular basis for substrate specificity and extremophilicity. In this work, we demonstrate the key role of the positive-subsite region in modulating substrate affinity and how the lack of a C-terminal extension impacts on oligomerization and structural stability of some GH39 members. The crystallographic and SAXS structures of a new GH39 member from the phytopathogen Xanthomonas citri support the importance of an extended C-terminal to promote oligomerization as a molecular strategy to enhance thermal stability. Comparative structural analysis along with site-directed mutagenesis showed that two residues located at the positive-subsite region, Lys166 and Asp167, are critical to substrate affinity and catalytic performance, by inducing local changes in the active site for substrate binding. These findings expand the molecular understanding of the mechanisms involved in substrate recognition and structural stability of the GH39 family, which might be instrumental for biological insights, rational enzyme engineering and utilization in biorefineries.

Analysis of Active Site Architecture and Reaction Product Linkage Chemistry Reveals a Conserved Cleavage Substrate for an Endo-alpha-mannanase within Diverse Yeast Mannans

Yeast α-mannan (YM) is a densely branched N-linked glycan that decorates the surface of yeast cell walls. Owing to the high degree of branching, cleavage of the backbone of YM appears to rely on the coupled action of side-chain-cleaving enzymes. Upon examining the genome sequences of bovine-adapted Bacteroides thetaiotaomicron strains, isolated for their ability to degrade YM, we have identified a tandem pair of genes inserted into an orphan pathway predicted to be involved in YM metabolism. Here, we investigated the activity of one of these enzymes, a predicted endo-mannanase from glycoside hydrolase (GH) family 76 (BtGH76-MD40). Purified recombinant BtGH76-MD40 displayed activity on structurally distinct YMs from Saccharomyces cerevisiae and Schizosaccharomyces pombe. Linkage analysis of released oligosaccharide products from S. cerevisiae and S. pombe mannan determined BtGH76-MD40 targets a specific linkage that is conserved in structurally diverse YM substrates. In addition, using two differential derivatization methods, we have shown that there is an absolute requirement for undecorated d-mannopyranose in the -1 subsite. Determination of the BtGH76-MD40 X-ray crystal structure and structural superimposition and molecular docking of a branched alpha-mannopentatose substrate supported these findings. In contrast, BtGH76-MD40 can accommodate extended side chains in the +1 and -2 subsites, highlighting that a single alpha-1,6-mannosyl residue is a prerequisite for activity, and cleavage occurs at the reducing end of the undecorated monosaccharide. Collectively these results demonstrate how acquisition of new enzymes within extant pathways contributes to the functional abilities of saccharolytic bacteria persisting in complex digestive ecosystems.

Glycoside Hydrolase Family 39 β-Xylosidases Exhibit β-1,2-Xylosidase Activity for Transformation of Notoginsenosides: A New EC Subsubclass

β-1,2-Xylosidase activity has not been recorded as an EC subsubclass. In this study, phylogenetic analysis and multiple sequence alignments revealed that characterized β-xylosidases of glycoside hydrolase family (GH) 39 were classified into the same subgroup with conserved amino acid residue positions participating in substrate recognition. Protein-ligand docking revealed that seven of these positions were probably essential to bind xylose-glucose, which is linked by a β-1,2-glycosidic bond. Amino acid residues in five of the seven positions are invariant, while those in two of the seven positions are variable with low frequency. Both the wild-type β-xylosidase rJB13GH39 and its mutants with mutation at the two positions exhibited β-1,2-xylosidase activity, as they hydrolyzed o-nitrophenyl-β-d-xylopyranoside and transformed notoginsenosides R₁ and R₂ to ginsenosides Rg₁ and Rh₁, respectively. The results suggest that all of these characterized GH 39 β-xylosidases probably show β-1,2-xylosidase activity, which should be assigned an EC number with these β-xylosidases as representatives.

Application of Transcriptomics to Compare the Carbohydrate Active Enzymes That Are Expressed by Diverse Genera of Anaerobic Fungi to Degrade Plant Cell Wall Carbohydrates

The efficiency with which the anaerobic fungi (phylum Neocallimastigomycota) degrade plant biomass is well-recognized and in recent years has received renewed interest. To further understand the biological mechanisms that are utilized by the rumen anaerobic fungi to break down lignocellulose, we have used a transcriptomic approach to examine carbohydrate digestion by Neocallimastix frontalis, Piromyces rhizinflata, Orpinomyces joyonii, and Anaeromyces mucronatus cultured on several carbon sources. The number of predicted unique transcripts ranged from 6,633 to 12,751. Pfam domains were identified in 62–70% of the fungal proteins and were linked to gene ontology terms to infer the biological function of the transcripts. Most of the predicted functions are consistent across species suggesting a similar overall strategy evolved for successful colonization of the rumen. However, the presence of differential profiles in enzyme classes suggests that there may be also be niche specialization. All fungal species were found to express an extensive array of transcripts encoding carbohydrate active enzymes (CAZymes) ranging from 8.3 to 11.3% of the transcriptome. CAZyme families involved in hemicellulose digestion were the most abundant across all four fungi. This study provides additional insight into how anaerobic fungi have evolved to become specialists at breaking down the plant cell wall in the complex and, strictly anaerobic rumen ecosystem.

SACCHARIS: an automated pipeline to streamline discovery of carbohydrate active enzyme activities within polyspecific families and de novo sequence datasets

BACKGROUND

Deposition of new genetic sequences in online databases is expanding at an unprecedented rate. As a result, sequence identification continues to outpace functional characterization of carbohydrate active enzymes (CAZymes). In this paradigm, the discovery of enzymes with novel functions is often hindered by high volumes of uncharacterized sequences particularly when the enzyme sequence belongs to a family that exhibits diverse functional specificities (i.e., polyspecificity). Therefore, to direct sequence-based discovery and characterization of new enzyme activities we have developed an automated in silico pipeline entitled: Sequence Analysis and Clustering of CarboHydrate Active enzymes for Rapid Informed prediction of Specificity (SACCHARIS). This pipeline streamlines the selection of uncharacterized sequences for discovery of new CAZyme or CBM specificity from families currently maintained on the CAZy website or within user-defined datasets.

RESULTS

SACCHARIS was used to generate a phylogenetic tree of a GH43, a CAZyme family with defined subfamily designations. This analysis confirmed that large datasets can be organized into sequence clusters of manageable sizes that possess related functions. Seeding this tree with a GH43 sequence from Bacteroides dorei DSM 17855 (BdGH43b, revealed it partitioned as a single sequence within the tree. This pattern was consistent with it possessing a unique enzyme activity for GH43 as BdGH43b is the first described α-glucanase described for this family. The capacity of SACCHARIS to extract and cluster characterized carbohydrate binding module sequences was demonstrated using family 6 CBMs (i.e., CBM6s). This CBM family displays a polyspecific ligand binding profile and contains many structurally determined members. Using SACCHARIS to identify a cluster of divergent sequences, a CBM6 sequence from a unique clade was demonstrated to bind yeast mannan, which represents the first description of an α-mannan binding CBM. Additionally, we have performed a CAZome analysis of an in-house sequenced bacterial genome and a comparative analysis of B. thetaiotaomicron VPI-5482 and B. thetaiotaomicron 7330, to demonstrate that SACCHARIS can generate "CAZome fingerprints", which differentiate between the saccharolytic potential of two related strains in silico.

CONCLUSIONS

Establishing sequence-function and sequence-structure relationships in polyspecific CAZyme families are promising approaches for streamlining enzyme discovery. SACCHARIS facilitates this process by embedding CAZyme and CBM family trees generated from biochemically to structurally characterized sequences, with protein sequences that have unknown functions. In addition, these trees can be integrated with user-defined datasets (e.g., genomics, metagenomics, and transcriptomics) to inform experimental characterization of new CAZymes or CBMs not currently curated, and for researchers to compare differential sequence patterns between entire CAZomes. In this light, SACCHARIS provides an in silico tool that can be tailored for enzyme bioprospecting in datasets of increasing complexity and for diverse applications in glycobiotechnology.