Skip to main content


Full-length genome and molecular characterization of dengue virus serotype 2 isolated from an imported patient from Myanmar



Dengue is the most common mosquito-borne infection worldwide and a serious threat to global public health. Sporadic dengue virus serotype 2 (DENV-2) imported cases from Myanmar have been documented almost every year in Yunnan Province of China since 2005. However, the complete genome sequences of DENV-2 isolates imported from Myanmar are not available.


The full-length genome of the DENV-2 strain (YNPE2), isolated from an imported case from Myanmar in 2013, was identified by the next-generation sequencing. The extreme ends of the viral genome were validated by 5′/3′ RACE and Sanger sequencing. Furthermore, phylogenetic, recombination and selection pressure analyses were conducted for the molecular characterization of YNPE2 strain.


Whole-genome sequencing revealed that the full-length sequence of YNPE2 strain was 10,724 bases, with an open reading frame encoding for 3391 amino acids. The YNPE2 strain had 99.0% nucleotide identity and 99.8% amino acid identity with two closely related strains, ThD2_0078_01 strain (DQ181797) and DENV-2/TH/BID-V2157/200 strain (FJ639832). The phylogenetic analysis suggested that the YNPE2 strain belonged to Asian I genotype and was likely derived from Thailand strain (DQ181797). Moreover, selection pressure analysis revealed two amino acid sites of the NS4B and NS5 proteins, with important evidence of positive selection.


This study revealed the first complete genome sequence and molecular characterization of a DENV-2 strain (YNPE2) isolated from an imported case from Myanmar, thus providing a valuable reference genome source for future surveillance, epidemiology and vaccine development of DENV-2 virus in Yunnan, China.


Dengue, a re-emerging mosquito-borne human viral disease, is caused by the dengue virus (DENV). The global prevalence of the disease has markedly increased in recent decades, with about half of the world′s population under threat [1]. More than 390 million people worldwide are infected by DENV every year, of which 96 million have clinical symptoms, ranging from mild dengue fever (DF) to dengue hemorrhagic fever (DHF) or Dengue shock syndrome (DSS) [2].

DENV is a single stranded and positive polarity RNA virus with a genome of about 11,000 bases in length [3], and belongs to the genus Flavivirus, family Flaviviridae [4]. The viral genome only has a long open reading frame (ORF) that encodes a polyprotein, and is flanked by 5′ untranslated region (UTR) and 3′ UTR. The polyprotein can be processed into three structural proteins (capsid, membrane and envelope) and seven non-structural (NS) proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) [3]. The DENV is mainly transmitted by Aedes aegypti mosquitoes in tropical and sub-tropical countries [5]. Serologically, DENV is classified into four serotypes (DENV-1 to DENV-4), which share 65–70% homology [6]. Due to its high mutation rate, DENV has generated different genotypes and lineages within each serotype [7]. Revealing the genotype of the circulating DENV in epidemic areas is critical to disease surveillance, epidemiology and vaccine development.

With the increase in population mobility, the DENV and mosquitoes can quickly travel by modern transportation across the globe. Thus, the endemic areas of dengue have expanded from 9 countries before 1970 to over 100 countries worldwide [1]. Among them, South-East Asia, America and the Western Pacific regions are seriously affected by dengue, due to their tropical and sub-tropical climate. In South-East Asia, all countries have reported outbreaks of dengue since 1950 [8,9,10,11,12,13,14]. In China, the first dengue epidemic was documented in Guangdong Province in 1978 [15]. In recent decades, the number of dengue patients have gradually increased in some provinces of China, including Guangdong, Yunnan, Fujian and Zhejiang [16,17,18,19]. Generally, numerous dengue outbreaks in China are closely related to imported cases from countries in the South-East Asia. Thus, dengue remains an imported disease in China [20].

Yunnan Province is located in the Southwest border of China, and shares a 4060-km border with three dengue-endemic countries: Myanmar, Laos and Vietnam. With the construction of the China-ASEAN Free Trade Area and the transport infrastructure in the Greater Mekong Subregion, Yunnan and South Asia and South-East Asia countries are experiencing an increase in trade and tourism. Hence, Yunnan Province faces an increasing need for dengue prevention and control due to numerous imported cases from the neighboring countries. The national dengue surveillance program was implemented in Yunnan province in 2005, and an increasing number of imported cases have been reported each year [21, 22]. In 2008, an outbreak of dengue was caused by imported cases from Myanmar [21]. Moreover, a large-scale dengue outbreak was reported in 2013 in Dehong and Xishuangbanna Prefectures of Yunnan Province. The outbreak was predominantly related to imported cases from Myanmar and Laos [23]. Therefore, imported cases from Myanmar and Laos were likely the main cause of the dengue outbreak in Yunnan Province. However, majority of the previous studies on the imported DENV-2 strain were based on partial sequencing data of the envelope (E) gene. The complete genome sequence of any DENV-2 strain imported from Myanmar is not available. In this study, we reported full-length genome sequence and molecular characterization of the YNPE2 strain isolated from an imported case from Myanmar. This study will likely facilitate future surveillance, epidemiology studies and vaccine development of DENV-2 virus in Yunnan, China.


Case presentation

In September 2013, a 32-year-old female patient, who came from Myanmar, was clinically diagnosed with dengue fever in Ruli, Yunnan, China. She had dengue-like symptoms including fever, joint pain, myalgia and headache. The patient′s acute serum was positive for NS1 antigen by One Step Dengue nonstructural protein (NS) 1 RapiDipTM InstaTest (Cortez). The patient’s acute serum (2–5 ml) was collected with informed consent and stored at -80 °C until testing.

Virus isolation, identification and serotyping

Acute-phase serum samples collected from the dengue patient were used for virus isolation. Briefly, the serum was diluted at a ratio of 1:10 using RPMI 1640 medium (Biological Industries, USA), followed by incubation for 1 h with Aedes albopictus C6/36 cells for virus adsorption in a 75 cm2 flask (Corning, USA). Then, RPMI 1640 medium supplemented with 2% fetal bovine serum (FBS; Gibco) was added to maintain cell survival. Cells were incubated at 28 °C in 5% CO2 and examined daily for cytopathic effects (CPE). When CPE was observed, the cell culture supernatant was collected, which served as virus stock. The virus stock was used for infecting C6/36 cells. Once CPE was observed, the cell culture supernatant was collected and analyzed by reverse transcription PCR (RT-PCR) to determine the serotype of the virus. The viral RNA was extracted from infected cell culture supernatant using TRIzol reagent (Tiangen, Beijing, China) according to the manufacturer′s protocol. Total RNA was reverse transcribed using GoScript™ Reverse Transcriptase (Promega, WI, USA) according to the manufacturer′s instructions. The reaction was conducted for 5 min at 25 °C, followed by 1 h at 42 °C and 15 min at 70 °C in a BioRad C1000 Cycler system. Then, PCR was performed using Taq polymerase (TaKaRa, Dalian, China) using a reference protocol with DENV universal primers (D1: 5′ -TGAATATGCTGAAACGCGCGAGAAACCG-3′ and D2: 5′ -TTGCACCAACAGTCAATGTCTTCAGGTTC-3′) and DENV-2 typing primer (D1 and TS2: 5′ -CGCCAGAAGGGCCATGAACAG-3′) [24]. Agarose gel electrophoresis was conducted to check the specificity of amplicons [24].

Whole-genome sequencing analysis

The viral RNA was extracted from infected cell culture supernatant using the TRIzol reagent (Invitrogen, CA, USA) according to the manufacturer′s protocol. The total RNA was depleted of rRNA with an epicentre Ribo-ZeroTM kit (Illumina, CA, USA). The rRNA-depleted RNA was purified and fragmented using a RNA fragmentation kit (Ambion, USA). The sequencing library was constructed as previously described [25]. The sequencing was performed on an Illumina Hiseq 4000 platform (LC Sciences, Hangzhou, China) according to the manufacturer′s standard protocol, and 2 × 150 bp paired-end reads (PE150) were generated. The joint and low quality reads were removed from the raw sequencing reads to obtain the clean reads. The viral genome was assembled from the clean reads by the SPAdes (v3.10.1) program [26].

5′ / 3′ rapid amplification of cDNA ends (RACE)

The viral RNA was extracted from infected cell culture supernatant using the TRIzol reagent (Invitrogen, CA, USA). Since the total RNA was extracted from a non-eukaryotic organism and lacked a polyadenylated tail, the 5′-first-strand cDNA was synthesized with random primers and the 3′-first-strand cDNA was synthesized after adding a poly(A) tail using Poly(A) Polymerase (TaKaRa, Dalian, China) according to the instruction of SMARTer 5′/3′ RACE kit (Clontech, CA, USA). The 5′/3′ RACE PCR was performed by using SeqAmp DNA Polymerase (TaKaRa, Dalian, China) with Universal Primer Mix and 5′ Gene-Specific Primer (5′-CGCCATTATGGTGAAGCCTGGATGTCTC-3′) / 3′ Gene-Specific Primer (5′-ATCTGGGAGGCCACAAACCATGGAAGC-3′) according to the instruction of SMARTer 5′/ 3′ RACE kit. The PCR reaction involved 40 cycles of 94 °C for 30 s, 60 °C for 30 s, and 72 °C for 2 min in a BioRad C1000 Cycler system. The PCR products of 5′ RACE and 3′ RACE were diluted at a ratio of 1:10 with nuclease-free water. The second round PCRs were performed in a reaction mixture of 50 μL containing 5 μL of 10 × La Taq Buffer (Mg2+ Plus), 8 μL of dNTPs (2.5 mM), 5 μL of 10 × Universal Primer Mix, 1 μL of 5′/3′ Gene-Specific Primer, 2.5 μL of diluted 5′/3′ RACE PCR products and 1 unit of high fidelity La Taq DAN polymerase (TaKaRa, Dalian, China). The second round PCRs were conducted under conditions of one cycle of 94 °C for 5 min followed by 40 cycles of 94 °C for 30 s, 60 °C for 30 s, and 72 °C for 2 min. The amplified products were purified from 1.5% agarose gel using the Wizard SV Gel and PCR Clean-Up System (Promega, WI, USA). The purified products were ligated to the T4 vector (Promega, WI, USA), transformed into E. coli, and positive clones were selected for Sanger sequencing (BGI, Shenzhen).

Accession number

The full-length genome sequence of the DENV-2 strain YNPE2 was deposited in GenBank (accession number MF459663.3).

Genomic characterization analysis

The viral genomic characterization and amino acid sequence was analyzed using the DNASTAR Lasergene 11.1 software package (DNASTAR Inc., USA). Multiple sequence alignments were performed using the ClustalX2.1 program [27]. The Mfold software package ( was utilized to predict the secondary structure of the UTRs of the viral genome with the default folding parameters [28].

Phylogenetic analysis

Phylogenetic analysis based on the DENV-2 coding sequence and the envelope gene were performed using the MEGA program (version 6.06) [29, 30]. The phylogenetic trees were deduced by employing the Maximum Likelihood (ML) method based on the General Time Reversible model [31]. The bootstrap replications was set at 1000 for estimating the Branch topology of the phylogenetic tree.

Recombination analysis

Recombination detection of the DENV-2 was conducted by seven methods, including RDP [32], Chimaera [33], BootScan [34], 3Seq [35], GENCONV [36], MaxChi [37] and SiScan [38] available in the Recombination Detection Program (RDP version 4.95), with p < 0.01. The recombinant event was identified if a sequence was detected by at least three methods under the multiple comparison correction setting option.

Selection pressure analysis

The selection pressure analysis for codons of the open reading frame (ORF) of the viruses (excluding potential recombinants) was performed by single-likelihood ancestor counting (SLAC) [39], fixed effects likelihood (FEL) [40], internal fixed effect likelihoods (IFEL) [41] and mixed effects model of evolution (MEME) [42] available in HyPhy (Hypothesis testing using phylogenies) open-source software package ( The ratio of non-synonymous to synonymous substitutions (dN/dS ratio, ω) were calculated in four methods. The sites were considered significant if they were detected by three or more methods.


Identification and serotyping of virus

To observe the proliferation properties of the virus, C6/36 cells were repeatedly infected with virus stock. At 8 days post-infection, typical CPE was observed in the C6/36 cells. Furthermore, the RT-PCR analysis showed two specific bands, 511 bp and 119 bp in length, which were specifically amplified by DENV universal primers and DENV-2 typing primers (Fig. 1), indicating that the virus was dengue virus serotype 2. The virus was named as YNPE2 strain.

Fig. 1

The new DENV-2 isolate identified by RT-PCR. M: DNA marker; 1: a 511-bp nucleotide fragment was obtained after amplification with DENV universal primers; 2: a 119-bp nucleotide fragment was obtained after amplification with DENV-2 type-specific primers

Full-length sequence and genomic characterization of YNPE2

The results of RNA-seq yielded 7,811,508 clean reads from 8,328,038 raw reads after removing unqualified reads. The full-length sequence of the YNPE2 strain was assembled by SPAdes software based on the RNA-seq data and the referenced genome of DENV-2 (DQ181797). The assembled full-length sequence of YNPE2 was 10,741 nucleotides (nts) with an ORF encoding a polyprotein of 3391 amino acids, flanked by 160 nts and 405 nts at the 5′ UTR and 3′ UTR, respectively. The RACE was considered a reliable method for confirming the extreme ends of RNA genome. Thus, the UTR sequences of viral genome were verified by 5′/3′ RACE and Sanger sequencing. Interestingly, the Sanger sequencing results of 5′ UTR and 3′ UTR were different from the Illumina′s next-generation sequencing (NGS). Sanger sequencing showed that the 5′ UTR and 3′ UTR were 96 nts and 452 nts, respectively (Fig. 2). Thus, based on 5′/3′ RACE and Sanger sequencing, the complete genome of YNPE2 was corrected to 10, 724 nts, flanked by 96 nts and 452 nts at the 5′ UTR and 3′ UTR, respectively. To determine the sequence identity between the YNPE2 strain and other DENV-2 strains, the complete genome sequence of YNPE2 was aligned with the prototype strain of DENV-2 (New Guinea C, AF038403) and 38 other globally representative DENV-2 strains (Additional file 1: Table S1). As shown in Table 1, the YNPE2 strain showed 89.5–99.0% sequence identity with the other strains. No insertion or deletion of nucleotides were observed in the protein-coding region when aligned with these representative strains. Moreover, the alignment results revealed that base substitutions were distributed over the whole stretch of the alignment. A single base substitution at position 54 in the 5′ UTR was observed when aligned with the prototype. As compared to the 3′ UTR of the prototype, the YNPE2 had 14 substitutions (position 10,275 [G → A], 10,279 [A → G], 10,289 [G → A], 10,321 [T → C], 10,387 [T → C], 10,389 [C → T], 10,390 [A → T], 10,413 [G → A], 10,418 [A → T], 10,553 [G → A], 10,555 [T → C], 10,570 [G → A], 10,575 [A → G] and 10,609 [A → G]), a single base (G) deletion at position 10,415 and a single nucleotide (A) insertion at position 10,455. The RNA secondary structure prediction results showed that these nucleotides substitutions in the 3′ UTR of YNPE2 led to minor change in the secondary stem loop structure as compared to the prototype (Fig. 3). The predicted RNA secondary structure of 3′ UTR of YNPE2 had fewer loops than the prototype.

Fig. 2

The 5′/3′ RACE PCR and Sanger sequencing results of YNPE2. a The gel shows the 5′ RACE and 3′ RACE amplifications of YNPE2. M: DNA marker; 5′ RACE: a 849-bp nucleotide fragment was obtained after amplification with Universal Primer Mix and 5′ Gene-Specific Primer; 3′ RACE: a 275-bp nucleotide fragment was obtained after amplification with Universal Primer Mix and 3′ Gene-Specific Primer. b the Sanger sequencing result of 5′ RACE, only part of the sequencing result is displayed; c the Sanger sequencing result of 3′ RACE

Table 1 Nucleotide (NT) and amino acid (AA) identity of YNPE2 (MF459663) complete genome compared to representative globally diverse full-length DENV-2 sequences
Fig. 3

The predicted secondary structure of 3′ UTRs of the prototype and YNPE2. Left: the predicted secondary structure of 3′ UTR of the prototype. Right: the predicted secondary structure of 3′ UTR of YNPE2

Amino acid sequence diversity

In order to characterize viral amino acid sequences, the deduced amino acid sequence of the YNPE2 strain was aligned with 39 complete genome sequences of DENV-2 representing the global diversity of DENV-2 (Additional file 1: Table S1). The sequence alignment revealed 96.7–99.8% identity among the sequences (Table 1) and amino acid substitutions throughout the complete coding region. The result indicated 66 amino acid replacements between YNPE2 and the prototype strain of DENV-2 (New Guinea C, AF038403) (Table 2). In addition, based on nucleotide and amino acid homology, two DENV-2 strains, ThD2_0078_01 strain (DQ181797) and DENV-2/TH/BID-V2157/200 strain (FJ639832), were the most closely related to YNPE2.

Table 2 Description of amino acid substitutions in YNPE2 compared to the prototype (New Guinea C)

Phylogenetic analysis

To determine the genetic relationship between YNPE2 and other DENV-2 strains isolated from geographically diverse areas, phylogenetic analysis of the whole coding sequence (Additional file 2: Table S2) and the envelope gene (Additional file 3: Table S3) were performed. The phylogenetic tree based on the complete coding sequences revealed that the DENV-2 strain belonged to the Asian I genotype along with multiple DENV-2 isolates from Vietnam, Thailand, and Cambodia. The YNPE2 strain was clustered into one group, together with Cambodia strain (FJ639718) and Thailand strain (DQ181797) (Fig. 4). The phylogenetic analysis suggested that the YNPE2 strain was likely derived from the Cambodia strain (FJ639718) or the Thailand strain (DQ181797). Moreover, the phylogenetic tree based on 73 envelope gene sequences isolated from geographically diverse areas showed that YNPE2 was clustered into a group with three Myanmar isolates (KJ470758, KJ470752 and KJ470762) and three Chinese isolates (KX262954, KX262957 and KY038916) circulating in 2013, which suggested their common origin (Fig. 5). In this study, the phylogenetic trees of the whole coding sequence and envelope gene were similar, with little diversity.

Fig. 4

Phylogenetic tree based on the complete coding sequences of 49 DENV-2 strains. Each strain was abbreviated as the GenBank accession number followed by the country and year of isolation. The numbers at the nodes revealed bootstrap support for that node (bootstrap replications was 1000). The DENV-2 strain sequenced in this study is marked with a solid circle (). The complete sequence of DENV-1 (US/Hawaii/1944, EU848545) was used to root the tree

Fig. 5

Phylogenetic tree based on 73 envelope gene sequences of DENV-2 strains. Each strain was abbreviated as the GenBank accession number followed by the country and year of isolation. The numbers at the nodes revealed bootstrap support for that node (bootstrap replications was 1000). The DENV-2 sequenced in this study is marked with a solid circle (), and the strains isolated from Yunnan Province of China in 2013 are marked with solid diamonds (). The envelope gene sequence of DENV-1 (US/Hawaii/1944, AB609588) was used to root the tree

Recombination analysis

The potential recombination events of the YNPE2 strain and other representative DENV-2 isolates (Additional file 1: Table S1) were analyzed by RDP4 software. The results of preliminary analysis revealed that four recombination events might have occurred among these DENV-2 strains (p < 0.01) (Additional file 4: Table S4). These recombination events seemed to be limited to some strains, accounting for 10% (4 out of 40 strains). However, no potential recombination event was detected for the YNPE2 strain within the tested isolates.

Selection pressure analysis

To determine the influence of the selection pressure acting on complete coding region of DENV-2 isolates (Additional file 1: Table S1), an analysis based on the selection coefficient ω (dN/dS ratio) at every codon was conducted with a cutoff value (ω > 1) of positive selection. The results showed that most of the codons in DENV-2 were under negative selection. However, a total of 13 sites were verified to be under positive selection pressure through at least three out of the four approaches (Table 3). Moreover, two of these sites, positioned on NS4B and NS5 respectively, indicated important evidence of episodic positive selection through at least two different methods with a low p-value (p < 0.1). Limited evidence of positive selection was confirmed by SLAC, FEL and IFEL .

Table 3 Selection pressure analysis of the ORF (3391 codons) of DENV-2 using MEME, SLAC, FEL and IFEL


For human viral diseases, knowledge of the full-length genome sequence, genotype and genomic characterization of a pathogen are essential for disease surveillance, epidemiology, and vaccine development. In this study, we used complete genome sequencing and comparative sequence analysis to determine the complete genome sequence and the molecular characterization of a DENV-2 strain (YNPE2) isolated from an imported case from Myanmar. It provides critical information on the whole genome sequence and genomic characterization of DENV-2 isolates imported from Myanmar to China.

Geographically, Yunnan Province in China borders with three dengue hyper-endemic countries, including Myanmar, Laos and Vietnam. With the increase in border trade between China and Myanmar, the risk of dengue import is also increasing. In 2008, a dengue outbreak caused by imported case from Myanmar to Ruili City was monitored [21]. Forty-nine dengue patients were clinically diagnosed based on the IgM antibody and viral nucleic acid detection, of which 48 were imported from Myanmar [21]. Moreover, a recent study on dengue epidemic characterization in Yunnan also demonstrated that the imported cases from Myanmar and Laos were mainly responsible for the dengue outbreak in Yunnan [43]. Thus, the investigation of complete genome and molecular characterization of different imported DENV isolates is essential for control and prevention of DENV circulating in Yunnan. In the present study, next-generation sequencing (NGS) Illumina Hiseq 4000 was employed for whole genome sequencing of the YNPE2 strain based on its advantages of high throughput, ultra-high speed, no mark and low cost [44, 45]. The UTR sequences were confirmed by 5′/3′ RACE, which is a credible method for determining the precise UTRs of RNA genome. The 5′/3′ RACE results were inconsistent with NGS, perhaps due to some artefacts of deep sequencing or assembly protocols. Therefore, 5′/3′ RACE and Sanger sequencing are essential for determining the precise viral genome. A highly reliable full-length genome sequence of YNPE2 was obtained using NGS, 5′/3′ RACE and Sanger sequencing. The conserved UTRs at both ends of the viral genome can form stable secondary structures, which play important roles in viral translation, replication and assembly [46, 47]. Interestingly, genomic characterization analysis from this study revealed some base substitutions in the 3′ UTR. The potential biological implications of the nucleotide substitution remains unknown, but their molecular configuration was slightly changed. Further studies are needed to determine the influence of the substitution on DENV virulence and replication in host cells.

Numerous amino acid substitutions were observed in the structural and non-structural proteins of the YNPE2 strain as compared to the prototype of DENV-2 (Table 2). These amino acid substitutions may potentially influence viral fitness and virulence in host cells. Interestingly, some amino acid substitutions, such as E71D, E126K and F402I observed in the envelope protein in this study, have been confirmed in previous studies [48, 49]. These amino acid substitutions enhanced the neurovirulence of the virus. Moreover, multiple amino acid substitutions in NS1, NS2A, NS2B and NS5 proteins of YNPE2 were observed. The potential influence of these substitutions needs further investigation.

Gene recombination plays a key role in the evolution and adaptation of DENV [50]. For molecular characterization of the YNPE2 strain, we selected 39 other representative globally diverse DENV-2 complete sequences to identify potential recombination events. No potential recombination event was found in the YNPE2 strain. Perhaps, more DENV-2 complete sequences are needed to uncover any potential recombination events of YNPE2. The analysis of site-specific selection pressures for DENV-2 indicated that strong negative selection was present, which was ascribed to the evolutionary constraints of arboviruses [7]. The positive selection from host-specific immunity contributes to viral protein mutations that allow the virus to escape host immune recognition [51, 52]. In this study, two sites with important evidence of positive selection were determined. One site in NS4B was reported in an earlier study of sylvatic DENV-2 and had an inhibitory effect on the interferon signaling pathway [53, 54]. Another site in NS5 was reported for the first time in our study, and its functional importance needs further investigation.

The phylogenetic analysis, based on the complete coding sequence, showed a close relationship between the YNPE2 strain and DQ181797 (Thailand, 2001), which suggested that the YNPE2 strain very likely originated in Thailand. This conclusion was further supported by the results of the sequence and amino acid identity. A phylogenetic tree based on the envelope gene revealed different subgenotype circulating in Myanmar. Moreover, the YNPE2 strain and seven Yunnan strains are closely related to isolates from Myanmar, suggesting that the dengue outbreaks in Yunnan Province are mainly caused by imported cases from Myanmar.


This study presented the first complete sequence of DENV-2 imported from Myanmar, and its molecular characterization. These results will stimulate further investigations of DENV infection, proliferation, pathogenicity, epidemiology and vaccine development.



Dengue virus serotype 2


Dengue fever


Dengue hemorrhagic fever


Dengue shock syndrome


Fetal bovine serum


Fixed effects likelihood


Internal fixed effect likelihoods


Mixed effects model of evolution


Maximum likelihood


Next-generation sequencing






Open reading frame


Polymerase chain reaction


Recombination Detection Program


Single-likelihood ancestor counting


Untranslated regions


  1. 1.

    WHO: World Health Organization. Dengue and Severe Dengue. Fact Sheet. Updated April 2017. 2017 Accessed 4 July 2017.

  2. 2.

    Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, Drake JM, Brownstein JS, Hoen AG, Sankoh O, et al. The global distribution and burden of dengue. Nature. 2013;496:504–7.

  3. 3.

    Chambers TJ, Hahn CS, Galler R, Rice CM. Flavivirus genome organization, expression, and replication. Annu Rev Microbiol. 1990;44:649–88.

  4. 4.

    Henchal EA, Putnak JR. The dengue viruses. Clin Microbiol Rev. 1990;3:376–96.

  5. 5.

    Salazar MI, Richardson JH, Sánchez-Vargas I, Olson KE, Beaty BJ. Dengue virus type 2: replication and tropisms in orally infected Aedes aegypti mosquitoes. BMC Microbiol. 2007;7:9.

  6. 6.

    Anoop M, Mathew AJ, Jayakumar B, Issac A, Nair S, Abraham R, Anupriya MG, Sreekumar E. Complete genome sequencing and evolutionary analysis of dengue virus serotype 1 isolates from an outbreak in Kerala, South India. Virus Genes. 2012;45:1–13.

  7. 7.

    Holmes E, Twiddy S. The origin, emergence and evolutionary genetics of dengue virus. Infect Genet Evol. 2003;3:19–28.

  8. 8.

    Khampapongpane B, Lewis HC, Ketmayoon P, Phonekeo D, Somoulay V, Khamsing A, Phengxay M, Sisouk T, Vongphrachanha P, Bryantcd JE. National dengue surveillance in the Lao People’s Democratic Republic,2006–2012: epidemiological and laboratory findings. Surveillance Report. 2014;5:7–13.

  9. 9.

    Cuong HQ, Vu NT, Cazelles B, Boni MF, Thai KTD, Rabaa MA, Quang LC, Simmons CP, Huu TN, Anders KL. Spatiotemporal dynamics of dengue epidemics, Southern Vietnam. Emerg Infect Dis. 2013;19:945–53.

  10. 10.

    Chen SP. Molecular evolution and epidemiology of four serotypes of dengue virus in Thailand from 1973 to 2007. Epidemiol Infect. 2013;141:419–24.

  11. 11.

    Kyaw AK, Ngwe Tun MM, Moi ML, Nabeshima T, Soe KT, Thwe SM, Myint AA, Maung KTT, Aung W, Hayasaka D, et al. Clinical, virological and epidemiological characterization of dengue outbreak in Myanmar, 2015. Epidemiol Infect. 2017;145:1886–97.

  12. 12.

    Massad E, Coutinho FA, Ma S, Burattini MN. A hypothesis for the 2007 dengue outbreak in Singapore. Epidemiol Infect. 2010;138:951–7.

  13. 13.

    Mohd-Zaki AH, Brett J, Ismail E, L’Azou M. Epidemiology of dengue disease in Malaysia (2000-2012): a systematic literature review. PLoS Negl Trop Dis. 2014;8:e3159.

  14. 14.

    Vong S, Khieu V, Glass O, Ly S, Duong V, Huy R, Ngan C, Wichmann O, Letson GW, Margolis HS, Buchy P. Dengue incidence in urban and rural Cambodia: results from population-based active fever surveillance, 2006–2008. PLoS Negl Trop Dis. 2010;4:e903.

  15. 15.

    Fan W, Yu S, Cosgriff TM. The reemergence of dengue in China. Rev Infect Dis. 1989;11:S847–53.

  16. 16.

    Guo R, Lin J, Li L, Ke C, He J, Zhong H, Zhou H, Peng Z, Yang F, Liang W. The prevalence and endemic nature of dengue infections in Guangdong, South China: an epidemiological, serological, and etiological study from 2005–2011. PLoS One. 2014;9:e85596.

  17. 17.

    Wang B, Li Y, Feng Y, Zhou H, Liang Y, Dai J, Qin W, Hu Y, Wang Y, Zhang L, et al. Phylogenetic analysis of dengue virus reveals the high relatedness between imported and local strains during the 2013 dengue outbreak in Yunnan, China: a retrospective analysis. BMC Infect Dis. 2015;15:142.

  18. 18.

    Wang J, Chen H, Huang M, Zhang Y, Xie J, Yan Y, Zheng K, Weng Y. Epidemiological and etiological investigation of dengue fever in the Fujian province of China during 2004–2014. Sci China Life Sci. 2017;60:72–80.

  19. 19.

    Sun J, Lin J, Yan J, Fan W, Lu L, Lv H, Hou J, Ling F, Fu T, Chen Z, et al. Dengue virus serotype 3 subtype III, Zhejiang Province, China. Emerg Infect Dis. 2011;17:321–3.

  20. 20.

    Sang S, Chen B, Wu H, Yang Z, Di B, Wang L, Tao X, Liu X, Liu Q. Dengue is still an imported disease in China: a case study in Guangzhou. Infect Genet Evol. 2015;32:178–90.

  21. 21.

    Zhang H, Fu S, Deng Z, Yuan J, Jiang H, Li M, Gao X, Wang J, Liu Y. 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. 2013;34:428–32.

  22. 22.

    Jinyong J, Xiaofang G, Yerong T, Mingdon Y, Chun W, Hongning Z. Surveillance and control of imported dengue cases in Yunnan from 2004 to 2014. Chin J Vec Biol Control. 2016;27:5–8.

  23. 23.

    Wang B, Yang H, Feng Y, Zhou H, Dai J, Hu Y, Zhang L, Wang Y, Baloch Z, Xia X. The distinct distribution and phylogenetic characteristics of dengue virus serotypes/genotypes during the 2013 outbreak in Yunnan, China phylogenetic characteristics of 2013 dengue outbreak in Yunnan, China. Infect Genet Evol. 2016;37:1–7.

  24. 24.

    Ministry of Health PsRoC. WS 216-2008 diagnostic criteria and principle of management of dengue fever. Beijing: China Standard Publishing House; 2008.

  25. 25.

    Riesgo A, Farrar N, Windsor PJ, Giribet G, Leys SP. The analysis of eight transcriptomes from all poriferan classes reveals surprising genetic complexity in sponges. Mol Biol Evol. 2014;31:1102–20.

  26. 26.

    Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.

  27. 27.

    Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.

  28. 28.

    Zuker M. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003;31:3406–15.

  29. 29.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

  30. 30.

    Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.

  31. 31.

    NMaK S. Molecular evolution and Phylogenetics: Oxford University Press; 2000.

  32. 32.

    Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1:vev003.

  33. 33.

    Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci U S A. 2001;98:13757–62.

  34. 34.

    Martin DP, Posada D, Crandall KA, Williamson C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retrovir. 2005;21:98–102.

  35. 35.

    Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007;176:1035–47.

  36. 36.

    Padidam M, Sawyer S, Fauquet CM. Possible emergence of new Geminiviruses by frequent recombination. Virology. 1999;265:218–25.

  37. 37.

    Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34:126–9.

  38. 38.

    Gibbs MJ, Gibbs AJ. Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000;16:573–82.

  39. 39.

    Pond SL, Frost SD. Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005;21:2531–3.

  40. 40.

    Kosakovsky Pond SL, Frost SD. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005;22:1208–22.

  41. 41.

    Pond SL, Frost SD, Grossman Z, Gravenor MB, Richman DD, Brown AJ. Adaptation to different human populations by HIV-1 revealed by codon-based analyses. PLoS Comput Biol. 2006;2:e62.

  42. 42.

    Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Kosakovsky Pond SL. Detecting individual sites subject to episodic diversifying selection. PLoS Genet. 2012;8:e1002764.

  43. 43.

    Hu T, Zhang H, Feng Y, Fan J, Tang T, Liu Y, Zhang L, Yin X, Chen G, Li H, et al. Epidemiological and molecular characteristics of emergent dengue virus in Yunnan Province near the China-MyanmarLaos border, 2013–2015. BMC Infect Dis. 2017;17:331.

  44. 44.

    Liu L, Li Y, Li S, Hu N, He Y, Pong R, Lin D, Lu L, Law M. Comparison of next-generation sequencing systems. J Biomed Biotechnol. 2012;2012:251364.

  45. 45.

    van Nimwegen KJM, van Soest RA, Veltman JA, Nelen MR, van der Wilt GJ, Vissers LELM, Grutters JPC. Is the $1000 genome as near as we think? A cost analysis of next-generation sequencing. Clin Chem. 2016;62:1458–64.

  46. 46.

    Chiu WW, Kinney RM, Dreher TW. Control of translation by the 5′- and 3′-terminal regions of the dengue virus genome. J Virol. 2005;79:8303–15.

  47. 47.

    Leardkamolkarn V, Sirigulpanit W, Kinney RM. Characterization of recombinant dengue-2 virus derived from a single nucleotide substitution in the 5′ noncoding region. J Biomed Biotechnol. 2010;2010:934694.

  48. 48.

    Bray M, Men R, Tokimatsu I, Lai C. Genetic determinants responsible for acquisition of dengue type 2 virus mouse neurovirulence. J Virol. 1998;72:1647–51.

  49. 49.

    Gualano R, Pryor M, Cauchi M, Wright P, Davidson A. Identification of a major determinant of mouse neurovirulence of dengue virus type 2 using stably cloned genomic-length cDNA. J Gen Virol. 1998;79:437–46.

  50. 50.

    Worobey M, Rambaut A, Holmes EC. Widespread intra-serotype recombination in natural populations of dengue virus. Proc Natl Acad Sci U S A. 1999;96:7352–7.

  51. 51.

    Suzuki Y. Natural selection on the influenza virus genome. Mol Biol Evol. 2006;23:1902–11.

  52. 52.

    Seoighe C, Ketwaroo F, Pillay V, Scheffler K, Wood N, Duffet R, Zvelebil M, Martinson N, McIntyre J, Morris L, Hide W. A model of directional selection applied to the evolution of drug resistance in HIV-1. Mol Biol Evol. 2007;24:1025–31.

  53. 53.

    Vasilakis N, Holmes EC, Fokam EB, Faye O, Diallo M, Sall AA, Weaver SC. Evolutionary processes among sylvatic dengue type 2 viruses. J Virol. 2007;81:9591–5.

  54. 54.

    Dash PK, Sharma S, Soni M, Agarwal A, Parida M, Rao PVL. Complete genome sequencing and evolutionary analysis of Indian isolates of dengue virus type 2. Biochem Biophys Res Commun. 2013;436:478–85.

Download references


We would like to thank Dr. Nigel W. Fraser of Department of Microbiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, USA, for proof-reading and reviewing our manuscript.


This work was supported by grants from the CAMS Innovation Fund for Medical Sciences (CIFMS) (grant no. 2017-I2M-3-022), the National Key Research and Development Program of the Ministry of Science and Technology of China (2016YFC1202300), the State Project for Essential Drug Researchand Development (SQ2018ZX091510), the National Natural Science Foundation of China (31500724), the Yunnan Applied Basic Research Projects (2017FB115, 2017ZF007), and the Peking Union Medical College Youth Fund (3332016113), and the Fundamental Research Funds for the Central Universities, China.

Availability of data and materials

The genome sequence of the DENV-2 (YNPE2 strain) has been deposited into GenBank database under the accession number MF459663.3.

Author information

JDS, and YZH designed the study. ZPZ, JDS, XFG, LM, NZH, JS and MNW were involved in cell culture, viral infection, and molecular biology experiments. ZPZ conducted the data analyses and wrote the paper with help from all the authors. All the authors read and approved the final form of the manuscript.

Correspondence to Hongning Zhou or Yunzhang Hu.

Ethics declarations

Ethics approval and consent to participate

Written consent was obtained from the patient. The study protocol was approved by the Institutional Ethics Committee of Institute of Medical Biology, Chinese Academy of Medical Sciences, and Peking Union Medical College.

Consent for publication

All the authors consent to publish.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Details of DENV-2 sequences of was used in multiple sequence alignments, recombination analysis and selection pressure analysis. (DOCX 21 kb)

Additional file 2:

Table S2. Details of the full-length sequences of DENV-2 was used in phylogenetic analysis. (DOCX 21 kb)

Additional file 3:

Table S3. Details of the envelope gene sequences of DENV-2 was used in phylogenetic analysis. (DOCX 23 kb)

Additional file 4:

Table S4. DENV-2 potential recombination events obtained by RDP4 package. (DOCX 18 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark


  • Dengue virus
  • YNPE2 strain
  • Complete genome
  • Molecular characterization
  • 5′/3′ RACE