Short report | Open | Published:
Molecular characterization of dengue virus reveals regional diversification of serotype 2 in Colombia
Virology Journalvolume 16, Article number: 62 (2019)
Dengue is hyperendemic in Colombia, where a cyclic behavior of serotype replacement leading to periodic epidemics has been observed for decades. This level of endemicity favors accumulation of dengue virus genetic diversity and could be linked to disease outcome. To assess the genetic diversity of dengue virus type 2 in Colombia, we sequenced the envelope gene of 24 virus isolates from acute cases of dengue or severe dengue fever during the period 2013–2016. The phylogenetic analysis revealed the circulation of the Asian-American genotype of dengue virus type 2 in Colombia during that period, the intra-genotype variability leading to divergence in two recently circulating lineages with differential geographic distribution, as well as the presence of nonsynonymous substitutions accompanying their emergence and diversification.
Dengue virus (DENV) is the etiological agent of dengue fever, one of the most important vector-borne viral diseases in terms of morbidity and mortality, according to the World Health Organization (WHO) . In tropical and subtropical regions, there are around 3.6 billion people susceptible to DENV infections. Annually, between 50 and 200 million people are infected worldwide, of which 500,000 progress to severe dengue (SD) and more than 20,000 cases are fatal . After DENV re-emergence in the 1970s and 1980s , Colombia has been considered a hyperendemic country with the presence of the four DENV serotypes, and a cyclic behavior of endemic/epidemic phases with peaks approximately every three to five years . The appearance of severe dengue in Colombia in 1989 coincided with the expansion of the Asian/American genotype of DENV-2 throughout the Americas and the displacement of the American genotype that had been circulating since the early 1970s . During the last two dengue epidemics (2010 and 2013), unprecedented numbers of dengue cases reached 157,152 and 127,219, respectively, followed by interepidemic years in which the number of cases significantly dropped . Intriguingly, the mortality rate of severe dengue cases in Colombia showed a gradual increase since 2007, which was only partially reduced during 2017 and 2018.
The determinants of DENV pathogenesis and disease outcome are multifactorial. The immunologic component as well as the lack of early medical attention have been considered the main factors associated with disease progression and case fatality. However, increasing in vitro, in vivo, and epidemiological evidence also suggests an important role of the viral genetic background in determining the virulence [7,8,9]. The epidemic behavior of the Asian/American genotype contributed to the accumulation of genetic variability conforming several intra-genotype lineages [10, 11], whose importance in explaining virulence differences has been demonstrated [12, 13]. The objective of this study was to determine the genotype and evaluate the genetic diversity and phylogenetic relationship of dengue virus type 2 isolates from patients with dengue and severe dengue in Colombia, during the period 2013–2016.
We performed a retrospective analysis of 1101 archived serum samples from patients with clinical presentations of dengue and severe dengue, collected during the period 2013–2016, according to the mandatory report format of the Program for Dengue Virus Surveillance of the National Institute of Health of Colombia. These samples had been confirmed for dengue infection and serotyped following standard methods as part of the surveillance program. The present study was approved by the Technical and Ethical Committee for Scientific Research (CTIN/CEIN 7–2014 and CTIN/CEIN 23–2014) at the National Institute of Health of Colombia. The final clinical classifications were adjusted according to the Epidemiologic Surveillance System of Colombia – Sivigila, following the WHO recommendations for dengue with warning signs, dengue without warning signs, and severe dengue . Serum samples were diluted 1/100 in Eagle’s Minimum Essential Medium, 200-μl aliquots were used for virus isolation in C6/36 cells and supernatants were collected after nine days post-inoculation or earlier if cytopathic effect was observed. A total of 45 samples were successfully isolated after the first or second passage and the serotype was confirmed by RT-PCR , 24 of which were analyzed in the present study, covering the different geographic regions of the country (Table 1). Nineteen of the selected viral strains were isolated from dengue fever patients, while the other 5 strains were isolated from severe dengue fever patients.
For RNA extraction the QIAamp Viral RNA Mini kit (Qiagen Inc., Chatsworth, CA, USA) was used by following the manufacturer instructions. Amplification of the DENV envelope gene, was performed with the serotype-specific oligonucleotides as described by Domingo et al. , which amplify a 1797 bp fragment. PCR products were purified through the QIAquick PCR purification kit (Qiagen®, Chatsworth, CA, USA) and processed for direct sequencing by using the BigDye® terminator cycle sequencing v3.1 (Applied Biosystems, Carlsbad, CA, USA) and the ABI 3130 Genetic Analyzer (Applied Biosystems, Carlsbad, CA, USA). The electropherograms were visualized, edited and assembled through the SeqMan module of LaserGene® v8.1 (DNASTAR Inc., Madison, WI, USA.).
The sequences obtained in the present study and fifty seven sequences representing the different genotypes of DENV-2 previously deposited in GenBank, mainly those covering the genetic variability within the Asian/American genotype, were aligned and used for phylogenetic reconstruction through Bayesian inference using the MrBayes software , and a total of four MCMCs (three cold, one hot) were evaluated at 1000000 generations with sampling frequency every 100 generations for a total of 10,000 trees. The consensus tree was visualized through FigTree v1.4.3 http://tree.bio.ed.ac.uk/software/figtree/ and was edited in MEGA 7.0 software .
Based on the phylogenetic tree, different well-supported lineages were defined into the Asian/American genotype of DENV-2. Overall mean, intra-lineage and inter-lineage genetic distances were estimated through the MEGA 7.0 software by using the best nucleotide substitution model. The nucleotide and protein alignments showing variable sites and non-synonymous substitutions through the different domains of the envelope protein are depicted (Additional file 2: Figure S2 and Fig. 2, respectively).
All Colombian DENV-2 strains included in the present study circulating during the period 1993–2016, belonged to the Asian/American genotype. Five well-supported intra-genotype lineages with marked spatial and temporal relationships were identified. Two of them (named Lineage 1 and Lineage 2) consisted of sequences from DENV-2 strains recently circulating in Colombia (Fig. 1). Lineages 1 and 2 were represented by sequences of strains circulating during the period 2000–2016 and 1998–2016, respectively. When estimating the global evolutionary divergence at the nucleotide level for the Asian/American genotype in the sequence alignment (using the Tamura-Nei nucleotide substitution model with proportion of invariant sites and gamma distribution with α shape = 2.9), an average of 0.031 substitutions per site was obtained between each pair of sequences of the Asian/American genotype, evidence of high intra-genotype diversity (Additional file 2: Figure S2). The estimated average evolutionary divergences within Asian/American Lineages 1 and 2 were 0.012 and 0.016 substitutions per site, respectively; while the average evolutionary divergence over sequence pairs between Lineages 1 and 2 was 0.031 substitutions per site, revealing the close relationship between strains belonging to each lineage and the marked within-country divergence of the epidemic DENV-2 strains belonging to these two lineages.
Lineage 1 was identified in the departments of Antioquia, Boyacá, Caldas, Cauca, Cundinamarca, Huila, Meta, Putumayo, Quindío, Santander, Valle del Cauca, Nariño and Tolima, that mainly encompass the Andean and Amazon regions in the Southwestern and Central portion of Colombia (Fig. 2a); while Lineage 2 was identified in the departments of Arauca, Casanare, Meta, Norte de Santander, Quindio and Santander, encompassing the Andean and mainly the Orinoquia natural regions in the East and Central portion of the country (Fig. 2a). From the analyzed dataset for the epidemic year 2013, lineages 1 and 2 co-circulated in the departments of Quindío and Santander. In the department of Meta, lineage 1 was identified in 2014 while lineage 2 was identified during 2012–2013. A recent study mainly including strains from the Santander department, allowed the identification of a single recently circulating lineage with a mean estimated time to the most recent common ancestor around 1987 and closely related to other isolates from Venezuela; however, the very low representation of sequences from the Andean and Amazon regions prevented the identification of what is denoted in the present study as Lineage 1 . Lineage 1 was found to be closely related to strains from Venezuela and Peru (Fig. 1), suggesting that DENV-2 circulation in these bordering countries is marked by importation and exportation of strains, and which is supported by the geographical proximity and commercial exchange between these regions. The third lineage included strains from Central America with evidence of introduction to Colombia in 2007, but there was no evidence of dispersion and diversification inside the country. The fourth lineage included strains from Bolivia, Brazil, and Peru during the period 2008–2014, without evidence of circulation in Colombia. The fifth lineage fell in an ancestral position in the phylogenetic tree and was conformed by strains that circulated during the period 1990–2007 in Colombia and other South American countries.
When the limited information related to the clinical classification of patients from the present and previous studies was mapped to the phylogenetic tree, dengue and severe dengue cases were associated with both recently circulating lineages belonging to the Asian/American genotype (Fig. 1). All Colombian sequences obtained in the present study contained the distinctive asparagine amino acid at position 390 of the envelope protein. Twenty-two nonsynonymous substitutions were observed when Colombian sequences of the Asian/American genotype were compared to the earliest Colombian sequence included in the dataset, isolated in 1993 (Fig. 2b). Most nonsynonymous substitutions (63.6%) occurred in the domain III (residues 296–394) which has been reported to directly interact with the cellular proteins during virus entry and constitutes a major target for neutralizing antibodies . An isoleucine to valine amino acid change (I312V) in the envelope protein was found to be exclusively present in one Colombian DENV-2 strain isolated from a severe dengue case in the present study. Further investigation will be needed to establish its role in viral pathogenesis.
Amino acid changes were mapped on the phylogenetic tree, enabling identification of a I462V substitution accompanying the emergence of the lineage 2 and being preserved in all descendants of the monophyletic group during the period 1998–2016. The T359I and V324I substitutions were present in lineage 2, in those strains diversifying during the period 2008–2016, as well as the K310R substitution in a subset of more recent strains belonging to a monophyletic group (2013–2016) (Fig. 1).
In contrast, the emergence of lineage 1 was not characterized by nonsynonymous changes at the envelope protein. Only a few amino acid substitutions appeared in subsets of strains. The I170T change accompanied the evolution of one strain isolated in 2011in Peru and six strains from Colombia covering the period 2013–2014. The independent occurrence of the I170T change within the lineage 1 and lineage 4 was evidence of convergent evolution. The reversion to the ancestral state (T170I and I324V) and the evidence of convergent evolution (K310R, T404I, E360G, I170T and T359I) are mainly due to amino acid changes located at domain III of the envelope protein and could be suggesting positive selection pressure acting at these sites that should be assessed in future studies.
The presence of Colombian isolates of DENV-2 through the whole branching of the highly diversified Asian/American genotype demonstrates the sustained transmission of the virus through time. The geographic and temporal segregation of the different lineages with strains from bordering countries are evidence of an intense dynamics determined by lineage extinction and a bi-directional flow of strains that could explain the drastic changes in the disease epidemiology [10, 12, 21].
The broad clinical spectrum of the disease ranging from asymptomatic to severe and fatal cases represents an opportunity for future clinical and virological studies attempting to demonstrate the existence of a viral genetic contribution to the disease outcome. The increase in the mortality rate of severe dengue cases during the last years (Additional file 1: Figure S1) suggests increased virulence of DENV strains through time. Nevertheless, unsolved difficulties in the clinical management and immunologic factors related to the hyperendemic circulation of the four serotypes or closely-related flaviviruses can also be contributing factors to the disease outcome. Under-reporting of the dengue cases to the National Surveillance System and a very low isolation rate from archived samples constitute limitations. Nevertheless, this study described the recent circulation of lineages 1 and 2 of the Asian/American genotype of DENV-2 in Colombia, the microevolution and differential geographic distribution at the national level.
Notwithstanding the growing epidemiologic and experimental data of the presence of determinants of virulence in the DENV genome [8, 22,23,24], further comparative analysis of full-length viral genomes and functional studies on the role of specific substitutions will be decisive for advancing the elucidation of its epidemiology and disease dynamics.
Reverse transcription-polymerase chain reaction
World Health Organization
World Health Organization. Global strategy for dengue prevention and control, 2012–2020. France: WHO Press; 2012. p. 43.
Murray MN, Quam MB, Wilder-Smith A. Epidemiology of dengue: past, present and future prospects. Clin Epidemiol. 2013;5:299–309.
Boshell J, Groot H, Gacharna MG, Márquez G, González M, Gaitán MO, et al. Dengue en Colombia. Biomédica. 1986;6:101–6.
Padilla JC, Rojas DP, Sáenz-Gómez R. Dengue en Colombia: Epidemiología de la reemergencia a la hiperendemia. 1st ed. Bogotá DC: Guías de Impresión Ltda; 2012. p. 246.
Méndez JA, Usme-Ciro JA, Domingo C, Rey GJ, Sánchez JA, Tenorio A, et al. Phylogenetic reconstruction of dengue virus type 2 in Colombia. Virol J. 2012;99(11):64.
Pan American Health Organization/World Health Organization. Dengue: Datos, mapas y estadísticas Washington: Pan American Health Organization; [Available from: https://www.paho.org/hq/index.php?option=com_topics&view=rdmore&cid=3274&Itemid=40734&lang=es.
Kyle JL, Harris E. Global spread and persistence of dengue. Annu Rev Microbiol. 2008;62:71–92.
Rico-Hesse R, Harrison LM, Salas RA, Tovar D, Nisalak A, Ramos C, et al. Origins of dengue type 2 viruses associated with increased pathogenicity in the Americas. Virology. 1997;230(2):244–51.
Rico-Hesse R. Dengue virus virulence and transmission determinants. Curr Top Microbiol Immunol. 2010;338:45–55.
Bennett SN, Holmes EC, Chirivella M, Rodriguez DM, Beltran M, Vorndam V, et al. Molecular evolution of dengue 2 virus in Puerto Rico: positive selection in the viral envelope accompanies clade reintroduction. J Gen Virol. 2006;87(Pt 4):885–93.
Drumond BP, Mondini A, Schmidt DJ, de Morais Bronzoni RV, Bosch I, Nogueira ML. Circulation of different lineages of dengue virus 2, genotype American/Asian in Brazil: dynamics and molecular and phylogenetic characterization. PLoS One. 2013;8(3):e59422.
OhAinle M, Balmaseda A, Macalalad AR, Tellez Y, Zody MC, Saborio S, et al. Dynamics of dengue disease severity determined by the interplay between viral genetics and serotype-specific immunity. Sci Transl Med. 2011;3(114):114ra28.
Williams M, Mayer SV, Johnson WL, Chen R, Volkova E, Vilcarromero S, et al. Lineage II of southeast Asian/American DENV-2 is associated with a severe dengue outbreak in the Peruvian Amazon. Am J Trop Med Hyg. 2014;91(3):611–20.
World Health Organization. Dengue: guidelines for diagnosis, treatment, prevention and control. New edition. Geneva: WHO Press; 2009. p. 10–2.
Usme-Ciro JA, Gómez-Castañeda AM. Gallego-Gómez JC. [Molecular detection and typing of dengue virus by RT-PCR and nested PCR using degenerated oligonucleotides. Salud Uninorte. 2012;28:1–15.
Domingo C, Niedrig M, Teichmann A, Kaiser M, Rumer L, Jarman RG, et al. 2nd international external quality control assessment for the molecular diagnosis of dengue infections. PLoS Negl Trop Dis. 2010;4(10):e833.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Jiménez-Silva CL, Carreño MF, Ortiz-Baez AS, Rey LA, Villabona-Arenas CJ, Ocazionez RE. Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia. PLoS One. 2018;13(8):e0203090.
Crill WD, Roehrig JT. Monoclonal antibodies that bind to domain III of dengue virus E glycoprotein are the most efficient blockers of virus adsorption to Vero cells. J Virol. 2001;75(16):7769–73.
Runge-Ranzinger S, Horstick O, Marx M, Kroeger A. What does dengue disease surveillance contribute to predicting and detecting outbreaks and describing trends? Tropical medicine & international health. Tropical Med Int Health. 2008;13(8):1022–41.
Cologna R, Armstrong PM, Rico-Hesse R. Selection for virulent dengue viruses occurs in humans and mosquitoes. J Virol. 2005;79(2):853–9.
Vaughn DW, Green S, Kalayanarooj S, Innis BL, Nimmannitya S, Suntayakorn S, et al. Dengue viremia titer, antibody response pattern, and virus serotype correlate with disease severity. J Infect Dis. 2000;181(1):2–9.
Cologna R, Rico-Hesse R. American genotype structures decrease dengue virus output from human monocytes and dendritic cells. J Virol. 2003;77(7):3929–38.
The authors thank the National Laboratory Network for routine virologic surveillance of dengue in Colombia, the Sequencing and Genomics Unit team at the Virology laboratory INS, and Dr. Yves Girerd-Chambaz from Sanofi Pasteur for providing the Colombian states and cities where CYD15 samples included in the phylogenetic analysis were collected.
This work was supported by COLCIENCIAS (Convocatoria 617–2013 and 210465740977) and the Colombian National Institute of Health (CTIN codes: 07–2014 and 23–2014).
Availability of data and materials
Envelope sequences generated in this study have been deposited in GenBank under the accession numbers KU878564 to KU878575, KY905139, KY905140, MK016290, MK016291 and MK016293 to MK016300.
Ethics approval and consent to participate
The study was approved by the Technical and Ethical Committee for Scientific Research (CTIN/CEIN 7–2014 and CTIN/CEIN 23–2014) at the National Institute of Health of Colombia. All samples analyzed were anonymized and used only for public health surveillance purposes.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Incidence of severe dengue and mortality rate in Colombia during the period 2007–2018. The Mortality rate of severe dengue cases was estimated as the number of fatal cases per hundred severe dengue cases. * Epidemiological Week 37. (PNG 50 kb)
Figure S2. Variation in the nucleotide sequence of the envelope gene of Colombian strains of DENV-2. Representative sequences of strains that have circulated in Colombia and those obtained in the present study were aligned and variable sites along the nucleotide sequence compared to the Asian/American genotype representative strain Jamaica N.1409. D: Dengue; SD: Severe dengue. (PDF 357 kb)