Jpn. J. Infect. Dis., 71, 350–353, 2018

Original Article Reclassification of spp. Isolated in South Korea Using Multilocus Sequence Typing Kyung-Hee Park1, Yeon-Joo Choi1,2, Jeoungyeon Kim1,2, Hye-Jin Park1,2, Dayoung Song1,2, and Won-Jong Jang1,2* 1Department of Microbiology and 2Institute of Glocal Disease Control, College of Medicine, Konkuk University, Seoul, Korea

SUMMARY: Here, we used multilocus sequence typing (MLST) to evaluate 3 intergenic genes (16S rRNA, ospA, and 5S-23S IGS) in Borrelia isolated from South Korea to analyze the relationships be- tween host, vector, and molecular background. We identified B. afzelii, B. yangtzensis, B. garinii, and B. bavariensis. This study is the first report for the identification of B. yangtzensis using MLST in South Korea.

related to B. tanukii and B. valaisiana (23). We postu- INTRODUCTION lated that such reclassifications can be explained by both is a common arthropod-borne epidemic the discovery and addition of new reference species as inflammatory disorder prevalent in Europe and North well as the development of experimental techniques. America (1,2). It is a multi-systemic zoonotic disorder Multilocus sequence typing (MLST) was originally caused by motile, helically-shaped Borrelia spp. proposed as a tool for the epidemiological study of infec- (3,4). Borrelia spp. is comprised of 2 major clades that tious diseases (24,25). MLST has been shown to be a form a monophyletic group in the Order Spirochetes (5). powerful and pragmatic molecular method for typing The Lyme borreliosis B. burgdorferi (sensu lato), trans- large numbers of microbial strains for population genetic mited by various species of Ixodid ticks, forms one studies, delineation of species, and assignment of strains monophyly. Relapsing fever Borrelia, mostly transmitted to defined bacterial species (26–31). Single-locus ap- by Argasid ticks, a few species of Ixodid ticks, and lice, proaches may result in distorted evolutionary patterns forms the other monophyly (6–11). There are several due to different evolutionary processes on different re- species in the genus Borrelia, all of which are parasites, gions of the genome. In contrast, MLST may more accu- either as specialists or generalists, on one or more verte- rately reflect evolutionary relationships (29,32), and has brate species, and are transmitted between these hosts by been advocated as an alternative method to determine a hematophagous arthropod (12). Borrelia spp. (26) because of its simplicity compared to Cases of Lyme disease have been widely reported in traditional DNA-DNA hybridization. A previous study North America, Europe, and East Asian countries includ- used MLST to define B. bavariensis as a new species ing South Korea (13–19). In South Korea, 4 isolates of that was sufficiently genetically distinct from B. garinii. Borrelia spp. were initially identified as both B. afzelii Another study reported that Japanese isolates of B. and B. garinii using polymerase chain reaction (PCR) ­valaisiana-like strains were closely related to Chinese analysis (20). Subsequently, 15 more Borrelia spp. isolates of B. yangtzensis, forming a monophyletic group strains were isolated in Korea, which included HN6, (33). HN7, HN8, HN9, HNM14, and HNM19 (21). All of However, MLST has not yet been used to study these strains were identified as B. garinii, except for ­Korean Borrelia spp. although it has been applied to HN9, which was identified as B. afzelii using PCR Campylobacter spp. and Staphylococcus aureus (34,35). ­analysis. However, when these isolates were classified In this study, we used an MLST scheme based on bac- by rRNA gene restriction fragment length polymor- terial housekeeping genes as described previously (26) phisms (RFLPs), HN6, HN8, HNM14, and HNM19 in addition to previously studied 16S rRNA, 5S-23S ­exhibited some degree of variation from the previously ­intergenic spacer (IGS), and the gene encoding the outer known reference strain, indicating that they may be a surface protein A (ospA) in the first examination of the new strain (22). A UPGMA phylogeny constructed based evolutionary position of Korean Borrelia strain isolates. on 16S rRNA further supported the view that HN6 and The data suggest the reclassification and possible mis- HN8 were clearly separated from B. garinii but closely identification of the previously studied strains.

Received March 27, 2018. Accepted May 11, 2018. MATERIALS AND METHODS J-STAGE Advance Publication May 31, 2018. DOI: 10.7883/yoken.JJID.2018.139 Strains and DNA extraction: Three samples (HN7, *Corresponding author: Mailing address: Department of HN9, and KW3) were isolated from Ixodes granulatus, ­Microbiology, College of Medicine, Konkuk University, and the remaining 4 samples (HN12, HN13, HNM17, 120 Neungdong-ro, Gwangjin-gu, Seoul 143-701, Korea. and KM4) were separated from Apodemus agrarius and Tel: +82-2-2030-7816, Fax: +82-2-2030-7845, E-mail: named according to the name of the locality where they ­[email protected] were found. Five Hae Nam strains (HN7, HN9, HN12,

350 Reclassification of Borrelia in South Korea

HN13, and HNM17), one Kangwon strain (KW3), and at the Borrelia MLST online under ID 2032–2036. for characterization. Genomic DNA was extracted and Sequence alignment and phylogenetic trees: Multi- purified as described previously (20,36). ple sequence alignment was generated with the ClustalW MLST and sequencing: The Borrelia MLST scheme algorithm and MegAlign software by manual inspection. suggested previously (26) is comprised of 8 housekeep- The alignment was made on the translated amino acid ing genes (clpA, clpX, pepX, pyrG, nifS, recG, rplB, and sequences and then back-translated to nucleotide se- uvrA). The PCR conditions for the housekeeping genes, quences to ensure in frame alignment. except for recG, were as follows: 95°C for 15 min, 94°C To construct the phylogenetic tree, the maximum like- for 30 s, annealing temperature from 55°C to 48°C for lihood method based on the Tamura-Nei model of nucle- 60 s, and an extension step of 72°C for 60 s. An addi- otide substitution was determined. An initial tree was tional 30 cycles were performed at 94°C for 30 s, anneal- ­obtained by applying the neighbor-joining and BioNJ ing temperature of 48°C for 60 s, and extension at 72°C ­algorithms to a matrix of pairwise distances estimated for 60 s. After a final extension step for 5 min at 72°C, using the maximum composite likelihood approach, and the samples were kept at 15°C until further analysis. For then selecting the topology with superior log likelihood recG, the PCR conditions for the first set of cycles were value. Evolutionary analyses were conducted in 95°C for 15 min; followed by 30 cycles at 94°C for 30 s, MEGA6, and bootstrap with 1,000 replications was 55°C for 30 s, 72°C for 30 s; and a final extension at ­performed. 72°C for 5 min. The PCR reactions were performed with a VertiTM 96 well Thermal Cycler (Applied Biosystems, RESULTS Foster City, CA, USA). PCR products of all housekeep- ing genes were sequenced in the forward and reverse Identification of Borrelia isolates: The partial se- ­directions (Sequencing was performed by Bioneer Co. quence of the 16S rRNA gene (1,482 bp) of the isolates Ltd., Daejeon, South Korea). The obtained sequences was identical to those of B. burgdorferi (HN13), B. val- were compared and aligned with reference Borrelia spp. aisiana (HN7, HN9, HN12, HNM17, and KM4) and B. gene sequences deposited in the GenBank database to garinii (KW3). Different and partial 5S-23S IGS gene identify published sequences with a high degree of simi- (350 bp) sequences were identical to those of B. valaisi- larity using a BLAST search. MLST sequences were ana (HN12), B. garinii (HN13 and KW3), and B. afzelii compared with the described alleles in the B. burgdorferi (HN7, HN9, HNM17, and KM4). In addition, the ospA MLST database , and new al- gene (402 bp) sequences were identical to those of B. leles were assigned arbitrary number designations. Al- ­valaisiana (HN12), B. garinii (HN13 and KW3), and B. lelic profiles already present in the B. burgdorferi MLST afzelii (HN7, HN9, HNM17, and KM4) (Table 1). database were assigned the corresponding sequence type MLST results: In this study, MLST analysis identi- (ST) number. Novel combinations of allele were as- fied 7 samples, 2 distinct STs, and 5 new ST profiles signed ST numbers, and the sequence data are available (Table 2). Briefly, HN7, HN9, HNM17, and KM4 were

Table 1. Summary of Borrelia genospecies by MLST and 3 genes typing Strain 16S rRNA1) Genospecies2) ospA1) Genospecies2) IGS1) Genospecies2) MLST type HN7 1 B. afzelii 3 B. afzelii 1 B. afzelii B. afzelii HN9 1 B. afzelii 3 B. afzelii 1 B. afzelii B. afzelii HN12 11 B. burgdorferi 23 B. valaisiana 31 B. valaisiana B. yangtzensis HN13 15 B. garinii 41 B. garinii 13 B. garinii B. garinii HNM17 11 B. burgdorferi 3 B. afzelii 1 B. afzelii B. afzelii KW3 5 B. garinii 11 B. garinii 15 B. garinii B. bavariensis KM4 11 B. burgdorferi 3 B. afzelii 1 B. afzelii B. afzelii 1): Alleles of 16S rRNA, ospA, and 5S–23S IGS were assigned numbers based on Borrelia database. (http://www.pubMLST.org/bburgdorferi/). 2): Result of genospecies based on BLAST based on uploaded references.

Table 2. Borrelia MLST sequence types (ST) and 8 allelic profiles Strain Source ST1) clpA clpX nifS pepX pyrG recG rplB uvrA HN7 I. granulatus 7362) 241 202 185 214 220 231 195 215 HN9 I. granulatus 7362) 241 202 185 214 220 231 195 215 HN12 A. agrarius 152 83 61 59 74 69 63 62 65 HN13 A. agrarius 7372) 242 203 29 215 221 232 196 216 HNM17 A. agrarius 7382) 243 204 47 216 222 233 197 215 KW3 I. persulcatus 122 60 45 64 54 50 46 41 45 KM4 A. agrarius 7392) 243 205 47 57 220 234 198 28 1): STs are based on eight housekeeping genes according to the MLST website, www.mlst.net/borrelia.mlst. 2): Novel combination with ST number assigned.

351 Fig. 1. Phylogenetic analyses of concatenated sequences of 8 housekeeping gene fragments. Concatenated se- quences of 8 housekeeping genes were aligned and analyzed in MEGA6. Phylogenetic tree was constructed us- ing the Maximum Likelihood method. The percentage of trees in which the associated taxa clustered together is shown next to branches. Bootstrap test was for 1,000 repetitions (†, Novel ST). identified as B. afzelii, HN12 was identified as B. yangt- as B. garinii. The 16S rRNA gene typing identified 3 zensis, HN13 was identified as B. garinii, and KW3 was isolates as B. burgdorferi, 2 as B. afzelii, and 2 as B. identified as B. bavariensis. Interestingly, HN12 classi- ­garinii. This classification was likely due to mislabeling fied as B. valaisiana in the ospA and 5S-23S IGS analy- of sequences in GenBank. Indeed, the MLST analysis ses were re-classified as B. yangtzensis in the MLST on the 16S rRNA, ospA, and 5S-23S IGS sequences analysis (Tables 1 and 2). ­revealed that these isolates belonged to 4 species: B. Phylogenetic tree of Borrelia strains: The 7 Borrelia ­afzelii, B. garinii, B. yangtzensis and B. bavariensis, as strains were classified into 4 genospecies: B. afzelii shown by the phylogenetic analysis of 8 housekeeping (HN7, HN9, HNM17, and KM4,), B. yangtzensis loci. Using single locus analyses, KW3 was identified (HN12), B. bavariensis (KW3), and B. garinii (HN13) as B. garinii. However, this isolate was identified as B. (Fig. 1). In addition to the experimental samples, other ­bavariensis—a species most closely related to B. STs were inserted into groups (ST330, ST364, ST411, ­garinii—in the MLST analysis. ST18, ST400, ST133, ST585, ST106, and ST424) and Furthermore, a previous study has reported a high ge- the tree was rooted with sequences of the B. parkeri netic heterogeneity of Borrelia within each genospecies, (ST670). including B. burgdorferi, B. garinii, and B. afzelii (37). B. afzelii reportedly has the highest haplotype diversity, which greatly contributes to increased genetic diversity. DISCUSSION Presently, B. afzelii remains a potent genospecies (57.1%, The 7 strains (5 HN, one KW, and one KM) that we n = 4), which displayed genetic heterogeneity with mul- re-analyzed using MLST for the 16S rRNA, ospA, and tiple branches (Fig. 1). Of these 4 strains, HN7 and HN9 5S-23S IGS sequences that have been previously classi- showed the same ST profile (241-202-185-214-220-231- fied through various analyses. For example, one study 195-215), and the rest were analyzed to have different identified HN7 and HN12 as B. lusitaniae and B. valaisi- ST profiles (Table 2). This genetic information may also ana, HN9 and KM4 strains as B. afzelii, and HN7, appear for other types of ST phenotypes according to HNM17, and KM4 as B. afzelii based on a groEL RFLP reservoir. For example, I. granulatus (HN7, HN9) analysis (21). showed the same phenotype as ST706, but A. agrarius To corroborate the affiliation of Borrelia isolates from (HNM17, and KM4) showed different ST phenotypes. South Korea, we performed MLST/16S rRNA, ospA, Originally found in Europe, B. valaisiana is mainly and 5S-23S IGS using conserved housekeeping loci transmitted by I. ricinus and birds (33). However, (­intra-species level) and inter-species levels, respec- slightly different B. valaisiana-related species were iden- tively. Typing using 5S-23S IGS and ospA genes identi- tified as B. yangtzensis in different reservoir hosts (38). fied 4 isolates as B. afzelii, one as B. valaisiana, and one As found presently, HN12 has been identified as B.

352 Reclassification of Borrelia in South Korea yangtzensis. This result is the first report of a Borrelia United States, 1992–2006. MMWR Surveill Summ. 2008;57:1-9. genospecies by MLST analysis in South Korea. Interest- 17. Masuzawa T. Terrestrial distribution of the Lyme borreliosis agent Borrelia burgdorferi sensu lato in East Asia. Jpn J Infect Dis. 2004; ingly, it has the same ST profile (ST152, 83-61-59-74- 57:229-35. 69-63-62-65) separated from reservoir host (A. agrarius) 18. Kim TH, Choi EH, Lee MG, et al. Serologically diagnosed Lyme from China. disease manifesting erythema migrans in Korea. J Korean Med Sci. B. yangtzensis has been reported (named valaisiana- 1999;14:85-8. 19. Ai CX, Wen YX, Zhang YG, et al. Clinical manifestations and epide- related) in China and Japan (17,33,38,39). This genospe- miological characteristics of Lyme disease in Hailin country, Hei- cies was detected only in small rodents rather than birds, longjiang Province, China. Ann NY Acad Sci. 1988;539:302-13. and had a similar ecological phenotype (40). It is pre- 20. Park KH, Chang WH, Schwan TG. Identification and characteriza- sumed that B. yangtzensis classified in this study will be tion of lyme disease Spirochetes, Borrelia burgdoferi sensu lato, iso- more closely related to Chinese strains than Japanese lated in Korea. J Clin Microbiol. 1993;31:1831-7. 21. Park HS, Lee JH, Jeong EJ, et al. Evaluation of groEL gene analysis strains. for identification of Borrelia burgdorferi sensu lato. J Clin Micro- The limitation of this study was that the number of biol. 2004;42:1270-3. identified Borrelia strains was small. Nevertheless, we 22. Park KH, Kim JH, Hahn MJ, et al. Classification of Borrelia burg- successfully surveyed and classify the species of Bor- dorferi sensu lato isolated in Korea by restriction fragment length polymorphism. J Korean Soc Microbiol. 1998;33:119-28. relia that are popular in South Korea. It is important to 23. Lee SH, Kim BJ, Kim JH, et al. Characterization of Borrelia identify and secure strains of predominant species burgdoferi strains isolated from Korea by 16S rDNA sequence analy- through the identification of Borrelia in Korea in the sis and PCR-RFLP analysis of rrf (5S)-rrl (23S) intergenic spacer ­future via the accumulation of analysis data on Borrelia amplicons. Int J Syst Evol Microbiol; 2000;50:857-63. genotypic diversity. 24. Hoen AG, Margos G, Bent SJ, et al. Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events. Proc Natl Acad Sci USA. 2009;106: Acknowledgments Funding for this work was provided by the 15013-8. Korea National Institute of Health, Korea Centers for Disease 25. Maiden MC. Multilocus sequence typing of bacteria. Annu Rev Control and Prevention, Osong, Chungcheongbuk-do (Province), ­Microbiol. 2006;60:561-88. Korea (Title: Novel emerging, variant pathogens and human- 26. Margos G, Vollmer SA, Cornet M, et al. A new Borrelia species derived biological resources bank, Project code No. 2014 ­defined by multilocus sequence analysis of housekeeping genes. ER470202). Appl Environ Microbiol. 2009;75:5410-6. 27. Bishop CJ, Aanensen DM, Jordan GE, et al. Assigning strains to ­bacterial species via the internet. BMC Biol. 2009;7:3. Conflict of interest None to declare. 28. Gevers D, Dawyndt P, Vandamme P, et al. Stepping stones towards a new prokaryotic . Philos. Trans R Soc Lond B Biol Sci. REFERENCES 2006;361:1911-6. 29. Maiden MC, Bygraves JA, Feil E, et al. Multilocus sequence typing: 1. Burgdorfer W, Barbour AG, Hayes SF, et al. Lyme diseases-A tick- A portable approach to the identification of clones within popula- borne Spirochetosis? Science. 1982;216:1317-9. tions of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998; 2. Steere AC. Lyme disease. N Engl J Med.1989;321:586-96. 95:3140-5. 3. Benach JL, Bosler EM, Hanrahan JP, et al. Spirochetes isolated from 30. Richter D, Postic D, Sertour N, et al. Delineation of Borrelia burg- the blood of two patients with Lyme disease. N Engl J Med. 1983; dorferi sensu lato species by multilocus sequence analysis and con- 308:740-2. firmation of the delineation of Borrelia spielmanii sp. nov. Int J Syst 4. Johnson RC, Schmid GP, Hyde FW, et al. Borrelia burgdorferi sp. Evol Microbiol. 2006;56:873-81. nov.: Etiologic agent of Lyme disease. Int J Syst Bacteriol. 1984;34: 31. Spratt BG. Multilocus sequence typing: molecular typing of bacterial 496-7. pathogens in an era of rapid DNA sequencing and the Internet. Curr 5. Barbour AG. Phylogeny of a relapsing fever Borrelia species trans- Opin Microbiol. 1999;2:312-6. mitted by the hard tick Ixodes scapularis. Infect Genet Evol. 2014; 32. Urwin R, Maiden MCJ. Multilocus sequence typing: A tool for 27:551-8. global epidemiology. Trends Microbiol. 2003;11:479-87. 6. Felsenfeld O. Borreliae, human relapsing fever, and parasite-vector- 33. Kawabata H, Takano A, Kadosaka T, et al. Multilocus Sequence host relationships. Bacteriol Rev. 1965;29:46-74. ­Typing and DNA Similarity Analysis Implicates that a Borrelia 7. Trees AJ. The transmission of Borrelia theileri by Boophilus annula- ­valaisiana–related sp. Isolated in Japan is Distinguishable from tus (Say, 1821). Trop Anim Health Prod. 1978;10:93-4. ­European B. valaisiana. J Vet Med Sci. 2013;75:1201-7. 8. Raoult D, Roux V. The body louse as a vector of reemerging human 34. Shin E, Hong H, Park J, et al. Characterization of Staphylococcus diseases. Clin Infect Dis. 1999;29:888-911. aureus faecal isolates associated with food-borne disease in Korea. J 9. Smith RD, Miranpuri GS, Adams JH, et al. Borrelia theileri: isola- Appl Microbiol. 2016;121:277-86. tion from ticks (Boophilus microplus) and tick-borne transmission 35. Shin E, Oh Y, Kim M, et al. Antimicrobial Resistance Patterns and between splenectomized calves. Am J Vet Res. 1985;46:1396-8. Corresponding Multilocus Sequence Types of the Campylobacter 10. Varela-Stokes AS. Transmission of bacterial agents from lone star ­jejuni Isolates from Human Diarrheal Samples. Microb Drug Res. ticks to white tailed deer. J Med Entomol. 2007;44:478-83. 2013;19:110-6. 11. Van Duijvendijk G, Coipan C, Wagemakers A, et al. Larvae of Ixodes 36. Park KH, Hahn MJ, Lee SH, et al. Isolation and antigenic identifica- ricinus transmit Borrelia afzelii and B. miyamotoi to vertebrate hosts. tion of Borrelia burgdoferi from ticks in Korea. J Korean Soc Micro- Parasit Vectors. 2016;9:97. biol. 1994;29:607-17. 12. Rollend L, Fish D, Childs JE. Transovarial transmission of Borrelia 37. Coipan EC, Jahfari S, Fonville M, et al. Imbalanced presence of spirochetes by Ixodes scapularis: a summary of the literature and ­Borrelia burgdorferi s.l. multilocus sequence types in clinical mani- ­recent observations. Ticks Tick Borne Dis. 2013;4:46-51. festations of Lyme borreliosis. Infect Genet Evol. 2016;42:66-76. 13. Walter KS, Carpi G, Caccone A, et al. Genomic insights into the 38. Chu CY, Liu W, Jian BG, et al. Novel genospecies of Borrelia burg- ­ancient spread of lyme disease across North America. Nat Ecol Evol. dorferi sensu lato from rodents and ticks in southwestern China. J 2017;1:1569-76. Clin Microbiol. 2008;46:3130-3. 14. Ni XB, Jia N, Jiang BG, et al. Lyme borreliosis caused by diverse 39. Masuzawa T, Takada N, Kudeken M, et al. Borrelia sinica sp. Nov., genospecies of Borrelia burgdorferi sensu lato in northeastern China. a lyme disease-related Borrelia species isolated in China. Int J Syst Clin Microbiol Infect. 2014;8:808-14. Evol Microbiol. 2001;51:1817-24. 15. Rizzoli A, Hauffe H, Carpi G, et al. Lyme borreliosis in Europe. Euro 40. Margos G, Chu CY, Takano A, et al. Borrelia yangtzensis sp. nov. Surveill. 2011;16:pii:19906. ­rodent-associated species in Asia, is related to Borrelia valaisiana. 16. Bacon RM, Kugeler KJ, Mead PS. Surveillance for Lyme disease - Int J Syst Evol Microbiol. 2015;65:3836-40.

353