Prevalence of tick-borne encephalitis virus in ixodid ticks collected from the republic of Korea during 2011-2012.

OBJECTIVES: In this study, we demonstrated that TBEV-infected ticks have been distributed in the ROK, combined with our previous results. These results suggest that TBEV may exist in the ROK, and H. longicornis, H. flava, and I. nipponensis may be potential vectors of TBEV. In addition, these results emphasize the need for further epidemiological research of TBEV.
METHODS: We examined for the presence of RNA of TBEV by reverse transcriptase-nested polymerase chain reaction (RT-nested PCR) using ixodid ticks captured in 25 localities of 10 provinces. Ticks were collected by the flagging and dragging method or using sentinel BG traps at forests, grass thickets, and grassland. A total of 13,053 ticks belonging to two genera and four species were collected and pooled (1292 pools), according to collection site, species of tick, and developmental stage.
RESULTS: Among 1292 pools, the envelope (E) protein gene of TBEV was detected using RT-nested PCR in 10 pools (3 pools of the 1,331 adult ticks and 7 pools of the 11,169 nymph ticks) collected from Gangwon-do province, Jeonrabuk-do province, and Jeju Island. The minimum infection rates for TBEV of Haemaphysalis longicornis, Haemaphysalis flava, and Ixodes nipponensis were 0.06%, 0.17%, and 2.38%, respectively. Phylogenetic analysis based on the partial E protein gene was performed to identify relationships between the TBEV strains. This showed that 10 Korean strains clustered with the Western subtype.
CONCLUSION: In this study, we investigated the prevalence of tick-borne encephalitis virus (TBEV) in ixodid ticks from various regions of the Republic of Korea (ROK) during 2011-2012 to identify whether TBEV is circulating and to determine the endemic regions of TBEV.

1. Introduction

Ixodid ticks transmit a number of zoonotic pathogens, such as Borrelia burgdorferi [1], Batonella [2], Ehrlichia and Anaplasma [3,4], Rickettsia [5], and tick-borne encephalitis virus (TBEV) [6] to mammalian hosts. As the causative agent of tick-borne encephalitis (TBE), TBEV is a member of the family Flaviviridae and genus Flavivirus, and one of the most important human infections of the central nervous system [7]. TBE occurs in endemic areas of the Eurasian continents, including Europe, Russia, and Far-Eastern Asia (China and Japan), and has a significant impact on public health in these endemic regions [8]. So far, TBEV has been subdivided into three subtypes: the Far-Eastern subtype, known as Russian spring summer encephalitis virus; the Western or European subtype, known as Central European encephalitis virus; and the Siberian subtype based on phylogenetic analysis [8,9]. The main tick vector of the Western subtype is Ixodes ricinus, and Ixodes persulcatus is the main vector for the other subtypes [10,11]. TBEVis transmitted by tick bite and is maintained in the zoonotic transmission cycle between ixodid ticks and wild vertebrate hosts, such as wild and domestic mammals, birds, and reptiles.

In the Republic of Korea (ROK), although TBE infections have not been reported among humans, we have identified molecular evidence of TBEV infections in infesting ticks of the wild animals or collected ticks such as Haemaphysalis longicornis, Haemaphysalis flava, Haemaphysalis japonica, and Ixodes niponensis, which previously had not been known as TBEV vectors, and have isolated TBEV from lung tissues of the wild rodent, Apodemus agrarius [12,13]. However, unlike our expectation, Korean isolates were identified as the Western subtype of TBEV by sequence and phylogenetic analyses compared with other TBEV strains from neighboring countries, including China, Japan, and northeastern Russia that belong to the Far-Eastern subtype [12,13].

In this study, we investigated the prevalence of TBEV in ixodid ticks from various regions of the ROK using sensitive reverse transcriptase-nested polymerase chain reaction (RT-nested PCR) method to identify whether TBEV is circulating and to determine the endemic regions of TBEV.

2. Materials and Methods

2.1. Collection of ixodid ticks

Ixodid ticks surveys were performed by the flagging and dragging method or using sentinel BG traps in various sites, including grass thicket, grassland, and broad-leaved and coniferous forests in 25 localities of 10 provinces of the ROK during 2011–2012. Figure 1 represents the geographical locations of the collection sites. After collection, ticks were placed in plastic tubes and transported to the medical entomology laboratory of Korea National Institute of Health, where they were identified according to their species and developmental stages under a dissecting microscope according to the classification method of Yamaguti et al [14]. Identified ticks were stored at 4 ℃ until further investigation.

Figure 1.

Geographical location of the collection sites in the ROK. Open circles show the sites at which positive tick pools for TBEV were detected [KorVec11-215, KorVec11-224, and KorVec11-226 (N 38°00’01.5’’ E 127°44’29.3’’ , Goseong-ri, Sabuk-myeon, Chuncheon-si), KorVec11-265 (N 38°00’51.1’’ E 126°46’36.2’’, Jugok-ri, Napo-myeon, Gunsan-si), KorVec12- 1131 (N 38°12’9.37’’ E 128°30’55.32’’, Nohak-dong, Sokchosi), KorVec12-1138 (N 38°12’9.11’’ E 128°30’52.84’’, Nohakdong, Sokcho-si), KorVec12-1152 (N 38°12’12.42’’ E 128°31’2.2’’, Nohak-dong, Sokcho-si), KorVec12-1195, KorVec12-1200, and KorVec12-1217 (N 33°25’27.3’’ E 126°33’14.8’’, Ara-dong, Jeju-si)]. CB = Chungcheongbukdo province; CN = Chungcheongnam-do province; GB = Gyeongsangbuk-do province; GG = Gyeonggi-do province; GN = Gyeongsangnam-do province; GW = Gangwon-do province; JB = Jeonllabuk-do province; JN = Jeonllanam-do province.

2.2. Tick processing and RNA extraction

A total of 13,053 identified ticks were pooled according to species, developmental stage, and locality. The size of pools ranged from 1 to 50 in larvae, from 1 to 30 in nymphs, and from 1 to 5 in adult males or females (Table 1). All the pooled ticks were homogenized using Precellys 24 homogenizer (Bertin Technologies, Montigny, Bretonneux, France) with a tissue lysis buffer in RNeasy Mini Kit (Qiagen GmbH, Hilden, Germany) and 2.8 mm stainless-steel beads. The homogenate was centrifuged at 10,000 rpm for 5 minutes, and supernatant was used for RNA extraction. RNA was extracted using the RNeasy Mini Kit (Qiagen GmbH) according to the manufacturer’s instructions.

Table 1.

The number of ixodid ticks classified by species and developmental stage, and MIR of TBEV in this study

Species	Developmental stages	Number of ticks	Number of pools	Positive pools	MIR (%)a
Haemaphysalis longicornis	Larvaeb	350	22	0	0
	Nymphc	10,436	460	6	0.06
	Adult (male)d	221	106	0	0
	Adult (female)d	849	329	1	0.12
	Subtotal	11,856	917	7	0.06
Haemaphysalis flava	Larvae	197	11	0	0
	Nymph	710	163	0	0
	Adult (male)	108	79	1	0.93
	Adult (female)	134	91	1	0.75
	Subtotal	1,149	344	2	0.17
Ixodes nipponensis	Larvae	0	0	0	0
	Nymph	23	9	1	4.35
	Adult (male)	8	8	0	0
	Adult (female)	11	11	0	0
	Subtotal	42	28	1	2.38
Ixodes persulcatus	Larvae	6	3	0	0
	Nymph	0	0	0	0
	Adult (male)	0	0	0	0
	Adult (female)	0	0	0	0
	Subtotal	6	3	0	0
Total	Larvae	553	36	0	0
	Nymph	11,169	632	7	0.06
	Adult (male)	337	193	1	0.30
	Adult (female)	994	431	2	0.20
	Total	13,053	1,292	10	0.08

aMIR: number of positive pools/total number of ticks assayed; . MIR = minimum infection rate; TBEV = tick-borne encephalitis virus.

2.3. Reverse transcriptase-nested polymerase chain reaction

To examine for the presence of TBEV envelope (E) gene, the one-step RT-PCR was carried out using a Maxime RT-PCR PreMix kit (iNtRON Biotechnology, Gyeonggi, Korea) with previously described primers, TBE-913F (5’-TGCACACAYYTGGAAAA CAGGGA-3’) and TBE-1738R (5’-TGGCCACTTTT CAGGTGGTACTTG-3’) [15]. The one-step RT-PCR reaction was performed in a PCR thermal cycler (GeneAmp PCR system 9700, Applied Biosystems, Foster City, CA, USA) under the following conditions: 30 minutes at 45 ℃ for reverse transcription and 5 minutes at 94 ℃ for denaturation as the initial step, followed by 25 cycles of 30 seconds at 94 ℃, 30 seconds at 52 ℃, and 1 minute at 72 ℃, and a final extension step of 5 minutes at 72 ℃. Nested PCR was carried out using an i-StarMaster Mix PCR kit (iNtRON Biotechnology) with previously designed primers, TBE- 1192F (5’-CAGAGTGATCGAGGCTGGGGYAA-3’) and TBE-1669R (5’-AACACTCCAGTCTGGTCTC CRAGGTTGTA-3’) [15]. The nested PCR reaction consisted of an initial denaturation step of 2 minutes at 94 ℃, followed by 30 cycles of 20 seconds at 94 ℃, 10 seconds at 62 ℃, and 20 seconds at 68 ℃, and a final extension step of 5 minutes at 72 ℃. To confirm the PCR products, we analyzed them by agarose gel electrophoresis, stained with SYBR safe DNA gel stain (Invitrogen, Carlsbad, CA, USA).

2.4. Sequencing and phylogenetic analysis

The TBEV-positive products of RT-nested PCR were purified using a QIAquick gel extraction kit (Qiagen GmbH) according to the manufacturer’s instructions and sequenced using ABI Prism BigDye terminator cycle sequencing kits and ABI 3730xl sequencer (Applied Biosystems) at Solgent Inc. (Daejeon, Korea). Results of sequencing were assembled using the SeqMan program implemented in DNASTAR software (version 5.0.6; DNASTAR Inc., Madison, WI, USA) to determine the consensus sequences. The TBEV strains used for phylogenetic analysis are listed in Table 2; multiple sequence alignments were performed using Clustal W implemented in MEGA software version 5 [16]. Phylogenetic analysis was performed using MEGA software version 5 by the maximum likelihood method. The sequences obtained from TBEV-positive products were submitted to GenBank (accession numbers JX282515–JX282524).

Table 2.

TBEV strains used in the phylogenetic analysis

Strain	Subtype	Source of virus	Geographical origin	Year of isolation	GenBank accession no.
Hypr	Western	Human blood	Czech Republic	1953	X75286
Neudoerfl	Western	Ixodes ricinus	Austria	1971	U27495
K23	Western	I. ricinus	Germany	1975	AF091010
Abesettarov	Western	Human blood	Russia	1951	AF091005
Kem I	Western	I. ricinus	Hungary	1952	AF091011
KrM 93	Western	A. agrarius	Republic of Korea	2006	EU276109
KrM 213	Western	A. agrarius	Republic of Korea	2006	EU276110
KrM 215	Western	A. agrarius	Republic of Korea	2006	EU276111
KrM 216	Western	A. agrarius	Republic of Korea	2006	EU276112
KrM 219	Western	A. agrarius	Republic of Korea	2006	EU276113
Oshima 5-10	Far-Eastern	Dog blood	Japan	1995	AB001026
Sofjin-HO	Far-Eastern	Human brain	Russia	1937	AB022703
Senzhang	Far-Eastern	Human brain	China	1953	AY174188
Oshima I-1	Far-Eastern	Ixodes ovatus	Japan	1996	AB022292
Oshima 5-11	Far-Eastern	Dog blood	Japan	1995	AB022290
KH98-2	Far-Eastern	Ixodes persulcatus	Russia	1998	AB022295
KH98-5	Far-Eastern	I. persulcatus	Russia	1998	AB022296
KH98-10	Far-Eastern	I. persulcatus	Russia	1998	AB022297
Crimea	Far-Eastern	I. persulcatus	Ukraine	1987	AF091008
D1283	Far-Eastern	Human brain	Russia	1998	AB049347
RK1424	Far-Eastern	I. persulcatus	Latvia	1977	AF091015
T-blood	Far-Eastern	Human blood	Russia	1939	AF091019
VL99-m11	Far-Eastern	I. persulcatus	Russia	1999	AB049345
N132	Far-Eastern	I. persulcatus	Russia	1979	AF091013
Vasilchenko	Siberian	Human blood	Russia	1969	M97659
Aina	Siberian	Human blood	Russia	1963	AF091006
IR99-1m1	Siberian	I. persulcatus	Russia	1999	AB049348
IR99-1m4	Siberian	I. persulcatus	Russia	1999	AB049349
IR99-2m3	Siberian	I. persulcatus	Russia	1999	AB049350
IR99-2m7	Siberian	I. persulcatus	Russia	1999	AB049351
IR99-2f7	Siberian	I. persulcatus	Russia	1999	AB049352
IR99-2f13	Siberian	I. persulcatus	Russia	1999	AB049353

TBEV = tick-borne encephalitis virus.

3. Results

3.1. Numbers and identification of ixodid ticks

A total of 13,053 ticks (553 larvae, 11,169 nymphs, 337 males, and 994 females) were collected from 2011 to 2012 in 25 localities of the ROK (Table 1 and Figure 1). The ixodid tick samples were identified to belong to four species in two genera: H. longicornis, H. flava, I. niponensis, and I. persulcatus. Of the identified ticks, H. longicornis (90.8%, 11,856/ 13,053) was the most abundant species in this study, followed by H. flava (8.8%, 1149/13,053), I. niponensis (0.3%, 42/13,053), and I. persulcatus (0.05%, 6/13,053).

3.2. Prevalence of TBEV in ixodid ticks

Of the 13,053 ixodid tick samples, 10 pools were positive for TBEV according to RT-nested PCR method. Among the 10 positive pools, three pools of the adult ticks and seven of the nymph ticks were collected from Gunsan in Jeollabuk-do province, Chuncheon or Sokcho in Gangwon-do province, and Jeju in Jeju Island. Regional prevalence of TBEV for each species was shown in Table 3. The minimum infection rate (MIR, calculated with the assumption that a positive pool contains one infected tick) of TBEV in H. longicornis, H. flava, and Ixodes nipponensis was 0.06%, 0.17%, and 2.38%, respectively (Table 1). Although H. longicornis was the most common species collected in this study, its positive rates of TBEV was lower than that of other species.

Table 3.

Regional prevalence of TBEV in ixodid ticks (minimum infection rate) collected in the Republic of Korea

Collection site	Larvae	Nymph	Male	Female	Total	Larvae	Nymph	Male	Female	Total
	Haemaphysalis flava: no. infected/no. examined (%)					Haemaphysalis flava: no. infected/no. examined (%)
Gyeonggi-do
Yangpyeong	–	0/12(0)	–	0/3(0)	0/15(0)	–	0/238(0)	0/4(0)	0/9(0)	0/251(0)
Gimpo	–	–	0/2(0)	0/1(0)	0/3(0)	–	0/4(0)	–	0/1(0)	0/5(0)
Incheon
Ganghwa	–	0/7(0)	–	0/1(0)	0/8(0)	–	0/523(0)	0/47(0)	0/43(0)	0/613(0)
Gangwon-do
Chuncheon	–	0/64(0)	1/6(16.67)	0/5(0)	1/75(1.33)	–	2/354(0.56)	–	0/4(0)	2/358(0.56)
Pyeongchang	0/17(0)	0/6(0)	0/11(0)	0/11(0)	0/45(0)	0/147(0)	0/1(0)	–	–	0/148(0)
Sokcho	–	0/125(0)	0/3(0)	1/6(16.67)	1/134(0.75)	–	1/122(0.82)	0/1(0)	0/5(0)	1/128(0.78)
Chungcheongbuk-do
Boeun	0/17(0)	0/217(0)	0/52(0)	0/59(0)	0/506(0)	0/190(0)	0/3598(0)	0/25(0)	0/353(0)	0/4166(0)
Chungcheongnam-do
Boryeong	–	0/18(0)	0/3(0)	0/5(0)	0/26(0)	–	0/436(0)	0/14(0)	0/27(0)	0/477(0)
Cheonan	–	0/1(0)	–	0/1(0)	0/2(0)	–	0/232(0)	–	0/41(0)	0/273(0)
Gongju	–	0/23(0)	0/2(0)	0/4(0)	0/29(0)	0/1(0)	0/359(0)	0/12(0)	0/19(0)	0/391(0)
Jeollabuk-do
Muju	–	0/14(0)	0/1(0)	0/3(0)	0/18(0)	–	0/29(0)	–	–	0/29(0)
Gunsan	–	0/42(0)	0/10(0)	0/7(0)	0/59(0)	–	1/516(0.19)	–	0/18(0)	1/534(0.19)
Jeonju	–	0/3(0)	–	–	0/3(0)	–	0/9(0)	–	0/2(0)	0/11(0)
Jeollanam-do
Boseong	–	0/11(0)	0/1(0)	0/2(0)	0/14(0)	–	0/269(0)	0/3(0)	0/18(0)	0/290(0)
Gokseong	–	0/5(0)	–	0/2(0)	0/7(0)	–	0/109(0)	0/1(0)	0/32(0)	0/142(0)
Wando	–	0/3(0)	0/2(0)	0/1(0)	0/6(0)	–	0/38(0)	0/28(0)	0/47(0)	0/113(0)
Gyeongsangbuk-do
Uiseong	–	0/21(0)	0/2(0)	0/2(0)	0/25(0)	–	0/76(0)	0/1(0)	0/2(0)	0/79(0)
Yeongju	–	0/57(0)	0/1(0)	0/5(0)	0/63(0)	–	0/83(0)	–	–	0/83(0)
Yeongyang	–	0/12(0)	–	–	0/12(0)	–	0/422(0)	0/14(0)	0/81(0)	0/517(0)
Uljin	–	0/15(0)	0/1(0)	0/3(0)	0/19(0)	–	0/574(0)	0/22(0)	0/69(0)	0/665(0)
Gyeongsangnam-do
Haman	–	0/1(0)	–	–	0/1(0)	–	0/52(0)	0/9(0)	0/33(0)	0/94(0)
Hamyang	–	0/3(0)	0/3(0)	0/2(0)	0/8(0)	–	0/1(0)	0/1(0)	–	0/2(0)
Miryang	–	0/6(0)	0/6(0)	0/7(0)	0/19(0)	–	0/80(0)	–	0/1(0)	0/81(0)
Jeju-do
Jeju	–	0/35(0)	–	–	0/35(0)	0/6(0)	2/1,418(0.14)	0/25(0)	1/26(3.85)	3/1,475(0.20)
Seogwipo	0/2(0)	0/9(0)	0/2(0)	0/4(0)	0/17(0)	0/6(0)	0/893(0)	0/14(0)	0/18(0)	0/931(0)
Total	0/197(0)	0/710(0)	0/710(0) 1/108(0.93)	1/134(0.75)	2/1,149(0.17)	0/350(0)	6/10,436(0.06)	0/221(0)	1/849(0.12)	7/11,856(0.06)
Gyeonggi-do
Yangpyeong	–	–	–	–	–	–	–	–	–	–
Gimpo	–	–	–	–	–	–	–	–	–	–
Incheon
Ganghwa	–	–	–	–	–	–	–	–	–	–
Gangwon-do
Chuncheon	–	0/1(0)	–	0/1(0)	0/2(0)	–	–	–	–	–
Pyeongchang	–	–	–	–	–	0/6(0)	–	–	–	0/6(0)
Sokcho	–	1/18(5.56)	0/1(0)	–	1/19(5.26)	–	–	–	–	–
Chungcheongbuk-do
Boeun		0/2(0)	0/4(0)	0/8(0)	0/14(0)	–	–	–	–	–
Chungcheongnam-do
Boryeong	–	–	0/1(0)	–	0/1(0)	–	–	–	–	–
Cheonan	–	–	–	–	–	–	–	–	–	–
Gongju	–	0/1(0)	-	–	0/1(0)	–	–	–	–	–
Jeollabuk-do
Muju		–	–	0/1(0)	0/1(0)	–	–	–	–	–
Gunsan		0/1(0)	–	–	0/1(0)	–	–	–	–	–
Jeonju		–	–	–	–	–	–	–	–	–
Jeollanam-do
Boseong	–	–	–	–	–	–	–	–	–	–
Gokseong	–	–	–	–	–	–	–	–	–	–
Wando	–	–	–	–	–	–	–	–	–	–
Gyeongsangbuk-do
Uiseong	–	–	0/1(0)	-	0/1(0)	–	–	–	–	–
Yeongju	–	–	0/1(0)	0/1(0)	0/2(0)	–	–	–	–	–
Yeongyang	–	–	–	–	–	–	–	–	–	–
Uljin	–	–	–	–	–	–	–	–	–	–
Gyeongsangnam-do
Haman	–	–	–	–	–	–	–	–	–	–
Hamyang	–	–	–	–	–	–	–	–	–	–
Miryang	–	–	–	–	–	–	–	–	–	–
Jeju-do
Jeju	–	–	–	–	–	–	–	–	–	–
Seogwipo	–	–	–	–	–	–	–	–	–	–
Total	-	1/23(4.35)	0/8(0)	0/11(0)	1/42(2.38)	0/6(0)	–	–	–	0/6(0)

ROK = Republic of Korea; TBEV = tick-borne encephalitis virus.

The highest MIR of TBEV infection was identified in I. nipponensis collected from Sokcho in Gangwon-do province (5.26%). In the other localities, the MIR of infected ticks varied from 0% to 1.33%. The highest rate of TBEV infection was identified in I. nipponensis nymph (4.35%). The infection rates in males (0.30%) were higher than those in females (0.20%) and nymph (0.06%). The mean TBEV prevalence in 25 collection sites was about 0.08%.

3.3. Sequence and phylogenetic analyses

Nucleotide sequence identities between the 10 Korean strains and the other 32 TBEV strains represented that the 10 Korean strains had high identity with the Western subtype strains with 97.2–99.6%, compared with the Far-Eastern subtype with 81.1–84.3% or the Siberian subtype with 83.9–85.7%.

Phylogenetic trees derived from nucleotide sequences of the E gene showed that the 10 Korean strains (Kor- Vec11-215, KorVec11-224, KorVec11-226, KorVec11-265, KorVec12-1131, KorVec12-1138, KorVec12-1152, KorVec12-1195, KorVec12-1200, and KorVec12-1217) identified in this study belong to the Western subtype. An envelope gene-based phylogenetic tree was divided into three distinct groups such as Far-Eastern, Siberian, and Western subtypes and it was also shown that the 10 Korean strains belonged to the same subgroup as the Western subtype, as shown in Figure 2.

Figure 2.

Phylogenetic analysis of TBEV strains based on the envelope region. Phylogenetic trees constructed using the ML method based on the Tamura-Nei model in MEGA version 5 (1000 bootstrap replicates). The ML tree were rooted with the sequences of the OHFV (Genbank accession no. X66694). The scale bar indicates the nucleotide substitutions per position. The 10 Korean strains identified in this study are marked with open circles. ML = maximum likelihood; OHFV = Omsk hemorrhagic fever virus.

4. Discussion

So far, TBEV vectors have been reported in 10 species of three genera (I. ricinus, I. persulcatus, Ixodes hexagonus, Ixodes arboricola, Ixodes ovatus, Haemaphysalis punctata, Haemaphysalis concinna, Haemaphysalis inermis, Dermacentor marginatus, and Dermacentor reticulatus) [17,18]. Among the known vectors, the total distribution of five species of TBEV vector (I. persulcatus, I. ovatus, Haemaphysalis concinna, D. marginatus, and D. reticulatus) in the ROK have been reported [19]. In our previous study, we suggested that TBE might exist in the ROK based on the following evidence. First, there was molecular evidence of TBEV infection in mammalian hosts and in potential vector ticks (H. longicornis, H. flava, H. japonica, and I. nipponensis) collected from the ROK [12,13,20], although these ticks had not been previously reported as TBEV vectors. Second, we reported on the isolation and identification of five TBEV strains from lung tissue homogenates of wild rodents captured in the ROK using in vitro and in vivo experiments [12]. Contrary to our expectations, the sequence comparisons and phylogenetic analyses based on the complete E gene or fullgenome of TBEV Korean isolates compared with other TBEV strains indicated that all the five Korean strains belonged to the Western subtype [21,22].

Based on this evidence, we carried out more extensive survey throughout the ROK for the prevalence of TBEV in ixodid ticks from various localities to identify whether TBEV is circulating and to determine the endemic areas of TBEV. A total of 13,053 ticks were divided into 1292 pools to identify for TBEV by RTnested PCR method, which detected 10 TBEV-positive pools. Sequence and phylogenetic analyses based on the envelope gene sequence revealed that the 10 Korean strains from infected ticks in this study belonged to the Western subtype, in accordance with our previous results.

The principal tick vector of the Western subtype of TBEV is I. ricinus [8], but in this study, it is reasonable to suppose that H. longicornis, H. flava, and I. nipponensis serve as the potential vectors of TBEV in the ROK. Until recently, the molecular evidence of TBEV infection was identified in ixodid ticks collected from Yangpyeong or Dongducheon in Gyeonggi-do province, Pyeongchang or Jeongseon in Gangwon-do province, and Jeju Island [13,20]. In this study, ixodid ticks such as H. longicornis, H. flava, and I. nipponensis collected from Gunsan in Jeollabuk-do province, Chuncheon or Sokcho in Gangwon-do province, and Jeju in Jeju Island were identified according to the molecular evidence of TBEV infection. These results indicated that TBEV may be endemic in these localities of theROKand H. longicornis, H. flava, and I. nipponensis may be potential vectors of the TBEV Western subtype, combined with previous findings. However, in order to prove these points, isolation of TBEV from infected vectors and characterization of TBEV isolates will be required.

In summary, we found that TBEV-infected ticks have been distributed in some localities of the ROK. These results emphasize the need for further epidemiological research of TBEV and preventive measures against the occurrence of TBE in the ROK.