- Open Access
Evolution and heterogeneity of multiple serotypes of Dengue virus in Pakistan, 2006–2011
Virology Journalvolume 10, Article number: 275 (2013)
Even though dengue has been recognized as one of the major public health threats in Pakistan, the understanding of its molecular epidemiology is still limited. The genotypic diversity of Dengue virus (DENV) serotypes involved in dengue outbreaks since 2005 in Pakistan is not well studied. Here, we investigated the origin, diversity, genetic relationships and geographic distribution of DENV to understand virus evolution during the recent expansion of dengue in Pakistan.
The study included 200 sera obtained from dengue-suspected patients from 2006 to 2011. DENV infection was confirmed in 94 (47%) sera by a polymerase chain reaction assay. These included 36 (38.3%) DENV-2, 57 DENV-3 (60.6%) and 1 DENV-4 (1.1%) cases. Sequences of 13 whole genomes (6 DENV-2, 6 DENV-3 and 1 DENV-4) and 49 envelope genes (26 DENV-2, 22 DENV-3 and 1 DENV-4) were analysed to determine the origin, phylogeny, diversity and selection pressure during virus evolution.
DENV-2, DENV-3 and DENV-4 in Pakistan from 2006 to 2011 shared 98.5-99.6% nucleotide and 99.3-99.9% amino acid similarity with those circulated in the Indian subcontinent during the last decade. Nevertheless, Pakistan DENV-2 and DENV-3 strains formed distinct clades characterized by amino acid signatures of NS2A-I116T + NS5-K861R and NS3-K590R + NS5-S895L respectively. Each clade consisted of a heterogenous virus population that circulated in Southern (2006–2009) and Northern Pakistan (2011).
DENV-2, DENV-3 and DENV-4 that circulated during 2006–2011 are likely to have first introduced via the southern route of Pakistan. Both DENV-2 and DENV-3 have undergone in-situ evolution to generate heterogenous populations, possibly driven by sustained local DENV transmission during 2006–2011 periods. While both DENV-2 and DENV-3 continued to circulate in Southern Pakistan until 2009, DENV-2 has spread in a Northern direction to establish in Punjab Province, which experienced a massive dengue outbreak in 2011.
Dengue is an acute febrile illness with a clinical spectrum of mild flu-like illness (dengue fever or DF) to fatal complications of severe dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). The disease is caused by the Dengue virus (DENV) complex that consists of 4 genetically and immunogenically distinct serotypes. DENV is a single stranded positive sense RNA virus with a 11.8 kb genome flanked by 5’ and 3’ untranslated regions and a single coding region for three structural proteins: capsid (C), the precursor of membrane (prM) and envelope (E); and seven nonstructural proteins: NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5 .
DENV is transmitted to humans through the bites of infected female mosquitoes of Aedes aegypti and Ae. albopictus. The transmission of DENV has increased in recent years in urban and semi-urban endemic settings, especially in Americas, South and South-east Asia and the Western Pacific. The magnitude, distribution, and clinical severity of dengue outbreaks have been alarmingly high in countries such as India , Sri Lanka , Nepal  Bangladesh  and Pakistan  in the Indian subcontinent during the last decade.
Pakistan is endemic to all four serotypes of DENV circulating throughout the year with a peak incidence during the post monsoon period between October and December [7, 8]. Factors such as crowded cities, unsafe drinking water, inadequate sanitation and large number of refugees facilitate the spread of dengue in different parts of the country, resulting in increased morbidity and mortality. It is believed that DENV was first introduced into Pakistan through the importation of tyres containing eggs of infected mosquitoes at Karachi sea port . Although serological evidence of DENV infections in Pakistan dates back to 1968 , the first confirmed outbreak associated with DHF occurred in the southern Pakistan city of Karachi in 1994 . Serological studies confirmed the circulation of both DENV-1 and DENV-2 during the 1994 outbreak . Since then, the disease has emerged as a major public health problem in the country . DENV-3 was first reported during the 2005–2006 outbreak in Karachi [14, 15]. By 2007, dengue started to emerge in the Northern Pakistan. DENV-2 and DENV-3 have been the dominant serotypes in Lahore from 2007 to 2009 . Dengue showed a resurgence in November 2010, especially in Sindh, Punjab and Khyber Pakhtunkhwa regions subsequent to massive floods in July same year . The outbreak escalated in Lahore, Punjab, in 2011 with 22,562 cases and 363 deaths due to severe DHF .
So far, there have been a few, comprehensive molecular epidemiological studies that describe DENV circulating in Pakistan. There is only one genome-wide analysis that describes Pakistani DENV-2 circulated during the 2011 outbreak . The studies of molecular epidemiology and evolutionary genetics are important to predict the origin and spread of viruses, to strengthen our understanding on the pathogenesis of disease, the cause of epidemics and the genetic basis of virulence. RNA viruses such as DENV evolve rapidly  and on rare occasions, certain mutations lead to phenotypic changes in the viruses that alter their potential to cause outbreaks associated with severe disease . Genotypic characterization has been a useful tool in determining the evolutionary origin of the DENV, identifying the circulating virus strains in an endemic area and detecting the introduction of new genotypes. In recent years, envelope (E) gene sequencing has been widely used to assess the phylogenetic relationships among DENV isolates . In Pakistan, sequencing of NS3 gene  and C-prM gene [14, 16] junction has mainly been used for the genetic analyses of DENV. Khan and co-workers recently performed a more comprehensive analysis of DENV2 strains associated with 2011 outbreak in Lahore . Nevertheless, the use of relatively short (up to 400 bp) and small number of sequences in available studies have limited the resolution of phylogenetic analyses aimed at determining the origin and evolution of DENV in Pakistan during the last decade.
Therefore, in the present study, we conducted a detailed molecular characterization of Pakistan DENV strains using the complete E gene and whole genome analyses. DENV isolates obtained from different parts of Pakistan between 2006 and 2011 were studied, in order to obtain a temporal evolutionary perspective of DENV in the country over a 6-year period. We describe the likely origin, evolution and geographic distribution of DENV that explain the recent changes of dengue epidemiology in Pakistan.
Results and discussion
Out of 200 sera analysed, 94 (47%) were positive for DENV by polymerase chain reaction (PCR). These included 36 DENV-2 (38.3%), 57 DENV-3 (60.6%) and 1 DENV-4 (1.1%) cases (Table 1). None of the samples was positive for DENV-1. DENV could be isolated only from 25 PCR-positive sera (26.6%). These included 12 DENV-2 (48.0%), 12 DENV-3 (48.0%) and 1 DENV-4 (4.0%) cases. The reduced success rate of virus isolation in PCR positive sera is due to collection of the majority of sera during the late viraemic phase and further reduction in virus titres during long-term storage, transfer and freeze-thawing of samples. Among the DENV-2 samples, 7 were classified as DF, 16 as DHF and 3 as DSS, whereas for DENV-3, there were 18 DF and 4 DHF cases. Only 1 DENV-4 sample classified as DF was included in this study (Additional file 1).
We detected only DENV-3 among sera collected during the 2006–2007 periods, whereas DENV-2, DENV-3 and DENV-4 were observed in sera collected after 2008. This is in agreement with previous findings that DENV-2 became dominant during the 2007–2009 periods . The absence of DENV-1 in our cohort was not surprising as it has not been documented in Pakistan between 2006 and 2009 [9, 16] and was a minor serotype in 2010 in Northern parts of the country . During the 2006–2009 periods, 87.7% (n=71) of analysed sera were collected in Karachi, the capital of Sindh province in Southern Pakistan. On the other hand, 69.2% (n=9) of cases in 2011 were from Lahore in Punjab province of the Northern Pakistan where a major dengue outbreak was reported in the same year. This case distribution pattern reflects the dengue epidemiology in Pakistan, where the disease was initially restricted to the southern part of the country, particularly the port city of Karachi, and subsequently expanded to the north, causing major epidemics. The spread of dengue to other parts of the country is believed to be either due to the movement of people carrying the virus or mosquito eggs carrying vertically-transmitted DENV from dengue-affected areas .
In order to understand the origin and evolution of different DENV serotypes during their spread in Pakistan, we performed a detailed analysis of 13 full genomes (6 DENV-2, 6 DENV-3 and 1 DENV-4) and 49 envelope (E) gene sequences detected in Southern and Northern parts during the 2006–2011 period. All E gene sequences were obtained directly from PCR-positive sera and comprised 26 DENV-2, 22 DENV-3 and 1 DENV-4 strains. Full genome sequencing was performed on viruses isolated from sera.
Characteristics of DENV-2 in Pakistan
Phylogenetic analysis of 26 DENV-2 E gene sequences obtained between 2008 and 2011 revealed that all viruses belonged to Indian subcontinent lineage of the cosmopolitan genotype. Although there were more DHF compared to DF among the samples collected, there was no unique signature in the E gene sequences that could differentiate between DF and DHF/DSS samples. With the exception of 1 isolate (D2-PS51) that was closely related to a previous strain from Saudi Arabia, the remaining 25 isolates formed a distinct monophyletic clade with 78% bootstrap support (Figure 1). In DENV-2 whole genome analysis, the same clade achieved 100% bootstrap support, strengthening the credibility of the clade further (data not shown). This monophyletic clade included DENV-2 isolates obtained from 14 sera collected in Southern Pakistan (Karachi) in 2008 and 2009 as well as 11 sera collected in Northern Pakistan in 2011. Among complete E gene sequences available in GenBank database, these 25 viruses shared the highest nucleotide (98.8-99.5%) and amino acid (99.5-100%) similarity with those reported from Sri Lanka [GenBank: GQ252676 and GQ252677] and India [GenBank: DQ448231] during the 2003–04 periods. The nucleotide (99.1-99.4%) and amino acid (99.7-99.9%) similarity was even higher at whole genome level between Pakistan and 2003–2004 Indian/Sri Lankan strains. Furthermore, whole genome analysis revealed a unique combination of 8 amino acid substitutions that was shared between the monophyletic clade of Pakistan DENV-2 and Indian/Sri Lankan isolates (Table 2). These findings indicated the close genetic relationship among DENV-2 strains circulated in the Indian subcontinent during the last decade. Even though the earliest DENV-2 isolate of our study cohort was in 2008, estimation of time to the most recent common ancestor (TMRCA) dated the origin of Pakistan DENV-2 clade to early 2006 [95% highest probability density (HPD) 2004.5 – 2007.4]. Given the fact that genetically similar strains existed in the region even in 2003, it is likely that the 2008–2011 clade of Pakistan DENV-2 viruses descended from an ancestral strain that was previously circulating in other parts of the Indian subcontinent region. Further evolution in subsequent years has deviated Pakistan DENV-2 from its common ancestor to form a monophyletic clade by 2008. During this process, Pakistan DENV-2 strains have acquired 2 additional fixed amino acid substitutions (NS2A-I116T and NS5-K861R), of which NS5-K861R is unique to Pakistan strains (Table 2). In DENV-2, residue NS5-861 is a strictly conserved Lysine. NS5-861 resides in the “thumb” domain of the C-terminal end of RNA dependent RNA polymerase (RdRp) of DENV . This residue sits outside structural elements of thumb domain that form the RNA template tunnel (residues 740–747 and 782–809) and a zinc binding site (residues 712, 714, 728 and 847). Therefore, even though the functional implications of the Lysine to Arginine substitution in our isolates are not clear, major structural changes are not expected as the K861R substitution does not alter the amino acid properties at the residue.
Furthermore, our findings indicated an in-situ spatio-temporal evolutionary pattern of Pakistan DENV-2 during the period from 2008 to 2011. This was evident by the step-ladder pattern topology of the 2008–2011 Pakistan clade (Figure 1), in which the 2009 viruses in Southern Pakistan were ancestral to the distinct sub-clade (99% bootstrap support) formed by 2011 viruses in the North. This 2011 sub-clade possessed a unique substitution (prM-S70A) that was shown to be under the diversifying selection (p=0.05, Table 2) by the mixed effect model of evolution (MEME) method. In contrast, single likelihhod ancestor counting (SLAC), fixed effects likelihood (FEL)  and internal fixed effects likelihood (IFEL)  methods showed the same residue to be under purifying selection. MEME is one of the methods available in Hyphy package to perform the codon-based selection pressure analysis . A simulation study using MEME has shown that MEME is preferable to other methods as it is capable of identifying not only pervasive positive selection at the level of an individual site, but also episodic (transient) positive selection affecting only a subset of lineages . Episodic selection is often missed by other methods that rely on the mean non-synonymous to synonymous ratio (dN/dS) over a period of time . Sequences from another recent study  have also shown prM-S70A to be fixed and unique to isolates from Lahore in 2011. Therefore, we postulate that amino acid residue prM-70 is likely to have been selected for during the DENV-2 outbreak in Northern Pakistan in 2011. As the support for diversifying selection in MEME analysis is marginally significant (p=0.05) and we obtained contrasting outputs from SLAC, FEL and IFEL methods, it is difficult to conclude whether the fixation of prM-S70A among Pakistani DENV-2 strains is either due to positive selection or genetic drift. Nevertheless, codon 70 constitutes the variable, middle amino acid residue of the N-glycosylation site sequon (N-S-T) of DENV-2 prM protein . N-glycosylation of prM protein is believed to play an important role in virion assembly and transport before release . The prM-70 residue is a strictly-conserved Serine in human DENV-2 isolates, except for two sylvatic viruses (N-A-T) reported from Malaysia in 2008  and 1970 . Therefore, whether this substitution from a hydrophilic Serine to a hydrophobic Alanine at prM-70 residue alters the conformation and thereby the functionality of the N-glycosylation site providing a selective advantage to DENV is of interest and necessitates further investigation.
As illustrated in Figure 1, virus populations from Southern and Northern areas are heterogenous, indicating the possibility of sustained virus transmission in both areas driving the viral diversity. While the majority of 2011 viruses belonged to a distinct clade, some remained within a 2009 virus cluster (Figure 1). Based on these observations, it is likely that DENV-2 first expanded in Southern Pakistan and subsequently spread to the Northern region where certain sub-population/s evolved further to become a distinct viral population in 2011. Therefore, the overall epidemiological understanding of South to North expansion of dengue in Pakistan is supported by our genetic analysis of DENV-2.
Characteristics of DENV-3 in Pakistan
The emergence of DENV-3 in Pakistan was first reported during the 2005 DHF outbreak in Karachi . Previous studies have revealed that DENV-3 sequences from the 2005 and 2006 outbreaks were most closely related to those circulating in India [14, 32] and belonged to genotype III . In support, all 22 DENV-3 E gene sequences from our study also belonged to the genotype III (Figure 2). Furthermore, all of our study isolates that were collected in southern Pakistan (Karachi and Hyderabad) from 2006 to 2009 shared the highest genetic similarity with two Indian DENV-3 strains [GenBank: AY770511 and FJ644564] and one Chinese strain [GenBank: GU363549]. Of these, one of the Indian strains [GenBank: AY770511] was reported in 2003  whereas the remaining Indian  and Chinese strains were reported in 2007 and 2009 respectively. Therefore, it was evident that ancestral strains of Pakistan DENV-3 were present in the Indian subcontinent before their emergence in Pakistan during the 2005–2006 periods. At E gene level, nucleotide and amino acid similarities between the Pakistan and the Indian strain in 2003 [GenBank: AY770511] were 98.8-99.2% and 99.5-100%, respectively. The respective whole genome similarities between Pakistan DENV-3 isolates (n=6) and the same Indian strain were 98.9-99.2% (nucleotide) and 99.6-99.8% (amino acid). Sharing of a unique combination of 4 amino acid substitutions (NS1-L178M, NS4B-M151I, NS5-T188A and NS5-Q562L) between Pakistan and Indian isolates (Table 3) further supported the previous notion that Pakistan DENV-3 isolates are genetically related to those circulating in India in preceding years. Of these, NS1-L178M was found to be a positively selected site (p: <0.01. Table 3) in the MEME-based selection pressure analysis, suggesting a lineage-specific selective advantage of DENV-3 that circulated in the Indian subcontinent before 2009. SLAC, FEL and IFEL-based slection pressure analyses predicted the same site to be neutral. The residue 178 sits in domain II of flavivirus NS1 protein , adjacent to one of the 12 strictly conserved Cysteine residues (DENV-3-Cys179) that form disulfide linkages . Disulfide bridges are believed to play an important role in NS1 dimerization . Interestingly, residue 178 appears to follow a lineage-specific variability as DENV-3 genotype III lineage in South-east Asia (Vietnam, Cambodia, Thailand) shows a fixed NS1-L178S substitution. This observation, together with our data, indicates the likelihood of episodic diversifying selection at DENV-3 NS1-178 residue, as predicted by the MEME analysis for positive selection pressure.
Overall, these findings suggested an Indian subcontinent ancestry of 2006–2009 Pakistan DENV-3. Interestingly, E gene-based TMRCA analysis dated the origin of 2006–2009 Pakistan DENV-3 to late 2002 (95% HPD 2001.6-2004.0). The available evidence of a genetically closest Indian DENV-3 isolate in 2003, therefore, confirms the circulation of similar virus strains in other parts of the Indian subcontinent at the time of emergence of the common ancestor of the 2006–2009 DENV-3 strains in Pakistan. However, despite their high similarity with Indian strains, 2006–2009 Pakistan DENV-3 formed a distinct E gene cluster with 76% bootstrap support (Figure 2). In DENV-3 whole genome phylogeny, the bootstrap support of the same clade improved to 100% (data not shown). Our whole genome analysis revealed that this genetic deviation of Pakistan strains from Indian sub continent isolates was driven by two additional fixed amino acid substitutions (NS3-K590R and NS5-S895L). Of these, NS3-K590R is unique to Pakistan DENV-3 isolates. NS3-K590R falls within a loop structure between alpha helices 9 and 10 of flavivirus helicase domain 3 . Even though domain 3 does not house major functional elements of helicase activity, it contributes to structural stability and the formation of RNA binding cleft [36–39]. Domain 3 is also believed to act as the RNA-dependent RNA polymerase binding site . As the residue 590 does not carry any known structural significance, the effect of NS3-K590R substitution on the domain 3 functionality of our study isolates is unclear. In summary, our analysis indicates that 2006 DENV-3 strains emerged from an ancestral strain that introduced into Southern Pakistan, possibly between 2001 and 2004. Pakistan DENV-3 strains evolved further to emerge as a distinct clade during the outbreak in Karachi in 2006.
Characteristics of DENV-4 in Pakistan
DENV-4 was first reported in Lahore during the 2008 dengue outbreak in Punjab, Pakistan . Phylogenetic analysis showed that DENV-4 from the Indian subcontinent formed a distinct group within genotype I (Figure 3). A DENV-4 isolate obtained from Karachi in 2009 during our study clustered with another Pakistan sequence reported previously in the same year . These Pakistan isolates shared the highest nucleotide (98.5%) and amino acid (99.3%) with a DENV-4 strain (GenBank: HM237349) from Andra Pradesh in Southren India in 2007 . Therefore, it is evident that DENV-4 strains closely related to those in Pakistan in 2009 have previously been circulating in the Indian subcontinent. Nevertheless, more sequence data from a wider geographical scale in the Indian subcontinent region is needed to confirm the spread and evolutionary history of DENV-4 in Pakistan.
DENV-2, DENV-3 and DENV-4 have co-circulated in Pakistan from 2008 to 2011. All three serotypes in Pakistan share an Indian subcontinent ancestry and are likely to have been introduced first into Southern Pakistan. DENV-2 and DENV-3 have undergone in situ evolution to emerge as distinct, heterogenous virus populations during the same period. While both DENV-2 and DENV-3 continued to circulate in Southern Pakistan until 2009, DENV-2 has subsequently spread in a Northern direction, a move that caused a massive outbreak in Punjab province in 2011. The study, therefore, highlights the implication of human migration and virus evolution in the expansion of dengue to areas where the disease was not previously a major public health threat, and the emergence of new strains that may be associated with severe outbreaks.
Aga Khan University Hospital (AKUH) is a 550 bed tertiary care hospital located in Karachi, Pakistan. The hospital is equipped with a state of the art clinical laboratory that caters not only for the in-patients but also reaches out to the community. It has a unique network of 200 sample collection units spread throughout the country. The laboratory receives samples from inpatients and outpatients of AKUH and from other clinics and hospitals in all major cities of the country. Karachi witnessed its first ever severe dengue outbreak in 2005–06 during which the clinical laboratory of AKUH was at the forefront of dengue diagnosis.
The patients presenting to the clinics and in-patients of the AKUH Karachi, with the history of fever, headache, joint pain, vomiting and diarrhoea were recruited based on presumptive clinical diagnosis of DENV infection. In addition, samples obtained from out patients at collection units were also included. After taking informed consent, a complete Performa about patient’s demographics and clinical presentation were filled. Serum samples of each patient were stored in three aliquots at −80°C. Two hundred dengue-suspected sera collected from 2006 to 2011 were retrieved from the serum bank at AKUH and were shipped to Environmental Health Institute, National Environment Agency, Singapore for the molecular epidemiological analysis. The work presented here was approved by the Ethical Review Committee, Aga Khan University (AKU), Pakistan.
DENV was isolated from sera using the Ae. albopictus clone C6/36 mosquito cell line (ATCC CRL-1660). Briefly, the monolayer of C6/36 was inoculated with serum in L-15 maintenance medium containing 3 % FCS at 33°C to allow for virus adsorption and replication. The infected fluid was harvested after 5–10 days. The presence of DENV in cell supernatants was confirmed by an immunofluoresent assay (IFA) during which cells were reacted with either DENV group-specific or serotype-specific monoclonal antibodies derived from hybridoma cultures (ATCC HB-46, HB-47, HB-48, and HB-49). The fluorescein isothiocyanate-conjugated goat anti-mouse antibody was used as the detector. A maximum of three passages was done for each sample based on the original virus load.
RNA extraction and conventional PCR for detection and serotyping of Dengue virus
Viral RNA was extracted either from sera or culture supernatants using the QIAGEN QIAamp viral RNA minikit (QIAGEN, Hilden, Germany) according to the manufacturer’s recommendations.
Conventional semi-nested PCR was performed using a modified procedure described by Lanciotti and colleagues . A one-step RT-PCR was performed using the AccessQuickTM RT-PCR System (Promega, Madison, WI, USA) in a 50 μl reaction volume containing 1X AccessQuick TM Master Mix, 2.5 units of AMV Reverse Transcriptase and 0.25 μM (each) of Primers D1 and D2 . The reactions were carried out in Tgradient PCR Thermocycler (Biometra, Gottingen, Germany) using the following programme; reverse transcription at 45°C for 30 min, inactivation at 94°C for 2 min, and PCR amplification of 40 cycles under the following conditions: 94°C for 30 sec, 55°C for 30 sec, 72°C for 1 min, and a final extension at 72°C for 5 min. Five microlitres of the 1:100 dilution of the first-round PCR product was further amplified in a 50 μl reaction mixture containing 1X Green GoTaq® Flexi Buffer, 1.5 mM MgCl2, 0.2 mM dNTPs, 1.25 units of GoTaq® DNA Polymerase (Promega, Madison, WI, USA) and 0.5 μM of each primer D1, TS1, TS2, TS3 and TS4 . After denaturation at 95°C for 2 min, the second round amplification was performed for 30 cycles (95°C for 30 sec, 55°C for 30 sec, 72°C for 1 min, and a final extension at 72°C for 10 min). Amplified products were visualized in 2% agarose gels stained with GelRed (Biotium Inc., USA).
PCR amplification of envelope gene for sequencing
The synthesis of complementary DNA (cDNA) from RNA extracted from sera was carried out using Superscript™ III First-strand Synthesis System (Invitrogen, Carlsbad, USA). The complete envelope (E) gene (~1.5 kb) was amplified by PCR using 0.5 μM of DENV serotype-specific primers (Additional file 2) and 1X Phusion™ Flash High-Fidelity PCR Master Mix (Finnzymes, Lafayette, CO). The amplification protocol is as follows: initial denaturation at 98°C for 5 sec, 35 cycles of denaturation at 98°C for 5 sec, annealing at 60°C for 8 sec, extension at 72°C for 25 sec and final extension at 72°C for 1 min. Semi-nested PCR (~1.4 kb, Additional file 2) was performed as above if no desired product yielded in the first-round amplification. Amplified products were visualized in 2% agarose gels stained with GelRed (Biotium Inc., USA).
Whole genome sequencing
A total of 13 DENV isolated from DENV-2 (n=6), DENV-3 (n=6) and DENV-4 (n=1) positive sera were selected for full genome sequencing. Isolates were selected to represent different groups of viruses based on E gene phylogeny. Amplification and sequencing primers were designed to cover the entire genome of different serotypes in overlapping fragments (Additional files 3, 4 and 5). Additional two sets of primers (Additional file 6) were used to capture 5' and 3' untranslated regions of complete genome for all four DENV serotypes. The PCR reaction was performed using the Phusion™ Flash High-Fidelity PCR Master Mix (Finnzymes, Lafayette, CO) as described for E gene amplification. Amplified products were visualized in 2% agarose gels stained with GelRed (Biotium Inc., USA).
Purification and sequencing of PCR products
PCR products were purified using the Qiaquick PCR Purification kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions. Sequencing of purified PCR products was perfromed by a commercial sequencing company (AITbiotech Pte Ltd, Singapore) according to the BigDye Terminator Cycle Sequencing kit (Applied Biosystems, USA) protocol. Seqeunces were deposited in Genbank database under the accession numbers KF041212 - KF041260.
Phylogenetic and mutation analyses of envelope genes and complete genomes
Nucleotide sequences were assembled using the Lasergene package version 8.0 (DNASTAR Inc., Madison, WI, USA). Contiguous sequences were aligned using ClustalX program  and compared with published sequences of DENV isolates in Genbank database. Phylogenetic analysis of E gene sequences was performed in MEGA5 program  using the maximum-likelihood method based on the general time reversible (GTR) model with gamma distribution and invariant sites. The robustness of the original tree was tested with 1000 bootstrap replications. Comparison of amino acid sequences deduced from nucleotide alignments was performed by using the BioEdit v7.0.5 software to identify unique mutations in study isolates.
Estimation of the time of viral emergence for DENV-2 and DENV-3
The time to the most recent common ancestor (TMRCA) was inferred using the Bayesian Markov Chain Monte Carlo (MCMC) method implemented in the BEAST package v1.7.4 [46, 47]. A GTR substitution model with gamma distribution of rate variation among sites and a proportion of invariable sites, an uncorrelated log-normal relaxed molecular clock model and a Bayesian skyline coalescent model were used for this analysis. One hundred million generations of MCMC chains were run with sampling at every 10,000 generations. The MCMC output was analysed using the Tracer Version 1.5 programme (http://tree.bio.ed.ac.uk/software/tracer/). The uncertainty in parameter estimates was expressed as 95% highest probability density (HPD). All runs of MCMC were repeated to ensure convergence with effective sample size (ESS) of >200.
Test for selection pressure
Genome-wide dN/dS ratios were computed in Hyphy 2.1.2 standalone version  using the single likelihood ancestor counting (SLAC), fixed effects likelihood (FEL), internal fixed effects likelihood (IFEL) and the mixed effect model of evolution (MEME) methods . The Muse and Gaut (MG94) codon-based substitution model and GTR-based maximum likelihood phylogeny with gamma (4) distributed invariant sites, constructed using MEGA 5 , were used and significance levels were set at p<0.05.
Dengue hemorrhagic fever
Dengue shock syndrome
Fixed effects likelihood
General time reversible
Highest probability density
Internal fixed effects likelihood
Monte Carlo Markov Chain
mixed effect model of evolution
Polymerase chain reaction
Single likelihood ancestor counting
Time to the most recent common ancestor.
Chambers TJ, Hahn CS, Galler R, Rice CM: Flavivirus genome organization, expression, and replication. Annu Rev Microbiol. 1990, 44: 649-688. 10.1146/annurev.mi.44.100190.003245.
Gupta B, Reddy BP: Fight against dengue in India: progresses and challenges. Parasitol Res. 2013, 112: 1367-1378. 10.1007/s00436-013-3342-2.
Dissanayake VH, Gunawardena ND, Gunasekara NC, Siriwardhana DR, Senarath N: Shift in the transmission pattern of dengue serotypes and concurrent infection with more than one dengue virus serotype. Ceylon Med J. 2011, 56: 176-178.
Dumre SP, Shakya G, Na-Bangchang K, Eursitthichai V, Grams HR, Upreti SR, Ghimire P, Kc K, Nisaluk A, Gibbons RV, Fernandez S: Dengue virus and Japanese Encephalitis virus epidemiological shifts in Nepal: a case of opposing trends. Am J Trop Med Hyg. 2013, 88: 677-680. 10.4269/ajtmh.12-0436.
Raheel U, Faheem M, Riaz MN, Kanwal N, Javed F, Zaidi N, Qadri I: Dengue fever in the Indian Subcontinent: an overview. J Infect Dev Ctries. 2011, 5: 239-247.
Rasheed SB, Butlin RK, Boots M: A review of dengue as an emerging disease in Pakistan. Public Health. 2013, 127: 11-17. 10.1016/j.puhe.2012.09.006.
Jahan F: Dengue Fever (DF) in Pakistan. Asia Pac Fam Med. 2011, 10: 1-10.1186/1447-056X-10-1.
Khan E, Kisat M, Khan N, Nasir A, Ayub S, Hasan R: Demographic and clinical features of dengue fever in Pakistan from 2003–2007: a retrospective cross-sectional study. PLoS One. 2010, 5: 12505-10.1371/journal.pone.0012505.
Idrees S, Ashfaq UA: A brief review on dengue molecular virology, diagnosis, treatment and prevalence in Pakistan. Genet Vaccines Ther. 2012, 10: 6-10.1186/1479-0556-10-6.
Hayes CG, Baqar S, Ahmed T, Chowdhry MA, Reisen WK: West Nile virus in Pakistan. 1. Sero-epidemiological studies in Punjab Province. Trans R Soc Trop Med Hyg. 1982, 76: 431-436. 10.1016/0035-9203(82)90130-4.
Chan YC, Salahuddin NI, Khan J, Tan HC, Seah CL, Li J, Chow VT: Dengue haemorrhagic fever outbreak in Karachi, Pakistan, 1994. Trans R Soc Trop Med Hyg. 1995, 89: 619-620. 10.1016/0035-9203(95)90412-3.
Akram DS, Igarashi A, Takasu T: Dengue virus infection among children with undifferentiated fever in Karachi. Indian J Pediatr. 1998, 65: 735-740. 10.1007/BF02731055.
Paul RE, Patel AY, Mirza S, Fisher-Hoch SP, Luby SP: Expansion of epidemic dengue viral infections to Pakistan. Int J Infect Dis. 1998, 2: 197-201. 10.1016/S1201-9712(98)90052-2.
Jamil B, Hasan R, Zafar A, Bewley K, Chamberlain J, Mioulet V, Rowlands M, Hewson R: Dengue virus serotype 3, Karachi, Pakistan. Emerg Infect Dis. 2007, 13: 182-183. 10.3201/eid1301.060376.
Khan E, Hasan R, Mehraj V, Nasir A, Siddiqui J, Hewson R: Co-circulations of two genotypes of dengue virus in 2006 out-break of dengue hemorrhagic fever in Karachi, Pakistan. J Clin Virol. 2008, 43: 176-179. 10.1016/j.jcv.2008.06.003.
Fatima Z, Idrees M, Bajwa MA, Tahir Z, Ullah O, Zia MQ, Hussain A, Akram M, Khubaib B, Afzal S, et al: Serotype and genotype analysis of dengue virus by sequencing followed by phylogenetic analysis using samples from three mini outbreaks-2007-2009 in Pakistan. BMC Microbiol. 2011, 11: 200-10.1186/1471-2180-11-200.
Khan MA, Ellis EM, Tissera HA, Alvi MY, Rahman FF, Masud F, Chow A, Howe S, Dhanasekaran V, Ellis BR, Gubler DJ: Emergence and diversification of dengue 2 cosmopolitan genotype in pakistan, 2011. PLoS One. 2013, 8: e56391-10.1371/journal.pone.0056391.
Steinhauer DA, Holland JJ: Rapid evolution of RNA viruses. Annu Rev Microbiol. 1987, 41: 409-433. 10.1146/annurev.mi.41.100187.002205.
Grenfell BT, Pybus OG, Gog JR, Wood JL, Daly JM, Mumford JA, Holmes EC: Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004, 303: 327-332. 10.1126/science.1090727.
Twiddy SS, Holmes EC, Rambaut A: Inferring the rate and time-scale of dengue virus evolution. Mol Biol Evol. 2003, 20: 122-129. 10.1093/molbev/msg010.
Chow VT, Yong R, Chan YC: Sequence analyses of NS3 genes of recent Pakistan and Singapore strains of dengue virus types 1 and 2. Res Virol. 1997, 148: 17-20. 10.1016/S0923-2516(97)81907-3.
Mahmood N, Rana MY, Qureshi Z, Mujtaba G, Shaukat U: Prevalence and molecular characterization of dengue viruses serotypes in 2010 epidemic. Am J Med Sci. 2012, 343: 61-64. 10.1097/MAJ.0b013e3182217001.
Yap TL, Xu T, Chen YL, Malet H, Egloff MP, Canard B, Vasudevan SG, Lescar J: Crystal structure of the dengue virus RNA-dependent RNA polymerase catalytic domain at 1.85-angstrom resolution. J Virol. 2007, 81: 4753-4765. 10.1128/JVI.02283-06.
Cardosa J, Ooi MH, Tio PH, Perera D, Holmes EC, Bibi K, Abdul Manap Z: Dengue virus serotype 2 from a sylvatic lineage isolated from a patient with dengue hemorrhagic fever. PLoS Negl Trop Dis. 2009, 3: e423-10.1371/journal.pntd.0000423.
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-9595. 10.1128/JVI.02776-06.
Blok J, McWilliam SM, Butler HC, Gibbs AJ, Weiller G, Herring BL, Hemsley AC, Aaskov JG, Yoksan S, Bhamarapravati N: Comparison of a dengue-2 virus and its candidate vaccine derivative: sequence relationships with the flaviviruses and other viruses. Virology. 1992, 187: 573-590. 10.1016/0042-6822(92)90460-7.
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-1222. 10.1093/molbev/msi105.
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-10.1371/journal.pcbi.0020062.
Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005, 21: 676-679. 10.1093/bioinformatics/bti079.
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-10.1371/journal.pgen.1002764.
Murray GG, Kosakovsky Pond SL, Obbard DJ: Suppressors of RNAi from plant viruses are subject to episodic positive selection. Proc Biol Sci. 2013, 280: 20130965-10.1098/rspb.2013.0965.
Dash PK, Parida MM, Saxena P, Abhyankar A, Singh CP, Tewari KN, Jana AM, Sekhar K, Rao PV: Reemergence of dengue virus type-3 (subtype-III) in India: implications for increased incidence of DHF & DSS. Virol J. 2013, 3: 55-
Sharma S, Dash PK, Agarwal S, Shukla J, Parida MM, Rao PV: Comparative complete genome analysis of dengue virus type 3 circulating in India between 2003 and. J Gen Virol. 2008, 92: 1595-1600.
Muller DA, Young PR: The flavivirus NS1 protein: molecular and structural biology, immunology, role in pathogenesis and application as a diagnostic biomarker. Antiviral Res. 2013, 98: 192-208. 10.1016/j.antiviral.2013.03.008.
Pryor MJ, Wright PJ: The effects of site-directed mutagenesis on the dimerization and secretion of the NS1 protein specified by dengue virus. Virology. 1993, 194: 769-780. 10.1006/viro.1993.1318.
Wu J, Bera AK, Kuhn RJ, Smith JL: Structure of the Flavivirus helicase: implications for catalytic activity, protein interactions, and proteolytic processing. J Virol. 2005, 79: 10268-10277. 10.1128/JVI.79.16.10268-10277.2005.
Mastrangelo E, Milani M, Bollati M, Selisko B, Peyrane F, Pandini V, Sorrentino G, Canard B, Konarev PV, Svergun DI, et al: Crystal structure and activity of Kunjin virus NS3 helicase; protease and helicase domain assembly in the full length NS3 protein. J Mol Biol. 2007, 372: 444-455. 10.1016/j.jmb.2007.06.055.
Lescar J, Luo D, Xu T, Sampath A, Lim SP, Canard B, Vasudevan SG: Towards the design of antiviral inhibitors against flaviviruses: the case for the multifunctional NS3 protein from Dengue virus as a target. Antiviral Res. 2008, 80: 94-101. 10.1016/j.antiviral.2008.07.001.
Luo D, Xu T, Watson RP, Scherer-Becker D, Sampath A, Jahnke W, Yeong SS, Wang CH, Lim SP, Strongin A, et al: Insights into RNA unwinding and ATP hydrolysis by the flavivirus NS3 protein. EMBO J. 2008, 27: 3209-3219. 10.1038/emboj.2008.232.
Humayoun MA, Waseem T, Jawa AA, Hashmi MS, Akram J: Multiple dengue serotypes and high frequency of dengue hemorrhagic fever at two tertiary care hospitals in Lahore during the 2008 dengue virus outbreak in Punjab, Pakistan. Int J Infect Dis. 2013, 14 (Suppl 3): 54-59.
Huang JH, Su CL, Yang CF, Liao TL, Hsu TC, Chang SF, Lin CC, Shu PY: Molecular characterization and phylogenetic analysis of dengue viruses imported into Taiwan during 2008–2010. Am J Trop Med Hyg. 87: 349-358.
Dash PK, Sharma S, Srivastava A, Santhosh SR, Parida MM, Neeraja M, Subbalaxmi MV, Lakshmi V, Rao PV: Emergence of dengue virus type 4 (genotype I) in India. Epidemiol Infect. 139: 857-861.
Lanciotti RS, Calisher CH, Gubler DJ, Chang GJ, Vorndam AV: Rapid detection and typing of dengue viruses from clinical samples by using reverse transcriptase-polymerase chain reaction. J Clin Microbiol. 1992, 30: 545-551.
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-2948. 10.1093/bioinformatics/btm404.
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-2739. 10.1093/molbev/msr121.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Drummond AJ, Suchard MA, Xie D, Rambaut A: Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012, 29: 1969-1973. 10.1093/molbev/mss075.
The study was funded by Aga Khan University Research Council Grant and National Environment Agency, Singapore. The funding sources of this study had no role in the study design, data collection, data analysis, data interpretation, writing of the report, or in the decision to submit the paper for publication.
Authors declare that they have no competing interests.
EK and NLC conceived the study. EK actively recruited patient and analyzed patient demographic and clinical data, AN executed PCR, ZH contributed intellectual input in optimization of the PCR assays at AKU. CK performed virus isolation and generated sequence data. CK, HCH and KSL performed sequence analyses. EK, CK, HCH and KSL wrote the manuscript. All authors read and approved the final manuscript.
Carmen Koo, Amna Nasir contributed equally to this work.