Intraspecific variation in antibiotic resistance potential within E. coli

Intraspecific genomic diversity brings the potential for an unreported and diverse reservoir of cryptic antibiotic resistance genes in pathogens, as cryptic resistance can occur without major mutations and horizontal transmission. Here, we predicted the differences in the types of antibiotics and genes that induce cryptic and latent resistance between micro-diverse Escherichia coli strains. For example, we hypothesize that known resistance genes will be the culprit of latent resistance within clinical strains. We used a modified functional metagenomics method to induce expression in eight E. coli strains. We found a total of 66 individual genes conferring phenotypic resistance to 11 out of 16 antibiotics. A total of 14 known antibiotic resistance genes comprised 21% of total identified genes, whereas the majority (52 genes) were unclassified cryptic resistance genes. Between the eight strains, 1.2% of core orthologous genes were positive (conferred resistance in at least one strain). Sixty-four percent of positive orthologous genes conferred resistance to only one strain, demonstrating high intraspecific variability of latent resistance genes. Cryptic resistance genes comprised most resistance genes among laboratory and clinical strains as well as natural, semisynthetic, and synthetic antibiotics. Known antibiotic resistance genes primarily conferred resistance to multiple antibiotics from varying origins and within multiple strains. Hence, it is uncommon for E. coli to develop cross-cryptic resistance to antibiotics from multiple origins or within multiple strains. We have uncovered prospective and previously unknown resistance genes as well as antibiotics that have the potential to trigger latent antibiotic resistance in E. coli strains from varying origins.IMPORTANCEIntraspecific genomic diversity may be a driving force in the emergence of adaptive antibiotic resistance. Adaptive antibiotic resistance enables sensitive bacterial cells to acquire temporary antibiotic resistance, creating an optimal window for the development of permanent mutational resistance. In this study, we investigate cryptic resistance, an adaptive resistance mechanism, and unveil novel (cryptic) antibiotic resistance genes that confer resistance when amplified within eight E. coli strains derived from clinical and laboratory origins. We identify the potential of cryptic resistance genes to confer cross-resistance to antibiotics from varying origins and within multiple strains. We discern antibiotic characteristics that promote latent resistance in multiple strains, considering intraspecific diversity. This study may help detect novel resistance genes and functional genes that could become responsible for cryptic resistance among diverse strains and antibiotics, thus also identifying potential novel antibiotic targets and mechanisms.

A novel family of sugar-specific phosphodiesterases that remove zwitterionic modifications of GlcNAc

The zwitterions phosphorylcholine (PC) and phosphoethanolamine (PE) are often found esterified to certain sugars in polysaccharides and glycoconjugates in a wide range of biological species. One such modification involves PC attachment to the 6-carbon of N-acetylglucosamine (GlcNAc-6-PC) in N-glycans and glycosphingolipids (GSLs) of parasitic nematodes, a modification that helps the parasite evade host immunity. Knowledge of enzymes involved in the synthesis and degradation of PC and PE modifications is limited. More detailed studies on such enzymes would contribute to a better understanding of the function of PC modifications and have potential application in the structural analysis of zwitterion-modified glycans. In this study, we used functional metagenomic screening to identify phosphodiesterases encoded in a human fecal DNA fosmid library that remove PC from GlcNAc-6-PC. A novel bacterial phosphodiesterase was identified and biochemically characterized. This enzyme (termed GlcNAc-PDase) shows remarkable substrate preference for GlcNAc-6-PC and GlcNAc-6-PE, with little or no activity on other zwitterion-modified hexoses. The identified GlcNAc-PDase protein sequence is a member of the large endonuclease/exonuclease/phosphatase superfamily where it defines a distinct subfamily of related sequences of previously unknown function, mostly from Clostridium bacteria species. Finally, we demonstrate use of GlcNAc-PDase to confirm the presence of GlcNAc-6-PC in N-glycans and GSLs of the parasitic nematode Brugia malayi in a glycoanalytical workflow.

Synthetic Biology Toolbox for Antarctic Pseudomonas sp. Strains: Toward a Psychrophilic Nonmodel Chassis for Function-Driven Metagenomics

One major limitation of function-driven metagenomics is the ability of the host to express the metagenomic DNA correctly. Differences in the transcriptional, translational, and post-translational machinery between the organism to which the DNA belongs and the host strain are all factors that influence the success of a functional screening. For this reason, the use of alternative hosts is an appropriate approach to favor the identification of enzymatic activities in function-driven metagenomics. To be implemented, appropriate tools should be designed to build the metagenomic libraries in those hosts. Moreover, discovery of new chassis and characterization of synthetic biology toolbox in nonmodel bacteria is an active field of research to expand the potential of these organisms in processes of industrial interest. Here, we assessed the suitability of two Antarctic psychrotolerant Pseudomonas strains as putative alternative hosts for function-driven metagenomics using pSEVA modular vectors as scaffold. We determined a set of synthetic biology tools suitable for these hosts and, as a proof of concept, we demonstrated their fitness for heterologous protein expression. These hosts represent a step forward for the prospection and identification of psychrophilic enzymes of biotechnological interest.

Differential response to prolonged amoxicillin treatment: long-term resilience of the microbiome versus long-lasting perturbations in the gut resistome

The collateral impact of antibiotics on the microbiome has attained increasing attention. However, the ecological consequences of long-term antibiotic exposure on the gut microbiome, including antibiotic resistance, are still limited. Here, we investigated long-term exposure effects to amoxicillin on the human gut microbiome and resistome. Fecal samples were collected from 20 patients receiving 3-months of amoxicillin or placebo treatment as part of a Norwegian multicenter clinical trial on chronic low back pain (AIM study). Samples were collected at baseline, last day of treatment, and 9 months after antibiotic cessation. The abundance and diversity of microbial and resistome composition were characterized using whole shotgun and functional metagenomic sequencing data. While the microbiome profiles of placebo subjects were stable over time, discernible changes in diversity and overall microbiome composition were observed after amoxicillin treatment. In particular, health-associated short-chain fatty acid producing species significantly decreased in proportion. However, these changes were short-lived as the microbiome showed overall recovery 9 months post-treatment. On the other hand, exposure to long-term amoxicillin was associated with an increase in total antimicrobial resistance gene load and diversity of antimicrobial resistance genes, with persistent changes even at 9 months post-treatment. Additionally, beta-lactam resistance was the most affected antibiotic class, suggesting a targeted response to amoxicillin, although changes at the gene level varied across individuals. Overall, our results suggest that the impact of prolonged amoxicillin exposure was more explicit and long-lasting in the fecal resistome than in microbiome composition. Such information is relevant for designing rational administration guidelines for antibiotic therapies.

Discovery of novel carbohydrate degrading enzymes from soda lakes through functional metagenomics

Extremophiles provide a one-of-a-kind source of enzymes with properties that allow them to endure the rigorous industrial conversion of lignocellulose biomass into fermentable sugars. However, the fact that most of these organisms fail to grow under typical culture conditions limits the accessibility to these enzymes. In this study, we employed a functional metagenomics approach to identify carbohydrate-degrading enzymes from Ethiopian soda lakes, which are extreme environments harboring a high microbial diversity. Out of 21,000 clones screened for the five carbohydrate hydrolyzing enzymes, 408 clones were found positive. Cellulase and amylase, gave high hit ratio of 1:75 and 1:280, respectively. A total of 378 genes involved in the degradation of complex carbohydrates were identified by combining high-throughput sequencing of 22 selected clones and bioinformatics analysis using a customized workflow. Around 41% of the annotated genes belonged to the Glycoside Hydrolases (GH). Multiple GHs were identified, indicating the potential to discover novel CAZymes useful for the enzymatic degradation of lignocellulose biomass from the Ethiopian soda Lakes. More than 73% of the annotated GH genes were linked to bacterial origins, with Halomonas as the most likely source. Biochemical characterization of the three enzymes from the selected clones (amylase, cellulase, and pectinase) showed that they are active in elevated temperatures, high pH, and high salt concentrations. These properties strongly indicate that the evaluated enzymes have the potential to be used for applications in various industrial processes, particularly in biorefinery for lignocellulose biomass conversion.

AMR-meta: a k-mer and metafeature approach to classify antimicrobial resistance from high-throughput short-read metagenomics data

BACKGROUND

Antimicrobial resistance (AMR) is a global health concern. High-throughput metagenomic sequencing of microbial samples enables profiling of AMR genes through comparison with curated AMR databases. However, the performance of current methods is often hampered by database incompleteness and the presence of homology/homoplasy with other non-AMR genes in sequenced samples.

RESULTS

We present AMR-meta, a database-free and alignment-free approach, based on k-mers, which combines algebraic matrix factorization into metafeatures with regularized regression. Metafeatures capture multi-level gene diversity across the main antibiotic classes. AMR-meta takes in reads from metagenomic shotgun sequencing and outputs predictions about whether those reads contribute to resistance against specific classes of antibiotics. In addition, AMR-meta uses an augmented training strategy that joins an AMR gene database with non-AMR genes (used as negative examples). We compare AMR-meta with AMRPlusPlus, DeepARG, and Meta-MARC, further testing their ensemble via a voting system. In cross-validation, AMR-meta has a median f-score of 0.7 (interquartile range, 0.2-0.9). On semi-synthetic metagenomic data-external test-on average AMR-meta yields a 1.3-fold hit rate increase over existing methods. In terms of run-time, AMR-meta is 3 times faster than DeepARG, 30 times faster than Meta-MARC, and as fast as AMRPlusPlus. Finally, we note that differences in AMR ontologies and observed variance of all tools in classification outputs call for further development on standardization of benchmarking data and protocols.

CONCLUSIONS

AMR-meta is a fast, accurate classifier that exploits non-AMR negative sets to improve sensitivity and specificity. The differences in AMR ontologies and the high variance of all tools in classification outputs call for the deployment of standard benchmarking data and protocols, to fairly compare AMR prediction tools.

Novel Viral DNA Polymerases From Metagenomes Suggest Genomic Sources of Strand-Displacing Biochemical Phenotypes

Viruses are the most abundant and diverse biological entities on the planet and constitute a significant proportion of Earth's genetic diversity. Most of this diversity is not represented by isolated viral-host systems and has only been observed through sequencing of viral metagenomes (viromes) from environmental samples. Viromes provide snapshots of viral genetic potential, and a wealth of information on viral community ecology. These data also provide opportunities for exploring the biochemistry of novel viral enzymes. The in vitro biochemical characteristics of novel viral DNA polymerases were explored, testing hypothesized differences in polymerase biochemistry according to protein sequence phylogeny. Forty-eight viral DNA Polymerase I (PolA) proteins from estuarine viromes, hot spring metagenomes, and reference viruses, encompassing a broad representation of currently known diversity, were synthesized, expressed, and purified. Novel functionality was shown in multiple PolAs. Intriguingly, some of the estuarine viral polymerases demonstrated moderate to strong innate DNA strand displacement activity at high enzyme concentration. Strand-displacing polymerases have important technological applications where isothermal reactions are desirable. Bioinformatic investigation of genes neighboring these strand displacing polymerases found associations with SNF2 helicase-associated proteins. The specific function of SNF2 family enzymes is unknown for prokaryotes and viruses. In eukaryotes, SNF2 enzymes have chromatin remodeling functions but do not separate nucleic acid strands. This suggests the strand separation function may be fulfilled by the DNA polymerase for viruses carrying SNF2 helicase-associated proteins. Biochemical data elucidated from this study expands understanding of the biology and ecological behavior of unknown viruses. Moreover, given the numerous biotechnological applications of viral DNA polymerases, novel viral polymerases discovered within viromes may be a rich source of biological material for further in vitro DNA amplification advancements.

Approaches to Unmask Functioning of the Uncultured Microbial Majority From Extreme Habitats on the Seafloor

Researchers have recognized the potential of enzymes and metabolic pathways hidden among the unseen majority of Earth’s microorganisms for decades now. Most of the microbes expected to colonize the seafloor and its subsurface are currently uncultured. Thus, their ability and contribution to element cycling remain enigmatic. Given that the seafloor covers ∼70% of our planet, this amounts to an uncalled potential of unrecognized metabolic properties and interconnections catalyzed by this microbial dark matter. Consequently, a tremendous black box awaits discovery of novel enzymes, catalytic abilities, and metabolic properties in one of the largest habitats on Earth. This mini review summarizes the current knowledge of cultivation-dependent and -independent techniques applied to seafloor habitats to unravel the role of the microbial dark matter. It highlights the great potential that combining microbiological and biogeochemical data from in situ experiments with molecular tools has for providing a holistic understanding of bio-geo-coupling in seafloor habitats and uses hydrothermal vent systems as a case example.

Identification of Glycoside Transporters From the Human Gut Microbiome

Transport is a crucial step in the metabolism of glycosides by bacteria, which is itself key for microbiota function and equilibrium. However, most transport proteins are function-unknown or only predicted, limiting our understanding of how bacteria utilize glycosides. Here, we present an activity-based screening method to identify functional glycoside transporters from microbiomes. The method is based on the co-expression in Escherichia coli of genes encoding transporters and carbohydrate-active enzymes (CAZymes) from metagenomic polysaccharide utilization loci (PULs) cloned in fosmids. To establish the proof of concept of the methodology, we used two different metagenomic libraries derived from human gut microbiota to select 18 E. coli clones whose metagenomic sequence contained at least one putative glycoside transporter and one functional CAZyme, identified by screening for various glycoside-hydrolase activities. Growth tests were performed on plant-derived glycosides, which are the target substrates of the CAZymes identified in each PUL. This led to the identification of 10 clones that are able to utilize oligosaccharides as sole carbon sources, thanks to the production of transporters from the PTS, ABC, MFS, and SusCD families. Six of the 10 hit clones contain only one transporter, providing direct experimental evidence that these transporters are functional. In the six cases where two transporters are present in the sequence of a clone, the transporters’ function can be predicted from the flanking CAZymes or from similarity with transporters characterized previously, which facilitates further functional characterization. The results expand the understanding of how glycosides are selectively metabolized by bacteria and offers a new approach to screening for glycoside-transporter specificity toward oligosaccharides with defined structures.

Discovery and Biotechnological Exploitation of Glycoside-Phosphorylases

Among carbohydrate active enzymes, glycoside phosphorylases (GPs) are valuable catalysts for white biotechnologies, due to their exquisite capacity to efficiently re-modulate oligo- and poly-saccharides, without the need for costly activated sugars as substrates. The reversibility of the phosphorolysis reaction, indeed, makes them attractive tools for glycodiversification. However, discovery of new GP functions is hindered by the difficulty in identifying them in sequence databases, and, rather, relies on extensive and tedious biochemical characterization studies. Nevertheless, recent advances in automated tools have led to major improvements in GP mining, activity predictions, and functional screening. Implementation of GPs into innovative in vitro and in cellulo bioproduction strategies has also made substantial advances. Herein, we propose to discuss the latest developments in the strategies employed to efficiently discover GPs and make the best use of their exceptional catalytic properties for glycoside bioproduction.

Identification of an Amylomaltase from the Halophilic Archaeon Haloquadratum walsbyi by Functional Metagenomics: Structural and Functional Insights

Amylomaltases are prokaryotic 4-α-glucanotransferases of the GH77 family. Thanks to the ability to modify starch, they constitute a group of enzymes of great interest for biotechnological applications. In this work we report the identification, by means of a functional metagenomics screening of the crystallization waters of the saltern of Margherita di Savoia (Italy), of an amylomaltase gene from the halophilic archaeon Haloquadratum walsbyi, and its expression in Escherichia coli cells. Sequence analysis indicated that the gene has specific insertions yet unknown in homologous genes in prokaryotes, and present only in amylomaltase genes identified in the genomes of other H. walsbyi strains. The gene is not part of any operon involved in the metabolism of maltooligosaccharides or glycogen, as it has been found in bacteria, making it impossible currently to assign a precise role to the encoded enzyme. Sequence analysis of the H. walsbyi amylomaltase and 3D modelling showed a common structure with homologous enzymes characterized in mesophilic and thermophilic bacteria. The recombinant H. walsbyi enzyme showed starch transglycosylation activity over a wide range of NaCl concentrations, with maltotriose as the best acceptor substrate compared to other maltooligosaccharides. This is the first study of an amylomaltase from a halophilic microorganism.

Secondary metabolite biosynthetic diversity in Arctic Ocean metagenomes

Polyketide synthases (PKSs) and non-ribosomal peptide synthetases (NRPSs) are mega enzymes responsible for the biosynthesis of a large fraction of natural products (NPs). Molecular markers for biosynthetic genes, such as the ketosynthase (KS) domain of PKSs, have been used to assess the diversity and distribution of biosynthetic genes in complex microbial communities. More recently, metagenomic studies have complemented and enhanced this approach by allowing the recovery of complete biosynthetic gene clusters (BGCs) from environmental DNA. In this study, the distribution and diversity of biosynthetic genes and clusters from Arctic Ocean samples (NICE-2015 expedition), was assessed using PCR-based strategies coupled with high-throughput sequencing and metagenomic analysis. In total, 149 KS domain OTU sequences were recovered, 36 % of which could not be assigned to any known BGC. In addition, 74 bacterial metagenome-assembled genomes were recovered, from which 179 BGCs were extracted. A network analysis identified potential new NP families, including non-ribosomal peptides and polyketides. Complete or near-complete BGCs were recovered, which will enable future heterologous expression efforts to uncover the respective NPs. Our study represents the first report of biosynthetic diversity assessed for Arctic Ocean metagenomes and highlights the potential of Arctic Ocean planktonic microbiomes for the discovery of novel secondary metabolites. The strategy employed in this study will enable future bioprospection, by identifying promising samples for bacterial isolation efforts, while providing also full-length BGCs for heterologous expression.

Predominant Biphenyl Dioxygenase From Legacy Polychlorinated Biphenyl (PCB)-Contaminated Soil Is a Part of Unusual Gene Cluster and Transforms Flavone and Flavanone

In this study, the diversity of bphA genes was assessed in a ¹³C-enriched metagenome upon stable isotope probing (SIP) of microbial populations in legacy PCB-contaminated soil with ¹³C-biphenyl (BP). In total, 13 bphA sequence variants (SVs) were identified in the final amplicon dataset. Of these, one SV comprised 59% of all sequences, and when it was translated into a protein sequence, it exhibited 87, 77.4, and 76.7% identity to its homologs from Pseudomonas furukawaii KF707, Cupriavidus sp. WS, and Pseudomonas alcaliphila B-367, respectively. This same BphA sequence also contained unusual amino acid residues, Alanine, Valine, and Serine in region III, which had been reported to be crucial for the substrate specificity of the corresponding biphenyl dioxygenase (BPDO), and was accordingly designated BphA_AVS. The DNA locus of 18 kbp containing the BphA_AVS-coding sequence retrieved from the metagenome was comprised of 16 ORFs and was most likely borne by Paraburkholderia sp. The BPDO corresponding to bphAE_AVS was cloned and heterologously expressed in E. coli, and its substrate specificity toward PCBs and a spectrum of flavonoids was assessed. Although depleting a rather narrow spectrum of PCB congeners, the efficient transformation of flavone and flavanone was demonstrated through dihydroxylation of the B-ring of the molecules. The homology-based functional assignment of the putative proteins encoded by the rest of ORFs in the AVS region suggests their potential involvement in the transformation of aromatic compounds, such as flavonoids. In conclusion, this study contributes to the body of information on the involvement of soil-borne BPDOs in the metabolism of flavonoid compounds, and our paper provides a more advanced context for understanding the interactions between plants, microbes and anthropogenic compounds in the soil.

Ultrahigh-Throughput Detection of Enzymatic Alcohol Dehydrogenase Activity in Microfluidic Droplets with a Direct Fluorogenic Assay

The exploration of large DNA libraries of metagenomic or synthetic origin is greatly facilitated by ultrahigh‐throughput assays that use monodisperse water‐in‐oil emulsion droplets as sequestered reaction compartments. Millions of samples can be generated and analysed in microfluidic devices at kHz speeds, requiring only micrograms of reagents. The scope of this powerful platform for the discovery of new sequence space is, however, hampered by the limited availability of assay substrates, restricting the functions and reaction types that can be investigated. Here, we broaden the scope of detectable biochemical transformations in droplet microfluidics by introducing the first fluorogenic assay for alcohol dehydrogenases (ADHs) in this format. We have synthesized substrates that release a pyranine fluorophore (8‐hydroxy‐1,3,6‐pyrenetrisulfonic acid, HPTS) when enzymatic turnover occurs. Pyranine is well retained in droplets for >6 weeks (i. e. 14‐times longer than fluorescein), avoiding product leakage and ensuring excellent assay sensitivity. Product concentrations as low as 100 nM were successfully detected, corresponding to less than one turnover per enzyme molecule on average. The potential of our substrate design was demonstrated by efficient recovery of a bona fide ADH with an >800‐fold enrichment. The repertoire of droplet screening is enlarged by this sensitive and direct fluorogenic assay to identify dehydrogenases for biocatalytic applications.

Portable and Rapid Sequencing Device with Microbial Community-Guided Culture Strategies for Precious Field and Environmental Samples

Culturing unculturable bacteria is a classic microbiology challenge; to successfully culture unculturable bacteria, microbiologists work hard to create hundreds of culture conditions. To improve the success rate and efficiency of culturing a broad spectrum of environmental microbes, it is helpful to know more about the microbial community composition. Shortening the amount of time spent sequencing, analyzing sequencing data, and predicting suitable culture conditions seems to be a critical step for improving knowledge of microbes in environmental samples and expanding culture collections. In this commentary, we introduce potential strategies using rapid and portable sequencing devices to help scientists design and plan for specific microbial culture media on their way back to the laboratory. Rapid metagenomic sequencing approaches and real-time analysis make it possible to choose microbes we are interested in and discover novel microbes in environments for cultivation.

A Titratable Cell Lysis-on-Demand System for Droplet-Compartmentalized Ultrahigh-Throughput Screening in Functional Metagenomics and Directed Evolution

Water-in-oil emulsion droplets are an attractive format for ultrahigh-throughput screening in functional metagenomics and directed evolution applications that allow libraries with more than 10⁷ members to be characterized in a day. Single library members are compartmentalized in droplets that are generated in microfluidic devices and tested for the presence of target biocatalysts. The target proteins can be produced intracellularly, for example, in bacterial hosts in-droplet cell lysis is therefore necessary to allow the enzymes to encounter the substrate to initiate an activity assay. Here, we present a titratable lysis-on-demand (LoD) system enabling the control of the cell lysis rate in Escherichia coli. We demonstrate that the rate of cell lysis can be controlled by adjusting the externally added inducer concentration. This LoD system is evaluated both at the population level (by optical density measurements) and at the single-cell level (on single-cell arrays and in alginate microbeads). Additionally, we validate the LoD system by droplet screening of a phosphotriesterase expressed from E. coli, with cell lysis triggered by inducer concentrations in the μM range. The LoD system yields sufficient release of the intracellularly produced enzymes to bring about a detectable quantity of product (measured by fluorescence in flow cytometry of double emulsions), while leaving viable cells for the downstream recovery of the genetic material.

Mosaic Ends Tagmentation (METa) Assembly for Highly Efficient Construction of Functional Metagenomic Libraries

Functional metagenomic libraries, physical bacterial libraries which allow the high-throughput capture and expression of microbiome genes, have been instrumental in the sequence-naive and cultivation-independent exploration of metagenomes. However, preparation of these libraries is often limited by their high DNA input requirement and their low cloning efficiency. Here, we describe a new method, mosaic ends tagmentation (METa) assembly, for highly efficient functional metagenomic library preparation. We applied tagmentation to metagenomic DNA from soil and gut microbiomes to prepare DNA inserts for high-throughput cloning into functional metagenomic libraries. The presence of mosaic end sequences in the resulting DNA fragments synergized with homology-based assembly cloning to result in a 300-fold increase in cloning efficiency compared to traditional blunt-cloning-based protocols. We show that compared to published libraries prepared by state-of-the-art protocols, METa assembly is on average ca. 20- to 200-fold more efficient and can prepare gigabase-sized libraries with as little as 200 ng of input DNA. We show the usefulness of METa assembly first by using a normative 5-μg mass of soil metagenomic DNA to prepare a 700-Gb library that allowed us to discover novel nourseothricin resistance genes and a potentially new mode of resistance to this antibiotic and second by using only 300 ng of goose fecal metagenomic DNA to prepare a 27-Gb library that captured numerous tetracycline and colistin resistance genes. METa assembly provides a streamlined, flexible, and efficient method for preparing functional metagenomic libraries, enabling new avenues of genetic and biochemical research into low-biomass or scarce microbiomes. IMPORTANCE Medically and industrially important genes can be recovered from microbial communities by high-throughput sequencing, but precise annotation is often limited to characterized genes and their relatives. Cloning a metagenome en masse into an expression host to produce a functional metagenomic library, directly connecting genes to functions, is a sequence-naive and cultivation-independent method to discover novel genes. The process of preparing these libraries is DNA greedy and inefficient, however. Here, we describe a library preparation method that is an order of magnitude more efficient and less DNA greedy. This method is consistently efficient across libraries prepared from cultures, a soil microbiome, and a goose fecal microbiome and allowed us to discover new antibiotic resistance genes and mechanisms. This library preparation method will potentially allow the functional metagenomic exploration of microbiomes that were previously off limits due to their rarity or low microbial biomass, such as biomedical swabs or exotic samples.

Discovery of New Carbonyl Reductases Using Functional Metagenomics and Applications in Biocatalysis

Enzyme discovery for use in the manufacture of chemicals, requiring high stereoselectivities, continues to be an important avenue of research. Here, a sequence directed metagenomics approach is described to identify short chain carbonyl reductases. PCR from a metagenomic template generated 37 enzymes, with an average 25% sequence identity, twelve of which showed interesting activities in initial screens. Six of the most productive enzymes were then tested against a panel of 21 substrates, including bulkier substrates that have been noted as challenging in biocatalytic reductions. Two enzymes were selected for further studies with the Wieland Miescher ketone. Notably, enzyme SDR-17, when co-expressed with a co-factor recycling system produced the anti-(4aR,5S) isomer in excellent isolated yields of 89% and 99% e.e. These results demonstrate the viability of a sequence directed metagenomics approach for the identification of multiple homologous sequences with low similarity, that can yield highly stereoselective enzymes with applicability in industrial biocatalysis.

Microbial Glycoside Hydrolases in the First Year of Life: An Analysis Review on Their Presence and Importance in Infant Gut

The first year of life is a crucial period during which the composition and functionality of the gut microbiota develop to stabilize and resemble that of adults. Throughout this process, the gut microbiota has been found to contribute to the maturation of the immune system, in gastrointestinal physiology, in cognitive advancement and in metabolic regulation. Breastfeeding, the “golden standard of infant nutrition,” is a cornerstone during this period, not only for its direct effect but also due to its indirect effect through the modulation of gut microbiota. Human milk is known to contain indigestible carbohydrates, termed human milk oligosaccharides (HMOs), that are utilized by intestinal microorganisms. Bacteria that degrade HMOs like Bifidobacterium longum subsp. infantis, Bifidobacterium bifidum, and Bifidobacterium breve dominate the infant gut microbiota during breastfeeding. A number of carbohydrate active enzymes have been found and identified in the infant gut, thus supporting the hypothesis that these bacteria are able to degrade HMOs. It is suggested that via resource-sharing and cross-feeding, the initial utilization of HMOs drives the interplay within the intestinal microbial communities. This is of pronounced importance since these communities promote healthy development and some of their species also persist in the adult microbiome. The emerging production and accessibility to metagenomic data make it increasingly possible to unravel the metabolic capacity of entire ecosystems. Such insights can increase understanding of how the gut microbiota in infants is assembled and makes it a possible target to support healthy growth. In this manuscript, we discuss the co-occurrence and function of carbohydrate active enzymes relevant to HMO utilization in the first year of life, based on publicly available metagenomic data. We compare the enzyme profiles of breastfed children throughout the first year of life to those of formula-fed infants.

Novel Soil-Derived Beta-Lactam, Chloramphenicol, Fosfomycin and Trimethoprim Resistance Genes Revealed by Functional Metagenomics

Antibiotic resistance genes (ARGs) in soil are considered to represent one of the largest environmental resistomes on our planet. As these genes can potentially be disseminated among microorganisms via horizontal gene transfer (HGT) and in some cases are acquired by clinical pathogens, knowledge about their diversity, mobility and encoded resistance spectra gained increasing public attention. This knowledge offers opportunities with respect to improved risk prediction and development of strategies to tackle antibiotic resistance, and might help to direct the design of novel antibiotics, before further resistances reach hospital settings or the animal sector. Here, metagenomic libraries, which comprise genes of cultivated microorganisms, but, importantly, also those carried by the uncultured microbial majority, were screened for novel ARGs from forest and grassland soils. We detected three new beta-lactam, a so far unknown chloramphenicol, a novel fosfomycin, as well as three previously undiscovered trimethoprim resistance genes. These ARGs were derived from phylogenetically diverse soil bacteria and predicted to encode antibiotic inactivation, antibiotic efflux, or alternative variants of target enzymes. Moreover, deduced gene products show a minimum identity of ~21% to reference database entries and confer high-level resistance. This highlights the vast potential of functional metagenomics for the discovery of novel ARGs from soil ecosystems.

Current Status of and Future Perspectives in Bacterial Degradation of Benzo[a]pyrene

Benzo[a]pyrene (BaP) is one the main pollutants belonging to the high-molecular-weight PAHs (HMW-PAHs) class and its degradation by microorganisms remains an important strategy for its removal from the environment. Extensive studies have been carried out on the isolation and characterisation of microorganisms that can actively degrade low-molecular-weight PAHs (LMW-PAHs), and to a certain extent, the HMW-PAH pyrene. However, so far, limited work has been carried out on BaP biodegradation. BaP consists of five fused aromatic rings, which confers this compound a high chemical stability, rendering it less amenable to biodegradation. The current review summarizes the emerging reports on BaP biodegradation. More specifically, work carried out on BaP bacterial degradation and current knowledge gaps that limit our understanding of BaP degradation are highlighted. Moreover, new avenues of research on BaP degradation are proposed, specifically in the context of the development of "omics" approaches.

The Antibacterial Agent Identified from Acidocella spp. in the Fluid of Nepenthes gracilis Against Multidrug-Resistant Klebsiella pneumoniae: A Functional Metagenomic Approach

Aims: The fluid of Nepenthes gracilis harbors diverse bacterial taxa that could serve as a gene pool for the discovery of the new genre of antimicrobial agents against multidrug-resistant Klebsiella pneumoniae. The aim of this study was to explore the presence of antibacterial genes in the fluids of N. gracilis growing in the wild. Methods: Using functional metagenomic approach, fosmid clones were isolated and screened for antibacterial activity against three strains of K. pneumoniae. A clone that exhibited the most potent antibacterial activity was sent for sequencing to identify the genes responsible for the observed activity. The secondary metabolites secreted by the selected clone was sequentially extracted using hexane, chloroform, and ethyl acetate. The chemical profiles of a clone (C6) hexane extract were determined by gas chromatography/mass spectrometry (GC-MS). Results: Fosmid clone C6 from the fluid of pitcher plant that exhibited antibacterial activity against three strains of K. pneumoniae was isolated using functional metagenome approach. A majority of the open reading frames detected from C6 were affiliated with the largely understudied Acidocella genus. Among them, the gene that encodes for coproporphyrinogen III oxidase in the heme biosynthesis pathway could be involved in the observed antibacterial activity. Based on the GC-MS analysis, the identities of the putative bioactive compounds were 2,5-di-tert-butylphenol and 1-ethyl-2-methyl cyclododecane. Conclusions: The gene that encodes for coproporphyrinogen III oxidase in the heme biosynthesis pathway as well as the secondary metabolites, namely 2,5-di-tert-butylphenol and 1-ethyl-2-methyl cyclododecane could be the potential antibacterial molecules responsible for the antibacterial activity of C6.

A New High-Throughput Screening Method to Detect Antimicrobial Volatiles from Metagenomic Clone Libraries

The ever-growing spread of resistance in medicine and agriculture highlights the need to identify new antimicrobials. Microbial volatile organic compounds (VOCs) are one of the most promising groups of chemicals to meet this need. These rarely exploited molecules exhibit antimicrobial activity and their high vapour pressure makes them ideal for application in surface sterilisation, and in particular, in biofumigation. Therefore, we adapted the previously developed Two Clamp VOCs Assay (TCVA) to a new high-throughput screening for the detection of novel antifungal VOCs from metagenomic clone libraries. As a proof of concept, we tested the new high-throughput TCVA (htTCVA) by sourcing a moss metagenomic library against Fusarium culmorum. This led to the identification of five clones that inhibited the growth of mycelium and spores in vitro by up to 8% and 30% and subsequently, to the identification of VOCs that are potentially, and in part responsible for the clones’ antifungal activity. For these VOCs, the in vitro effect of the pure compounds was as high as 100%. These results demonstrate the robustness and feasibility of the htTCVA, which provides access to completely new and unexplored biosynthetic pathways and their secondary metabolites.

Microbiome and Metagenome Analyses of a Closed Habitat during Human Occupation

This study provides the first assessment of monitoring cultivable and viable microorganisms on surfaces within a submerged, closed, analog habitat. The results of the analyses presented herein suggest that the surface material plays a role in microbial community structure, as the microbial populations differed between LDP and metal/glass surfaces. The metal/glass surfaces had less-complex community, lower bioburden, and more closely resembled the controls. These results indicated that material choice is crucial when building closed habitats, even if they are simply analogs. Finally, while a few species were associated with previously cultivated isolates from the International Space Station and MIR spacecraft, the majority of the microbial ecology of the submerged analog habitat differs greatly from that of previously studied analog habitats.

Microbial contamination during long-term confinements of space exploration presents potential risks for both crew members and spacecraft life support systems. A novel swab kit was used to sample various surfaces from a submerged, closed, analog habitat to characterize the microbial populations. Samples were collected from various locations across the habitat which were constructed from various surface materials (linoleum, dry wall, particle board, glass, and metal), and microbial populations were examined by culture, quantitative PCR (qPCR), microbiome 16S rRNA gene sequencing, and shotgun metagenomics. Propidium monoazide (PMA)-treated samples identified the viable/intact microbial population of the habitat. The cultivable microbial population ranged from below the detection limit to 10⁶ CFU/sample, and their identity was characterized using Sanger sequencing. Both 16S rRNA amplicon and shotgun sequencing were used to characterize the microbial dynamics, community profiles, and functional attributes (metabolism, virulence, and antimicrobial resistance). The 16S rRNA amplicon sequencing revealed abundance of viable (after PMA treatment) Actinobacteria (Brevibacterium, Nesternkonia, Mycobacterium, Pseudonocardia, and Corynebacterium), Firmicutes (Virgibacillus, Staphylococcus, and Oceanobacillus), and Proteobacteria (especially Acinetobacter) on linoleum, dry wall, and particle board (LDP) surfaces, while members of Firmicutes (Leuconostocaceae) and Proteobacteria (Enterobacteriaceae) were high on the glass/metal surfaces. Nonmetric multidimensional scaling determined from both 16S rRNA and metagenomic analyses revealed differential microbial species on LDP surfaces and glass/metal surfaces. The shotgun metagenomic sequencing of samples after PMA treatment showed bacterial predominance of viable Brevibacterium (53.6%), Brachybacterium (7.8%), Pseudonocardia (9.9%), Mycobacterium (3.7%), and Staphylococcus (2.1%), while fungal analyses revealed Aspergillus and Penicillium dominance.

IMPORTANCE This study provides the first assessment of monitoring cultivable and viable microorganisms on surfaces within a submerged, closed, analog habitat. The results of the analyses presented herein suggest that the surface material plays a role in microbial community structure, as the microbial populations differed between LDP and metal/glass surfaces. The metal/glass surfaces had less-complex community, lower bioburden, and more closely resembled the controls. These results indicated that material choice is crucial when building closed habitats, even if they are simply analogs. Finally, while a few species were associated with previously cultivated isolates from the International Space Station and MIR spacecraft, the majority of the microbial ecology of the submerged analog habitat differs greatly from that of previously studied analog habitats.

Personalized Mapping of Drug Metabolism by the Human Gut Microbiome

The human gut microbiome harbors hundreds of bacterial species with diverse biochemical capabilities. Dozens of drugs have been shown to be metabolized by single isolates from the gut microbiome, but the extent of this phenomenon is rarely explored in the context of microbial communities. Here, we develop a quantitative experimental framework for mapping the ability of the human gut microbiome to metabolize small molecule drugs: Microbiome-Derived Metabolism (MDM)-Screen. Included are a batch culturing system for sustained growth of subject-specific gut microbial communities, an ex vivo drug metabolism screen, and targeted and untargeted functional metagenomic screens to identify microbiome-encoded genes responsible for specific metabolic events. Our framework identifies novel drug-microbiome interactions that vary between individuals and demonstrates how the gut microbiome might be used in drug development and personalized medicine.

Novel Genes Involved in Resistance to Both Ultraviolet Radiation and Perchlorate From the Metagenomes of Hypersaline Environments

Microorganisms that thrive in hypersaline environments on the surface of our planet are exposed to the harmful effects of ultraviolet radiation. Therefore, for their protection, they have sunscreen pigments and highly efficient DNA repair and protection systems. The present study aimed to identify new genes involved in UV radiation resistance from these microorganisms, many of which cannot be cultured in the laboratory. Thus, a functional metagenomic approach was used and for this, small-insert libraries were constructed with DNA isolated from microorganisms of high-altitude Andean hypersaline lakes in Argentina (Diamante and Ojo Seco lakes, 4,589 and 3,200 m, respectively) and from the Es Trenc solar saltern in Spain. The libraries were hosted in a UV radiation-sensitive strain of Escherichia coli (recA mutant) and they were exposed to UVB. The resistant colonies were analyzed and as a result, four clones were identified with environmental DNA fragments containing five genes that conferred resistance to UV radiation in E. coli. One gene encoded a RecA-like protein, complementing the mutation in recA that makes the E. coli host strain more sensitive to UV radiation. Two other genes from the same DNA fragment encoded a TATA-box binding protein and an unknown protein, both responsible for UV resistance. Interestingly, two other genes from different and remote environments, the Ojo Seco Andean lake and the Es Trenc saltern, encoded two hypothetical proteins that can be considered homologous based on their significant amino acid similarity (49%). All of these genes also conferred resistance to 4-nitroquinoline 1-oxide (4-NQO), a compound that mimics the effect of UV radiation on DNA, and also to perchlorate, a powerful oxidant that can induce DNA damage. Furthermore, the hypothetical protein from the Es Trenc salterns was localized as discrete foci possibly associated with damaged sites in the DNA in cells treated with 4-NQO, so it could be involved in the repair of damaged DNA. In summary, novel genes involved in resistance to UV radiation, 4-NQO and perchlorate have been identified in this work and two of them encoding hypothetical proteins that could be involved in DNA damage repair activities not previously described.

A Novel, Integron-Regulated, Class C β-Lactamase

AmpC-type β-lactamases severely impair treatment of many bacterial infections, due to their broad spectrum (they hydrolyze virtually all β-lactams, except fourth-generation cephalosporins and carbapenems) and the increasing incidence of plasmid-mediated versions. The original chromosomal AmpCs are often tightly regulated, and their expression is induced in response to exposure to β-lactams. Regulation of mobile ampC expression is in many cases less controlled, giving rise to constitutively resistant strains with increased potential for development or acquisition of additional resistances. We present here the identification of two integron-encoded ampC genes, bla_IDC-1 and bla_IDC-2 (integron-derived cephalosporinase), with less than 85% amino acid sequence identity to any previously annotated AmpC. While their resistance pattern identifies them as class C β-lactamases, their low isoelectric point (pI) values make differentiation from other β-lactamases by isoelectric focusing impossible. To the best of our knowledge, this is the first evidence of an ampC gene cassette within a class 1 integron, providing a mobile context with profound potential for transfer and spread into clinics. It also allows bacteria to adapt expression levels, and thus reduce fitness costs, e.g., by cassette-reshuffling. Analyses of public metagenomes, including sewage metagenomes, show that the discovered ampCs are primarily found in Asian countries.

Dairy Products and Dairy-Processing Environments as a Reservoir of Antibiotic Resistance and Quorum-Quenching Determinants as Revealed through Functional Metagenomics

Here, the role of the dairy-processing chain as a reservoir of antimicrobial resistance (AR) determinants and a source of novel biocontrol quorum-sensing inhibitors is assessed through a functional metagenomics approach. A metagenomic library comprising ∼22,000 recombinant clones was built from DNA isolated from raw milk, raw milk cheeses, and cheese-processing environment swab samples. The high-throughput sequencing of 9,216 recombinant clones showed that lactic acid bacteria (LAB) dominated the microbial communities of raw milk cheese, while Gram-negative microorganisms of animal or soil origin dominated the microbiota of raw milk and cheese-processing environments. Although functional screening of the metagenomic library did not recover potential quorum-sensing inhibitors, in silico analysis using an in-house database built specifically for this study identified homologues to several genes encoding proteins with predicted quorum-quenching activity, among which, the QsdH hydrolase was the most abundant. In silico screening of the library identified LAB, and especially Lactococcus lactis, as a relevant reservoir of AR determinants in cheese. Functional screening of the library allowed the isolation of 13 recombinant clones showing an increased resistance toward ampicillin, which in all cases was accompanied by a reduced susceptibility to a wide range of β-lactam antibiotics. This study shows that the dairy-processing environment is a rich reservoir of AR determinants, which vary by sample source, and suggests that combining next-generation sequencing with functional metagenomics can be of use in overcoming the limitations of both approaches.IMPORTANCE The study shows the potential of functional metagenomics analyses to uncover the diversity of functions in microbial communities prevailing in dairy products and their processing environments, evidencing that lactic acid bacteria (LAB) dominate the cheese microbiota, whereas Gram-negative microorganisms of animal or soil origin dominate the microbiota of milk and cheese-processing environments. The functional and in silico screening of the library allowed the identification of LAB, and especially Lactococcus lactis, as a relevant reservoir of antimicrobial resistance (AR) determinants in cheese. Quorum-quenching (QQ) determinants were not recovered through the execution of wet-lab function-based screenings but were detected through in silico sequencing-based analyses.

Metagenomics of Atacama Lithobiontic Extremophile Life Unveils Highlights on Fungal Communities, Biogeochemical Cycles and Carbohydrate-Active Enzymes

Halites, which are typically found in various Atacama locations, are evaporitic rocks that are considered as micro-scaled salterns. Both structural and functional metagenomic analyses of halite nodules were performed. Structural analyses indicated that the halite microbiota is mainly composed of NaCl-adapted microorganisms. In addition, halites appear to harbor a limited diversity of fungal families together with a biodiverse collection of protozoa. Functional analysis indicated that the halite microbiome possesses the capacity to make an extensive contribution to carbon, nitrogen, and sulfur cycles, but possess a limited capacity to fix nitrogen. The halite metagenome also contains a vast repertory of carbohydrate active enzymes (CAZY) with glycosyl transferases being the most abundant class present, followed by glycosyl hydrolases (GH). Amylases were also present in high abundance, with GH also being identified. Thus, the halite microbiota is a potential useful source of novel enzymes that could have biotechnological applicability. This is the first metagenomic report of fungi and protozoa as endolithobionts of halite nodules, as well as the first attempt to describe the repertoire of CAZY in this community. In addition, we present a comprehensive functional metagenomic analysis of the metabolic capacities of the halite microbiota, providing evidence for the first time on the sulfur cycle in Atacama halites.

Functional metagenomics-guided discovery of potent Cas9 inhibitors in the human microbiome

CRISPR-Cas systems protect bacteria and archaea from phages and other mobile genetic elements, which use small anti-CRISPR (Acr) proteins to overcome CRISPR-Cas immunity. Because Acrs are challenging to identify, their natural diversity and impact on microbial ecosystems are underappreciated. To overcome this discovery bottleneck, we developed a high-throughput functional selection to isolate ten DNA fragments from human oral and fecal metagenomes that inhibit Streptococcus pyogenes Cas9 (SpyCas9) in Escherichia coli. The most potent Acr from this set, AcrIIA11, was recovered from a Lachnospiraceae phage. We found that AcrIIA11 inhibits SpyCas9 in bacteria and in human cells. AcrIIA11 homologs are distributed across diverse bacteria; many distantly-related homologs inhibit both SpyCas9 and a divergent Cas9 from Treponema denticola. We find that AcrIIA11 antagonizes SpyCas9 using a different mechanism than other previously characterized Type II-A Acrs. Our study highlights the power of functional selection to uncover widespread Cas9 inhibitors within diverse microbiomes.

Clinical Resistome Screening of 1,110 Escherichia coli Isolates Efficiently Recovers Diagnostically Relevant Antibiotic Resistance Biomarkers and Potential Novel Resistance Mechanisms

Multidrug-resistant pathogens represent one of the biggest global healthcare challenges. Molecular diagnostics can guide effective antibiotics therapy but relies on validated, predictive biomarkers. Here we present a novel, universally applicable workflow for rapid identification of antimicrobial resistance (AMR) biomarkers from clinical Escherichia coli isolates and quantitatively evaluate the potential to recover causal biomarkers for observed resistance phenotypes. For this, a metagenomic plasmid library from 1,110 clinical E. coli isolates was created and used for high-throughput screening to identify biomarker candidates against Tobramycin (TOB), Ciprofloxacin (CIP), and Trimethoprim-Sulfamethoxazole (TMP-SMX). Identified candidates were further validated in vitro and also evaluated in silico for their diagnostic performance based on matched genotype-phenotype data. AMR biomarkers recovered by the metagenomics screening approach mechanistically explained 77% of observed resistance phenotypes for Tobramycin, 76% for Trimethoprim-Sulfamethoxazole, and 20% Ciprofloxacin. Sensitivity for Ciprofloxacin resistance detection could be improved to 97% by complementing results with AMR biomarkers that are undiscoverable due to intrinsic limitations of the workflow. Additionally, when combined in a multiplex diagnostic in silico panel, the identified AMR biomarkers reached promising positive and negative predictive values of up to 97 and 99%, respectively. Finally, we demonstrate that the developed workflow can be used to identify potential novel resistance mechanisms.

Sex-Specific Linkages Between Taxonomic and Functional Profiles of Tick Gut Microbiomes

Ticks transmit the most diverse array of disease agents and harbor one of the most diverse microbial communities. Major progress has been made in the characterization of the taxonomic profiles of tick microbiota. However, the functional profiles of tick microbiome have been comparatively less studied. In this proof of concept we used state-of-the-art functional metagenomics analytical tools to explore previously reported datasets of bacteria found in male and female Ixodes ovatus, Ixodes persulcatus, and Amblyomma variegatum. Results showed that both taxonomic and functional profiles have differences between sexes of the same species. KEGG pathway analysis revealed that male and female of the same species had major differences in the abundance of genes involved in different metabolic pathways including vitamin B, amino acids, carbohydrates, nucleotides, and antibiotics among others. Partial reconstruction of metabolic pathways using KEGG enzymes suggests that tick microbiome form a complex metabolic network that may increase microbial community resilience and adaptability. Linkage analysis between taxonomic and functional profiles showed that among the KEGG enzymes with differential abundance in male and female ticks only 12% were present in single bacterial genera. The rest of these enzymes were found in more than two bacterial genera, and 27% of them were found in five up to ten bacterial genera. Comparison of bacterial genera contributing to the differences in the taxonomic and functional profiles of males and females revealed that while a small group of bacteria has a dual-role, most of the bacteria contribute only to functional or taxonomic differentiation between sexes. Results suggest that the different life styles of male and female ticks exert sex-specific evolutionary pressures that act independently on the phenomes (set of phenotypes) and genomes of bacteria in tick gut microbiota. We conclude that functional redundancy is a fundamental property of male and female tick microbiota and propose that functional metagenomics should be combined with taxonomic profiling of microbiota because both analyses are complementary.

Investigating Host Microbiota Relationships Through Functional Metagenomics

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

Development and Application of a High-Throughput Functional Metagenomic Screen for Glycoside Phosphorylases

Glycoside phosphorylases (GPs) catalyze the reversible phosphorolysis of glycosidic bonds, releasing sugar 1-phosphates. To identify a greater range of these under-appreciated enzymes, we have developed a high-throughput functional screening method based on molybdenum blue formation. In a proof-of-principle screen focused on cellulose-degrading GPs we interrogated ∼23,000 large insert (fosmid) clones sourced from microbial communities inhabiting two separate environments and identified seven novel GPs from carbohydrate active enzyme family GH94 and one from GH149. Characterization identified cellobiose phosphorylases, cellodextrin phosphorylases, laminaribiose phosphorylases, and a β-1,3-glucan phosphorylase. To demonstrate the versatility of the screening method, varying substrate combinations were used to identify GP activity from families GH13, GH65, GH112, and GH130 in addition to GH94 and GH149. These pilot screen and substrate versatility results provide a screening paradigm platform for recovering diverse GPs from uncultivated microbial communities acting on different substrates with considerable potential to unravel previously unknown degradative pathways within microbiomes.

Expanding the Toolbox of Broad Host-Range Transcriptional Terminators for Proteobacteria through Metagenomics

As the field of synthetic biology moves toward the utilization of novel bacterial chassis, there is a growing need for biological parts with enhanced performance in a wide number of hosts. Is not unusual that biological parts (such as promoters and terminators), initially characterized in the model bacterium Escherichia coli, do not perform well when implemented in alternative hosts, such as Pseudomonas, therefore limiting the construction of synthetic circuits in industrially relevant bacteria, for instance Pseudomonas putida. In order to address this limitation, we present here the mining of transcriptional terminators through functional metagenomics to identify novel parts with broad host-range activity. Using a GFP-based terminator trap strategy and a broad host-range plasmid, we identified 20 clones with potential terminator activity in P. putida. Further characterization allowed the identification of 4 unique sequences ranging from 58 to 181 bp long that efficiently terminate transcription in P. putida, E. coli, Burkholderia phymatum, and two Pseudomonas strains isolated from Antarctica. Therefore, this work presents a new set of biological parts useful for the engineering of synthetic circuits in Proteobacteria.

Discovery of Novel Antibiotic Resistance Determinants in Forest and Grassland Soil Metagenomes

Soil represents a significant reservoir of antibiotic resistance genes (ARGs), which can potentially spread across distinct ecosystems and be acquired by pathogens threatening human as well as animal health. Currently, information on the identity and diversity of these genes, enabling anticipation of possible future resistance development in clinical environments and the livestock sector, is lacking. In this study, we applied functional metagenomics to discover novel sulfonamide as well as tetracycline resistance genes in soils derived from forest and grassland. Screening of soil metagenomic libraries revealed a total of eight so far unknown ARGs. The recovered genes originate from phylogenetically diverse soil bacteria (e.g., Actinobacteria, Chloroflexi, or Proteobacteria) and encode proteins with a minimum identity of 46% to other antibiotic resistance determinants. In particular forest soil ecosystems have so far been neglected in studies focusing on antibiotic resistance. Here, we detected for the first time non-mobile dihydropteroate synthase (DHPS) genes conferring resistance to sulfonamides in forest soil with no history of exposure to these synthetic drugs. In total, three sulfonamide resistant DHPSs, differing in taxonomic origin, were discovered in beech or pine forest soil. This indicates that sulfonamide resistance naturally occurs in forest-resident soil bacterial communities. Besides forest soil-derived sulfonamide resistance proteins, we also identified a DHPS affiliated to Chloroflexi in grassland soil. This enzyme and the other recovered DHPSs confer reduced susceptibility toward sulfamethazine, which is widely used in food animal production. With respect to tetracycline resistance, four efflux proteins affiliated to the major facilitator superfamily (MFS) were identified. Noteworthy, one of these proteins also conferred reduced susceptibility toward lincomycin.

Nasal Resistome Development in Infants With Cystic Fibrosis in the First Year of Life

Polymicrobial infections of the respiratory tract due to antibiotic resistant bacteria are a great concern in patients with cystic fibrosis (CF). We therefore aimed at establishing a functional metagenomic method to analyze the nasal resistome in infants with CF within the first year of life. We included samples from patients before antibiotic treatment, which allowed obtaining information regarding natural status of the resistome. In total, we analyzed 130 nasal swabs from 26 infants with CF and screened for β-lactams (ampicillin, amoxicillin-clavulanic acid, and cefuroxime) and other classes of antibiotic resistances (tetracycline, chloramphenicol and trimethoprim-sulfamethoxazole). For 69 swabs (53% of total), we found at least one non-susceptible phenotype. Analyses of the inserts recovered from non-susceptible clones by nanopore MinION sequencing revealed a large reservoir of resistance genes including mobile elements within the antibiotic naïve samples. Comparing the data of the resistome with the microbiota composition showed that the bacterial phyla and operational taxonomic units (OTUs) of the microbiota rather than the antibiotic treatment were associated with the majority of non-susceptible phenotypes in the resistome. Future studies will reveal if characterization of the resistome can help in clinical decision-making in patients with CF.

Discovery and Characterization of Cas9 Inhibitors Disseminated across Seven Bacterial Phyla

CRISPR-Cas systems in bacteria and archaea provide immunity against bacteriophages and plasmids. To overcome CRISPR immunity, phages have acquired anti-CRISPR genes that reduce CRISPR-Cas activity. Using a synthetic genetic circuit, we developed a high-throughput approach to discover anti-CRISPR genes from metagenomic libraries based on their functional activity rather than sequence homology or genetic context. We identified 11 DNA fragments from soil, animal, and human metagenomes that circumvent Streptococcus pyogenes Cas9 activity in our selection strain. Further in vivo and in vitro characterization of a subset of these hits validated the activity of four anti-CRISPRs. Notably, homologs of some of these anti-CRISPRs were detected in seven different phyla, namely Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, Cyanobacteria, Spirochaetes, and Balneolaeota, and have high sequence identity suggesting recent horizontal gene transfer. Thus, anti-CRISPRs against type II-A CRISPR-Cas systems are widely distributed across bacterial phyla, suggesting a more complex ecological role than previously appreciated.

Functional Metagenomics Reveals an Overlooked Diversity and Novel Features of Soil-Derived Bacterial Phosphatases and Phytases

Phosphatases, including phytases, play a major role in cell metabolism, phosphorus cycle, biotechnology, and pathogenic processes. Nevertheless, their discovery by functional metagenomics is challenging. Here, soil metagenomic libraries were successfully screened for genes encoding phosphatase activity. In this context, we report the largest number and diversity of phosphatase genes derived from functional metagenome analysis. Two of the detected gene products carry domains which have never been associated with phosphatase activity before. One of these domains, the SNARE-associated domain DedA, harbors a so-far-overlooked motif present in numerous bacterial SNARE-associated proteins. Our analysis revealed a previously unreported phytase activity of the alkaline phosphatase and sulfatase superfamily (cl23718) and of purple acid phosphatases from nonvegetal origin. This suggests that the classical concept comprising four classes of phytases should be modified and indicates high performance of our screening method for retrieving novel types of phosphatases/phytases hidden in metagenomes of complex environments.IMPORTANCE Phosphorus (P) is a key element involved in numerous cellular processes and essential to meet global food demand. Phosphatases play a major role in cell metabolism and contribute to control the release of P from phosphorylated organic compounds, including phytate. Apart from the relationship with pathogenesis and the enormous economic relevance, phosphatases/phytases are also important for reduction of phosphorus pollution. Almost all known functional phosphatases/phytases are derived from cultured individual microorganisms. We demonstrate here for the first time the potential of functional metagenomics to exploit the phosphatase/phytase pools hidden in environmental soil samples. The recovered diversity of phosphatases/phytases comprises new types and proteins exhibiting largely unknown characteristics, demonstrating the potential of the screening method for retrieving novel target enzymes. The insights gained into the unknown diversity of genes involved in the P cycle highlight the power of function-based metagenomic screening strategies to study Earth's phosphatase pools.

iVikodak-A Platform and Standard Workflow for Inferring, Analyzing, Comparing, and Visualizing the Functional Potential of Microbial Communities

Background: The objectives of any metagenomic study typically include identification of resident microbes and their relative proportions (taxonomic analysis), profiling functional diversity (functional analysis), and comparing the identified microbes and functions with available metadata (comparative metagenomics). Given the advantage of cost-effectiveness and convenient data-size, amplicon-based sequencing has remained the technology of choice for exploring phylogenetic diversity of an environment. A recent school of thought, employing the existing genome annotation information for inferring functional capacity of an identified microbiome community, has given a promising alternative to Whole Genome Shotgun sequencing for functional analysis. Although a handful of tools are currently available for function inference, their scope, functionality and utility has essentially remained limited. Need for a comprehensive framework that expands upon the existing scope and enables a standardized workflow for function inference, analysis, and visualization, is therefore felt.

Methods: We present iVikodak, a multi-modular web-platform that hosts a logically inter-connected repertoire of functional inference and analysis tools, coupled with a comprehensive visualization interface. iVikodak is equipped with microbial co-inhabitance pattern driven published algorithms along with multiple updated databases of various curated microbe-function maps. It also features an advanced task management and result sharing system through introduction of personalized and portable dashboards.

Results: In addition to inferring functions from 16S rRNA gene data, iVikodak enables (a) an in-depth analysis of specific functions of interest (b) identification of microbes contributing to various functions (c) microbial interaction patterns through function-driven correlation networks, and (d) simultaneous functional comparison between multiple microbial communities. We have bench-marked iVikodak through multiple case studies and comparisons with existing state of art. We also introduce the concept of a public repository which provides a first of its kind community-driven framework for scientific data analytics, collaboration and sharing in this area of microbiome research.

Conclusion: Developed using modern design and task management practices, iVikodak provides a multi-modular, yet inter-operable, one-stop framework, that intends to simplify the entire approach toward inferred function analysis. It is anticipated to serve as a significant value addition to the existing space of functional metagenomics.

iVikodak web-server may be freely accessed at https://web.rniapps.net/iVikodak/.

Novel Ethanol- and 5-Hydroxymethyl Furfural-Stimulated β-Glucosidase Retrieved From a Brazilian Secondary Atlantic Forest Soil Metagenome

Beta-glucosidases are key enzymes involved in lignocellulosic biomass degradation for bioethanol production, which complete the final step during cellulose hydrolysis by converting cellobiose into glucose. Currently, industry requires enzymes with improved catalytic performance or tolerance to process-specific parameters. In this sense, metagenomics has become a powerful tool for accessing and exploring the biochemical biodiversity present in different natural environments. Here, we report the identification of a novel β-glucosidase from metagenomic DNA isolated from soil samples enriched with decaying plant matter from a Secondary Atlantic Forest region. For this, we employed a functional screening approach using an optimized and synthetic broad host-range vector for library production. The novel β-glucosidase – named Lfa2 – displays three GH3-family conserved domains and conserved catalytic amino acids D283 and E487. The purified enzyme was most active in pH 5.5 and at 50°C, and showed hydrolytic activity toward several pNP synthetic substrates containing β-glucose, β-galactose, β-xylose, β-fucose, and α-arabinopyranose, as well as toward cellobiose. Lfa2 showed considerable glucose tolerance, exhibiting an IC₅₀ of 300 mM glucose and 30% of remaining activity in 600 mM glucose. In addition, Lfa2 retained full or slightly enhanced activity in the presence of several metal ions. Further, β-glucosidase activity was increased by 1.7-fold in the presence of 10% (v/v) ethanol, a concentration that can be reached in conventional fermentation processes. Similarly, Lfa2 showed 1.7-fold enhanced activity at high concentrations of 5-hydroxymethyl furfural, one of the most important cellulase inhibitors in pretreated sugarcane bagasse hydrolysates. Moreover, the synergistic effect of Lfa2 on Bacillus subtilis GH5-CBM3 endoglucanase activity was demonstrated by the increased production of glucose (1.6-fold). Together, these results indicate that β-glucosidase Lfa2 is a promissory enzyme candidate for utilization in diverse industrial applications, such as cellulosic biomass degradation or flavor enhancement in winemaking and grape processing.

Functional metagenomics identifies an exosialidase with an inverting catalytic mechanism that defines a new glycoside hydrolase family (GH156)

Exosialidases are glycoside hydrolases that remove a single terminal sialic acid residue from oligosaccharides. They are widely distributed in biology, having been found in prokaryotes, eukaryotes, and certain viruses. Most characterized prokaryotic sialidases are from organisms that are pathogenic or commensal with mammals. However, in this study, we used functional metagenomic screening to seek microbial sialidases encoded by environmental DNA isolated from an extreme ecological niche, a thermal spring. Using recombinant expression of potential exosialidase candidates and a fluorogenic sialidase substrate, we discovered an exosialidase having no homology to known sialidases. Phylogenetic analysis indicated that this protein is a member of a small family of bacterial proteins of previously unknown function. Proton NMR revealed that this enzyme functions via an inverting catalytic mechanism, a biochemical property that is distinct from those of known exosialidases. This unique inverting exosialidase defines a new CAZy glycoside hydrolase family we have designated GH156.

Distinctive Archaeal Composition of an Artisanal Crystallizer Pond and Functional Insights Into Salt-Saturated Hypersaline Environment Adaptation

Hypersaline environments represent some of the most challenging settings for life on Earth. Extremely halophilic microorganisms have been selected to colonize and thrive in these extreme environments by virtue of a broad spectrum of adaptations to counter high salinity and osmotic stress. Although there is substantial data on microbial taxonomic diversity in these challenging ecosystems and their primary osmoadaptation mechanisms, less is known about how hypersaline environments shape the genomes of microbial inhabitants at the functional level. In this study, we analyzed the microbial communities in five ponds along the discontinuous salinity gradient from brackish to salt-saturated environments and sequenced the metagenome of the salt (halite) precipitation pond in the artisanal Cáhuil Solar Saltern system. We combined field measurements with spectrophotometric pigment analysis and flow cytometry to characterize the microbial ecology of the pond ecosystems, including primary producers and applied metagenomic sequencing for analysis of archaeal and bacterial taxonomic diversity of the salt crystallizer harvest pond. Comparative metagenomic analysis of the Cáhuil salt crystallizer pond against microbial communities from other salt-saturated aquatic environments revealed a dominance of the archaeal genus Halorubrum and showed an unexpectedly low abundance of Haloquadratum in the Cáhuil system. Functional comparison of 26 hypersaline microbial metagenomes revealed a high proportion of sequences associated with nucleotide excision repair, helicases, replication and restriction-methylation systems in all of them. Moreover, we found distinctive functional signatures between the microbial communities from salt-saturated (>30% [w/v] total salinity) compared to sub-saturated hypersaline environments mainly due to a higher representation of sequences related to replication, recombination and DNA repair in the former. The current study expands our understanding of the diversity and distribution of halophilic microbial populations inhabiting salt-saturated habitats and the functional attributes that sustain them.

A Mechanism-Based Approach to Screening Metagenomic Libraries for Discovery of Unconventional Glycosidases

Functional metagenomics has opened new opportunities for enzyme discovery. To exploit the full potential of this new tool, the design of selective screens is essential, especially when searching for rare enzymes. To identify novel glycosidases that employ cleavage strategies other than the conventional Koshland mechanisms, a suitable screen was needed. Focusing on the unsaturated glucuronidases (UGLs), it was found that use of simple aryl glycoside substrates did not allow sufficient discrimination against β-glucuronidases, which are widespread in bacteria. While conventional glycosidases cannot generally hydrolyze thioglycosides efficiently, UGLs follow a distinct mechanism that allows them to do so. Thus, fluorogenic thioglycoside substrates featuring thiol-based self-immolative linkers were synthesized and assessed as selective substrates. The generality of the approach was validated with another family of unconventional glycosidases, the GH4 enzymes. Finally, the utility of these substrates was tested by screening a small metagenomic library.

Mining Novel Constitutive Promoter Elements in Soil Metagenomic Libraries in Escherichia coli

Although functional metagenomics has been widely employed for the discovery of genes relevant to biotechnology and biomedicine, its potential for assessing the diversity of transcriptional regulatory elements of microbial communities has remained poorly explored. Here, we experimentally mined novel constitutive promoter sequences in metagenomic libraries by combining a bi-directional reporter vector, high-throughput fluorescence assays and predictive computational methods. Through the expression profiling of fluorescent clones from two independent soil sample libraries, we have analyzed the regulatory dynamics of 260 clones with candidate promoters as a set of active metagenomic promoters in the host Escherichia coli. Through an in-depth analysis of selected clones, we were able to further explore the architecture of metagenomic fragments and to report the presence of multiple promoters per fragment with a dominant promoter driving the expression profile. These approaches resulted in the identification of 33 novel active promoters from metagenomic DNA originated from very diverse phylogenetic groups. The in silico and in vivo analysis of these individual promoters allowed the generation of a constitutive promoter consensus for exogenous sequences recognizable by E. coli in metagenomic studies. The results presented here demonstrates the potential of functional metagenomics for exploring environmental bacterial communities as a source of novel regulatory genetic parts to expand the toolbox for microbial engineering.

Highly Promiscuous Oxidases Discovered in the Bovine Rumen Microbiome

The bovine rumen hosts a diverse microbiota, which is highly specialized in the degradation of lignocellulose. Ruminal bacteria, in particular, are well equipped to deconstruct plant cell wall polysaccharides. Nevertheless, their potential role in the breakdown of the lignin network has never been investigated. In this study, we used functional metagenomics to identify bacterial redox enzymes acting on polyaromatic compounds. A new methodology was developed to explore the potential of uncultured microbes to degrade lignin derivatives, namely kraft lignin and lignosulfonate. From a fosmid library covering 0.7 Gb of metagenomic DNA, three hit clones were identified, producing enzymes able to oxidize a wide variety of polyaromatic compounds without the need for the addition of copper, manganese, or mediators. These promiscuous redox enzymes could thus be of potential interest both in plant biomass refining and dye remediation. The enzymes were derived from uncultured Clostridia, and belong to complex gene clusters involving proteins of different functional types, including hemicellulases, which likely work in synergy to produce substrate degradation.

Developing a Bacteroides System for Function-Based Screening of DNA from the Human Gut Microbiome

Human gut microbiome research has been supported by advances in DNA sequencing that make it possible to obtain gigabases of sequence data from metagenomes but is limited by a lack of knowledge of gene function that leads to incomplete annotation of these data sets. There is a need for the development of methods that can provide experimental data regarding microbial gene function. Functional metagenomics is one such method, but functional screens are often carried out using hosts that may not be able to express the bulk of the environmental DNA being screened. We expand the range of current screening hosts and demonstrate that human gut-derived metagenomic libraries can be introduced into the gut microbe Bacteroides thetaiotaomicron to identify genes based on activity screening. Our results support the continuing development of genetically tractable systems to obtain information about gene function.

Functional metagenomics is a powerful method that allows the isolation of genes whose role may not have been predicted from DNA sequence. In this approach, first, environmental DNA is cloned to generate metagenomic libraries that are maintained in Escherichia coli, and second, the cloned DNA is screened for activities of interest. Typically, functional screens are carried out using E. coli as a surrogate host, although there likely exist barriers to gene expression, such as lack of recognition of native promoters. Here, we describe efforts to develop Bacteroides thetaiotaomicron as a surrogate host for screening metagenomic DNA from the human gut. We construct a B. thetaiotaomicron-compatible fosmid cloning vector, generate a fosmid clone library using DNA from the human gut, and show successful functional complementation of a B. thetaiotaomicron glycan utilization mutant. Though we were unable to retrieve the physical fosmid after complementation, we used genome sequencing to identify the complementing genes derived from the human gut microbiome. Our results demonstrate that the use of B. thetaiotaomicron to express metagenomic DNA is promising, but they also exemplify the challenges that can be encountered in the development of new surrogate hosts for functional screening.

IMPORTANCE Human gut microbiome research has been supported by advances in DNA sequencing that make it possible to obtain gigabases of sequence data from metagenomes but is limited by a lack of knowledge of gene function that leads to incomplete annotation of these data sets. There is a need for the development of methods that can provide experimental data regarding microbial gene function. Functional metagenomics is one such method, but functional screens are often carried out using hosts that may not be able to express the bulk of the environmental DNA being screened. We expand the range of current screening hosts and demonstrate that human gut-derived metagenomic libraries can be introduced into the gut microbe Bacteroides thetaiotaomicron to identify genes based on activity screening. Our results support the continuing development of genetically tractable systems to obtain information about gene function.

New Bacterial Phytase through Metagenomic Prospection

Alkaline phytases from uncultured microorganisms, which hydrolyze phytate to less phosphorylated myo-inositols and inorganic phosphate, have great potential as additives in agricultural industry. The development of metagenomics has stemmed from the ineluctable evidence that as-yet-uncultured microorganisms represent the vast majority of organisms in most environments on earth. In this study, a gene encoding a phytase was cloned from red rice crop residues and castor bean cake using a metagenomics strategy. The amino acid identity between this gene and its closest published counterparts is lower than 60%. The phytase was named PhyRC001 and was biochemically characterized. This recombinant protein showed activity on sodium phytate, indicating that PhyRC001 is a hydrolase enzyme. The enzymatic activity was optimal at a pH of 7.0 and at a temperature of 35 °C. β-propeller phytases possess great potential as feed additives because they are the only type of phytase with high activity at neutral pH. Therefore, to explore and exploit the underlying mechanism for β-propeller phytase functions could be of great benefit to biotechnology.

Acid-resistant genes of oral plaque microbiome from the functional metagenomics

Acid resistance is one of key properties assisting the survival of cariogenic bacteria in a dental caries environment, but only a few genes conferring acid resistance have been identified to data. Functional metagenomics provides a systematic method for investigating commensal DNA to identify genes that encode target functions. Here, the host strain Escherichia coli DH10B and a constructed bidirectional transcription vector pSKII⁺-lacZ contributed to the construction of a metagenomic library, and 46.6 Mb of metagenomic DNA was cloned from carious supragingival plaque of 8children along with screening for lethal functionality. The screen identified 2 positive clones that exhibited a similar aciduric phenotype to that of the positive controls. Bioinformatic analysis revealed that these two genes encoded an ATP/GTP-binding protein and a malate dehydrogenase. Moreover, we also performed functional screening of Streptococcus mutans, since it is one of the predominant cariogenic strains but was not identified in our initial screening. Five positive clones were retrieved. In conclusion, our improved functional metagenomics screening method helped in the identification of important acid resistance genes, thereby providing new insights into the mechanism underlying caries formation as well as in the prevention and treatment of early childhood caries (ECC).