Wound Biology Group,1 Department of Oral Surgery, Medicine and Pathology,2 Wound Healing Research Unit,4 Department of Surgery, and Cardiff Institute of Tissue Engineering and Repair, University of Wales College of Medicine, Cardiff CF14 4XY, United Kingdom3
Received 27 June 2003/ Returned for modification 19 September 2003/ Accepted 7 May 2004
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
There have been significant advances in recent years in the characterization and identification of bacteria by molecular techniques. This has facilitated a shift from conventional phenotypic methods for the detection of pathogens in clinical samples to the increased use of molecular approaches for both taxonomic and phylogenetic analysis (13). Conventional analysis of chronic wound communities has relied on traditional bacteriological culture methods (5, 6, 11, 12, 16, 18), thereby completely missing the unculturable population that may be present (53) and thus strengthening the need for the application of such molecular techniques to wound ecology studies.
We applied molecular methods to the analysis of the microflora of a single chronic venous leg ulcer wound (23). Direct sequencing of PCR-amplified 16S rRNA genes demonstrated significantly greater bacterial diversity than that revealed by culture alone. Furthermore, sequences which may possibly represent novel species of bacteria, were retrieved. Hence, the study revealed that a far more comprehensive analysis of the wound microflora was possible by using molecular techniques. However, due to the labor-intensive nature of direct cloning and sequencing to identify the unculturable portion of the microflora, the study was limited to the analysis of a single clinical sample.
In environmental microbiology, denaturing gradient gel electrophoresis (DGGE) has been used as a tool for profiling complex microbial populations without the biases of cultural analysis for many years (30, 32, 51, 52), and the technique has now been applied to the study of a limited number of human microbial populations (17, 21, 28, 45, 55). The advantage of this approach is that it creates a genetic fingerprint or profile of total community diversity by separating mixed 16S rRNA PCR amplification products on the basis of their sequence melting behavior. Subtractive analysis can then be used to identify specific bands of interest such as those sequences that are present only by molecular means (the so-called uncultured fraction) by screening sequences amplified from clinical wound samples alongside those derived from bacteria cultured from these same wounds. These bands can then be cloned and sequenced, without the need to sequence amplification products from cultured organisms as well. DGGE also has other advantages in that direct comparisons can be made between wound samples run on the same gels.
In the present study, we describe the use of DGGE to determine the extent of bacterial diversity within 18 chronic venous leg ulcers, with particular emphasis on the unculturable microflora. Comparison of healing and nonhealing wounds was undertaken to determine whether specific bacterial species were important in the nonhealing phenotype of this painful and often debilitating condition.
| MATERIALS AND METHODS |
|---|
|
|
|---|
0.8. Patients were also screened to exclude the presence of systemic disease (including diabetes), and patients who had received systemic or topical antimicrobial therapy in the previous month were also excluded. All patients in the study were treated with a standard regimen of compression therapy. Patients were followed up for more than 6 months postbiopsy to determine healing rates over this time. Ulcers with 100% closure of the wound were considered healed, but only wounds which did not reoccur within 1 month of wound closure were included in the present study. The study group selected for molecular analysis included 10 healers and 10 nonhealers, although two patients in the healing group were later excluded because the wounds recurred within 1 month of complete healing.
Swab and tissue sample collection. The wound surface was irrigated with 5 ml of sterile saline (0.9% [wt/vol] NaCl), and the surface microflora was sampled by swabbing a 1-cm2 area of the wound by using a cotton tipped swab. The swab was immediately inoculated into 2 ml of transport medium (TM) (3). This medium had been prereduced by incubation in an anaerobic environment (10% CO2, 10% H2, and 80% N2) at 37°C for 24 h. After the swabbing step, local anesthetic, 5 ml of Lidocaine Hydrochloride BP (Phoenix Pharma, Ltd., Gloucester, United Kingdom) was injected subcutaneously. After 5 to 12 min, a specimen from the central region of the ulcer bed (where no healing was observed) was obtained by using a 6-mm disposable sterile punch biopsy (Stiefel Laboratories, Ltd., Sligo, Ireland) and placed immediately into a second vial containing 2 ml of TM. Both samples were transported to the laboratory at room temperature within 2 h of collection for processing. The tissue sample was weighed and bisected with one portion used for traditional cultural analysis and the second was stored at 20°C for subsequent molecular analysis.
Cultural analysis of the wound microflora. The swab sample was vortex mixed in TM for 5 min. The tissue for cultural analysis was cut up finely with a sterile scalpel and also vortex mixed in TM for 5 min. Serial dilutions in 0.85% saline were performed to obtain quantitative culture counts by plating onto the following media (LabM, Bury, United Kingdom): blood agar and Fastidious Anaerobe Agar (FAA), both supplemented with 5% (vol/vol) horse blood, MacConkey no. 3, and Sabouraud medium. All plates except FAA were incubated aerobically at 37°C. FAA plates were incubated under anaerobic conditions at 37°C. Prolonged incubation of the macerated tissue was undertaken in Fastidious Anaerobe Broth for 7 days prior to plating on FAA to allow recovery of fastidious and slow-growing anaerobic species. Primary isolation plates were initially examined after 48 h and then incubated for at least 10 days. Identification of bacteria followed standard microbiological schemes by examination of a range of phenotypic properties (staining reactions, colonial morphology, and carbohydrate fermentation patterns) and, where appropriate, by using commercial identification kits. In this way it was possible to identify each isolate to the genus or species level. Isolates were stored at 80°C until needed.
DNA extraction and manipulation. Molecular analysis was undertaken on the tissue samples only. DNA was extracted by standard techniques (24), with the addition of a freeze-thaw step that was repeated three times (2 min in liquid nitrogen and 2 min at 65°C). DNA was extracted directly from the tissue sample and also from individual microbial species cultured from the same wound.
PCR amplification. The V3 variable region of the 16S rRNA gene was amplified by PCR with primers 341f and 518r (Escherichia coli numbering) as described by Muyzer et al. (32). An improved 40-nucleotide GC clamp was attached to the 5' end of the forward primer, 5'-CGCCCGCCGCGCGCGGCGGGCGGGGCGGGGGCACGGGGGG-3' (15). The previous use of such universal primers is well documented in the literature for amplification of a broad range of bacteria, including their application in a clinical setting (28, 45, 50). The expected product size for the PCR amplification with these primers was a single band 193 bp in length. The PCR mixture (50 µl) contained 50 mM KCl, 10 mM Tris-HCl (pH 9.0), 0.1% Triton X-100, 25 pmol of each primer, 0.2 mM concentrations of each deoxynucleoside triphosphate (Roche, Lewes, United Kingdom), 1 µg of DNA, 1.5 mM MgCl2 (Promega, Southampton, United Kingdom), and 1.25 U of Taq DNA polymerase (Promega). Touchdown PCR was performed with an initial denaturation step of 95°C for 5 min; followed by 20 cycles of 94°C for 1 min, 65°C (0.5°C cycle1) for 45 s, and 72°C for 1 min; followed by 10 cycles of 94°C for 1 min, 55°C for 45 s, and 72°C for 1 min; followed finally by a final elongation step of 72°C for 10 min. Water controls in place of target DNA were carried out for all PCRs. All PCRs were carried out under a laminar flow hood by using filter-sealed tips, ultrapure water, and sterile, clean-room-produced tips and PCR tubes. Molecular weight markers (1 Kb Plus; Invitrogen, Paisley, United Kingdom) were run in the outer lanes of each gel. Gels were stained with 0.5 µg of ethidium bromide (Sigma, Poole, United Kingdom) ml1 and then visualized by using the Molecular Analyst image analysis system (Bio-Rad Laboratories, Ltd., Hemel Hempstead, United Kingdom).
DGGE. DGGE was done by using the Bio-Rad D-CODE system. Four bacterial wound isolates (Peptostreptococcus sp., Proteus mirabilis, Pseudomonas aeruginosa, and Micrococcus sp.) chosen on the basis of their widely differing G+C content (27 to 37%, 38 to 45%, 39 to 59%, and 58 to 70%, respectively) (29) were used to optimize the denaturant gradient on a perpendicular gel, and the isolates were then used as marker strains for all subsequent DGGE gels. PCRs for the "in-house" molecular weight markers were done in bulk and completely separately from all other PCRs. The perpendicular gel contained a denaturing gradient of 0 to 100% and was formed according to the manufacturer's instructions. DGGE optimization demonstrated that the optimal range of denaturants for the control strains was between 30 and 60%.
Parallel DGGE was performed essentially as described previously (31, 32). PCR fragments were separated by using 10% (wt/vol) polyacrylamide (acrylamide-bisacrylamide [37.5:1]; Sigma) containing a 30 to 60% linear gradient of denaturants (urea and formamide) increasing in the direction of electrophoresis. Gradients were formed by using a Bio-Rad Gradient Former model 385. PCR samples were applied to gels in aliquots of 20 µl per lane, with 10 µl of DGGE loading buffer and 10 µl of agarose loading buffer. The running buffer used was 1x TAE buffer (40 mM Tris, 20 mM acetic acid, 50 mM EDTA [pH 8.0]). Electrophoresis was performed at 56°C at 70 V for 10 min, followed by 170 V for 3 h and 50 min. Gels were stained with SYBR Green I (Sigma) diluted in 1x TAE buffer (1:10,000) for 15 min and visualized by using the Gel Doc system (Bio-Rad). The PCR products from all cultured isolates from one wound were electrophoresed alongside those amplified directly from tissue, allowing for direct comparison between cultured and amplified products from the same wound.
Excision and sequencing of DGGE fragments. Bands on the DGGE gel derived from tissue samples that did not align to any bands from bacteria cultured from the same wound were considered unculturable and were excised and sequenced. Bands were cut from the DGGE gels by using sterile scalpel blades and purified by using the "crush-and-soak" method (44). DGGE bands from three cultured wound isolates were also excised and sequenced as controls. The purified DNA was then reamplified by using primers 341f and 534r as described previously. PCR was followed by ligation into the pCR2.1 TOPO vector (Invitrogen) and transformation into Top10-competent Escherichia coli cells (Invitrogen). Blue-white screening of transformants was done on Luria-Bertani agar (44) containing 50 mg of ampicillin (Sigma) ml1 and top spread with 40 ml of X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactopyranoside; Sigma) (20 mg ml1). Clones were screened for inserts of the correct size by M13 amplification. DNA for sequencing was prepared from clones by using Wizard Plus SV Minipreps (Promega). Two clones from each excised band were sequenced, and both strands were sequenced on an automated laser fluorescence sequencer (ABI Prism 377; Applied Biosystems, Warrington, United Kingdom) by using primers M13f and M13r, giving 100% double coverage of the 193-bp 16S ribosomal DNA (rDNA) product for phylogenetic analysis. Sequences obtained were compared to those in the EMBL database (release 71) (47) by using FASTA3 (35, 36) at the European Bioinformatics Institute and with those of the Ribosomal Database Project (RDP) by using SEQUENCE MATCH (10) to identify closely related gene sequences. CHIMERA CHECK (10) at the RDP was used to detect possible chimeric sequence structures, along with manual inspection of the alignments.
Nucleotide sequence accession numbers. Sequences for selected clones in this study have been submitted to EMBL/GenBank under accession numbers AJ605727 to AJ605739.
| RESULTS |
|---|
|
|
|---|
|
|
DGGE.
The DGGE profiles of healing and nonhealing wounds amplified directly from chronic wound tissue were compared in Fig. 1. The individual profiles proved to be complex, with the mean number of bands for the healers and nonhealers being 8 and 10, respectively. Although patients had some bands in common, each patient produced a unique banding pattern that was reproducible between successive PCRs (results not shown). Significantly, however, an apparently identical band was found to be present in all of the samples for both healing and nonhealing wounds (highlighted with an arrow in Fig. 1) and was identified from sequence analysis of seven of these bands as being most closely related to Pseudomonas sp. despite the fact that this species was not cultured from the wounds of all patients (i.e., patients 6, 12, 15, 19, 24, and 28). The shortness of the sequence information obtained for these bands was insufficient to allow for species identification, although all but one band 48a (belonging to the Pseudomonas aeruginosa subgroup) had exactly the same identity to up to 13 species in the Pseudomonas putida and Pseudomonas tolaasii subgroups of the
-proteobacteria. In some cases, the PCR signal for this common band was very weak (i.e., for patients 2, 5, 15, 19, 26, 48, 49, and 53), even though for most of these patients (i.e., patients 2, 5, 26, 48, 49, and 53) high numbers of pseudomonads had been cultured (results not shown). The intensity of the PCR-derived signal for this common band was in general markedly stronger for the nonhealers than for the healers. For a few patients the signal was very faint, although still clearly visible. Also, migration of this fragment in patient 15, was marginally less than for the common bands in all of the other patients. If we assume that the ratio of target human to bacterial DNA does not change significantly between samples and since standard concentrations of target DNA (1 µg) were used in all PCRs, the PCR signal could be regarded as semiquantitative; the strength of signal therefore suggests that the bacterial numbers for this organism may be significantly higher in the nonhealers compared to the healers. From the quantitative cultural analysis of these wounds (results not shown), a large range of Pseudomonas bacterial counts was obtained (8 x 102 to 3.6 x 108 CFU), and these counts revealed no significant differences between healing and nonhealing wounds. Therefore, the PCR results suggest that there may have been an underestimation of bacterial numbers (for Pseudomonas sp. at least) by cultural analysis alone.
|
|
| DISCUSSION |
|---|
|
|
|---|
One of the limitations of DGGE is the length of sequence that can be separated effectively using this technique. The PCR products for DGGE analysis must be <500 bp (33) to allow for efficient resolution and analysis, thereby undoubtedly limiting the amount of sequence information that can be obtained. Although the sequences derived from the DGGE-PCR cloning of products of such small size is insufficient to allow great precision in the construction of phylogenetic trees, it does allow for a presumptive identification at least to genus level (54). The primers used in the present study yielded a very small product (<200 bp), and this was insufficient in many cases to give unambiguous species identification. Better discrimination, however, may be possible with primer products of increased length or from another variable region of the 16S rRNA. In general, DGGE will only separate the predominant species present in a community. In addition, comigration can be a problem for retrieving clean sequences from individual bands. It is also clear from other studies with DGGE that 16S rRNA gene sequences affiliated to specific bacterial species can be found in more than one position in DGGE gels. This is due to the presence of multiple rrN operons. Paenibacillus polymyxa produces more than one DGGE band due to slight sequence heterogeneities between operons (34), and as many as nine different positions were found among 10 Corynebacterium-affiliated sequences (45). There is evidence of multiple operons in the in-house DNA ladder used in the present study in Fig. 1 and 2, particularly since both of these gels were overexposed to enhance the detail of tissue-amplified DNAs. Despite these shortcomings, however, DGGE remains a valuable technique, providing substantial information about complex mixed bacterial populations.
PCR has previously been shown to be more sensitive than culture for the detection of bacteria in clinical samples (40, 45). In the present study, DGGE allowed the identification of a number of strains not detected by cultural means, with ca. 40% of the DGGE fragments sequenced representing organisms not cultured from the wound from which they had been amplified. This highlights the fact that a significant proportion of the resident microflora was not amenable to analysis by culture. In keeping with other studies of the human microflora, it was not that these "molecular isolates" represented unculturable species. In fact, the majority of these amplification products represented strains which, although not cultured from the wound in question, had >95% sequence identity to typical wound microflora organisms, and similar strains were cultured from other wounds in this set of patients. Rather, they were uncultured, perhaps as a result of being in a viable but unculturable state. Moreover, most of these "uncultured" bacteria (i.e., those detected by DGGE alone) were isolated from samples containing large numbers of related bacteria. Hence, they may not have been detected during the initial cultural screening due to competition with other more numerous species and overgrowth on the selective/nonselective media by related microflora. Only one sequence showed <93% sequence identity to any database sequences but this was determined to be a chimeric clone, 48b (S. pyogenes).
Four organisms were identified by sequencing which have not traditionally been associated with chronic wounds. Paenibacillus spp. have been isolated from acute wounds (7), whereas Gemella spp. are closely related to Staphylococcus spp., which were the predominant organisms found in all of the wounds. Sphingomonas and Afipia spp. are closely related to each other, both belonging to the
-division of the proteobacteria, and are species that have not, to our knowledge, been isolated from chronic wounds. Sphingomonas sp. is ubiquitous in the natural environment, has been implicated in community-acquired and nosocomial infections, and was isolated in the present study from both healed and nonhealed patient wounds. Afipia felis is a facultatively intracellular pathogen, although the role of other Afipia sp. as pathogens remains speculative (29).
The contrast between healing and nonhealing wounds was only apparent when looking at the carriage of Micrococcus and Streptococcus spp. in the wounds, both having a higher incidence in the nonhealer group than in the healers. However, although these percentages appear to be potentially significant in this relatively small group of patients, in a larger study in which our group evaluated 66 patients by cultural analysis of the wound microfloras (unpublished results), no significant associations were found between specific bacterial groups and the healing or nonhealing phenotype. Streptococcus and Micrococcus spp. are commensal organisms, but they are also opportunistic pathogens. Both species have been implicated in a range of human infections with Streptococcus sp. in particular initiating (well-documented) soft tissue and skin infections (2) and Micrococcus sp. being implicated in infections in certain immunocompromised patients (29).
Comparison of the amplification signal of the DGGE bands in the nonhealing group of patients to those of the healing group, suggested a possible discrepancy with pseudomonad counts by culture. PCR bias has been widely reported as a reason to view the quantitative interpretation of PCR results with caution (1, 38, 49), although a recent study (8) gives some validity to the quantitative nature of 16S rDNA PCR-DGGE profiles of bacterial communities generated by using this same primer pair (32). Moreover, the particular ratios of human to bacterial DNA are unknown in these samples, and it is these that will significantly influence any PCR biases observed (49). Interestingly, however, bacterial communities in human neonates (17) showed a similar anomaly, with a predominance of Ruminococcus sp. over many months being revealed by DGGE, although these species had been shown by previous cultural analysis to be present in only moderate or low numbers. Recently, real-time PCR was used for the rapid quantification of Pseudomonas aeruginosa in acute wound biopsy samples (37). The application of techniques such as real-time or quantitative PCR (20) to wound DNA would be useful to quantify the pseudomonad (and other species) populations in these wounds to a high degree of accuracy and hence elucidate whether they have a particular significance in nonhealing. Moreover, the sequences generated in the present study can be used directly to aid the design of new primers and probes specifically for such applications. The advantage of such molecular analysis of clinical samples is to further study the relevance and association with healing of selected phylotypes in other samples by DNA probing or PCR analysis by targeting phylotype-specific sequences.
Although a number of strategies such as DNase treatment or restriction endonuclease digestion have been tested for decontaminating PCR reagents, it would seem that none has been wholly successful (9, 25). Hence, it still remains a problem when trying to detect low concentrations of bacterial 16S rDNAs with broad-range PCR primers (22). The fact that we amplified two fragments that appeared to arise from native purified Taq polymerase suggests that a pre-PCR treatment would have been appropriate in the present study to prevent universal primers from amplifying contaminating bacterial DNA. Hence, additional strategies, along with those currently used (i.e., clean-room produced plastics, filter tips, dedicated PCR room, and ultrapure DNase- and RNase-treated water), such as DNase treatment of Taq polymerase (and perhaps other PCR reagents), should be recommended for amplification of bacterial DNA from such clinical samples to improve the validity of universal 16S rRNA gene PCR results by prevention of false positives (22).
The use of multiple primers in conjunction with PCR and sequencing to analyze clinical samples has been highlighted previously in the study of dental (14, 26) and gastrointestinal (48) microflora and is also now beginning to be applied to the field of wound microbiology (41). In a previous study of a single chronic venous leg ulcer (23), we found that using multiple primer sets in conjunction with the sequencing of large numbers of clones, a far more detailed picture of the wound microflora was revealed. Indeed, a far greater diversity of organisms was identified in this way.
The results of the present study together with our previous work (23) have demonstrated the greater diversity of the wound microflora as assessed by a molecular approach in comparison to culture alone. We have also shown that, by the use of DGGE profiling, the molecular analysis of the human microflora need not be limited to small numbers of samples as would be the case for the cloning and sequencing approach. By extension, this labor-saving technique offers the possibility of new clinical approaches, for example, monitoring the status of wounds by serial DGGE profiling. Although the additional complexity documented by a molecular approach could be seen to further complicate the study of these chronic wounds, it will only be by the comprehensive microbial analysis of clinical specimens and the thorough investigation of the biases of both cultural and molecular approaches that associations between specific bacterial types and healing will be established or refuted.
| ACKNOWLEDGMENTS |
|---|
We thank Gareth Lewis for sequencing.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|