Skip to main content
Advertisement
  • Loading metrics

Molecular genomic characterization of tick- and human-derived severe fever with thrombocytopenia syndrome virus isolates from South Korea

  • Seok-Min Yun ,

    Contributed equally to this work with: Seok-Min Yun, Su-Jin Park

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing

    Current address: Division of Emerging Infectious Disease and Vector Research, Center for Infectious Diseases Research, National Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

    Affiliations Division of Arboviruses, National Research Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea, College of Medicine and Medical Research Institute, Chungbuk National University, Cheongju-si, Republic of Korea

  • Su-Jin Park ,

    Contributed equally to this work with: Seok-Min Yun, Su-Jin Park

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Validation, Writing – original draft, Writing – review & editing

    Affiliation College of Medicine and Medical Research Institute, Chungbuk National University, Cheongju-si, Republic of Korea

  • Sun-Whan Park,

    Roles Methodology, Resources, Writing – review & editing

    Current address: Division of Emerging Infectious Disease and Vector Research, Center for Infectious Diseases Research, National Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

    Affiliation Division of Arboviruses, National Research Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

  • WooYoung Choi,

    Roles Resources, Writing – review & editing

    Current address: Division of Viral Diseases, Center for Laboratory Control of Infectious Diseases, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

    Affiliation Division of Arboviruses, National Research Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

  • Hye Won Jeong,

    Roles Resources, Writing – review & editing

    Affiliation College of Medicine and Medical Research Institute, Chungbuk National University, Cheongju-si, Republic of Korea

  • Young-Ki Choi ,

    Roles Conceptualization, Data curation, Funding acquisition, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

    choiki55@chungbuk.ac.kr (YKC); leewonja@gmail.com (WJL)

    Affiliation College of Medicine and Medical Research Institute, Chungbuk National University, Cheongju-si, Republic of Korea

  • Won-Ja Lee

    Roles Conceptualization, Data curation, Funding acquisition, Project administration, Supervision, Validation, Writing – review & editing

    choiki55@chungbuk.ac.kr (YKC); leewonja@gmail.com (WJL)

    Current address: Division of Emerging Infectious Disease and Vector Research, Center for Infectious Diseases Research, National Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

    Affiliation Division of Arboviruses, National Research Institute of Health, Korea Centers for Disease Control and Prevention, Cheongju-si, Republic of Korea

Abstract

Background

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging tick-borne viral disease caused by the SFTS virus (SFTSV) from Bunyaviridae that is endemic in East Asia. However, the genetic and evolutionary characteristics shared between tick- and human-derived Korean SFTSV strains are still limited.

Methodology/Principal findings

In this study we identify, for the first time, the genome sequence of a tick (Haemaphysalis longicornis)-derived Korean SFTSV strain (designated as KAGWT) and compare this virus with recent human SFTSV isolates to identify the genetic variations and relationships among SFTSV strains. The genome of the KAGWT strain is consistent with the described genome of other members of the genus Phlebovirus with 6,368 nucleotides (nt), 3,378 nt, and 1,746 nt in the Large (L), Medium (M) and Small (S) segments, respectively. Compared with other completely sequenced human-derived Korean SFTSV strains, the KAGWT strain had highest sequence identities at the nucleotide and deduced amino acid level in each segment with the KAGWH3 strain which was isolated from SFTS patient within the same region, although there is one unique amino acid substitution in the Gn protein (A66S). Phylogenetic analyses of complete genome sequences revealed that at least four different genotypes of SFTSV are co-circulating in South Korea, and that the tick- and human-derived Korean SFTSV strains (genotype B) are closely related to one another. Although we could not detect reassortant, which are commonly observed in segmented viruses, further large-scale surveillance and detailed genomic analysis studies are needed to better understand the molecular epidemiology, genetic diversity, and evolution of SFTSV.

Conclusions/Significance

Full-length sequence analysis revealed a clear association between the genetic origins of tick- and human-derived SFTSV strains. While the most prevalent Korean SFTSV is genotype B, at least four different genotypes of SFTSV strains are co-circulating in South Korea. These findings provide information regarding the molecular epidemiology, genetic diversity, and evolution of SFTSV in East Asia.

Author summary

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging tick-borne viral disease caused by the SFTS virus (SFTSV). During entomological surveillance of SFTSV infection in Korean ticks collected from SFTS outbreak areas, we isolated a single SFTSV strain which we designated KAGWT. In addition, we isolated three SFTSVs from human patients with typical SFTS symptoms. In this study, we report the genomic sequences of each of these isolates and compare the genetic and evolutionary characteristics between tick- and human-derived Korean SFTSV isolates. Genetic and phylogenetic analyses of these sequences revealed that the tick-derived Korean SFTSV strain is clustered into genotype B, the most prevalent genotype in South Korea, and was closely related to other SFTSV in the same group. Furthermore, our results show that at least four different genotypes of SFTSV strains are co-circulating in South Korea.

Introduction

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging tick-borne viral disease characterized by fever, gastrointestinal symptoms, leukopenia, and thrombocytopenia. It was first reported in China in 2010 [1] and was subsequently identified in South Korea and Japan in 2013 [24]. The causative agent, SFTS virus (SFTSV), belongs to the genus Phlebovirus in the family Bunyaviridae [1]. Other novel tick-borne phleboviruses, including Heartland virus (HRTV), Malsoor virus (MV) and Hunter Island Group virus (HIGV), which are genetically related to but distinctly different from SFTSV, have been isolated from humans and ticks in the United States [5, 6], bats in India [7] and ticks in Australia [8], respectively.

Like other members of the genus Phlebovirus, SFTSV contains a tripartite RNA genome consisting of three single-stranded RNA segments of negative polarity, designated large (L), medium (M), and small (S). The L, M, and S segments encode the RNA-dependent RNA polymerase (RdRp), the viral envelope glycoproteins (Gn and Gc) and both a nucleoprotein (NP) and a nonstructural protein (NSs) in an ambisense orientation, respectively [1, 9].

Although human-to-human transmission of SFTSV through contact with blood and/or body secretions of patients has been reported [1012], the virus is generally transmitted to humans by tick bites. Several studies have reported SFTSV isolation or detection from tick species including Haemaphysalis longicornis, Rhipicephalus microplus, H. flava, H. concinna, Amblyomma testudinarium, and Ixodes nipponensis [1, 1318]. Further, this virus has also been isolated or detected from domestic animals (e.g., cattle, goats, dogs, chickens and cats), small mammals such as rodents and shrews [1922], and reptiles [23].

Since the first reported fatal case in 2012 in Gangwon Province in South Korea [2], concern regarding SFTSV has grown as the numbers of SFTS patients has increased annually with 36 cases reported in 2013, 55 cases in 2014, 79 cases in 2015, and 165 cases in 2016 [24]. Moreover, the mean mortality rate of these cases was approximately 21.8%. During our previous survey of carrier ticks from affected areas, a single SFTSV strain (designated KAGWT) was initially isolated from H. longicornis nymphs collected from the Samcheok-si, Gangwon Province [15]. In addition, we also isolated SFTSV from two recovered- and one fatal-human cases that presented with typical SFTS symptoms. In this study, we analyzed the whole genome sequence of the first tick-derived Korean SFTSV and human SFTSV strains to compare the genetic and evolutionary characteristics between tick- and human-derived Korean SFTSV isolates. Genetic characterization revealed that the tick-derived SFTSV is closely associated with recent KAGWH3 Korean human isolates and that at least four different genotypes of SFTSV strain are co-circulating in South Korea. Taken together, our results suggest that more intensive and continuous surveillance of SFTSV is essential for better understanding of the molecular epidemiology, genetic diversity, and evolution of this virus in East Asia.

Methods

Ethics statement

Chungbuk National University Hospital received written consent for sample collection from each patient with SFTSV infection. All participants were adults and this study was approved by the institutional review board (IRB) of Chungbuk National University Hospital (IRB no. 2015-08-009-001).

Virus propagation, cell lines, RNA extraction and cDNA synthesis

The KAGWT strain was isolated from H. longicornis ticks collected from the Samcheok-si, Gangwon Province in South Korea as described previously [15]. The CB1, CB2, and CB3 strains were isolated from the sera of SFTS patients who were hospitalized with typical SFTS symptoms at Chungbuk National University Hospital.

For virus propagation, the virus was passaged five times on confluent monolayers of Vero E6 cells (ATCC No. CRL-1586; American Type Culture Collection, Manassas, VA) in Dulbecco’s Modified Eagle Medium (DMEM; Gibco, Grand Island, NY) containing 8% fetal bovine serum (FBS; Gibco) with penicillin (100IU/mL) and streptomycin (100μg/mL; P/S, Gibco) placed in 37°C incubator supplemented with 5% CO2. Cell culture supernatant was collected after seven days and stored at -80°C as the working virus stock for whole genome sequencing.

Viral RNA was extracted from 140 μL of the viral stock using QIAamp Viral RNA Mini Kits (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. Single strand cDNA was synthesized using viral RNA and primers specific for each segment using a cDNA synthesis kit (Cosmogenetech co, Ltd., Seoul, Korea) according to the manufacturer’s instructions.

Whole genome sequencing by Sanger method

For the whole genome sequencing, one to six overlapping PCR fragments covering the SFTSV genome containing full-length L, M, and S segments were amplified by PCR from cDNA using SP-Taq DNA polymerase (Cosmogenetech, Seoul, Korea). PCR for each segment was performed at 95°C for 5 minutes, followed by 45 cycles of amplification consisting of 95°C for 30 seconds, annealing at 53 or 60°C for 30 seconds according to the primer sets for each segment, and 72°C for 2 minutes, with a final extension at 72°C for 5 minutes. The non-coding 5′ and 3′ ends of the viral genome were determined by rapid amplification of cDNA ends (RACE) method as described previously [25]. PCR products were then purified using a QIAquick Gel Extraction Kit (Qiagen) according to the manufacturer’s instructions and direct sequenced using ABI Prism BigDye Terminator Cycle Sequencing Kits (Applied Biosystems, Foster City, CA) and an ABI 3730x1 sequencer (Applied Biosystems) at Cosmogenetech Co, Ltd. The nucleotide sequences obtained from this study were assembled using the SeqMan program in the DNASTAR software (version 5.0.6; DNASTAR Inc., Madison, WI) to determine the complete genomic sequence.

Genetic and phylogenetic analyses

Genetic and phylogenetic analyses were conducted by aligning published full-length sequences of SFTSV obtained from China, Japan and Korea isolates, which are available in GenBank, together with the closely related SFTSV sequences obtained from the basic local alignment search tool results. A total of 41 full-length sequences of the L, M, and S segments, including the strains isolated in this study (Table 1), were included in the analyses.

thumbnail
Table 1. Sequence information for SFTSV strains identified in this study.

https://doi.org/10.1371/journal.pntd.0005893.t001

Multiple sequence alignments were performed using the Clustal W algorithm in DNASTAR version 5.0.6 or MEGA version 6.0 [27]. The aligned nucleotide and deduced amino acid sequences were analyzed using the MegAlign program of DNASTAR to compare the sequence homologies and amino acid substitutions.

Phylogenetic analyses were performed based on the full-length L, M, and S segments of SFTSVs. Phylogenetic trees were constructed with MEGA version 6.0 software using the Maximum Likelihood (ML) method based on the Kimura 2-parameter model. The reliability of the ML tree was evaluated by the bootstrap test with 1,000 replications.

Accession numbers

The full-length L, M, and S segment sequences of tick- and human-derived strains determined in this study have been deposited in GenBank under the following accession numbers: KAGWT (KY273136 to KY273138), CB1 (KY789433, KY789436, and KY789439), CB2 (KY789434, KY789437, and KY789440), and CB3 (KY789435, KY789438, and KY789441) strains.

Results

Genome organization of tick- and human-derived Korean SFTSV strains

The genomes of tick- and human-derived Korean SFTSV strains analyzed in this study were consistent with 6,368 nucleotide (nt), 3,378 nt and 1,746 nt present in the L, M and S segments respectively, consistent with what has been reported for other members of the genus Phlebovirus [9]. The L segment encodes a 6,255 nt long ORF [2,084 amino acids (aa)] for the RNA-dependent RNA polymerase gene, the M segment comprises a 3,222 nt precursor of the glycoprotein gene, coding for Gn (516 aa) and Gc (511 aa) proteins, while the S segment contains 882 nt and 738 nt long ORFs, which translate into a nonstructural protein (293 aa) and a nucleoprotein (245 aa), respectively. The non-coding regions (NCRs) of the L, M and S segments at the 5′ termini are 16 nt, 18 nt and 43 nt, respectively; and at the 3′ termini are 97 nt, 138 nt and 29 nt, respectively. As shown in Fig 1, complementary sequences within the 5′ and 3′ NCRs of the three segments are highly conserved between tick- and human-derived Korean SFTSV strains.

thumbnail
Fig 1. Nucleotide sequence alignment of complementary sequences within the 5′ and 3′-noncoding regions (NCRs) of L (blue boxes), M (yellow boxes), and S (green boxes) segments of 18 completely sequenced SFTSV strains, including tick- and human-derived Korean strains (KAGWT, CB1, CB2, and CB3).

Compared with the sequence of the KAGWT strain, the different, and deleted sequences are indicated by open boxes, and hyphens, respectively. The closed boxes represent conserved complementary sequences within the 5′ and 3′ NCRs of the three segments.

https://doi.org/10.1371/journal.pntd.0005893.g001

Genetic analyses based on the complete L, M, S segment sequences of SFTSV strains

Alignment and pairwise comparisons of each segment between the tick-derived KAGWT strain and other fully sequenced human-derived Korean SFTSV strains showed nucleotide (and deduced amino acid) sequence homology ranging from 95.9 (99.3) to 99.9 (100), 94.0 (98.4) to 99.8 (99.8), and 94.9 (98.3) to 99.9 (100)% for L, M and S segments, respectively (Table 2). This high homology between the tick-derived KAGWT and other human-derived SFTSV strains reflects the close level of relatedness of these strains to each other.

thumbnail
Table 2. Comparison of the full-length nucleotide sequence (ORF region) and deduced amino acid sequence identities between the tick-derived KAGWT strain and other human-derived Korean SFTSV strains.

https://doi.org/10.1371/journal.pntd.0005893.t002

In particular, the KAGWT strain showed the highest sequence identities at the nucleotide and deduced amino acid levels with the KAGWH3 strain isolated in 2014 from the serum of a patient from Gangwon Province, South Korea who experienced typical SFTS symptoms [25]. These results indicate that tick- and human-derived Korean SFTSV strains are most closely related with one another. Moreover, the deduced amino acid sequence of KAGWT revealed two amino acid variations (A66S, I89V) in the Gn protein compared with the KAGWH3 strain. In particular, one unique amino acid substitution in the Gn protein, at position 66 (alanine to serine), was found only in the KAGWT strain compared with other human SFTSVs. According to previous reports, a change from phenylalanine to serine at position 330 in the M segment polyprotein (F330S) occurred in cell culture-adapted SFTSV and resulted in the large-focus phenotype [28]. However, all SFTSV strains analyzed in this study had a conserved amino acid sequence (F) at position 330 of the M segment.

Compared with the other fully sequenced SFTSV strains, KAGWT showed higher identity with the genotype B SFTSV strains than with the other genotypes at each segment as shown in S1 Table. The nucleotide sequence identities between the tick-derived KAGWT and other SFTSV strains from China, Japan, and South Korea belonging to the genotype B were 96.4 to 99.9% (L), 95.8 to 99.8% (M), 94.9 to 99.9% (NP) and 95.8 to 99.8% (NSs) similar, respectively. However, KAGWT showed relatively low nucleotide sequence identity ranging from 95.7 to 96.3% (L genes), 93.6 to 94.4% (M genes), 95.1 to 96.2% (NP genes) and 94.8 to 95.9% (NSs genes), respectively with other genotypes of SFTSV (Table 2 and S1 Table).

In comparison with the other fully sequenced SFTSV strains, CB2 was identified as the first genotype A SFTSV strain out of 14 full-length sequenced Korean SFTSV isolates. Genetic comparison results showed that the CB2 strain had nucleotide sequence identity ranging from 97.9% (L), 97.9 to 98.0% (M), 98% (NP) and 97.2% (NSs), respectively with the other genotype A SFTSV strains circulating in China. However, the CB2 strain showed relatively low nucleotide sequence identity ranging from 96.1 to 96.7% (L genes), 93.5 to 96.2% (M genes), 95.3 to 97.6% (NP genes) and 93.9 to 96.8% (NSs genes), respectively with other SFTSV strains (Table 2 and S1 Table).

Phylogenetic analyses and genotype distribution

To investigate the genetic evolutionary origins and relationship between the tick- and human-derived Korean SFTSV strains, phylogenetic analyses were conducted with full-length complete sequences of previous SFTSV isolates from China, Japan, and South Korea. Phylogenetic analyses of complete genome sequences (L, M, and S segments) indicated that the tick-derived Korean SFTSV strain (KAGWT) was clustered with the genotype B SFTSV strains circulating in humans in China, Japan, and South Korea (Figs 24), although KAGWT also exhibited a close relationship with the KAGWH3 human isolate [25].

thumbnail
Fig 2. Phylogenetic analysis based on the complete nucleotide sequences of L segment of SFTSV strains using the Maximum Likelihood (ML) method based on the Kimura 2-parameter model.

The numbers on the branches indicate bootstrap percentages based on 1,000 replications, and the scale bar indicates the nucleotide substitutions per site. The phylogenetic branches were supported with greater than 70% bootstrap values. The red, blue, brown, green, purple, and grey branches were designated as SFTSV strains belonging to the genotypes A, B, C, D, E, and F, respectively. The tick- and human-derived Korean SFTSV strains analyzed in this study are marked with red and blue closed circles, respectively.

https://doi.org/10.1371/journal.pntd.0005893.g002

thumbnail
Fig 3. Phylogenetic analysis based on the complete nucleotide sequences of M segment of SFTSV strains using the Maximum Likelihood (ML) method based on the Kimura 2-parameter model.

The numbers on the branches indicate bootstrap percentages based on 1,000 replications, and the scale bar indicates the nucleotide substitutions per site. The phylogenetic branches were supported with greater than 70% bootstrap values. The red, blue, brown, green, purple, and grey branches were designated as SFTSV strains belonging to the genotypes A, B, C, D, E, and F, respectively. The tick- and human-derived Korean SFTSV strains analyzed in this study are marked with red and blue closed circles, respectively.

https://doi.org/10.1371/journal.pntd.0005893.g003

thumbnail
Fig 4. Phylogenetic analysis based on the complete nucleotide sequences of S segment of SFTSV strains using the Maximum Likelihood (ML) method based on the Kimura 2-parameter model.

The numbers on the branches indicate bootstrap percentages based on 1,000 replications, and the scale bar indicates the nucleotide substitutions per site. The phylogenetic branches were supported with greater than 70% bootstrap values. The red, blue, brown, green, purple, and grey branches were designated as SFTSV strains belonging to the genotypes A, B, C, D, E, and F, respectively. The tick- and human-derived Korean SFTSV strains analyzed in this study are marked with red and blue closed circles, respectively.

https://doi.org/10.1371/journal.pntd.0005893.g004

Furthermore, the CB1 and CB3 human isolates also clustered with genotype B viruses although they are separated into a different sub-node from tick-driven Korean SFTSV strains (KAGWT) and are closely related with recent Korean SFTSV human isolates, KADGH and KACNH3. In addition, the CB2 human isolate was clustered together with recent China genotype A viruses.

Overall the phylogenetic trees revealed that Korean SFTSV strains can be classified with high branch support into four distinct genotypes (designated A, B, D, and F) out of the six genotypes described previously [26]. To investigate the prevalence of each genotype of SFTSV, we analyzed the SFTSV sequences available in the GenBank database. As shown in Table 3, most SFTSV strains from South Korea belong to genotype B (11 out of 14 isolates) and only one isolate each was reported from genotype A, D, and F. It is noteworthy that the CB2 is the first report case of a genotype A virus in South Korea. Isolates from China were diverse with all six genotypes being represented, although the most prevalent was genotype F followed by genotypes A and D (Table 3). In contrast, only three SFTSV strain genotypes (B, C, and D) observed in Japan and the most prevalent genotype was B.

thumbnail
Table 3. Genotype distribution of SFTSV strains isolated from China, Japan, and South Korea.

https://doi.org/10.1371/journal.pntd.0005893.t003

Discussion

SFTS is an emerging tick-borne infectious disease caused by a novel Phlebovirus, SFTSV that is highly endemic in China, Japan, and South Korea [1, 2, 4].

In this study, we determined the whole genome sequences of the first tick and patient-derived Korean SFTSV strains and compared them with other available whole genomic SFTSV sequences. To the best of our knowledge, this is the first report of the whole genomic sequence of an SFTSV strain isolated from ticks and of a genotype A SFTSV strain collected from South Korea. Comparison of the whole genome sequence of tick- and human-derived strains revealed that all SFTSVs consist of three segments with 6,368 nt in the L segment, 3,378 nt in the M segment, and 1,746 nt in the S segment. Thus, the genome organization of Korean strains are consistent with the known genome organization of SFTSV belonging to the Phlebovirus genus [9]. Complementary sequences within the 5′ and 3′ NCRs of the three segments analyzed in this study are highly conserved, as in other SFTSV strains (Fig 1). According to previous studies of Bunyaviruses, these complementary sequences form panhandle-like structures and may have important roles in transcription initiation, viral RNA replication, viral RNA encapsidation, and viral genome packaging [29, 30]. All Korean SFTSV strains have 53 unique amino acid variations in the L, Gn, Gc, N, and NSs proteins compared with other SFTSV strains. Among them, one unique amino acid substitution in the Gn protein at position 66 (alanine to serine) was found only in the tick-driven KAGWT strain. Since little information is known about SFTSV, additional research regarding the association of amino acid variations in each gene product with the pathogenesis of SFTSV infection and detailed molecular studies utilizing reverse genetics will be required to further explain the functional role of these unique substitutions and in particular, the Gn protein substitution (66S).

Genetic and phylogenetic analysis revealed that the tick-driven KAGWT strain isolated from H. longicornis nymphs in 2013 [15] has high sequence identity in each segment and is closely related with the KAGWH3 strain isolated in 2014 from the serum of an SFTS patient who experienced high fever, vomiting, diarrhea, and fatigue. This is not surprising given that both viruses were from the Samcheok-si, Gangwon Province [25], and further suggests that tick- and human-derived Korean SFTSV strains are closely related to each other.

Although several studies about SFTSV classification have been reported [26, 3134], the uniform classification of SFTSVs has not yet to be established. Therefore, to use unified nomenclature of SFTSV genotypes, more mutual understanding and discussion are needed. In this study, we adapted the nomenclature as previously described by Fu et al. [26] which described results of most large numbers (205 SFTSV strains) and complete sequences. That paper showed SFTSV strains can be divided into six distinct genotypes. Through full-length sequence analysis of human-driven SFTSV isolates we detected the first genotype A strain (CB2) from a patient exhibiting severe SFTSV-like symptoms in the Chungbuk Province of South Korea in 2015. The Korean strains belong to four of these (genotypes A, B, D, and F) of which genotype B strains appear to be predominant. It should be noted that in China all six genotypes of SFTSV have been reported, while only one isolate was reported for each genotype A, D and F in Korea. Thus, further detailed surveillance of tick- and human-derived SFTSVs are needed to understand the actual prevalence and possible transmission of each genotype of SFTSV strain into South Korea. Due to the segmented nature of the Phlebovirus, genetic reassortment is known to be an important process resulting in the sequence diversity necessary for viral evolution [35]. Previous studies have shown reassortment occurs in Phleboviruses including Hantavirus [36, 37], Rift Valley fever virus (RVFV) [38], and SFTSV [26, 3134]. Although SPL087A and AHL/China/2011 SFTSV strains isolated in Japan and China, respectively, were identified as reassortants which contained different genotype segments in the same SFTSV strain [26, 3134], no reassortant was identified in the Korean SFTSV strains tested in this study. Therefore, continuous analysis consisting of in-depth surveillance and genome sequencing will be needed in South Korea to obtain more detailed information of the molecular epidemiology, genetic diversity, and evolution of SFTSV.

In conclusion, we report here the first whole genomic sequence of the tick-derived SFTSV KAGWT strain isolated from the ticks in the Gangwon province of South Korea and the comparison of the genetic characteristics of this virus with recent human SFTSV isolates. Genetic and phylogenetic analyses with full-length genome sequences revealed that tick- and human-derived Korean SFTSV strains are closely related to one another and that at least four different genotypes of SFTSV are co-circulating in South Korea. These results provide insight into the genetic origins of human of SFTSV strains as well as shed light on the molecular epidemiology, genetic diversity, and evolution of SFTSV. Furthermore, these associations will have important implications for the design of diagnostic procedures and vaccines against SFTSV.

Supporting information

S1 Table. Nucleotide sequence (ORF region) and deduced amino acid sequence identities of (A) L, (B) M, and (C and D) S segments of SFTSV strains belonging to the different genotypes.

Nucleotide identities (%) are shown above the diagonal and the deduced amino acid identities (%) are shown below the diagonal. 1, CB2 strain; 2, HL/Injected strain; 3, HL/Egg/G2 strain; 4, HL/Larvae/G2 strain; 5, HL/Nymph/G2 strain; 6, HL/Adult/G2 strain; 7, KAGWT strain; 8, KAGWH3 strain; 9, KAGBH5 strain; 10, KAGBH6 strain; 11, KACNH3 strain; 12, KADGH strain; 13, CB1 strain; 14, CB3 strain; 15, KAGNH strain; 16, KAGNH4 strain; 17, KAJNH2 strain; 18, YG1 strain; 19, SPL003A strain; 20, SPL010A strain; 21, SPL035A strain; 22, Zhao strain; 23, ZJZHSH-FDE/China/06/2012 strain; 24, ZJZHSH-LWL/China/08/2014 strain; 25, ZJZHSH-WRF/China/08/2014 strain; 26, SPL087A strain; 27, AHL/China/2011 strain; 28, KASJH strain; 29, HB154/China/2011 strain; 30, HB155/China/2011 strain; 31, HB156/China/2011 strain; 32, SD24 strain; 33, HB29 strain; 34, JS4 strain; 35, SD4 strain; 36, Gangwon/Korea/2012 strain; 37, SDLZtick12/2010 strain; 38; JS2012-tick01 strain; 39, AH12 strain; 40, QD7 strain; 41, AH-YTY/China/05/2012 strain.

https://doi.org/10.1371/journal.pntd.0005893.s001

(DOCX)

Acknowledgments

We thank our collaborators for tick surveillance: Ye-Ji Lee at the Division of Arboviruses, National Research Institute of Health, Korea Centers for Disease Control and Prevention; Bong Gu Song, Won Il Park, Jong Yul Roh, and Young-Ran Ju at the Division of Medical Entomology, National Research Institute of Health, Korea Centers for Disease Control and Prevention.

References

  1. 1. Yu XJ, Liang MF, Zhang SY, Liu Y, Li JD, Sun YL, et al. Fever with thrombocytopenia associated with a novel bunyavirus in China. N Engl J Med. 2011;364(16):1523–32. Epub 2011/03/18. pmid:21410387.
  2. 2. Kim KH, Yi J, Kim G, Choi SJ, Jun KI, Kim NH, et al. Severe fever with thrombocytopenia syndrome, South Korea, 2012. Emerg Infect Dis. 2013;19(11):1892–4. Epub 2013/11/12. pmid:24206586.
  3. 3. Park SW, Han MG, Yun SM, Park C, Lee WJ, Ryou J. Severe fever with thrombocytopenia syndrome virus, South Korea, 2013. Emerg Infect Dis. 2014;20(11):1880–2. Epub 2014/10/24. pmid:25341085.
  4. 4. Takahashi T, Maeda K, Suzuki T, Ishido A, Shigeoka T, Tominaga T, et al. The first identification and retrospective study of Severe Fever with Thrombocytopenia Syndrome in Japan. J Infect Dis. 2014;209(6):816–27. Epub 2013/11/16. pmid:24231186.
  5. 5. McMullan LK, Folk SM, Kelly AJ, MacNeil A, Goldsmith CS, Metcalfe MG, et al. A new phlebovirus associated with severe febrile illness in Missouri. N Engl J Med. 2012;367(9):834–41. Epub 2012/08/31. pmid:22931317.
  6. 6. Savage HM, Godsey MS Jr., Lambert A, Panella NA, Burkhalter KL, Harmon JR, et al. First detection of heartland virus (Bunyaviridae: Phlebovirus) from field collected arthropods. Am J Trop Med Hyg. 2013;89(3):445–52. Epub 2013/07/24. pmid:23878186.
  7. 7. Mourya DT, Yadav PD, Basu A, Shete A, Patil DY, Zawar D, et al. Malsoor virus, a novel bat phlebovirus, is closely related to severe fever with thrombocytopenia syndrome virus and heartland virus. J Virol. 2014;88(6):3605–9. Epub 2014/01/07. pmid:24390329.
  8. 8. Wang J, Selleck P, Yu M, Ha W, Rootes C, Gales R, et al. Novel phlebovirus with zoonotic potential isolated from ticks, Australia. Emerg Infect Dis. 2014;20(6):1040–3. Epub 2014/05/27. pmid:24856477.
  9. 9. Liu Q, He B, Huang SY, Wei F, Zhu XQ. Severe fever with thrombocytopenia syndrome, an emerging tick-borne zoonosis. Lancet Infect Dis. 2014;14(8):763–72. Epub 2014/05/20. pmid:24837566.
  10. 10. Gai Z, Liang M, Zhang Y, Zhang S, Jin C, Wang SW, et al. Person-to-person transmission of severe fever with thrombocytopenia syndrome bunyavirus through blood contact. Clin Infect Dis. 2012;54(2):249–52. Epub 2011/11/19. pmid:22095565.
  11. 11. Liu Y, Li Q, Hu W, Wu J, Wang Y, Mei L, et al. Person-to-person transmission of severe fever with thrombocytopenia syndrome virus. Vector Borne Zoonotic Dis. 2012;12(2):156–60. Epub 2011/10/01. pmid:21955213.
  12. 12. Kim WY, Choi W, Park SW, Wang EB, Lee WJ, Jee Y, et al. Nosocomial transmission of severe fever with thrombocytopenia syndrome in Korea. Clin Infect Dis. 2015;60(11):1681–3. Epub 2015/02/20. pmid:25694652.
  13. 13. Zhang YZ, Zhou DJ, Qin XC, Tian JH, Xiong Y, Wang JB, et al. The ecology, genetic diversity, and phylogeny of Huaiyangshan virus in China. J Virol. 2012;86(5):2864–8. Epub 2011/12/23. pmid:22190717.
  14. 14. Meng K, Sun W, Cheng Z, Guo H, Liu J, Chai T. First detection of severe fever with thrombocytopenia syndrome virus in the tick species Haemaphysalis concinna in Shandong Province, China. Parasitol Res. 2015;114(12):4703–7. Epub 2015/09/10. pmid:26350381.
  15. 15. Yun SM, Song BG, Choi W, Roh JY, Lee YJ, Park WI, et al. First Isolation of Severe Fever with Thrombocytopenia Syndrome Virus from Haemaphysalis longicornis Ticks Collected in Severe Fever with Thrombocytopenia Syndrome Outbreak Areas in the Republic of Korea. Vector Borne Zoonotic Dis. 2016;16(1):66–70. Epub 2016/01/09. pmid:26745758.
  16. 16. Yun SM, Lee WG, Ryou J, Yang SC, Park SW, Roh JY, et al. Severe fever with thrombocytopenia syndrome virus in ticks collected from humans, South Korea, 2013. Emerg Infect Dis. 2014;20(8):1358–61. Epub 2014/07/26. pmid:25061851.
  17. 17. Park SW, Song BG, Shin EH, Yun SM, Han MG, Park MY, et al. Prevalence of severe fever with thrombocytopenia syndrome virus in Haemaphysalis longicornis ticks in South Korea. Ticks Tick Borne Dis. 2014;5(6):975–7. Epub 2014/08/29. pmid:25164614.
  18. 18. Yun SM, Lee YJ, Choi W, Kim HC, Chong ST, Chang KS, et al. Molecular detection of severe fever with thrombocytopenia syndrome and tick-borne encephalitis viruses in ixodid ticks collected from vegetation, Republic of Korea, 2014. Ticks Tick Borne Dis. 2016;7(5):970–8. Epub 2016/05/24. pmid:27211914.
  19. 19. Liu JW, Wen HL, Fang LZ, Zhang ZT, He ST, Xue ZF, et al. Prevalence of SFTSV among Asian house shrews and rodents, China, January-August 2013. Emerg Infect Dis. 2014;20(12):2126–8. Epub 2014/11/25. pmid:25418111.
  20. 20. Zhao L, Zhai S, Wen H, Cui F, Chi Y, Wang L, et al. Severe fever with thrombocytopenia syndrome virus, Shandong Province, China. Emerg Infect Dis. 2012;18(6):963–5. Epub 2012/05/23. pmid:22608264.
  21. 21. Niu G, Li J, Liang M, Jiang X, Jiang M, Yin H, et al. Severe fever with thrombocytopenia syndrome virus among domesticated animals, China. Emerg Infect Dis. 2013;19(5):756–63. Epub 2013/05/08. pmid:23648209.
  22. 22. Cui F, Cao HX, Wang L, Zhang SF, Ding SJ, Yu XJ, et al. Clinical and epidemiological study on severe fever with thrombocytopenia syndrome in Yiyuan County, Shandong Province, China. Am J Trop Med Hyg. 2013;88(3):510–2. Epub 2013/01/23. pmid:23339197.
  23. 23. Suh JH, Kim HC, Yun SM, Lim JW, Kim JH, Chong ST, et al. Detection of SFTS Virus in Ixodes nipponensis and Amblyomma testudinarium (Ixodida: Ixodidae) Collected From Reptiles in the Republic of Korea. J Med Entomol. 2016. Epub 2016/03/10. pmid:26957392.
  24. 24. Korea CDC. Disease web statistics system. 2017. http://is.cdc.go.kr/dstat. Accessed 19 April 2017. Korean.
  25. 25. Yun MR, Park SW, Kwon T, Lee S, Yoo WG, Choi W, et al. Full-genome sequences of severe Fever with thrombocytopenia syndrome virus, isolated from South Korea in 2014. Genome Announc. 2015;3(2). Epub 2015/04/18. pmid:25883276.
  26. 26. Fu Y, Li S, Zhang Z, Man S, Li X, Zhang W, et al. Phylogeographic analysis of severe fever with thrombocytopenia syndrome virus from Zhoushan Islands, China: implication for transmission across the ocean. Sci Rep. 2016;6:19563. Epub 2016/01/26. pmid:26806841.
  27. 27. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. Epub 2013/10/18. pmid:24132122.
  28. 28. Brennan B, Li P, Zhang S, Li A, Liang M, Li D, et al. Reverse genetics system for severe fever with thrombocytopenia syndrome virus. J Virol. 2015;89(6):3026–37. Epub 2015/01/02. pmid:25552716.
  29. 29. Guu TS, Zheng W, Tao YJ. Bunyavirus: structure and replication. Adv Exp Med Biol. 2012;726:245–66. Epub 2012/02/03. pmid:22297517.
  30. 30. Lei XY, Liu MM, Yu XJ. Severe fever with thrombocytopenia syndrome and its pathogen SFTSV. Microbes Infect. 2015;17(2):149–54. Epub 2014/12/17. pmid:25498868.
  31. 31. Ding NZ, Luo ZF, Niu DD, Ji W, Kang XH, Cai SS, et al. Identification of two severe fever with thrombocytopenia syndrome virus strains originating from reassortment. Virus Res. 2013;178(2):543–6. Epub 2013/09/24. pmid:24055465.
  32. 32. Liu JW, Zhao L, Luo LM, Liu MM, Sun Y, Su X, et al. Molecular Evolution and Spatial Transmission of Severe Fever with Thrombocytopenia Syndrome Virus Based on Complete Genome Sequences. PLoS One. 2016;11(3):e0151677. Epub 2016/03/22. pmid:26999664.
  33. 33. Shi J, Hu S, Liu X, Yang J, Liu D, Wu L, et al. Migration, recombination, and reassortment are involved in the evolution of severe fever with thrombocytopenia syndrome bunyavirus. Infect Genet Evol. 2017;47:109–17. Epub 2016/11/26. pmid:27884653.
  34. 34. Lv Q, Zhang H, Tian L, Zhang R, Zhang Z, Li J, et al. Novel sub-lineages, recombinants and reassortants of severe fever with thrombocytopenia syndrome virus. Ticks Tick Borne Dis. 2017;8(3):385–90. Epub 2017/01/25. pmid:28117273.
  35. 35. Perez-Losada M, Arenas M, Galan JC, Palero F, Gonzalez-Candelas F. Recombination in viruses: mechanisms, methods of study, and evolutionary consequences. Infect Genet Evol. 2015;30:296–307. Epub 2014/12/30. pmid:25541518.
  36. 36. Li D, Schmaljohn AL, Anderson K, Schmaljohn CS. Complete nucleotide sequences of the M and S segments of two hantavirus isolates from California: evidence for reassortment in nature among viruses related to hantavirus pulmonary syndrome. Virology. 1995;206(2):973–83. Epub 1995/02/01. pmid:7856108.
  37. 37. Rodriguez LL, Owens JH, Peters CJ, Nichol ST. Genetic reassortment among viruses causing hantavirus pulmonary syndrome. Virology. 1998;242(1):99–106. Epub 1998/03/17. pmid:9501041.
  38. 38. Sall AA, Zanotto PM, Sene OK, Zeller HG, Digoutte JP, Thiongane Y, et al. Genetic reassortment of Rift Valley fever virus in nature. J Virol. 1999;73(10):8196–200. Epub 1999/09/11. pmid:10482570.