ResultsWe obtained 110 Candida strains, which included 40 Candida albicans (36.36%), 37 C. parapsilosis (33.64%), 21 C. tropicalis (19.09%), 9 C. glabrata (8.18%), 2 C. rugose (1.82%), and 1 C. haemulonii (0.91%) isolates. At a limiting point of 0.80, C. albicans isolates could be grouped into five clusters, C. parapsilosis and C. tropicalis isolates into seven clusters, and C. glabrata isolates into only one cluster comprising six strains by RAPD typing. Antifungal susceptibility testing revealed that the isolates showed the greatest overall resistance against fluconazole (6.36%), followed by voriconazole (4.55%). All C. albicans and C. parapsilosis isolates exhibited 100% susceptibility to echinocandins (i.e., anidulafungin, caspofungin, and micafungin), whereas one C. glabrata strain was resistant to echinocandins. The most common amino acid substitutions noted in our study was 132aa (Y132H, Y132F) in the azole-resistant strains. No missense mutation was identified in the hotpot regions of FKS1. Comparison of the selected virulence factors detectable in a laboratory environment, such as biofilm, caseinase, and hemolysin production, revealed that most Candida isolates were caseinase and hemolysin producers with a strong activity (Pz < 0.69). Furthermore, C. parapsilosis had greater total biofilm biomass (average Abs620 = 0.712) than C. albicans (average Abs620 = 0.214, p < 0.01) or C. tropicalis (average Abs620 = 0.450, p < 0.05), although all C. glabrata strains were either low- or no-biofilm producers. The virulence level of the isolates from different specimen sources or clusters showed no obvious correlation. Interesting, 75% of the C. albicans from cluster F demonstrated azole resistance, whereas two azole-resistant C. tropicalis strains belonged to the cluster Y.