Phylogenetic studies reveal existence of multiple lineages of a single genotype of DENV-1 (genotype III) in India during 1956–2007
- Himani Kukreti1,
- Paban Kumar Dash2,
- Manmohan Parida2,
- Artee Chaudhary1,
- Parag Saxena2,
- RS Rautela1,
- Veena Mittal3,
- Mala Chhabra3,
- D Bhattacharya3,
- Shiv Lal1, 3,
- PV Lakshmana Rao2 and
- Arvind Rai1Email author
© Kukreti et al. 2009
Received: 11 August 2008
Accepted: 06 January 2009
Published: 06 January 2009
Dengue virus type 1 (DENV-1) have been mostly circulating silently with dominant serotypes DENV-2 and DENV-3 in India. However recent times have marked an increase in DENV-1 circulation in yearly outbreaks. Many studies have not been carried out on this virus type, leaving a lacunae pertaining to the circulating genotypes, since its earliest report in India. In the present study, we sequenced CprM gene junction of 13 DENV-1 isolated from Delhi and Gwalior (North India) between 2001–2007 and one 1956 Vellore isolate as reference. For comparison, we retrieved 11 other Indian and 70 global reference sequences from NCBI database, making sure that Indian and global isolates from all decades are available for comparative analysis.
The region was found to be AT rich with no insertion or deletion. Majority of the nucleotide substitutions were silent, except 3 non-conservative amino acid changes (I → T, A → T and L → S at amino acid positions 59,114 and 155 respectively) in the Indian DENV-1 sequences, sequenced in this study. Except two 1997–98 Delhi isolates, which group in genotype I; all other Indian isolates group in genotype III. All Indian genotype III DENV-1 exhibited diversity among them, giving rise to at least 4 distinct lineages (India 1–4) showing proximity to isolates from diverse geographic locations.
The extensive phylogenetic analysis revealed consistent existence of multiple lineages of DENV-1 genotype III during the last 5 decades in India.
Dengue fever (DF) is one of the most important arboviral diseases of humans in tropic and sub-tropics [1, 2]. In South-East Asia, with a total population of 1.5 billion, approximately 1.3 billion people live at risk of acquiring DF or DHF [3, 4]. Its etiological agent, Dengue virus belongs to family Flaviviridae, genus Flavivirus; and exists in 4 antigenically distinct serotypes, Dengue virus type 1–4 (DENV-1 to 4) . Although history of dengue virus in India dates back to 1946, however the first major outbreak was reported in 1963 in Calcutta. Since then many outbreaks have been reported from all over the country [6–9]. Although all 4 dengue serotypes have been reported to circulate in the country  but only DENV-2 and DENV-3 have been implicated in major DF/dengue hemorrhagic fever (DHF) outbreaks [11–14]. We and some other workers have earlier reported CprM gene based genotyping of DENV-2 and DENV-3 which has proved to be useful for carrying out molecular epidemiology of these viruses [12, 14–16]. Lately there has been a rise in DENV-1 associated cases, which account for around 30% of the total cases in most recent 2006 DF outbreak, when it co-existed with pre-dominant DENV-3 . Earlier DENV-1 has been reported from south India (Vellore) in 1956 and 1962–64; and from north India (Delhi and Gwalior) in 1970, 1982, 1997–98 and 2002–2006. Inspite of all these reports, no attempt has been made to study the phylogeny of this virus since its emergence. With increase in DENV-1 cases, the need to understand the genetic nature of circulating DENV-1 has become even more necessary as it will harbors significant information regarding the genotypes of the DENV-1 circulating in the country for so long. In addition it will also indicate whether these viruses were pre-existing or an importation. With this view, present study was undertaken to understand the genetic nature of circulating DENV-1 in India and to trace their evolution during the last 50 years, by sequencing CprM gene junction of 13 DEN-1 isolated in India during 2001–2007 and comparing 354 bp of this region with 11 other Indian DENV-1 sequences reported till date with at least one sequence from each decade, since it was first reported in the country. Seventy other global reference sequences were also retrieved from NCBI nucleotide database for comparison. This study shall help to fill the lacunae pertaining to knowledge of circulating DENV-1 genotypes in India since its first report.
On comparison of sequences, we observed that one 2002 Gwalior and two 2001 Delhi CprM sequences exhibited close nucleotide identity of 99.7% to a 2004 Reunion Island sequence. It was also found that the most recent Indian isolates from 2006 and 2007 show mean sequence identities of 96.45%, 95.35%, 97.5% and 96.75% with 1956 Vellore, 1970 Delhi, 1982 Delhi and 2004 Gwalior sequences; and 98.6–99.2% (mean 98.9%) nucleotide identity among themselves. On analyzing the Indian reference sequences taken for comparison, it was observed that 2005 Delhi sequences exhibited 95.8–99.2% (mean 97.5%) identity with sequences from South American isolates sampled between 1989–2001. Three sequences from south India (Vellore) sampled in 1962–64 show a nucleotide identity of 96.6–96.9% (mean 96.75%) with a 1956 sequence from the same area. Two Delhi sequences from years 1997 and 1998 show more nucleotide identity i.e. 95.8–96.6% (mean 96.2%) with Japan Mochizuki sequence isolated in 1943; than to the other Indian sequences, with whom theses two sequences show a nucleotide identity of only 93%.
Genotype I in which two of the Indian isolates from Delhi (1997–98) cluster; consists of viruses from Asia except 1 African isolate from Djibouti isolated in 1998. Although, most Asian strains in this genotype are from South East Asian countries viz. Cambodia, Vietnam, Thailand and Singapore; 2 recent isolates from China and 1 very old strain from Japan (Japan 43 Mochizuki strain) are the only 3 strains from North Asia that lie in this genotype.
None of the Indian strains was classified as genotype II. This genotype was more cosmopolitan in nature with strains from East Africa, Asia and South Pacific. East African strains include recent viruses from East African islands of Seychelles and Reunion isolated during 2003–04; South Pacific strains include viruses from Hawaii, Nauru Island and Micronesia isolated between 74-04; and Asian strains include viruses from China, Philippines and Indonesia isolated during 1995–2004. Although African strains are clustering together, making a distinct grouping; strains from South Pacific and Asia do not show any such geographical clustering within this genotype.
The large scale climatic and demographic changes have resulted in distribution of Aedes species in hitherto unknown areas. This has led to appearance of dengue infection in many newer areas with an estimated 3 billion people living at risk of this infection around the world. Considering the high rates of Aedes infestation and the presence of high-risk susceptible population, the probability of dengue epidemic always remain very high in India . The pathogenesis of dengue is not yet fully understood. Though re-infection with a heterologous dengue serotype remain a major risk factor for DHF, however, occurrence of DHF even in primary infection makes the role of virus in pathogenesis all the more important. The incursion of new genotype into an area is also being attributed to the severe form of the disease . All these facts drive the scientific community to pay more attention towards the genetic nature of dengue viruses and their spread in the population.
Like other RNA viruses, the dengue viruses also revealed strong genetic diversity. Dengue serotypes are further classified into 4–5 'genotypes' based on their genetic diversity . However in recent reports many researchers have categorized DENV-1 in three distinct genotypes [21–23]. Researchers have also reported clustering of sequences below the genotype level that correlate with the geographical origin and/or time of isolation and defined as lineages/clades [21, 23, 24]. Different regions of dengue genome like Envelope, E-NS1, C-prM and complete genome has been utilized for the genotyping [15, 20, 21, 24]. We have utilized the sequence information of CprM gene junction in this study. The genotyping based on the CprM gene junction has been adopted by several researchers including us in the recent past [12, 14, 15, 25]. This results in faster and economical genotyping due to utilization of a single set of primer pair for both amplification and sequencing .
The involvement of DENV-1 in major dengue outbreaks in India was not recorded. Most of the major dengue outbreaks are attributed to DENV-2 and DENV-3 viruses. However, DENV-1 were isolated from different parts of India at regular intervals, since its first isolation from southern India (Vellore) in 1956. It is increasingly implicated as a minor serotype during recent outbreaks in India, including the Delhi outbreak in 2006 [17, 27]. The lack of sequence analysis of Indian DENV-1 also affects the effective monitoring of circulating genotypes in India. In this study, we determined the nucleotide sequence of 13 recent DENV-1 directly from clinical samples to avoid selection and sequencing of mutant viral RNA which may occur during isolation/propagation in mouse brain/cell culture . The comparison of deduced amino acid sequence revealed 3 important unique substitutions in 2001–02 and some of the 2006–07 Indian isolates. These involve substitutions of hydrophobic Isoleucine and Leucine by hydrophilic Threonine and Serine respectively in some 2006–07 isolates; and hydrophobic Alanine by hydrophilic Threonine in Delhi 2001 and Gwalior 2002 isolates. The hydrophilic amino acids are found mostly on the surface of the protein and are involved in the immunological interactions. The implications of these unique changes in recent Indian DENV-1 warrants further studies to understand their virulence and epidemic potential.
The phylogenetic analysis revealed the Indian DENV-1 collected over last 5 decades fall into two genotypes (I and III). However, majority of the viruses are grouped within genotype III, two isolates sampled from Delhi during 1997–98 group in genotype I. These genotype I viruses from Delhi were found closely related to the 1943 Japanese Mochizuki strain. Their proximity was however difficult to explain. Similar type of genotype I viruses were also found circulating predominantly in other Asian countries including Thailand and China around the same time. Presence of an African isolate from Djibouti in this group probably signifies the introduction of these Asian viruses in Djibouti or vice-versa. Thailand strains belonging to this genotype were isolated over four decades (1964–2001); indicating persistence of this genotype in Thailand for a long time. However, the restricted circulation of genotype I viruses in Delhi for just 2 years is very surprising.
The phylogenetic analysis clearly revealed continuous circulation of genotype III viruses over the last 50 years in India. Interestingly, majority of the Indian genotype III viruses are found to be phylogenetically quite distinct. It is surprising to find that the recent viruses recovered in quick succession in 2001–02, 2004 and 2006–07 from northern India also belong to separate lineages. In contrast, this type of distinct lineage pattern was not observed among Indian DENV-2 and 3 [12, 14].
Domingo et al.,  has earlier established 2 distinct lineages of DENV-1 in India and designated them as India-1 and India-2. While India-1 has shown closeness to South American strains, India-2 exhibited close proximity to a 1993 isolate from Singapore. In our earlier study we have substantiated the presence of India-1 and India-2 lineages during 2005 and 2006 respectively and also reported a third lineage that was present in south India (Vellore) between 1962–64, whose present existence could not be established . In this study we now report another lineage, which consists of viruses from Delhi and Gwalior sampled in 2001 and 2002 respectively; that we designated as India-4. During the same period presence of India-2 lineage was also reported , there by suggesting the co-circulation of India-2 and India-4 lineage in India during 2001. As evident from the phylogenetic tree (Figure 3) Indian DENV-1 isolated in 1956, 1970, 1982 and 2004 could not be designated in a separate lineage but their relatedness could be deciphered by clustering together of these isolates in the phylogenetic tree.
In the cosmopolitan genotype (genotype II) with viruses from East Africa, Asia and South Pacific, only African isolates clustered together indicating independent evolution of these viruses in this area.
Thus we conclude that during the last 5 decades there has been a persistence of genotype III of DENV-1 in India with genotype I being present only during 2 years i.e. 1997–98 in Delhi. Most of the Indian DENV-1, grouping in genotype III were quite diverse, giving rise to distinct lineages. Thus, there has been persistence of multiple lineages of DENV-1 genotype III during the last 50 years in India.
This study clearly identified the genotypes of DENV-1 circulating in India since 1956. It also confirmed the utility of CprM gene junction for rapid and economical genotyping in endemic areas. The prevalence of multiple lineages of DENV-1 in India warrants sustained monitoring of the circulating viruses, to implement effective control measures at the earliest.
Materials and methods
Serum samples from febrile patients suspected for dengue infection were collected from Delhi and Gwalior during 2001–2007. Informed consent from all the patients and/or their parents (in minors) was obtained, before collection of clinical samples. Approval of the ethical committee of both institutions was obtained to carry out the present study. Thirteen serum samples found positive for DENV-1 RNA were included in this study.
Dengue virus serotype 1 (P-23086) obtained from the National Institute of Virology (NIV), Pune, India was used as reference strain in this study. This is the prototype Indian DENV-1 virus isolated from Vellore, India in 1956.
Extraction of viral RNA
Viral RNA was extracted from 140 μl of serum samples and P-23086 infected C6/36 supernatant by using QIAamp viral RNA mini kit (Qiagen, Germany) in accordance with the manufacturer's instructions. Finally RNA was eluted in 50 μl of nuclease free water and stored at -80°C until use.
Reverse transcription- Polymerase chain reaction (RT-PCR)
The RT-PCR was carried out in a 25 μl reaction volume using the access quick one-step RT-PCR kit (Promega, USA) containing PCR master mix, AMV-RT, and respective sense and antisense primers [(D1: 5'-TCAATATGCTGAAACGCGCGAGAAACCG-3') (D2: 5'-TTGCACCAACAGTCAATGTCTTCAGGTTC-3')]  in a thermal cycler (BioRad, USA). The thermal profile of the RT-PCR reaction was- RT step at 42°C for 45 min, followed by PCR step of initial denaturation at 95°C for 2 min, followed by 35 cycles of denaturation at 95°C for 1 min, annealing at 55°C for 1 min, extension at 72°C for 2 min and final extension at 72°C for 10 min. The PCR products were gel purified from 1.2% agarose gel using the QIAquick PCR purification kit (Qiagen, Germany) and used as template in sequencing reaction.
Both strands of the purified amplicons were sequenced employing Big dye terminator cycle sequencing ready reaction kit (Applied Biosystems, USA) following the standard protocol described by us earlier . The cycle sequenced product was purified by precipitation with 75% isopropanol and was vaccum dried. The DNA pellet was resuspended in 10 μl of template suppression reagent (TSR), heated at 95°C for 2 min and loaded on the ABI 310 automated DNA sequencer (Applied Biosystems, USA).
Sequence alignment and phylogenetic analysis
Details of all the Indian DENV-1 isolates included in this study.
Year of Isolation
India 56 vellore
India 62 623996/vellore
India 63 631286/1963/vellore
India 64 64411/1964/vellore
India 70 703180/1970/delhi
India 82 826891/1982/delhi
India 97 1021/1997/delhi
India 98 1412/1998/delhi
India 01 D1/1CprM/Del01
India 01 D1/2CprM/Del01
India 02 GWL14
India 04 GWL19
India 05 del/05/1474/D1
India 05 del/05/902/D1
India 05 AARS01
India 05 AARS02-2005
India 06 07/1/del2006
India 06 08/1/del2006
India 06 09/1/del2006
India 06 10/1/del2006
India 06 11/1/del2006
India 06 12/1/del2006
India 06 13/1/del2006
India 07 D1/1CprM/Del07
India 07 D1/2CprM/Del07
The authors are thankful to the Director, DRDE, Gwalior and Director, NICD, Delhi for their keen interest and constant support in this study. The authors are also thankful to the Superintendent, JA Hospital, Gwalior and Dr. Ashok Rawat, Deputy Health Officer, Municipal Corporation of Delhi, for providing the clinical samples. HK also acknowledges partial financial support form Council of Scientific and Industrial Research (CSIR), Delhi and WHO country office, Delhi, India.
- Guzmán MG, Kourí G: Dengue: an update. Lancet Infect Dis 2002, 2:33–42.View ArticlePubMed
- Gubler DJ: Dengue/dengue haemorrhagic fever: history and current status. Novartis Found Symp 2006, 277:3–16.View ArticlePubMed
- Setiati TE, Wagenaar JFP, de Kruif MD, Mairuhu ATA, Gorp ECM, Soemantri A: Changing epidemiology of changing haemorrhagic fever in Indonesia. Dengue Bull 2006, 30:1–14.
- World Health Organization: Dengue and dengue haemorrhagic fever. 2002.
- Lindenbach BD, Rice CM: Molecular biology of flaviviruses. Adv Virus Res 2003, 59:23–61.View ArticlePubMed
- Ram S, Khurana S, Kaushal V, Gupta R, Khurana SB: Incidence of dengue fever in relation to climatic factors in Ludhiana, Punjab. Indian J Med Res 1998, 108:128–33.PubMed
- Dar L, Broor S, Sengupta S, Xess I, Seth P: The first major outbreak of dengue haemorrhagic fever in Delhi, India. Emerg Infect Dis 1999, 5:589–90.View ArticlePubMed
- Joshi PT, Pandya AP, Anjan JK: Epidemiological and entomological investigation in dengue outbreak area of Ahmedabad district. J Commun Dis 2000, 32:22–7.PubMed
- Gore MM: Need for constant monitoring of dengue infections. Indian J Med Res 2005,121(1):9–12.PubMed
- Dar L, Gupta E, Narang P, Broor S: Cocirculation of dengue serotypes, Delhi, India, 2003. Emerg Infect Dis 2006, 12:352–3.PubMed
- Singh UB, Maitra A, Broor S, Rai A, Pasha ST, Seth P: Partial nucleotide sequencing and molecular evolution of epidemic causing dengue 2 strains. J Infect Dis 1999, 180:959–65.View ArticlePubMed
- Dash PK, Parida MM, Saxena P, Kumar M, Rai A, Pasha ST, Jana AM: Emergence and continued circulation of Dengue-2 (Genotype IV) virus strains in northern India. J Med Virol 2004, 74:314–22.View ArticlePubMed
- Dash PK, Saxena P, Abhyankar A, Bhargava R, Jana AM: Emergence of dengue virus type-3 in northern India. Southeast Asian J Trop Med Public Health 2005,36(2):370–377.PubMed
- Dash PK, Parida MM, Saxena P, Abhyankar A, Singh CP, Tewari KN, Jana AM, Sekhar K, Rao PVL: Reemergence of dengue virus type-3 (subtype-III) in India: implications for increased incidence of DHF & DSS. Virol J 2006, 3:55.View ArticlePubMed
- Avilés G, Rowe J, Meissner J, Manzur Caffarena JC, Enria D, St Jeor S: Phylogenetic relationships of dengue-1 viruses from Argentina and Paraguay. Arch Virol 2002, 147:2075–87.View ArticlePubMed
- Singh UB, Seth P: Use of nucleotide sequencing of the genomic cDNA fragments of the capsid/premembrane junction region for molecular epidemiology of dengue type 2 viruses. Southeast Asian J Trop Med Public Health 2001, 32:326–35.PubMed
- Kukreti H, Chaudhary A, Rautela RS, Anand R, Mittal V, Chhabra M, Bhattacharya D, Lal S, Rai A: Emergence of an independent lineage of dengue virus type 1 (DENV-1) and its co-circulation with predominant DENV-3 during the 2006 dengue fever outbreak in Delhi. Int J Infect Dis 2008,12(5):542–549.View ArticlePubMed
- Pandya G: Prevalence of dengue infections in India. Def Sci J 1982, 32:359–70.
- Rico-Hesse R, Harrison LM, Salas RA, Tovar D, Nisalak A, Ramos C, Boshell J, de Mesa MT, Nogueira RMR, da Rosa AT: Origins of dengue type-2 viruses associated with increased pathogenicity in the Americas. Virology 1997,230(2):244–251.View ArticlePubMed
- Rico-Hesse R: Molecular evolution and distribution of dengue viruses type 1 and 2 in nature. Virology 1990,174(2):479–493.View ArticlePubMed
- Zhang C, Mammen PM, 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:15123–30.View ArticlePubMed
- A-Nuegoonpipat A, Berlioz-Arthaud A, Chow V, Endy T, Lowry K, Mai le 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.PubMed
- Thu HM, Lowry K, Jian L, Hlaing T, Holmes EC, Aaskov J: Lineage extinction and replacement in dengue type 1 virus populations is due to stochastic events rather than natural selection. Virology 2005,336(2):163–172.View Article
- Domingo C, Palacios G, Jabado O, Reyes N, Niedrig M, Gasco'n J: Use of a short fragment of the C-terminal E gene for detection and characterization of two new lineages of dengue virus 1 in India. J Clin Microbiol 2006, 44:1519–29.View ArticlePubMed
- Kumar M, Pasha ST, Mittal V, Rawat DS, Arya SC, Agarwal N, Bhattacharya D, Lal S, Rai A: Unusual emergence of Guate98-like molecular subtype of DEN-3 during 2003 dengue outbreak in Delhi. Dengue Bull 2004, 28:161–7.
- 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–51.PubMed
- Bharaj P, Chahar HS, Pandey A, Diddi K, Dar L, Guleria R, Kabra SK, Broor S: Concurrent infections by all four dengue virus serotypes during an outbreak of dengue in 2006 in Delhi, India. Virology J 2008, 5:1.View Article
- Leitmeyer KC, Vaughn DW, Watts DM, Salas R, Chacon de IV, Ramos C, Rico-Hesse R: Dengue virus structural differences that correlate with pathogenesis. J Virol 1999, 73:4738–4747.PubMed
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform 2004, 5:150–63.View ArticlePubMed
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987, 4:406–25.PubMed
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.