Isolation of Candida
The distribution of various types of Candida strains identified in the hospital from 2016 to 2019 was analyzed. As shown in Fig. 1, a total of 2872 strains of Candida were identified. After the Cochran Q test, the composition ratio (i.e., the percentage of all detected C. tropicalis species) of various Candida species in each year showed statistical significance (Q = 1101.094 - 1904.945; all P < 0.001). For each year, C. albicans was the major strain, accounting for more than 60.00% of the total. C. tropicalis ranked third; there was no significant variation in the proportion of isolates between 2016, 2018 and 2019. Additionally, the composition ratio of each kind of Candida was compared between these four years. Our results showed statistical significance only for C. albicans (χ2 = 12.620; P = 0.006) and C. tropicalis (χ2 = 20.410; P < 0.001) (Additional file 1:Fig. S1).
Sources of C. tropicalis samples
The 319 strains of C. tropicalis collected from 2016-2019 were classified according to the sample sources, and the composition ratios were also analyzed. Our results showed that significant differences were observed in the proportion of samples from which C. tropicalis was detected in each year (Q = 123.949 - 194.898; P < 0.001) (Fig. 2). Moreover, our results from the trend test showed that, except for the proportion in the stool samples (χ2 = 3.550; P = 0.060), the proportions in the airway secretions and urine samples increased over the years, but the proportions decreased for the other sample types year by year (all P < 0.050) (Additional file 2:Fig. S2).
Department distribution of C. tropicalis
The department sources of C. tropicalis from 2016-2019 were analyzed. Our results showed that , there were significant differences in the distribution of C. tropicalis in these departments from 2016 to 2019 (Q = 40.746 - 92.691; all P < 0.001). The departments with relatively high detection rates were the Departments of Respiratory Medicine, ICU, and Geriatrics (Fig. 3). However, in comparing the C. tropicalis amount detected in each department in each of the four years, statistically significant differences were observed in the amount of C. tropicalis in the Department of Geriatrics (χ2 = 8.623; P = 0.035) and the ICU (χ2 = 27.148; P < 0.001) (Additional file 3:Fig. S3).
Susceptibility of C. tropicalis to antifungal drugs
The susceptibility to azole antifungal agents including fluconazole, itraconazole and voriconazole was analyzed in the 319 samples of C. tropicalis collected from 2016 to 2019. Our results showed that C. tropicalis had high resistance rates to fluconazole, itraconazole, and voriconazole and even exhibited cross-resistance. In 2018, the resistance rate of C. tropicalis to fluconazole reached 39.51% (32/81) and thus this year had the most resistant C. tropicalis. There were no statistically significant differences in the resistance rates of C. tropicalis to azole antifungal drugs over these four years (χ2 = 1.156; P = 0.979). However, statistical analysis found that there were significant differences in the resistance rates of C. tropicalis to fluconazole (χ2 = 10.455; P = 0.015) and voriconazole (χ2 = 9.154; P = 0.027) from 2016 to 2019 (Table 1).
Table 1 The resistance rate of C. tropicalis to antifungal drugs from 2016 to 2019(%)
Antifungal drug
|
2016(n=98)
|
2017(n=59)
|
2018(n=81)
|
2019(n=81)
|
Fluconazole
|
18.37(18/98)
|
23.73(14/59)
|
39.51(32/81)
|
28.40(23/81)
|
Itraconazole
|
18.37(18/98)
|
28.81(17/59)
|
30.86(25/81)
|
23.46(19/81)
|
Voriconazole
|
17.35(17/98)
|
23.73(14/59)
|
37.04(30/81)
|
25.93(21/81)
|
n is the number of strains.
Relative mRNA expression levels of ERG11 and UPC2 in C. tropicalis
To further study the resistance-related genes of C. tropicalis, the mRNA expression levels of ERG11 and UPC2 were detected with quantitative real-time PCR in 20 fluconazole-susceptible and 30 fluconazole-resistant strains extracted from 319 strains. Our results showed that the relative mRNA expression levels of ERG11 and UPC2 genes in C. tropicalis from the fluconazole-resistant group (27/30 strains fully resistant to fluconazole, itraconazole, and voriconazole) and fluconazole- susceptible group were normally distributed. In fluconazole-resistant group, 20 strains had ERG11 overexpression and 18 strains had upregulated UPC2 expression. However, the expression level of ERG11 and UPC2 had no significant difference in the fluconazole-susceptible group. Based on the independent samples t-tests, the relative mRNA expression level of ERG11 in the drug-resistant group was 1.579 ± 0.896, while the relative mRNA expression of ERG11 in the susceptible group was 0.483 ± 0.259, and the difference was statistically significant (t = 4.511; P <0.001) (Fig. 4A). On the other hand, the relative mRNA expression level of UPC2 in the resistant group was 1.400±0.919, while the relative mRNA expression level of UPC2 in the susceptible group was 0.448 ± 0.272, with a statistically significant difference (t = 3.970; P <0.001) (Fig. 4B). These results suggest that the resistance of C. tropicalis to fluconazole is related to the expression levels of ERG11 and UPC2.
Correlation analysis of UPC2 and ERG11 mRNA expression in C. tropicalis
The mRNA expression levels of resistance-related genes (i.e., ERG11 and UPC2) in C. tropicalis were detected, and the correlation between the gene expression levels was further analyzed. Our results from Spearman correlation analysis showed that there was no linear correlation between the expression levels of UPC2 and ERG11 in the susceptible group (r = - 0.074; P = 0.757) (Fig. 5A). However, the UPC2 and ERG11 expression levels were positively correlated in the drug-resistance group (r = 0.571; P = 0.001) (Fig. 5B).