Quantifying fluorescent glycan uptake to elucidate strain-level variability in foraging behaviors of rumen bacteria.

Gut microbiomes, such as the microbial community that colonizes the rumen, have vast catabolic potential and play a vital role in host health and nutrition. By expanding our understanding of metabolic pathways in these ecosystems, we will garner foundational information for manipulating microbiome structure and function to influence host physiology. Currently, our knowledge of metabolic pathways relies heavily on inferences derived from metagenomics or culturing bacteria in vitro. However, novel approaches targeting specific cell physiologies can illuminate the functional potential encoded within microbial (meta)genomes to provide accurate assessments of metabolic abilities. Using fluorescently labeled polysaccharides, we visualized carbohydrate metabolism performed by single bacterial cells in a complex rumen sample, enabling a rapid assessment of their metabolic phenotype. Specifically, we identified bovine-adapted strains of Bacteroides thetaiotaomicron that metabolized yeast mannan in the rumen microbiome ex vivo and discerned the mechanistic differences between two distinct carbohydrate foraging behaviors, referred to as "medium grower" and "high grower." Using comparative whole-genome sequencing, RNA-seq, and carbohydrate-active enzyme fingerprinting, we could elucidate the strain-level variability in carbohydrate utilization systems of the two foraging behaviors to help predict individual strategies of nutrient acquisition. Here, we present a multi-faceted study using complimentary next-generation physiology and "omics" approaches to characterize microbial adaptation to a prebiotic in the rumen ecosystem. Video abstract.

Molecular basis of an agarose metabolic pathway acquired by a human intestinal symbiont.

In red algae, the most abundant principal cell wall polysaccharides are mixed galactan agars, of which agarose is a common component. While bioconversion of agarose is predominantly catalyzed by bacteria that live in the oceans, agarases have been discovered in microorganisms that inhabit diverse terrestrial ecosystems, including human intestines. Here we comprehensively define the structure-function relationship of the agarolytic pathway from the human intestinal bacterium Bacteroides uniformis (Bu) NP1. Using recombinant agarases from Bu NP1 to completely depolymerize agarose, we demonstrate that a non-agarolytic Bu strain can grow on GAL released from agarose. This relationship underscores that rare nutrient utilization by intestinal bacteria is facilitated by the acquisition of highly specific enzymes that unlock inaccessible carbohydrate resources contained within unusual polysaccharides. Intriguingly, the agarolytic pathway is differentially distributed throughout geographically distinct human microbiomes, reflecting a complex historical context for agarose consumption by human beings.

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.

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

Enzyme activities that improve digestion of recalcitrant plant cell wall polysaccharides may offer solutions for sustainable industries. To this end, anaerobic fungi in the rumen have been identified as a promising source of novel carbohydrate active enzymes (CAZymes) that modify plant cell wall polysaccharides and other complex glycans. Many CAZymes share insufficient sequence identity to characterized proteins from other microbial ecosystems to infer their function; thus presenting challenges to their identification. In this study, four rumen fungal genes (nf2152, nf2215, nf2523, and pr2455) were identified that encode family 39 glycoside hydrolases (GH39s), and have conserved structural features with GH51s. Two recombinant proteins, NF2152 and NF2523, were characterized using a variety of biochemical and structural techniques, and were determined to have distinct catalytic activities. NF2152 releases a single product, ß1,2-arabinobiose (Ara2) from sugar beet arabinan (SBA), and ß1,2-Ara2 and α-1,2-galactoarabinose (Gal-Ara) from rye arabinoxylan (RAX). NF2523 exclusively releases α-1,2-Gal-Ara from RAX, which represents the first description of a galacto-(α-1,2)-arabinosidase. Both ß-1,2-Ara2 and α-1,2-Gal-Ara are disaccharides not previously described within SBA and RAX. In this regard, the enzymes studied here may represent valuable new biocatalytic tools for investigating the structures of rare arabinosyl-containing glycans, and potentially for facilitating their modification in industrial applications.

Comparative variation within the genome of Campylobacter jejuni NCTC 11168 in human and murine hosts.

Campylobacteriosis incited by C. jejuni is a significant enteric disease of human beings. A person working with two reference strains of C. jejuni National Collection of Type Cultures (NCTC) 11168 developed symptoms of severe enteritis including bloody diarrhea. The worker was determined to be infected by C. jejuni. In excess of 50 isolates were recovered from the worker's stool. All of the recovered isolates and the two reference strains were indistinguishable from each other based on comparative genomic fingerprint subtyping. Whole genome sequence analysis indicated that the worker was infected with a C. jejuni NCTC 11168 obtained from the American Type Culture Collection; this strain (NCTC 11168-GSv) is the genome sequence reference. After passage through the human host, major genetic changes including indel mutations within twelve contingency loci conferring phase variations were detected in the genome of C. jejuni. Specific and robust single nucleotide polymorphism (SNP) changes in the human host were also observed in two loci (Cj0144c, Cj1564). In mice inoculated with an isolate of C. jejuni NCTC 11168-GSv from the infected person, the isolate underwent further genetic variation. At nine loci, mutations specific to inoculated mice including five SNP changes were observed. The two predominant SNPs observed in the human host reverted in mice. Genetic variations occurring in the genome of C. jejuni in mice corresponded to increased densities of C. jejuni cells associated with cecal mucosa. In conclusion, C. jejuni NCTC 11168-GSv was found to be highly virulent in a human being inciting severe enteritis. Host-specific mutations in the person with enteritis occurred/were selected for in the genome of C. jejuni, and many were not maintained in mice. Information obtained in the current study provides new information on host-specific genetic adaptation by C. jejuni.

An ancestral member of the polysaccharide lyase family 2 displays endolytic activity and magnesium dependence.

Polysaccharide lyases (PLs) are enzymes that cleave glycosidic linkages in hexuronate polysaccharides, such as homogalacturonan (HG), using a ß-elimination mechanism. Traditionally, PL activities on HG have been associated with catalytic calcium cofactors, unusually high pH optima, and arginine Brønstead bases. Recently, however, PL families that harness transition metal cofactors, utilize lysine and histidine Brønstead bases, and display more neutral pH optima have been described. One such family is PL2, which has members found primarily in phytopathogenic (e.g., Dickeya spp. and Pectobacterium spp.) or enteropathogenic (e.g., Yersinia spp.) bacterial species. PL2 is divided into two major subfamilies that are correlated with either an endolytic or exolytic activity. This study has focused on the activity of a PL2 member, which is not classified within either subfamily and helps to illuminate the origin of enzyme activities within the family. In addition, the role of Mg(2+) as a preferential catalytic metal for an intracellular PL2 (PaePL2) is described. The implications for the relationship between catalytic metal selectivity and the cellular location of pectate lyase-mediated catalysis are discussed.

Molecular methods to measure intestinal bacteria: a review.

The intestine is an exceptionally rich ecosystem encompassing a complex interaction among microorganisms, influenced by host factors, ingested food, and liquid. Characterizing the intestinal microbiota is currently an active area of research. Various molecular-based methods are available to characterize the intestinal microbiota, but all methods possess relative strengths, as well as salient weaknesses. It is important that researchers are cognizant of the limitations of these methods, and that they take the appropriate steps to mitigate weaknesses. Here, we discuss methodologies used to monitor intestinal bacteria including: (i) traditional clone libraries; (ii) direct sequencing using next-generation parallel sequencing technology; (iii) denaturing gradient gel electrophoresis and temperature gradient gel electrophoresis; (iv) terminal restriction fragment length polymorphism analysis; (v) fluorescent in situ hybridization; and (vi) quantitative PCR. In addition, we also discuss experimental design, sample collection and storage, DNA extraction, gene targets, PCR bias, and methods to reduce PCR bias.

spyder, a new method for in silico design and assessment of 16S rRNA gene primers for molecular microbial ecology.

Molecular microbial ecology studies are heavily reliant on 'Universal' 16S rRNA gene primers for elucidating microbial community structure and composition, and yet primer design and optimization is often overlooked. Primers that exhibit minor biases due to primer-template mismatches can substantially alter the pool of amplicons from a community DNA sample, resulting in inaccurate conclusions. As a result, it is important that primers are critically evaluated against the most comprehensive data sets available before commencing molecular microbial community studies. We present a user-friendly, multi-platform (e.g. Windows, Linux, Mac) method named spyder for the in silico design and assessment of 16S rRNA gene primers. The method utilizes the Ribosomal Database Project's Probe Match feature coupled with a compact program (available at http://people.uleth.ca/~selibl/Spyder/Spyder.html) that aligns and identifies mismatches between primers and templates. To demonstrate the value of spyder, we assessed commonly used 'Universal' and phyla-specific primers and identified primer modifications that improved the coverage of target organisms by 5-42% as well as removed excessive degeneracies.

Sequence variant discovery in DNA repair genes from radiosensitive and radiotolerant prostate brachytherapy patients.

PURPOSE: The presence of intrinsic radiosensitivity within prostate cancer patients may be an important factor contributing to development of radiation toxicity. We investigated whether variants in genes responsible for detecting and repairing DNA damage independently contribute to toxicity following prostate brachytherapy. EXPERIMENTAL DESIGN: Genomic DNA was extracted from blood samples of 41 prostate brachytherapy patients, 21 with high and 20 with low late toxicity scores. For each patient, 242 PCR amplicons were generated containing 173 exons of eight candidate genes: ATM, BRCA1, ERCC2, H2AFX, LIG4, MDC1, MRE11A, and RAD50. These amplicons were sequenced and all sequence variants were subjected to statistical analysis to identify those associated with late radiation toxicity. RESULTS: Across 41 patients, 239 sites differed from the human genome reference sequence; 170 of these corresponded to known polymorphisms. Sixty variants, 14 of them novel, affected protein coding regions and 43 of these were missense mutations. In our patient population, the high toxicity group was enriched for individuals with at least one LIG4 coding variant (P = 0.028). One synonymous variant in MDC1, rs28986317, was associated with increased radiosensitivity (P = 0.048). A missense variant in ATM, rs1800057, associated with increased prostate cancer risk, was found exclusively in two high toxicity patients but did not reach statistical significance for association with radiosensitivity (P = 0.488). CONCLUSIONS: Our data revealed new germ-line sequence variants, indicating that existing sequence databases do not fully represent the full extent of sequence variation. Variants in three DNA repair genes were linked to increased radiosensitivity but require validation in larger populations.