Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2008, p. 723-730, Vol. 74, No. 3
0099-2240/08/$08.00+0 doi:10.1128/AEM.01911-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

He Kyoung Lim,3,
Jin-Cheol Kim,3
Gyung Ja Choi,3
Eun Jin Park,1
Myung Hwan Lee,4
Young Ryun Chung,2 and
Seon-Woo Lee1,4*
Division of Applied Biology, College of Natural Resources and Life Science, Dong-A University, Busan 604-714, Republic of Korea,1 Division of Applied Life Sciences, Gyeongsang National University, Jinju 660-701, Republic of Korea,2 Biological Function Research Team, Korea Research Institute of Chemical Technology, 100 Jang-Dong, Yuseong-Gu, Daejeon 305-343, Republic of Korea,3 Department of Medical Bioscience, Dong-A University, Busan 604-714, Republic of Korea4
Received 20 August 2007/ Accepted 26 November 2007
|
|
|---|
|
|
|---|
The recognition of microbial diversity and the importance of unculturable bacteria in microbial ecosystems have initiated the screening of unculturable bacterial resources. This requires a genetic approach to exploring microbial resources from both aquatic and terrestrial ecosystems (4, 45). The metagenome approach includes cloning the total microbial genome isolated directly from microbial ecosystems into a culturable surrogate host and attempting to screen for valuable microbial resources (23, 39, 49). The metagenome consists primarily of the unculturable bacterial genome and a minor fraction of the culturable bacterial genome. This is the case because, in most microbial habitats, unculturable bacteria are the dominant members of the total microbial community (57). It has frequently been estimated that 0.001 to 1% of bacteria in soils and aquatic habitats are culturable and characterized as a pure culture (1). Similarly, previous studies with 16S rRNA clone analysis showed that the most abundant species in a forest was from the phylum Acidobacteria, which has few cultured representatives (19, 30, 32).
Compared to the successful discovery of novel enzymes from metagenomes (11, 18, 27, 28, 31, 45, 53), the discovery of antimicrobial compounds is a difficult task, since multiple genes and entire gene clusters involved in the production of the target compounds should be cloned together with resistance and regulatory genes in one contiguous segment (52). This implies that relatively large metagenomic DNA fragments should be cloned and expressed in a surrogate host. The genes encoding the modular enzyme type I polyketide synthase (PKS) frequently encompass DNA fragments of more than 50 kb (23). In many cases, expression of large DNA fragments would be a barrier in a host bacterium such as Escherichia coli. Nevertheless, several papers have reported the discovery of antimicrobial compounds and their respective gene clusters by a metagenomic approach from terrestrial and symbiotic microbial habitats (5, 41, 42). Alternatively, the selection of metagenomic clones based on DNA sequence homology can be conducted either by direct PCR amplification or by screening metagenome libraries with DNA probes. Soil DNA and metagenome clones encoding type I and type II PKSs have been detected from soils by sequence-based screening (16, 50, 55).
In this study, Saccharomyces cerevisiae was used as a target organism to search for antifungal activity from a forest soil metagenome. This fungus is fast growing and is suitable for screening the metagenome library using the double-agar-layer method. Corran et al. (9) reported that 67% of tested fungicides exhibited S. cerevisiae growth inhibition at 5- to 100-µg/ml concentrations and proposed that S. cerevisiae may be an excellent indicator organism to search for candidate fungicidal compounds. Therefore, it may be feasible to use S. cerevisiae as a model fungus to screen for antifungal activity from metagenomes. Here, we describe the isolation of a metagenome clone with antifungal activity using S. cerevisiae as a target fungus. The selected clone was characterized to define the gene cluster involved in antifungal activity expression.
|
|
|---|
Recombinant DNA technology, plasmids, and sequencing clone pEAF66.
Plasmid preparation, restriction endonuclease digestion, DNA ligation, plasmid DNA transformation, agarose gel electrophoresis, and other standard recombinant DNA techniques were carried out following standard methods (46). DNA sequencing and primer synthesis were performed commercially at the DNA sequencing facility of GenoTech Corp. (Daejeon, South Korea). Complete DNA sequencing of the pEAF66 clone was performed by the shotgun sequencing method using plasmid DNA purified from a 500-ml culture of E. coli carrying pEAF66. A shotgun library was constructed by subcloning mechanically sheared DNA (1.5 to 4 kb) into pUC118 digested with SmaI. Insert DNAs were sequenced with a ABI 3730xl DNA analyzer (Applied Biosystems, Foster City, CA). Chromatograms were processed by the phred/phrap/consed software package (http://www.phrap.org). The sequencing gaps were closed by primer walking. The sequence data provided approximately sixfold coverage of pEAF66 DNA. DNA sequences were analyzed with the BLAST program provided by the National Center for Biotechnology Information. Open reading frames (ORFs) were analyzed with the ORF Finder software at the National Center for Biotechnology Information. The following criteria were used to identify the potential ORFs among the potential ORFs detected. Nonoverlapping sequences longer than 50 amino acids were retained. For overlapping ORFs in different reading frames, we selected those with BLAST homologues. The start codon was determined by the presence of conserved Shine-Dalgarno sequences 4 to 10 bp ahead of ATG codon for the individual ORF.
Two genes encoding response regulators were subcloned in high-copy-number plasmids such as pUC119 (59) and pUC129 (26). A 2.58-kb SphI fragment of pEAF66 was subcloned into pUC119 to generate pEJ101. A 1.25-kb XhoI fragment of pEAF66 was subcloned into pUC129 to generate pEJ102. Plasmids pEJ101 and pEJ102 carried genes encoding LysR and IclR family response regulators under the control of the lac promoter, respectively (Fig. 1).
![]() View larger version (16K): [in a new window] |
FIG. 1. Map of a 40-kb fragment of pEAF66. (A) Map of pEAF66 carrying 39 ORFs encoding bacterial type II PKS components (gray arrows) and putative response regulators (dark arrows). Circles at the top of the map indicate the 17 transposon insertion sites causing loss of antifungal activity of the original pEAF66. (B) Enlarged map of orf20 to orf29. orf20 and orf28 (dark arrows) represent potential response regulators in the IclR family and the LysR family, respectively. Two subclones, pEJ101 and pEJ102, from among a high copy number of plasmids are indicated, and the small dark arrow in the subclones represents the lac-operon promoter. Sp and Xh, SphI and XhoI sites for subcloning, respectively.
|
Selection of antifungal metagenome clones.
To select active antifungal clones from the previously constructed and stored forest soil metagenome library (32), the library pool stocks were diluted in buffer (per liter: NaCl, 8.5 g; KH2PO4, 0.3 g; Na2HPO4, 0.6 g; MgSO4, 0.2 g; gelatin, 0.1 g), and metagenomic clones of E. coli were cultured on LB agar supplemented with chloramphenicol for 3 days at 37°C. The number of clones per plate was adjusted to approximately 500 by dilution of the library stock, and at least fivefold the number per initial stock was bioassayed using the procedure described by Rondon et al. (45). Briefly, the S. cerevisiae Y-139 culture suspension was adjusted to an optical density at 600 nm of 0.2 to 0.3 and mixed with LB soft agar at 42°C just before the medium was poured in order to maintain the viability of the yeast. The precultured E. coli clones were overlaid with 5 ml of LB soft agar containing 0.5 ml of S. cerevisiae Y-139 culture. The plates were incubated overnight at 28°C and scored for activity by looking for zones of growth inhibition in the yeast lawn. One antifungal active clone was selected and designated pEAF66. We further analyzed this clone by sequencing the whole insert DNA and transposon insertion mutagenesis. G+C content (%) was determined by using software available at http://tim.saraogtim.com/molbio/gccontent.php.
To investigate the effect of transcriptional regulators for antifungal activity expression of pEAF66, various plasmid constructs were introduced into E. coli EPI-100 carrying pEAF66. E. coli EPI-100 cells carrying pEAF66 and various subclones were grown for 2 days at 37°C on LB medium supplemented with 0.15 mM of isopropylthio-β-D-galactoside, and LB soft agar containing S. cerevisiae was overlaid as described above. Yeast growth inhibition by E. coli EPI-100 carrying various constructs was compared 2 days after incubation.
Transposon mutagenesis of the selected clone.
To define the genes involved in the antifungal activity production from pEAF66, in vitro transposon mutagenesis was carried out following the mutagenesis protocol provided with the GPS-mutagenesis system (New England BioLabs Inc., Beverly, MA). Transposon insertion mutants were selected on LB supplemented with chloramphenicol and kanamycin. The mutants were subjected to antifungal activity analysis with the S. cerevisiae Y-139 strain by the double-agar-layer method. Mutants deficient in antifungal activity were selected, and a transposon insertion site with flanking DNA fragments was subcloned into pUC119. Transposon insertion sites were determined by DNA sequencing and restriction analysis of the mutated subclones.
Antifungal activity analysis.
E. coli cells carrying pEAF66 were grown in 500 ml LB broth supplemented with chloramphenicol in a shaking incubator at 37°C for 2 days. Although the cell-free culture supernatant did not exhibit clear inhibition of S. cerevisiae by double-agar-layer assay, the culture supernatant was extracted twice serially with an equal volume of n-hexane, ethyl acetate, and butanol. The residual aqueous fraction was concentrated to dryness. The organic solvent extracts were also concentrated to dryness. All four residues were dissolved in an appropriate solvent and tested for antifungal activity against S. cerevisiae Y-139. The E. coli cells carrying pEAF66 were also grown in LB agar medium at 37°C for 2 days. The solid cultures were chopped into small pieces and submerged overnight in twice the volume of acetone. Acetone extracts were evaporated to remove the acetone. The residual aqueous fraction was serially extracted with solvents such as n-hexane, ethyl acetate, and butanol. The residual aqueous fraction was concentrated to dryness. The organic solvent extracts were also concentrated to dryness. Residues were dissolved in an appropriate solvent and tested for antifungal activity against S. cerevisiae.
Nucleotide sequence accession number.
The nucleotide sequence of the DNA insert of pEAF66 was deposited in the GenBank database under accession number EU099626.
|
|
|---|
Characterization of the antifungal clone.
Since efforts to isolate and identify active antifungal components produced from clone pEAF66 were not successful, we analyzed the genes responsible for antifungal activity from pEAF66 by transposon insertional mutagenesis and subsequent bioassays of the insertion mutant clones. Among 120 mutant clones tested, 17 were deficient in activity against S. cerevisiae. All 17 mutants also lacked activity against B. subtilis, whereas the original pEAF66 clone inhibited B. subtilis growth. This result indicated that the same locus might be responsible for both antifungal and antibacterial activity. Partial DNA sequences of the transposon insertion sites were determined as described in Materials and Methods, and they were similar to the genes of bacterial type II PKSs or bacterial fatty acid synthases. Therefore, DNA sequences of the full-length insert of pEAF66 were determined and analyzed. The overall G+C content of pEAF66 was 63.7%. We defined the potential ORFs and their organization by a BLAST search (Fig. 1A). The size of the insert DNA is 40,225 bp, and 39 ORFs, including two partial ORFs at both ends of the insert, were identified. The 17 transposon insertion sites were clustered in a potential single transcriptional unit from orf29 to orf35 (Fig. 1A). Two transposon insertions occurred in the promoter area of the antifungal gene cluster. In fact, in silico promoter sequence analysis revealed that the two transposon insertions were exactly in the middle of promoter consensus sequences of the transcriptional unit (data not shown). Therefore, the gene cluster from orf29 to orf35 appeared to be directly responsible for producing the antifungal component. The close relatives of the pEAF66 genes and the identities of deduced amino acid sequences are summarized in Table 1. Since there are no phylogenetic marker genes in the insert DNA, the pEAF66 organism of origin was not known. However, most of the ORFs were highly similar to genes of Proteobacteria.
|
View this table: [in a new window] |
TABLE 1. List of ORFs from antifungal clone pEAF66, gene function, and similar genes in GenBank
|
and KSβ (chain length factor) and the acyl carrier protein (ACP) (34). Homologs of all three core enzyme genes were identified in our gene cluster (Table 1 and Fig. 1A). Our antifungal gene cluster comprised a potential operon with 11 ORFs, since most start codons either overlapped a preceding stop codon or were several bases from the preceding stop codon. The cluster consisted of a gene for ACP (orf29) and four similar genes encoding ketoacyl synthases (from orf30 to orf33) as core enzymes of type II PKS. In addition, genes similar to those encoding aminotransferase, ACP reductase, and ketoreductase (from orf34 to orf36) were identified in the downstream region. Although four genes encoding ketoacyl synthase were detected, they were not highly similar to each other. The highest identity of deduced amino acid sequences among the four genes was 34% between orf30 and orf32 (data not shown). Ketoreductases are probably involved in the hydroxylation of the poly-keto-acyl chain before cyclization, although transposon insertion at orf36 abolishing antifungal activity was not identified. No transposon insertion mutant defective in antifungal activity was identified from orf36 to orf39. We could not find any cyclase or dehydratase homologs in our gene cluster. An identical organization from orf29 to orf39 was found in the operon (PP2777 to 2788) of the P. putida KT2440 genome, which is probably related to the biosynthesis of secondary metabolites such as antibiotics (38). The identity of deduced amino acid sequences between genes of pEAF66 and those of P. putida KT2440 ranged from 31% to 53% (Table 1). A fosmid library of P. putida KT2440 was constructed and subsequently screened for antifungal activity. Screening of more than 20,000 clones did not select any active antifungal clone. It is interesting that a similar gene organization of the region upstream of the antifungal gene cluster was not found in the genome of Gammaproteobacteria such as P. putida KT2440. Instead, the upstream part of the gene organization (from orf1 to orf28) was found in the genomes of Betaproteobacteria, such as Ralstonia solanacearum and Ralstonia metallidurans (Table 1). Whereas the average G+C content of pEAF66 was 63.7%, there were some regions with a relatively AT-rich (50 to 55% G+C) satellite. These were located mostly in the intergenic area, which may be the promoter region (Fig. 2). A long intergenic space (804 bp) found in front of orf29 with a G+C content of 50.6% may be the promoter area for the antifungal gene cluster. The relatively high G+C content (67.8%) in the antifungal gene cluster (orf29 to orf35) was recognized in the area from 27 kb to 36.5 kb (Fig. 2).
![]() View larger version (16K): [in a new window] |
FIG. 2. G+C percentage of a 40-kb insert DNA of pEAF66. Note that the G+C content in the antifungal gene cluster (from 27 to 36.5 kb) is relatively high and the potential promoter region is relatively AT rich.
|
Positive regulation for antifungal activity.
DNA sequence analysis revealed that two genes encoding transcriptional regulators were located upstream of the antifungal gene cluster. orf20 and orf28 encode IclR and LysR family response regulators, respectively. Deduced amino acid sequences of both genes exhibited a conserved DNA binding domain such as the helix-turn-helix motif (data not shown). We subcloned each gene in the high-copy-number plasmid, where the gene expression is driven by the lac promoter of the pUC plasmid. The subclone pEJ101 carried orf28, encoding the LysR family response regulator in pUC119, and pEJ102 carried orf20, encoding the IclR family response regulator in pUC129. Introduction of either pEJ101 or pEJ102 into E. coli carrying pEAF66 apparently increased the yeast growth inhibition zone by the double-agar-layer method compared to E. coli carrying pEAF66 and pUC119 or pUC129 (Fig. 3). The negative control with E. coli carrying pEPI-FOS5 and pEJ102 did not show any yeast growth inhibition (Fig. 3). Similarly, E. coli carrying pEPI-FOS5 and pEJ101 did not exhibit any yeast growth inhibition (data not shown). Repeated experiments with E. coli carrying the various subclones exhibited the consistent results.
|
View larger version (37K): [in a new window] |
FIG. 3. Enhanced antifungal activity of E. coli EPI-100 carrying pEAF66 and genes for transcriptional regulators. S. cerevisiae growth inhibition by E. coli carrying pEAF66 plus pUC129 (A), pEJ102 (B), pUC119 (C), pEJ101 (D), or pEPI-FOS5 and pEJ102 (E) was assessed.
|
|
|
|---|
The active antifungal clone pEAF66 exhibited clear S. cerevisiae growth inhibition. However, the intensive effort to isolate the active component failed when E. coli carrying pEAF66 or E. coli carrying pEAF66 and pEJ101 or pEJ102, which gave increased antifungal activity, was used. We speculated that the antifungal component of pEAF66 might be unstable during our extraction process. Instead of characterizing the component conferring antifungal activity, the DNA sequence analysis of pEAF66 indicated that a gene cluster encoding bacterial type II PKS homologs is involved in the antifungal activity expression in E. coli. Microbial PKSs are involved in the biosynthesis of many natural products of industrial interest and are classified as one of two structurally different types of enzymes (25). Type I PKSs are modular multifunctional enzymes, while the bacterial type II PKSs are multienzyme complexes composed of three or more separate mono- or bifunctional enzymes, which act iteratively during the synthesis of antibiotics. In this study, we identified genes encoding type II PKS homologs, suggesting that the active antifungal component from pEAF66 is probably related to the polyketides. So far, we have been unable to find an extraction procedure that can separate the active antifungal component. Although we do not know the chemical nature of the antifungal component from pEAF66, the gene cluster could serve as a genetic resource for polyketide engineering from unculturable microorganisms (7, 24, 40, 54). Regardless of the difficulty in purifying the antifungal component, it is also possible to speculate that the antifungal component could be the result of a hybrid synthesis, in which some of the enzymes originated from E. coli chromosome, while others were from our clone. Due to the presence of genes encoding three core enzymes for iteratively acting type II PKS, it is likely that E. coli carrying pEAF66 may produce polyketide antifungal compounds. However, without the chemical identity of the active substance, it is difficult to predict the structure of the potential polyketide compounds. Because of four ketoacyl synthases (from orf30 to orf33) and aminotransferase (orf34), the final product will be the result of repeated addition of two carbons to condense malonyl-ACP to a certain chain length of polyketide and subsequent addition of an amine group to the growing acyl chain of the polyketide (3, 24, 25, 40). However, the identification of the final active antifungal components will certainly provide the answer regarding the detailed process of biosynthesis.
Gene organization similar to that of the antifungal gene cluster of pEAF66 was identified from several genome sequences of Gammaproteobacteria such as P. putida KT2440. Thus, it is logical to expect similar antifungal activity against S. cerevisiae from the gene cluster of P. putida KT2440. The P. putida KT2440 operon (PP2777 to PP2788) may be related to the production of antibiotics, although the active component was not isolated nor predicted from strain KT2440 (38). Our screening of more than 20,000 fosmid clones of P. putida KT2440 did not yield any active antifungal clone, suggesting that either the operon (PP2777 to PP2788) of P. putida KT2440 was not expressed in E. coli or the expressed product of the operon does not have antifungal activity. In fact, when precultured P. putida KT2440 was overlaid with LB soft agar containing S. cerevisiae Y-139 culture, clear growth inhibition of the yeast was not detected. This result suggested that even though our antifungal gene cluster showed gene organization similar to that of the P. putida operon, it is likely that the P. putida operon encodes a different activity.
Interestingly, while the downstream part of pEAF66 involved in antifungal activity has a gene organization similar to that of the genome of Gammaproteobacteria, such as P. putida KT2440, the upstream part of the gene organization was similar to genes from the genome of Betaproteobacteria, such as R. solanacearum. This finding raised the question of whether pEAF66 is a chimeric metagenome clone that resulted from the simultaneous ligation of two random DNA fragments originating from Betaproteobacteria and Gammaproteobacteria. Through G+C content analysis, we identified a relatively AT-rich promoter region in front of the antifungal gene cluster of pEAF66. However, G+C content analysis did not clearly indicate whether two different DNA fragments were coligated to generate pEAF66. Interestingly, the IclR family response regulator (orf20) and the LysR family response regulator (orf28) appeared to play roles as positive regulators in the expression of the antifungal gene cluster. Since orf20 is somewhat distal from and not linked to the antifungal gene cluster but is found in betaproteobacterial gene organization, the role of the IclR regulator in positive regulation of antifungal gene cluster expression was surprising. The similar IclR family response regulator and the LysR family response regulator were not found in the upstream part of the operon (PP2777 to PP2788) of P. putida KT2440. The regulators of antifungal gene expression are located in the upstream part, which is similar to the betaproteobacterial genome, while structural genes for antifungal activity are located in the downstream part. Thus, it is likely that the insert DNA of pEAF66 comes from one bacterium.
The IclR family response regulators are widespread among prokaryotes and act as repressors, activators, and proteins with a dual function (36). The IclR family response regulators always have an N-terminal DNA binding domain and a C-terminal effector binding site. Since a clear consensus sequence for the IclR family response regulator has yet to be defined, we could not determine if there is a consensus sequence for IclR protein in the promoter area of the antifungal gene cluster. The LysR family response regulator is one of the most frequently found positive transcriptional regulators in prokaryotes (48). The LysR regulators contained a helix-turn-helix motif and a coinducer recognition domain. Binding of small molecules to the LysR activates transcription of target genes. At this stage, it is not clear whether IclR and LysR directly regulate the expression of the antifungal gene cluster in a cooperative manner or a hierarchical manner. A transposon insertion mutation in either orf20 or orf28 abolishing antifungal activity of pEAF66 was not identified. It is likely that the LysR and IclR families of transcriptional regulators interact with certain chemicals either by an antifungal compound precursor or by foreign signal molecules and subsequently regulate expression of the downstream antifungal gene cluster.
The flanking genes of antifungal clusters of pEAF66 were similar to those for membrane proteins and ABC transporters probably involved in the secretion of antifungal compounds or uptake of primary metabolites for the biosynthesis of the antifungal component. Since transposon insertion mutants deficient in antifungal activity were not obtained from either locus, it is not clear whether they are directly involved in the secretion of antifungal compounds by pEAF66. Although they are involved in the transport process regarding antifungal activity, if E. coli transporters can restore the transport of antifungal activity despite the mutation of orf37 to orf39, the mutation of transporters of pEAF66 would not disrupt the expression of antifungal activity.
S. cerevisiae is a fast-growing fungus that is sensitive to several fungicides. This yeast could be used as an indicator strain to search for candidate fungicides (9). In this study, we assessed the utility of S. cerevisiae in selecting clones with antifungal activity from a metagenome, and we succeeded in screening clones for such activity. On the other hand, the active antifungal component was not isolated and identified from E. coli carrying the antifungal clone. However, molecular characterization of the gene cluster indicated the presence of a novel antifungal gene cluster which may produce a polyketide family component. Further effort is necessary to isolate and identify the antifungal component encoded by the pEAF66 gene cluster. If the antifungal substance is identified as a novel compound, the substance could be used as an antifungal lead compound for novel fungicide development. To our knowledge, this is the first report of an antifungal gene cluster obtained from a soil metagenome.
Published ahead of print on 7 December 2007. ![]()
E. J. Chung and H. K. Lim contributed equally to these studies. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»