DOI QR코드

DOI QR Code

Multi-locus sequence typing을 이용한 한국에서 분리한 Candida glabrata 임상균주의 유전자 유형 분석

Genetic Variations of Candida glabrata Clinical Isolates from Korea using Multi-locus Sequence Typing

  • 강민지 (충남대학교 의과대학 미생물학교실) ;
  • 이경은 (부산가톨릭대학교 보건과학대학 임상병리학과) ;
  • 진현우 (부산가톨릭대학교 보건과학대학 임상병리학과)
  • Kang, Min Ji (Department of Microbiology, College of Medicine, Chungnam National University) ;
  • Lee, Kyung Eun (Department of Clinical Laboratory Science, College of Health Science, Catholic University of Pusan) ;
  • Jin, Hyunwoo (Department of Clinical Laboratory Science, College of Health Science, Catholic University of Pusan)
  • 투고 : 2019.10.23
  • 심사 : 2020.01.31
  • 발행 : 2020.02.28

초록

Candida albicans가 칸디다 혈증의 주요 원인 균으로 알려져 있지만, 최근에는 non-albicans Candida (NAC)에 의한 심각한 감염이 증가하고 있다. NAC 중에서 C. glabrata는 두 번째로 병원성을 나타내는 원인 균이지만 C. glabrata의 구조, 역학 등의 기본적인 연구는 거의 없는 실정이다. 따라서 본 연구에서는 다양한 유형의 임상 표본으로부터 분리된 총 102개의 C. glabrata균주로 multi-locus sequence typing (MLST)을 수행하였다. FKS, LEU2, NMT1, TRP1, UGP1 및 URA3을 포함한 6개의 하우스키핑 유전자를 증폭하여 염기서열을 분석하였다. 총 3,345개의 DNA 염기서열 중 49개의 가변 염기서열 부위가 발견되었으며, 그 결과 102개의 균주에서 총 12개의 상이한 서열 유형을 확인하였고, 알려지지 않은 서열 유형(USTs) 중 UST1이 가장 많이 나타났다. 이 연구의 결과는 국내 C. glabrata 항생제 처방에 도움이 될 것으로 사료되며, 한국에서 유행하는 C. glabrata에 대한 추가 연구를 위한 기본 역학자료로 사용될 것으로 기대한다.

Although Candida albicans is the major fungal pathogen of candidemia, severe infections by non-albicans Candida (NAC) spp. have been increasing in recent years. Among NAC spp., C. glabrata has emerged as the second most common pathogen. However, few studies have been conducted to investigate its structure, epidemiology, and basic biology. In the present study, multi-locus sequence typing (MLST) was performed with a total of 102 C. glabrata clinical isolates that were isolated from various types of clinical specimen. For MLST, six housekeeping genes-FKS, LEU2, NMT1, TRP1, UGP1, and URA3-were amplified and sequenced. The results were analyzed using the C. glabrata database. Out of a total of 3,345 base-pair DNA sequences, 49 variable nucleotide sites were found, and the results showed that 12 different sequence types (ST) were identified from the 102 clinical isolates. The data also demonstrated that the undetermined ST1 was the most predominant ST in Korea. Further, seven undetermined STs (UST) containing UST2-8 were classified at specific loci. The data from this study may provide a fundamental database for further studies on C. glabrata, including its epidemiology and evolution. The data may also contribute to the development of novel antifungal agents and diagnostic tests.

키워드

Introduction

Candida spp. belong to the normal flora of the vaginal tract, the gastrointestinal tract, and the oral cavity in human [8,18]. However, rarely, serious infections, ranging from mucosal infections to systemic infection have been caused by Candida spp. [8,18]. Fungal infections caused by Candida spp. have increased significantly, especially in acquired immune deficiency syndrome (AIDS) and immunocompromised individuals, including intensive care and, elderly patients [3,8]. Also, candidemia is associated with a high mortality rate approximately 30 to 40% in hospitalized patients and is difficult to treat, thus increasing the cost of medical care [3,8].

C. glabrata has emerged as the second or third most common Candida pathogen after C. albicans in the United States, depending on the site [3, 7, 21]. Despite its increased prevalence, there have been relatively few studies on the population structure, epidemiology, and basic biology of C. glabrata compared to those conducted on other Candida spp.[6-8].

As mentioned above, C. albicans has considered to be the major fungal pathogen of candidemia in the past [6,23]. However, as the number of severe infections caused by non-albicans Candida spp. (NAC) have increased, studies have shifted from C. albicans to NAC such as C. glabrata in recent years [11,21]. Furthermore, since C. glabrata infections are often resistant to azole antifungal drugs, especially fluconazole, it is important to distinguish NAC from C. albicans to ensure the appropriate antifungal therapy and clinical management [8, 20, 21]. Thus, the discrimination of subtypes in these species are required for investigating their epidemiology and evolutionary biology [15, 23, 24].

In recent years, there has been substantial progress in the development of several molecular methods for typing subspecies and strains of fungi [26]. For instance, pulsed-field gel electrophoresis (PFGE) compares total DNA band patterns with or without restriction enzyme digestion, while multilocus variable-number tandem-repeat (VNTR) analysis examines length variations in six to nine PCR-amplified loci that contain polymorphic tandem repeats. Further, the random amplification of polymorphic DNA compares banding patterns following PCR with a nonspecific primer. Finally, multilocus enzyme electrophoresis, studies the different electrophoretic mobility of multiple core metabolic enzymes. These four approaches have some limitations, such as a lower reproducibility and portability [10,22], and the results obtained in different laboratories are difficult to compare [2,13].

Among these genotyping methods, multilocus sequence typing (MLST) is a useful tool to assign single nucleotide polymorphisms as allele numbers, which are stored in a database on line (PubMLST) and determine the differences from between closely related isolates by their geographical origins, sources, and other properties [23]. Also, it is possible that database are accessed by laboratories worldwide [13,24].

In the present study, MLST was performed with a total of 102 C. glabrata clinical isolates from various clinical specimens such as blood, urine, and other body fluids in Korea. The results were analyzed by using the C. glabrata MLST database (http://pubmlst.org/cglabrata/). The aim of the study is to discriminate STs in the same C. glabrata spp. by using common MLST and investigate the most prevalent ST from the C. glabrata in Korea.

Materials and Methods

Clinical strains

A total of 102 C. glabrata clinical isolates were provided from Korean Culture Collection of Medical Fungi (KCMF) and those isolates were collected from tertiary hospitals in Korea. Clinical isolates were isolated from a wide variety of clinical samples, including blood, catheterized urine, bile and other body fluids (Table 1).

Table 1. Details of diverse yeast isolates used in this study

SMGHBM_2020_v30n2_122_t0001.png 이미지

1)Others: ascitic fluid, joint fluid, pleural fluid, tissue etc.

2)No clinical information.

Genomic DNA extraction from fungal isolates

Genomic DNA of C. glabrata clinical isolates was extracted using a I-genomic BYF DNA Extraction Mini Kit (iNtRON Inc., Seongnam, Korea) according to the manufacturer's instructions [5]. The concentration and purity of the genomic DNA were checked by 260/280 optical density using a Nanodrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The extracted genomic DNA was stored at 4℃ until use.

Polymerase chain reaction and sequence analysis of the fungal ITS region for precise identification of C. glabrata clinical isolates

The fungal internal transcribed spacer (ITS) region, the conserved region between the 18S and 28S ribosomal RNA (rRNA) was amplified and sequenced using each primer pairs (Table 2). Target amplification was carried out in 20 μl reaction mixture containing 10 μl Prime Taq Primix (Genet Bio Inc., Daejeon, Korea), 5 μl of distilled ultra-pure water, 1 μl of each primer (10 pmol/μl), and 3 μl of genomic DNA template. The PCR condition was: an initial denaturation at 94℃ for 1 min, 30 cycles including subsequent denaturation at 94℃ for 30 sec, annealing at 57℃ for 30 sec and extension at 72℃ for 45 sec followed by final extension at 72℃ for 7 min and holding at 4℃. The amplified products were visualized by gel electrophoresis to confirm the presence of desired product.

Table 2. Sequences of primer pairs for fungal species identification

SMGHBM_2020_v30n2_122_t0002.png 이미지

The resulting amplicon was purified and sequenced by Macrogen Inc. (Daejeon, Korea). All sequences with lowquality bases in the chromatogram were re-sequenced for the high-quality results. The obtained sequences were aligned with reference sequences in the Genebank database using the basic local alignment search tool (BLAST) at the National Center for Biotechnology Information (NCBI), and percent homology scores were generated to precise identification of  C. glabrata clinical isolates.

MLST analysis for identifying sequence type of C. glabrata clinical isolates

Table 3 shows primers for the amplification and sequence analysis of C. glabrata six housekeeping gene fragments including FKS, LEU2, NMT1, TRP1, UGP1, and URA3. For the PCR amplification, 20 μl of final mixture contained 10 μl of Prime Taq Premix, 5 μl of distilled ultra-pure water, 1 μl of each forward and reverse primer (10 pmol/μl), and 3 μl of genomic DNA template. To amplify each six gene, the PCR reaction conditions were as follows: 7 min at 94℃, 30 cycles of 1 min at the relevant annealing temperature (Table 3), and 1 min at 74℃, followed by 10 min at 74℃.

Table 3. List of gene fragments and primer sequences for C. glabrata MLST analysis

SMGHBM_2020_v30n2_122_t0003.png 이미지

C. glabrata clinical isolates.

The resulting sample was analyzed by gel electrophoresis. The PCR product of all loci were purified and sequenced using reverse sequence primer at Macrogen Inc.. The obtained sequences were analyzed by using the C. glabrata MLST database (http://pubmlst.org/cglabrata/). Each unique sequence at a locus defined an allele number, and unique combinations of alleles assigned as a ST.

Data analysis

The alignment of combined six target gene sequences and loci (3,345 bp) was performed using the Molecular Evolutionary Genetics Analysis (MEGA) v. 7.0 software [12]. For relatedness of the same species, the phylogenetic tree was drawn with the Unweighted Pair Group Method using Arithmetic algorithm (UPGMA) with randomized 1,000 bootstrapping. And then the eBURST package (http://eburst.mlst.net/) was used to determine that all related isolates were grouped into clonal complexes.

Results

PCR and sequence analysis of the fungal ITS region for species identification

1.5% TBE agarose gel DNA electrophoresis data showed that the size of amplified fungal ITS region was 978 bp, and amplicons have shown one clear band (data not shown). As a analysis result of comparison by Genebank BLAST tool for verifying the amplified PCR products, all clinical isolates used in this study were identified as C. glabrata with high concordance rate (97.8%±2.9).

PCR, sequence analysis, obtaining allele number of six housekeeping genes for the MLST analysis

In order to perform the MLST analysis, six housekeeping genes of 102 C. glabrata clinical isolates were amplified by PCR. The size of amplified fragments of FKS, LEU2, NMT1, TRP1, UGP1, and URA3 were 589 bp, 512 bp, 607 bp, 419 bp, 616 bp, and 602 bp, respectively as the expected size and they represented clear band (data not shown).

Amplified six housekeeping gene fragments were sequenced and then sequenced data was trimmed as each length manually by using Chromas software. Trimmed sequence of six housekeeping gene fragments was analyzed at the C. glabrata PubMLST database then allele number of each gene was obtained. Obtained allele number of FKS gene were 3, 5, 7, 8, and 10. Those of LEU2 gene were 5, 6, 16, and 17. Those of NMT1 gene were 3, 4, 6, 8, 11, 14, 19, and 22. Those of TRP1 gene were 2, 3, 5, 7, 10, 12, 19, and 50. Those of UGP1 gene were 1, 3, 5, and 51. Those of URA3 gene were 2, 4, 6, 8, 9, 17, and 20 (Table 4). A total of 49 (1.5%) polymorphic sites were found among six housekeeping gene fragments. The number of polymorphic site in each gene fragment were as follows: FKS (5, 0.9%); LEU2 (4, 1.2%); NMT1 (8, 2.0%); TRP1 (8, 3.1%); UGP1 (4, 0.7%); and URA3 (7, 1.5%), as shown in Table 4. Additionally, insertions, deletions, or heterozygosity were only detected in the NMT1 fragments.

Table 4. The number of polymorphic site and different alleles, and sequence variable position of six housekeeping gene fragments in 102 C. glabrata clinical isolates

SMGHBM_2020_v30n2_122_t0004.png 이미지

Sequence type and cluster of C. glabrata clinical isolates

The MLST scheme revealed a high diversity of C. glabrata isolates with a total of 12 STs, 8 of which were identified as undetermined STs (USTs) that were not discovered in the previous studies.

The data demonstrates that the UST1 was the most predominant ST in this study as a total of 54 clinical isolates (52.9%) were contained in this ST, and the following most predominant ST was the ST63 as a total of 24 clinical isolates (23.5%) were contained in this ST. In addition, this study obtained the ST55, ST22, and ST43 were as a total of 3 (2.9%), 6 (5.9%), and 2 (1.2%) clinical isolates were contained in respective ST and the UST2 was identified in 6 isolates (5.9%). The UST3 was identified in 2 isolates (1.2%). The remaining 5 USTs (UST4, -5, -6, -7, and -8) were classified only once each (1%) (Table 5).

Table 5. Unique sequence types determined with a combination of six loci

SMGHBM_2020_v30n2_122_t0005.png 이미지

Combined sequence (3,345 bp) of six housekeeping genes was used for the phylogenetic tree analysis. Apart from 3 outliers, the isolates were divided into 2 major clusters: cluster 1 and cluster 2 (Fig. 1). Cluster 1 consisted of the UST1 and cluster 2 consisted of the ST63 (Fig. 2, Fig. 3).

SMGHBM_2020_v30n2_122_f0001.png 이미지

Fig. 1. Evolutionary relationships of distinct sequence types. UPGMA dendrogram used to determine the pairwise differences of the concatenated sequence of C. glabrata clinical isolates.

SMGHBM_2020_v30n2_122_f0002.png 이미지

Fig. 2. Phylogenetic tree generated from each isolate. The MLST allele profile, sequence type (ST), and eBURST group (Cluster 1) are listed.

SMGHBM_2020_v30n2_122_f0003.png 이미지

Fig. 3. Phylogenetic tree generated from each isolate. The MLST allele profile, sequence type (ST), and eBURST group (Cluster 2) are listed.

Discussion

C. glabrata is a highly opportunistic pathogen of the urogenital tract and the bloodstream in humans [17]. It is especially prevalent in the elderly and within the human immunodeficiency virus positive population [19]. Although candidiasis is frequently treated with azole antifungal agents, treatment failure has become a serious concern with azole-resistant clinical isolates due to widespread and long-term use of these agents. Nevertheless, few studies have been conducted on the structure, epidemiology, and basic biology of C. glabrata.

Healthcare-associated infections may be endogenous in origin or nosocomially transmitted, and the only way to distinguish them is through strain typing. Recently, MLST directly investigated the DNA sequence variations in a set of housekeeping genes and characterized the strains by their unique allelic profiles. The principle of MLST is simple, involving PCR amplification followed by DNA sequence analysis. Nucleotide differences between strains can be verified at a variable number of genes depending on the desired degree of discrimination. MLST schemes now exist for several important bacterial pathogens including Neisseria meningitidis, Streptococcus pneumoniae, Staphylococcus aureus, Streptococcus pyogenes, and Campylobacter jejuni. The technique has also been used to assess genetic relatedness among strains of Candida spp. including C. glabrata, C. albicans, C. tropicalis, and C. krusei. However, MLST scheme for C. glabrata has only been used by Dodgson, A. R. (2003) and Katiyar, S.(2016).

Hence, in this study, the first MLST analysis with the yeast pathogen C. glabrata was performed and evaluated in Korea. 6 loci were selected for this study, as recommended by previous studies. While a ST3 was defined as prevalent ST in Dodgson, A. R. (2003), the data in this study demonstrates that the UST1 was the most predominant ST. Additionally, the data defined a total of 12 STs among the 102 clinical isolates, and found 8 USTs as a result.

In conclusion, prevalent and novel C. glabrata STs were found in the present study. The data might provide a fundamental database for further studies on C. glabrata, including its epidemiology and evolution. Furthermore, these data might also contribute to the development of novel antifungal agents and diagnostic tests. It might even be possible to discover the virulence factors associated with disease, which population genetic studies currently struggle to monitor.

Acknowledgements

This study was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2017R1C1B5076822) and Brain Busan 21 Plus project.

The Conflict of Interest Statement

The authors declare that they have no conflicts of interest with the contents of this article.

참고문헌

  1. Berila, N. and Subik, J. 2010. Molecular analysis of Candida glabrata clinical isolates. Mycopathologia 170, 99-105. https://doi.org/10.1007/s11046-010-9298-1
  2. Bidet, P., Lalande, V., Salauze, B., Burghoffer, B., Avesani, V., Delmee, M., Rossier, A., Barbut, F. and Petit, J. C. 2000. Comparison of PCR-ribotyping, arbitrarily primed PCR, and pulsed-field gel electrophoresis for typing Clostridium difficile. J. Clin. Microbiol. 38, 2484-2487. https://doi.org/10.1128/JCM.38.7.2484-2487.2000
  3. Cheng, M. F., Yang, Y. L., Yao, T. J., Lin, C. Y., Liu, J. S., Tang, R. B., Yu, K. W., Fan, Y. H., Hsieh, K. S., Ho, M. and Lo, H. J. 2005. Risk factors for fatal candidemia caused by Candida albicans and non-albicans Candida species. BMC Infect. Dis. 5, 22-26. https://doi.org/10.1186/1471-2334-5-22
  4. Da Matta, D. A., Melo, A. S., Guimaraes, T., Frade, J. P., Lott, T. J. and Colombo, A. L. 2010. Multilocus sequence typing of sequential Candida albicans isolates from patients with persistent or recurrent fungemia. Med. Mycol. 48, 757-762. https://doi.org/10.3109/13693780903501689
  5. Da Silva-Rocha, W. P., Lemos, V. L., Svidizisnki, T. I., Milan, E. P. and Chaves, G. M. 2014. Candida species distribution, genotyping and virulence factors of Candida albicans isolated from the oral cavity of kidney transplant recipients of two geographic regions of Brazil. BMC Oral Health 14, 20-28. https://doi.org/10.1186/1472-6831-14-20
  6. Dodgson, A. R., Pujol, C., Denning, D. W., Soll, D. R. and Fox, A. J. 2003. Multilocus sequence typing of Candida glabrata reveals geographically enriched clades. J. Clin. Microbiol. 41, 5709-5717. https://doi.org/10.1128/JCM.41.12.5709-5717.2003
  7. Dodgson, A. R., Pujol, C., Pfaller, M. A., Denning, D. W. and Soll, D. R. 2005. Evidence for recombination in Candida glabrata. Fungal Genet. Biol. 42, 233-243. https://doi.org/10.1016/j.fgb.2004.11.010
  8. Fidel, P. L. Jr., Vazquez, J. A. and Sobel, J. D. 1999. Candida glabrata: review of epidemiology, pathogenesis, and clinical disease with comparison to C. albicans. Clin. Microbiol. Rev. 12, 80-96. https://doi.org/10.1128/CMR.12.1.80
  9. Hesham, A. E. L., Wambui, V., JO, H. O. and Maina, J. M. 2014. Phylogenetic analysis of isolated biofuel yeasts based on 5.8 S-ITS rDNA and D1/D2 26S rDNA sequences. J. Genet. Engineer. Biotechnol. 12, 37-43. https://doi.org/10.1016/j.jgeb.2014.01.001
  10. Katiyar, S., Shiffrin, E., Shelton, C., Healey, K., Vermitsky, J. P. and Edlind, T. 2016. Evaluation of polymorphic locus sequence typing for Candida glabrata epidemiology. J. Clin. Microbiol. 54, 1042-1050. https://doi.org/10.1128/JCM.03106-15
  11. Krcmery, V. and Barnes, A. J. 2002. Non-albicans Candida spp. causing fungaemia: pathogenicity and antifungal resistance. J. Hosp. Infect. 50, 243-260. https://doi.org/10.1053/jhin.2001.1151
  12. Kumar, S., Stecher, G. and Tamura, K. 2016. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870-1874. https://doi.org/10.1093/molbev/msw054
  13. Maiden, M. C., Bygraves, J. A., Feil, E., Morelli, G., Russell, J. E., Urwin, R., Zhang, Q., Zhou, J., Zurth, K., Caugant, D. A., Feavers, I. M., Achtman, M. and Spratt, B. G. 1998. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc. Natl. Acad. Sci. USA. 95, 3140-3145. https://doi.org/10.1073/pnas.95.6.3140
  14. Meats, E., Feil, E. J., Stringer, S., Cody, A. J., Goldstein, R., Kroll, J. S., Popovic, T. and Spratt, B. G. 2003. Characterization of encapsulated and noncapsulated Haemophilus influenzae and determination of phylogenetic relationships by multilocus sequence typing. J. Clin. Microbiol. 41, 1623-1636. https://doi.org/10.1128/JCM.41.4.1623-1636.2003
  15. Odds, F. C. and Jacobsen, M. D. 2008. Multilocus sequence typing of pathogenic Candida species. Eukaryot. Cell 7, 1075-1084. https://doi.org/10.1128/EC.00062-08
  16. Pryce, T. M., Palladino, S., Kay, I. D. and Coombs, G. W. 2003. Rapid identification of fungi by sequencing the ITS1 and ITS2 regions using an automated capillary electrophoresis system. Med. Mycol. 41, 369-381. https://doi.org/10.1080/13693780310001600435
  17. Rahimkhani, M., Saberian, M., Mordadi, A., Varmazyar, S. and Tavakoli, A. 2015. Urinary tract infection with Candida glabrata in a patient with spinal cord injury. Acta. Med. Iran. 53, 516-517.
  18. Rodrigues, C. F., Silva, S. and Henriques, M. 2013. Candida glabrata: a review of its features and resistance. Eur. J. Clin. Microbiol. Infect. Dis. 33, 673-688. https://doi.org/10.1007/s10096-013-2009-3
  19. Sanson, G. F. and Briones, M. R. 2000. Typing of Candida glabrata in clinical isolates by comparative sequence analysis of the cytochrome c oxidase subunit 2 gene distinguishes two clusters of strains associated with geographical sequence polymorphisms. J. Clin. Microbiol. 38, 227-235. https://doi.org/10.1128/JCM.38.1.227-235.2000
  20. Silva, S., Negri, M., Henriques, M., Oliveira, R., Williams, D. W. and Azeredo, J. 2012. Candida glabrata, Candida parapsilosis and Candida tropicalis: biology, epidemiology, pathogenicity and antifungal resistance. FEMS Microbiol. Rev. 36, 288-305. https://doi.org/10.1111/j.1574-6976.2011.00278.x
  21. Spampinato, C. and Leonardi, D. 2013. Molecular fingerprints to identify Candida species. Biomed. Res. Int. 2013, 923742.
  22. Tatay-Dualde, J., Prats-van der Ham, M., de la Fe, C., Gomez-Martin, A., Paterna, A., Corrales, J. C., Contreras, A. and Sanchez, A. 2016. Multilocus sequence typing of Mycoplasma mycoides subsp. capri to assess its genetic variability in a contagious agalactia endemic area. Vet. Microbiol. 191, 60-64. https://doi.org/10.1016/j.vetmic.2016.06.002
  23. Tavanti, A., Davidson, A. D., Johnson, E. M., Maiden, M. C., Shaw, D. J., Gow, N. A. and Odds, F. C. 2005. Multilocus sequence typing for differentiation of strains of Candida tropicalis. J. Clin. Microbiol. 43, 5593-5600. https://doi.org/10.1128/JCM.43.11.5593-5600.2005
  24. Urwin, R. and Maiden, M. C. 2003. Multi-locus sequence typing: a tool for global epidemiology. Trends Microbiol. 11, 479-487. https://doi.org/10.1016/j.tim.2003.08.006
  25. Wang, Y., Shi, C., Liu, J. Y., Li, W. J., Zhao, Y. and Xiang, M. J. 2016. Multilocus sequence typing of Candida tropicalis shows clonal cluster enrichment in azole-resistant isolates from patients in Shanghai, China. Infect. Genet. Evol. 44, 418-424. https://doi.org/10.1016/j.meegid.2016.07.026
  26. Xu, J. 2006. Fundamentals of fungal molecular population genetic analyses. Curr. Issues Mol. Biol. 8, 75-89.