Screening for novel enzymes from metagenome and SIGEX,as a way to improve it

Metagenomics has been successfully applied to isolate novel biocatalysts from the uncultured microbiota in the environment. Two types of screening have been used to identify clones carrying desired traits from metagenomic libraries: function-based screening, and sequence-based screening. Both function- and sequence- based screening have individual advantages and disadvantages, and they have been applied successfully to discover biocatalysts from metagenome. However, both strategies are laborious and tedious because of the low frequency of screening hits. A recent paper introduced a high throughput screening strategy, termed substrate-induced gene-expression screening (SIGEX). SIGEX is designed to select the clones harboring catabolic genes induced by various substrates in concert with fluorescence activated cell sorting (FACS). This method was applied successfully to isolate aromatic hydrocarbon-induced genes from a metagenomic library. Although SIGEX has many limitations, it is expected to provide economic advantages, especially to industry.     

Substrate-induced gene-expression screening of environmental metagenome libraries for isolation of catabolic genes

Recent awareness that most microorganisms in the environment are resistant to cultivation has prompted scientists to directly clone useful genes from environmental metagenomes. Two screening methods are currently available for the metagenome approach, namely, nucleotide sequence-based screening and enzyme activity-based screening. Here we have introduced and optimized a third option for the isolation of novel catabolic operons, that is, substrate-induced gene expression screening (SIGEX). This method is based on the knowledge that catabolic-gene expression is generally induced by relevant substrates and, in many cases, controlled by regulatory elements situated proximate to catabolic genes. For SIGEX to be high throughput, we constructed an operon-trap gfp-expression vector available for shotgun cloning that allows for the selection of positive clones in liquid cultures by fluorescence-activated cell sorting. The utility of SIGEX was demonstrated by the cloning of aromatic hydrocarbon-induced genes from a groundwater metagenome library and subsequent genome-informatics analysis.     

Substrate-induced gene expression (SIGEX) screening of metagenome libraries

Substrate-induced gene-expression screening (SIGEX) has been developed for isolating novel catabolic genes from environmental metagenomes, particularly genes that are difficult to obtain using conventional gene-cloning methods. In SIGEX, restriction enzyme-digested metagenome fragments are ligated into an operon-trap vector (e.g., p18GFP), and a library is constructed in a liquid culture by transforming a cloning host (e.g., Escherichia coli). The library is subjected to a substrate-dependent gene-induction assay, and positive cells are selected by detecting activity of a co-expressed marker (e.g., GFP) encoded in the vector. High-throughput screening is possible if fluorescence-activated cell sorting (FACS) is used to select GFP-expressing cells. The abovementioned SIGEX procedure requires approximately 17 d. In this protocol, a widely applicable SIGEX scheme is presented along with typical experimental results.     

Insight into Dominant Cellulolytic Bacteria from Two Biogas Digesters and Their Glycoside Hydrolase Genes

Diverse cellulolytic bacteria are essential for maintaining high lignocellulose degradation ability in biogas digesters. However, little was known about functional genes and gene clusters of dominant cellulolytic bacteria in biogas digesters. This is the foundation to understand lignocellulose degradation mechanisms of biogas digesters and apply these gene resource for optimizing biofuel production. A combination of metagenomic and 16S rRNA gene clone library methods was used to investigate the dominant cellulolytic bacteria and their glycoside hydrolase (GH) genes in two biogas digesters. The 16S rRNA gene analysis revealed that the dominant cellulolytic bacteria were strains closely related to Clostridium straminisolvens and an uncultured cellulolytic bacterium designated BG-1. To recover GH genes from cellulolytic bacteria in general, and BG-1 in particular, a refined assembly approach developed in this study was used to assemble GH genes from metagenomic reads; 163 GH-containing contigs ≥ 1 kb in length were obtained. Six recovered GH5 genes that were expressed in E. coli demonstrated multiple lignocellulase activities and one had high mannanase activity (1255 U/mg). Eleven fosmid clones harboring the recovered GH-containing contigs were sequenced and assembled into 10 fosmid contigs. The composition of GH genes in the 163 assembled metagenomic contigs and 10 fosmid contigs indicated that diverse GHs and lignocellulose degradation mechanisms were present in the biogas digesters. In particular, a small portion of BG-1 genome information was recovered by PhyloPythiaS analysis. The lignocellulase gene clusters in BG-1 suggested that it might use a possible novel lignocellulose degradation mechanism to efficiently degrade lignocellulose. Dominant cellulolytic bacteria of biogas digester possess diverse GH genes, not only in sequences but also in their functions, which may be applied for production of biofuel in the future.     

A novel bioinformatics strategy for searching industrially useful genome resources from metagenomic sequence libraries

Although remarkable progress in metagenomic sequencing of various environmental samples has been made, large numbers of fragment sequences have been registered in the international DNA databanks, primarily without information on gene function and phylotype, and thus with limited usefulness. Industrial useful biological activity is often carried out by a set of genes, such as those constituting an operon. In this connection, metagenomic approaches have a weakness because sets of the genes are usually split up, since the sequences obtained by metagenome analyses are fragmented into 1-kb or much shorter segments. Therefore, even when a set of genes responsible for an industrially useful function is found in one metagenome library, it is usually difficult to know whether a single genome harbors the entire gene set or whether different genomes have individual genes. By modifying Self-Organizing Map (SOM), we previously developed BLSOM for oligonucleotide composition, which allowed classification (self-organization) of sequence fragments according to genomes. Because BLSOM could reassociate genomic fragments according to genomes, BLSOM may ameliorate the abovementioned weakness of metagenome analyses. Here, we have developed a strategy for clustering of metagenomic sequences according to phylotypes and genomes, by testing a gene set contributing to environment preservation.     

Accurate phylogenetic classification of variable-length DNA fragments

Metagenome studies have retrieved vast amounts of sequence data from a variety of environments leading to new discoveries and insights into the uncultured microbial world. Except for very simple communities, the encountered diversity has made fragment assembly and the subsequent analysis a challenging problem. A taxonomic characterization of metagenomic fragments is required for a deeper understanding of shotgun-sequenced microbial communities, but success has mostly been limited to sequences containing phylogenetic marker genes. Here we present PhyloPythia, a composition-based classifier that combines higher-level generic clades from a set of 340 completed genomes with sample-derived population models. Extensive analyses on synthetic and real metagenome data sets showed that PhyloPythia allows the accurate classification of most sequence fragments across all considered taxonomic ranks, even for unknown organisms. The method requires no more than 100 kb of training sequence for the creation of accurate models of sample-specific populations and can assign fragments >or=1 kb with high specificity.     

Metagenomic Characterization of Chesapeake Bay Virioplankton

Viruses are ubiquitous and abundant throughout the biosphere. In marine systems, virus-mediated processes can have significant impacts on microbial diversity and on global biogeocehmical cycling. However, viral genetic diversity remains poorly characterized. To address this shortcoming, a metagenomic library was constructed from Chesapeake Bay virioplankton. The resulting sequences constitute the largest collection of long-read double-stranded DNA (dsDNA) viral metagenome data reported to date. BLAST homology comparisons showed that Chesapeake Bay virioplankton contained a high proportion of unknown (homologous only to environmental sequences) and novel (no significant homolog) sequences. This analysis suggests that dsDNA viruses are likely one of the largest reservoirs of unknown genetic diversity in the biosphere. The taxonomic origin of BLAST homologs to viral library sequences agreed well with reported abundances of cooccurring bacterial subphyla within the estuary and indicated that cyanophages were abundant. However, the low proportion of Siphophage homologs contradicts a previous assertion that this family comprises most bacteriophage diversity. Identification and analyses of cyanobacterial homologs of the psbA gene illustrated the value of metagenomic studies of virioplankton. The phylogeny of inferred PsbA protein sequences suggested that Chesapeake Bay cyanophage strains are endemic in that environment. The ratio of psbA homologous sequences to total cyanophage sequences in the metagenome indicated that the psbA gene may be nearly universal in Chesapeake Bay cyanophage genomes. Furthermore, the low frequency of psbD homologs in the library supports the prediction that Chesapeake Bay cyanophage populations are dominated by Podoviridae.     

Phylogenetic screening of a bacterial,metagenomic library using homing endonuclease restriction and marker insertion

Metagenomics provides access to the uncultured majority of the microbial world. The approaches employed in this field have, however, had limited success in linking functional genes to the taxonomic or phylogenetic origin of the organism they belong to. Here we present an efficient strategy to recover environmental DNA fragments that contain phylogenetic marker genes from metagenomic libraries. Our method involves the cleavage of 23S ribsosmal RNA (rRNA) genes within pooled library clones by the homing endonuclease I-CeuI followed by the insertion and selection of an antibiotic resistance cassette. This approach was applied to screen a library of 6500 fosmid clones derived from the microbial community associated with the sponge Cymbastela concentrica. Several fosmid clones were recovered after the screen and detailed phylogenetic and taxonomic assignment based on the rRNA gene showed that they belong to previously unknown organisms. In addition, compositional features of these fosmid clones were used to classify and taxonomically assign a dataset of environmental shotgun sequences. Our approach represents a valuable tool for the analysis of rapidly increasing, environmental DNA sequencing information.     

Individual genome assembly from complex community short-read metagenomic datasets

Assembling individual genomes from complex community metagenomic data remains a challenging issue for environmental studies. We evaluated the quality of genome assemblies from community short read data (Illumina 100 bp pair-ended sequences) using datasets recovered from freshwater and soil microbial communities as well as in silico simulations. Our analyses revealed that the genome of a single genotype (or species) can be accurately assembled from a complex metagenome when it shows at least about 20 × coverage. At lower coverage, however, the derived assemblies contained a substantial fraction of non-target sequences (chimeras), which explains, at least in part, the higher number of hypothetical genes recovered in metagenomic relative to genomic projects. We also provide examples of how to detect intrapopulation structure in metagenomic datasets and estimate the type and frequency of errors in assembled genes and contigs from datasets of varied species complexity.     

Assembly of a Marine Viral Metagenome after Physical Fractionation

Metagenomic analyses of marine viruses generate an overview of viral genes present in a sample, but the percentage of the resulting sequence fragments that can be reassembled is low and the phenotype of the virus from which a given sequence derives is usually unknown. In this study, we employed physical fractionation to characterize the morphological and genomic traits of a subset of uncultivated viruses from a natural marine assemblage. Viruses from Kāne‘ohe Bay, Hawai‘i were fractionated by equilibrium buoyant density centrifugation in a cesium chloride (CsCl) gradient, and one fraction from the CsCl gradient was then further fractionated by strong anion-exchange chromatography. One of the fractions resulting from this two-dimensional separation appeared to be dominated by only a few virus types based on genome sizes and morphology. Sequences generated from a shotgun clone library of the viruses in this fraction were assembled into significantly more numerous contigs than have been generated with previous metagenomic investigations of whole DNA viral assemblages with comparable sequencing effort. Analysis of the longer contigs (up to 6.5 kb) assembled from our metagenome allowed us to assess gene arrangement in this subset of marine viruses. Our results demonstrate the potential for physical fractionation to facilitate sequence assembly from viral metagenomes and permit linking of morphological and genomic data for uncultivated viruses.     

Construction, Analysis, and β-Glucanase Screening of a Bacterial Artificial Chromosome Library from the Large-Bowel Microbiota of Mice

A metagenomic (community genomic) library consisting of 5,760 bacterial artificial chromosome clones was prepared in Escherichia coli DH10B from DNA extracted from the large-bowel microbiota of BALB/c mice. DNA inserts detected in 61 randomly chosen clones averaged 55 kbp (range, 8 to 150 kbp) in size. A functional screen of the library for β-glucanase activity was conducted using lichenin agar plates and Congo red solution. Three clones with β-glucanase activity were detected. The inserts of these three clones were sequenced and annotated. Open reading frames (ORF) that encoded putative proteins with identity to glucanolytic enzymes (lichenases and laminarinases) were detected by reference to databases. Other putative genes were detected, some of which might have a role in environmental sensing, nutrient acquisition, or coaggregation. The insert DNA from two clones probably originated from uncultivated bacteria because the ORF had low sequence identity with database entries, but the genes associated with the remaining clone resembled sequences reported in Bacteroides species.     

Semi-Automatic In Silico Gap Closure Enabled De Novo Assembly of Two Dehalobacter Genomes from Metagenomic Data

Typically, the assembly and closure of a complete bacterial genome requires substantial additional effort spent in a wet lab for gap resolution and genome polishing. Assembly is further confounded by subspecies polymorphism when starting from metagenome sequence data. In this paper, we describe an in silico gap-resolution strategy that can substantially improve assembly. This strategy resolves assembly gaps in scaffolds using pre-assembled contigs, followed by verification with read mapping. It is capable of resolving assembly gaps caused by repetitive elements and subspecies polymorphisms. Using this strategy, we realized the de novo assembly of the first two Dehalobacter genomes from the metagenomes of two anaerobic mixed microbial cultures capable of reductive dechlorination of chlorinated ethanes and chloroform. Only four additional PCR reactions were required even though the initial assembly with Newbler v. 2.5 produced 101 contigs within 9 scaffolds belonging to two Dehalobacter strains. By applying this strategy to the re-assembly of a recently published genome of Bacteroides, we demonstrate its potential utility for other sequencing projects, both metagenomic and genomic.     

Construction of a bacterial artificial chromosome (BAC) library for citrus and identification of BAC contigs containing resistance gene candidates

A BAC library was constructed from the genomic DNA of an intergeneric Citrus and Poncirus hybrid. The library consists of 24,576 clones with an average insert size of 115 kb, representing approximately seven haploid genome equivalents and is able to give a greater than 99% probability of isolating single-copy citrus DNA sequences from this library. High-density colony hybridization-based library screening was performed using DNA markers linked to the citrus tristeza virus (CTV) resistance gene and citrus disease resistance gene candidate (RGC) sequences. Between four and eight clones were isolated with each of the CTV resistance gene-linked markers, which agrees with the library’s predicted genome coverage. Three hundred and twenty-two clones were identified using 13 previously cloned citrus RGC sequences as probes in library screening. One to four fragments in each BAC were shown to hybridize with RGC sequences. One hundred and nine of the RGC BAC clones were fingerprinted using a sequencing gel-based procedure. From the fingerprints, 25 contigs were assembled, each having a size of 120–250 kb and consisting of 2–11 clones. These results indicate that the library is a useful resource for BAC contig construction and molecular isolation of disease resistance genes. Received: 22 May 2000 / Accepted: 25 September 2000     

Suppressive subtractive hybridization as a tool for identifying genetic diversity in an environmental metagenome: the rumen as a model

Molecular techniques previously used for genome comparisons of closely related bacterial species could prove extremely valuable for comparisons of complex microbial communities, or metagenomes. Our study aimed to determine the breadth and value of suppressive subtractive hybridization (SSH) in a pilot-scale analysis of metagenomic DNA from communities of microorganisms in the rumen. Suppressive subtractive hybridization was performed using total genomic DNA isolated from rumen fluid samples of two hay-fed steers, arbitrarily designated as tester or driver. Ninety-six subtraction DNA fragments from the tester metagenome were amplified, cloned and the DNA sequences were determined. Verification of the isolation of DNA fragments unique to the tester metagenome was accomplished through dot blot and Southern blot hybridizations. Tester-specific SSH fragments were found in 95 of 96 randomly selected clones. DNA sequences of subtraction fragments were analysed by computer assisted DNA and amino acid comparisons. Putative translations of 26 (32.1%) subtractive hybridization fragments exhibited significant similarity to Bacterial proteins, whereas 15 (18.5%) distinctive subtracted fragments had significant similarity to proteins from Archaea. The remainder of the subtractive hybridization fragments displayed no similarity to GenBank sequences. This metagenomic approach has exposed an unexpectedly large difference in Archaeal community structure between the rumen microbial populations of two steers fed identical diets and housed together. 16S rRNA dot blot hybridizations revealed similar proportions of Bacteria and Archaea in both rumen samples and suggest that the differences uncovered by SSH are the result of varying community structural composition. Our study demonstrates a novel approach to comparative analyses of environmental microbial communities through the use of SSH.     

Novel Nickel Resistance Genes from the Rhizosphere Metagenome of Plants Adapted to Acid Mine Drainage

Metal resistance determinants have traditionally been found in cultivated bacteria. To search for genes involved in nickel resistance, we analyzed the bacterial community of the rhizosphere of Erica andevalensis, an endemic heather which grows at the banks of the Tinto River, a naturally metal-enriched and extremely acidic environment in southwestern Spain. 16S rRNA gene sequence analysis of rhizosphere DNA revealed the presence of members of five phylogenetic groups of Bacteria and the two main groups of Archaea mostly associated with sites impacted by acid mine drainage (AMD). The diversity observed and the presence of heavy metals in the rhizosphere led us to construct and screen five different metagenomic libraries hosted in Escherichia coli for searching novel nickel resistance determinants. A total of 13 positive clones were detected and analyzed. Insights about their possible mechanisms of resistance were obtained from cellular nickel content and sequence similarities. Two clones encoded putative ABC transporter components, and a novel mechanism of metal efflux is suggested. In addition, a nickel hyperaccumulation mechanism is proposed for a clone encoding a serine O-acetyltransferase. Five clones encoded proteins similar to well-characterized proteins but not previously reported to be related to nickel resistance, and the remaining six clones encoded hypothetical or conserved hypothetical proteins of uncertain functions. This is the first report documenting nickel resistance genes recovered from the metagenome of an AMD environment.     

A cosmid library specific for human Chromosome regions 6p21.3 and 6q27

We have explored the potential of irradiation-fusion gene transfer (IFGT) hybrids as a source of well-defined human chromosome fragments from which probes can be derived. Extensive characterization of the IFGT hybrid 4J4 with a full panel of markers from Chromosome (Chr) 6 showed that the human DNA content derives largely from 6p21.3 and 6q27. A cosmid library has been constructed from 4J4 DNA, and 370 recombinants containing human DNA have been isolated and overlapping clones ordered into 20 contigs. Regional localization of representative clones from each contig, determined by fluorescent in situ hybridization (FISH), places 13 contigs in 6q27 and 6 in 6p21.3. Preliminary screening of cDNA libraries with selected cosmids has identified two expressed sequences. Since there are a number of medically important genes in both these regions of human Chr 6 with several disease loci linked to the HLA-A region in 6p21.3 and various tumor suppressor genes to 6q27, this library will provide a valuable resource to aid the isolation of candidate genes for these diseases. In addition, unique markers for detailed physical and genetic mapping of these regions of human Chr 6 can be easily obtained.     

Phylogenetic and Functional Analysis of Gut Microbiota of a Fungus-Growing Higher Termite: Bacteroidetes from Higher Termites Are a Rich Source of β-Glucosidase Genes

Fungus-growing termites, their symbiotic fungi, and microbiota inhibiting their intestinal tract comprise a highly efficient cellulose-hydrolyzing system; however, little is known about the role of gut microbiota in this system. Twelve fosmid clones with β-glucosidase activity were previously obtained by functionally screening a metagenomic library of a fungus-growing termite, Macrotermes annandalei. Ten contigs containing putative β-glucosidase genes (bgl1–10) were assembled by sequencing data of these fosmid clones. All these contigs were binned to Bacteroidetes, and all these β-glucosidase genes were phylogenetically closed to those from Bacteroides or Dysgonomonas. Six out of 10 β-glucosidase genes had predicted signal peptides, indicating a transmembrane capability of these enzymes to mediate cellulose hydrolysis within the gut of the termites. To confirm the activities of these β-glucosidase genes, three genes (bgl5, bgl7, and bgl9) were successfully expressed and purified. The optimal temperature and pH of these enzymes largely resembled the environment of the host’s gut. The gut microbiota composition of the fungus-growing termite was also determined by 454 pyrosequencing, showing that Bacteroidetes was the most dominant phylum. The diversity and the enzyme properties of β-glucosidases revealed in this study suggested that Bacteroidetes as the major member in fungus-growing termites contributed to cello-oligomer degradation in cellulose-hydrolyzing process and represented a rich source for β-glucosidase genes.     

Generation of polymerase chain reaction-specific probes for library screening using single degenerate primers

Degenerate oligonucleotide primers were made to peptide sequences from hydroxylamine oxidoreductase (HAO) from Nitrosomonas europaea. The primers were used singly in PCR reactions to amplify portions of the gene for HAO from genomic DNA. Southern hybridizations using fragments amplified with each primer showed that they labeled the same genomic DNA fragments. The PCR-amplified fragments were successfully used to screen a gene library for clones containing the HAO gene. The method of isolating genes by PCR with single primers has general utility.