Complete Genome Characterization of the 2017 Dengue Outbreak in Xishuangbanna, a Border City of China, Burma and Laos
Songjiao Wen1,2,3† 
Dehong Ma4† 
Yao Lin1,2,3 
Lihua Li4 Shan Hong1,2,5 Xiaoman Li6 
Xiaodan Wang1,2,3 
Juemin Xi1,2,3 Lijuan Qiu6 Yue Pan1,2,3 Junying Chen1,2,3 
Xiyun Shan4* 
Qiangming Sun1,2,3*- 1Institute of Medical Biology, Peking Union Medical College, Chinese Academy of Medical Sciences, Kunming, China
- 2Yunnan Key Laboratory of Vaccine Research and Development on Severe Infectious Diseases, Kunming, China
- 3Yunnan Key Laboratory of Vector-borne Infectious Disease, Kunming, China
- 4Xishuangbanna Dai Autonomous Prefecture People's Hospital, Xishuangbanna, China
- 5School of Basic Medicine, Kunming Medical University, Kunming, China
- 6Institute of Pediatric Disease Research, The Affiliated Children's Hospital of Kunming Medical University, Kunming, China
A dengue outbreak abruptly occurred at the border of China, Myanmar, and Laos in June 2017. By November 3rd 2017, 1184 infected individuals were confirmed as NS1-positivein Xishuangbanna, a city located at the border. To verify the causative agent, complete genome information was obtained through PCR and sequencing based on the viral RNAs extracted from patient samples. Phylogenetic trees were constructed by the maximum likelihood method (MEGA 6.0). Nucleotide and amino acid substitutions were analyzed by BioEdit, followed by RNA secondary structure prediction of untranslated regions (UTRs) and protein secondary structure prediction in coding sequences (CDSs). Strains YN2, YN17741, and YN176272 were isolated from local residents. Stains MY21 and MY22 were isolated from Burmese travelers. The complete genome sequences of the five isolates were 10,735 nucleotides in length. Phylogenetic analysis classified all five isolates as genotype I of DENV-1, while isolates of local residents and Burmese travelers belonged to different branches. The three locally isolates were most similar to the Dongguan strain in 2011, and the other two isolates from Burmese travelers were most similar to the Laos strain in 2008. Twenty-four amino acid substitutions were important in eight evolutionary tree branches. Comparison with DENV-1SS revealed 658 base substitutions in the local isolates, except for two mutations exclusive to YN17741, resulting in 87 synonymous mutations. Compared with the local isolates, 52 amino acid mutations occurred in the CDS of two isolates from Burmese travelers. Comparing MY21 with MY22, 17 amino acid mutations were observed, all these mutations occurred in the CDS of non-structured proteins (two in NS1, 10 in NS2, two in NS3, three in NS5). Secondary structure prediction revealed 46 changes in the potential nucleotide and protein binding sites of the CDSs in local isolates. RNA secondary structure prediction also showed base changes in the 3′UTR of local isolates, leading to two significant changes in the RNA secondary structure. To our knowledge, this study is the first complete genome analysis of isolates from the 2017 dengue outbreak that occurred at the border areas of China, Burma, and Laos.
Introduction
Dengue virus (DENV) is a member of the family Flaviviridae and genus Flavivirus, which has four closely related serotypes (DENV-1 to DENV-4) (Anez et al., 2017). Among these serotypes, DENV-1 has five distinct genotypes (I–V) (Sasmono et al., 2015). DENV is a single-stranded positive-sense RNA virus approximately 11 kb in length, comprising an open reading frame (ORF) and 5′ and 3′untranslated regions (UTRs). Translation of the ORF encodes three structural proteins (capsid [C], membrane [PrM/M], and envelope [E]) and seven non-structural proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) (Nukui et al., 2006; Parameswaran et al., 2017).
Dengue is a highly prevalent disease spread by the mosquito species Aedes albopictus and Ae. aegypti (Lustig et al., 2017). This virus is endemic to more than 100 countries worldwide, especially in cities and semi-urban areas of tropical and subtropical countries (Imrie et al., 2010; Anez et al., 2017). A recent study estimated that 390 million cases of dengue infection occur worldwide each year, representing a 30-fold increase in the past 50 years (Anez et al., 2017). Dengue virus has become a worldwide public problem (Santos-Sanz et al., 2014). The World Health Organization (WHO) classified dengue disease as dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) in 1997 (Narvaez et al., 2011). However, in 2008, the WHO proposed a new classification system, classifying dengue into dengue without warning signs (DWWS), dengue with warning signs (DWS), and severe dengue (SD) (Lima et al., 2013).
The first dengue outbreak in China was reported in Foshan, Guangdong Province, in 1978. Since then, dengue epidemics have been reported in Guangdong, Jiangsu, Hainan, Fujian, Guangxi, Zhejiang, and Henan (Huang et al., 2014; Liu et al., 2014; Xiong and Chen, 2014; Chen and Liu, 2015; Ren et al., 2015). In 2008, an unexpected dengue outbreak occurred in Ruili, a city on the border of China and Burma (Zhang et al., 2013). In 2013 and 2015, large-scale dengue outbreaks occurred in Xishuangbanna, Yunnan Province, a Chinese city on the border of China, Burma, and Laos. In 2014 and 2016, sporadic cases of dengue fever were also reported. In June 2017, a dengue outbreak developed abruptly at the borders of China, Burma, and Laos. By 03 November 2017, 1,184 infected individuals in Xishuangbanna had been confirmed as dengue NS1-positive.
Materials and Methods
Ethics Statement
Each participant was informed of the purpose of the study. Written informed consent was obtained from all participants. The study protocol was approved by the Institutional Ethics Committee (Institute of Medical Biology, Chinese Academy of Medical Sciences, and Peking Union Medical College). The study was conducted in accordance with the Declaration of Helsinki for Human Research of 1974 (last modified in 2000).
Dengue Virus RNA Extraction and Identification
Serum samples were collected from DENV NS1-positive patients at Xishuangbanna Dai Autonomous Prefecture People's Hospital. The serum samples were separated from the collected blood, followed by viral RNA extraction. Viral RNA was extracted from 140 μl of dengue virus-infected serum using the QIAamp viral RNA mini kit (Qiagen, Hilden, Germany; No. 52906). The RNA was eluted in 60 μl of nuclease-free water and stored at −80°C. Viral RNA was reverse transcribed into cDNA using the PrimeScript™ II 1st Strand cDNA Synthesis Kit (Takara Bio, Shiga, Japan; No. 6210A). The cDNA was amplified by using dengue universal, 1, 2, 3, 4-type primers (Supplementary Table 1) for specific identification. The reaction conditions in each 25-μl volume were denaturation at 94°C for 5 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, and elongation at 72°C for 30 s; with a final elongation step at 72°C for 7 min. The PCR products were then subjected to gel electrophoresis.
Primer Design
A total of 21 synthetic oligonucleotide primer pairs (Supplementary Table 2) were designed to amplify overlapping fragments with sizes between 200 and 900 nucleotides (nt) spanning the entire viral genome of DENV-1. All primers were selected online from Primer-BLAST at the NCBI website, based on the strain WestPac (GenBank accession no. U88535.1). All primers were synthesized and purified by Shuo GE Biotechnology Co., Ltd. (Kunming, China).
Genome Synthesis for Amplification and Sequencing
PCR was performed with the following protocol (50 μl volume): denaturation at 98°C for 2 min, followed by 40 cycles of denaturation at 98°C for 10 s, annealing at 56°C for 10 s, and elongation at 72°C for 30 s; with a final elongation step at 72°C for 1 min. The PCR products were confirmed by agarose gel electrophoresis and sequenced at Shuo GE Biotechnology Co. Ltd. Forward and reverse sequencing were done.
Analysis of Molecular Characteristics
A total of 21 nucleotide sequences were assembled using BioEdit7.1.3 (http://www.mbio.ncsu.edu/bioedit/bioedit.html). Complete genome sequences of five DENV-1 isolates were submitted to the NCBI GenBank database with accession numbers of MF681692 (for YN2), MF683116 (YN1774), MF681693 (YN176272), MG679800 (MY21), and MG679801 (MY22). Nucleotide sequences and translated amino acid sequences were analyzed by BioEdit7.1.3. The secondary structure of isolates and reference strains CDS and UTR was predicted with the Predict Protein server (https://www.predictprotein.org/) and RNA fold web server (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi), respectively.
Phylogenetic Analysis
The sequences were assembled using BioEdit 7.1.3. Phylogenetic analysis based on the whole genome genes was conducted using the Molecular Evolutionary Genetics Analysis (MEGA) software version 6.0 (Maximum Likelihood method (ML), Bootstrap1000, Jones-Taylor-Thornton (JTT) model). The reference viral sequences used to construct the distinct phylogenetic branches were collected from the GenBank sequence database under the following country accession numbers: China (JQ048541, KP7752, KU365900, DQ193572, FJ176780, AB608788, KU094071, and KT827374), Singapore (GU370049, KJ806945, M87512, JN544411, and GQ357692), Malaysia (KX452068 and KX452065), South Korea (KP406802), Indonesia (KC762646), Japan (AB178040), Thailand (HM469966 and AF180817), Myanmar (AY726554 and AY726553), Laos (KC172829, KC172835, and KC172834), Cambodia (HQ624984), Viet Nam (JQ045660), Hawaii (EU848545, DQ672564, and DQ672562), Western Pacific (U88535), South Pacific (JQ915076), Brazil (AF226685), Haiti (KT279761), USA (FJ390379), India (KF289072), and standard strains (M29095, M93130, and AF326573).
Results
Geographic Analysis of the Outbreak
Xishuangbanna has become the central area of dengue fever epidemics occurring at the borders of China, Myanmar, and Laos. The latest was an abrupt onset outbreak that began in June 2017. As of November 3rd 2017, 1,184 infected individuals had been confirmed as dengue NS1-positive as of November 3rd in Xishuangbanna, Yunnan, China. These cases included 71 cases involving non-residents of Xishuangbanna and 1,113 local cases.
Base Sequence Analysis of the Complete Genome Sequences of Five DENV Isolates
The complete genome sequences of five isolates were obtained by amplifying 21 overlapping fragments. Two of the five complete genome sequences were obtained from travelers from Burma. The remaining three genome sequences were obtained from local residents. Sequencing and clustering, All isolates were 10,735 nt in length. The 5′ and3′UTRs of the five isolates were 94 and 465 nt in length, respectively. The ORF of all five isolates was located between nt 95 and 10270, and encoded 3,392 amino acids.
Phylogenetic Analysis
The whole genome of the five isolates was aligned and compared with 40 sequences from standard strains of four serotypes and some typical DENV-1 strains of various geographical origins obtained from the NCBI GenBank database. The phylogenetic analysis was constructed with the ML method. The five isolates belonged to DENV-1 and were further classified as genotype. Three local isolates were located in one cluster of the ML tree, with the closest relationship to the Dongguan strain of China (2011, JQ048541). Additionally, three isolates were also closely related to the following strains: China (Hubei 2014 KP772252, Taiwan 2014 KU365900, and Fujian 2005 DQ193572), Singapore (2008, GU370049), and Malaysia (2014 KX452068 and KX452065). However, two isolates obtained from Burmese travelers were located in another cluster of the ML tree and were most similar to the Laos strain in 2008 (KC172834). DENV-3 (M93130), DENV-4 (AF326573), and DENV-2 (M29095) were outgroups (Figure 1A).
Figure 1. Phylogenetic tree data. (A) Phylogenetic tree of DENV-1 epidemic strains in Xishuangbanna, Yunnan, China, in 2017. The phylogenetic tree was constructed with the Maximum Likelihood method with 1,000 bootstrap values. • isolates from Xishuangbanna, Yunnan, China; ▴ standard strains of DENV1-4 serotypes. (B) Relationship between evolutionary tree branching and amino acid mutations in genotype I of DENV-1.
We further analyzed the relationship between the evolutionary tree branching and amino acid substitution of genotype ? of DENV-1. A total of 24 amino acid substitutions were detected that had important roles in eight evolutionary treebranches. Particularly, A1843S, I2392V, D2921E, S3372L, and E3386D resulted in branch 1; N75S, K1345R, R2871K, K3044R, V3280I, and E3326G resulted in branch 2; I781V, K945R, M1172I, C1184S, V2264D, A2738V, and Y3142H resulted in branch 3; and S3059A, I1937V/A2863T, N90S, H2620Y, and I2628M resulted in branches 4–8, respectively (Figures 1B, 2).
Figure 2. Identification of variations among genotype I of DENV-1. Variations are presented in accordance with their location in the viral genome. The various colors indicate the different variations described in Figure 1B.
Base Substitution and Amino Acid Mutation Analysis of Coding DNA Sequences (CDSS) From the Isolated Strains
Base substitution and amino acid mutation analyses were performed by comparison with the DENV-1 standard strain (DENV-1SS) (Hawaii, 1944, NO. EU848545). The results revealed 187 base substitutions in the structural protein (C/PrM/E) in three local isolates, of which 155 substitutions were non-synonymous. These were 471 base substitutionsin non-structural protein regions in the three local isolates, leading to 416 non-synonymous mutations (Supplementary Figure 1). The genome sequence of the three local isolates was almost the same, except for nt 3786(C → T) and 10630(A → G) in isolate YN17741. All the isolates harbored a base deletion (T) at nt 79 of the 5′UTR, and a few base substitutions occurred in the 3′UTRs in all the isolates.
The two isolates obtained from Burmese travelers were different from the three local isolates. Compared with local isolates, a total of 52 amino acid mutations occurred in the CDS of the two isolates (MY21 and 22) obtained from Burmese travelers. There were also differences between MY21 and 22. Comparison of MY21with MY22 revealed 17 amino acid mutations, which all occurred in CDS of non-structural proteins. Of these, two amino acid mutations (R874Q, D995V) occurred in the non-structural protein NS1, three amino acid mutations (A1317T, V1331L, G1340K) occurred in the non-structural protein NS2A, seven amino acid mutations (D1350N, G1355A, G1364S, A1365S, E1368K, G1369N, G1374A) occurred in the non-structural protein NS2B, two amino acid mutations (D1960S, I2096V) occurred in the non-structural protein NS3, and three amino acid mutations (H2778Y, S3122L, S3239G) occurred in the non-structural protein NS5 (Figure 3).
Figure 3. Amino acid mutations between CDS of locally isolated strains and two strains isolated from Myanmar (2017).
Compared to the closest Dongguan strain (China 2011, JQ048541), a total of 14 amino acid mutations occurred in the CDS of the local isolates (Figure 4). Three mutations leading to changes in amino acid polarity occurred in non-structural protein regions. At amino acid 1031in L1031Y, the amino acid changed from a non-polar to a polar residue, while the change was from polar to non-polar residues in positions S1183L and V1615L).
Figure 4. Amino acid mutations between CDS of Dongguan strain (2011) and Yunnan strains (2017).
Prediction of the RNA Secondary Structure Based on the Untranslated Regions of the Standard and Local Isolated Strains
All sequences wereconserved in the UTRs. Compared to the standard strain, only one base deletion occurred, and thus the predicted RNA structure of the 5′UTR was only marginally changed (Supplementary Figures 2A,B). Comparison with the 3′UTR of the standard strains revealed that all local isolates had eight substitutions (nt 10273 [A → G], nt 10276 [T → C], nt 10286 [C → T], nt 10293 [C → T], nt 10408 [G → A], nt 10467 [C → T], nt 10536 [A → G], nt 10621 [C → T]) and 3 inserts (nt 10455 [G], nt 10479 [A], and nt 10730 [C]). In addition to these mutations, an additional mutation at nt 10630 [A → G] was observed in YN17741. RNA secondary structure prediction results showed that some base changes in the 3′UTR of the local isolates lead to two significant changes in the RNA secondary structure. In contrast to the long interior loop in the standard strain, a 47 bp stem downstream of nt 10277 was formed in the isolates. Long stem and loop structures were formed before nt 10734 in the standard strain, while the stem and loops were shorter in the isolates before nt 10719 (Supplementary Figures 2C,D).
Potential Secondary Structure of the CDS Regions
The potential secondary structures of the CDS regions of the local isolates were predicted and compared with that of the standard strain. In structural proteins C/prM/E, there were 13 changes inputative nucleotide/protein binding sites among the 775 amino acids. Among these, six putative protein-binding sites (sites 8, 68, 87, 161, 335, and 436) were detected only in the isolates. However, five potential protein-binding sites (sites 1, 7, 21, 86, and 429) and two nt-binding sites (sites 6 and 7) were lost compared to the standard strain. Approximately eight distinct changes appeared in the helix and strand regions. Two distinct changes were found in the exposed and buried regions and the disordered regions. The helical transmembrane regions in the local isolates were highly consistent with those in the standard strain (Supplementary Figure 3). In the non-structural protein NS1, comparison with the standard strain revealed only one potential nucleotidesite (site 118) and protein binding site (site 99) in the isolates. Approximately four distinct changes were found in the helix and strand regions (Supplementary Figure 3B). In the non-structural protein NS2A, three changes were observed in putative protein binding sites among the 218 amino acids. Among these changes, four potential protein-binding sites were increased, and two other sites (118 and 180) were lost in the isolated strains (Supplementary Figure 3C). In the non-structural protein NS2B, comparison with reference strain revealed only one change (site 49) in the putative protein binding sites of the isolates (Supplementary Figure 3D). In the non-structural protein NS3, 15 changes were observed in the putative nucleotide and protein binding sites among the 619 amino acids. Among these changes, five potential protein-binding sites (sites 132, 134, 143, 338, and 471) were increased in the isolates, while 10 putative protein-binding sites (sites 107, 161, 172, 304, 461, 481, 483, 531, 539, and 593) were lost compared with the standard strain (Supplementary Figure 3E). In the non-structural protein NS4A, only two distinct changes were found in the exposed and buried regions and disordered regions (Supplementary Figure 3F). In the non-structural protein NS4B, the isolates lost potential nucleotide (site 249) and protein binding sites (site 243) compared to the standard strain (Supplementary Figure 3G). In the non-structural protein NS5, 11 changes were detected in putative nucleotide and protein binding sites among the 899 amino acids. Four putative protein-binding sites (sites 51, 289, 558, and 691) and two nucleotide-binding sites (sites 598, 853) were increased in the isolated strains, while five putative protein-binding sites (sites 350, 465, 499, 502, and 895) were lost. In addition, three distinct changes were found in the disordered regions (Supplementary Figure 3H).
Discussion
Yunnan Province is located on the southwest border of China. The northwest corner is closely attached to the Tibet autonomous region. The west borders Myanmar and the south is connected with Laos and Vietnam. Yunnan is located in the vast expanse of the Asian continent, south of the Indian Ocean and the Pacific Ocean. The region is influenced by monsoons that originate to the southeast and southwest, and also by the Tibetan plateau area. These climatic influences have produced a complex and diverse natural geographical environment. The region is a popular touristic destination. These and other factors have made this area important for the prevention and control of dengue fever in Asia.
A 2006 study revealed the existence of DENV infection in the Dushan and Xingyi areas of the Yunnan-Guizhou Plateau, China (Liu et al., 2006). In 2008, 48 imported cases and one native case were reported in Yunnan Ruili. This was the first report of a dengue outbreak in Yunnan. Of the 49 confirmed cases, 18 involved people from Mujie City, Myanmar, 30 cases involved Chinese residents who had returned to Myanmar from elsewhere, and one native case. This outbreak was caused by DENV-1 and DENV-3 (Zhang et al., 2013). A DENV-2 isolate was obtained from a patient with fever who returned from Laos in 2009 (Guo et al., 2013). The first aboriginal outbreak occurred in Xishuangbanna Dai Autonomous Prefecture, Yunnan Province, China, in 2013. This outbreak was caused by DENV-3 and 1,331 cases were confirmed. A mixed outbreak of DENV-1 and DENV-2 appeared in Ruili in Dehong prefecture, with 246 dengue fever cases recorded (Zhang et al., 2014; Guo et al., 2015; Wang et al., 2017). From June to December of 2014, a total of 292 cases of dengue fever were diagnosed in Ruili City. These included 139 local cases and 153 imported cases caused by DENV-1, DENV-2, and DENV-4. In 2015, 240 DENV-1 cases were reported in Lincang, and 1,067 DENV-2 cases were reported in Xishuangbanna (Guo et al., 2016; Zhao et al., 2016). In 2016, 25 isolates were obtained from the sera of dengue fever patients. These included 15 local cases and 10 imported cases from Burma. These cases were caused by DENV-1, DENV-2, DENV-3, and DENV-4. A dengue fever outbreak caused by DENV-1occurred in Xishuangbanna, Yunnan Province, in 2017. This was the first known outbreak due to DENV-1in this area. According to government statistics, as of 03 November 2017, a total of 1,184 infected cases have been verified, including 71 imported cases and 1,113 local cases. Interestingly, large-scale dengue outbreak shave occurred in Xishuangbanna every 2 years since 2013, involving different serum types of DENV.
In the present study, viral RNA was extracted from the sera of patients from China and Myanmar, followed by amplification of overlapping fragments to obtain five complete genomes (YN2, YN17741, YN176272, MY21, and MY22). Two isolates (MY21 and MY22) were obtained from Burmese travelers. The five isolates were 10,735 bp in length. YN2 and YN176272 displayed the same genome sequence. Compared with YN2, YN17741 possessed has two base changes at nt 3786 (C → T) and nt 10630 (A → G). These base changes were non-synonymous. However, the two isolates from Burmese travelers were different from three isolated strains from local residents. Compared with locally isolated strains, a total of 52 amino acid mutations occurred in the CDS of the two isolates from Burmese travelers. These two isolates differed from each other. Comparison of MY21 and MY22 revealed 17 amino acid mutations. All occurred in CDS of non-structured proteins. These results indicate that the 2017 dengue outbreak in Xishuangbanna was not caused by the imported DENV-1, but rather was a local outbreak. The 2017 dengue outbreak in the border areas around Xishuangbanna was not caused by only one DENV-1 source strain, but belong to a single genotype.
The 5′ and 3′ UTRs of dengue virus play an important role in the genome replication of RNA (Wei et al., 2008; Vashist et al., 2009). According to different conservative levels, the 3′ UTR is divided into three regions: (1) a variable area directly following the ORF; (2) a highly conserved cyclization sequence motif (CS1) and stable stem-loop structure; and (3) the region between the variable and 3-terminal regions, which have modest conservative levels and contain several hairpin motifs (Bryant et al., 2005). The 3′UTR is dispensable for the replication of dengue type 1 virus in vitro (Tajima et al., 2006). The variable area in the 3' UTR of Japanese encephalitis virus (JEV) reportedly did not affect its growth in vitro or its pathogenicity in mice (Kato et al., 2012). However, a previous report showed that the variable area of DENV-2 can enhance viral replication in baby hamster kidney cells, although this area is dispensable in mosquito cells (Alvarez et al., 2005). Another study showed that VR mutations of DENV-1 might not affect the translation process (Tajima et al., 2007). In the present study, mutations were observed in the 3′UTRs of the isolates compared to those of the reference strains. Nucleotide changes in 3′UTRs of the isolates lead to two significant changes in the RNA secondary structure. In contrast to the long interior loop in the standard strain, a 47-bp stem downstream of nt 10277 was formed in the isolated strains. Long stem and loop structures were formed before nt 10734 in the standard strain, while the stem and loops were shorter in the isolated strains before nt 10719, especially in the variable area. The effect of these changes on dengue type 1 replication needs further clarification.
In recent years, several studies reported that the mutation of amino acids in some special sites has a significant impact on the virulence, replication, and infectivity of flaviviruses. An alanine-to-valine amino acid substitution at residue 188 of the NS1 protein promotes the infectivity of Zika virus (ZIKV) Asian genotype in mosquitoes (Delatorre et al., 2017; Liu et al., 2017). A single mutation (S139N) in the PrM protein of Zika virus promotes fetal microcephaly (Yuan et al., 2017). The single mutation G66A in the NS2A gene of Japanese encephalitis-live vaccine virus reduces neurovirulence and neuroinvasiveness in mice(Ye et al., 2012). The mutation of putative N-glycosylation sites (Asn-58 and Asn-62) affects the biological function of DENV NS4B in the viral replication complex (Naik and Wu, 2015). A DENV-2 vaccine strain (PDK53) has three mutations that are critical for the attenuated state of PDK53 (Butrapet et al., 2000). The present results indicate that the region of the viral structural protein is relatively conserved, with less amino acid substitutions, with mutations being more prevalent in the non-structural protein regions in genotype I of DENV-1. Based on phylogenetic analysis, we identified a total of 24 amino acid substitutions that play an important role in evolutionary tree branching. These substitutions included two important amino acid substitutions (N90S and S3059A) in Chinese strains (DQ193572, KU365900, KP772252, JQ048541, and MF681692), and an amino acid substitution (I2628M) in the local isolated strains and the Dongguan strain of China (JQ048541). The effects of these amino acid substitutions on virulence, replication, infection, and other functions of the virus are unclear.
In 2017, there have been extensive outbreaks of dengue fever in many areas including Yunnan, Guangdong, Zhejiang, Shandong, Taiwan, Xianggang, Macao, and other areas (Chen et al., 2017). Xishuangbanna is geographically near Myanmar, Laos, and Vietnam. The present study is the first report on the complete genome characterization of the 2017 Xishuangbanna DENV outbreak. The present study should serve as a reference to follow-up studies of DENV outbreaks in Southeast Asia in 2017. The findings are a reference for future research on the virulence, replication, infection, pathogenicity, and vaccine development of DENV.
Author Contributions
XS and QS contributed to the design of the study. SW and DM were involved in data acquisition and analysis. YL, LL, SH, and XL provided software technology support. XW, JX, YP, LQ, and JC provided experiment technology support. All the authors agree to the final version of the submission.
Funding
This research was supported by the Foundation of the CAMS Initiative for Innovative Medicine (CAMS-I2M) (grant no. 2016-I2M-1-19), the Natural Science Foundation of Yunnan Province (2016FA029), Yunnan Medical Discipline Leader Project (D-201605), Kunming Science and Technology Innovation Team Project (grant no. 2015-1-R-00956), the National Twelfth Five Year major new drug discovery technology major project (grant no. 2012ZX09104-302), the Jointly Supported Foundation of the National Project in Yunnan Province (2013GA018), and the National Natural Science Foundation of China (81171946).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank Dr. Penghua Wang of Department of Microbiology and Immunology, School of Medicine, New York Medical College for proofreading and reviewing our manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2018.00148/full#supplementary-material
References
Alvarez, D. E., De Lella Ezcurra, A. L., Fucito, S., and Gamarnik, A. V. (2005). Role of RNA structures present at the 3'UTR of dengue virus on translation, RNA synthesis, and viral replication. Virology 339, 200–212. doi: 10.1016/j.virol.2005.06.009
Añez, G., Volkova, E., Jiang, Z., Heisey, D. A. R., Chancey, C., and Fares, R. C. G. (2017). Collaborative study to establish World Health Organization international reference reagents for dengue virus Types 1 to 4 RNA for use in nucleic acid testing. Transfusion 57, 1977–1987. doi: 10.1111/trf.14130
Bryant, J. E., Vasconcelos, P. F., Rijnbrand, R. C., Mutebi, J. P., Higgs, S., and Barrett, A. D. (2005). Size heterogeneity in the 3' noncoding region of South American isolates of yellow fever virus. J. Virol. 79, 3807–3821. doi: 10.1128/JVI.79.6.3807-3821.2005
Butrapet, S., Huang, C. Y., Pierro, D. J., Bhamarapravati, N., Gubler, D. J., and Kinney, R. M. (2000). Attenuation markers of a candidate dengue type 2 vaccine virus, strain 16681 (PDK-53), are defined by mutations in the 5' noncoding region and nonstructural proteins 1 and 3. J. Virol. 74, 3011–3019. doi: 10.1128/JVI.74.7.3011-3019.2000
Chen, B., and Liu, Q. (2015). Dengue fever in China. Lancet 385, 1621–1622. doi: 10.1016/S0140-6736(15)60793-0
Chen, Q., Song, W., Mou, D., Li, Y., Yin, W., Li, Z., et al. (2017). Analysis on epidemiological characteristics of dengue fever in China, as of 31th August,2017. Dis. surveill. 10, 801–804. doi: 10.3784/j.issn.1003-9961.2017.10/11.005
Delatorre, E., Mir, D., and Bello, G. (2017). Tracing the origin of the NS1 A188V substitution responsible for recent enhancement of Zika virus Asian genotype infectivity. Mem. Inst. Oswaldo Cruz 112, 793–795. doi: 10.1590/0074-02760170299
Guo, X., Yang, H., Wu, C., Jiang, J., Fan, J., Li, H., et al. (2015). Molecular characterization and viral origin of the first dengue outbreak in Xishuangbanna, Yunnan Province, China, 2013. Am. J. Trop. Med. Hyg. 93, 390–393. doi: 10.4269/ajtmh.14-0044
Guo, X., Yang, M., Jiang, J., Li, H., Zhu, C., Gui, Q., et al. (2016). [Molecular characteristics of dengue virus outbreak in China-Myanmar border region, Yunnan province, 2015]. Zhonghua Liu Xing Bing Xue Za Zhi 37, 398–401. doi: 10.3760/cma.j.issn.0254-6450.2016.03.022
Guo, X., Zhao, Q., Wu, C., Zuo, S., Zhang, X., Jia, N., et al. (2013). First isolation of dengue virus from Lao PDR in a Chinese traveler. Virol. J. 10:70. doi: 10.1186/1743-422X-10-70
Huang, X. Y., Ma, H. X., Wang, H. F., Du, Y. H., Su, J., Li, X. L., et al. (2014). Outbreak of dengue Fever in central China, 2013. Biomed. Environ. Sci. 27, 894–897. doi: 10.1186/s12916-015-0336-1
Imrie, A., Roche, C., Zhao, Z., Bennett, S., Laille, M., Effler, P., et al. (2010). Homology of complete genome sequences for dengue virus type-1, from dengue-fever- and dengue-haemorrhagic-fever-associated epidemics in Hawaii and French Polynesia. Ann. Trop. Med. Parasitol. 104, 225–235. doi: 10.1179/136485910X12647085215570
Kato, F., Kotaki, A., Yamaguchi, Y., Shiba, H., Hosono, K., Harada, S., et al. (2012). Identification and characterization of the short variable region of the Japanese encephalitis virus 3' NTR. Virus Genes 44, 191–197. doi: 10.1007/s11262-011-0685-6
Lima, F. R., Croda, M. G., Muniz, D. A., Gomes, I. T., Soares, K. R., Cardoso, M. R., et al. (2013). Evaluation of the traditional and revised World Health Organization classifications of dengue cases in Brazil. Clinics 68, 1299–1304. doi: 10.6061/clinics/2013(10)02
Liu, C., Liu, Q., Lin, H., Xin, B., and Nie, J. (2014). Spatial analysis of dengue fever in Guangdong Province, China, 2001-2006. Asia Pac. J. Public Health 26, 58–66. doi: 10.1177/1010539512472356
Liu, W., Zuo, L., and Zhou, Y. (2006). The distribution of DEN infected people in Dushan and Xingyi area of Yunnan-Guizhou Plateau, China. Cell. Mol. Immunol. 3, 473–476.
Liu, Y., Liu, J., Du, S., Shan, C., Nie, K., Zhang, R., et al. (2017). Evolutionary enhancement of Zika virus infectivity in Aedes aegypti mosquitoes. Nature 545, 482–486. doi: 10.1038/nature22365
Lustig, Y., Wolf, D., Halutz, O., and Schwartz, E. (2017). An outbreak of dengue virus (DENV) type 2 Cosmopolitan genotype in Israeli travellers returning from the Seychelles, April 2017. Euro Surveill. 22:30563. doi: 10.2807/1560-7917.ES.2017.22.26.30563
Naik, N. G., and Wu, H. N. (2015). Mutation of putative N-glycosylation sites on dengue virus NS4B decreases RNA replication. J. Virol. 89, 6746–6760. doi: 10.1128/JVI.00423-15
Narvaez, F., Gutierrez, G., Pérez, M. A., Elizondo, D., Nuñez, A., Balmaseda, A., et al. (2011). Evaluation of the traditional and revised WHO classifications of Dengue disease severity. PLoS Negl. Trop. Dis. 5:e1397. doi: 10.1371/journal.pntd.0001397
Nukui, Y., Tajima, S., Kotaki, A., Ito, M., Takasaki, T., Koike, K., et al. (2006). Novel dengue virus type 1 from travelers to Yap State, Micronesia. Emerging Infect. Dis. 12, 343–346. doi: 10.3201/eid1202.050733
Parameswaran, P., Wang, C., Trivedi, S. B., Eswarappa, M., Montoya, M., Balmaseda, A., et al. (2017). Intrahost selection pressures drive rapid dengue virus microevolution in acute human infections. Cell Host Microbe 22, 400–410.e405. doi: 10.1016/j.chom.2017.08.003
Ren, H., Ning, W., Lu, L., Zhuang, D., and Liu, Q. (2015). Characterization of dengue epidemics in mainland China over the past decade. J. Infect. Dev. Ctries. 9, 970–976. doi: 10.3855/jidc.5998
Santos-Sanz, S., Sierra-Moros, M. J., Oliva-Iñiguez, L., Sanchez-Gómez, A., Suarez-Rodriguez, B., Simón-Soria, F., et al. (2014). [Possible introduction and autochthonous transmission of dengue virus in Spain]. Rev. Esp. Salud Publica 88, 555–567. doi: 10.4321/S1135-57272014000500002
Sasmono, R. T., Wahid, I., Trimarsanto, H., Yohan, B., Wahyuni, S., Hertanto, M., et al. (2015). Genomic analysis and growth characteristic of dengue viruses from Makassar, Indonesia. Infect. Genet. Evol. 32, 165–177. doi: 10.1016/j.meegid.2015.03.006
Tajima, S., Nukui, Y., Ito, M., Takasaki, T., and Kurane, I. (2006). Nineteen nucleotides in the variable region of 3' non-translated region are dispensable for the replication of dengue type 1 virus in vitro. Virus Res. 116, 38–44. doi: 10.1016/j.virusres.2005.08.015
Tajima, S., Nukui, Y., Takasaki, T., and Kurane, I. (2007). Characterization of the variable region in the 3' non-translated region of dengue type 1 virus. J. Gen. Virol. 88, 2214–2222. doi: 10.1099/vir.0.82661-0
Vashist, S., Anantpadma, M., Sharma, H., and Vrati, S. (2009). La protein binds the predicted loop structures in the 3' non-coding region of Japanese encephalitis virus genome: role in virus replication. J. Gen. Virol. 90, 1343–1352. doi: 10.1099/vir.0.010850-0
Wang, X., Ma, D., Huang, X., Li, L., Li, D., Zhao, Y., et al. (2017). Complete genome analysis of dengue virus type 3 isolated from the 2013 dengue outbreak in Yunnan, China. Virus Res. 238, 164–170. doi: 10.1016/j.virusres.2017.06.015
Wei, Y., Jiang, T., Li, X., Zhao, H., Liu, Z., Deng, Y., et al. (2008). [Effects of RNA elements within 3'untranslated region on dengue virus translation]. Wei Sheng Wu Xue Bao 48, 583–588. doi: 10.3321/j.issn:0001-6209.2008.05.004
Xiong, Y., and Chen, Q. (2014). [Epidemiology of dengue fever in China since 1978]. Nan Fang Yi Ke Da Xue Xue Bao 34, 1822–1825. doi: 10.3969/j.issn.1673-4254.2014.12.24
Ye, Q., Li, X. F., Zhao, H., Li, S. H., Deng, Y. Q., Cao, R. Y., et al. (2012). A single nucleotide mutation in NS2A of Japanese encephalitis-live vaccine virus (SA14-14-2) ablates NS1' formation and contributes to attenuation. J. Gen. Virol. 93, 1959–1964. doi: 10.1099/vir.0.043844-0
Yuan, L., Huang, X. Y., Liu, Z. Y., Zhang, F., Zhu, X. L., Yu, J. Y., et al. (2017). A single mutation in the prM protein of Zika virus contributes to fetal microcephaly. Science 358, 933–936. doi: 10.1126/science.aam7120
Zhang, F. C., Zhao, H., Li, L. H., Jiang, T., Hong, W. X., Wang, J., et al. (2014). Severe dengue outbreak in Yunnan, China, 2013. Int. J. Infect. Dis. 27, 4–6. doi: 10.1016/j.ijid.2014.03.1392
Zhang, H. L., Fu, S. H., Deng, Z., Yuan, J., Jiang, H. Y., Li, M. H., et al. (2013). [An outbreak of imported dengue fever from Myanmar to the border of China, with its viral molecular epidemiological features]. Zhonghua Liu Xing Bing Xue Za Zhi 34, 428–432. doi: 10.3760/cma.j.issn.02546450.2013.05.004
Keywords: dengue outbreak, complete genome, phylogenetic trees, RNA secondary structure, protein secondary structure
Citation: Wen S, Ma D, Lin Y, Li L, Hong S, Li X, Wang X, Xi J, Qiu L, Pan Y, Chen J, Shan X and Sun Q (2018) Complete Genome Characterization of the 2017 Dengue Outbreak in Xishuangbanna, a Border City of China, Burma and Laos. Front. Cell. Infect. Microbiol. 8:148. doi: 10.3389/fcimb.2018.00148
Received: 23 December 2017; Accepted: 20 April 2018;
Published: 08 May 2018.
Edited by:
Jianming Qiu, University of Kansas Medical Center, United StatesReviewed by:
Xiao-Guang Chen, Southern Medical University, ChinaHaijian Zhou, Chinese Center for Disease Control and Prevention, China
Copyright © 2018 Wen, Ma, Lin, Li, Hong, Li, Wang, Xi, Qiu, Pan, Chen, Shan and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiyun Shan, sgtrhh@126.com
Qiangming Sun, qsun@imbcams.com.cn
†These authors have contributed equally to this work.