Invasive candidiasis (IC) bears has high risk of morbidity and mortality in the intensive care units (ICU). The Candidia spp was the main cause of hospitalized bloodstream infections and it is associated with mortality highly rate. ICU receive huge number of patients exposed to multi-accidents types such asterrorist attacks; Cerebrovascular accidents(CVA), etc., most of these cases lead to penetration of Candidaspp to blood stream. The aim of the study was prospective, accurate and rapped identification molecular methods for detection of Candida from blood stream and oral infections in intensive care units in Iraq. 400 samples clinically diagnosed candidemia were collected over a period from November / 2015 to April / 2016 from ICU of Margan Teaching Hospital and Hilla Teaching Hospital in the province of Babylon, Iraq, the Middle East. Out of 264 blood samples (53) were culture positive (32n=male,21n=female) ,while 71 out of 136 oral swabs were culture positive. Out of which 63.8% isolates were C. parapsilosis20.34% were C.albicans12.7%,C.membrenifaciens 2.97% and 0.25% for C.sake. Our conclusion: This study gave more attention of the risk factor of candidemia, and showed many Candidaspp were penetrate blood stream and showed unique genetic polymorphism patterns of Candida in Iraq. The new recorded for the C. parapsilosis,C.sake and C.membrenifaciens from blood samples for the first time in Iraq based on sequencing analysis of the whole ITS region and the ITS2 of the rDNA emphasizes the low precise of conventional identification methods.
Keywords : Invasive candidiasis, ITS and ITS2 sequencing analysis, multiple alignment sequence.
Candidemia and others types of Candidiasis are generally common causes of nosocomial infections through the worldwide, the incidence ,prevalence and identification pictures are attention in many countries such as USA ,European and some of Asian countries and ICU have been notarized for C.albicans,C.sake, C.krusei, C.parapsilosis and C.glabrata strain  .
In Intensive Care Units(ICU) Candidia spp. was the 3rd generality common cause of hospitalized bloodstream infections and is associated with death rate of 47% .C.albicansis considered as the more serious Candida spp and main cause of candidiasis, oral thrush, Candiduria and vaginal candidiasis  .Although the C.albicans has long been the most common types isolated during blood stream infection, recent reports has shown a shift to Non Candida Albicans(NCA) like C.sake, C.krusei, C.parapsilosisand C.glabrata [6-8], especially in ICU and patients  . Many studies showed the C.albicans and C.glabrata representative a major threat to the patients in hospitals worldwide  . NCAcandidemia has developed a rising significant infection and it to replace C.albicans especially uncommon species such C.sake and C.mambranifaciens in most clinical sites any bloodstream infections  .
Huge number of patients exposed to multi-accidents types lead to penetration of Candidaspp to blood stream as in related to inquiries in car accident and war or terrorist attacks. Rodriguez et al., studied the risk factors associated with invasive fungal infections in combat trauma. Blyth et al.  identified many Candidaspp colonization and infection of Combat-Related Injured Patients of USA from Iraq and Afghanistan.
A conventional diagnostic methods such as CHROMagar,fermentation and other phenotypic testwere consuming time and gave bias diagnoses of Candida isolates to the species level [16-18] reports that the CHROMagar unreliable test for distinguish between C.glabrata and C.parapsilosis, both showed white-pink also the CHROMagar showed confuse color to distinguish between many species like C.kefyr, C.utilis, C.robusta, C.famata, C.rugosa, C.guilliermondiiand C.pelliculosa. All make the same kind of glossy pink colonies as C.glabrata, consequently misidentification can happen  .
Advance molecular methods especially PCR and sequencing analysis is being increasingly used for the rapid diagnose of Candida isolates to the species level . The most widespread targets of PCR amplification are rDNA genes included typing of ITS regions and ITS2 regions of fungi were beneficial for the rapid identification of clinically important fungi . DNA sequencing is reinforcement new discovery which revolt the conceptual foundations of numerous fields .
The main Candidaspp was identified in nosocomial candidemia in Suadia Arabia like C.parapsilosis ,C,tropicalis ,C.krusei [23-25]. Shokohi et al.,  showed many Candidia spp. likeC. parapsilosis, C. glabrata, C. albicans, C. tropicalis,C.kruseandC.guilliermondiiin Iran. Also Vijayakumar et al., was identified C. parapsilosis was the common yeast isolated from IC patients.C. glabrata, C. albicans, C. tropicalis, and C. kefyr, were identified from candidemia patient in India based on PCR and RFLP-PCR .but few attention in the Middle East .
Unfortunately ,the picture of candidemiais no clear till now in Iraq due to absence of any attention and lack of real studies on candidemia except the short study of Sabeeh et al,  diagnosed the C.albicans based on conventional methods as causes agent of candidemia in leukemia in one province of Iraq. Risan  isolated C.glabrata only from blood samples of acute leukemia in Baghdad province. On the other hand, the intensive care units patients of the most hospitals in Iraq were receiving injured persons due to car accident or war or terrorist attacks. The variant human accidents lead to penetration of Candida spp. to blood stream and may be reflective of a population at higher risk for patients severity and death.
The aim of this study was prospective, accurate and rapped identification molecular methods for detection of Candida spp. from blood stream and oral infections in intensive care units in Iraq .
2. Material and methods
Four hundreds clinical samples (246 blood, 77 swab immune patients in cancer 23 swab Hematology Unit, 36 swab patients in ICU; were collected based on standard methods  .The samples were collected from patient. These samples collected with different age groups ranging from (4 to 97 years). The diseases associated with immune deficiency included attenuated and hospitalized Margan Teaching Hospital and Hilla Teaching Hospital in the province of Babylon, the survey covering the period from November / 2015 to April / 2016 was performed.
Two milliliters of blood were collected from each patient, blood put in EDTA tube for keeping blood from clotting, then incubated blood cultures at 37Co for 24h to 96h on Sabouraud’s dextrose agar ( SDA).Blood and swab were streaking on SDA (12),single colonies were subjected to species identification basedon character’s color production on CHROMagar medium .
2.1. PCR assay
The phenotypic results were confirmed by simple PCR by specific primer pair for Candida. One μLof DNA(20µg/ml)from each of 20 Candida isolates were mixed with PCR mixture ( final reaction volume 25 μL) consisted of 12 μL of 20x Master Mix (Promega),2 μL of primers (10 pmole) and restmolecular-gradewater. The PCR conditions for primer pairs ITS5/ITS4,ITS3/ITS4 primers were: 95 ºC for 3 min followed by30 cycles 94 ºC for 1min.55ºC for 40sec. and 72ºC for 1 min.and final extension 72ºC for 5 min.The PCR mixture was amplified by thermal cycler PCR System (Labnet, USA) .
The PCR products for each target region were run on 1.2% agarose gel (Bio Basic Canada Inc.).Electrophoresis performed at 100 V. in TBE buffer. The gel was pre-stained with 0.05% ethidium bromide. The DNA bands were detected by Desktop Gel imager scope  ultraviolet transilluminator (Korea Com.). The purified products for 17 isolates were subjected to sequencing. Sequence analysis was performed in MacrogeneLab., USA. The sequence alignment of Candida spp was compared with the BLAST data base and were aligned with sequences from the BLAST data base derived from the following reference strains. The phylogenetic tree(UPGM) based on sequencing were constructed employing the Mega 6 software, pairwise and multiple sequence alignment based on BioEdit software was performed .
All isolates of Candida spp. were preliminary classified based on color on CHROMagar, later only 24 isolates were subjected for PCR assays amplified whole ITS region ITS2 region by primer pair ITS5/ITS4 and ITS3/ITS4 respectively (Figure 1).Only seventeen of the PCR products for both ITS regions of selectively isolates based on similarity and difference in amplicons sizes were sequenced (Table 1).
[Figure is omitted from this preview]
Figure (1): Agarose gel electrophoresis of PCR products for Candida spp isolates amplified by primer pair (ITS5-ITS4) and primer pair (ITS3-ITS4).
3.1. Sequence analysis
Five species of Candida out of 17 sequence were identified based on pairwis alignment with reference isolates in the gene bank ,they are C.sake ,C.prapsilosis ,C.membrenifaciens, C.albicans and the telomorphe of C.kruse is Pichiakudriavzevii for both target region (ITS and ITS2 respectively(Table 1).C.sake and C.membrenifaciens were recorded for the first time in Iraq.
3.2. Multiple alignment of 17 sequences for both ITS and ITS2 regions
The results of Multiple alignment analysis based on BioEdit software was performed for 17 isolates of Candida spp. (5n for C.sake,2n for C.albicans,4n for both C.parapsilosis and C.membrenifaciens and one isolate for Pichiakudriavzevii (Figures 2-3) . Each set of isolates for Candidaspp was showed high similarity with leading sequences with some mutation or substitutions, these sequence variation were indicated to microevolution in each set of isolates.
[Figure is omitted from this preview]
Figure (2): The Sequence Multiple-alignment of sequence analysis of ITS2 region amplified by ITS5/ITS4.
High similarity of multiple alignment partial sequence of ITS and ITS2 regions amplified by ITS5/ITS4 and ITS3/ITS4 primer pairs respectively. Seventeen and 16 isolates of Candida spp. (6n and 5n for C.sake,2n for C.albicans,4n for both C.parapsilosis and C.membrenifaciens and one isolate for Pichiakudriavzevii for to regions respectively. All isolates of C.sake showed highly similarity with leading sequence with low variation in some sequences sites. These variations may be Single Nucleotide Polymorphism( SNP) or bases substitution ,more variation was observed at the end of sequences of C.sake (Figures2,3).Other species showed intra or inter- species variations, sometimes similarity in all species may due to the uniformity of belonging to same Candidagenus.
The phenotypic identification assays showed fluctuation results compare with molecular identification assays were showed reliable coincidence results, many isolates of Candida spp. were identified as C.kefyr based on their reaction color on CHROMagar while identified as C.membranifaciens based on whole ITS region and ITS2 in rDNA region, and also with some isolates of C.albicans, and C.albicans Ssp.dubliniesis (Table 1).
[Figure and figure are omitted from this preview]
Figure (3): Multiple-alignment of sequence analysis of ITS2 region amplified by ITS3/ITS4.
Table (1): Comparison among of phenotypic identification and molecular identification based on sequence analysis data.
3.3. Phylogeny tree
The phylogenetic tree (UPGM) for16 Candidaspp was constructed based on sequences of ITS region. The isolates of Candidaspp were isolated from cancer patients showed closed related intra-isolates groups as in Figure (4), results 5 clusters of Candida spp. the observation show C.sake in cluster 1, C.parapsilosis in cluster2, C.albicans in cluster 3, C.membranifaciens in cluster 4 and Pichiakudriavzevii in the cluster 5.
[Figure is omitted from this preview]
Figure (4): Phylogeny tree (UPGM) based on sequence of 16 Candidasppisolates.
The results of bioinformative pairwise sequence alignment data for two each of Candida species were showed distinct of each species based on the data of the score and the Query cover percent as in Figures (5-8). The pairwise sequence alignment data of two C, albicans isolates showed high score 558 and 91% of query cover, this means belonging for same species ,while the same bioinformative between the others species showed them as a distinct species for example: The pairwise sequence alignment data of two , C.membrenifaciens and C.albicans showed high score 159 and 37% of query cover, this means belonging distinct species (Figure 5 )
[Figures are omitted from this preview]
Figure (5): The pairwise sequence alignment data of two, C.albicans and C.albicans showed high score 558 and 91% of query cover, this means belonging to the same species.red cycles mismatches ,blue cycles gaps .
Figure (6): The pairwise sequence alignment data of two, C.membrenifaciensand C.albicans showed high score 159 and 37% of query cover, this means belonging distinct species.
Figure ( 7): The pairwise sequence alignment data of two , C.sake with C.albicans showed high score 252 and 48% of query cover, this means belonging distinct species.red cycles mismatches ,blue cycles gaps .
Figure (8): The pairwise sequence alignment data of two , C.sake with C.parapsilosis showed high score 265 and 43% of query cover, this means belonging distinct species . red cycles mismatches ,blue cycles gaps .
Our results of isolation and identification of C.sake and C.membrenifaciens and other species in the blood samples of ICU patients were coincidence with recently studies .The final identification of 17 Candida spp.(Table 1) showed that molecular markers more accurate for identification based on [33-35]. The missing identification based on phenotypic markers being fluctuation,CHROMagar considered as preliminary identification test but not accurate forever, this interpretation was reviewed in  .
The variation of PCR product for both region showed coincidence. The identification of Candidaspp was followed(Fujita et al.36 andBellemainet al. . These sequences variations may raised due to micro evolutionary trends in any species groups,Imran and Ali18 mention to the same sequences variation in C.glabrata and C.parapsilosis.
Due to the fluctuation of the identification of Candidaspp based on CHROMagar for example: some of isolates were identified as C.albicans based on green color on CHROMagar but identified as C.sake based on sequences analysis, the same case of C.albicans Ssp.dublinies is identified as C.albicans based on Sequencing methods (Table 1). These confusion in identification between Candida spp. Wassolved the confuse by pairwise sequence alignment between two each species, and the results of paiwise alignment showed genetic variation and microevolution degree for each alignment (Figures 5-8).
Our finding, with regards to the lacking previous data on Candidemia in Iraq, the finding C.parapsilosis, C.sake C.membrenifaciens and C.albicans being the most frequent recovered species in this study is indicative of unique trait specific to the epidemiology of candidemia in ICU. Moreover, similar to our finding, many studies were reported the isolation most of the Candida spp. from neighboring countries from Sudia Arabia , Gulf region  and Iran  and other countries in Europe and USA .
Most of previous identification studies were based on small sample sizes ranged from 30-50 clinical isolates such as Vijayakumaret al.,  based on  consecutive clinical isolates patients from various ICUs .
Our study suffered from several limitations. First, because we used blood culture positivity as the criterion for diagnosing candidemia instead of direct PCR detection, this estimation could be less than the actual frequency of candidemia. Second, small samples number makes our results less accurate to estimate the actual prevalence of candidemia in Iraq. Finally, we faced difficulties in the following up our patient groups for along period due to the critical patient's cases,such as death of the patients and built our explanation on the results of culture positive diagnosis.
On the other hand the misidentification of blood septicemia ,which may be due to the systematic candidiasis causes problems in ICU for patients, may be indistinguishable from bacterial septicemia, and these infection are severe difficult to diagnose may lead to severely disease from 1-year to elderly patients. The major risk factors for developing IC include prior antibiotic use, central venous catheterization, a recent major surgery, use of steroids, dialysis and immunosuppression . According to a recent survey, 50-80% of the critically ill patients who were admission to ICU had already been exposed to risk factors for IC, 5-15% had Candida colonization on admission and 5-30% actually had IC (Leroy et al.,2009 ) .
The results of this study were agreed with review of Yapar40that only 15 out of 150 Candida spp isolated from patients as infectious agents , it has been determined that in 95% of infections, involved C. albicans, C. glabrata, C. parapsilosis, C. tropicalis, and C. krusie.The contribution of this study showed the foggy picture of candidemia in Iraq ,and emphasis the requirement for futures studies about candidemia and IC incidence compare with huge studies in development countries were performed .
Our conclusion,in spite ofthe ICU receiving multi-hundreds of patients in each year, undergo any one of multi-accidents types such as car accidents, wars accidents ,terrorist attacks and Cerebrovascular accidents(CVA) cases and the patients were exposed to different levels of injuring cases and contaminated with fungi and bacteria . Unfortunately the studies in Iraq gave low attention of identification species caused candidemia based and depending on simple conventional methods for identification compare with finding our study, due to these reasons the value of this study still low and the picture of IC is ambiguous because the absence of pervious interest studies in Iraq ,also the local epidemiological data on fungal infections in ICU staff hospitals are lacking . Our perspective of the incidence and the impact of IC in our ICUs should be further clarified in future investigations.
The first author supervised on the work and performed the molecular analysis and final editing of this article, while the second author performed the practical part and written and pre-editing this article. The works were conducted in the Biotechnology Laboratory, The Department of Biology, All Women College of Science, Babylon University, Iraq. The authors would like to thank ICU staff of Margan Teaching Hospital and Hilla Teaching Hospital for samples collection and thank hospital for their logistic support and their help in sampling.
 Hedderwick, S. A. ; Lyons, M. J. ; Liu, M. ; Vazquez, J. A. and Kauffman CA. (2000). Epidemiology of yeast colonization in the intensive care unit.Eur J ClinMicrobiol Infect Dis. 19(9):663-70.
 Yang, Z.T. ; Wu, L. ; Xiao-Ying Liu, X. Y. ; Min Zhou, M. ; Jie Li, J. ; Wu, J. Y. ; Yong Cai, Y. ; En-Qiang Mao, E. ; Chen, E. and Lortholary, O. (2014). Epidemiology, species distribution and outcome of nosocomial Candida spp. bloodstream infection in Shanghai.Yang et al. BMC Infectious Diseases. 241(14):1-10.
 Malani, A.; Hmoud, J.;Chiu, L.; Peggy, L.; Carver ; Bielaczyc, A. and Carol A. K. (2005). Candida glabrata Fungemia: Experience in a Tertiary care center.CID.41:975-981.
 Imran, Z . K. and Abuad, S. H. (2015). Genetic Diagnosis and Prevalence of Urinary Tract Fungal Pathogen with Antifungal Susceptibility Pattern in Iraq BJMMR 7(5): 410-418 .
 Imran, Z.K, Al.Shukry, H.N. (2014). Molecular diagnosis of vaginal candidiasis by polymerase chain reaction (PCR) and random amplification polymorphism DNA(RAPD-PCR)in Babylon Province, Iraq. Afr. J.Microb. Reach. 8(6):496-502.
 Chang, A. ;Neofytos, D. and Horn, D. (2008).Candidemia in the 21st century. Future Microbiol. 3(4):463–472.
 Bassetti, M. ;Taramasso, L. ; Nicco, E. ; Molinari, M. P. ; Mussap, M. and Viscoli, C. (2011). Epidemiology, species distribution, antifungal susceptibility and outcome of nosocomial candidemia in a tertiary care hospital in Italy.J PLoS One. 6(9):e24198 :1-6.
 Nucci, M.;Queiroz-Telles, F.; Alvarado-Matute, T.;Tiraboschi, I.N.; Cortes, J. et al.(2013) Epidemiology of candidemia in Latin-America: a laboratory-based survey. PLos One 8: e59374.
 Zirkel, J.; Klinker, H.; Kuhn, A. ; Abele-Horn, M.; Tappe, D. ; Turnwald, D.; Einsele, H. and Heinz, W. J. (2012). Epidemiology of Candida blood stream infections in patients with hematological malignancies or solid tumors. Med Mycol. 50:50–5.
 Montagna, M. T. ; Caggiano, G. ; Lovero, G. ; De Giglio, O. ; Coretti, C. ; Cuna, T. ; Iatta, R. ; Giglio, M. ; Dalfino, L. ; Bruno, F. and Puntillo, F. (2013). Epidemiology of invasive fungal infections in the intensive care unit: results of a multicenter Italian survey (AURORA Project) Infection. PMID. 41(3):645–653.
 Wisplinghoff, H.; Bischoff, T.; Tallent, S. M.; Seifert, H.; Wenzel, R. P. and Edmond, M. B. (2004). Nosocomial bloodstream infections in US hospitals: analysis of 24,179 cases from a prospective nationwide surveillance study. Clin Infect Dis. 39(3):309-17.
 Sabeeh, S. ; Al-Attraqhchi, A. A. F. and Al-Aswad, E. (2013). PCR in Comparison with Culture Methods for The Diagnosis of Candida albicans Responsible for Candidemia in Leukemic Patient.Diyala Journal of Medicine. 5(2):29-35.
 Aghili, S. R. ;Shokohi, T. ; Boroumand. M. A. ;Fesharaki, S. H. and Salmanian, B. (2015). Intravenous Catheter-Associated Candidemia due to Candida membranaefaciens: The First Iranian Case J TehUniv Heart Ctr. 10(2):101-105.
 Rodriguez, C.J.;Weintrob, A.C.; Shah, J. et al.(2014). Risk factors associated with invasive fungal infections in combat trauma. Surg Infect (Larchmt) 15:521–6.
 Blyth, D. M.;Mende, K.;Weintrob ,A. C.;Beckius, M.L.;Zera, W. C.; Bradley, W.; Lu, D.;Tribble, D. R. and Murray, C. K. (2014).Resistance Patterns and Clinical Significance of Candida Colonization and Infection in Combat-Related Injured Patients From Iraq and Afghanistan.Open Forum Infect Dis. 1(3):1-10.
 Imran ,Z. and Al- asadi, Y. (2014). Multiple molecular markers for diagnosis of conjunctivitis caused by Candida spp. in Iraq, Afr. J. Microbiol. Res. 8(38): 3482-3488.
 Imran, Z. K.(2015).Candida albicans Ssp. dubliniensis stat.et comb. nov., a new combination for Candida dubliniensis based on genetic criteria . Afr. J. Microbiol. Res.,9(17):1205-1214.
 Imran, Z . K. and Ali, E. K. (2015).Molecular identification of Candida glabrata and C. parapsilosis based on sequencing analysis of rDNA .International J. of Med. Sci. and Clin. Inventions .2(12):1490-1497.
 Baumgartner, C. ;Freydiere, A.M. and Gille, Y. (1996). Direct identification and recognition of yeast species from clinical material by using Albicans ID and CHROMagar Candida plates.J ClinMicrobiol. 34:454-456.
 Freydiere, A.M. ;Buchaille, L. andGille, Y. (1997). Comparison of three media for direct identification and discrimination of Candida species in clinical specimens.Eur. J ClinMicrobiol. Infect. Dis. 16:464-467.
 Turenne, C. Y.;Sanche, S. E.; Hoban, D. J.;Karlowsky, J. A.;Kabani, A. M.(1999).Rapid identification of fungi by using the internal transcribed spacer 2 genetic region and an automated fluorescent capillary electrophoresis system. J ClinMicrobiol. 37:1846–1851.
 Franca, L. T. C. ;Carrilho, E. and Kist, T. B. L. (2002). A review of DNA sequencing techniques.Quarterly Reviews of Biophysics. 35(2):169–200.
 Bukharie, H.A.(2000).Nosocomial Candidemia in a tertiary care hospital in Sudia Arabia .Mycopathologia 154(4):195-198.
 Mokaddas, E. M. ; Ramadan, S. A. ; Abo el Maaty, S. H. and Sanyal, S. C. (2000). Candidemia in pediatric surgery patients.J Chemother. 12:332-8.
 Al-Tawfiq, J. A. (2007). Distribution and epidemiology of Candida species causing fungemia at a Saudi Arabian hospital, 1996–2004.Int J Infect Dis. 11:239–44.
 Shokohi,T.;HashemiSoteh,M.B.;Pouri,Z.S.;Hedayati,M.T.;Mayahi,S.(2010).Identification speciesusing PCR-RFLPin cancer patients inIran.IndianJ.Med.Microbiol.28:147-51.
 Vijayakumar,R.;Giri,S. and Kindo, A. J. ( 2012).Molecular species identification of Candida from Blood Samples of Intensive Care Unit Patients by Polymerase Chain Reaction – Restricted Fragment Length Polymorphism.J Lab Physicians. 4(1): 1–4.
 Taj-Aldeen, S. J. ; Kolecka, A.;Boesten, R. ; Alolaqi, A.; Almaslamani, M.; Chandra, P.; Meis, J. F. and Boekhout, T. (2014). Epidemiology of candidemia in Qatar, the Middle East: performance of MALDI-TOF MS for the identification of Candida species, species distribution, outcome, andsusceptibility pattern. Clinical And Epidemiological Study. Springer-Verlag Berlin Heidelberg. Volume 42(2):393–404.
 Risan,M.H.2016.Molecular identification of yeast Candida glabrata from candidemia patients in Iraq.IraqJ.of Science.57(2):808-813.
 Rao, P. K. (2012). Oral Candidiasis – A Review Scholarly Journals International Scholarly Journal of Medicine. 2(2): 26-30.
 Nadeem, S.G.; Hakim, S.T.;Kazm, S.U. (2010).Use chromoagar candida medium for the presumptive identification of Candida species directly from clinical specimens in resource –limited setting. Libyan J. Med. 5:1-6.
 Kurzmann, C. P. ; Fell, J. W. and Boekhout, T. (2011). The yeasts: a taxonomic study. 5th ed. Amsterdam: Elsevier.
 Peng, C. F. ; Lee, K. M. and Lee, S. H. (2007). Characterization of Two Chromogenic Media of Candida ID2 and CHROMagar Candida for Preliminary Identification of Yeasts.J Biomed Lab.
 Bellemain, E. ;Carlsen, T. ; Brochmann, C. ; Coissac, E. ; Taberlet, P. and HåvardKauserud, H. (2010). ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiology. 189 (10): 1-9.
 Fujita, S.I;Senda, Y.;Nakaguchi, S.; Hashimoto, T.( 2001). Multiplex PCR using internal transcribed spacer 1 and 2 regions for rapid detection and identification of yeast strains. J ClinMicrobiol 39:3617–3622.
 Bellemain, E.;Carlsen, T.;Brochmann, C.;Coissac, E.;Taberlet, P.;Kauserud, H.(2010). ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiology. 10:189
 Osorio, J.J.; Roman, A.R.; Torre-Cisneros, J. (2007).Spectrum and risk factors of invasive fungal infection]. EnfermedadesInfecciosas y MicrobiologiaClinica 25: 467–476.
 Leroy, O.;Gangneux, J.P.;Montravers,P.;Mira, J.P.;Gouin, F.;Sollet, J.P, et al.( 2009 ) Epidemiology, management, and risk factors for death of invasive Candida infections in critical care: a multicenter, prospective, observational study in France (2005–2006). Crit Care Med 37: 1612–1618.
 Yapar, N.(2014).Epidemiology and risk factors for invasive candidiasis.TherClinRisk Manag.10:95-105.