Open Access

Evolution and heterogeneity of multiple serotypes of Dengue virus in Pakistan, 2006–2011

  • Carmen Koo1,
  • Amna Nasir2,
  • Hapuarachchige Chanditha Hapuarachchi1,
  • Kim-Sung Lee1Email author,
  • Zahra Hasan2,
  • Lee-Ching Ng1 and
  • Erum Khan2Email author
Contributed equally
Virology Journal201310:275

https://doi.org/10.1186/1743-422X-10-275

Received: 16 May 2013

Accepted: 30 August 2013

Published: 4 September 2013

Abstract

Background

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.

Methods

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.

Results

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).

Conclusions

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.

Keywords

Dengue Evolution Diversity Molecular epidemiology Phylogenetics Genome Pakistan

Background

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 [1].

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 [2], Sri Lanka [3], Nepal [4] Bangladesh [5] and Pakistan [6] 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 [9]. Although serological evidence of DENV infections in Pakistan dates back to 1968 [10], the first confirmed outbreak associated with DHF occurred in the southern Pakistan city of Karachi in 1994 [11]. Serological studies confirmed the circulation of both DENV-1 and DENV-2 during the 1994 outbreak [12]. Since then, the disease has emerged as a major public health problem in the country [13]. 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 [16]. Dengue showed a resurgence in November 2010, especially in Sindh, Punjab and Khyber Pakhtunkhwa regions subsequent to massive floods in July same year [6]. The outbreak escalated in Lahore, Punjab, in 2011 with 22,562 cases and 363 deaths due to severe DHF [6].

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 [17]. 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 [18] and on rare occasions, certain mutations lead to phenotypic changes in the viruses that alter their potential to cause outbreaks associated with severe disease [19]. 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 [20]. In Pakistan, sequencing of NS3 gene [21] 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 [17]. 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).
Table 1

Number of dengue positive samples in study cohort during 2006-2011

Collection year

No. of sera analysed

No. of positive samples

DENV Serotype

DENV-1

DENV-2

DENV-3

DENV-4

2006

97

46

0

0

46

0

2007

3

3

0

0

3

0

2008

8

4

0

2

2

0

2009

61

28

0

21

6

1

2011

30

13

0

13

0

0

Unknown

1

0

0

0

0

0

Total

200

94

0

36

57

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 [16]. 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 [22]. 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 [6].

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 [23]. 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.
Figure 1

Phylogenetic tree of DENV-2. The maximum-likelihood tree was constructed based on the complete envelope gene sequences generated during this study and retrieved from GenBank. Sequences from Pakistan are highlighted in red. Figures on branches are bootstrap percentages. Only bootstrap values more than 75% are shown on the major nodes.

Table 2

Genome-wide mutation and selection pressure analysis of Pakistan DENV-2 isolates

Polyprotein position

Gene

Gene position

Wild type

Post-2008 Pakistan isolates (n=5)

Other reports of similar substitutions

Mixed effect model of evolution (MEME) analysis

α

β-

q-

β+

P value

Log L

130

prM

16

R

K

Sri Lanka2003-04§ & Pakistan [17]

1.87

0

0.97

91.05

0.15

−35.52

184

prM

70

S

A

Pakistan £ [17] & Malaysia (sylvatic) [24, 25]

3.61

0

1

134.67

0.05

−41.97

862

NS1

87

T

S

Sri Lanka2003-04§

1.25

0

0.97

6.49

0.74

−27.76

1160

NS2A

33

I

V

India2001-06§ & Sri Lanka 1996-2004§

0

1

0.06

0.33

0.66

−13.61

1243

NS2A

116

I

T

Thailand 1984 [26]

0

0

0.6

2.7

0.42

−40.58

2891

NS5

400

T

A

Sri Lanka 2003-04§

0

0

0.93

2.85

0.48

−10.8

2953

NS5

462

L

I

India2001-06§ & Sri Lanka 1996-04§

1.12

0

0.99

513.74

0.20

−38.62

3107

NS5

616

E

G

India2001-06§ & Sri Lanka 1996-2004§

0

1

0.62

0.85

0.44

−9.39

3122

NS5

631

S

G

India2001-06§ & Sri Lanka 1996-2004§

0.35

0.45

0.16

1.44

0.23

−32.15

3352

NS5

861

K

R

Unique to Pakistan £

0

0

0.98

20.75

0.26

−11.06

3367

NS5

876

N

S

India2001-06§ & Sri Lanka 1996-04§

0

0

0.98

19.81

0.19

−11.44

§Directly obtained from the GenBank database.

£Reported in this study.

Mutations under positive selection and those that are unique to Pakistan DENV-2 isolates are shown in bold.

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) [27] and internal fixed effects likelihood (IFEL) [28] 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 [29]. 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 [30]. 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 [31]. Sequences from another recent study [17] 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 [1]. N-glycosylation of prM protein is believed to play an important role in virion assembly and transport before release [1]. 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 [24] and 1970 [25]. 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 [14]. 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 [16]. 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 [32] whereas the remaining Indian [33] 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 [34], adjacent to one of the 12 strictly conserved Cysteine residues (DENV-3-Cys179) that form disulfide linkages [1]. Disulfide bridges are believed to play an important role in NS1 dimerization [35]. 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.
Figure 2

Phylogenetic tree of DENV-3. The maximum-likelihood tree was constructed based on the complete envelope gene sequences generated during this study and retrieved from GenBank. Sequences from Pakistan are highlighted in red. Figures on branches are bootstrap percentages. Only bootstrap values more than 75% are shown on the major nodes.

Table 3

Genome-wide mutation and selection pressure analysis of Pakistan DENV-3 isolates

Polyprotein position

Gene

Gene position

Wild type

Pakistan isolates (n=6)

Other reports of similar substitutions

Mixed effect model of evolution (MEME) analysis

α

β-

q-

β+

P value

Log L

951

NS1

178

L

M

India 2003–07 [32, 33] & China 2009 §

1.83

0.36

0.99

1844.02

<0.01

−53.95

2063

NS3

590

K

R

Unique to Pakistan £

1.28

0

0.98

42.42

0.58

−17.20

2393

NS4B

151

M

I

India 2003–07 [32, 33] & China 2009§

1.35

0.16

1.00

0.31

1.00

−9.87

2678

NS5

188

T

A

India 2003–07 [32, 33], China 2009§, Venezuela 2001-07§, Puerto Rico 2006§

0.85

0.38

0.99

707.57

0.21

−35.65

3052

NS5

562

Q

L

India 2003–07 [32, 33], China 2009§, Colombia 2004§

1.74

1.00

0

2.65

0.64

−33.19

3385

NS5

895

S

L

Puerto-Rico 2002-03§, Colombia 2003§, Vietnam 2006§

2.00

0

0.99

35.8

0.47

−21.96

§Directly obtained from the GenBank database.

£Reported in this study.

Mutations under positive selection and those that are unique to Pakistan DENV-3 isolates are shown in bold.

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 [36]. 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 [3639]. Domain 3 is also believed to act as the RNA-dependent RNA polymerase binding site [36]. 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 [40]. 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 [41]. 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 [42]. 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.
Figure 3

Phylogenetic tree of DENV-4. The maximum-likelihood tree was constructed based on the complete envelope gene sequences generated during this study and retrieved from GenBank. Sequences from Pakistan are highlighted in red. Figures on branches are bootstrap percentages. Only bootstrap values more than 75% are shown on the major nodes.

Conclusions

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.

Methods

Sample collection

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.

Virus isolation

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 [43]. 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 [43]. 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 [43]. 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 [44] and compared with published sequences of DENV isolates in Genbank database. Phylogenetic analysis of E gene sequences was performed in MEGA5 program [45] 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 [29] 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 [30]. 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 [45], were used and significance levels were set at p<0.05.

Notes

Abbreviations

DENV: 

Dengue virus

DHF: 

Dengue hemorrhagic fever

DSS: 

Dengue shock syndrome

E: 

Envelope

FEL: 

Fixed effects likelihood

GTR: 

General time reversible

HPD: 

Highest probability density

IFEL: 

Internal fixed effects likelihood

MCMC: 

Monte Carlo Markov Chain

MEME: 

mixed effect model of evolution

PCR: 

Polymerase chain reaction

SLAC: 

Single likelihood ancestor counting

TMRCA: 

Time to the most recent common ancestor.

Declarations

Acknowledgements

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’ Affiliations

(1)
Environmental Health Institute, National Environment Agency, [as part of its work as a] WHO Collaborating Centre for Reference and Research of Arbovirus and their Associated Vectors
(2)
Department of Pathology Microbiology, Aga Khan University

References

  1. 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.PubMedView ArticleGoogle Scholar
  2. Gupta B, Reddy BP: Fight against dengue in India: progresses and challenges. Parasitol Res. 2013, 112: 1367-1378. 10.1007/s00436-013-3342-2.PubMedView ArticleGoogle Scholar
  3. 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.PubMedGoogle Scholar
  4. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  5. 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.PubMedGoogle Scholar
  6. 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.PubMedView ArticleGoogle Scholar
  7. Jahan F: Dengue Fever (DF) in Pakistan. Asia Pac Fam Med. 2011, 10: 1-10.1186/1447-056X-10-1.PubMedPubMed CentralView ArticleGoogle Scholar
  8. 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.View ArticleGoogle Scholar
  9. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  10. 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.PubMedView ArticleGoogle Scholar
  11. 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.PubMedView ArticleGoogle Scholar
  12. 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.PubMedView ArticleGoogle Scholar
  13. 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.PubMedView ArticleGoogle Scholar
  14. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  15. 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.PubMedView ArticleGoogle Scholar
  16. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  17. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  18. Steinhauer DA, Holland JJ: Rapid evolution of RNA viruses. Annu Rev Microbiol. 1987, 41: 409-433. 10.1146/annurev.mi.41.100187.002205.PubMedView ArticleGoogle Scholar
  19. 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.PubMedView ArticleGoogle Scholar
  20. 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.PubMedView ArticleGoogle Scholar
  21. 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.PubMedView ArticleGoogle Scholar
  22. 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.PubMedView ArticleGoogle Scholar
  23. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  24. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  25. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  26. 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.PubMedView ArticleGoogle Scholar
  27. 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.PubMedView ArticleGoogle Scholar
  28. 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.PubMedView ArticleGoogle Scholar
  29. Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005, 21: 676-679. 10.1093/bioinformatics/bti079.PubMedView ArticleGoogle Scholar
  30. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  31. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  32. 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-View ArticleGoogle Scholar
  33. 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.View ArticleGoogle Scholar
  34. 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.PubMedView ArticleGoogle Scholar
  35. 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.PubMedView ArticleGoogle Scholar
  36. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  37. 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.PubMedView ArticleGoogle Scholar
  38. 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.PubMedView ArticleGoogle Scholar
  39. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  40. 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.Google Scholar
  41. 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.Google Scholar
  42. 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.Google Scholar
  43. 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.PubMedPubMed CentralGoogle Scholar
  44. 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.PubMedView ArticleGoogle Scholar
  45. 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.PubMedPubMed CentralView ArticleGoogle Scholar
  46. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMedPubMed CentralView ArticleGoogle Scholar
  47. 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.PubMedPubMed CentralView ArticleGoogle Scholar

Copyright

© Koo et al.; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.