- Research article
- Open Access
Dengue virus type 1 clade replacement in recurring homotypic outbreaks
BMC Evolutionary Biologyvolume 13, Article number: 213 (2013)
Recurring dengue outbreaks occur in cyclical pattern in most endemic countries. The recurrences of dengue virus (DENV) infection predispose the population to increased risk of contracting the severe forms of dengue. Understanding the DENV evolutionary mechanism underlying the recurring dengue outbreaks has important implications for epidemic prediction and disease control.
We used a set of viral envelope (E) gene to reconstruct the phylogeny of DENV-1 isolated between the periods of 1987–2011 in Malaysia. Phylogenetic analysis of DENV-1 E gene revealed that genotype I virus clade replacements were associated with the cyclical pattern of major DENV-1 outbreaks in Malaysia. A total of 9 non-conservative amino acid substitutions in the DENV-1 E gene consensus were identified; 4 in domain I, 3 in domain II and 2 in domain III. Selection pressure analyses did not reveal any positively selected codon site within the full length E gene sequences (1485 nt, 495 codons). A total of 183 (mean dN/dS = 0.0413) negatively selected sites were found within the Malaysian isolates; neither positive nor negative selection was noted for the remaining 312 codons. All the viruses were cross-neutralized by the respective patient sera suggesting no strong support for immunological advantage of any of the amino acid substitutions.
DENV-1 clade replacement is associated with recurrences of major DENV-1 outbreaks in Malaysia. Our findings are consistent with those of other studies that the DENV-1 clade replacement is a stochastic event independent of positive selection.
Dengue virus (DENV) is a virus of the Flaviviridae family. It is an enveloped positive-sense single-stranded RNA arbovirus with a genome of approximately 11 kb . There are four antigenically distinct DENV serotypes; DENV-1, DENV-2, DENV-3 and DENV-4  and each serotype shows phylogenetically distinct genotypes . The virus is transmitted to susceptible hosts through bites of infected mosquitoes. The virus is maintained in sylvatic nonhuman primate/sylvatic mosquitoes and endemic human/urban/peridomestic mosquitoes cycles. All four DENV serotypes are believed to have independently evolved from separate sylvatic ancestral lineages through either peridomestic/urban mosquitoes or human hosts 100–1,500 years ago . Currently, an estimated 3.6 billion persons living in dengue-endemic countries are at risk of contracting dengue. The number of countries reporting dengue has in recent years escalated to more than 125 suggesting successful adaptation and dissemination of the virus .
In dengue endemic regions, heterotypic and homotypic major dengue outbreaks occur in cyclical patterns of approximately every 3–5 years and 7–10 years, respectively [6–10]. The major concern associated with recurring dengue outbreak in endemic countries is the risk of contracting the severe forms of dengue especially following second infection with a heterotypic virus . Antibody-dependent enhancement [12, 13], original antigenic sin [14, 15], cytokine storm , and autoimmune responses [17, 18] are the possible mechanisms contributing to the manifestation of severe dengue. It has been reported that infection with one DENV serotype confers lifelong protection against homotypic reinfection but only temporary cross-protection against heterotypic infection . The presence of sub-neutralizing and cross-reacting antibodies is suggested to play important role in the manifestation of the severe dengue [20, 21]. Reports of repeated infection with dengue is however, not uncommon in dengue endemic regions [22, 23]. Understanding the factors contributing to the recurrence of dengue outbreaks has important implications for our understanding of dengue epidemiology. Knowledge gained from this understanding could help improve dengue surveillance and outbreak prediction and preparation. It could also help to facilitate selection of better dengue vaccine candidates.
Earlier studies have suggested that DENV clade replacement is linked to the recurring and cyclical pattern of dengue outbreaks in many endemic countries [6, 24–27]. From these studies, it is suggested that clade replacement is associated with positive selection due to the differences in viral fitness between clades; new virus with a higher viremia level in human  or enhanced infectivity to mosquito [28–31] could be positively selected to replace the old virus which was less fit. In contrary, several other phylogenetic studies of DENV have suggested that the clade replacement is solely a stochastic event due to the virus population bottleneck effects [24, 25]. Although the potential mechanisms of DENV evolution underlying the clade replacement have been investigated, correlation studies involving the host immunological factors have not been adequately addressed. Zhang et al. proposed that the DENV-1 clade replacement is associated with the cross-protective immunity accorded by DENV-4 based on a longitudinal dengue epidemiological study in Bangkok . Whereas Adams et al. employed a mathematical model to demonstrate that the degree of interserotypic cross-protective immunity could account for the cyclical pattern of heterotypic outbreak in Bangkok .
In the present study, we used the recurring DENV-1 outbreaks in Malaysia which occurred in 1987, 1997 and 2004 [8, 33], as a study model. The availability of serially collected DENV-1 since 1987 within a single locality (Klang Valley) provides us with an opportunity to explore the temporal phylogenetic evolution that shapes the virus clade replacement in recurring DENV-1 outbreaks in Malaysia. We used DENV-1 full length envelope (E) gene sequence to reconstruct the phylogeny and investigated the presence of homotypic cross-neutralizing antibody of patients with primary DENV-1 infection.
Results and discussion
In the study, a total of 335 Malaysian DENV-1 isolates collected from 1987 to 2011 were used. All viruses were isolated from dengue patients living in single locality, Klang Valley. Between the periods, three major DENV-1 outbreaks occurred in 1987, 1997 and 2004 with a cyclical pattern of ~8 years interval (Figure 1) [8, 33]. Preliminary phylogenetic analysis performed using the full E gene sequence revealed that the virus collection comprised of DENV-1 from genotype I, II, III and sylvatic group with frequency of 80.6% (n = 270), 16.1% (n = 54), 3.0% (n = 10) and 0.3% (n = 1), respectively (data not shown), suggesting that DENV-1 outbreaks in Malaysia over the last 3 decades were mainly caused by the genotype I and II viruses. Figure 2 summarized the relative percentages of DENV-1 genotypes isolated during the three major outbreaks. Due to the large number of DENV-1 available for each outbreak, the viral sequences with >98% nucleotide identity within each monophyletic group were reduced to 3–6 sequences. A total of 44 representative isolates that covered all possible viral genetic diversity over time were presented in this study.
The Bayesian phylogenetic tree constructed using the full E gene sequence (1485 nt) revealed the existence of at least six distinct genotypes (Figure 3): three ancestral genotypes (Japan/Hawaii 1943–1945, Thailand 1954–1964 and Malaysia/sylvatic 1972–2005) and three endemic/epidemic genotypes (genotype I, II and III) in concordance to the previous classification method [33, 34]. DENV-1 genotype I viruses were involved in outbreaks in 1987, 1997 and 2004, while genotype II viruses were only involved in the 1997 outbreak.
Genotype I comprised mainly of the Asian strains with the exception of an African strain (D1.Djibouti/98). The genotype I viruses from the 1987, 1997 and 2004 outbreaks grouped into clade Ia, Ib and Ic, respectively in a ladder-like topology (Figure 3) suggesting temporal clade replacement occurred in each of the major DENV-1 outbreak. All the genotype I clades consisted of Thai viruses with approximately similar isolation years (±7 years), suggesting that the same pool of DENV-1 strains circulated within Thailand . Viruses in clade Ia and Ib shared clear ancestral lineages with Thai isolates whereas those forming clade Ic could have evolved from another Malaysian isolate seen in 1996 (D1.Malaysia.04834/96). The clade Ic viruses caused the largest major DENV-1 outbreak in Malaysia in 2004 and the viruses shared high sequence similarity to isolates recovered from a major outbreak in Singapore in 2005 . The clade Ic viruses were recovered much earlier in Thailand (1997), China (1998), Cambodia (2001), Vietnam (2002) and Myanmar (2002) before it caused outbreak in Malaysia and Singapore. Though all the clade Ic viruses shared common ancestral lineage to the Malaysian 1996 isolate, it is possible that the Malaysian 1996 virus was endemic in Malaysia but remained in the background until it achieved fitness to cause outbreak later in 2004 . This fit the expected homotypic cycle of major dengue outbreak in Malaysia at approximately every 8 years. Since the 2004 outbreak, clade Ic viruses were continuously isolated in Malaysia, Singapore, Thailand, Vietnam, Cambodia and China until 2011 and no clade replacement event of DENV-1 was observed then. The magnitude of 2004 outbreak involving many Asian countries could probably extent the time interval between the occurrences of major homotypic DENV-1 outbreaks in the future. To date, the clade Ic viruses have been imported into Japan (2001), South Korea (2005), Taiwan (2008), Sri Lanka (2009), Indonesia (2010) and Ireland (2010) [36–39] suggesting a possible global spread of the virus to other parts of the world. This could represent the most successful distribution of DENV-1 to date akin to that observed for the Cosmopolitan DENV-2 . On the other hand, the persistence of a ladder-like phylogenetic tree topology predicts possible emergence of a clade Id DENV-1 in the next major DENV-1 outbreak in Malaysia perhaps in or around 2019.
The genotype II viruses obtained in the study shared high sequence similarities to viruses from a wider geographical distribution including Asia and Pacific Ocean regions with the exception of an Indian Ocean strain (D1.Seychelles.NIID41/03). In addition to clade Ib, genotype II viruses (clade IIa and IIb) were also recovered during the 1997 outbreak. Phylogenetic analysis suggests that the clade IIa viruses could have evolved from another Malaysian isolate seen in 1991 (D1.Malaysia.94946/91). The clade IIa viruses are likely the indigenous DENV-1 of Malaysia as no virus from other country was found within the clade. In contrast, clade IIb viruses grouped with isolates from Indonesia and South Pacific islands including Micronesia and Polynesia. All the clade IIb viruses shared common ancestral lineage to the Indonesia 1988 isolate. The presence of the Indonesian isolates during the 1997 outbreak could be due to importation of the viruses through influx of migrant workers as there was a surge of workers from Indonesia to Malaysia during that time . In 2004, two isolates, D1.Malaysia.32858/04 and D1.Malaysia.32581/04, which shared high sequence similarities to the clade IIb viruses were recovered. This could be either another incidence of importation of the Indonesian isolate or an in situ evolution of the clade IIb viruses. The importations of DENV-1 strains from the neighboring countries would increase the genetic diversity of DENV-1 in Malaysia.
Genotype III consisted of mainly DENV-1 isolates from Latin America. Asian and African strains of the genotype III, however, have been occasionally isolated. The genotype III viruses from Malaysia, D1.Malaysia.59266/87 and D1.Malaysia.11738/97, grouped with isolates from Saudi Arabia, Brunei, India, Singapore and China and shared common ancestral lineage to the Thai 1980 isolate. In our samples, only 10 genotype III isolates (6 in 1987, 2 in 1995, 1 in 1996 and 1 in 1997) were recovered from patients during the period 1987–1997. The virus has not been isolated since 1998 which suggests its possible extinction from Malaysia. Similar findings were observed in Thailand and Myanmar in 1983 and 1998, respectively. The lineage extinction of genotype III was hypothesized to be due to a stochastic event attributable to possible low rate of virus transmission during an inter-epidemic period [6, 42]. However, the emergence of genotype III virus in Singapore in 2007 suggests possible reintroduction of genotype III virus to Malaysia in the future .
Using the ancestral sylvatic virus E amino acid sequence (D1.Malaysia.Sylvatic.P72_1244/72) as the basal sequence, the consensus amino acid sequences of all the Malaysian isolates were found highly conserved, with identities ranging from 96 to 98% between the genotypes (Figure 4). Although the genotype Japan/Hawaii 1943–1945 viruses served as the common ancestral lineage for all the genotype I viruses (Figure 3), the identities of consensus amino acid sequences between sylvatic virus and genotype I viruses were found to be higher (97%) than those between the Japan/Hawaii 1943–1945 viruses and genotype I viruses (ranging from 95 to 97%) (data not shown). Three conserved yet functionally important regions of DENV-1 E were observed: the twelve disulfide bond-forming cysteine residues; the glycine-rich fusion domain (residues 98–111); and two glycosylated asparagines at position Asn-67 and Asn-153 (data not shown) [44–46]. Out of 27 amino acid substitutions observed, only 9 (highlighted in grey) resulted in polarity changes (non-conservative amino acid substitutions). These non-conservative substitutions were found mostly in domain I (residues 37, 52, 157 and 161), followed by domain II (residues 88, 120 and 272) and domain III (residues 305 and 369), suggesting that domain I may be less functionally critical than domain II and domain III . The individual functional effect of each conservative amino acid substitution may be unnoticeable or relatively smaller than those of non-conservative substitutions, but the cumulative effects of two or more neighboring conservative substitutions at the three-dimensional space could be important for the proper protein function.
In this study, we noted that the non-conservative amino acid substitutions E120K (domain II) and G157E/D (domain I) differentiated the three endemic genotypes from the ancestral sylvatic group. Several other amino acid substitutions gave rise to additional three endemic genotypes (Figure 4). Within genotype I, two non-conservative amino acid substitutions P305S (domain III) and N390S (domain III) were identified in the 1997 viruses when compared against the 1987 viruses. No non-conservative amino acid substitution was observed between the viruses from 1997 and 2004 outbreaks. The decreasing number of amino acid substitution among the genotype I viruses from 1987 to 2004 could suggest that the viruses may have reached a point where all of the non-conservative amino acid substitutions were deleterious.
Within genotype II, two non-conservative amino acid substitutions D157E and T161I in domain I further differentiated the viruses into the clade IIa and IIb. Temporal analysis of amino acid substitution was not performed for the clade IIa and IIb viruses due to their similar emergence time during the 1997 outbreak. Whereas within genotype III, two non-conservative amino acid substitutions L305S (domain III) and M272T (domain II) were found between clade IIIa and IIIb. The M272T substitution located at the “kl” β-hairpin region (residues 270 to 279). This “kl” β-hairpin is a pH-dependent hinge region involved in the fusion-activating conformational change . Therefore, any non-synonymous mutation in the “kl” β-hairpin region could affect virus replication by altering the fusion pH threshold resulting in inefficient virus replication [48, 49]. This could help explain the rapid extinction of the genotype III virus as the virus could not be efficiently propagated.
Selection pressure analyses of viral E gene collectively found 183 (mean dN/dS = 0.0413), 136 (mean dN/dS = 0.0852), 245 (mean dN/dS = 0.0559) and 196 (mean dN/dS = 0.0704) negatively selected sites within the Malaysian isolates, the genotype I, II and III isolates, respectively. No positively selected codon site was noted within the Malaysian isolates and within each of the genotypes. All observed amino acid substitutions in Malaysian DENV-1 E gene consensuses (Figure 4) were found to be under negative selection except at the codon 8, 52, 161, 272, 339 and 390 which had no evidence of either negative or positive selection but probably happened through neutral genetic drift. Based on the Bayesian inference, the evolutionary rates of DENV-1 E gene were estimated to be 3.39 × 10-4-14.64 × 10-4 substitutions/site/year. The evolutionary rates of DENV-1 E gene did not differ substantially from those of other DENV serotypes (5.42 × 10-4-11.58 × 10-4 substitutions/site/year) . These rates, however, were relatively lower than those of non-vector-borne RNA viruses (~10-3 substitutions/site/year) such as influenza viruses and HIV . The DENV-1 evolution in Malaysia was probably constrained by purifying selection .
The cross-neutralizing capacity of convalescent patient serum samples from 1997 (Ib, IIa and IIb) and 2004 (Ic) outbreaks were evaluated against viruses from seven clades: Ia (D1.Malaysia.59280/87), Ib (D1.Malaysia.11177/97), Ic (D1.Malaysia.36000/05), IIa (D1.Malaysia.11708/97), IIb (D1.Malaysia.12796/97), IIIa (D1.Malaysia.59266/87) and sylvatic group (D1.Malaysia.36046/05) (Table 1). The neutralization titers of Ib sera against genotype I viruses (FRNT80 = 320–1280) were higher than those of Ic, IIa and IIb sera (FRNT80 = 80–320). All sera had FRNT80 = 80 against genotype II and IIIa viruses except Ic sera, which showed FRNT80 = 320 against the genotype IIb and IIIa viruses. The neutralization titers of the IIa sera against sylvatic virus (FRNT80 = 80) were lower than those of Ib, Ic and IIb sera (FRNT80 = 320). Overall, all viruses were neutralized by pooled sera from the 1997 and 2004 outbreaks with neutralization titers ranging from FRNT80 = 80–1280. Interserotypic interaction of DENV-1 patient’s sera with other DENV serotypes, however, is not available in this study to rule out any possible cross-reactivity among them. Results presented here suggested that there was no apparent immunological advantage accorded by the new amino acid substitutions for the virus in countering the homotypic human herd immunity. Therefore, the availability of a new susceptible human population may be obligatory for DENV-1 to initiate a new outbreak.
The human herd immunity developed after the DENV-1 outbreak could suppress the horizontal transmission of virus to human population by infected mosquito . This could restrict the virus to continue to exist only in mosquito population through vertical transmission [53, 54]. The population dynamics of mosquito over time, however, are mainly fluctuated by the stochastic environmental factors and the intermittent vector control measures [55, 56]. This could repeatedly cause virus population bottlenecks (small population size) which favor the amino acid substitutions by genetic drift over those by natural selection . Even if the positively selected substitutions are advantageous, it could be randomly lost over time due to the bottleneck effect. By contrast, the neutral substitutions could be easily fixed in small virus population over time by chance . In our study, six codons with observed amino acid substitutions in DENV-1 E gene consensuses were under neutral genetic drift.
Our results suggest that the DENV-1 evolution in recurring outbreak in Malaysia is likely to be a stochastic phenomenon, possibly driven by both negative selection and genetic drift [6, 24, 25, 52]. Whether the observed amino acid substitutions would cause enhanced virus replication in mosquito or human requires further investigation. During the inter-epidemic period, the increased viral fitness, however, may have no potential epidemic impact due to the human herd immunity that broadly cross-neutralized among different DENV-1 genotypes. Phylogenetic analysis in this study revealed that the emergences of clade IIb (1997 outbreak) and clade Ic viruses (2004 outbreak) were detected as early as in 1991 (D1.Malaysia.94946/91) and 1996 (D1.Malaysia.04834/96), respectively. The reason why these viruses did not cause outbreak when they were first encountered could be due to the presence in the population of cross-neutralizing immunity accorded by other DENV-1 genotypes. The cyclical pattern of recurring DENV-1 outbreaks in Malaysia suggests that a lag period of 7–10 years after the old DENV-1 outbreak would be required for a newly evolved virus to initiate another homotypic outbreak. During this period, an increasing number of new susceptible human hosts may gradually break down the human herd immunity below the protective threshold, thus allowing the restoration of virus horizontal transmission cycle in human population. This suggests the possibility that it is not always the fittest, but the fortuitous DENV-1 would be selected to be the clade replacement candidate which would initiate the next outbreak. And this could occur when the ecological situation favors the virus transmission. The appearances of rare isolates such as D1.Malaysia.92844/91, D1.Malaysia.10837/97, D1.Malaysia.32581/04, D1.Malaysia.32858/04 and D1.Malaysia.296236/11 (Figure 3, marked with *) could represent the potential clade replacement candidates for the future major outbreaks.
The present study is limited by working only with the viral E gene, such that evolutionary pressures on other gene regions will be missed. Also, during the inter-epidemic periods between the homotypic outbreaks, the available sample set selected for sequencing is limited. This leaves gaps in understanding the extinction and emergence of the different clades. Thirdly, the interserotypic interaction between the DENV-1 with other co-circulating virus serotypes is not addressed in this study. This is important as the interserotypic immune reaction could account for the alternating epidemic pattern of DENV serotypes in endemic regions. The cross-neutralization assay in this study is within the limitation where the possible cross-reactivity between different DENV serotypes was not examined. This was due to the unavailability of sufficient convalescent sera from dengue patients since the study was using retrospectively collected sera. As such further studies which will include early convalescent serum from DENV-2, DENV-3 or DENV-4 infection is desirable. Another limitation of the present study is that the immune responses of patients were averaged in the pooled sera; the neutralizing capacity of each serum might vary depending on the time of blood collection and genetic background of patients. Notwithstanding the limitation of the study, findings from the study provide opportunities to understand the possible mechanisms driving to cyclical pattern of major dengue outbreaks in endemic regions.
DENV-1 clade replacement is associated with the recurring major DENV-1 outbreaks in Malaysia. Our findings from selection pressure analyses and neutralization assays are consistent with earlier studies [24, 25] suggesting that virus clade replacement is stochastic rather than driven by positive selection.
The DENV-1 isolates used in this study were isolated from Klang Valley during the period of time spanning 25 years (1987–2011). Almost all the isolates came from patients who seeked medical treatment at the University Malaya Medical Centre (UMMC), a major referral hospital serving an estimated of 8.1 million people, mostly living within 25 km of the hospital. All the viruses were archived at the UMMC Diagnostic Virology Laboratory. The isolates at passage 1 were used to inoculate C6/36 (Aedes albopictus) mosquito cells for one week and virus RNA was extracted from this infected cell culture supernatant.
RNA extraction, RT-PCR and nucleotide sequencing
Viral RNA was extracted using QIAamp Viral RNA Mini Kit (Qiagen, Germany) strictly following the manufacturer’s instructions. One-step RT-PCR was performed to amplify DENV-1 E gene using the primers as previously described . The amplified DNA fragments were purified using QIAquick Gel Extraction Kit (Qiagen, Germany) and sequenced using the BigDye Terminator v3.1 Cycle Sequencing Kit on an automated capillary DNA sequencer 3730xl DNA Analyzer (Applied Biosystems, USA).
The DENV-1 E gene dataset for phylogenetic analysis comprised of 44 isolates from Malaysia and 121 global sequences representing the 6 distinct DENV-1 genotypes as previously reported [33, 34]. The maximum clade credibility tree (MCC) was inferred by using the Bayesian Markov Chain Monte Carlo (MCMC) method implemented in BEAST version 1.6.2 . The best-fit model of nucleotide substitution was selected by Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) as implemented in jModelTest 0.1.1 . Both AIC and BIC approaches found the TN93 + G4 + I model as the best-fit model for the dataset. Two independent MCMC analyses (100 million steps) were performed utilizing relaxed uncorrelated lognormal molecular clock with Bayesian skyline as coalescent prior. The convergence of the chain was evaluated by using Tracer 1.5 . The effective sample size (ESS) values of >200 indicated sufficient level of sampling. The MCC tree was generated and visualized by using TreeAnnotator program and FigTree 1.2.3, respectively. The degree of uncertainty in each parameter estimate is provided by the 95% highest posterior density (HPD) values.
Selection pressure analysis
The selection pressure analysis was performed by using the online facility of Datamonkey web server (http://www.datamonkey.org) based on the neighbor joining tree . Four datasets were included, one containing only Malaysian isolates (n = 44) and the other three containing isolates belonged to the three different virus groups, genotype I (n = 63), genotype II (n = 36) and genotype III (n = 57). Five likelihood-based algorithms: SLAC, FEL, IFEL, REL and PARRIS methods were used to identify the existence of positive selection pressure at individual codon sites in the DENV-1 E gene. Sites were considered to be under positive selection if the ratio of non-synonymous (dN) to synonymous (dS) substitutions per site (ratio dN/dS) indicated with high statistical significance (P <0.1 / Bayes factor >50).
The study obtained ethics approval from the UMMC Medical Ethics Committee (Ethics Committee/IRB Reference Number: 860.24 and 908.9). Paired (acute and convalescent) DENV-1 patient serum samples from 1997 (Ib, IIa and IIb) and 2004 (Ic) outbreaks were obtained from UMMC Diagnostic Virology Laboratory. Informed consents were not obtained from the patients as this study was using retrospectively collected sera. The acute and convalescent serum samples were taken at a mean of 4.3 days (range, 1–6 days) and 12.5 days (range, 5–24 days) post-onset, respectively. Acute DENV-1 infection was confirmed by virus isolation and a seroconversion of anti-DENV IgM in paired serum samples. Serum samples with optical density (OD) absorbance value ≥2.0 in the in-house IgM capture enzyme-linked immunosorbent assay (ELISA) were detected as IgM-positive sera . Serum samples were classified as either primary or secondary infection sera using the hemagglutination inhibition (HI) test . Convalescent serum samples with an HI titer <1:1280 were classified as primary infection sera, whilst those with an HI titer ≥1:2560 were classified as secondary infection sera .
Foci reduction neutralization test (FRNT)
Serum of convalescing dengue patient with primary infection with the respective DENV-1 genotype I and genotype II were used to determine their cross-neutralizing capacities against viruses from other subgroups. Due to the very limited amount of serum, serum samples (2–3 serum samples) from each of the different genotype group were pooled. The pooled serum was heat-inactivated at 56°C for 30 min and diluted four-fold (1:10 to 1:10240) in FBS-free EMEM medium. Serum samples (100 μl) were incubated with 100 FFU of virus (100 μl) at 37°C for 1 h. Subsequently, the mixtures (200 μl) were transferred to the C6/36 cell monolayer in 24-well plate and incubated at room temperature for 1 h. The mixtures were replaced with overlay EMEM medium [supplemented with 2% FBS and 1.5% carboxymethylcellulose (CMC)] and the plate was incubated at 28°C for 4 days. Foci of infected cells were visualized by using peroxidase-based foci staining assay as described previously [62, 65]. The neutralizing antibody titer was expressed as the maximum serum dilution yielding 80% reduction in foci formed (FRNT80).
Henchal EA, Putnak JR: The dengue viruses. Clin Microbiol Rev. 1990, 3 (4): 376-396.
Russell PK, Nisalak A: Dengue virus identification by the plaque reduction neutralization test. J Immunol. 1967, 99 (2): 291-296.
Holmes EC, Burch SS: The causes and consequences of genetic variation in dengue virus. Trends Microbiol. 2000, 8 (2): 74-77. 10.1016/S0966-842X(99)01669-8.
Wang E, Ni H, Xu R, Barrett AD, Watowich SJ, Gubler DJ, Weaver SC: Evolutionary relationships of endemic/epidemic and sylvatic dengue viruses. J Virol. 2000, 74 (7): 3227-3234. 10.1128/JVI.74.7.3227-3234.2000.
Gubler DJ: The economic burden of dengue. Am J Trop Med Hyg. 2012, 86 (5): 743-744. 10.4269/ajtmh.2012.12-0157.
Zhang C, Mammen MP, Chinnawirotpisan P, Klungthong C, Rodpradit P, Monkongdee P, Nimmannitya S, Kalayanarooj S, Holmes EC: Clade replacements in dengue virus serotypes 1 and 3 are associated with changing serotype prevalence. J Virol. 2005, 79 (24): 15123-15130. 10.1128/JVI.79.24.15123-15130.2005.
Nisalak A, Endy TP, Nimmannitya S, Kalayanarooj S, Thisayakorn U, Scott RM, Burke DS, Hoke CH, Innis BL, Vaughn DW: Serotype-specific dengue virus circulation and dengue disease in Bangkok, Thailand from 1973 to 1999. Am J Trop Med Hyg. 2003, 68 (2): 191-202.
Abubakar S, Shafee N: Outlook of dengue in Malaysia: a century later. Malays J Pathol. 2002, 24 (1): 23-27.
Vu TT, Holmes EC, Duong V, Nguyen TQ, Tran TH, Quail M, Churcher C, Parkhill J, Cardosa J, Farrar J, Wills B, Lennon NJ, Birren BW, Buchy P, Henn MR, Simmons CP: Emergence of the Asian 1 genotype of dengue virus serotype 2 in viet nam: in vivo fitness advantage and lineage replacement in South-East Asia. PLoS Negl Trop Dis. 2010, 4 (7): e757-10.1371/journal.pntd.0000757.
Bennett SN, Drummond AJ, Kapan DD, Suchard MA, Munoz-Jordan JL, Pybus OG, Holmes EC, Gubler DJ: Epidemic dynamics revealed in dengue evolution. Mol Biol Evol. 2010, 27 (4): 811-818. 10.1093/molbev/msp285.
Sam SS, Syed Omar SF, Teoh BT, Abd-Jamil J, AbuBakar S: Review of dengue hemorrhagic fever fatal cases seen among adults: a retrospective study. PLoS Negl Trop Dis. 2013, 7 (5): e2194-10.1371/journal.pntd.0002194.
Halstead SB: Observations related to pathogensis of dengue hemorrhagic fever. VI. Hypotheses and discussion. Yale J Biol Med. 1970, 42 (5): 350-362.
Kouri GP, Guzman MG, Bravo JR, Triana C: Dengue haemorrhagic fever/dengue shock syndrome: lessons from the Cuban epidemic, 1981. Bull World Health Organ. 1989, 67 (4): 375-380.
Halstead SB, Rojanasuphot S, Sangkawibha N: Original antigenic sin in dengue. Am J Trop Med Hyg. 1983, 32 (1): 154-156.
Mongkolsapaya J, Dejnirattisai W, Xu XN, Vasanawathana S, Tangthawornchaikul N, Chairunsri A, Sawasdivorn S, Duangchinda T, Dong T, Rowland-Jones S, Yenchitsomanus PT, McMichael A, Malasit P, Screaton G: Original antigenic sin and apoptosis in the pathogenesis of dengue hemorrhagic fever. Nat Med. 2003, 9 (7): 921-927. 10.1038/nm887.
Pang T, Cardosa MJ, Guzman MG: Of cascades and perfect storms: the immunopathogenesis of dengue haemorrhagic fever-dengue shock syndrome (DHF/DSS). Immunol Cell Biol. 2007, 85 (1): 43-45. 10.1038/sj.icb.7100008.
Markoff LJ, Innis BL, Houghten R, Henchal LS: Development of cross-reactive antibodies to plasminogen during the immune response to dengue virus infection. J Infect Dis. 1991, 164 (2): 294-301. 10.1093/infdis/164.2.294.
Falconar AK: The dengue virus nonstructural-1 protein (NS1) generates antibodies to common epitopes on human blood clotting, integrin/adhesin proteins and binds to human endothelial cells: potential implications in haemorrhagic fever pathogenesis. Arch Virol. 1997, 142 (5): 897-916. 10.1007/s007050050127.
Kurane I, Ennis FE: Immunity and immunopathology in dengue virus infections. Semin Immunol. 1992, 4 (2): 121-127.
Teerasut C, Lei HY, Natthakorn K, Jittmittraphap A, Thammapalo S, Witayathawornwong P, Lin YT, Jarman RG, Sawanpanyalert P, Jampangern W, Limkittikul K: The serotype-independent but concentration-dependent enhancing antibodies among Thai dengue patients. Southeast Asian J Trop Med Public Health. 2012, 43 (3): 624-633.
Dejnirattisai W, Jumnainsong A, Onsirisakul N, Fitton P, Vasanawathana S, Limpitikul W, Puttikhunt C, Edwards C, Duangchinda T, Supasa S, Chawansuntati K, Malasit P, Mongkolsapaya J, Screaton G: Cross-reacting antibodies enhance dengue virus infection in humans. Science. 2010, 328 (5979): 745-748. 10.1126/science.1185181.
Guzman MG, Kouri G, Valdes L, Bravo J, Alvarez M, Vazques S, Delgado I, Halstead SB: Epidemiologic studies on Dengue in Santiago de Cuba, 1997. Am J Epidemiol. 2000, 152 (9): 793-799. 10.1093/aje/152.9.793. discussion 804
Sangkawibha N, Rojanasuphot S, Ahandrik S, Viriyapongse S, Jatanasen S, Salitul V, Phanthumachinda B, Halstead SB: Risk factors in dengue shock syndrome: a prospective epidemiologic study in Rayong, Thailand. I. The 1980 outbreak. Am J Epidemiol. 1984, 120 (5): 653-669.
Myat Thu H, Lowry K, Jiang L, Hlaing T, Holmes EC, Aaskov J: Lineage extinction and replacement in dengue type 1 virus populations are due to stochastic events rather than to natural selection. Virology. 2005, 336 (2): 163-172. 10.1016/j.virol.2005.03.018.
Wittke V, Robb TE, Thu HM, Nisalak A, Nimmannitya S, Kalayanrooj S, Vaughn DW, Endy TP, Holmes EC, Aaskov JG: Extinction and rapid emergence of strains of dengue 3 virus during an interepidemic period. Virology. 2002, 301 (1): 148-156. 10.1006/viro.2002.1549.
Bennett SN, Holmes EC, Chirivella M, Rodriguez DM, Beltran M, Vorndam V, Gubler DJ, McMillan WO: Selection-driven evolution of emergent dengue virus. Mol Biol Evol. 2003, 20 (10): 1650-1658. 10.1093/molbev/msg182.
Carrillo-Valenzo E, Danis-Lozano R, Velasco-Hernandez JX, Sanchez-Burgos G, Alpuche C, Lopez I, Rosales C, Baronti C, de Lamballerie X, Holmes EC, Ramos-Castaneda J: Evolution of dengue virus in Mexico is characterized by frequent lineage replacement. Arch Virol. 2010, 155 (9): 1401-1412. 10.1007/s00705-010-0721-1.
Anderson JR, Rico-Hesse R: Aedes aegypti vectorial capacity is determined by the infecting genotype of dengue virus. Am J Trop Med Hyg. 2006, 75 (5): 886-892.
Armstrong PM, Rico-Hesse R: Efficiency of dengue serotype 2 virus strains to infect and disseminate in Aedes aegypti. Am J Trop Med Hyg. 2003, 68 (5): 539-544.
Hanley KA, Nelson JT, Schirtzinger EE, Whitehead SS, Hanson CT: Superior infectivity for mosquito vectors contributes to competitive displacement among strains of dengue virus. BMC Ecol. 2008, 8: 1-10.1186/1472-6785-8-1.
Lambrechts L, Fansiri T, Pongsiri A, Thaisomboonsuk B, Klungthong C, Richardson JH, Ponlawat A, Jarman RG, Scott TW: Dengue-1 virus clade replacement in Thailand associated with enhanced mosquito transmission. J Virol. 2012, 86 (3): 1853-1861. 10.1128/JVI.06458-11.
Adams B, Holmes EC, Zhang C, Mammen MP, Nimmannitya S, Kalayanarooj S, Boots M: Cross-protective immunity can account for the alternating epidemic pattern of dengue virus serotypes circulating in Bangkok. Proc Natl Acad Sci U S A. 2006, 103 (38): 14234-14239. 10.1073/pnas.0602768103.
Teoh BT, Sam SS, Abd-Jamil J, AbuBakar S: Isolation of ancestral sylvatic dengue virus type 1, Malaysia. Emerg Infect Dis. 2010, 16 (11): 1783-1785. 10.3201/eid1611.100721.
A-Nuegoonpipat A, Berlioz-Arthaud A, Chow V, Endy T, Lowry K, le Mai Q, Ninh TU, Pyke A, Reid M, Reynes JM, Su Yun ST, Thu HM, Wong SS, Holmes EC, Aaskov J: Sustained transmission of dengue virus type 1 in the Pacific due to repeated introductions of different Asian strains. Virology. 2004, 329 (2): 505-512. 10.1016/j.virol.2004.08.029.
Schreiber MJ, Holmes EC, Ong SH, Soh HS, Liu W, Tanner L, Aw PP, Tan HC, Ng LC, Leo YS, Low JG, Ong A, Ooi EE, Vasudevan SG, Hibberd ML: Genomic epidemiology of a dengue virus epidemic in urban Singapore. J Virol. 2009, 83 (9): 4163-4173. 10.1128/JVI.02445-08.
Ito M, Yamada K, Takasaki T, Pandey B, Nerome R, Tajima S, Morita K, Kurane I: Phylogenetic analysis of dengue viruses isolated from imported dengue patients: possible aid for determining the countries where infections occurred. J Travel Med. 2007, 14 (4): 233-244. 10.1111/j.1708-8305.2007.00130.x.
Yamanaka A, Mulyatno KC, Susilowati H, Hendrianto E, Ginting AP, Sary DD, Rantam FA, Soegijanto S, Konishi E: Displacement of the predominant dengue virus from type 2 to type 1 with a subsequent genotype shift from IV to I in Surabaya, Indonesia 2008–2010. PLoS One. 2011, 6 (11): e27322-10.1371/journal.pone.0027322.
Tissera HA, Ooi EE, Gubler DJ, Tan Y, Logendra B, Wahala WM, de Silva AM, Abeysinghe MR, Palihawadana P, Gunasena S, Tam CC, Amarasinghe A, Letson GW, Margolis HS, De Silva AD: New dengue virus type 1 genotype in Colombo, Sri Lanka. Emerg Infect Dis. 2011, 17 (11): 2053-2055.
Jeong YE, Kim YH, Cho JE, Han MG, Ju YR: Identification of dengue type 1 virus (DENV-1) in Koreans traveling abroad. Public Health Res Perspect. 2011, 2 (1): 34-40. 10.1016/j.phrp.2011.04.002.
Twiddy SS, Farrar JJ, Vinh Chau N, Wills B, Gould EA, Gritsun T, Lloyd G, Holmes EC: Phylogenetic relationships and differential selection pressures among genotypes of dengue-2 virus. Virology. 2002, 298 (1): 63-72. 10.1006/viro.2002.1447.
Migration issues in the Asia Pacific: Issues paper from Malaysia. http://www.unesco.org/most/apmrnwp9.htm,
Sittisombut N, Sistayanarain A, Cardosa MJ, Salminen M, Damrongdachakul S, Kalayanarooj S, Rojanasuphot S, Supawadee J, Maneekarn N: Possible occurrence of a genetic bottleneck in dengue serotype 2 viruses between the 1980 and 1987 epidemic seasons in Bangkok. Thailand. Am J Trop Med Hyg. 1997, 57 (1): 100-108.
Lee KS, Lo S, Tan SS, Chua R, Tan LK, Xu H, Ng LC: Dengue virus surveillance in Singapore reveals high viral diversity through multiple introductions and in situ evolution. Infect Genet Evol. 2012, 12 (1): 77-85. 10.1016/j.meegid.2011.10.012.
Modis Y, Ogata S, Clements D, Harrison SC: A ligand-binding pocket in the dengue virus envelope glycoprotein. Proc Natl Acad Sci U S A. 2003, 100 (12): 6986-6991. 10.1073/pnas.0832193100.
Goncalvez AP, Escalante AA, Pujol FH, Ludert JE, Tovar D, Salas RA, Liprandi F: Diversity and evolution of the envelope gene of dengue virus type 1. Virology. 2002, 303 (1): 110-119. 10.1006/viro.2002.1686.
Modis Y, Ogata S, Clements D, Harrison SC: Structure of the dengue virus envelope protein after membrane fusion. Nature. 2004, 427 (6972): 313-319. 10.1038/nature02165.
Wahala WM, Silva AM: The human antibody response to dengue virus infection. Viruses. 2011, 3 (12): 2374-2395. 10.3390/v3122374.
Lee E, Weir RC, Dalgarno L: Changes in the dengue virus major envelope protein on passaging and their localization on the three-dimensional structure of the protein. Virology. 1997, 232 (2): 281-290. 10.1006/viro.1997.8570.
Butrapet S, Childers T, Moss KJ, Erb SM, Luy BE, Calvert AE, Blair CD, Roehrig JT, Huang CY: Amino acid changes within the E protein hinge region that affect dengue virus type 2 infectivity and fusion. Virology. 2011, 413 (1): 118-127. 10.1016/j.virol.2011.01.030.
Twiddy SS, Holmes EC, Rambaut A: Inferring the rate and time-scale of dengue virus evolution. Mol Biol Evol. 2003, 20 (1): 122-129. 10.1093/molbev/msg010.
Jenkins GM, Rambaut A, Pybus OG, Holmes EC: Rates of molecular evolution in RNA viruses: a quantitative phylogenetic analysis. J Mol Evol. 2002, 54 (2): 156-165. 10.1007/s00239-001-0064-3.
Holmes EC: Patterns of intra- and interhost nonsynonymous variation reveal strong purifying selection in dengue virus. J Virol. 2003, 77 (20): 11296-11298. 10.1128/JVI.77.20.11296-11298.2003.
Martins VE, Alencar CH, Kamimura MT, de Carvalho Araujo FM, De Simone SG, Dutra RF, Guedes MI: Occurrence of natural vertical transmission of dengue-2 and dengue-3 viruses in Aedes aegypti and Aedes albopictus in Fortaleza, Ceara. Brazil. PLoS One. 2012, 7 (7): e41386-10.1371/journal.pone.0041386.
Mulyatno KC, Yamanaka A, Yotopranoto S, Konishi E: Vertical transmission of dengue virus in Aedes aegypti collected in Surabaya, Indonesia, during 2008–2011. Jpn J Infect Dis. 2012, 65 (3): 274-276. 10.7883/yoken.65.274.
Rohani A, Suzilah I, Malinda M, Anuar I, Mohd Mazlan I, Salmah Maszaitun M, Topek O, Tanrang Y, Ooi SC, Rozilawati H, Lee HL: Aedes larval population dynamics and risk for dengue epidemics in Malaysia. Trop Biomed. 2011, 28 (2): 237-248.
Yusoff N, Budin H, Ismail S: Simulation of population dynamics of Aedes aegypti using climate dependent model. World Acad Sci Eng Technol. 2012, 62: 477-482.
Masel J: Genetic drift. Curr Biol. 2011, 21 (20): R837-R838. 10.1016/j.cub.2011.08.007.
Ohta T: Evolution by nearly-neutral mutations. Genetica. 1998, 102–103 (1–6): 83-90.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25 (7): 1253-1256. 10.1093/molbev/msn083.
Pond SL, Frost SD: Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. Bioinformatics. 2005, 21 (10): 2531-2533. 10.1093/bioinformatics/bti320.
Wong SS, Abd-Jamil J, Abubakar S: Antibody neutralization and viral virulence in recurring dengue virus type 2 outbreaks. Viral Immunol. 2007, 20 (3): 359-368. 10.1089/vim.2006.0100.
Clarke DH, Casals J: Techniques for hemagglutination and hemagglutination-inhibition with arthropod-borne viruses. Am J Trop Med Hyg. 1958, 7 (5): 561-573.
World Health Organization: Laboratory diagnosis. Dengue haemorrhagic fever: diagnosis, treatment and control. 1997, Geneva: World Health Organization, 34-47. 2
Okuno Y, Fukunaga T, Srisupaluck S, Fukai K: A modified PAP (peroxidase-anti-peroxidase) staining technique using sera from patients with dengue hemorrhagic fever (DHF): 4 step PAP staining technique. Biken J. 1979, 22 (4): 131-135.
This study was supported by University Malaya Postgraduate Research Fund (PS412-2010A), the Ministry of Higher Education Long-Term Research Grant Scheme (LRGS/TD/2011/UM/Penyakit-Berjangkit) and the Ministry of Science, Technology, and Innovation (Malaysia Genome Institute initiative grant 07-05-MGI-GMB015), Malaysia.
The authors declare that they have no competing interests.
BTT and SSS analyzed and interpreted the data, performed the experiments, and wrote the manuscript. KKT, JJ, MHS, MBD, JAJ, and NAMR performed the virus propagation, amplification and sequencing of the viral envelop gene. NMM helped to analyze and interpret the data. SAB conceived and designed the study, coordinated the experiments, analyzed and interpreted the data, and wrote the manuscript. All authors have read and approved the final manuscript.