Background
==========

Dengue fever (DF) and two more severe syndromes, dengue shock syndrome (DSS) and dengue hemorrhagic fever (DHF) are important mosquito-borne diseases in tropical and subtropical regions \[[@B1],[@B2]\]. Since the first documented outbreak in 1779 in Jakarta, Indonesia, outbreaks have been documented in tropical and subtropical regions. It has been the maximum public health burden in South-East Asia countries \[[@B3]\]. Dengue epidemiology changes fast among epidemic countries, and keeps on expanding to the non-epidemic area \[[@B1]\]. Since the pathogens were first discovered by Japanese scientists in 1943, dengue viruses (DENV) were isolated from almost all South-East Asia countries including Thailand, Cambodia, Vietnam, Laos, Myanmar, Malaysia, Philippines, and Indonesia \[[@B3]\]. The dengue virus belongs to Flaviviridae family, and has four closely related but different serotypes (DENV-1 through -4) in nature that are circulating or co-circulating \[[@B1],[@B2]\]. Their genetic diversity remains wide, due to 1) absence of a proof-reading capacity in RNA-dependent RNA polymerases \[[@B4]\], 2) emergence of different lineages or clades during epidemic \[[@B5],[@B6]\], 3) increasing natural recombination \[[@B7]-[@B10]\], and 4) co-circulation of more than one serotypes in a locality \[[@B11]-[@B14]\].

Since imported DF epidemic had been reported in Hankou, Hangzhou, Shanghai and Guangzhou in 1920s and 1940s \[[@B15]\], there was no DF case reported in China till the outbreak occurred in 1978 in Foshan, Guangdong Province, with DENV-4 infection \[[@B16]\]. And then DF was prevalent in Guangdong, Guangxi and Hainan province. A DENV-1 epidemic occurred in 1979 in Zhongshan, Guangdong Province. This serotype of virus continued causing outbreaks over 2-3 years intervals. It become the dominating serotype and caused the latest outbreak in 2006 \[[@B17],[@B18]\]. DENV-3 epidemic was only recorded once in 1980 at Zhan County, Hainan Island, where in 1985 DENV-2 caused an epidemic during which the first DHF case was reported \[[@B19]\]. After this DENV-2 epidemic, DENV-2 continued to be transmitted into Guangdong, Guangxi and Hainan until 2001, including three outbreaks in Foshan in 1993, 1998 and in Jiangmen in 2001 \[[@B20]\]. The second DENV-4 outbreak occurred in Guangzhou city in 1990 \[[@B20]\]. It was estimated more than 700,540 hospitalized cases with 513 deaths from 1978 to 2007 \[[@B15]\]. Although the first isolate was sampled 30 years ago, dengue epidemic in southern China keeps increasing \[[@B15]\]. Lack of longitudinal research on dengue epidemics has hampered our understanding of dengue molecular origin and evolution in China

In this study, we determined the complete genome sequences of 14 dengue isolates collected in southern China during each epidemic year between 1978 and 2006 and E gene sequences from six patients of Guangdong, 2006. In combination with those published sequences in GenBank, we conducted an extensive molecular epidemiological analysis, aiming to determine where the DENV isolates in China originally came from, and what shaped their evolution.

Materials and methods
=====================

Ethics statement
----------------

This research was approved by the Review Board of Guangzhou Medical Research Institute, the Review Board of Shenzhen Second people\'s Hospital, the Review Board of Beijing Institute of Genomics, the Review Board of Beijing Genomics Institute in Shenzhen and the Review Board of Southwest University. Written informed consent was obtained from the patient for publication of this case report and accompanying images. A copy of the written consent is available for review by the Editor-in-Chief of this journal.

Sera
----

Six dengue fever patients\' sera were collected in Guangzhou city during the epidemic of DENV-1 in 2006 by Guangzhou Medicine Institute. These sera were collected within 7 days after onset of symptoms and stored at -20°C. All sera were tested positive for DEN-1 IgM by indirect immunofluorescence assay at the Guangzhou Medicine Institute.

Viruses
-------

Fourteen DENV strains were obtained from patients sera in DF epidemic in China during 1978 to 2006 (Table [1](#T1){ref-type="table"}). Sucking BALB/c mice and C6/36 cells were used to isolate viruses from clinical specimens. Serotypes of those isolates were confirmed by indirect immunofluorescence using anti-DENV monoclonal IgG (Guangzhou Medicine Institute, China). These isolates included eight DENV-1 isolates collected from Guangzhou and Chaozhou City of Guangdong province in 1991, 1995, 1997, 1999, 2003 and 2006 epidemic. Four DENV-2 isolates were collected from Fusan and Jiangmen City of Guangdong Province in 1993, 1998 and 2001 epidemic. Two DENV-4 isolates were collected in Foshan City during 1978 epidemic and in Guangzhou city during 1990 epidemic respectively.

###### 

Description of DENV isolates sequenced in this study

  Serotype   **Isolate**^***a***^   **Passage history (host, no. of passages)**^***b***^   Locality    Year
  ---------- ---------------------- ------------------------------------------------------ ----------- ------
             GD03/91                SMB, 3; C6/36, 1                                       Guangzhou   1991
             GD95/95                SMB, 2; C6/36, 1                                       Guangzhou   1995
             GD01/97                SMB, 1; C6/36, 1                                       Chaozhou    1997
  DENV-1     GD99/99                SMB, 1; C6/36, 1                                       Chaozhou    1999
             GD54/03                SMB, 1; C6/36, 1                                       Guangzhou   2003
             GD66/03                SMB, 1; C6/36, 1                                       Guangzhou   2003
             GD01/06                SMB, 1                                                 Guangzhou   2006
             GD02/06                SMB, 1                                                 Guangzhou   2006
                                                                                                       
             GD01/93                SMB, 2; C6/36, 1                                       Foshan      1993
  DENV-2     GD06/93                SMB, 2; C6/36, 1                                       Foshan      1993
             GD08/98                SMB, 1; C6/36, 1                                       Foshan      1998
             GD01/01                SMB, 1; C6/36, 1                                       Jiangmen    2001
                                                                                                       
  DENV-4     GD07/78                SMB, 2                                                 Foshan      1978
             GD09/90                SMB, 2; C6/36, 1                                       Guangzhou   1990

^*a*^All isolates were isolated from patients with primary infection manifesting classic DF symptoms at Guangdong Province, People\'s Republic of China.

^*b*^SMB, sucking mouse brain; C6/36, *Aedes albopictus*cell line.

Viral RNA extraction and RT-PCR methods
---------------------------------------

Viral RNA was extracted using TRIzol reagent (Invitrogen Corp., California, USA) according to the manufacturer\'s instruction. The cDNA was synthesized by reverse transcription from 10 μL of RNA at 42°C for 50 min in a 50 μL solution containing 50 mM Tris-HCl (pH 8.3), 75 mM KCl, 3 mM MgCl2, 10 mM DTT, 100 ng of the random hexamer primers, 200 U of Superscript II reverse transcriptase (Invitrogen), 25 U of RNasin (Invitrogen) and 0.5 mM dNTPs. The PCR was performed in a 25 μl mixture containing 2 μl of cDNA, 10 mM Tris/HCl (pH 8.4), 50 mM KCl, 2.5 mM MgCl2, 100 μM each dNTP, 1 U Taq DNA Polymerase (Shanghai Promage Corp., Shanghai, China), 0.25 μM each primer (Additional file [1](#S1){ref-type="supplementary-material"}, Table S1-S3). The amplification reactions consisted of an initial denature step of 3 min at 95°C, followed by 35 cycles of 30 s at 94°C, 30 s at 55°C, 60 s at 72°C.

Viral genome sequencing and E gene sequencing
---------------------------------------------

The complete genomes of the viruses were sequenced by a PCR-by-PCR strategy. At the same time, Envelope (E) gene from six patients\' sera was also amplified. All primers were synthesized by BGI (Beijing Genomics Institute, Beijing, China) according to the reference sequence data of DENV-1 Singapore strain S275/90 \[GenBank: [M87512](M87512)\], DENV-2 Jamaica strain 1409/88 \[GenBank: [M20558](M20558)\], and DENV-4 Dominica strain 814669/81 \[GenBank: [AF326573](AF326573)\]. Amplified products were detected by agarose gel electrophoresis and sequenced using ABI 3730xl DNA Analyzers (Applied Biosystems, Foster City, CA).

Assembly of Genome Sequences and Sequence Analysis
--------------------------------------------------

Genome assembly was performed independently by distinct operators using Phred-Phrap-Consed \[[@B21]\]. The consensus sequence was yielded exactly with BLASTN and custom script which was written by Perl for all strains. A single contig was obtained for each of the 14 isolates. The single contig was aligned to reference by BLASTN with the E-value of 1e-5. The computation of substitution tables were performed using a custom program written by Perl.

Neighbor-joining trees were constructed using MEGA version 3.1 with the Kimura-2 parameter corrections of multiple substitutions \[[@B22]\]. Reliability of nodes was assessed by bootstrap resampling with 1,000 replicates. One hundred and twenty-five E gene sequences and 60 complete genome sequences published in the GenBank were also involved.

RDP3 was used for recombination analysis. We employed Bootscan, Chimera, GENECONV, MaxChi, and RDP methods incorporated in RDP3 Beta 27 program \[[@B23]-[@B27]\]. General recombination settings for all methods were as follows: Sequences were considered as linear, and the highest acceptable P-value was set to 0.01, and event detected by two or more methods was taken into consideration.

Results
=======

Nucleotide sequence accession numbers
-------------------------------------

The determined DENV nucleotide sequences were deposited in GenBank database under the accession numbers [FJ196841](FJ196841)-[FJ196860](FJ196860).

Phylogenetic tree of DENV-1, -2, and -4 generated by E gene sequences
---------------------------------------------------------------------

To determine the DENV molecular epidemiology in southern China, 145 E gene sequences were used to infer a ML phylogenetic tree (Figure [1](#F1){ref-type="fig"}), including 1) 20 E gene sequences determined in the study, 2) 47 sequences with China localities available from GenBank, and 3) 78 representative sequences from GenBank referring to diverse genotype and geographic localities.

![**Neighbor-joining phylogenetic trees of 145 E gene sequences were conducted using Mega3.1 with the Kimura-2 parameter corrections of multiple substitutions**. Bootstrap values (in percentage \>65%) on each node were generated by using 1000 replications. The branch lengths are proportional to the number of nucleotide changes. All the Chinese isolates are in bold for clarity. A: phylogenetic tree of serotype DENV-1, DENV-2 and DENV-4 virus; B: phylogenetic tree of the DENV-1 isolates; C: phylogenetic tree of the DENV-2 isolates; D: phylogenetic tree of the DENV-4 isolates.](1743-422X-8-322-1){#F1}

The DENV-1 isolates from Guangdong belong to genotype I and IV (Figure [1B](#F1){ref-type="fig"}). Within genotype I and IV, several sub-genotypes can be identified. Two isolates (CN.GD01/97 and CN.GD99/99) determined in the study clustered together as a single clade with a 96% bootstrapping support, which was named as a new Asia-4 clade. DENV-1 strains isolated in 2001 and 2006 belong to Asia 3 clade in genotype I. However, the stains of 2006 could further be divided into three different sub-clades. Eleven strains of 2006 positioned at the end of the Asia 3 clade and the other four strains positioned at two different sub-clades. Even the eight strains (two isolates and six patients\' sera) collected in Guangzhou city during the epidemic of DENV-1 in 2006 also clustered into 2 subclades in Asia 3 clade. All eight strains collected during 2002 and 2003 period at Guangzhou were identical and clustered closely as genotype IV type.

Chinese DENV-2 isolates were characterized into four genotypes except for American genotype (Figure [1C](#F1){ref-type="fig"}). Broadly, those collected in pre-1990s fell into Asia/American and Asia II genotypes, while the post-1990s belonged to Asia I and Cosmopolitan genotypes. Two isolates collected in 1998 in Foshan, which belonged to Asia I genotype, were identical to the strain from Thailand in 1993. Three isolates of Jiangmen in 2001 epidemic showed closely identical to the strain from Philippines in 2000 epidemic.

Two DENV-4 epidemics occurred in China in the year of 1978 and 1990. Five isolates collected from the two epidemic years belonged to genotype I and genotype II (Figure [1D](#F1){ref-type="fig"}). Two isolates collected during 1990 epidemic belonged to genotype I and clustered closely with a strain isolated in 1984 epidemic of Philippines (PH. 12123/84) \[GenBank: [U18435](U18435)\]. The three other DENV-4 isolates from 1978 epidemic clustered together and fell into genotype II. They resided in phylogenetic tree much closer to the Indonesia isolate 30153 (ID. 30153/73) \[GenBank: [U18428](U18428)\] collected in 1973 when there was a severe dengue epidemic.

Phylogenetic tree of DENV-1, -2, and -4 generated by ORF sequences
------------------------------------------------------------------

Based on the 74 complete coding regions (14 determined in the study and 60 from GenBank), we constructed maximum likelihood (ML) phylogenetic tree (Figure [2](#F2){ref-type="fig"}).

![**Neighbor-joining phylogenetic trees of 74 complete coding regions were conducted using Mega3.1 with the Kimura-2 parameter corrections of multiple substitutions**. Bootstrap values (in percentage \>65%) on each node were generated by using 1000 replications. The branch lengths are proportional to the number of nucleotide changes. All the Chinese isolates are in bold for clarity.](1743-422X-8-322-2){#F2}

In the study, the whole genome and the single E gene phylogenetic tree are in good coincidence with only little differences. Two DENV-1 isolates, CN.GD01/97 and CN.GD99/99, were assumed to be new members of a defined Asia-4 clade in E gene tree (Figure [1B](#F1){ref-type="fig"}), while in the coding region phylogenetic tree; they clustered into the Asia-3 clade (Figure [2](#F2){ref-type="fig"}).

Recombination detection
-----------------------

RDP3 software was employed to examine the potential recombination events of all the Chinese isolates (43 DENV-1,20 DENV-2,2 DENV-3,5 DENV-4) on E gene. Three recombinant events were detected on the E gene of DENV-1 with low P-value (P \< 0.01) (Table [2](#T2){ref-type="table"}), which was confirmed by phylogenetic construction (Figure [3](#F3){ref-type="fig"}). Neither DENV-2 nor DENV-4 has any evidence of such recombination across their counterparts. The potential recombination events were also analyzed across the 32 complete ORF of all the Chinese isolates by RDP3 software. No inter-serotype recombination was detected, but 22 intra-serotype potential recombination events (17 for DENV-1, 5 for DENV-2) were confirmed by more than two methods (Table [3](#T3){ref-type="table"}).

###### 

Detected recombination events on the E gene region of Chinese DENV-1\*

  ------------------------------------------------------------------------------------------------------
  Event   Daughter                 Major Parent             Minor Parent             **Region**^**†**^
  ------- ------------------------ ------------------------ ------------------------ -------------------
  1       D1.GD23/95\              D1.GD14/97\              D1.GD95/95\              1-326
          ([AY373427](AY373427))   ([AY376737](AY376737))   ([FJ196844](FJ196844))   

                                                                                     

  2       D1.GD01/97\              D1.NHD98039/98\          D1.GZ02/03\              60-760
          ([FJ196847](FJ196847))   ([EF508198](EF508198))   ([DQ211349](DQ211349))   

          D1.GD99/99\                                                                71-760
          ([FJ196848](FJ196848))                                                     

                                                                                     

  3       D1.GZ01/04\              D1.STD06098\             D1.71/02\                6-1299
          ([EF032589](EF032589))   ([EF508207](EF508207))   ([EF025110](EF025110))   
  ------------------------------------------------------------------------------------------------------

\*Recombinant region and parental sequences were identified by five procedures (Bootscan, Chimarea, GENECONV, MaxChi, and RDP) integrated in RDP3 Package with a multiple comparison corrected P-value \< 0.01.

^†^Region coordinates are nucleotide positions of detected breakpoints in recombinant E gene sequence (1485bp).

![**Phylogenetic analysis based on non-recombination regions (Left) and recombination regions (Right) of event 1 (A), 2 (B), and 3 (C) constructed by Mega 3.1 using Neighbor-joining methods**. The event numbers are in accordance with table 2. Black circle, square, and triangle indicate the sequences of recombinant, major parent and minor parent, respectively. The trees are mid-rooted for clarity.](1743-422X-8-322-3){#F3}

###### 

Recombination events detected by RDP3^†^

  Recombinant Event Number   Breakpoint Positions in Recombinant Sequence   Recombinant Sequence(s)   Minor Parental Sequence   Major Parental Sequence   
  -------------------------- ---------------------------------------------- ------------------------- ------------------------- ------------------------- ------------
  1                          677                                            10104\*                   D2.GD06/93                D2.GD01/93                D2.GD08/98
  2                          5514                                           6548\*                    D1.GD05/99                D1.GD95/95                D1.GD01/06
  3                          1\*                                            7487                      D1.GD01/97                D1.GD14/97                D1.GD95/95
                                                                                                      D1.GD99/99                D1.ZJ01/04                D1.GZ01/95
  4                          900                                            1605                      D1.GD01/97                D1.GZ71/02                D1.GD14/97
                                                                                                      D1.GD99/99                D1.GZ01/95                D1.ZJ01/04
  5                          5763                                           6401                      D1.GD99/99                D1.GD66/03                D1.GD01/06
  6                          1\*                                            466                       D1.GD05/99                D1.GD23/95                D1.GD99/99
  7                          5585                                           6219                      D1.GD01/97                D1.GD66/03                D1.GD14/97
  8                          749                                            1167                      D1.GD23/95                D1.GD14/97                D1.GD95/95
  9                          8240                                           9343                      D2.GD08/98                D2.GD01/93                D2.GZ44/89
  10                         5457                                           6201\*                    D1.GZ71/02                D1.GZ01/95                D1.GD54/03
  11                         1606\*                                         2555                      D1.GD14/97                D1.GD99/99                D1.GD01/97
  12                         3711                                           4913                      D2.GD01/01                D2.GD06/93                D2.ZS01/01
  13                         2588\*                                         2898                      D1.GD23/95                D1.GD14/97                D1.GZ01/95
  14                         5171                                           5506                      D1.GD23/95                D1.GD14/97                D1.GD95/95
  15                         7890                                           8795                      D1.GD01/97                D1.GD99/99                Unknown
  16                         2708                                           3121                      D1.GD01/06                D1.GD54/03                D1.ZJ01/04
  17                         4889                                           5368                      D2.GD08/98                D2.ZS01/01                D2.GZ44/89
  18                         9645                                           10146\*                   D1.GD95/95                D1.GD03/91                D1.GZ01/95
  19                         7890\*                                         8385                      D1.GD99/99                D1.GZ71/02                Unknown
  20                         8820                                           9106                      D1.GD99/99                D1.GD23/95                D1.GD05/99
  21                         2335                                           2398\*                    D1.GD95/95                Unknown                   D1.GD23/95
  22                         1043\*                                         10104\*                   D2.GD06/93                D2.GD01/93                D2.ZS01/01

^†^= 22 final recombination events resulting from at least two methods confirmation integrated in RDP3 package with a multiple comparison corrected P-value \< 0.01 and one-by-one manual checking. Only the closest major and minor parental sequences were shown in the table. Recombinants of DENV-1 were bolded for clarity.

\* = The actual breakpoint position is undetermined (it was most likely overprinted by a subsequent recombination event).

Minor Parent = Parent contributing the smaller fraction of sequence.

Major Parent = Parent contributing the larger fraction of sequence.

Unknown = Only one parent and a recombinant need be in the alignment for a recombination event to be detectable.

Discussion
==========

Since the first well documented DF outbreak in Foshan of Guangdong province in 1978, dengue has been reported periodically in southern China, with the latest epidemic at Guangdong province in 2006. In this study, 14 Dengue isolates from 1978 to 2006 were sequenced and analyzed for molecular evolution.

Dengue virus exists as four antigenically distinct viruses designated as serotypes (DENV-1 through DENV-4), belonging to genus *Flavivirus*of family *Flaviviridae*. It has a positive-sense RNA genome that is translated as a single polyprotein and posttranslationally cleaved into three structural proteins and seven nonstructural proteins \[[@B28]\]. The envelope protein (E) is considered to be the immunodominant protein \[[@B29]\]. Dengue viruses also could be divided into different genotypes by the E gene \[[@B30],[@B31]\]. So there is correlation between serotypes and genotypes of Dengue virus. Phylogenetic trees of DENV-1, -2 and -4 generated by ORF and E gene sequences showed that there was a wide genetic diversity of DENVs isolated in southern China (Figure [1](#F1){ref-type="fig"}, Figure [2](#F2){ref-type="fig"}). DENV-1 strains exist in almost all of the clades of genotype I and IV but Asia 1 clade of genotype I (Figure [1B](#F1){ref-type="fig"}, Figure [2](#F2){ref-type="fig"}). DENV-2 strains are grouped into four genotypes except American genotype (Figure [1C](#F1){ref-type="fig"}, Figure [2](#F2){ref-type="fig"}). DENV-4 strains are grouped into 2 genotypes (I and II) (Figure [1D](#F1){ref-type="fig"}, Figure [2](#F2){ref-type="fig"}). It also showed that different serotypes and genotypes epidemic prevalence exist in the same city and even during the same epidemic (Figure [1](#F1){ref-type="fig"}, Figure [2](#F2){ref-type="fig"}). It indicated that these dengue viruses maybe have different origination.

Some of isolates, especially the first emerging of certain dengue serotype strains in southern China were closely related to those strains which were isolated in prior epidemics from neighboring Southeast Asia countries. It indicated that these dengue epidemics may be imported into China from the neighbor countries. For DENV-2, there were two outbreaks in Foshan in 1993 and 1998. Isolates from 1998 epidemic was closely related to the isolates from 1993 epidemic in Thailand (Figure [1C](#F1){ref-type="fig"}). The same situation was showed in DENV-4 isolates from 1978 and 1990 epidemics. Isolates from 1978 epidemic was closely related to the isolates from 1973 epidemic in Indonesia (Figure [1D](#F1){ref-type="fig"}). Isolates form 1990 epidemic was closely related to the isolates from 1984 epidemic in Philippines (Figure [1D](#F1){ref-type="fig"}). It indicated that these epidemics maybe imported into China from neighbor countries.

Although some data showed the dengue fever were imported from neighbor countries, those epidemics with the same serotype of DENV during continues years showed that endemic infection of dengue circulating locally may be also the important cause of Dengue epidemic in southern China. In the study, two isolates (CN.GD01/97 and CN.GD99/99) clustered together as a new Asia-4 clade (Figure [1B](#F1){ref-type="fig"}), suggesting that the isolates of 1999 maybe originated from the isolates of 1997. Sequence analysis showed that the DENV-1 strains isolated during the 2003 epidemic were closely related to the strains isolated during the 2002 epidemic of Guangzhou, suggesting that the 2002 strains had evolved in local and eventually cause the epidemic of 2003 (Figure [1B](#F1){ref-type="fig"}). Similar case was also showed in Guangzhou outbreak of DF in 2006. The DENV-1 strain isolated during the 2006 epidemic is closely related to the strains isolated during the 2001 epidemic (Figure [1B](#F1){ref-type="fig"}). It suggested that DENV probably circulated in local and caused the epidemics of 2001 and 2006. Moreover, there were several clear divisions within the Chinese 2006 isolates in the tip of the tree, even the eight strains collected in Guangzhou city by us during the epidemic of DENV-1 in 2006 also clustered into 2 subclades in Asia 3 clade, showing viral evolution or increase of genetic complexity even during a single epidemic (Figure [1B](#F1){ref-type="fig"}).

Recombination plays a role in dengue virus evolution and adaptation \[[@B7]\]. We analyzed the potential recombination events of all the Chinese isolates (43 DENV-1,20 DENV-2,2 DENV-3,5 DENV-4) on E gene. Three recombinant events were detected on the E gene of DENV-1 with low P-value (P \< 0.01) (Table [2](#T2){ref-type="table"}). The incongruence of phylogenetic trees constructed separately by recombinant regions and non-recombinant regions confirmed the recombinant events (Figure [3](#F3){ref-type="fig"}). The recombination event 1 showed recombination between present GD14/97 isolate \[GenBank: [AY376737](AY376737)\] as the major parent and the GD95/95 isolate \[GenBank: [FJ196846](FJ196846)\] as the minor parent, which led to a recombinant isolate of GD23/95 \[GenBank: [AY373427](AY373427)\](Figure [3A](#F3){ref-type="fig"}). The recombination event 2 showed recombination between present NHD98039/98 isolate \[GenBank: [EF508198](EF508198)\] as the major parent and the GZ02/03 isolate \[GenBank: [DQ211349](DQ211349)\] as the minor parent, which led to two recombinant isolates of GD01/97 \[GenBank: [FJ196847](FJ196847)\] and GD99/99 \[GenBank: [FJ196848](FJ196848)\](Figure [3B](#F3){ref-type="fig"}). The recombination event 3 showed recombination between present STD06098/06 isolate \[GenBank: [EF508206](EF508206)\] as the major parent and the 71/02 isolate \[GenBank: [EF025110](EF025110)\] as the minor parent, which led to a recombinant isolate of GZ01/04 \[GenBank: [EF032589](EF032589)\](Figure [3C](#F3){ref-type="fig"}). These studies of recombination provided an important context to understand the evolution of these viruses and sequence diversity generated in viruses in China.

We also analyzed the potential recombination across the 32 complete ORF of all the Chinese isolates, and found no inter-serotype recombination occurred. The results are coincided with previous data \[[@B8],[@B10],[@B32]\]. However, we found 22 intra-serotype potential recombination events (17 for DENV-1, 5 for DENV-2) (Table [3](#T3){ref-type="table"}). Interestingly, for DENV-1, 13 of 17 potential recombination events happened in the strains isolated in 1995, 1997 and 1999 epidemics (Table [3](#T3){ref-type="table"}), prior to which, in 1995, the severest dengue outbreak ever reported, indicating a possibility of viruses escaping human immunity by recombination following epidemic. However, whether the presumption makes sense or just a coincidence remains an open question.

Conclusions
===========

Our study suggested that Dengue fever epidemic in Southern China over the past 30 years presents two important modes, 1) imported-cases-induced endemic prevalence; 2) endogenous epidemic outbreak with natural epidemic focuses. Recombination analysis showed that recombination plays an important role in dengue virus evolution and adaptation.

Conflict of interests declaration
=================================

The authors declare that they have no competing interests.

Authors\' contributions
=======================

WLW, HQZ, ZT and LCL carried out the genes sequencing and phlogenetics analysis. ZJB and MYF carried out the virus\'s culture and serotype identification. BHT and JHL participated in the samples detection. JWL and WJC designed the study, perform the data analysis and wrote the manuscript. All authors read and approved the final manuscript.

Supplementary Material
======================

###### Additional file 1

**Table S1-S3**. Overlapping PCR primers for DENV-1, DENV-2 and DENV-4 genome sequences

###### 

Click here for file

Acknowledgements
================

We thank Dr. Lin Wu for critical reading of the manuscript. This work was sponsored by Infectious Diseases Special Project, Minister of Health of China (2008ZX10004-001, 2009ZX10004-103).
