Skip to main content


  • Methodology article
  • Open Access

Reliable differentiation of Meyerozyma guilliermondii from Meyerozyma caribbica by internal transcribed spacer restriction fingerprinting

  • 1, 2,
  • 1, 2,
  • 2 and
  • 1Email author
BMC Microbiology201414:52

  • Received: 30 July 2013
  • Accepted: 24 February 2014
  • Published:



Meyerozyma guilliermondii (anamorph Candida guilliermondii) and Meyerozyma caribbica (anamorph Candida fermentati) are closely related species of the genetically heterogenous M. guilliermondii complex. Conventional phenotypic methods frequently misidentify the species within this complex and also with other species of the Saccharomycotina CTG clade. Even the long-established sequencing of large subunit (LSU) rRNA gene remains ambiguous. We also faced similar problem during identification of yeast isolates of M. guilliermondii complex from indigenous bamboo shoot fermentation in North East India. There is a need for development of reliable and accurate identification methods for these closely related species because of their increasing importance as emerging infectious yeasts and associated biotechnological attributes.


We targeted the highly variable internal transcribed spacer (ITS) region (ITS1-5.8S-ITS2) and identified seven restriction enzymes through in silico analysis for differentiating M. guilliermondii from M. caribbica. Fifty five isolates of M. guilliermondii complex which could not be delineated into species-specific taxonomic ranks by API 20 C AUX and LSU rRNA gene D1/D2 sequencing were subjected to ITS-restriction fragment length polymorphism (ITS-RFLP) analysis. Taq I ITS-RFLP distinctly differentiated the isolates into M. guilliermondii (47 isolates) and M. caribbica (08 isolates) with reproducible species-specific patterns similar to the in silico prediction. The reliability of this method was validated by ITS1-5.8S-ITS2 sequencing, mitochondrial DNA RFLP and electrophoretic karyotyping.


We herein described a reliable ITS-RFLP method for distinct differentiation of frequently misidentified M. guilliermondii from M. caribbica. Even though in silico analysis differentiated other closely related species of M. guilliermondii complex from the above two species, it is yet to be confirmed by in vitro analysis using reference strains. This method can be used as a reliable tool for rapid and accurate identification of closely related species of M. guilliermondii complex and for differentiating emerging infectious yeasts of the Saccharomycotina CTG clade.


  • Meyerozyma guilliermondii
  • Emerging infectious yeasts
  • mtDNA-RFLP
  • PFGE karyotyping
  • Soibum


Meyerozyma guilliermondii is a genetically heterogenous complex belonging to the Saccharomycotina CTG clade [1]. This complex consists of phenotypically indistinguishable and closely related species namely Meyerozyma guilliermondii (anamorph Candida guilliermondii), Meyerozyma caribbica (anamorph Candida fermentati), Candida carpophila, Candida smithsonii, Candida athensensis, Candida elateridarum and Candida glucosophila[26]. Apart from its presence in healthy human [7, 8], M. guilliermondii also exists in clinical [3, 9] and environmental samples [10]. This organism is widely studied in various aspects due to its clinical importance, biotechnological applications and biological control potential [11].

C. guilliermondii is regarded as an emerging infectious yeast of the non-albicans Candida (NAC) species group which accounts for 1 – 5% of nosocomial blood stream infections worldwide [9, 12, 13]. However, in certain geographical regions such as Brazil, India and Italy, over 10% of all the candidaemia cases are caused by this species [14]. The threat posed by this organism is ever increasing due to the decreased susceptibility and emergence of strains resistant to antifungal drugs like polyene (amphotericin B) and azoles (fluconazole and itraconazole), leading to mortality in candidaemia patients [9, 12, 15]. C. fermentati has been rarely found to be associated with candidaemia [16, 17]. But due to the poor discernability of C. fermentati from C. guilliermondii, they are commonly misidentified in clinical laboratories. Apart from being organisms of clinical importance, M. guilliermondii and M. caribbica are often linked with fermented foods [1820]. M. guilliermondii is known for the production of flavour compounds in fermented food products [21]. Further, in a study with soybean paste fermentation, M. guilliermondii and M. caribbica have been claimed for the efficient production of isoflavone aglycone which is a widely known bioactive compound for its various health promoting functions [22]. M. guilliermondii is a flavinogenic yeast which is known for the overproduction of vitamin B2 (riboflavin) [23]. Moreover, isolates of M. guilliermondii and M. caribbica have exhibited great potential in the biological control of fungi responsible for postharvest spoilage of fruits and vegetables [2426]. These yeast species with enhanced biological control efficacy have emerged as a potential alternative to the conventional fungicide treatment. Considering the various importance and applications of the two species, there is a need for the development of accurate and reliable method to identify and distinctly discriminate the closely related species.

Current methods of yeast identification, mostly in clinical practice, are mainly based on the conventional and rapidly evolving commercial phenotypic and biochemical methods. However, such methods are often unreliable for accurate identification of closely related yeast species [13, 27]. According to recent studies, M. guilliermondii and M. caribbica are extremely difficult to differentiate by the phenotypic methods [2831]. We also faced similar problem during differentiation of yeast isolates from soibum, an indigenous fermented bamboo shoot product of North East India (Additional file 1: Table S1). The widely used API 20 C AUX yeast identification system and sequencing of large subunit (LSU) rRNA gene D1/D2 domain failed to give proper species-level taxonomic assignment to these isolates (Additional file 1: Tables S2 and S3). Moreover, the phylogenetic tree reconstructed from the publicly available D1/D2 sequences of different strains of M. guilliermondii and M. caribbica failed to discriminate the two species (Additional file 2: Figure S1).

Several attempts have been made using molecular approaches such as DNA base composition, electrophoretic karyotyping [6, 32], multi locus sequence typing (MLST) [3], multi locus enzyme electrophoresis (MLEE), randomly amplified polymorphic DNA (RAPD) [4], sequencing of internal transcribed spacer (ITS) [28, 30], intergenic spacer restriction fragment length polymorphism (IGS-RFLP) [29] and RFLP of housekeeping genes such as riboflavin synthetase gene RIBO[17] in order to resolve the misidentification. Some recent studies have claimed that the matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF-MS) is advantageous over previous approaches for reliable identification of clinically important NAC and non-Candida yeast species [28, 31, 33, 34]. Unfortunately, MALDI-TOF-MS requires reference spectra of accurately identified closely related strains otherwise the results may be erroneous. On the other hand, the sequence-based studies have considered the ITS1-5.8S-ITS2 region as universal DNA barcode for yeast identification [35] and the RFLP of ITS1-5.8S-ITS2 region has successfully separated the closely related species in the genera Candida and Pichia[36, 37]. Therefore, in this study, we targeted the ITS1-5.8S-ITS2 region to develop a simple RFLP method for accurate taxonomic assignment of M. guilliermondii and M. caribbica.

With this background, the aim of the present study was (i) to perform in silico prediction of restriction enzymes to discriminate M. guilliermondii and M. caribbica using the publicly available ITS1-5.8S-ITS2 sequences, (ii) to evaluate the selected enzymes by in vitro ITS-RFLP analysis of ambiguously identified 55 yeast isolates for species-specific taxonomic assignment, and (iii) to validate the taxonomic assignment by ITS1-5.8S-ITS2 sequencing, mitochondrial DNA (mtDNA)-RFLP and pulsed field gel electrophoresis (PFGE) karyotyping.


Yeast isolates and strains

The yeast isolates used in the present study are listed in Additional file 1: Table S1. These isolates were obtained from samples collected at different stages of indigenous bamboo shoot fermentation for the production of soibum in Manipur state of North East India [38]. The sample (10 g) was homogenized in 90 mL of sterile physiological saline (1 g/L bacteriological peptone, 8.5 g/L NaCl, pH 6.1) using Stomacher® 400 Circulator (Seward, Worthing, West Sussex) at 250 rpm for 3 min. The yeasts were isolated by serial dilution spread-plating of the above homogenate on yeast extract peptone dextrose (YEPD) agar medium (pH 6.5) (HiMedia, Mumbai, India) containing 100 μg/mL each of filter-sterilized ampicillin and tetracycline (Sigma-Aldrich, Bangalore, India), followed by incubation at 30°C for 48 − 72 h under aerobic conditions. All the isolates were purified by sub-culturing twice on the same agar medium and preserved at −80°C in YEPD broth containing 10% (v/v) sterile glycerol (Sigma-Aldrich). For short term storage, the cultures were maintained at 4°C on YEPD agar. The type strain C. guilliermondii ATCC 6260 used for comparison was obtained from American Type Culture Collection.

Phenotypic characterization and morphological observation

Phenotypic identification of the yeast isolates was carried out using the API 20 C AUX yeast identification system (bioMérieux, New Delhi, India) following manufacturer’s instructions. Colony and cell morphology of the isolates were studied using SZ-PT stereo binocular microscope (Olympus, Japan) and BX61 phase contrast microscope (Olympus).

In silico analysis and restriction enzyme selection

The full length ITS1-5.8S-ITS2 sequences of M. guilliermondii and M. caribbica were retrieved from NCBI ( and Centraalbureau voor Schimmelcultures (CBS-KNAW) yeast nucleotide databases ( Type strain sequences of the two species, C. guilliermondii ATCC 6260 [GenBank: AY939792.1] and M. caribbica CBS 9966 ( were subjected to in silico PCR amplification using primers ITS1 (5′-TCCGTAGGTGAACCTGCGG-3′) and ITS4 (5′-TCCTCCGCTTATTGATATGC-3′) [39] to trim off the untargeted regions on both 5′ and 3′ ends of the sequences using the online Sequence Manipulation Suite ( Using NEBcutter, version 2.0 (, the ITS amplicon sequences were subjected to in silico restriction digestion using the commercially available type-II restriction endonucleases listed in the REBASE database ( [40] to select the enzymes which cut the two sequences differently at not more than 5 cleavage sites. Multiple sequence alignment of 10 additional ITS1-5.8S-ITS2 sequences of different strains from different ecological niches for each species was performed using Clustal X, version 2.0 ( and BioEdit, version 7.2.0 ( to confirm the taxa-specificity of the selected restriction enzymes.

DNA extraction

DNA was extracted from pure cultures as cell-free DNA lysate using lyticase-heat lysis method. Briefly, a single colony of 24 − 48 h old culture from YEPD agar was inoculated to 5 mL of YEPD broth supplemented with antibiotics, and incubated for 18 h at 30°C with shaking at 200 rpm. Cells were harvested from 1 mL of the culture broth at 5,000 g for 5 min (FA-45-24-11, Centrifuge 5424, Eppendorf, Hamburg, Germany). The cell pellet was washed twice with 1 mL sterile 0.5 M NaCl followed by sterile deionized water (Milli Q, Millipore, Molsheim, France). The cells were finally resuspended in 500 μL of 1× TE buffer (10 mM Tris-Cl, 1 mM EDTA, pH 8.0) containing 10 μL of lyticase (5U/μL) (Sigma-Aldrich) and incubated at 37°C for 1 h. After the incubation, the spheroplasts were lysed by heating at 95°C for 20 min. The crude cell-free lysate was collected by centrifugation at 10,000 g for 10 min at 4°C and the DNA was quantified spectrophotometrically (Nanodrop ND-1000, NanoDrop Technologies, Inc., Rockland, USA). The cell-free lysate with absorbance ratio (A260/280) of 1.8 − 2.2 was used for PCR analysis and stored at −20°C until required.


ITS1-5.8S-ITS2 was amplified from the cell-free DNA lysate using primers ITS1 and ITS4 mentioned elsewhere. The amplification was carried out in a 25 μL final reaction volume containing 50 ng of the genomic DNA as previously described [41]. The amplified ITS fragment was analyzed by 2.0% (w/v) agarose gel electrophoresis at 80 V in 0.5× TBE (45 mM Tris-borate, 1 mM EDTA, pH 8.0) buffer to check its intactness and absence of non-specific amplification. The PCR product (4 μL) was digested with 5 U of Taq I (Promega, Madison, USA) in a 10 μL reaction volume at 65°C as per manufacturer’s instructions. The restriction patterns were analyzed by electrophoresis of the 10 μL reaction volume on 2.0% (w/v) agarose gel in parallel with PCR 100 bp Low DNA ladder (Sigma-Aldrich) as molecular size standard. The electrophoresis was run at 80 V for 2 h in 0.5× TBE buffer. The gel was then stained in 0.5 μg/mL ethidium bromide solution for 30 min with rocking at 15 rpm on a platform rocker (Tarsons, Kolkata, India). After destaining for 30 min in sterile deionized water, the gel was photographed using ChemiDoc MP gel documentation system (Bio Rad, Hercules, USA). The restriction fingerprints were analysed for the absence or presence of discriminating fragments using GelCompar II software, version 6.5 (Applied Maths, Sint-Martens-Latem, Belgium).


A single colony of 24 − 48 h old culture from YEPD agar was inoculated to 5 mL of YEPD broth supplemented with antibiotics, and incubated for 18 h at 30°C with shaking at 200 rpm. The grown culture was inoculated into 50 mL of fresh YEPD broth (initial OD600 =0.1) and incubated in the above conditions till mid-logarithmic growth phase (final OD600=0.4 − 0.8). Cells of 20 OD600 were harvested at 1,800 g for 5 min at 4°C (A-4-81, Centrifuge 5810R, Eppendorf). The mtDNA was extracted as previously described [42] with some modifications. The cells were resuspended and washed with 5 mL of yeast resuspension buffer (50 mM Tris-Cl, 20 mM EDTA, pH 8.0) and stored at −20°C for 10 min. Lyticase (50 U) (Sigma-Aldrich) was used to produce spheroplast and 15 μL of 1 mg/mL RNase A solution (Sigma-Aldrich) was added during cell lysis. The total DNA was precipitated at −20°C for 1 h. After quantifying the DNA content spectrophotometrically, the DNA was freeze dried, re-dissolved in sterile deionized water to a final concentration of 1 μg/μL and stored at −20°C till further use. Restriction digestion was carried out on 10 μg of the DNA in a 20 μL reaction volume using 10 U each of Hae III and Hinf I (Promega) according to manufacturer’s instructions. The restriction patterns were generated by 1.0% (w/v) agarose gel electrophoresis of the 20 μL reaction volume at 80 V in 0.5× TBE buffer for 4 h in parallel with 1 kb DNA ladder (Promega). After staining and documentation, the restriction fingerprints were subjected to cluster analysis using unweighted pair group method with arithmetic mean (UPGMA) algorithm on Jaccard similarity coefficients using GelCompar II. Composite data set of the restriction digestion profiles was generated with 1.0% position tolerance to generate the clustering. Bootstrap analysis with 1,000 replicates was performed to indicate the branch quality.

Electrophoretic karyotyping

Intact chromosomal DNA for electrophoretic karyotyping using PFGE was prepared as previously described [32]. The electrophoresis was carried out in 1.0% (w/v) PFGE-grade agarose gel (Sigma-Aldrich) and 0.5× TBE buffer at 13 − 14°C and 150 V in contour-clamped homogeneous electric field electrophoresis apparatus (Gene Navigator, Amersham Biosciences, Uppsala, Sweden). The gel was run for 22 h with a switch interval of 90 s for 8 h followed by 105 s for 6 h and finally 120 s for 8 h in parallel with PFGE marker (225 − 22,000 kb) from Saccharomyces cerevisiae strain YPH80 (Sigma-Aldrich). Staining and documentation were performed as mentioned elsewhere.

ITS and D1/D2 sequencing and sequence analysis

The representative isolates from each ITS-RFLP genotype group were randomly selected for sequencing ITS1-5.8S-ITS2 and LSU rRNA gene D1/D2 domain. ITS1-5.8S-ITS2 was PCR amplified as mentioned elsewhere. Amplification of D1/D2 region was carried out using primers NL1 (5′-GCATATCAATAAGCGGAGGAAAAG-3′) and NL4 (5′-GGTCCGTGTTTCAAGACGG-3′) as previously described [43]. The amplified products were purified using NucleoSpin® Extract II gel extraction kit (Machery-Nagel, Düren, Germany) following manufacturer’s instructions. The PCR products were sequenced using ABI 3100 Genetic Analyser (Merck, Bangalore, India) with the same primers used for the amplification. The sequence reads were validated by analysing the electropherogram data using ChromasLITE software, version 2.01 ( To identify the closest known relatives, the sequences were queried with NCBI and CBS yeast nucleotide databases. The sequences obtained from both sequencing and nucleotide databases were aligned using Clustal X algorithm and a neighbour joining tree was constructed by Kimura’s evolutionary distance matrix obtained from the multiple sequence alignment using MEGA4 phylogenetic software. Bootstrap values for 1000 replicates were shown at the node of cluster branch. The sequences were deposited to NCBI GenBank under the following accession numbers: JF439366 − JF439369 and KF268351 − KF268354.


In silico selection of differentiating restriction enzymes

The full length ITS1-5.8S-ITS2 sequences of M. guilliermondii and M. caribbica were retrieved from NCBI and CBS yeast nucleotide databases and subjected to multiple sequence alignment followed by in silico restriction digestion. Three variable regions differentiating M. guilliermondii from M. caribbica were identified. Seven restriction enzymes (Ars I, Bfa I, Bsr I, Hpy 188I, Hpy CH4III, Mme I and Taq I) which cut the variable regions differently were identified (Figure 1A and Additional file 2: Figure S2). Considering the length and the number of polymorphic fragments with sizes greater than 100 bp (for easy analysis in normal agarose gel), Bfa I, Mme I and Taq I were found appropriate. Notably, commonly available Taq I gave distinct species-specific differentiation between the two species (Additional file 1: Table S4 and Additional file 2: Figure S3). We also tested the selected three restriction enzymes (Bfa I, Mme I and Taq I) for differentiating M. guilliermondii and M. caribbica from other closely related members of M. guilliermondii complex (Additional file 1: Table S4). Except C. carpophila and M. caribbica, all other members of M. guilliermondii complex were distinctly differentiated during the analysis.
Figure 1
Figure 1

Differentiation of M. guilliermondii and M. caribbica by Taq I digestion of ITS1-5.8S-ITS2. A: Multiple sequence alignment of representative ITS1-5.8S-ITS2 sequences of the two species obtained from NCBI GenBank and CBS yeast database showing position of Taq I recognition site (highlighted) which distinctly differentiated the two species. B: Taq I restriction digestion profile of ITS1-5.8S-ITS2 amplicons obtained from some of the representative isolates. Lane 1: C. guilliermondii ATCC 6260; Lane 2 − 12: isolates of M. guilliermondii genotype group MG (A1S10Y1, A2S10Y1, A3S9Y1, A2S9Y1, A3S11Y1, A3S2Y1, A3S6Y1, A2S6Y1, A1S9Y1, Kw3S3Y1 and Kw2S11Y2); Lane 13 – 20: isolates of M. caribbica genotype group MC (A1S10Y2a, A1S10Y3, A1S10Y5, Kw3S2Y1, Kw2S3Y1, Kw3S3Y3, Kw3S3Y4 and Kw1S7Y2); Lane M: PCR 100 bp Low DNA ladder (Sigma-Aldrich).

Evaluation of in silico selected restriction enzymes by in vitro ITS-RFLP

To validate the above in silico selection, the 55 yeast isolates of M. guilliermondii complex (which were not differentiated by phenotypic characterization and D1/D2 sequencing) were analysed by ITS-RFLP using the selected Taq I restriction enzyme in comparison with the type strain C. guilliermondii ATCC 6260. All the tested isolates and the type strain gave a single PCR amplicon of molecular size of 607 bp. As predicted by the in silico analysis, Taq I ITS-RFLP distinctly differentiated the isolates into two genotype groups. Forty seven isolates produced M. guilliermondii-specific pattern (MG), while the remaining eight isolates generated M. caribbica-specific pattern (MC) (Table 1). Examples of Taq I ITS-RFLP profiles differentiating the above two species are shown in Figure 1B.
Table 1

Differentiation of ambiguous 55 yeast isolates obtained from soibum into Meyerozyma guilliermondii and Meyerozyma caribbica

Group (Number of isolates)

Representative strains

Taxonomic designation








MG (47)

A1S10Y1, Kw2S11Y2

M. guilliermondii

M. guilliermondii

M. guilliermondii/M. caribbica (JF439368, JF439369)

M. guilliermondii (KF268351, KF268352)

M. guilliermondii

M. guilliermondii

MC (08)

Kw1S7Y2, Kw3S2Y1

M. guilliermondii

M. caribbica

M. guilliermondii/M. caribbica (JF439366, JF439367)

M. caribbica (KF268353, KF268354)

M. caribbica

M. caribbica

Type strain

ATCC 6260

M. guilliermondii

M. guilliermondii

M. guilliermondii/M. caribbica (AJ508562.1)

M. guilliermondii (AY939792.1)

M. guilliermondii

M. guilliermondii

*No identification data for M. caribbica is included in the database.

GenBank accession numbers.

Validation of ITS-RFLP assay

The above ITS-RFLP based discrimination of M. guilliermondii and M. caribbica was further confirmed by ITS1-5.8S-ITS2 sequencing, mtDNA-RFLP fingerprinting and PFGE karyotyping (Table 1). The ITS sequences of the isolates in each genotype group MG and MC matched with the sequences of the type strains C. guilliermondii ATCC 6260 and M. caribbica CBS 9966 with 99.6% and 99.8% similarity respectively. The sequences between the two groups were 99% identical showing only 5 nucleotide differences which were the same as shown by the above type strain sequences. Unlike D1/D2 region, the ITS sequences formed distinct cluster of M. guilliermondii and M. caribbica during phylogenetic analysis (Figure 2). The ITS sequences of M. guilliermondii strains PX-PAT, CanR-56 and SD 337; M. caribbica strains UCLM 44A and IWBT-Y836 showed 1 – 3 nucleotide differences which was reflected as divergence within the clusters. The mean evolutionary divergence of 0.0131 between the two clusters was 6 times more than the divergence within each cluster.
Figure 2
Figure 2

Neighbour-joining (NJ) phylogenetic tree showing taxa-specific separation of M. guilliermondii from M. caribbica. The tree was constructed based on the evolutionary distance calculated using Kimura-2 parameter from the nucleotide sequence of ITS1-5.8S-ITS2. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches for values >40%. The bar represents 1% sequence divergence. GenBank accession numbers are mentioned within the parentheses. S. cerevisiae was the outgroup in the analysis. T =Type strain.

The mtDNA-RFLP using Hae III and Hinf I distinctly segregated the yeast isolates into M. guilliermondii and M. caribbica. mtDNA-RFLP profile-based dendrogram formed two clusters (Figure 3) similar to the ITS-RFLP groups. Between the two enzymes used, Hinf I showed higher polymorphism than Hae III. Electrophoretic karyotyping also distinctly discriminated the above two species (Figure 4). The species-specific mtDNA-RFLP pattern suggested that the isolates of each group belonged to only one strain (Figure 3). Whereas electrophoretic karyotyping brought out strain level diversity in both the groups which confirmed that multiple strains of M. guilliermondii and M. caribbica were involved in the indigenous bamboo shoot fermentation (Figure 4 and Additional file 2: Figure S4).
Figure 3
Figure 3

mtDNA-RFLP based dendrogram showing distinct clustering of M. guilliermondii and M. caribbica . The dendrogram was constructed using UPGMA algorithm on Jaccard similarity coefficients generated from Hae III and Hinf I restriction digestion profile of mtDNA of some of the representative isolates. Value at each branch node indicates the branch quality with 1000 bootstrap replications. The scale represents the similarity.

Figure 4
Figure 4

PFGE karyotype patterns of isolates belonging to M. guilliermondii and M. caribbica genotype groups. Lane 1: C. guilliermondii ATCC 6260; Lane 2 − 3: M. guilliermondii isolates A1S10Y1 and Kw2S11Y2; Lane 4 − 11: M. caribbica isolates A1S10Y2a, A1S10Y3, A1S10Y5, Kw3S2Y1, Kw2S3Y1, Kw3S3Y3, Kw3S3Y4 and Kw1S7Y2; Lane M: S. cerevisiae PFGE marker (Sigma-Aldrich). Right arrow indicates the co-migrating chromosomal doublets showing strain level diversity.


In recent times, the frequency of emerging infectious diseases caused by the opportunistic yeast species of NAC and non-Candida groups has increased in immunosuppressed patients [12, 44]. This is linked with the indiscriminate use of broad-spectrum antifungal drugs and global climate change [4547]. Most of these closely related yeast species are often misidentified by the conventional phenotypic, biochemical and antibiotic susceptibility methods. Thus, accurate identification of these species has become essential for clinical management and effective antifungal therapy.

In the present study, we showed that the developed ITS-RFLP method was reliable and consistent for distinct differentiation of closely related M. guilliermondii from M. caribbica for which phenotypic methods and D1/D2 sequencing were inconclusive. Our results also indicated that sequencing of both D1/D2 and ITS regions will increase the resolution of species identification which can be further improved by multigene sequence-based phylogenetic approach [3, 48, 49]. However, the presence of incorrectly identified, insufficiently annotated and non-updated entries in the public nucleotide databases may underestimate the resolving power of these taxonomic markers [50]. Out of the 29 sequences of LSU rRNA gene for M. guilliermondii available in NCBI GenBank, 17 sequences (58.62%) clustered with M. caribbica type strain CBS 9966 [GenBank: EU348786] (Additional file 2: Figure S1).

The choice of appropriate restriction endonucleases is critical for RFLP experiments. The commonly used Cfo I, Hae III and Hinf I [37, 41] failed to segregate M. guilliermondii from other species of the same genus during in silico and in vitro ITS-RFLP analysis. Our results indicated that in silico selection of restriction enzymes using the publicly available sequences from various strains of the target species is a better approach than randomly selecting the previously described and commonly used enzymes. This approach has been proven to be highly effective and reproducible [36, 5153], and many online resources have been developed for this purpose [5457]. Clinical isolates of Candida famata and Candida palmioleophila were also frequently misidentified as M. guilliermondii[30, 31]. In silico analysis confirmed that the developed ITS-RFLP method can also discriminate these species (data not shown). This in silico selection approach can be effectively applied to other cryptic yeast species of clinical importance for the development of RFLP based diagnostic tools.

The developed method of ITS-RFLP using Taq I differentiated M. guilliermondii and M. caribbica at species level. This method is simple, rapid and reliable in comparison to the commonly used sequencing methods. The entire analysis starting from DNA extraction to ITS-RFLP profiling can be completed within 8 h. Further studies using higher number of strains of these two species from different clinical sources are required to confirm the robustness of this method for diagnostic applications. Though the combination of ITS-RFLP profiles generated by Taq I, Bfa I and Mme I differentiated other closely related species of the M. guilliermondii complex from M. guilliermondii and M. caribbica during in silico analysis, it is yet to be confirmed through in vitro analysis using reference strains.

In our present study the use of either Hae III or Hinf I for mtDNA-RFLP was sufficient to differentiate the two species, but failed to give strain level differentiation even when both enzymes were used. This demonstrated that mtDNA-RFLP can also be used for distinct differentiation of closely related species. The Hinf I mtDNA-RFLP pattern of our M. guilliermondii isolates was similar with the mtDNA restriction pattern ‘E’ of M. guilliermondii strains isolated from wineries in Alentejo, Portugal [58]. This genotype was linked with the production of flavour compound, 4-ethylphenol in wine. The major phenolic flavour compound (4-methylphenol) detected from fermented bamboo shoot product, soibum (Singh NR: unpublished observations) might also have originated from M. guilliermondii. Future study is required to characterize the flavour compound producing strain for starter culture development. Though fresh bamboo shoots are highly perishable, the fermented bamboo shoot can be preserved up to one year after fermentation without any deterioration or change in its organoleptic character. This long term preservation may be linked with the dominant presence of M. guilliermondii which has been reported as an efficient biological control agent [24, 25]. Being an emerging infectious yeast, the presence of M. guilliermondii in fermented food is a great concern regarding the safety of its consumption. Further study in strain level is required to unravel the pathogenic potential of M. guilliermondii associated with soibum fermentation.


In this study, we described an ITS-RFLP method developed through an integrated approach of in silico selection of restriction enzymes and in vitro validations for distinct differentiation of frequently misidentified M. guilliermondii from M. caribbica, which can be used as an alternative or an adjunct to ITS sequencing. This method may be used for rapid and accurate identification of emerging infectious yeasts of the Saccharomycotina CTG clade. This approach can also be used for other closely related species complex when phenotypic methods and D1/D2 sequencing are ambiguous.



The research was supported by the Department of Biotechnology (DBT), Govt. of India funded project (BT/PR-9268/FNS/20/342/2007). Wahengbam Romi is a recipient of Senior Research Fellowship from Council of Scientific and Industrial Research (CSIR), Govt. of India (112417/2 K10/1). We are grateful to Prof. N. Rajmuhon Singh for providing the data of flavour compounds associated with soibum. The authors would like to thank the indigenous producers of soibum in Andro and Kwatha villages, Manipur, India for their support during sample collection.

Authors’ Affiliations

Microbial Resources Division, Institute of Bioresources and Sustainable Development (IBSD), Takyelpat Institutional Area, Imphal, 795001, Manipur, India
Department of Biotechnology, Gauhati University, Guwahati, 781014, Assam, India


  1. Dujon B: Yeast evolutionary genomics. Nat Rev Genet. 2010, 11: 512-524. 10.1038/nrg2811.View ArticlePubMedGoogle Scholar
  2. Lachance MA, Boekhout T, Scorzetti G, Fell JW, Kurtzman CP: Candida Berkhout (1923). The Yeasts: A Taxonomic Study, Volume 2. Edited by: Kurtzman CP, Fell JW, Boekhout T. 2011, San Diego: Elsevier, 987-1278. 5View ArticleGoogle Scholar
  3. Lan L, Xu J: Multiple gene genealogical analyses suggest divergence and recent clonal dispersal in the opportunistic human pathogen Candida guilliermondii. Microbiology. 2006, 152: 1539-1549. 10.1099/mic.0.28626-0.View ArticlePubMedGoogle Scholar
  4. San Millan RM, Wu LC, Salkin IF, Lehmann PF: Clinical isolates of Candida guilliermondii include Candida fermentati. Int J Syst Bacteriol. 1997, 47: 385-393. 10.1099/00207713-47-2-385.View ArticlePubMedGoogle Scholar
  5. Suh SO, Blackwell M: Three new beetle-associated yeast species in the Pichia guilliermondii clade. FEMS Yeast Res. 2004, 5: 87-95. 10.1016/j.femsyr.2004.06.001.View ArticlePubMedGoogle Scholar
  6. Vaughan-Martini A, Kurtzman CP, Meyer SA, O’Neill EB: Two new species in the Pichia guilliermondii clade: Pichia caribbica sp. nov., the ascosporic state of Candida fermentati, and Candida carpophila comb. nov. FEMS Yeast Res. 2005, 5: 463-469. 10.1016/j.femsyr.2004.10.008.View ArticlePubMedGoogle Scholar
  7. Kam AP, Xu J: Diversity of commensal yeasts within and among healthy hosts. Diagn Microbiol Infect Dis. 2002, 43: 19-28. 10.1016/S0732-8893(02)00364-4.View ArticlePubMedGoogle Scholar
  8. Xu J, Mitchell TG: Geographical differences in human oral yeast flora. Clin Infect Dis. 2003, 36: 221-224. 10.1086/345672.View ArticlePubMedGoogle Scholar
  9. Krcmery V, Barnes AJ: Non-albicans Candida spp. causing fungaemia: pathogenicity and antifungal resistance. J Hosp Infect. 2002, 50: 243-260. 10.1053/jhin.2001.1151.View ArticlePubMedGoogle Scholar
  10. Savini V, Catavitello C, Onofrillo D, Masciarelli G, Astolfi D, Balbinot A, Febbo F, D’Amario C, D’Antonio D: What do we know about Candida guilliermondii? A voyage throughout past and current literature about this emerging yeast. Mycoses. 2011, 54: 434-441. 10.1111/j.1439-0507.2010.01960.x.View ArticlePubMedGoogle Scholar
  11. Papon N, Savini V, Lanoue A, Simkin AJ, Creche J, Giglioli-Guivarc’h N, Clastre M, Courdavault V, Sibirny AA: Candida guilliermondii: biotechnological applications, perspectives for biological control, emerging clinical importance and recent advances in genetics. Curr Genet. 2013, (in press) (doi:10.1007/s00294-013-0391-0)Google Scholar
  12. Miceli MH, Diaz JA, Lee SA: Emerging opportunistic yeast infections. Lancet Infect Dis. 2011, 11: 142-151. 10.1016/S1473-3099(10)70218-8.View ArticlePubMedGoogle Scholar
  13. Neppelenbroek K, Seo R, Urban V, Silva S, Dovigo L, Jorge J, Campanha N: Identification of Candida species in the clinical laboratory: a review of conventional, commercial, and molecular techniques. Oral Dis. 2013, (in press) (doi:10.1111/odi.12123)Google Scholar
  14. Sandven P: Epidemiology of candidemia. Rev Iberoam Micol. 2000, 17: 73-81.PubMedGoogle Scholar
  15. Pfaller MA, Diekema DJ, Gibbs DL, Newell VA, Ellis D, Tullio V, Rodloff A, Fu W, Ling TA: Results from the ARTEMIS DISK Global Antifungal Surveillance Study, 1997 to 2007: a 10.5-year analysis of susceptibilities of Candida species to fluconazole and voriconazole as determined by CLSI standardized disk diffusion. J Clin Microbiol. 2010, 48: 1366-1377. 10.1128/JCM.02117-09.PubMed CentralView ArticlePubMedGoogle Scholar
  16. Chen CY, Huang SY, Tang JL, Tsay W, Yao M, Ko BS, Chou WC, Tien HF, Hsueh PR: Clinical features of patients with infections caused by Candida guilliermondii and Candida fermentati and antifungal susceptibility of the isolates at a medical centre in Taiwan, 2001-10. J Antimicrob Chemother. 2013, (in press) (doi:10.1093/jac/dkt214)Google Scholar
  17. Lockhart SR, Messer SA, Pfaller MA, Diekema DJ: Identification and susceptibility profile of Candida fermentati from a worldwide collection of Candida guilliermondii clinical isolates. J Clin Microbiol. 2009, 47: 242-244. 10.1128/JCM.01889-08.PubMed CentralView ArticlePubMedGoogle Scholar
  18. Daniel HM, Vrancken G, Takrama JF, Camu N, de Vos P, de Vuyst L: Yeast diversity of Ghanaian cocoa bean heap fermentations. FEMS Yeast Res. 2009, 9: 774-783. 10.1111/j.1567-1364.2009.00520.x.View ArticlePubMedGoogle Scholar
  19. Hidalgo C, Mateo E, Mas A, Torija MJ: Identification of yeast and acetic acid bacteria isolated from the fermentation and acetification of persimmon (Diospyros kaki). Food Microbiol. 2012, 30: 98-104. 10.1016/ ArticlePubMedGoogle Scholar
  20. Nova MX, Schuler AR, Brasileiro BT, Morais MA: Yeast species involved in artisanal cachaca fermentation in three stills with different technological levels in Pernambuco, Brazil. Food Microbiol. 2009, 26: 460-466. 10.1016/ ArticlePubMedGoogle Scholar
  21. Wah TT, Walaisri S, Assavanig A, Niamsiri N, Lertsiri S: Co-culturing of Pichia guilliermondii enhanced volatile flavor compound formation by Zygosaccharomyces rouxii in the model system of Thai soy sauce fermentation. Int J Food Microbiol. 2013, 160: 282-289. 10.1016/j.ijfoodmicro.2012.10.022.View ArticlePubMedGoogle Scholar
  22. Kim WC, So JH, Kim SI, Shin JH, Song KS, Yu CB, Kho YH, Rhee IK: Isolation, identification, and characterization of Pichia guilliermondii K123-1 and Candida fermentati SI, producing isoflavone β-glycosidase to hydrolyze isoflavone glycoside efficiently, from the korean traditional soybean paste. J Appl Biol Chem. 2009, 52: 163-169. 10.3839/jabc.2009.028.View ArticleGoogle Scholar
  23. Boretsky YR, Pynyaha YV, Boretsky VY, Fedorovych DV, Fayura LR, Protchenko O, Philpott CC, Sibirny AA: Identification of the genes affecting the regulation of riboflavin synthesis in the flavinogenic yeast Pichia guilliermondii using insertion mutagenesis. FEMS Yeast Res. 2011, 11: 307-314. 10.1111/j.1567-1364.2011.00720.x.PubMed CentralView ArticlePubMedGoogle Scholar
  24. Wilson CL, Chalutz E: Pichia guilliermondii (anamorph Candida guilliermondii) useful for the biological control of postharvest rots in fruits. US Patent. 1991, 8: 5041384-Google Scholar
  25. Lima JR, Gondime DMF, Oliveirab JTA, Oliveirad FSA, Gonçalvesc LRB, Viana FMP: Use of killer yeast in the management of postharvest papaya anthracnose. Postharvest Biol Technol. 2013, 83: 58-64.View ArticleGoogle Scholar
  26. Cao J, Zhang H, Yang Q, Ren R: Efficacy of Pichia caribbica in controlling blue mold rot and patulin degradation in apples. Int J Food Microbiol. 2013, 162: 167-173. 10.1016/j.ijfoodmicro.2013.01.007.View ArticlePubMedGoogle Scholar
  27. Pincus DH, Orenga S, Chatellier S: Yeast identification−past, present, and future methods. Med Mycol. 2007, 45: 97-121. 10.1080/13693780601059936.View ArticlePubMedGoogle Scholar
  28. Castanheira M, Woosley LN, Diekema DJ, Jones RN, Pfaller MA: Candida guilliermondii and other species of Candida misidentified as Candida famata: assessment by Vitek 2, DNA sequencing analysis, and matrix-assisted laser desorption ionization-time of flight mass spectrometry in two global antifungal surveillance programs. J Clin Microbiol. 2013, 51: 117-124. 10.1128/JCM.01686-12.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Cornet M, Sendid B, Fradin C, Gaillardin C, Poulain D, Nguyen HV: Molecular identification of closely related Candida species using two ribosomal intergenic spacer fingerprinting methods. J Mol Diagn. 2011, 13: 12-22. 10.1016/j.jmoldx.2010.11.014.PubMed CentralView ArticlePubMedGoogle Scholar
  30. Desnos-Ollivier M, Ragon M, Robert V, Raoux D, Gantier JC, Dromer F: Debaryomyces hansenii (Candida famata), a rare human fungal pathogen often misidentified as Pichia guilliermondii (Candida guilliermondii). J Clin Microbiol. 2008, 46: 3237-3242. 10.1128/JCM.01451-08.PubMed CentralView ArticlePubMedGoogle Scholar
  31. Jensen RH, Arendrup MC: Candida palmioleophila: characterization of a previously overlooked pathogen and its unique susceptibility profile in comparison with five related species. J Clin Microbiol. 2011, 49: 549-556. 10.1128/JCM.02071-10.PubMed CentralView ArticlePubMedGoogle Scholar
  32. Bai FY, Liang HY, Jia JH: Taxonomic relationships among the taxa in the Candida guilliermondii complex, as revealed by comparative electrophoretic karyotyping. Int J Syst Evol Microbiol. 2000, 50: 417-422. 10.1099/00207713-50-1-417.View ArticlePubMedGoogle Scholar
  33. Marklein G, Josten M, Klanke U, Muller E, Horre R, Maier T, Wenzel T, Kostrzewa M, Bierbaum G, Hoerauf A, Sahl HG: Matrix-assisted laser desorption ionization-time of flight mass spectrometry for fast and reliable identification of clinical yeast isolates. J Clin Microbiol. 2009, 47: 2912-2917. 10.1128/JCM.00389-09.PubMed CentralView ArticlePubMedGoogle Scholar
  34. Spanu T, Posteraro B, Fiori B, D’Inzeo T, Campoli S, Ruggeri A, Tumbarello M, Canu G, Trecarichi EM, Parisi G, Tronci M, Sanguinetti M, Fadda G: Direct MALDI-TOF mass spectrometry assay of blood culture broths for rapid identification of Candida species causing bloodstream infections: an observational study in two large microbiology laboratories. J Clin Microbiol. 2012, 50: 176-179. 10.1128/JCM.05742-11.PubMed CentralView ArticlePubMedGoogle Scholar
  35. Schoch CL, Seifert KA, Huhndorf S, Robert V, Spouge JL, Levesque CA, Chen W: Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for fungi. Proc Natl Acad Sci U S A. 2012, 109: 6241-6246. 10.1073/pnas.1117018109.PubMed CentralView ArticlePubMedGoogle Scholar
  36. Trost A, Graf B, Eucker J, Sezer O, Possinger K, Gobel UB, Adam T: Identification of clinically relevant yeasts by PCR/RFLP. J Microbiol Methods. 2004, 56: 201-211. 10.1016/j.mimet.2003.10.007.View ArticlePubMedGoogle Scholar
  37. Villa-Carvajal M, Querol A, Belloch C: Identification of species in the genus Pichia by restriction of the internal transcribed spacers (ITS1 and ITS2) and the 5.8S ribosomal DNA gene. Antonie Van Leeuwenhoek. 2006, 90: 171-181. 10.1007/s10482-006-9071-0.View ArticlePubMedGoogle Scholar
  38. Jeyaram K, Singh TA, Romi W, Devi AR, Singh WM, Dayanidhi H, Singh NR, Tamang JP: Traditional fermented foods of Manipur. Indian J Tradit Knowl. 2009, 8: 115-121.Google Scholar
  39. White TJ, Bruns T, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protocols: A Guide to Methods and Applications. Edited by: Innis MA, Gelfand DH, Sninsky JJ, White TJ. 1990, New York: Academic Press Inc, 315-322.Google Scholar
  40. Roberts RJ, Vincze T, Posfai J, Macelis D: REBASE−a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 2010, 38: 234-236. 10.1093/nar/gkp874.View ArticleGoogle Scholar
  41. Jeyaram K, Singh WM, Capece A, Romano P: Molecular identification of yeast species associated with ‘Hamei’−a traditional starter used for rice wine production in Manipur, India. Int J Food Microbiol. 2008, 124: 115-125. 10.1016/j.ijfoodmicro.2008.02.029.View ArticlePubMedGoogle Scholar
  42. Querol A, Barrio E, Huerta T, Ramon D: Molecular monitoring of wine fermentations conducted by active dry yeast strains. Appl Environ Microbiol. 1992, 58: 2948-2953.PubMed CentralPubMedGoogle Scholar
  43. Kurtzman CP, Robnett CJ: Identification and phylogeny of ascomycetous yeasts from analysis of nuclear large subunit (26S) ribosomal DNA partial sequences. Antonie Van Leeuwenhoek. 1998, 73: 331-371. 10.1023/A:1001761008817.View ArticlePubMedGoogle Scholar
  44. Roden MM, Zaoutis TE, Buchanan WL, Knudsen TA, Sarkisova TA, Schaufele RL, Sein M, Sein T, Chiou CC, Chu JH, Kontoyiannis DP, Walsh TJ: Epidemiology and outcome of zygomycosis: a review of 929 reported cases. Clin Infect Dis. 2005, 41: 634-653. 10.1086/432579.View ArticlePubMedGoogle Scholar
  45. Birrenbach T, Bertschy S, Aebersold F, Mueller NJ, Achermann Y, Muehlethaler K, Zimmerli S: Emergence of Blastoschizomyces capitatus yeast infections, Central Europe. Emerg Infect Dis. 2012, 18: 98-101. 10.3201/eid1801.111192.PubMed CentralView ArticlePubMedGoogle Scholar
  46. Garcia-Solache MA, Casadevall A: Global warming will bring new fungal diseases for mammals. mBio. 2010, 1: e00061-10.PubMed CentralView ArticlePubMedGoogle Scholar
  47. Raffa RB, Eltoukhy NS, Raffa KF: Implications of climate change (global warming) for the healthcare system. J Clin Pharm Ther. 2012, 37: 502-504. 10.1111/j.1365-2710.2012.01355.x.View ArticlePubMedGoogle Scholar
  48. Tavanti A, Davidson AD, Gow NA, Maiden MC, Odds FC: Candida orthopsilosis and Candida metapsilosis spp. nov. to replace Candida parapsilosis groups II and III. J Clin Microbiol. 2005, 43: 284-292. 10.1128/JCM.43.1.284-292.2005.PubMed CentralView ArticlePubMedGoogle Scholar
  49. Tsui CK, Daniel HM, Robert V, Meyer W: Re-examining the phylogeny of clinically relevant Candida species and allied genera based on multigene analyses. FEMS Yeast Res. 2008, 8: 651-659. 10.1111/j.1567-1364.2007.00342.x.View ArticlePubMedGoogle Scholar
  50. Nilsson RH, Ryberg M, Kristiansson E, Abarenkov K, Larsson KH, Koljalg U: Taxonomic reliability of DNA sequences in public sequence databases: a fungal perspective. PLoS One. 2006, 1: e59-10.1371/journal.pone.0000059.PubMed CentralView ArticlePubMedGoogle Scholar
  51. Brugger SD, Frei L, Frey PM, Aebi S, Muhlemann K, Hilty M: 16S rRNA terminal restriction fragment length polymorphism for the characterization of the nasopharyngeal microbiota. PLoS One. 2012, 7: e52241-10.1371/journal.pone.0052241.PubMed CentralView ArticlePubMedGoogle Scholar
  52. Jeyaram K, Romi W, Singh TA, Adewumi GA, Basanti K, Oguntoyinbo FA: Distinct differentiation of closely related species of Bacillus subtilis group with industrial importance. J Microbiol Methods. 2011, 87: 161-164. 10.1016/j.mimet.2011.08.011.View ArticlePubMedGoogle Scholar
  53. Mirhendi H, Bruun B, Schonheyder HC, Christensen JJ, Fuursted K, Gahrn-Hansen B, Johansen HK, Nielsen L, Knudsen JD, Arendrup MC: Molecular screening for Candida orthopsilosis and Candida metapsilosis among Danish Candida parapsilosis group blood culture isolates: proposal of a new RFLP profile for differentiation. J Med Microbiol. 2010, 59: 414-420. 10.1099/jmm.0.017293-0.View ArticlePubMedGoogle Scholar
  54. Bikandi J, San Millan R, Rementeria A, Garaizar J: In silico analysis of complete bacterial genomes: PCR, AFLP-PCR and endonuclease restriction. Bioinformatics. 2004, 20: 798-799. 10.1093/bioinformatics/btg491.View ArticlePubMedGoogle Scholar
  55. Collins RE, Rocap G: REPK: an analytical web server to select restriction endonucleases for terminal restriction fragment length polymorphism analysis. Nucleic Acids Res. 2007, 35: W58-W62. 10.1093/nar/gkm384.PubMed CentralView ArticlePubMedGoogle Scholar
  56. Junier P, Junier T, Witzel KP: TRiFLe, a program for in silico terminal restriction fragment length polymorphism analysis with user-defined sequence sets. Appl Environ Microbiol. 2008, 74: 6452-6456. 10.1128/AEM.01394-08.PubMed CentralView ArticlePubMedGoogle Scholar
  57. Vincze T, Posfai J, Roberts RJ: NEBcutter: A program to cleave DNA with restriction enzymes. Nucleic Acids Res. 2003, 31: 3688-3691. 10.1093/nar/gkg526.PubMed CentralView ArticlePubMedGoogle Scholar
  58. Martorell P, Barata A, Malfeito-Ferreira M, Fernandez-Espinar MT, Loureiro V, Querol A: Molecular typing of the yeast species Dekkera bruxellensis and Pichia guilliermondii recovered from wine related sources. Int J Food Microbiol. 2006, 106: 79-84. 10.1016/j.ijfoodmicro.2005.05.014.View ArticlePubMedGoogle Scholar


© Romi et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.