|
|
||||||||
1,2Department of Biomedical Sciences1 and Gynaecology and Obstetrics Unit2, Faculty of Medicine and Health Sciences, University Putra Malaysia, Selangor, Malaysia 3Department of Medical Microbiology, Faculty of Medicine, University of Malaya, 59100 Kuala Lumpur, Malaysia
Correspondence Pei Pei Chong cpp{at}medic.upm.edu.my
Received May 20, 2002
Accepted January 28, 2003
| Abstract |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Recurrent vaginal candidiasis (RVC) can be defined as the occurrence of at least four mycologically proven symptomatic episodes within 12 months (Sobel, 1992), or at least three episodes unrelated to antibiotic therapy within 1 year (Ringdahl, 2000). Two main hypotheses have been proposed for the cause of recurrent vaginitis: (i) reinfection through sexual transmission or from reservoirs in other organ systems (e.g. gastrointestinal or urinary tracts) (Fidel & Sobel, 1996); (ii) relapse due to incomplete eradication of Candida species that have increased virulence or are drug-resistant (Sobel, 1993; Faro, 1996). Decreased cell-mediated local immunity is also a risk factor for RVC (Fidel & Sobel, 1996).
In a study on 18 patients with RVC, Lockhart et al. (1996) used three probes to evaluate the genetic relatedness of Candida strains; they reported that the predominant scenario was strain maintenance with minor genetic variations, a process that they termed substrain shuffling'. Vazquez et al. (1994) found that eight of ten patients maintained a single strain of C. albicans through sequential infectious episodes but, in two patients, the strain had undergone replacement by another strain. In a similar study, Mercure et al. (1993) demonstrated that the majority of patients had maintenance of a single strain, whereas a few patients had strain replacement. These studies were focused only on C. albicans, whereas patients with non-albicans infections were excluded. Therefore, the present study aimed to identify the different Candida species that give rise to recurrent infections in a representative Malaysian population and to assess their genetic relatedness by using random amplification of polymorphic DNA (RAPD).
| METHODS |
|---|
|
|
|---|
Vaginal secretions were cultured for 27 days at 37 °C on Sabouraud's dextrose agar (Difco) alone, and also supplemented with chloramphenicol (final concentration, 50 µg ml-1). Morphologically distinct colonies from each culture were subcultured and stored on Sabouraud's dextrose slant for species differentiation, DNA fingerprinting and storage. C. albicans was identified by using the germ-tube test (Rippon, 1974). Candida species were also identified by using the biochemical sugar assimilation test (Rippon, 1974) as well as the API 20C AUX yeast identification kit (bioMérieux).
DNA extraction and RAPD-PCR.
Genomic DNA was isolated from each Candida strain for genotyping. A loopful of a single colony of Candida was added to 2.0 ml Sabouraud's dextrose broth and incubated with shaking at 37 °C for 48 h. Yeast cells were collected by centrifugation, resuspended in 1.0 ml PBS and washed twice with PBS. The cell pellet was resuspended in 500 µl lysis buffer [50 mM Tris/HCl (pH 7.5), 10 mM EDTA, 0.5 % ß-mercaptoethanol and 0.2 mg lyticase ml-1 (Sigma)] and incubated at 37 °C for 30 min on a rotary shaker (150 r.p.m.). Lysed cells were treated with 25 µl 20 % (v/v) SDS and 2.5 µl proteinase K (20 mg ml-1), incubated at 56 °C for 60 min and then boiled for 5 min. DNA was extracted with an equal volume of phenol/chloroform/isoamyl alcohol (25 : 24 : 1, pH 8.0; Pierce Biotechnology). DNA was then precipitated with 0.1 vol 3 M sodium acetate (pH 5.2) and an equal volume of 2-propanol; the resultant pellet was washed with cold 70 % (v/v) ethanol and finally resuspended in 20 µl sterile distilled water.
An arbitrary oligonucleotide primer, PA03 (5'-AGTCAGCCAC-3'; custom-made by Operon Technologies) was used for DNA fingerprinting, using the conditions described by Xu et al. (1999) with slight modifications. DNA amplifications were performed in 25 µl volumes that contained 1040 ng DNA template, 2.5 µl 10x PCR buffer, 1.5 µl 25 mM MgCl2, 20 pmol PA03 primer, 2.5 µl 0.5 mM dNTP mixture and 0.5 U SuperTherm DNA polymerase (Bertec Enterprise). For the negative control, sterile distilled water was used as substitute for the DNA template. RAPD-PCR was carried out by using a thermal cycler (Biometra) with an initial denaturation step of 97 °C for 3 min, followed by 45 cycles of 60 s at 93 °C, 60 s at 36 °C and 120 s at 72 °C, and a final extension cycle of 5 min at 72 °C. Amplicons were separated by electrophoresis in 1.5 % agarose gel with added ethidium bromide (final concn, 0.5 µg ml-1) in 1x TBE buffer for 4 h at 2 V cm-1, using the Maxicell Primo horizontal gel electrophoresis system (EC Apparatus). Bands were visualized by using a UV-light transilluminator and pictures were captured digitally with the Gene Genius BioImaging system (Syngene).
Analysis of RAPD profiles and statistical test.
The RAPD technique produced a mixture of fragments from each DNA sample that was compared. When these fragments were separated by gel electrophoresis in an immobile gel, a characteristic pattern of bands was obtained from each sample. When fragment mixtures from two or more samples were compared, the total number of unique bands produced by a particular DNA sample compared to all other samples was determined. To compare DNA samples that were run in different gels or in non-adjacent lanes of the same gel, the molecular sizes of all bands were computed by using the GeneTool software (Syngene); bands with highly similar or identical sizes were deemed to be homologous. Subsequently, the presence or absence of every band in each pattern was recorded in a chart. Relatedness of the DNA samples was assessed using the RAPDistance program (Armstrong et al., 1994). The Dice metric (Dice, 1945) was the statistical method chosen to calculate the similarity coefficients (SABs) of sample pairs. SAB was calculated as follows: SAB = 2n11/(2n11+n01+n10), where n11 represents the number of bands shared by samples A and B, n01 represents the number of bands present in sample A but not in sample B and n10 represents the number of bands present in sample B but absent in sample A.
Dendrograms based on a pairwise distance (1-SAB) matrix of the DNA samples (which determines within- and between-group similarities) were generated by UPGMA. The MEGA software package, version 2.1 (Kumar et al., 2001) was used to generate the UPGMA phenograms. A paired t-test was used to compare the mean SAB values among defined groups.
In order to compare the distribution of species among different groups of patients, species diversity for each group was determined. Species diversity was calculated as (1-
ps2), where ps represents the frequency of a particular species (Xu et al., 1999). The species diversity ranges from 0, where all isolates are of the same species, to 1, where every isolate is of a different species.
| RESULTS |
|---|
|
|
|---|
|
|
|
Interspecies RAPD profiles
For 87 isolates from 71 patients, 61 polymorphic PCR fragments of distinct sizes were detected and scored. However, only 80 RAPD profiles were obtained (seven of the isolates failed to produce any PCR amplicons, despite repeated attempts). Fig. 1 parts (a), (b) and (c) are representative of the electrophoresis of PCR products obtained from isolates of patients in groups A, B and C, respectively.
|
Substantial interspecies differences in RAPD profiles were observed when isolates of C. albicans, C. glabrata, C. lusitaniae, C. famata, C. krusei and C. parapsilosis were compared, as evident from Fig. 1. The number and size of DNA fragments generated were fairly characteristic for each species, especially C. glabrata, C. lusitaniae and C. krusei. However, C.albicans, C. famata and C. parapsilosis had similar patterns to each other. Despite the relatively similar array of species-specific DNA fragments, distinct profiles generated by RAPD analysis were observed for C. albicans isolates, as shown by that of HRY2 compared with that of ZKH1 (Fig. 1a). The former generated bands that ranged from 200 to 1400 bp, whereas the latter produced bands that ranged from 200 to 1700 bp. Most C. albicans isolates produced bands that resembled the pattern of HRY2. However, some C. albicans isolates only produced between one and three bands. C. krusei isolates also produced two different profiles, as seen in those of GCL2 and TVM.
Genetic similarity of Candida species within patients in group A
Among the 12 patients in group A, five sequential Candida isolates were obtained from one patient (pTHN; THN15), three isolates from another patient (pLML; LML13) and two isolates from each of the remaining ten patients. To assess the genetic similarity of sequential isolates within an individual, SAB values between each pair of isolates were computed and tabulated in Table 4.
|
As a quantitative measure of the level of unrelatedness of isolates, a collection of unrelated isolates (each obtained from different individuals) must be used as a control collection for determination of the threshold of unrelatedness (Marco et al., 1999). For the present study, 25 unrelated C. albicans isolates (each collected from a different patient in group C) were used as the control collection. The mean SAB value for this control collection was 0.60 ± 0.20, which represents an estimate of unrelatedness for vaginal isolates that were used in this study. Based on the empirical estimation of unrelatedness derived from the control group, an SAB value of 1.00 represents identical strains, an SAB value of 0.800.99 represents highly similar (but non-identical) strains and usually reflects microevolution of a single strain, and an SAB value of < 0.80 represents unrelated strains.
Pairwise similarities of only C. albicans isolates (n = 17) in group A patients were compared. It was found that 12.5 % of C. albicans isolates were genotypically identical, 56.25 % were highly similar (but non-identical) while 31.25 % of the isolates were unrelated. Fig. 2 presents the genetic similarity of all 28 isolates of different Candida species from group A patients in a UPGMA phenogram. This figure shows that isolates of the same species clustered together, with one exception: isolate ZKH1, which was identified morphologically and biochemically as C. albicans, was more closely related to C. glabrata and C. krusei in the phenogram.
|
The three isolates obtained from patient pLML were genotypically identical to each other, with an SAB value of 1.00. The isolates were also clustered together on the phylogenetic tree (Fig. 2). The most closely related but non-identical strains from two different isolates of a single patient were from patients pELZ, pLSL and pMHS, as both of their strains of Candida species clustered together on the phylogenetic tree. The two different isolates of patients pNFZ, pLSY and pNML were also related to each other, although not as closely related as the isolates from pELZ, pLSL and pMHS.
For patient pTHN, who had five consecutive isolates, the first four isolates (THN1, THN2, THN3 and THN4) were fairly closely related to each other, whereas the fifth isolate, THN5, was totally unrelated to the first four isolates (SAB < 0.8). This was because this patient was colonized by a different Candida species during the fifth episode of vaginitis compared with the first four episodes of infection. For patient pGCL, the second isolate obtained was of a different Candida species (C. krusei) from the first isolate (C. albicans).
Although conventional methods for the identification of Candida species showed that some of the sequential isolates were phenotypically identical, RAPD analysis revealed them to be genotypically different, as observed for ZKH1 vs ZKH2, HRY1 vs HRY2 and MHR1 vs MHR2. Analysis of C. albicans isolates gave rise to the observation of two distinct groups in the phenogram. The first group consisted of ZKH1, the only isolate that was on its own and did not cluster with any other isolate. The second group comprised two subgroups, with MHR1 and MHR2 in one subgroup and all remaining C. albicans isolates in the other.
In summary, for five patients (pLML, pELZ, pMHS, pLSL and pMHR), subsequent infection was caused by identical or very similar strains, whereas, for the remaining seven patients (pTHN, pZKH, pGCL, pNML, pLSY, pNFZ and pHRY), a different Candida strain or species caused subsequent infections.
Comparison of between-group and within-group genetic relatedness of isolates
Table 5 summarizes the genetic similarity between pairs of isolates for the within- and between-group comparisons. A composite UPGMA dendrogram for all 87 isolates from patients in groups A, B and C was generated, as shown in Fig. 3. Scrutiny of the dendrogram revealed that all isolates that belong to the same Candida species were grouped together within the same cluster, except for the C. famata isolates; there is a greater tendency for group B and group C isolates to cluster within their own groups than to cluster with the other groups.
|
|
Analysis of the results also indicated that the pregnancy status of a patient had no correlation with higher genetic similarity of the Candida isolates.
(i) Comparison of isolates from patients in groups A and C. For comparison of Candida isolates between patients with recurrent vaginitis (groups A and B) or single-episode vaginitis (group C), only isolates from the species C. albicans were used, so that an additional variable factor of different species would not be introduced. The SAB values of isolates from these groups were calculated; the mean SAB was found to be 0.56 ± 0.21, which represented the unrelatedness between the strains from patients in groups A and C. There were no genotypically identical pairs of strains when these two groups were compared.
For the within-group (intragroup) genetic relatedness of C. albicans isolates in the two groups, the mean SAB for group A was 0.532 ± 0.25, whereas that for group C was 0.636 ± 0.206 (Table 5). A paired t-test showed that this difference was significant (P = 0.036). Thus, it can be deduced that the C. albicans isolates in group C were more closely related to each other than to the C. albicans isolates in group A.
(ii) Comparison of isolates from patients in groups B and C. When the C. albicans isolates of patients in group B and group C were compared, two pairs of isolates were found to be genotypically identical. Isolate KTK from group B was genotypically identical to PTI and RKH from group C (SAB = 1.00). Of all possible pairs of strains that were compared, only 8.9 % of the strain-pairs from groups B and C were closely related but non-identical to each other. The rest of the strain-pairs (91.1 %) were unrelated to each other. This finding was consistent with the mean SAB of 0.55 ± 0.17, derived from pairwise comparison of group B isolates versus group C isolates, which reflected unrelatedness between strains from these two groups. This suggested that intragroup relatedness was higher than intergroup relatedness. The mean SAB value for group B, calculated from C. albicans isolates (n = 14), was 0.561 ± 0.173. This contrasts with the mean SAB value for group C, which was 0.636 ± 0.206 for C. albicans isolates only. However, this difference was not statistically significant (P = 0.325).
| DISCUSSION |
|---|
|
|
|---|
In addition, the results of this study showed that a wider spectrum of Candida species is involved in causing recurrent vaginitis than one-episode vaginitis. The lower species diversity in patients in group C, compared with those in groups A and B, was because the majority of cases of one-episode vaginitis were caused by C. albicans, except for a few isolated cases that were attributed to C. glabrata, C. lusitaniae, C. parapsilosis and C. famata. Likewise, the higher number of non-albicans species contributed to the higher values obtained for the species diversities of patients with recurrent vaginitis.
RAPD profiles of C. albicans isolates
It was interesting to note that isolates that were typed as C. albicans by using conventional methods and the API 20C AUX kit exhibited three main types of distinct RAPD profiles. Moreover, there exist at least three groups of clusters among the C. albicans isolates when all isolates from groups A, B and C were looked at, as shown in the dendrogram in Fig. 3. This could be due to several causes. Firstly, it could be that some of these isolates were not actually C. albicans strains but belonged to Candida dubliniensis, a newly delineated species that is phenotypically highly similar to C. albicans. However, this possibility could be ruled out as the API 20C AUX kit has the capacity to distinguish between C. albicans and C. dubliniensis, based on differences in growth on xylose and methyl
-D-glucose. The second possibility is that some of these isolates were of type I Candida stellatoidea, a species that shares high DNA homology with C. albicans. However, opinions still differ among taxonomists on whether or not to classify C. stellatoidea as a separate species from C. albicans. A recent finding by Biswas et al. (2001) reported the successful use of cytochrome b gene sequence analysis to differentiate between C. stellatoidea and C. albicans; hence, future work could include the use of this technique to genotype ambiguous strains. Nonetheless, there exists the third possibility that the C. albicans isolates that showed distinct genotypes were in fact due to subgroups or variants within the species, as it has been noted that a genotyping method that used RFLP analysis of rRNA genes could differentiate four genotypic groups of C. albicans (McCullough et al., 1999).
Genetic relatedness of sequential Candida isolates
Results from previous studies have led to the suggestion of three basic scenarios for the genetic relatedness of isolates obtained from patients with sequential episodes of RVC, best described by Lockhart et al. (1996). These are: (i) strain maintenance with no genetic variation; (ii) strain maintenance with minor genetic variation; and (iii) strain replacement by a different strain. Overall, it was observed from this study that all three scenarios were occurring, differing from patient to patient and from one episode to another. However, the majority of patients had recurrent infections caused by highly similar strains with minor genetic variations that can be ascribed to microevolution. This result is in agreement with the study of Lockhart et al. (1996), which reported that the main scenarios for the relatedness of sequential isolates from 18 patients were strain maintenance with and without minor variation.
From our study, an interesting observation was noted for patient pTHN, from whom five sequential isolates were obtained within a period of 1 year. The genetic similarity (measured by SAB) between the first and subsequent isolates seemed to be in descending order compared with the similarity between the first and second isolates, as more episodes of vaginitis occurred (shown by progressive decrease of computed SAB values from 0.857 to 0.8, 0.632 and finally to 0.000, for comparison of THN1 to THN2, THN3, THN4 and THN5). This seemed to suggest that a particular strain has undergone gradual genetic variation or microevolution each time it recolonized patient pTHN, probably due to drug pressure selection; by the fifth episode, it was replaced by a different species. In contrast, patient pLML seemed to exhibit strain maintenance without genetic variation, in which all three sequential isolates were genotypically identical with an SAB value of 1. For patient pGCL, the third scenario of strain replacement seemed to have occurred, as the second isolate was of a different species from the first. The remaining patients appeared to display strain maintenance with minor genetic variation.
Prevalence of strain replacement
The present study reports 14 % occurrence of strain replacement by an unrelated strain or species (4 of 28). Two patients (pHRY and pZKH) were subsequently colonized by a strain of C. albicans that was unrelated from the first episode. Patients pGCL and pTHN had subsequent vaginitis episodes caused by a totally different Candida species. Nonetheless, Lockhart et al. (1996) did not find any instances of strain replacement among the sequential isolates of 18 patients studied by using DNA fingerprinting probes. Mercure et al. (1993) found that 14 % of recurrent cases were due to colonization by a different strain, by using a Southern hybridization method. By using karyotyping resolved by PFGE, Vazquez et al. (1994) found that as many as 20 % of sequential isolates involved strain replacement.
Intragroup versus intergroup relatedness
Mean SAB values derived from statistical analysis suggested that intragroup relatedness for group C isolates was higher than groups B vs C and A vs C intergroup relatedness. In contrast, the intragroup relatedness of both group A and group B isolates did not differ significantly from the intergroup comparison with group C isolates. This suggests that strains that cause recurrent vaginitis (from groups A and B) do not possess any significant difference in their genotypes from strains that cause one-off episodes of vaginitis (group C).
From the outset, we had hoped to identify specific genotypes that correlated with recurrent vaginitis infections, based on our premise that strains that cause RVC are more virulent and would display a common genotype not found in the strains that cause occasional infections. However, we found that even strains of the same species isolated from the same patient did not share an identical phylogenetic profile. Hence, RVC in these patients could be due to a myriad of other factors, such as non-compliance with therapy, reinfection with a new strain of endogenous or exogenous origin or infection with a drug-resistant strain.
Reproducibility of RAPD analysis and comparison of results to earlier reports
To determine the reproducibility of the RAPD analysis procedure, RAPD-PCR was repeated at least three times for each DNA sample, so that the amplification profiles for each attempt could be compared. We found the profiles for a given strain to be highly reproducible, with very little variation from one RAPD analysis to another. On average, the mean SAB values reported here for each group are much lower (0.50.6) than those reported by Xu et al. (1999), which ranged from 0.6 to 0.9. However, the study by Xu et al. (1999) compared pairs of isolates obtained from different body sites within the same patients, whereas this study mainly compared pairs of vaginal isolates from different patients from either the same or different groups. In addition to using the same primer as in the current study (PA03), the study by Xu et al. (1999) also used another primer (the M13 core sequence primer); this could explain the discrepancy between the computed SAB values. Another reason could be the different algorithms used by those authors and the current study. A total of 51 polymorphic PCR fragments were obtained in the study by Xu et al. (1999), whereas the total number of polymorphic fragments obtained in the present study was 61.
Lockhart et al. (1996) employed a different method from this study, namely Southern probing with Ca3, C1 and CARE2 probes for fingerprinting. As well as the band positions of the fragments, signal intensities were also considered in computation of SAB values by Lockhart et al. (1996). Band intensities were not taken into account in this study, as PCR is known to produce fragments of varying intensities. As the PA03 primer is a random arbitrary oligonucleotide, it generated a vast range of SAB values, from as low as 0 to the maximum value of 1. In contrast, the range of SAB values obtained by Lockhart et al. (1996) was 0.691.00. When a different algorithm, the simple matching (Apostol) metric with the formula SAB = (n11+n00)/n, was used to calculate SAB values in the present study, higher values were obtained than those obtained by using the Dice metric (Dice, 1945).
Another recent study in our laboratory involved genotyping of the Candida isolates from this study by using different primer sets, namely primers ERIC1 and ERIC2 (S. R. Abdul Hadi, P. P. Chong, B. C. Tan & K. P. Ng, unpublished data). The resultant genotypic clusters and subgroups correlated well with results obtained by using the PA03 primer, although only 68 unique RAPD profiles and 39 distinct polymorphic fragments were generated, which is comparatively fewer than the number obtained with the PA03 primer. Metzgar et al. (1998) also used these ERIC primers to perform RAPD analysis of Candida strains isolated from HIV-infected patients.
In summary, results from this study demonstrated that genotyping of Candida species by using the PA03 primer and RAPD-PCR is feasible for differentiation of various species and strains. The results also strongly suggest that recurrent vaginal candidiasis patients in Malaysia are most frequently colonized by a single persistent Candida strain over an extended period of time, with minor microevolution due to accumulated genetic changes.
| Acknowledgements |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
X. Zhou, R. Westman, R. Hickey, M. A. Hansmann, C. Kennedy, T. W. Osborn, and L. J. Forney Vaginal Microbiota of Women with Frequent Vulvovaginal Candidiasis Infect. Immun., September 1, 2009; 77(9): 4130 - 4135. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. C. Odds and M. D. Jacobsen Multilocus Sequence Typing of Pathogenic Candida Species Eukaryot. Cell, July 1, 2008; 7(7): 1075 - 1084. [Full Text] [PDF] |
||||
![]() |
P. P. Chong, D. C.-S. Chieng, L. Y. Low, A. Hafeez, M. N. Shamsudin, H. F. Seow, and K. P. Ng Recurrent candidaemia in a neonate with Hirschsprung's disease: fluconazole resistance and genetic relatedness of eight Candida tropicalis isolates. J. Med. Microbiol., April 1, 2006; 55(Pt 4): 423 - 428. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. C. Baeza, M. T. Matsumoto, A. M. F. Almeida, and M. J. S. Mendes-Giannini Strain differentiation of Trichophyton rubrum by randomly amplified polymorphic DNA and analysis of rDNA nontranscribed spacer. J. Med. Microbiol., April 1, 2006; 55(Pt 4): 429 - 436. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Sampaio, L. Gusmao, A. Correia, C. Alves, A. G. Rodrigues, C. Pina-Vaz, A. Amorim, and C. Pais New Microsatellite Multiplex PCR for Candida albicans Strain Typing Reveals Microevolutionary Changes J. Clin. Microbiol., August 1, 2005; 43(8): 3869 - 3876. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |