Retrieval of glycoside hydrolase family 9 cellulase genes from environmental DNA by metagenomic gene specific multi-primer PCR

A new method, termed metagenomic gene specific multi-primer PCR (MGSM-PCR), is presented that uses multiple gene specific primers derived from an isolated gene from a constructed metagenomic library rather than degenerate primers designed based on a known enzyme family. The utility of MGSM-PCR was shown by applying it to search for homologues of the glycoside hydrolase family 9 cellulase in metagenomic DNA. The success of the multiplex PCR was verified by visualizing products on an agarose gel following gel electrophoresis. A total of 127 homologous genes were amplified with combinatorial multi-primer reactions from 34 soil DNA samples. Multiple alignments revealed extensive sequence diversity among these captured sequences with sequence identity varying from 26 to 99.7%. These results indicated that significantly diverse homologous genes were indeed readily accessible when using multiple metagenomic gene specific primers.     

Analysis of shuffled gene libraries

In vitro recombination of homologous genes (family shuffling) has been proposed as an effective search strategy for laboratory evolution of genes and proteins. Few data are available, however, on the composition of shuffled gene libraries, from which one could assess the efficiency of recombination and optimize protocols. Here, probe hybridization is used in a macroarray format to analyze chimeric DNA libraries created by DNA shuffling. Characterization of hundreds of shuffled genes encoding dioxygenases has elucidated important biases in the shuffling reaction. As expected, crossovers are favored in regions of high sequence identity. A sequence-based model of homologous recombination that captures this observed bias was formulated using the experimental results. The chimeric genes were found to show biases in the incorporation of sequences from certain parents, even before selection. Statistically different patterns of parental incorporation in genes expressing functional proteins can help to identify key sequence-function relationships.     

Discovery of a novel esterase subfamily sharing an identified arm sequence (ArmEst) by gene-specific metagenomic PCR

A gene-specific, metagenomic PCR method has led to the discovery of a novel esterase subfamily consisting of five homologous members. Sequence analysis of this esterase subfamily, named the ArmEst subfamily, revealed a unique conserved pattern with a significant variable interior sequence flanked by two symmetric and identical long arm sequences. The two homologous long arm sequences had 100 % sequence identity and symmetry at both ends between the five members of this esterase class, but only 17–58 % identity was shared for the internal sequence. The biochemical properties of two of the ArmEst esterases definitively demonstrated that they are true active esterases rather than pseudogenes. This is the first report presenting an esterase subfamily containing a unique arm sequence, indicating a rare homologous recombination occurring in the coding area of a functional gene to generate their functional diversity.     

Cloning and identification of cellulase genes from uncultured microorganisms in pulp sediments from paper mill effluent]

The metagenomic DNA of pulp sediments from paper mill effluent was extracted and purified. The 16S rDNA was amplified using the purified metagenomic DNA as template and a 16S rDNA library was prepared. Sequence analysis of 16S rDNA clones showed that diverse of uncultured bacteria inhabit in this environment, which can be classified into 4 clusters as Spirochaetes, Proteobacteria, Bacteroidetes and Firmicutes. A metagenomic library containing 10000 clones was constructed into cosmid vector, and the capacity of inserted DNA of which was 3.53 x 10(8) bp. Functional screening of the library resulted in isolation of two independent clones expressing endoglucanase activity, three independent clones expressing exoglucanase activity and two independent clones expressing beta-glucosidase activity. One clone expressing strongest enzyme activity from each activity category was chosen to be further analyzed. Three novel cellulase genes designated as umcel5L, umcel5M and umbgl3D were identified by subcloning, sequencing and expression. The umcel5L encodes an endoglucanase belonging to glycosyl hydrolase family 5, which is most related to an endoglucanase from Bradyrhizobium japonicum at 43% identity and 59% similarity. The umcel5M encodes a cellodextrinase belonging to glycosyl hydrolase family 5, which is most similar to a cellodextrinase from Fibrobacter succinogenes at 48% identity and 69% similarity. The umbgl3D encodes a putative beta-glucosidase belonging to glycosyl hydrolase family 3, which shares highest homology with a beta-glucosidase from Thermotoga maritima at 46% identity and 61% similarity. It is the first time to reveal the bacterial diversity of pulp sediments from paper mill effluent and clone novel cellulase genes from the bacteria by culture-independent method.     

Chimeric gene library construction by a simple and highly versatile method using recombination-dependent exponential amplification

A simple and efficient method for the construction of chimeric gene libraries termed RDA-PCR (recombination-dependent exponential amplification polymerase chain reaction) was developed by modifying polymerase chain reaction. A chimeric gene library is generated from homologous parental genes with additional primer-annealing sequences at their "heads" and "tails". Two primers ("skew primers") are designed to exclusively anneal to either the heads of maternal genes or the tails of paternal genes. During the RDA-PCR, short annealing/extension periods facilitate homologous recombination. The chimeric sequences can be exponentially amplified to form the chimeric gene library, whereas parental sequences without crossovers are not amplified. As a model, we constructed a chimeric gene library of yellow and green fluorescent protein (yfp and gfp, respectively). The crossover point profile of RDA-PCR clones was compared with those obtained by (modified) family shuffling. PCR restriction fragment polymorphism (PCR-RFLP) analysis of the RDA-PCR clones showed a high content of chimeric genes in the library, whereas family shuffling required the modification using skew primers for selective enrichment of chimeric sequences. PCR-RFLP analysis also indicated that the crossover points of RDA-PCR chimeras were distributed over the entire protein-coding region. Moreover, as few as 2 bp of the continual identity of nucleotides were found at the crossover points at high frequency (30% of the tested clones), suggesting that RDA-PCR resulted in a higher diversity in crossover points than family shuffling.     

DNA shuffling of subgenomic sequences of subtilisin.

DNA family shuffling of 26 protease genes was used to create a library of chimeric proteases that was screened for four distinct enzymatic properties. Multiple clones were identified that were significantly improved over any of the parental enzymes for each individual property. Family shuffling, also known as molecular breeding, efficiently created all of the combinations of parental properties, producing a great diversity of property combinations in the progeny enzymes. Thus, molecular breeding, like classical breeding, is a powerful tool for recombining existing diversity to tailor biological systems for multiple functional parameters.     

Extending the diversity of cytochrome P450 enzymes by DNA family shuffling

The cytochrome P450 enzymes involved in xenobiotic metabolism are an excellent starting point for the directed evolution of novel biocatalysts due to their wide substrate specificity. A shuffled library of three highly homologous mammalian genes (for P450 2C9, P450 2C11 and P450 2C19) was constructed by applying a modified DNA family shuffling procedure. The modifications made to the traditional DNA shuffling protocols involved non-random digestion via the use of different combinations of restriction enzymes (REs) followed by isolation of fragments under 300 bp by size-selective filtration. Shuffled cytochrome P450 mutants were co-expressed in Escherichia coli with their redox partner, NADPH-cytochrome P450 reductase (NPR). We report here how non-random fragmentation may help in chimeragenesis within the areas of low sequence similarity such as substrate recognition sites (SRSs) that are generally underrepresented in recombination using the random fragmentation process. Size-selective filtration was used to limit recovery of incompletely digested fragments and consequently minimize the chances for contamination of the shuffled library with parental forms. No parental forms could be detected in the shuffled library using restriction fragment length polymorphism (RFLP) analysis, suggesting the library was free of parental contamination. Sequencing of randomly selected mutants demonstrated a high level of chimeragenesis with on average of 8.0+/-2.2 crossovers and a low level of mutagenesis with 5.2+/-2.8 spontaneous mutations per approximately 1.5 kbp of the full-length P450 sequence. The proportion of properly folded protein as indicated by the observation of characteristic Fe(II).CO vs. Fe(II) difference spectra was 15% (4/27) of analysed mutants. Screening of the shuffled library for indole oxidation revealed four clones with similar or higher levels of indigo pigment production to those of the parental P450s and two clones with elevated P450 expression. In this paper we present a method for the effective family shuffling of cytochrome P450 enzymes, applicable to the creation of mutant libraries with expanded metabolic diversity and with a significant proportion of functional clones.     

Identification of a new subfamily of salt-tolerant esterases from a metagenomic library of tidal flat sediment

To search for novel lipolytic enzymes, a metagenomic library was constructed from the tidal flat sediment of Ganghwa Island in South Korea. By functional screening using tributyrin agar plates, 3 clones were selected from among the 80,050 clones of the fosmid library. The sequence analysis revealed that those clones contained different open reading frames, which showed 50–57% amino acid identity with putative lipolytic enzymes in the database. Based on the phylogenetic analysis, they were identified to encode novel members, which form a distinct and new subfamily in the family IV of bacterial lipolytic enzymes. The consensus sequence, GT(S)SA(G)G, encompassing the active site serine of the enzymes was different from the GDSAG motif, conserved in the other subfamily. The genes were expressed in Escherichia coli and recombinant proteins were purified as active soluble forms. The enzymes showed the highest activity toward p-nitrophenyl valerate (C5) and exhibited optimum activities at mesophilic temperature ranges and slightly alkaline pH. In particular, the enzymes displayed salt tolerance with over 50% of the maximum activity remained in the presence of 3 M NaCl (or KCl). In this study, we demonstrated that the metagenomic approach using marine tidal flat sediment as a DNA source expanded the diversity of lipolytic enzyme-encoding genes.     

Evolution of function in protein superfamilies, from a structural perspective

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.     

Construction of DNA-shuffled and incrementally truncated libraries by a mutagenic and unidirectional reassembly method: changing from a substrate specificity of phospholipase to that of lipase

A method of mutagenic and unidirectional reassembly (MURA) that can generate libraries of DNA-shuffled and randomly truncated proteins was developed. The method involved fragmenting the template gene(s) randomly by DNase I and reassembling the small fragments with a unidirectional primer by PCR. The MURA products were treated with T4 DNA polymerase and subsequently with a restriction enzyme whose site was located on the region of the MURA primer. The N-terminal-truncated and DNA-shuffled library of a Serratia sp. phospholipase A(1) prepared by this method had an essentially random variation of truncated size and also showed point mutations associated with DNA shuffling. After high-throughput screening on triglyceride-emulsified plates, several mutants exhibiting absolute lipase activity (NPL variants) were obtained. The sequence analysis and the lipase activity assay on the NPL variants revealed that N-terminal truncations at a region beginning with amino acids 61 to 71, together with amino acid substitutions, resulted in the change of substrate specificity from a phospholipase to a lipase. We therefore suggest that the MURA method, which combines incremental truncation with DNA shuffling, can contribute to expanding the searchable sequence space in directed evolution experiments.     

Screening for novel lipolytic enzymes from uncultured soil microorganisms

The construction and screening of metagenomic libraries constitute a valuable resource for obtaining novel biocatalysts. In this work, we present the construction of a metagenomic library in Escherichia coli using fosmid and microbial DNA directly isolated from forest topsoil and screened for lipolytic enzymes. The library consisted of 33,700 clones with an average DNA insert size of 35 kb. Eight unique lipolytic active clones were obtained from the metagenomic library on the basis of tributyrin hydrolysis. Subsequently, secondary libraries in a high-copy-number plasmid were generated to select lipolytic subclones and to characterize the individual genes responsible for the lipolytic activity. DNA sequence analysis of six genes revealed that the enzymes encoded by the metagenomic genes for lipolytic activity were novel with 34–48% similarity to known enzymes. They had conserved sequences similar to those in the hormone-sensitive lipase family. Based on their deduced amino acid similarity, the six genes encoding lipolytic enzymes were further divided into three subgroups, the identities among which ranged from 33% to 45%. The six predicted gene products were successfully expressed in E. coli and secreted into the culture broth. Most of the secreted enzymes showed a catalytic activity for hydrolysis of p-nitrophenyl butyrate (C₄) but not p-nitrophenyl palmitate (C₁₆).     

Identification and characterization of metagenomic fragments from tidal flat sediment

Phylogenetic surveys based on cultivation-independent methods have revealed that tidal flat sediments are environments with extensive microbial diversity. Since most of prokaryotes in nature cannot be easily cultivated under general laboratory conditions, our knowledge on prokaryotic dwellers in tidal flat sediment is mainly based on the analysis of metagenomes. Microbial community analysis based on the 16S rRNA gene and other phylogenetic markers has been widely used to provide important information on the role of microorganisms, but it is basically an indirect means, compared with direct sequencing of metagenomic DNAs. In this study, we applied a sequence-based metagenomic approach to characterize uncultivated prokaryotes from tidal flat sediment. Two large-insert genomic libraries based on fosmid were constructed from tidal flat metagenomic DNA. A survey based on end-sequencing of selected fosmid clones resulted in the identification of clones containing 274 bacterial and 16 archaeal homologs in which majority were of proteobacterial origins. Two fosmid clones containing large metagenomic DNAs were completely sequenced using the shotgun method. Both DNA inserts contained more than 20 genes encoding putative proteins which implied their ecological roles in tidal flat sediment. Phylogenetic analyses of evolutionary conserved proteins indicate that these clones are not closely related to known prokaryotes whose genome sequence is known, and genes in tidal flat may be subjected to extensive lateral gene transfer, notably between domains Bacteria and Archaea. This is the first report demonstrating that direct sequencing of metagenomic gene library is useful in underpinning the genetic makeup and functional roles of prokaryotes in tidal flat sediments.     

Construction, analysis, and beta-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 beta-glucanase activity was conducted using lichenin agar plates and Congo red solution. Three clones with beta-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.     

Enhanced crossover SCRATCHY: construction and high-throughput screening of a combinatorial library containing multiple non-homologous crossovers

SCRATCHY is a methodology for the construction of libraries of chimeras between genes that display low sequence homology. We have developed a strategy for library creation termed enhanced crossover SCRATCHY, that significantly increases the number of clones containing multiple crossovers. Complementary chimeric gene libraries generated by incremental truncation (ITCHY) of two distinct parental sequences are created, and are then divided into arbitrarily defined sections. The respective sections are amplified by skewed sets of primers (i.e. a combination of gene A specific forward primer and gene B specific reverse primer, etc.) allowing DNA fragments containing non-homologous crossover points to be amplified. The amplified chimeric sections are then subjected to a DNA shuffling process generating an enhanced crossover SCRATCHY library. We have constructed such a library using the rat theta 2 glutathione transferase (rGSTT2) and the human theta 1 glutathione transferase (hGSTT1) genes (63% DNA sequence identity). DNA sequencing analysis of unselected library members revealed a greater diversity than that obtained by canonical family shuffling or with conventional SCRATCHY. Expression and high-throughput flow cytometric screening of the chimeric GST library identified several chimeric progeny that retained rat-like parental substrate specificity.     

Thionin genes specifically expressed in barley leaves

Complementary-DNA (cDNA) clones encoding thionin were identified as one of the most frequent types of clones in a cDNA library constructed from total polyadenylated RNA from young barley leaf cells. One full-length clone codes for a precursor protein that starts with a signal peptide (28 amino acids) followed by the mature thionin (46 amino acids) and terminated by a long acidic extension (63 amino acids). The amino-acid sequence of the leaf thionin is 52% homologous to thionins from barley endosperm and in the C-terminal extension the homology decreases to 41%. In contrast, the leaf thionin is 72% homologous to viscotoxin from mistletoe leaves. Leaf thionin is coded by a multigene family with an estimated nine to eleven genes and analysis of the cDNA clones showed that at least two extremely homologous genes are expressed. Northern hybridization experiments indicate that the leaf thionin genes are not expressed in endosperm and roots. In leaves, the expression of the thionin genes is strongly repressed by light.Abbreviations cDNA complementary DNA - poly(A)RNA polyadenylated RNA     

Sequence homology of polymorphic AFLP markers in garlic (Allium sativum L.).

Linkage mapping and genetic diversity studies with DNA markers in plant species assume that comigrating bands are identical, or at least that they have homologous sequences. To test this assumption in a plant with a large genome, sequence identities of 7 polymorphic amplified fragment length polymorphism (AFLP) markers of garlic, previously used to estimate similarity in genetic diversity studies, were characterized. Among 37 diverse garlic clones, 87 bands from these 7 polymorphisms were excised, amplicons were cloned, and 2 to 6 colonies were sequenced from each band, to yield a total of 191 DNA amplicons. Of these 87 bands, 83 bands (95.4%) contained AFLP amplicons that were identical or highly homologous to the typical marker of that band; only 4 bands contained amplicons with little homology to the same-sized amplicons of other garlic clones. Of these 83 bands, 64 (73.6%) contained only highly homologous amplicons (>90% sequence identity), whereas 19 (21.8%) contained both homologous and nonhomologous amplicons, with sequence identities less than 60%. Of the 37 nonhomologous amplicons identified, 25 (67.5%) differed in length from other amplicons in the band. Sequence conservation of AFLP amplicons followed patterns similar to phylogenetic relationships among garlic clones, making them useful for developing simple PCR-based markers in genetic mapping and diversity assessment.     

Isolation of novel lipolytic genes from uncultured bacteria of pond water

Metagenomic libraries give access to gene pool of bacteria present in environmental samples avoiding the culture bias. A metagenomic library of pond water microbial assemblage in plasmid vector containing about 532 Mb of community DNA was prepared. Screening of a part of the unamplified library resulted in isolation of 11 unique lipolytic clones with an ability to hydrolyze tributyrin. DNA sequence of the lipolytic genes varied in G+C composition from 57% to 75%. Twelve lipolytic genes encoding proteins with 25-70% amino acid identity with proteins in the databases were identified. Ten of the encoded proteins belonged to seven known lipolytic protein families. One of the proteins was similar to recently identified esterase BioH. A lipolytic protein with high similarity to yet uncharacterized alpha/beta hydrolase protein family abh_upf0017 was identified from one of the clones. Conserved motif for lipolytic enzymes GXSXG, conserved aspartic and histidine residues were identified in this encoded protein.