Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2004, p. 6619-6627, Vol. 70, No. 11
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.11.6619-6627.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Institute of Biological Sciences, University of Wales, Aberystwyth, Ceredigion, Wales,1 Department of Chemistry, University of Manchester Institute of Sciences and Technology, Manchester, United Kingdom2
Received 25 February 2004/ Accepted 29 June 2004
|
|
|---|
|
|
|---|
The ability to measure bacterial diversity is prerequisite for the systematic study of bacterial biogeography and community assembly. It is therefore central to the ecology of surface waters, the oceans and soils, waste treatment, agriculture, and global elemental cycles. However, the experimental definition of bacterial diversity has never been undertaken for any naturally occurring bacterial community, and the extent of prokaryotic diversity is widely held to be beyond practical calculation (7). While the gold standard environmental diversity metric is currently based on 16S rDNA sequences (49), it is widely recognized that polyphasic taxonomic approaches are better, since they take into account all available phenotypic and genotypic data and integrate them in a consensus type of classification (44).
The successful isolation and enumeration of actinomycetes from the environment is usually achieved by dilution plate techniques with a medium containing selective nutrients and certain antibiotics (16, 40). The rapid characterization and identification of microbiological isolates are prerequisites for an effective industrial screening program. Detection of duplicates is necessary to reduce redundancy in screening assays, a process referred to as dereplication (3), and the selection of unusual strains can help in the discovery of novel compounds. Phylogenetic characterization of isolates on the basis of 16S rDNA and rRNA sequences and DNA-DNA homology are considered by many microbial taxonomists to be the most robust approaches to microbial identification at the species level. Based on 16S rRNA sequence analysis, Stackebrandt et al. (42) suggested the reclassification at the family level of actinomycetes and found that there were taxon-specific 16S rRNA and rDNA signature nucleotides.
The application of both traditional techniques, such as macro- and micromorphological, chemosystematic, and numerical taxonomy of phenotypic measurements, and the modern molecular systematic techniques were instrumental in clarifying the relationships within the actinomycetes (2, 42), although this process is still unfinished (11). However, 16S rDNA sequencing is time-consuming, and its capacity to differentiate at the subspecies level is often insufficient. As a result, a variety of alternative approaches have been developed for characterizing actinomycetes. These include fatty acid methyl esters (19), fatty acid composition by gas chromatography (9), pyrolysis mass spectrometry of cells (4, 41), random amplified polymorphic DNA method (1), and rep-PCR (51).
Amplified fragment length polymorphism (AFLP) is a relatively new genomic fingerprinting technique that is based on the selective amplification of a subset of restriction fragments (22, 46). There is plenty of evidence to show that AFLP can be applied in bacterial taxonomy for the differentiation of closely related taxa and for diagnostic, epidemiological, and source-tracking applications. Indeed, its particular power lies in the fact that it can discriminate bacteria down to the strain level (5, 13, 21). In addition, fluorescent AFLP analysis can be easily automated, digitized, and standardized for long-term database establishment and even exchange of data between laboratories.
Fourier transform infrared spectroscopy (FT-IR) is a rapid, nondestructive spectroscopic approach for whole-organism fingerprinting (10, 18, 33, 38, 39). FT-IR measures vibrations of functional groups and highly polar bonds (14). Thus, these (bio)chemical fingerprints are made up of the vibrational features of all the sample components, and for microbial samples these will include DNA and RNA, proteins, and membrane and cell wall components. FT-IR spectrometers record the interaction of IR radiation with samples, measuring the frequencies at which the sample absorbs the radiation and the intensities of these absorptions. The resolving power of FT-IR is excellent, and it has even been shown to detect differences between yeast mutants (39) and in the metabolic footprint of microbial cells differing in the presence of single genes (25). Our established protocol for FT-IR (10, 35, 43, 47, 48) has the major advantages that it is nondestructive, reproducible, and very rapid both for a single sample (1 to 10 s) and with respect to the automated high throughput of samples in batches of 384.
In this study, the main aim was to investigate the use of FT-IR spectroscopy to characterize 20 coastal sediment isolates putatively identified by morphological and cultural characteristics as Micromonospora. The same isolates were also analyzed by 16S rDNA sequencing and AFLP, and the phylogenetic results of all three methods are compared.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Isolation of actinomycetes from marine sediments
|
Nearly complete (
1,400 nucleotides) 16S rDNA sequences of all isolates were aligned with Streptomyces coelicolor A(3)2 (designated Y00411 [which was retrieved from the EMBL database]) and a selection of Micromonospora strains (also retrieved from EMBL) with ClustalW (20). After automatic alignment, the sequences were manually checked. TreeCon version 1.3b software was used to estimate evolutionary distances (Jukes and Cantor, 1969), and to construct phylogenetic trees by the neighbor-joining method (45).
AFLP analysis.
Restriction, ligation, and amplification were performed as described previously (22). Approximately 250 ng of chromosomal DNA was digested with BamHI and MspI (Promega) and subsequently ligated with T4 DNA ligase (Promega) to double-stranded restriction-half-site-specific adaptors (1 µM).
The double-stranded restriction-half-site-specific adaptors were made as follows. For the BamHI adaptor, 100 µl each of the complementary oligonucleotides 5'-TCGTAGACTGCGTACA-3' and 5'-GATCTGTACGCAGTCTACGA-3', each at a concentration of 100 µM, were used, and for the MspI adaptor, the complementary oligonucleotides 5'-CGTCAGGACTCATCGTC-3' and 5'-GACGATGAGTCCTGA-3' were employed at the same concentration and details as for the BamHI adaptor. The two BamHI adaptor complementary single-stranded oligonucleotides were mixed with 14 µl of 5 M NaCl, heated at 95°C for 3 min, then cooled down at room temperature, and incubated at 25°C overnight. Then 20 µl of 3 M sodium acetate (pH 5.2) and 200 µl of pure 2-isopropanol were added to precipitate the constructed adaptor DNA. TE (10 mM Tris-HCl, 1 mM EDTA, pH 8.0) was then used to dissolve the above DNA (100 µM). The final concentrations of MspI adaptor and BamHI adaptor for ligation of restriction fragment were 6.6 and 2.3 µM, respectively.
For PCR amplification of the ligated fragments, the BamHI+T primer (5' modified with fluorescence 7',8'-benzo-5'-fluoro-2',4,7-trichloro-5-carboxyfluorescein, 5'-TCGTAGACTGCGTACAGATCCT-3') and MspI+0 primer (5'-GACGATGAGTCCTGACGG-3') were used. Amplified fragments were separated on an ABI 3100 DNA sequencer (Applied Biosystems) with a 36-cm array, filter set D, and 400HD (Rox) standard as the internal molecular size marker.
The AFLP results were captured with ABI Prism Genescan 3.7 software (Applied Biosystems), and the fragment data were tabulated by size and fluorescent intensity with the ABI Prism Genotyper 3.7 software (Applied Biosystems). All electropherograms were visually inspected for polymorphous peaks before the final fragment table was produced. Only fragments in the range from 50 to 500 bp were considered further because the resolving capability of the sequencing gel is generally good in this range, when with the Rox internal ladder. Analysis of fluorescent AFLP data was performed with the TreeCon version 1.3b software (45), and an unrooted phylogenetic tree was produced with simple matching and unweighted pair group mean average (UPGMA).
FT-IR spectroscopy.
The 20 Micromonospora isolates were inoculated onto the surface of nylon membranes placed on top of nonsporulation agar medium (NSM; Casamino Acids 2%, soluble starch 2%, yeast extract 0.4%, agar 1.2%, pH 7.2) and incubated at 25°C for 7 days. Each isolate was cultured six times, giving six biological replicates for analysis. After incubation, the biomass was collected carefully and suspended in 0.9% NaCl and then homogenized by 10 s of sonication. Homogenized biomass was loaded into 400-well aluminum IR sample carriers and subsequently dried at 50°C for 30 min. Each biological replicate was analyzed three times (so-called machine replicates).
The FT-IR instrument used was a Bruker IFS28 FT-IR spectrometer (Bruker Spectrospin Ltd., Banner Lane, Coventry, United Kingdom) equipped with an MCT (mercury-cadmium-telluride) detector cooled with liquid N2. The aluminum plate was then loaded onto the motorized stage of a reflectance thin-layer chromatography accessory and analyzed as detailed (10, 35, 43, 47). The IBM-compatible personal computer used to control the IFS28 was also programmed (with Opus version 2.1 software running under IBM O/S2 Warp provided by the manufacturers) to collect spectra over the wave number range 4,000 cm1 to 500 cm1. Spectra were acquired at a rate of 20 s1. The spectral resolution used was 4 cm1.
To improve the signal-to-noise ratio, 256 scans were coadded and averaged. All spectra were collected in reflectance mode and displayed in terms of absorbance as calculated from the reflectance-absorbance spectra with the Opus software, which is based on the Kubelka-Munk theory (14). Typical FT-IR spectra are shown in Fig. 1. Ranges of wave numbers are informative and can be associated with special chemical bonds (17, 36); 3,050 to 2,800 cm1 is the so-called fatty acid region, where peaks show the vibration of CH2 and CH3 groups of fatty acids. The amide section, 1,750 to 1,500 cm1, is where protein and peptide bands dominate. The range from 1,500 to 1,200 cm1 is a mixed region containing vibrations of fatty acids, proteins, and polysaccharide. Polysaccharide dominates the region from 1,200 to 900 cm1. The fingerprint region is from 900 to 700 cm1. This contains bands which are most characteristic at the species level, but only a few peaks can be assigned to the vibrations of special substances.
![]() View larger version (25K): [in a new window] |
FIG. 1. Typical FT-IR spectra from isolates N0072 and N0073. The spectra are offset on the ordinate for clarity.
|
Since the six biological replicates from each isolate clustered on top of one another (data not shown), the means of these were used for dendrogram construction. The Euclidean distance between a priori group centers in discriminant function analysis space from the biological replicates was used to construct a similarity measure, with the Gower general similarity coefficient SG (12), and these distance measures were then processed by an agglomerative clustering algorithm to construct a dendrogram (32).
Nucleotide sequence accession numbers.
The 16S rDNA sequences for the 20 Micromonospora isolates have been deposited in GenBank, and the sequence accession numbers are from AY221484 to AY221503, as detailed in Table 2.
|
View this table: [in a new window] |
TABLE 2. Micromonospora isolates studied and top three 16S rDNA sequence hits against the RDPII database
|
|
|
|---|
Many colonies with different morphological and cultural characteristics were picked from the 1- to 4-week-old isolation plates and transferred to starch casein agar medium for purification and morphology observation. Most of the isolates were tentatively assigned as belonging to the genera Streptomyces and Micromonospora according to their macro- and micromorphological and cultural characteristics. Streptomyces isolates had rich aerial mycelia and matured to form short, long, spiral, or straight spore chains. Most of their colonies showed fuzzy surface because of sporulation. Compared with the Streptomyces isolates, Micromonospora spp. have smaller colonies with yellow, pink, or purple vegetative mycelium and sporulate after longer incubation. Some rare actinomycetes, such as Nocardia spp. with branched substrate mycelium, were also tentatively identified.
While an investigation of the exact community diversity and structure of this marine sediment would be useful, this was not our current objective. It would not be practical to pursue this in the present study, since only cultured organisms are considered and sampling errors will be manifest in taking such a small subsample from such a large beach area. We therefore considered whether FT-IR and AFLP can be used as complementary tools to 16S rDNA sequencing for the characterization of a subset of these actinomycetes.
Twenty Micromonospora isolates were selected for analysis by a polyphasic taxonomic approach. These were coded as N0072, N0073, N0074, N0075, N0087, N0088, N0093, N0097, N0098, N0152, N0153, N0155, N0156, N0157, N0158, N0161, N0162, N0163, N0168, and N0179. All 20 Micromonospora isolates gave rise to relatively small colonies on starch casein agar, and the substrate mycelia were orange or purplish. On further incubation, the colony surface of some of them became black as they started to sporulate. A typical scanning electron micrograph of one of the isolates is shown in Fig. 2, where it can be seen that spores were generated singly from the substrate mycelium.
![]() View larger version (162K): [in a new window] |
FIG. 2. Scanning electron micrograph of isolate N0072 grown on starch casein agar for 3 weeks at 25°C. Single spores were produced from the substrate mycelium.
|
Nearly complete 16S rDNA sequences from the 20 isolates were amplified and sequenced. The sequences were compared with all the 16S rDNA data stored in the Ribosomal Database Project II (RDP II) with a fast Sequence Match search program (31). All 20 isolates were placed in the genus Micromonospora according to the hit lists of similar sequences found in RDP II (Table 2).
A selection of Micromonospora spp. from RDP II were downloaded and aligned with the 20 isolates from this study along with the 16S rDNA sequence from S. coelicolor A(3)2 Y00411 (the outgroup). All of these Micromonospora strains exhibited a very high level of 16S rDNA sequence similarity (>97%, data not shown). A phylogenetic tree was generated with the neighbor-joining method, and the results are shown in Fig. 3.
![]() View larger version (20K): [in a new window] |
FIG. 3. Neighbor-joining phylogenetic tree constructed with 20 Micromonospora isolates from this study (denoted N072 to N0179) plus Micromonospora from the RDPII based on their 16S rDNA sequences. S. coelicolor A(3)2 Y00411 was used as the outgroup. Bootstrap values greater than 50% are indicated at the nodes (100 replications). The bar stands for 0.1 substitution per nucleotide position.
|
In order to simplify the comparison of this phylogeny based on 16S rDNA sequences with those generated from AFLP and FT-IR spectroscopy, phylogenetic analysis of just the 20 isolates, with Y00411 as the outgroup, was undertaken. This led to the production of the tree shown in Fig. 4, which clearly shows that the strains fall into two subgroups (labeled A and B).
![]() View larger version (12K): [in a new window] |
FIG. 4. Neighbor-joining phylogenetic tree constructed with 20 Micromonospora isolates based on their 16S rDNA sequences. Streptomyces coelicolor A(3)2 Y00411 was used as the outgroup. Bootstrap values greater than 50% are indicated at the nodes (100 replications). The bar stands for 0.1 substitution per nucleotide position.
|
100 fragments which ranged from 50 to 500 bp and the fragments were compiled into a data matrix in Microsoft Excel for visual inspection prior to cluster analysis. The frequency of fluorescent AFLP bands shared between the 20 Micromonospora isolates analyzed was assessed, and the results are shown in Table 3. It can be seen in this table that no single band was present in all 20 isolates, and moreover, the majority of the isolates had few common bands. That many of the isolates had unique fluorescent AFLP barcodes with few shared polymorphisms will have implications for phylogenetic analyses (i.e., the branches on the trees will be long and the bootstrap values are likely to be low).
|
View this table: [in a new window] |
TABLE 3. Frequency of fluorescent AFLP bands shared among the 20 Micromonospora isolates analyzed
|
![]() View larger version (47K): [in a new window] |
FIG. 5. Cluster analysis of the fluorescent AFLP data from the 20 Micromonospora isolates with simple matching and UPGMA. Also shown are the AFLP profiles for these isolates.
|
FT-IR spectroscopy was performed on the six biological replicates from each of the 20 isolates. Therefore, 120 samples were cultivated, and each of these was analyzed three times (the machine replicates). The cluster analysis of these data by the discriminant function analysis algorithm used a priori group knowledge encoded at the machine replicate level; that is to say, 120 groups were used rather than 20, allowing an unbiased analysis and allowing us to ascertain how reproducible growth and sample preparation were. Considering that 360 objects were clustered in this discriminant function analysis plot (data not shown), it is necessarily cluttered; however, each of the six biological replicates clustered synonymously for each of the 20 Micromonospora isolates. Therefore, we were satisfied that the growth and preparation of these actinomycetes were conducted in a reproducible manner.
The dendrogram was therefore constructed with the group means in the discriminant function analysis plot from the biological replicates, as detailed above, and is shown in Fig. 6. The clustering of the Micromonospora isolates from their FT-IR spectra (Fig. 6) was in good agreement with 16S rDNA sequencing analysis. The two main groups, A and B, were the same, although the subclustering within these two groups was slightly different. Isolates N0163 and N0168 have 100% 16S rDNA similarity and are clustered together based on FT-IR spectral information. This shows that they are duplicates, at least on the species level. The nearly complete (
1,400 nt) 16S rDNA sequence similarity between N0073 and N0075, N0152 and N0155, N0072 and N0161, N0097 and N0098, and N0093 and N0153 is 99.6%, 99.6%, 99.4%, 99.9%, 99.4%, respectively, and they are grouped together in FT-IR clustering analysis. Isolates N0087, N0088, and N0157 have identical 16S rDNA sequences and N0153 has 99% similarity with N0087 and N0157, which means that they are the same species. The FT-IR dendrogram also support this grouping, that those five isolates are in the same group, A. All the isolates from the subgroups in groups A and B have more than 98.2% 16S rDNA similarity. That means that FT-IR can detect duplicates on the species level. This result is highly encouraging when one considers that FT-IR is a phenotypic typing method and measures the total biochemical makeup of the sample rather than being based on a single gene.
![]() View larger version (16K): [in a new window] |
FIG. 6. Dendrogram constructed from the FT-IR spectra from the 20 Micromonospora isolates with hierarchical cluster analysis as detailed in the text.
|
Concluding remarks.
FT-IR spectroscopy is becoming established for the rapid identification of microorganisms (for a review see reference 33). The interest in FT-IR for the characterization of actinomycetes is growing, and strains from several genera have been analyzed, including Streptomyces and Micromonospora (15) and Corynebacterium and Rhodococcus (37). One of the exciting prospects is the generation of validated comprehensive FT-IR spectral reference databases, and this has been demonstrated in food-borne yeasts (27) and some selected coryneform bacteria and related taxa from the suborders Micrococcineae and Corynebacterineae (Actinomycetales and Actinobacteria) (38).
Oberrueuter and colleagues (37) compared their FT-IR results with 16S rDNA sequence analysis to assay the intraspecific diversity of three different actinomycetes species, Brevibacterium linens, Corynebacterium glutamicum, and Rhodococcus erythropolis. No correlation was found between 16S rDNA sequence similarity and FT-IR spectral distance for the three species below the species level. In our study, the characterization of a group of 20 marine sediment isolates belonging to Micromonospora by FT-IR was compared to two molecular methods, the established 16S rDNA sequence analysis and AFLP, which is itself a relatively new molecular typing approach. All methods gave very similar dendrograms, and FT-IR and 16S rDNA were found to be the most congruent, at least on the species level. We also studied 15 Streptomyces isolates with both the FT-IR and 16S rDNA sequence analysis methods (unpublished data), which showed that the FT-IR spectral dendrogram correlated to the phylogeny tree very well. Thus, FT-IR could be used as preliminary dereplication for screening novel isolates.
All three methods in this study characterize bacteria in different ways. 16S rDNA analysis focuses on a small highly conserved region of the genome, which, while useful for inferring deep phylogenies, suffers in that only a small part of the genome is analyzed. By contrast, AFLP analyzes DNA fragments from restriction analysis of genomic DNA and so more of the genotype is interrogated. FT-IR spectroscopy also effectively analyzes more of the organism, although at the phenotypic level, by measuring the vibrational properties of molecules.
FT-IR is a very rapid method, typically taking 1 to 10 s per sample, and samples can be loaded in batches of 96 or 384. It is very easy to use and inexpensive with respect to laboratory consumables (although an initial capital outlay for the equipment is needed). By contrast, both 16S rDNA sequencing and AFLP are slower (taking many hours if not days), require skilled operators, and are expensive at both the consumables and instrument levels. It is worth noting that a DNA sequencer is equivalent in price to an FT-IR instrument. Since FT-IR spectroscopy gives a measure of the organism's phenotype, standard growth parameters such as medium, temperature, and incubation time need to be controlled stringently in order to get reproducible results, especially for construction and comparison of spectral databases.
In conclusion, although currently FT-IR cannot be used alone for microbial classification because extensive microbial databases do not currently exist, we believe that FT-IR will become a very valuable first step in microbial screening programs. In particular, it can be used as a reliable method for the dereplication of environmental isolates prior to more expensive and time-consuming characterization of organisms producing potentially important and novel pharmacophores.
We are grateful to Steve Wade for scanning electron microscopic analysis.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»