Abstract
Tick-borne encephalitis virus (TBEV) has been classified into three subtypes, namely the European (Eu-TBEV), Far Eastern (FE-TBEV), and Siberian (Sib-TBEV). In this study, we discovered a new subtype of TBEV in wild rodent Marmota himalayana in Qinghai-Tibet Plateau in China, proposed as subtype Himalayan (Him-TBEV). Two complete genomes of TBEV were obtained from respiratory samples of 200 marmots. The phylogenetic analysis using the E protein and polyprotein demonstrated that the two strains of Him-TBEV formed an independent branch, separated from Eu-TBEV, Sib-TBEV, and FE-TBEV. The nomenclature of Him-TBEV as a new subtype was also supported by comparative analysis using nucleotide and amino acid sequences of E protein and polyprotein. For E protein, The Him-TBEV showed 82.6–84.6% nucleotide identities and 92.7–95.0% amino acid identities with other three subtypes. For polyprotein, the Him-TBEV showed 83.5–85.2% nucleotide identities and 92.6–94.2% amino acids identities with other three subtypes. Furthermore, of 69 amino acid substitutions profiles detected in complete polyprotein of 112 strains of TBEV, Him-TBEV subtype displayed unique amino acids in the 36 positions. Notably, for the subtype-specific amino acid position 206 of E protein, Him-TBEV shared the Val with Eu-TBEV, but differed from FE-TBEV and Sib-TBEV. The evolutionary analysis with BEAST suggested that Him-TBEV diverged from other subtypes of eastern TBEV group about 2469 years ago. It should be mentioned that Qinghai-Tibet Plateau in China is the plague endemic region where Marmota himalayana is the primary host. The public health significance of discovery of Him-TBEV in Marmota himalayana must be carefully evaluated.
Introduction
Tick-borne encephalitis (TBE) is a severe central nervous system infection caused by tick-borne encephalitis virus (TBEV). Over 10,000 human infections have been reported annually in endemic regions of Eurasian continent, especially in Russia.Citation1–Citation3 The causative agent TBEV is a member of genus Flavivirus, family Flaviviridae. The virus genome is a positive single-stranded RNA molecule, approximately 11,000 bases in length, which encodes three structural (capsid, C; membrane, M, which is expressed as precursor prM; envelope, E) and seven non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5).Citation4 Based on phylogenetic analysis of E protein, the TBEV has been classified into three subtypes, namely European (Eu-TBEV), Far-Eastern (FE-TBEV), and Siberian (Sib-TBEV).Citation5, Citation6 The Sib-TBEV and FE-TBEV were also called eastern TBEV. Recently, 886-84-like strains were divided from Sib-TBEV and proposed as a new subtype named Baikalian (Bkl-TBEV).Citation7 Members of these three subtypes showed difference in geographical distribution, virulence, clinical severity, and the difference in the nucleotide (14.6–16.5%) and amino acid (5.0–6.8%) sequences of their polyprotein.Citation8 In China, the first human infection was traced back to 1943, and the strains of TBEV were first isolated from patients and ticks in 1952.Citation9, Citation10 Recently, the number of human infections in China has increased notably. A total of 2117 human infections were reported in the period from 2006 to 2013. Two subtypes of TBEV had been identified in China, namely FE-TBEV and Sib-TBEV. The FE-TBEV was widely detected in several provinces, such as Heilongjiang, Jilin, Liaoning, Xinjiang, Yunnan, and Tibet. The Sib-TBEV was only detected in Xinjiang.Citation11–Citation14 Furthermore, TBEV has been detected in previously unaffected areas,Citation3, Citation12, Citation15 making the disease a growing public health threat. In this study, we discovered two new strains of TBEV in wild rodent Marmota himalayana in Qinghai-Tibet Plateau in China, which represented a new subtype that had been concealed for hundred years.
Results
The respiratory virome of Marmota himalayana
Only 66 contigs with a length range of 125–7419 nt were assembled from RNA-seq sequencing showing similarity to RNA viral sequences. These contigs could be classified into six RNA virus families and one RNA virus genus: Flaviviridae (n = 27), Phenuiviridae (n = 8), Picornaviridae (n = 3), Picobirnaviridae (n = 20), Dicistroviridae (n = 1), Narnaviridae (n = 3) and Sobemovirus (n = 4). Importantly, 3 contigs with lengths of 2965 nt, 4647 nt and 7419 nt showed 99% amino acid identity to Marmota Himalayana hepatovirus (MHHAV) within family Picornaviridae which has been reported by Yu.Citation16 Contigs with a length range of 125–775 nt showed 89.1–100% amino acid identities to different strains of tick-borne encephalitis virus (TBEV) within family Flaviviridae.
Discovery of genomes of TBEV from Marmota himalayana
The existence of TBEV in Marmota himalayana was revealed by RNA-seq sequencing. Two TBEV-positive low respiratory tissues were detected by RT-PCR method. The complete genome sequence data for two strains of TBEV from Marmota himalayana in the present study were obtained by combination method of RNA-seq sequencing and RT-PCR products amplified by a set of primers covering the whole genome region. The two genomes shared the same length of 10,761 nucleotides, with 99.5% identity. The complete genome sequences were submitted to the GenBank with accession numbers MG599476 and MG599477. The single open reading frame (ORF) of the genome with 10,245-nt in length that encodes three structural proteins (C, prM, and E) and 7 nonstructural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) and was flanked by two non-coding regions (NCR), with a 132-nt NCR at 5′ end and a 384-nt NCR at 3′ end. The genome organization and sequence analysis identified the two genomes as members of TBEV, which were named as Himalaya-1 and Himalaya-2, respectively.
Genome analysis
The nucleotide sequences of the E protein and polyprotein were used to reveal the phylogenetic relationship of the Himalaya-1 and Himalaya-2 with the members of other TBEV lineages. The complete nucleotide sequences of E protein (1488 nt) and polyprotein (10,245 nt) of 23 virus genomes were included in the study, including 7 strains of FE-TBEV, 4 strains of Sib-TBEV, 3 strains of Eu-TBEV, 1 strain of Baikalian subtype and TBEV strain 178–79. In addition, 4 strains from animals in European countries were also included, such as Turkish sheep encephalitis virus (TSEV), Greek goat encephalitis virus (GGEV), Louping ill virus (LIV) and Spanish sheep encephalitis virus (SSEV), which closely related to Eu-TBEV.
Phylogenetic trees were constructed by the neighbor-joining (NJ) method with 1000 bootstrap replications. For phylogenetic tree of E protein (Fig. ), Omsk hemorrhagic fever virus (OHFV) was used as outgroup, the TBEV strains were formed two clades. One clade included Eu-TBEV, TSEV, GGEV, LIV, and SSEV, which discovered from western European countries. Another clade included eastern TBEV group (FE-TBEV, Sib-TBEV and Bkl-TBEV), Himalaya-1 and Himalaya-2. Clearly, the two TBEV strains from Marmota himalayana formed a separate lineage, from where the FE-TBEV and Bkl-TBEV were further diverged.
a The phylogenetic tree of E protein. b The phylogenetic tree of polyprotein. Strains of Him-TBEV were indicated in boldface. Phylogenetic trees were constructed using MEGA 6.0 with neighbor-joining method (1000 bootstrap replications). Bootstrap values (>70%) are shown at the branches. Scale bar below indicates the nucleotide substitutions per site
When the complete nucleotide sequences of polyprotein (10,245 nt) from the same set of strains were used to construct the phylogenetic tree with same method, the same topology was obtained (Fig. ). These results strongly suggested that the Himalaya-1 and Himalaya-2 are indeed representing a new subtype of TBEV, which named Himalayan subtype (Him-TBEV), differed from FE-TBEV, Sib-TBEV, and Eu-TBEV as well.
We then conducted comparative analysis of E protein and polyprotein of between and within the subtypes. The nucleotide and amino acid identities of E protein of Himalaya-1 and Himalaya-2 were both of 99.6%. The nucleotide and amino acid identities of polyprotein between Himalaya-1 and Himalaya-2 were 99.5 and 99.8%, respectively (Supplementary Table S1). For polyprotein, the Him-TBEV showed 83.5–85.2% nucleotide identities and 92.6–94.2% amino acid identities with other three subtypes; the identity level of nucleotide and amino acid among Eu-TBEV, Sib-TBEV, and FE-TBEV subtypes were 83.2%–85.5% and 93.0%–95.2%, respectively. For E protein, the Him-TBEV shared 82.6–84.6% nucleotide identities and 92.7–95.0% amino acid identities with other three subtypes; the identity level of nucleotide and amino acid among the Eu-TBEV, Sib-TBEV, and FE-TBEV subtypes were 82.9%–85.7% and 95.2%–98.0%, respectively (Supplementary Table S1).
The unique signature amino acid sequence profile of Him-TBEV
The published genomes from a total of 112 strains of TBEV currently available in GenBank were downloaded for analysis, including 60 strains of FE-TBEV, 22 strains of Sib-TBEV and 30 strains of Eu-TBEV (Supplementary Table S2).
Notably, the analysis with deduced amino acid sequences of Himalaya-1 and Himalaya-2 revealed 69 amino acid substitutions in the complete polyprotein (Table ). Most of these substitutions were located at E, NS1, NS2A, and NS4B proteins. No substitutions were found in the well-known features of E protein, such as fusion peptide, tick-specific peptide EHLPTA, 12 cysteine residues.Citation17 Three substitutions (positions 16, 143, and 151) were located in domain I, two substitutions (positions 130 and 201) were located in domain II and three substitutions (positions 331, 346 and 349) were located in domain III of E protein.Citation18 The substitution of amino acid position 201 of E protein (E → A) would increase the net positive charge of the virus surface, as with mutant (E → K) described previouslyCitation19.
Amino acid substitutions between Himalayan subtype and other three subtypesa
Of 69 amino acid substitutions observed in the complete polyprotein, Him-TBEV has 36 unique substitutions. At these 36 positions, the amino acids were conserved in all other three subtypes, showing identical amino acids. However, Him-TBEV displayed unique amino acids in the 36 positions (Table ). These data strongly supported the hypothesis that the Him-TBEV is indeed a new subtype of TBEV. Interestingly, for the so-called subtype-specific amino acid position 206 of E protein, Him-TBEV shared the Val with that of Eu-TBEV, but differed from that of FE-TBEV and Sib-TBEV.Citation5
Of the 17 pathogenicity associated amino acid residues, the Him-TBEV shared 9 substitutions that are specific to pathogenic strains, and 5 substitutions that are specific to strains from patients with the subclinical presentation. These pathogenicity associated amino acid residues were proposed by analyzing the polyprotein of pathogenic strains isolated from the brains of dead patients with the encephalitis, from the blood of patients with febrile, and from patients with subclinical symptoms.Citation20, Citation21 In short, the profile of pathogenic associated amino acid substitution of Him-TBEV is similar to low virulence strain Oshima 5–10 (Table ).
17 key amino acid positions associated with virulence of TBEV strains
Estimated divergence time of Him-TBEV
To understand the temporal constraints on the origin and dispersal of the TBEV group, BEAST was used to estimate the time of the most recent common ancestor (tMRCA) of each lineage, based on nucleotide sequences of E protein (1488 nt). Maximum clade credibility (MCC) tree was constructed with 40 representative strains of TBEV and related flaviviruses, including 5strains of Eu-TBEV, 5 strains of Sib-TBEV, 21 strains of FE-TBEV, 1 strain of Baikalian subtype and TBEV strain 178–79. The resulting tree is shown in Fig. . The most recent common ancestor of the Him-TBEV and the other subtypes of TBEV was estimated about 3119 (95% HPD = 1123–7695) years ago, when the western TBEV group (Eu-TBEV) and the eastern TBEV group were diverged. The Him-TBEV was first diverged with the eastern TBEV group about 2469 (95% HPD = 956–6376) years ago. The Sib-TBEV was formed about 800 years ago. The Baikalian subtype was diverged about 1820 years ago. Most nodes had 100% bootstrap support in the neighbor-joining tree. It should be mentioned that Chinese isolates belong to two subtypes. The strain of Sib-TBEV (Sib-XJ-X5) emerged about 480 years ago. The strains of FE-TBEV in china emerged about 210 years ago, which widely distributed the northeastern region of China (Fig. ).
The year of sampling, strain name and accession number are on the tip labels. Node labels indicate the most recent common ancestor (tMRCA) and the posterior probabilities. Strains of Him-TBEV were indicated in boldface. The TBEV strains detected in China were marked with black triangle
Discussion
We have discovered two novel strains (Himalaya-1 and Himalaya-2) of TBEV in Marmota himalayana for the first time in Qinghai-Tibet Plateau, China. Phylogenetic analysis, which has been widely used to study the genetic population of TBEV, clearly suggested that the Him-TBEV separated from other three subtypes (FE-TBEV, Sib-TBEV, and Eu-TBEV). The topology of phylogenetic trees demonstrated that the Him-TBEV could be considered as a new member of eastern TBEV group, together with FE-TBEV and Sib-TBEV (Fig. ). Classification of the Him-TBEV was also supported by comparative analysis of E protein and polyprotein. At amino acid level, the diversity of the E protein was less than 2.2% within the subtype and 3.6–5.6% between the subtypes, which is in the range of variation reported for other flaviviruses (3–10%).Citation5, Citation22–Citation25 Him-TBEV differed by 5.0–7.3% from the other subtypes. The same was observed for the polyprotein where 4.8–7.4% divergence calculated between Him-TBEV and the other subtypes was almost identical to the value 5.0–6.8%Citation8, the range of distances found between conventional three subtypes (Supplementary Table S1).
Classification of Him-TBEV as a new subtype was further supported by signature amino acid analysis. Notably, 69 amino acid substitutions in complete polyprotein were found in comparison with 112 strains of the three subtypes, of which 36 substitutions were unique for Him-TBEV, but the amino acids substitution in the same position were conserved in the three subtypes (Table ). For the so-called subtype-specific amino acid position 206 of E protein, Him-TBEV shared the Val with that of Eu-TBEV, but differed from that of FE-TBEV and Sib-TBEV, which suggested that the amino acid at the 206 position were not unique between subtypes of TBEV.Citation5
Although the classification of TBEV is mainly based on phylogenetic analysis and differences in nucleotide and amino acid sequences of E protein and polyprotein, isolation of a strain is required for full recognition of Him-TBEV as a new subtype.
With the linear and nonlinear regression analysis of genetic versus geographic distance combined with BEAST analysis, Heinze et al. reported that the most recent common ancestor of the TBEV group was present at 3119 year agoCitation26. When Him-TBEV was included and analyzed with similar methodology, it was observed that the Him-TBEV and the eastern TBEV group were diverged about 2469 years ago. After the divergence, the Sib-TBEV was separated from the FE-TBEV subtype about 2176 years ago (Fig. ).
The topology of our MCC tree was almost consistent with that reported by Heinze et al.Citation26. Those results suggested that Him-TBEV is a very ancient lineage of eastern TBEV (Figs. , ). When geographic site was considered, it is possible that after divergence with western TBEV, the Him-TBEV was separated with rest eastern TBEV subtypes about 2469 years ago, in Himalaya region (Fig. ). The FE-TBEV and Sib-TBEV were then progressed eastwards to Russian, and northeastern China.
The maximum clade credibility tree is plotted onto geographic location. The time of estimated tMRCA of each phylogenetic linage are indicated at the major nodes in the tree
It should be mentioned that the Qinghai-Tibet Plateau is the highest and one of the most extensive plateaus in the world, covering an area of 2.5 × 106 km2 with an average elevation of more than 4000 m. Qinghai-Tibetan Plateau of China is known to be the plague endemic region where marmot Marmota himalayana is the primary hostCitation27. Most human-inhabited areas at risk of exposure to enzootic plague are distributed in the east and south of the Plateau. There were 27 species of ticks detected in this region. The tick vectors and possible human infections of Him-TBEV should be urgently evaluatedCitation28.
Materials and methods
Sample collection
The animals of Marmota himalayana were sampled in Haixi prefecture (with an altitude of 2994 m above sea level [a.s.l]), in Qinghai province of China in July 2013. A total of 200 marmots were captured by cages in the field and sampled in the laboratory of the Centre for Disease Control and prevention of Haixi prefecture. The sampling was performed in accordance with the protocol for national plague surveillance program in animals. The low respiratory tissue and contents of each animal were separately collected in 2 ml sterile tubes with viral transport medium, which were then stored at −20 ℃ freezer in local laboratory immediately, transported in the same freezer to the laboratory in Beijing. The samples were transferred and storied in −80 °C freezer until use. The study has been reviewed and approved by the ethic committee of National Institute for Communicable Diseases Control and Prevention, China CDC.
RNA extraction and RNA-seq sequencing
Total RNA of low respiratory sample was extracted using RNeasy Fibrous Tissue Mini Kit (Qiagen, Hilden, Germany).The RNA was eluted in 50 μl RNase-Free Water (Qiagen, Hilden, Germany), which was used for RNA-seq library construction and as template for reverse transcriptase PCR (RT-PCR). The RNAs of low respiratory samples from 100 marmots were randomly chose and merged into 2 pools for RNA-seq library construction and sequencing. After DNase I digestion, total RNA was subjected to an rRNA removal step by using Ribo-Zero Magnetic Gold Kit (Human/Rat/Mouse) and Ribo-Zero Magnetic Gold Kit (Bacterial) (Epicentre, Madison, USA). The remaining RNA was used to construct RNA-seq library according to protocol provided. Briefly, mRNA was purified and enriched using oligo (dT) magnetic beads. The mRNA then was fragmented, and the cDNA was synthesized from the RNA fragments using reverse transcriptase and random primers. After the synthesis of the cDNA, ends repair was performed, followed by adenylation of the 3′ end, ligation of sequencing adapter and quantification using Aglient 2100 Bioanalyzer and ABI StepOnePlus Real-time PCR system. Pair-end (125 bp) sequencing was performed on Hiseq 2500 platform (Illumina, Sandiego, USA). The library construction and sequencing procedures were performed in BGI Tech (Shenzhen, China).
The resulting sequencing reads were trimmed and assembled de novo into contigs using the Trinity softwareCitation29. The assembled contigs were translated and compared to reference protein sequences of all RNA viruses by using local Blastx with an E-value of 1e−5. The targeted contigs were extracted by a perl script and compared to the entire non-redundant protein database to exclude non-viral sequences using online Blastx. The resulting viral sequences were merged by using SeqMan software within the Lasergene software package version 7.1 (DNAstar, USA) by identifying unassembled overlaps between neighboring contigs.Citation30 The newly generated viral sequences were performed a second round of online Blastx and the highly scoring hit was considered as the closet homolog.
RT-PCR for TBEV detection
For the detection of TBEV between the 200 marmots, a semi-nested RT-PCR method was performed by amplifying a 597-bp fragment spanning the prM-E-coding junction of the virus using specific primers designed based on nucleotide sequences obtained by RNA-seq sequencing.
The first round of semi-nested RT-PCR was performed using PrimeScript One Step RT-PCR kit (TaKaRa, Japan). The RT-PCR mixture (50 μl) contained 2 μl of total RNA, 2 μl PrimeScript one Step Enzyme Mix and 1 μl (20 μM) of each outer primers (F1: 5′-GACTCACTGTCCTATGAGTG-3′, R1:5′-GACCTCCATGACCACTGTGTCAT-3′). After Reverse transcription at 50 °C for 30 min and initial denaturation at 94 °C for 2 min, the mixtures were amplified with 30 cycles of 94 °C for 30 s, 50 °C for 81 s, and 72 °C for 30 s and a final extension at 72 °C for 10 min. Subsequently, 5 μl of the first-round products was used as the template of the second round of PCR amplification. In this round, 20 μl of PCR mixture included 2.5U ExTaq DNA polymerase (TaKaRa, Japan) and 1 μl (20 μM) of each inner primers (F2: 5′-TATGATGCCAACAAGATCGT-3′, R1). After initial denaturation at 94 °C for 5 min, the mixtures were amplified with 30 cycles of 94 °C for 30 s, 50 °C for 36 s, and 72 °C for 30 s and a final extension at 72 °C for 10 min.
Products of second-round PCR were gel purified using a QIAquick gel extraction kit (Qiagen, Hilden, Germany) and sequenced using the ABI prism 3700 DNA Analyzer (Applied Biosystems, Foster City, CA. USA) by Sanger method. The sequences of PCR products were compared with the known sequences of TBEV in the GenBank database.
Full genome sequencing
The complete genomes of TBEV Himalaya-1 and Himalaya-2 were amplified and sequenced using the RNA extracted from the original sample as templates. The RNA was reverse transcripted and amplified by 11 sets of specific primers (Supplementary Table S3) designed based on multiple alignment of contigs obtained by RNA-seq sequencing and genome sequences of TBEV strains used in Supplementary Table S1. Each amplified segments were sequenced for three times. Genomic Sequences were assembled using SeqMan software (DNAstar, USA).
Genome analysis
The nucleotide sequences of the genomes and the deduced amino acid sequences of the open reading frame (ORF) were compared to those of other TBEV strains with complete genome. Complete nucleotide sequences of polyprotein and E protein were aligned with clustalW within the BioEdit software (version 7.1). Phylogenetic analysis was conducted using the neighbor-joining method using MEGA 6.0 with 1000 bootstrap replications. Virus strains used in this study were listed in Supplementary Table S2.
Evolutionary analysis
Complete nucleotide sequences of E protein with known year of collection were used to estimate the mean time of the most recent common ancestor (tMRCA) using the Bayesian Markov Chain Monte Carlo (MCMC) approach employed by BEAST package Version.Citation31 The jModeltest software 2.1.7 was usedCitation32 to estimate the best-fit nucleotide substitution model according to the Akaike information criterion (AIC), with GTR + I + G as the best substitution model. Analyses were conducted under the best-fit nucleotide substitution model and using a relaxed (uncorrelated lognormal) molecular clock model. The MCMC analysis was performed with 20 million generations and was sampled every 1000 generations with 10% burn-in. Convergence of parameters was assessed on the basis of the ESS reaching values >200 by using Tracer software version 1.5. Maximum clade credibility (MCC) trees were subsequently generated after 10% burn-in using Tree Annotator and viewed by FigTree.
Supplementary material for revised manuscript
Download MS Word (250.5 KB)Acknowledgements
This work was supported by The National Key Research and Development Plan, 2016YFC1201903), National Natural Science Foundation of China (81290340 and 81290345), and the Ministry of Science and Technology, People′s Republic of China (Mega Project of Research on The Prevention and Control of HIV/AIDS, Viral Hepatitis Infectious Diseases, 2013ZX10004-101).
Conflict of interest
The authors declare that they have no conflict of interest.
Electronic supplementary material
Supplementary Information accompanies this paper at (10.1038/s41426-018-0081-6).
Related Research Data
References
- LindquistLVapalahtiOTick-borne encephalitisLancet2008371 1861 187110.1016/S0140-6736(08)60800-4
- GritsunTSLashkevichVAGouldEATick-borne encephalitisAntivir. Res.20035712914610.1016/S0166-3542(02)00206-1
- BogovicPStrleFTick-borne encephalitis: a review of epidemiology, clinical characteristics, and managementWorld J. Clin. Cases2015343044110.12998/wjcc.v3.i5.4304419106
- ChambersTJFlavivirus genome organization, expression, and replicationAnnu. Rev. Microbiol.19904464968810.1146/annurev.mi.44.100190.003245
- EckerMSequence analysis and genetic classification of tick-borne encephalitis viruses from Europe and AsiaJ. Gen. Virol.19998017918510.1099/0022-1317-80-1-179
- GrardGGenetic characterization of tick-borne flaviviruses: new insights into evolution, pathogenetic determinants and taxonomyVirology2007361809210.1016/j.virol.2006.09.015
- KovalevSYMukhachevaTAReconsidering the classification of tick-borne encephalitis virus within the Siberian subtype gives new insights into its evolutionary historyInfect. Genet. Evol.20175515916510.1016/j.meegid.2017.09.014
- DeminaTVGenotyping and characterization of the geographical distribution of tick-borne encephalitis virus variants with a set of molecular probesJ. Med. Virol.20108296597610.1002/jmv.21765
- LiuYJEpidemiological characteristics of tick borne encephalitis natural foci in JilinBull. Acad. Mil. Med. Sci. [Zhongguo Ren. Min. Jie Fang. Jun. Jun. Shi Yi Xue Ke Xue Yuan Yuan Kan.]19793109120
- YinDMLiuRZReview on the control of forest encephalitis in the forest areas of north-east ChinaChin. J. Epidemiol. [Zhonghua Liu Xing Bing. Xue Za Zhi]200021387389
- LuZBrokerMLiangGTick-borne encephalitis in mainland ChinaVector Borne Zoonotic Dis.2008871372010.1089/vbz.2008.0028
- SunRXMapping the distribution of tick-borne encephalitis in mainland ChinaTicks Tick Borne Dis.2017863163910.1016/j.ttbdis.2017.04.009
- LiuRDetection of the Siberian Tick-borne Encephalitis Virus in the Xinjiang Uygur Autonomous Region, northwestern ChinaChin. J. Virol. [Bing. du xue bao]2016322631
- ZhangGLInvestigation on the endemic foci of new emerged tick-borne encephalitis in Charles Hilary, XinjiangChin. J. Epidemiol. [Zhonghua liu xing Bing. xue za zhi]201334438442
- HeinzFXEmergence of tick-borne encephalitis in new endemic areas in Austria: 42 years of surveillanceEur. Surveill.20152091610.2807/1560-7917.ES2015.20.13.21077
- YuJMA novel hepatovirus identified in wild woodchuck Marmota himalayanaSci. Rep.2016610.1038/srep223614770319
- GritsunTSHolmesECGouldEAAnalysis of flavivirus envelope proteins reveals variable domains that reflect their antigenicity and may determine their pathogenesisVirus Res.19953530732110.1016/0168-1702(94)00090-Y
- ReyFAThe envelope glycoprotein from tick-borne encephalitis virus at 2 A resolutionNature199537529129810.1038/375291a0
- MandlCWAdaptation of tick-borne encephalitis virus to BHK-21 cells results in the formation of multiple heparan sulfate binding sites in the envelope protein and attenuation in vivoJ. Virol.2001755627563710.1128/JVI.75.12.5627-5637.2001114275
- BelikovSIThe relationship between the structure of the tick-borne encephalitis virus strains and their pathogenic propertiesPLoS One20149e9494610.1371/journal.pone.00949463989262
- LeonovaGNComprehensive assessment of the genetics and virulence of tick-borne encephalitis virus strains isolated from patients with inapparent and clinical forms of the infection in the Russian Far EastVirology2013443899810.1016/j.virol.2013.04.029
- ChuMCO’RourkeEJTrentDWGenetic relatedness among structural protein genes of dengue 1 virus strainsJ. Gen. Virol.1989701701171210.1099/0022-1317-70-7-1701
- NiHBarrettADNucleotide and deduced amino acid sequence of the structural protein genes of Japanese encephalitis viruses from different geographical locationsJ. Gen. Virol.19957640140710.1099/0022-1317-76-2-401
- LepiniecLGeographic distribution and evolution of yellow fever viruses based on direct sequencing of genomic cDNA fragmentsJ. Gen. Virol.19947541742310.1099/0022-1317-75-2-417
- LewisJAPhylogenetic relationships of dengue-2 virusesVirology199319721622410.1006/viro.1993.1582
- HeinzeDMGouldEAForresterNLRevisiting the clinal concept of evolution and dispersal for the tick-borne flaviviruses by using phylogenetic and biogeographic analysesJ. Virol.2012868663867110.1128/JVI.01013-123421761
- QianQMapping risk of plague in Qinghai-Tibetan Plateau, ChinaBmc Infect. Dis.20141410.1186/1471-2334-14-3824227279
- YangYSInvestigation on kinds and nature geographic distribution of ticks in Qinghai provinceChin. J. Hyg. Equip. [Zhonghua Wei Sheng Sha Chong Yao Xie]200814201203
- GrabherrMGFull-length transcriptome assembly from RNA-Seq data without a reference genomeNat. Biotechnol.20112964465210.1038/nbt.18833571712
- LiCXUnprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruseseLife20154e053784384744
- DrummondAJRambautABEAST: Bayesian evolutionary analysis by sampling treesBmc Evol. Biol.2007721410.1186/1471-2148-7-2142247476
- PosadaDjModelTest: phylogenetic model averagingMol. Biol. Evol.2008251253125610.1093/molbev/msn083