Mohammad Asadzadeh
Suhail Ahmad
Noura Al-Sweih
Ziauddin Khan- Department of Microbiology, College of Medicine, Kuwait University, Jabriya, Kuwait
Backgrounds: Candida glabrata is a frequently isolated non-albicans Candida species and invasive C. glabrata infections in older patients are associated with high mortality rates. Opportunistic Candida infections in critically ill patients may be either endogenous or nosocomial in origin and this distinction is critical for effective intervention strategies. This study performed multi-locus sequence typing (MLST) to study genotypic relatedness among clinical C. glabrata isolates in Kuwait.
Methods: Candida glabrata isolates (n = 91) cultured from 91 patients were analyzed by MLST. Repeat isolates (n = 16) from 9 patients were also used. Antifungal susceptibility testing for fluconazole, voriconazole, caspofungin and amphotericin B (AMB) was determined by Etest. Genetic relatedness was determined by constructing phylogenetic tree and minimum spanning tree by using BioNumerics software.
Results: Resistance to fluconazole, voriconazole and AMB was detected in 7, 2 and 10 C. glabrata isolates, respectively. MLST identified 28 sequence types (STs), including 12 new STs. ST46 (n = 33), ST3 (n = 8), ST7 (n = 6) and ST55 (n = 6) were prevalent in ≥4 hospitals. Repeat isolates obtained from same or different site yielded identical ST. No association of ST46 with source of isolation or resistance to antifungals was apparent. Microevolution and cross-transmission of infection was indicated in two hospitals that yielded majority (57 of 91, 67%) of C. glabrata.
Conclusion: Our data suggest that C. glabrata undergoes microevolution in hospital environment and can be nosocomially transmitted to other susceptible patients. Thus, proper infection control practices during routine procedures on C. glabrata-infected patients may prevent transmission of this pathogen to other hospitalized patients.
1. Introduction
Invasive fungal infections (IFIs) are mostly caused by opportunistic yeast pathogens in nosocomial settings and are associated with high mortality rates in critically ill patients (1–4). The incidence of IFIs is increasing globally due to increasing use of antifungal prophylaxis in an expanding population of immunocompromised/immunosuppressed individuals and other susceptible patients with multiple comorbidities, particularly in intensive care unit (ICU) settings (1, 2, 5, 6). Several novel yeasts with reduced susceptibility or resistance to antifungal drugs have emerged in recent years (7, 8). The increasing incidence of drug resistance in commonly encountered yeasts is also causing major shifts in the epidemiology of IFIs (5, 9–14). Candida spp. are also recognized as an important cause of late onset septicemia in low/very low-birth-weight neonates (15). Although Candida albicans is the most pathogenic and most frequently isolated species from clinical specimens, infections caused by non-albicans Candida species (NACS) have increased considerably in the last two decades and more than 50% of all Candida infections are now caused by NACS (5, 13, 16, 17). Widespread use of azole antifungal drugs and improvements in diagnostic methods have likely contributed to the increased proportion of NACS infections in recent years (5, 13, 18). Candida glabrata, along with two other closely related species; Candida bracarensis and Candida nivariensis, now form part of the Nakaseomyces clade, and hence, have recently been transferred to a new genus, Nakaseomyces (19). The emergence of C. glabrata (N. glabrata) as a common yeast pathogen is due to the intrinsic and acquired resistance of this species to azoles, the most widely used class of antifungal agents (20). Candida glabrata is the second or third most frequently isolated yeast species from blood and other clinical specimens from immunocompromised, cancer and diabetic patients, its isolation frequency is higher from older (>60 years) patients and those with indwelling medical devices and invasive infections are associated with higher mortality rates compared to infections with other Candida spp. (5, 13, 21, 22).
Opportunistic Candida infections in critically ill patients may be either endogenous or nosocomial in origin and this distinction is critical for effective intervention strategies (23–25). Several molecular methods have been used previously to determine genotypic similarities/differences among clinical C. glabrata strains (26–30). However, only few robust and/or portable methods namely, multi-locus sequence typing (MLST) based on sequence analyses of 5–7 conserved housekeeping loci, multi-locus microsatellite typing (MLMT), polymorphic locus sequence typing (PLST) and whole-genome sequencing (WGS) have been used in recent years (24, 25, 31). Although WGS provides the most comprehensive information and resolution, it is still costly for routine use in clinical settings and requires sophisticated data analyses while non-coding region-based methods (PLST) do not fulfill the evolutionary criteria which are exhibited by the conserved housekeeping loci (24, 26, 27). In this study, we performed MLST of C. glabrata isolates collected from hospitalized patients to determine its genetic structure, microevolution and nosocomial transmission of infection in Kuwait.
2. Materials and methods
2.1. Reference strains, clinical isolates, growth conditions and phenotypic and molecular identification
The C. glabrata (ATCC 90030 and CBS 135), C. nivariensis (CBS 9983), C. bracarensis (CBS 10154), C. albicans (ATCC 76615) and Candida parapsilosis (ATCC 22019) were used as reference strains. A total of 91 clinical C. glabrata isolates cultured from various clinical specimens and collected during 2007 to 2015 from 91 patients hospitalized in 9 different major government hospitals in Kuwait were used. These isolates were recovered from urine (n = 29), oral/respiratory specimens (n = 27 including 14 sputum, 9 tracheal aspirate, 3 oral swabs and 1 bronchoalveolar lavage samples), blood (n = 21), wound swab (n = 5), rectal swab (n = 4), pleural fluid/abdominal fluid/gastric aspirate (n = 3), axilla swab (n = 1) and vaginal swab (n = 1). Repeat isolates (n = 16) collected from the same or different anatomic sites of 9 patients were also used to ascertain the robustness of genotyping data.
The clinical specimens yielding C. glabrata were collected from the hospitalized patients after obtaining informed verbal consent only as part of routine diagnostic work-up for the isolation of bacterial and fungal pathogens. The data on deidentified samples are reported in this study. The clinical specimens were cultured by using BACTEC 9240 (Becton Dickinson, Sparks, MD, United States), BacT/Alert 3D (bioMérieux, Marcy-l’Étoile, France) and/or Versa TREK™240 (Thermo Fisher Scientific, United States) for the isolation of yeasts, as described previously (14, 17). All growth-positive cultures were subcultured on Sabouraud dextrose agar and Mast ID-CHROMagar Candida (Mast Diagnostics, Merseyside, United Kingdom) for phenotypic colony characteristics, as described previously (32). Species-specific identification was achieved by assimilation profiles obtained by commercial Vitek2 yeast identification system and/or by protein profiles by MALDI-TOF MS (bioMérieux, Marcy-l’Étoile, France), as described previously (33, 34). For molecular identification, the genomic DNA was extracted by using the Gentra Puregene Yeast DNA extraction kit (Qiagen, Hilden, Germany) according to kit instructions or by the rapid method using Chelex-100 as described previously (35). Molecular identification was performed by a multiplex PCR assay that discriminates C. glabrata sensu stricto from C. nivariensis and C. bracarensis among C. glabrata sensu lato strains, as described previously (36). The identity of 70 randomly selected isolates was further confirmed by DNA sequencing of internal transcribed spacer (ITS) region of rDNA by using panfungal ITS1 and ITS4 primers, performed as described previously (37). The study was approved by the Health Sciences Center Ethics Committee (vide approval No. VDR/EC/3724), Kuwait University and all the methods and investigations reported in this study were carried out according to their guidelines.
2.2. Antifungal susceptibility testing
In vitro antifungal susceptibility testing (AST) was conducted by Etest (bioMérieux SA, Marcy-l’Étoile, France) in accordance with the manufacturer’s instructions and as previously described (38). C. parapsilosis (ATCC22019) and C. albicans (ATCC90028) reference strains were employed for quality control. The minimum inhibitory concentration (MIC) values determined at the point of 80% growth inhibition for azoles and complete inhibition for amphotericin B (AMB) were read after 24 h. Interpretive susceptibility criteria for fluconazole (FLU), voriconazole (VOR) and caspofungin (CFG) were those recommended by Clinical and Laboratory Standards Institute (CLSI) (FLU: ≤32 μg/mL, susceptible dose-dependent, ≥64 μg/mL, resistant, VOR: ≤2 μg/mL, susceptible, >2 μg/mL, resistant, CFG: ≤0.12 μg/mL, susceptible, 0.25 μg/mL to <0.5 μg/mL, intermediate, ≥0.5 μg/mL, resistant) (39, 40) and for AMB were: ≤1 μg/mL, susceptible, >1 μg/mL, resistant (11).
2.3. Multi-locus sequence typing of Candida glabrata isolates
A total of six (FKS2, LEU2, NMT1, TRP1, UGP1, and URA3) housekeeping gene fragments accepted as an international standard for C. glabrata MLST1 were used (26). PCR was used to amplify each gene fragment from each isolate and PCR amplicons were purified by using PCR product purification kit (Qiagen Inc., Valencia, CA, United States) according to kit instructions. Both strands were sequenced using the BigDye terminator DNA sequencing kit (Applied Biosystems) with the same primers used for PCR amplification. The sequencing products were processed and loaded onto the DNA sequencer’s capillary as instructed by the manufacturer (Applied Biosystems) and as described previously (41). Clustal omega2 was used to verify forward and reverse allele profiles. The reproducibility of the MLST protocol was ensured by including the reference C. glabrata strain (ATCC 90030) in each run. The C. glabrata database was used to determine the allelic profile (allele number) for each gene. The allelic combinations for all the six loci for each isolate (sequence type, ST) were assigned or new allele and new STs were provided by the curator, Prof. Oliver Bader of the C. glabrata MLST database. The new allelic profiles and STs obtained in this study were submitted to the MLST database. The isolates belonged to the same ST if they contained the same allelic combinations for all the six loci. Based on the allele number for each isolate’s six gene fragments, a dendrogram was created using BioNumerics v7.5 (Applied Maths, Sint-Martens-Latem, Belgium) using the standard unweighted pair group method with arithmetic mean (UPGMA) parameters as described previously (41). A cluster was defined as a group of two or more C. glabrata isolates belonging to the same ST. The genetic relationship among the isolates was determined by constructing a minimum spanning tree using BioNumerics software which predicts putative relationships among the isolates and records two isolates as more closely related when five of six loci are identical (34, 41). Genotypic relationship between Kuwaiti and Asian C. glabrata isolates was also investigated.
2.4. Statistical analyses
The distribution of STs with respect to different clinical specimens or hospitals was compared by using Fisher’s exact test and probability levels <0.05 by the two-tailed test were considered as significant. Statistical analyses were performed by using WinPepi software ver. 11.65 (PEPI for Windows, Microsoft Inc., Redmond, WA, United States).
3. Results
3.1. Molecular identification of Candida glabrata isolates
A total of 1,149 C. glabrata isolates were cultured from 1,149 hospitalized patients in different government hospitals across Kuwait during 2007–2015. Repeat isolates cultured from the same and/or different anatomical sites were also available from some patients. More than 400 of these isolates were characterized in detail by various molecular methods including PCR-sequencing of rDNA (36). A total of 91 isolates from 91 patients and 16 repeat isolates from 9 of the same 91 patients were selected from this pool for fingerprinting studies by MLST (Supplementary Table S1). The largest number of isolates came from Hospital A (n = 33) followed by Hospital B (n = 24), Hospital C (n = 12), Hospital D (n = 7), Hospital E (n = 6), Hospital F (n = 4) and Hospital G (n = 3) while only one isolate each was included from Hospital H and Hospital I (Supplementary Table S1). The vast majority (74 of 91, 81.3%) of the isolates analyzed in this study were obtained in 2008 (n = 9), 2009 (n = 14), 2010 (n = 36) and 2011 (n = 15) and represented 20.4% (74 of 362) of all the isolates collected during this period. Repeat isolates (n = 9) from five (P26, P55, P56, P58, and P60) patients were from the same anatomical site as the first isolate while repeat isolates (n = 7) from four (P57, P59, P79 and P80) patients included isolates from different anatomical sites (Supplementary Table S1). All isolates were identified as C. glabrata sensu lato by Vitek2 yeast identification system and as C. glabrata sensu stricto by the MALDI-TOF MS and/or multiplex PCR. The identification was further confirmed by PCR-sequencing of 70 selected isolates as their ITS region of rDNA sequences showed <1% difference with the corresponding sequence from reference C. glabrata strains CBS135 and ATCC 90030, as expected.
3.2. Antifungal susceptibility testing of Candida glabrata isolates
The results of AST for C. glabrata isolates to four antifungal agents are shown in Table 1. The isolates exhibited MIC50 and MIC90 values of 0.11 and 1.5 μg/mL for AMB, 8 and 32 μg/mL for fluconazole, and 0.25 and 0.75 μg/mL for voriconazole, and 0.125 and 0.19 μg/mL for caspofungin, respectively. Majority (84 of 91, 92.3%) of C. glabrata isolates were susceptible dose-dependent while 7 isolates were resistant to fluconazole. Resistance to voriconazole was rare as only 2 isolates were resistant while the remaining 89 isolates were susceptible (Table 1). Eighty-one (89%) C. glabrata isolates were susceptible to AMB while 10 isolates were resistant. Eighty-three isolates were susceptible while 8 isolates were intermediate for caspofungin. Collectively, two (Kw1018/12 and Kw129/12) isolates were resistant to three drugs (AMB, fluconazole and voriconazole).
Table 1. Minimum inhibitory concentration (MIC) ranges and susceptibility data for 91 Candida glabrata isolates by Etest.
3.3. Genotyping of Candida glabrata isolates by MLST
All six housekeeping gene fragments were successfully amplified by using the corresponding gene-specific forward and reverse primers. The DNA sequences were used to obtain the allelic profiles for each gene fragment and allele combinations to obtain STs from the C. glabrata MLST database. TRP1 had the best typing efficiency (14 distinct genotypes/alleles), whereas URA3 yielded the lowest typing efficiency (8 distinct genotypes/alleles). The polymorphic alleles combined to form 28 STs among 91 C. glabrata isolates from 91 patients yielding a ratio of number of STs to the number of C. glabrata isolates of 3.25 (one ST for every 3.25 isolates). The most prevalent ST was ST46 (33/91, 36.3%), followed by ST3 (8/91, 8.8%) and ST7 and ST55, each with a frequency of 6 (6.6%). Of the 28 STs, only 17 STs from 79 isolates were present in the C. glabrata MLST database while 11 STs (ST155 to ST165) from 12 isolates were new and were added to the C. glabrata MLST database. Sixteen (57.1%) of STs were represented by a single isolate which also represented 17.6% of the total 91 isolates. The allelic profile and final ST for each isolate are provided in Table 2. Repeat isolates of each patient collected from the same or different sites yielded the same allelic profile and ST as the first isolate. These included 5 bloodstream isolates from one (P60) patient and blood and non-blood isolates from three (P57, P79 and P80) patients. While the difference in the distribution of ST46, ST7 and ST55 among urine, oral/respiratory and blood specimens was not significant, the distribution of ST3 was significantly different in blood (4 of 21) compared to oral/respiratory specimens (0 of 27) (p = 0.031).
Table 2. Summary of fingerprinting data for 91 C. glabrata isolates using MLST.
The distribution of various STs among C. glabrata isolates collected from patients hospitalized in 9 hospitals are shown in Supplementary Table S1. The most common genotype, ST46 was found in eight hospitals while ST3, ST7, ST55, ST14 and ST122 were found in six, four, four, three and three hospitals, respectively (Supplementary Table S1). The distribution of four major genotypes (ST46, ST3, ST7 and ST55) among C. glabrata cultured at Hospital A, Hospital B and Hospital C which yielded the majority of the isolates was nearly same. The frequency of genotypic heterogeneity among C. glabrata isolates was nearly same (48.5% versus 54.2%) among isolates from the two major hospitals which contributed 63% of the total 91 isolates as 16 STs were detected among 33 isolates from Hospital A while 13 STs were found among 24 isolates from Hospital B (Supplementary Table S1). The frequency of genotypic heterogeneity was slightly higher among isolates from Hospital C (8 STs among 12 isolates, 66.7%), Hospital D (5 STs among 7 isolates, 71.4%) and Hospital E (4 STs among 6 isolates, 66.7%) (Supplementary Table S1). No specific association was found between a specimen type and ST as the six most common and previously reported STs (ST46, ST3, ST7, ST55, ST15 and ST122) were found among urine, bloodstream and oral/respiratory isolates. Similarly, there was no specific association of an ST with resistance of the isolate to an antifungal drug.
The genetic association between the STs of C. glabrata isolates from Kuwait was determined by construction of an unrooted phylogenetic tree based on MLST data. The dendrogram (Figure 1) showed that 16 of 91 (17.6%) isolates were dispersed as unrelated singletons belonging to a single unique ST while the remaining 75 isolates typed into 13 clusters with the largest cluster (ST46) containing 33 of 75 (44%) of all cluster isolates. No obvious relationship was found between specific STs with patient’s nationality.
Figure 1. UPGMA tree based on MLST sequence data from 91 Candida glabrata isolates. The columns on left show MLST based sequence type (ST). Similarity is presented in percentages using the scale bar in the upper left corner. The columns from left to right include, patient number, age, sex, nationality, hospital name, isolate number, clinical source and MLST-based sequence type (ST).
The microevolution and nosocomial transmission of C. glabrata isolates among hospitalized patients in Kuwait was also investigated by using the minimal spanning tree (MST) algorithm of the BioNumerics (v.7.6.1 software) software and the MST is depicted in Figure 2. The microevolution and nosocomial transmission of infection was more evident among C. glabrata isolates collected during 2008 to 2011 as majority (74 of 91, 81.3%) of the tested isolates were from this period and 57 of 91 (62.6%) isolates were cultured at Hospital A and Hospital B which are also located close (within 500 meters) to each other. Four specific cases of microevolution and nosocomial transmission of C. glabrata infection were recognized at Hospital A (Cases 1 to 3) and Hospital B (Case 4) which occurred within the same time-frame (±1 year) and are summarized in Table 3.
Figure 2. Minimum spanning tree showing genotypic relationship among C. glabrata isolated in different hospitals in Kuwait. Each circle corresponds to a unique genotype, and lines between circles represent relative distance between isolates. The sizes of the circles correspond to the number of isolates of the same genotype (ST). Connecting lines correspond to the number of allele differences between genotypes, with a solid thick line connecting genotypes that differ in one locus, a solid thin line connecting genotypes that differ in two-three loci, a dashed line connecting genotypes that differ in four loci, and a dotted line connecting genotypes that differ in more than four loci.
Table 3. Microevolution and cross transmission of C. glabrata among hospitalized patients in two major hospitals in Kuwait.
Case 1 in Hospital A involved microevolution and nosocomial transmission of a C. glabrata isolate (Kw2450/10) obtained from wound and belonging to ST163 which evolved from an ST3 isolate obtained from urine (Kw2826/10) (Table 3). Case 2 involved a C. glabrata isolate (Kw1134/10) obtained from tracheal aspirate and belonging to ST165 which evolved from an ST161 urine isolate (Kw1073/10) while Case 3 involved two urine isolates, an ST160 isolate (Kw3025/09) which evolved from an ST162 isolate (Kw1659/10) (Table 3). Case 4 in Hospital B involved an ST157 C. glabrata urine isolate (Kw3368/09) which evolved from a bloodstream ST7 isolate (Kw2820/11) (Table 3). Interestingly, in all four cases, the two isolates contained identical sequence of the polymorphic ITS region of rDNA (Figure 3).
Figure 3. UPGMA tree based on internal transcribed spacer (ITS) region of rDNA sequence data from C. glabrata isolates. ITS sequences from four cases of C. glabrata (n = 10 isolates) with epidemiological linkage indicating microevolution and nosocomial transmission in the two hospitals in Kuwait.
The population structure of C. glabrata isolates from Kuwait was also investigated by determining their genetic relationship with isolates from other Asian countries in the central data library of the C. glabrata MLST database by using the MST and the data are depicted in Figure 4. Most of the isolates from Kuwait clustered with isolates from Iran while few isolates also clustered with C. glabrata isolates from India, Japan, and China (Figure 4).
Figure 4. Minimum spanning tree showing relationship of C. glabrata from Kuwait with 498 isolates from other Asian countries available from the MLST website as of June 19, 2023. Each circle corresponds to a unique genotype, and lines between circles represent relative distance between isolates. The sizes of the circles correspond to the number of isolates of the same genotype (ST). Connecting lines correspond to the number of allele differences between genotypes, with a solid thick line connecting genotypes that differ in one locus, a solid thin line connecting genotypes that differ in two-three loci, a dashed line connecting genotypes that differ in four loci, and a dotted line connecting genotypes that differ in more than four loci.
4. Discussion
The IFIs are increasing worldwide and are mostly caused by opportunistic Candida and other yeast species in nosocomial settings (1, 2, 5–8). Opportunistic Candida infections may be either endogenous or nosocomial in origin and this distinction is critical for effective intervention strategies. Molecular fingerprinting studies are crucial in determining the source of infection and nosocomial transmission of yeast pathogens in healthcare settings (25, 31, 34, 41, 42). Fingerprinting studies of Candida spp. isolates have mostly been carried out by employing MLST or MLMT-based methods (28–31, 34, 41–43). Although MLST schemes were devised for C. albicans (44) and C. glabrata (26) at nearly the same time, extensive studies have been carried out on clinical C. albicans isolates while only few studies have explored the population structure of C. glabrata (28, 31, 45–53).
C. glabrata is among the most frequently isolated NACS from clinical specimens in Kuwait (17, 18). A total of 91 C. glabrata isolates collected during 2007 to 2015 were selected for MLST analyses. Although these represented only ~8% of the total 1,149 isolates collected during 2007–2015, significantly higher number of isolates were analyzed from the years 2008 to 2011 as 74 of 362 (20.4%) isolates were included for MLST studies. MLST identified 28 STs among 91 isolates giving a ratio (number of STs to the number of isolates) of 3.25 which is nearly same as that reported from Australia (2.88, 18 STs among 52 isolates) (50), Iran (3.12, 16 STs among 50 isolates) (51) and Tanzania (3.6, 13 STs among 47 isolates) (52). On the contrary, a higher ratio of 5.3 (15 STs among 80 isolates) has been reported from Taiwan (46), 7.41 (31 STs among 230 isolates) (47) and 9.13 (22 STs among 201 isolates) (48) among C. glabrata isolates from the United States, 6.97 (30 STs among 209 isolates) from Republic of Korea (49) and 7.82 (17 STs among 133 isolates) (28) and 11.4 (35 STs among 411 isolates) (53) from China. The higher genetic diversity of C. glabrata from Kuwait, an Arabian Gulf country much smaller than Australia, Iran or Tanzania, is likely due to the large (nearly three times of Kuwaiti nationals) and highly diverse expatriate population (54, 55). This is also supported by the relatively higher frequency (16 of 91, 17.6%) of C. glabrata isolates from Kuwait represented by a single and unique ST. Similar values have also been reported from Australia (8 of 52, 15.3%) (50) and Iran (9 of 50, 18%) (51) while lower frequencies have been reported from Taiwan (6 of 80, 7.5%) (46), Republic of Korea (20 of 209, 9.6%) (49), China (23 of 411, 5.6%) (28) and (7 of 133, 5.2%) (53) and from the United States (16 of 230, 6.9%) (47) and (7 of 201, 3.5%) (48).
Although high genetic diversity was evident among 91 C. glabrata isolates collected over an extended period-of-time (2007 to 2015) from Kuwait, ST46 was the dominant genotype detected in 33 of 91 (36.3%) isolates and this genotype was present in nearly all hospitals and all major specimen types. Interestingly, ST46 has been detected previously in other countries, however, it was not found as the dominant genotype previously, including studies from other (Iran and Qatar) Middle Eastern countries (31, 51). On the contrary, ST7 is the predominant genotype in East Asian countries as this genotype was detected in 23 of 80 (28.7%) C. glabrata isolates from Taiwan (46), 100 of 209 (47.8%) isolates from Republic of Korea (49), and 273 of 411 (66.4%) (28) and 81 of 133 (60.9%) (53) isolates from China. The second most common genotype from Kuwait, ST3, has been reported as the most common genotype among C. glabrata isolates from Australia (8 of 51, 15.7%) (50) and United States (38 of 201, 18.1%) (48), it is the second most common genotype among isolates from Taiwan (13 of 80, 16.2%) (46), Republic of Korea (47 of 209, 22.5%) (49) and China (39 of 411, 9.5%) (28) and has also been detected among the isolates from Iran (4 of 50, 8%) (51). The ST3 is also a common genotype among C. glabrata isolates from Belgium, France, Germany, Spain and the United Kingdom (26, 50), reflecting the universal spreading of this genotype. Other common genotypes (ST55, ST15 and ST10) detected in this study also appear to have a worldwide distribution as they have been detected previously in several other countries, including, Australia, China, France, Japan, Iran, Qatar, Spain, Taiwan, Tanzania, United Kingdom and the United States (26, 28, 31, 47–53).
The distribution of three major genotypes (ST46, ST7 and ST55) among urine, oral/respiratory and blood specimens was not significantly different, however, ST3 was not detected among oral/respiratory specimens. There was no association of a particular ST with resistance of C. glabrata isolates in Kuwait to an antifungal drug. Only few previous studies have shown an association of specific STs with resistance of C. glabrata to antifungal drugs. Fluconazole resistant C. glabrata isolates from China were associated with ST7 and/or ST3 (28, 53) while ST18 isolates from Tanzania exhibited lower mean fluconazole susceptibility than isolates belonging to other STs (52).
Repeat C. glabrata isolates from the same patient obtained from the same site in Kuwait yielded the same ST. Only few studies have performed MLST on repeat isolates from the same patient. The only repeat isolate in the study from Australia belonged to the same genotype (ST8) as the first isolate (50). Lin et al. (46), analyzed multiple isolates from 25 patients and their data showed that repeat isolates from 19 patients belonged to the same ST or very closely related genotype. However, multiple isolates from 6 patients belonged to distinct STs (46). Our study also showed that bloodstream C. glabrata isolates from three patients belonged to the same ST as non-blood isolates obtained from the catheter tip (patient P79) or mucosal sites (oral cavity and rectum) (patients P57 and P80, respectively). These findings suggest translocation of C. glabrata from mucosal sites to bloodstream causing candidemia in these three patients. Similar findings have also been recorded recently among patients colonized with multidrug-resistant Candida auris who later developed candidemia (56, 57). A total of 5 isolates were analyzed from a patient (P60) with persistent bloodstream infection and all isolates belonged to ST46. A recent study showed that multiple C. glabrata colonies from candidemia patients exhibited heterogenous susceptibility results but were indistinguishable by MLST as they all belonged to ST3 (58). Unfortunately, whole genome sequencing (WGS) was not performed on C. glabrata isolates described by Knoll et al. (58) as well as on the ST46 isolates from P60 in our study. However, a more recent study employing WGS has shown that C. glabrata isolates from candidemia patients contain mixed populations of clonal but genetically diverse strains which cannot be differentiated by MLST alone (59). These findings suggest that MLST likely records events which happened several months/years ago while WGS can detect changes occurring in the same patient over days/weeks.
An interesting observation of our study is the detection of microevolution and nosocomial transmission of C. glabrata in the two major hospitals which yielded most of the isolates analyzed in this study. Three specific cases of microevolution and nosocomial transmission of C. glabrata were detected in Hospital A and one in Hospital B. In each case, the two isolates contained identical sequence of the ITS region of rDNA which, unlike other Candida spp., is often variable among clinical C. glabrata sensu stricto isolates (36, 43, 60). Although clinical details of patients, disease severity and outcome were not available, there was no apparent relationship among microevolution-related C. glabrata isolates with respect to the age or gender of patients. Nosocomial transmission of C. glabrata, initially indicated by PLST (27), was recently demonstrated for cluster P and cluster N isolates in a recent study by including additional analyses of another highly polymorphic (CgMT-C) locus (25). Although not specifically stated, microevolution of C. glabrata isolate belonging to ST182 from an isolate of ST26 and its nosocomial transmission was also evident among the two strains (L255 and L256) isolated from the two patients in LH Hospital in a recent study from China (53).
Our study has a few limitations. (1) Only 91 of the total 1,149 (7.9%) C. glabrata isolates collected during 2007 to 2015 in Kuwait were analyzed by MLST. These isolates were selected from >400 of C. glabrata isolates that were recently characterized in detail by various molecular methods including PCR-sequencing of rDNA. (2) The distribution of the isolates for some years and from some hospitals was too small to draw any meaningful conclusions for these settings. (3) The microevolution and nosocomial transmission of C. glabrata suggested by MLST analyses was not confirmed by WGS studies. (4) The clinical details of patients, disease severity and outcome were not available. Since MLST is labor intensive and time consuming, we intend to extend these observations by performing WGS on C. glabrata isolates collected more recently.
5. Conclusion
We have performed molecular fingerprinting of C. glabrata isolates by multilocus sequence typing for the first time in Kuwait and the Arabian Peninsula. Our results have identified a dominant genotype (ST46) among C. glabrata isolates in Kuwait and have also suggested that C. glabrata undergoes microevolution in the hospital environment and is nosocomially transmitted to other susceptible patients. Further investigations by WGS are needed to confirm these findings. If WGS confirms microevolution and nosocomial transmission of infection, proper infection control practices during routine procedures on C. glabrata-infected patients will greatly help in preventing transmission of this pathogen to other hospitalized patients.
Data availability statement
All relevant data described in the paper are present within the manuscript and in the supplementary material. The rDNA sequence data have been submitted to GenBank under accession no. FN652301, LS398116, LT837722, LT837724 and OR226607-OR226612 and are freely available. The sequence data for different alleles comprising new MLST STs reported in this study have also been submitted to the C. glabrata MLST database (https://pubmlst.org/bigsdb?db=pubmlst_cglabrata_seqdef) which are also freely available.
Ethics statement
The studies involving humans were approved by Ethics Committee of the Health Sciences Center, Kuwait University. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired as a by-product of routine patient care. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
MA, SA, NA-S, and ZK conceptualized the research and wrote the manuscript. MA, SA, and ZK supervised the study, designed the experiments, and interpreted the data. All authors contributed to the article and approved the submitted version.
Funding
The study was supported by the Department of Microbiology, College of Medicine, Kuwait University and Research Core Facility grant SRUL02/13.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2023.1242622/full#supplementary-material
Footnotes
References
1. Brown, GD, Denning, DW, Gow, NA, Levitz, SM, Netea, MG, and White, TC. Hidden killers: human fungal infections. Sci Transl Med. (2012) 4:165rv13. doi: 10.1126/scitranslmed.3004404
2. McCarty, TP, and Pappas, PG. Invasive candidiasis. Infect Dis Clin N Am. (2016) 30:103–24. doi: 10.1016/j.idc.2015.10.013
3. Leroy, O, Gangneux, JP, Montravers, P, Mira, JP, Gouin, F, Sollet, JP, et al. 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. (2009) 37:1612–8. doi: 10.1097/CCM.0b013e31819efac0
4. Pfaller, MA, Andes, DR, Diekema, DJ, Horn, DL, Reboli, AC, Rotstein, C, et al. Epidemiology and outcomes of invasive candidiasis due to non-albicans species of Candida in 2496 patients: data from the prospective antifungal therapy (PATH) registry 2004–2008. PLoS One. (2014) 9:e101510. doi: 10.1371/journal.pone.0101510
5. Lamoth, F, Lockhart, SR, Berkow, EL, and Calandra, T. Changes in the epidemiological landscape of invasive candidiasis. J Antimicrob Chemother. (2018) 73:i4–i135. doi: 10.1093/jac/dkx444
6. Jenks, JD, Cornely, OA, Chen, SC, Thompson, GR, and Hoenigl, M. Breakthrough invasive fungal infections: who is at risk? Mycoses. (2020) 63:1021–32. doi: 10.1111/myc.13148
7. Ahmad, S, and Alfouzan, W. Candida auris: epidemiology, diagnosis, pathogenesis, antifungal susceptibility, and infection control measures to combat the spread of infections in healthcare facilities. Microorganisms. (2021) 9:807. doi: 10.3390/microorganisms9040807
8. Wiederhold, NP . Emerging fungal infections: new species, new names, and antifungal resistance. Clin Chem. (2021) 68:83–90. doi: 10.1093/clinchem/hvab217
9. Ostrosky-Zeichner, L . Candida glabrata and FKS mutations: witnessing the emergence of the true multidrug-resistant Candida. Clin Infect Dis. (2013) 56:1733–4. doi: 10.1093/cid/cit140
10. Al-Baqsami, Z, Ahmad, S, and Khan, Z. Antifungal drug susceptibility, molecular basis of resistance to echinocandins and molecular epidemiology of fluconazole resistance among clinical Candida glabrata isolates in Kuwait. Sci Rep. (2020) 10:6238. doi: 10.1038/s41598-020-63240-z
11. Ahmad, S, Joseph, L, Parker, JE, Asadzadeh, M, Kelly, SL, Meis, JF, et al. ERG6 and ERG2 are major targets conferring reduced susceptibility to amphotericin B in clinical Candida glabrata isolates in Kuwait. Antimicrob Agents Chemother. (2019) 63:e01900–18. doi: 10.1128/AAC.01900-18
12. Khan, Z, Ahmad, S, Mokadda, E, Meis, JF, Joseph, L, Abdullah, A, et al. Development of echinocandin resistance in Candida tropicalis following short-term exposure to caspofungin for empiric therapy. Antimicrob Agents Chemother. (2018) 62:e01926–17. doi: 10.1128/AAC.01926-17
13. Pfaller, MA, Diekema, DJ, Turnidge, JD, Castanheira, M, and Jones, RN. Twenty years of the SENTRY antifungal surveillance program: results for Candida species from 1997–2016. Open Forum Infect Dis. (2019) 6:S79–94. doi: 10.1093/ofid/ofy358
14. Alobaid, K, Ahmad, S, Asadzadeh, M, Mokaddas, E, Al-Sweih, N, Albenwan, K, et al. Epidemiology of candidemia in Kuwait: a nationwide, population-based study. J Fungi. (2021) 7:673. doi: 10.3390/jof7080673
15. Greenberg, RG, Kandefer, S, Do, BT, Smith, PB, Stoll, BJ, Bell, EF, et al. Late-onset sepsis in extremely premature infants: 2000–2011. Pediatr Infect Dis J. (2017) 36:774–9. doi: 10.1097/INF.0000000000001570
16. Basmaciyan, L, Bon, F, Paradis, T, Lapaquette, P, and Dalle, F. Candida albicans interactions with the host: crossing the intestinal epithelial barrier. Tissue Barriers. (2019) 7:1612661. doi: 10.1080/21688370.2019.1612661
17. Khan, Z, Ahmad, S, Al-Sweih, N, Mokaddas, E, Al-Banwan, K, Alfouzan, W, et al. Changing trends in epidemiology and antifungal susceptibility patterns of six bloodstream Candida species isolates over a 12-year period in Kuwait. PLoS One. (2019) 14:e0216250. doi: 10.1371/journal.pone.0216250
18. Khan, Z, Ahmad, S, Al-Sweih, N, Mokaddas, E, Al-Banwan, K, Alfouzan, W, et al. Increasing trends of reduced susceptibility to antifungal drugs among clinical Candida glabrata isolates in Kuwait. Microb Drug Resist. (2020) 26:982–90. doi: 10.1089/mdr.2019.0437
19. Kidd, SE, Abdolrasouli, A, and Hagen, F. Fungal nomenclature: managing change is the name of the game. Open Forum Infect Dis. (2023) 10:ofac559. doi: 10.1093/ofid/ofad395
20. Pfaller, MA . Antifungal drug resistance: mechanisms, epidemiology, and consequences for treatment. Am J Med. (2012) 125:S3–S13. doi: 10.1016/j.amjmed.2011.11.001
21. Logan, C, Martin-Loeches, I, and Bicanic, T. Invasive candidiasis in critical care: challenges and future directions. Intensive Care Med. (2020) 46:2001–14. doi: 10.1007/s00134-020-06240-x
22. Hassan, Y, Chew, SY, and Than, LTL. Candida glabrata: pathogenicity and resistance mechanisms for adaptation and survival. J Fungi. (2021) 7:667. doi: 10.3390/jof7080667
23. Bouza, E, and Muñoz, P. Epidemiology of candidemia in intensive care units. Int J Antimicrob Agents. (2008) 32:S87–91. doi: 10.1016/S0924-8579(08)70006-2
24. Gabaldón, T, Gómez-Molero, E, and Bader, O. Molecular typing of Candida glabrata. Mycopathologia. (2020) 185:755–64. doi: 10.1007/s11046-019-00388-x
25. Katiyar, S, and Edlind, T. New locus for Candida glabrata sequence-based strain typing provides evidence for nosocomial transmission. J Clin Microbiol. (2021) 59:e02933–20. doi: 10.1128/JCM.02933-20
26. Dodgson, AR, Pujol, C, Denning, DW, Soll, DR, and Fox, AJ. Multilocus sequence typing of Candida glabrata reveals geographically enriched clades. J Clin Microbiol. (2003) 41:5709–17. doi: 10.1128/JCM.41.12.5709-5717.2003
27. Katiyar, S, Shiffrin, E, Shelton, C, Healey, K, Vermitsky, JP, and Edlind, T. Evaluation of polymorphic locus sequence typing for Candida glabrata epidemiology. J Clin Microbiol. (2016) 54:1042–50. doi: 10.1128/JCM.03106-15
28. Hou, X, Xiao, M, Chen, SC-A, Kong, F, Wang, H, Chu, Y-Z, et al. Molecular epidemiology and antifungal susceptibility of Candida glabrata in China (August 2009 to July 2014): a multi-center study. Front Microbiol. (2017) 8:880. doi: 10.3389/fmicb.2017.00880
29. Arastehfar, A, Daneshnia, F, Zomorodian, K, Najafzadeh, MJ, Khodavaisy, S, Zarrinfar, H, et al. Low level of antifungal resistance in Iranian isolates of Candida glabrata recovered from blood samples in a multicenter study from 2015 to 2018 and potential prognostic values of genotyping and sequencing of PDR1. Antimicrob Agents Chemother. (2019) 63:e02503–18. doi: 10.1128/AAC.02503-18
30. Arastehfar, A, Daneshnia, F, Salehi, MR, Yaşar, M, Hoşbul, T, Ilkit, M, et al. Low level of antifungal resistance of Candida glabrata blood isolates in Turkey: fluconazole minimum inhibitory concentration and FKS mutations can predict therapeutic failure. Mycoses. (2020) 63:911–20. doi: 10.1111/myc.13104
31. Arastehfar, A, Marcet-Houben, M, Daneshnia, F, Taj-Aldeen, SJ, Batra, D, Lockhart, SR, et al. Comparative genomic analysis of clinical Candida glabrata isolates identifies multiple polymorphic loci that can improve existing multilocus sequence typing strategy. Stud Mycol. (2021) 100:100133. doi: 10.1016/j.simyco.2021.100133
32. Khan, Z, Ahmad, S, Al-Sweih, N, Joseph, L, Alfouzan, W, and Asadzadeh, M. Increasing prevalence, molecular characterization and antifungal drug susceptibility of serial Candida auris isolates in Kuwait. PLoS One. (2018) 13:e0195743. doi: 10.1371/journal.pone.0208951
33. Jamal, WY, Ahmad, S, Khan, ZU, and Rotimi, VO. Comparative evaluation of two matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) systems for the identification of clinically significant yeasts. Int J Infect Dis. (2014) 26:167–70. doi: 10.1016/j.ijid.2014.05.031
34. Al-Obaid, K, Asadzadeh, M, Ahmad, S, and Khan, Z. Population structure and molecular genetic characterization of clinical Candida tropicalis isolates from a tertiary-care hospital in Kuwait reveal infections with unique strains. PLoS One. (2017, 2017) 12:e0182292. doi: 10.1371/journal.pone.0182292
35. Asadzadeh, M, Ahmad, S, Hagen, F, Meis, JF, Al-Sweih, N, and Khan, Z. Simple, low-cost detection of Candida parapsilosis complex isolates and molecular fingerprinting of Candida orthopsilosis strains in Kuwait by ITS region sequencing and amplified fragment length polymorphism analysis. PLoS One. (2015) 10:e0142880. doi: 10.1371/journal.pone.0142880
36. Asadzadeh, M, Alanazi, AF, Ahmad, S, Al-Sweih, N, and Khan, Z. Lack of detection of Candida nivariensis and Candida bracarensis among 440 clinical Candida glabrata sensu lato isolates in Kuwait. PLoS One. (2019) 14:e0223920. doi: 10.1371/journal.pone.0223920
37. Khan, ZU, Ahmad, S, Hagen, F, Fell, JW, Kowshik, T, Chandy, R, et al. Cryptococcus randhawai sp. nov, a novel anamorphic basidiomycetous yeast isolated from tree trunk hollow of Ficus religiosa (peepal tree) from New Delhi, India. Antonie Van Leeuwenhoek. (2010) 97:253–9. doi: 10.1007/s10482-009-9406-8
38. Asadzadeh, M, Al-Sweih, NA, Ahmad, S, and Khan, ZU. Antifungal susceptibility of clinical Candida parapsilosis isolates in Kuwait. Mycoses. (2008) 51:318–23. doi: 10.1111/j.1439-0507.2008.01492.x
39. Clinical and Laboratory Standards Institute . Fourth informational supplement M27-S4 In: Reference method for broth dilution antifungal susceptibility testing of yeasts. 4th ed. Wayne, PA: Clinical Laboratory Standards Institute (2012)
40. Bourgeois, N, Laurens, C, Bertout, S, Balard, Y, Krasteva, D, Rispail, P, et al. Assessment of caspofungin susceptibility of Candida glabrata by the Etest, CLSI, and EUCAST methods, and detection of FKS1 and FKS2 mutations. Eur J Clin Microbiol Infect Dis. (2014) 33:247–52. doi: 10.1007/s10096-014-2069-z
41. Asadzadeh, M, Ahmad, S, Al-Sweih, N, and Khan, Z. Molecular fingerprinting studies do not support intrahospital transmission of Candida albicans among candidemia patients in Kuwait. Front Microbiol. (2017) 8:247. doi: 10.3389/fmicb.2017.00247
42. Asadzadeh, M, Ahmad, S, Al-Sweih, N, Hagen, F, Meis, JF, and Khan, Z. High-resolution fingerprinting of Candida parapsilosis isolates suggests persistence and transmission of infections among neonatal intensive care unit patients in Kuwait. Sci Rep. (2019) 9:1340. doi: 10.1038/s41598-018-37855-2
43. Asadzadeh, M, Ahmad, S, Al-Sweih, N, and Khan, Z. Population structure and molecular genetic characterization of 5-flucytosine-susceptible and-resistant clinical Candida dubliniensis isolates from Kuwait. PLoS One. (2017) 12:e0175269. doi: 10.1371/journal.pone.0175269
44. Tavanti, A, Gow, NA, Senesi, S, Maiden, MC, and Odds, FC. Optimization and validation of multilocus sequence typing for Candida albicans. J Clin Microbiol. (2003) 41:3765–76. doi: 10.1128/JCM.41.8.3765-3776.2003
45. McManus, BA, and Coleman, DC. Molecular epidemiology, phylogeny and evolution of Candida albicans. Infect Genet Evol. (2014) 21:166–78. doi: 10.1016/j.meegid.2013.11.008
46. Lin, CY, Chen, YC, Lo, HJ, Chen, KW, and Li, SY. Assessment of Candida glabrata strain relatedness by pulsed-field gel electrophoresis and multilocus sequence typing. J Clin Microbiol. (2007) 45:2452–9. doi: 10.1128/JCM.00699-07
47. Lott, TJ, Frade, JP, and Lockhart, SR. Multilocus sequence type analysis reveals both clonality and recombination in populations of Candida glabrata bloodstream isolates from U.S. surveillance studies. Eukaryot Cell. (2010) 9:619–25. doi: 10.1128/EC.00002-10
48. Lott, TJ, Frade, JP, Lyon, GM, Iqbal, N, and Lockhart, SR. Bloodstream and non-invasive isolates of Candida glabrata have similar population structures and fluconazole susceptibilities. Med Mycol. (2012) 50:136–42. doi: 10.3109/13693786.2011.592153
49. Byun, SA, Won, EJ, Kim, MN, Lee, WG, Lee, K, Lee, HS, et al. Multilocus sequence typing (MLST) genotypes of Candida glabrata bloodstream isolates in Korea: association with antifungal resistance, mutations in mismatch repair gene (Msh 2), and clinical outcomes. Front Microbiol. (2018) 13:1523. doi: 10.3389/fmicb.2018.01523
50. Biswas, C, Marcelino, VR, Van Hal, S, Halliday, C, Martinez, E, Wang, Q, et al. Whole genome sequencing of Australian Candida glabrata isolates reveals genetic diversity and novel sequence types. Front Microbiol. (2018) 9:2946. doi: 10.3389/fmicb.2018.02946
51. Amanloo, S, Shams-Ghahfarokhi, M, Ghahri, M, and Razzaghi-Abyaneh, M. Genotyping of clinical isolates of Candida glabrata from Iran by multilocus sequence typing and determination of population structure and drug resistance profile. Med Mycol. (2018) 56:207–15. doi: 10.1093/mmy/myx030
52. Mushi, MF, Gross, U, Mshana, SE, and Bader, O. High diversity of Candida glabrata in a tertiary hospital-Mwanza, Tanzania. Med Mycol. (2019) 57:914–7. doi: 10.1093/mmy/myy151
53. Chen, Y, Wu, Y, Lulou, K, Yao, D, and Ying, C. Multilocus sequence typing and antifungal susceptibility of vaginal and non-vaginal Candida glabrata isolates from China. Front Microbiol. (2022) 13:808890. doi: 10.3389/fmicb.2022.1085000
54. Ahmad, S, Mokaddas, E, and Al-Mutairi, NM. Prevalence of tuberculosis and multidrug resistant tuberculosis in the Middle East Region. Expert Rev Anti-Infect Ther. (2018) 16:709–21. doi: 10.1080/14787210.2018.1519393
55. Al-Awadhi, M, Ahmad, S, and Iqbal, J. Current status and the epidemiology of malaria in the Middle East Region and beyond. Microorganisms. (2021) 9:338. doi: 10.3390/microorganisms9020338
56. Al-Obaid, I, Asadzadeh, M, Ahmad, S, Alobaid, K, Alfouzan, W, Bafna, R, et al. Fatal breakthrough candidemia in an immunocompromised patient in Kuwait due to Candida auris exhibiting reduced susceptibility to echinocandins and carrying a novel mutation in hotspot-1 of FKS1. J Fungi. (2022) 8:267. doi: 10.3390/jof8030267
57. Asadzadeh, M, Mokaddas, E, Ahmad, S, Abdullah, AA, de Groot, T, Meis, JF, et al. Molecular characterisation of Candida auris isolates from immunocompromised patients in a tertiary-care hospital in Kuwait reveals a novel mutation in FKS1 conferring reduced susceptibility to echinocandins. Mycoses. (2022) 65:331–43. doi: 10.1111/myc.13419
58. Knoll, MA, Lackner, N, Ulmer, H, Samardzic, E, Steinmann, J, Krause, R, et al. Multiple colony antifungal susceptibility testing detects polyresistance in clinical Candida cultures: a European Confederation of Medical Mycology excellence centers study. Clin Microbiol Infect. (2022) 28:e1–1288.e7. doi: 10.1016/j.cmi.2022.04.014
59. Badrane, H, Cheng, S, Dupont, C, Hao, B, Driscoll, E, Morder, K, et al. Genotypic diversity and unrecognized antifungal resistance among populations of Candida glabrata from positive blood cultures. Res Sq. (2023):rs3.rs-2706400. doi: 10.21203/rs.3.rs-2706400/v1
Keywords: Candida glabrata , MLST, microevolution, nosocomial transmission, fingerprinting
Citation: Asadzadeh M, Ahmad S, Al-Sweih N and Khan Z (2023) Molecular fingerprinting by multi-locus sequence typing identifies microevolution and nosocomial transmission of Candida glabrata in Kuwait. Front. Public Health. 11:1242622. doi: 10.3389/fpubh.2023.1242622
Edited by:
Jian Wu, Suzhou Municipal Hospital, ChinaReviewed by:
Aylin Döğen, Mersin University, TürkiyeHéctor Armando Rincón León, Mexican Social Security Institute (IMSS), Mexico
Copyright © 2023 Asadzadeh, Ahmad, Al-Sweih and Khan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Suhail Ahmad, suhail.ahmad@ku.edu.kw