# Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia

Cinthy Jiménez-Silva, María Fernanda Carreño, Ayda Susana Ortiz-Baez, Luz

Aida Rey, Christian Julián Villabona-Arenas, Raquel Ocazionez

► To cite this version:

Cinthy Jiménez-Silva, María Fernanda Carreño, Ayda Susana Ortiz-Baez, Luz Aida Rey, Christian Julián Villabona-Arenas, et al.. Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia. PLoS ONE, 2018, 13 (8), pp.e0203090. 10.1371/journal.pone.0203090. lirmm-03113139

HAL Id: lirmm-03113139

https://hal-lirmm.ccsd.cnrs.fr/lirmm-03113139

Submitted on 18 Jan 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia

Cinthy L. Jiménez-Silva¹, María Fernanda Carreño¹, Ayda Susana Ortiz-Baez¹, Luz Aida Rey², Christian Julián Villabona-Arenas³,⁴*, Raquel E. Ocazionez¹

1. Laboratorio de Arbovirus, Centro de Investigaciones en Enfermedades Tropicales (CINTROP), Universidad Industrial de Santander, Bucaramanga, Colombia,
2. Unidad de Epidemiología Clínica, Escuela de Medicina, Universidad Industrial de Santander, Bucaramanga, Colombia,
3. Laboratoire TransVIHMI, Institut de Recherche pour le Développement (IRD), Université de Montpellier, Montpellier, France,
4. Laboratoire d'Informatique, de Robotique et de Microélectronique de Montpellier (LIRMM), Institut de Biologie Computationnelle (IBC), Université de Montpellier, Montpellier, France

* christian-julian.villabona-arenas@ird.fr

<figure><img src="image_2.png" /></figure>

## Abstract

Dengue is a prevalent disease in Colombia and all dengue virus serotypes (DENV-1 to -4) co-circulate in the country since 2001. However, the relative impact of gene flow and local diversification on epidemic dynamics is unknown due to heterogeneous sampling and lack of sufficient genetic data. The region of Santander is one of the areas with the highest incidence of dengue in Colombia. To provide a better understanding of the epidemiology of dengue, we inferred DENV population dynamics using samples collected between 1998 and 2015. We used Bayesian phylogenetic analysis and included 143 new envelope gene sequences from Colombia, mainly from the region of Santander, and 235 published sequences from representative countries in the Americas. We documented one single genotype for each serotype but multiple introductions. Whereas the majority of DENV-1, DENV-2, and DENV-4 strains fell into one single lineage, DENV-3 strains fell into two distinct lineages that co-circulated. The inferred times to the most recent common ancestors for the most recent clades of DENV-1, DENV-2, and DENV-4 fell between 1977 and 1987, and for DENV-3 was around 1995. Demographic reconstructions suggested a gradual increase in viral diversity over time. A phylogeographical analysis underscored that Colombia mainly receives viral lineages and a significant diffusion route between Colombia and Venezuela. Our findings contribute to a better understanding of the viral diversity and dengue epidemiology in Colombia.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Jiménez-Silva CL, Carreño MF, Ortiz-Baez AS, Rey LA, Villabona-Arenas CJ, Ocazionez RE (2018) Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia. PLoS ONE 13(8): e0203090.
https://doi.org/10.1371/journal.pone.0203090

**Editor:** Humberto Lanz-Mendoza, Instituto Nacional de Salud Pública, MEXICO

**Received:** July 1, 2018

**Accepted:** August 14, 2018

**Published:** August 29, 2018

**Copyright:** © 2018 Jiménez-Silva et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This study was supported by the Colombian Government: Science, Technology and Innovation Fund of Sistema General de Regalias, Santander, Casanare, Valle del Cauca (BPIN 2013000100011).

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Dengue disease is highly prevalent in tropical countries due to climate, population growth, unplanned rapid urbanization and increased travel and trade [1,2]. Consequently, the global burden of dengue disease is high: a total of 58.4 million (23.6 million–121.9 million) apparent
cases in 2013 and 1.14 million (0.73 million–1.98 million) disability-adjusted life-years [3]. Dengue is caused by four closely related viruses referred to as serotypes (DENV-1, DENV-2, DENV-3, and DENV-4), which are further subdivided into genotypes [4]. A higher risk of severe dengue is attributed to the co-circulation of multiple serotypes due to antibody-dependent enhancement of infection [2] and to some particular strains [5]. Therefore, documenting serotype prevalence and the dynamics of genetic variants help forecast epidemic impact, ensuring epidemic management and preparedness.

Colombia is a tropical country in the northwest of South America with high dengue incidence. For instance, a total of 752,429 cases were reported during 1980–2007 [6–9]: the country-wide 2001/2002, 2010, and 2013/2014 epidemics had 143,820, 147,257, and 247,075 notified dengue cases, respectively. DENV-2, DENV-1, and DENV-4 were reported for the first time in the country in 1971, 1977, and 1982, respectively [10]. DENV-3 was reported only once in 1975 and re-emerged again in 2001 [10,11]; the four DENV serotypes have been frequently documented in Colombia from that year onwards [9,11,12].

A total of six Colombian states (out of 32) accumulate more than half of the cases in the country: Valle del Cauca (19%), Santander (12%), Antioquia (8%), Norte de Santander (7%), Huila (6%), and Tolima (6%) [7,9,12]. The region of Santander is located in the northeast of Colombia and comprises the states of Santander and Norte de Santander. Previous studies reported changes in the predominance of serotypes over time in this region [11,13–15]. The argument over whether there is a genetic basis or there are random fluctuations that explain the temporal distribution of serotypes has not been settled yet.

Previous studies [16–19] that focused on the molecular epidemiology of dengue in Colombia included sequence data collected up to 2008 and did not include DENV-4. Consequently, many aspects of the local epidemiology remain unclear. In the present study, we obtained 143 new envelope gene sequences from serum samples collected in the region of Santander and six sequences from serum samples collected in other regions between 1998 and 2015. We provide insights into the evolution and population dynamics of all four dengue serotypes in the region of Santander, which help understanding dengue epidemiology in Colombia and might be relevant for future control programs, including vaccination.

# Materials and methods

## Study sites

Serum samples were collected from the metropolitan area of Bucaramanga in the state of Santander and the province of Ocaña in the state of Norte de Santander; these constitute the region of Santander that is located in the central northern part of Colombia [Fig 1]. Bucaramanga is the capital city of Santander and together with three nearby municipalities constitutes the seventh largest metropolitan area of the country (around 1.4 million inhabitants). It has an average altitude of ~1,200 meters above the sea, an annual mean temperature of 25°C and an average amount of annual precipitation of 1,041 mm. The city's urban mass transportation system moves over 100,000 passengers daily and extends to the Metropolitan area [20]. Moreover, Bucaramanga is a transportation hub for the northeast of Colombia and it has a bus terminal where 1.4 million passengers commute to diverse regions of the country on a yearly basis [21]. Bucaramanga's airport has national and international traffic operations, which mobilized more than 1.2 million people in 2010. At least 94% of dengue cases from Santander were reported in the metropolitan area of Bucaramanga [7,9], and it was one of the participant cities in the dengue vaccine CYD-TDV clinical trial [22].

The province of Ocaña is a conurbation of 10 municipalities in an area of 2,065 meters above the sea in the state of Norte de Santander. Ocaña is its capital and its third largest
<figure><img src="image_4.png" /><figcaption>**Fig 1. Geographical origin of the samples from Colombia and dengue incidence.** The region of Santander (Santander and Norte de Santander states) is colored in red.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g001

municipality (around 98,229 inhabitants in 2014). This city has an annual mean temperature of 22°C and an average amount of annual precipitation of 1,000 mm. Cucuta is the capital city of Norte de Santander, it is bordered by Venezuela which makes it an important commercial city. Ocaña is located at a distance of 197 km (122 mi) from Cucuta. Ocaña has a small airport with only regular flights to Cucuta [23].

## Viral strains

DENV strains were obtained from the collection of the Laboratorio de Arbovirus, Centro de Investigaciones en Enfermedades Tropicales (CINTROP), Universidad Industrial de Santander, Bucaramanga. Viruses were isolated by culturing in C6/36 mosquito cells from patient sera collected in previous studies [11,14,15,24,25]. Sera were collected either for routine dengue laboratory diagnosis at medical institutions—for which an Institutional Ethics Committee approval was not required—or collected from patients enrolled in cross-sectional clinical trials, which were approved by the Research Ethics Committee of the Universidad Industrial de Santander. In the latter, an informed consent from each patient was obtained. All virus samples were analyzed anonymously. We selected samples to represent years in which each virus serotype was recorded between 1998 and 2015: 137 samples from the region of Santander (60 DENV-1, 33 DENV-2, 39 DENV-3, and 11 DENV-4) and six samples that were available from the Valle del Cauca, Bolivar, and Cesar states (southwestern and northeastern regions).

## Full-length envelope gene (E-gene) sequencing

Viral RNA was extracted directly from supernatants of DENV-infected C6/36 cells using the QIAamp Viral RNA mini kit (Qiagen). The RNA was transcribed to cDNA using RevertAid reverse transcriptase (Thermo Scientific, USA) and random hexamer primers (Thermo Scientific). The full E-gene (~1485 pb) was amplified by PCR using DENV serotype specific primers [S1 Table] and Thermo Scientific Taq DNA Polymerase according to manufacturer
instructions. Each reaction of DNA amplification produced overlapping amplicons of 863-1520 nucleotides in length. Amplicons were sequenced using a Sanger sequencing commercial service (Macrogen DNA Sequencing Service, Seoul Korea). Sequence assembly was performed with CLC Genomics Workbench 4.5 (CLC Bio, Denmark) and sequences were deposited in GenBank [S2 Table].

## Data selection

Available full-length E-gene accessions with known location and sampling date were retrieved from GenBank (last accessed in March 2017). The total dataset (n = 9.669) was aggregated by serotype (3676 of DENV-1, 2925 of DENV-2, 1746 of DENV-3 and 1322 of DENV-4) and analyzed separately. Identical sequences from the same country and year were removed using the UCLUST algorithm in the USEARCH v.10.0.240 software [26]. The filtered datasets (3452 of DENV-1, 2789 of DENV-2, 1645 of DENV-3 and 1100 of DENV 4 Sequences) were combined with novel Colombian isolates, aligned using MUSCLE v.3.8 [27] and visualized using Bioedit v.7.2.5 software [28]. Preliminary phylogenetic analyses were done using Maximum Likelihood (ML) methods with a non-parametric bootstrap using PhyML v.3.1 software [29] in order to identify genotypes and to focus the analysis in the genotypes that were observed in Colombia [S2 Fig].

Sequences from the same year were downsampled when they fell within a monophyletic group by country and were overrepresented. The resulting datasets (DENV-1 = 118, DENV-2 = 118, DENV-3 = 66 and DENV-4 = 77) had only sequences from the region of the Americas because the phylogenetic trees showed a single introduction of each genotype into this region from non-American countries.

A Likelihood Mapping Analysis in TREE-PUZZLE v.5.3 [30] software and a substitution saturation analysis in DAMBE v.6 [31] showed that there was enough phylogenetic information in each dataset [S1 Fig]. All datasets were free of recombination signal following a pairwise homoplastic index test using SplitsTree v.4.12.6 [32] (pairwise homoplastic index ≤ 0.5 for each dataset).

## Data analysis

Haplotype diversity and nucleotide diversity were calculated using the software DNAsp v.5 [33]. For each dataset, the best-fit model of nucleotide substitution was selected based on Akaike Information Criterion (AIC) using the APE package v.3.4.2 [34] and the R v.3.3.0 software [35]; The Tamura and Nei plus invariant and discrete gamma models were the best-fit nucleotide substitution models for DENV-1 and -2, and the generalized time reversible plus invariant and discrete gamma models were the best-fit nucleotide substitution models for DENV-3 and -4. Regression of root-to-tip genetic distance against sampling time [S3 Table] using TempEst v 1.5 [36], showed that there was sufficient temporal signal in each dataset to proceed with phylogenetic molecular clock analyses. The best-fit clock model (strict vs. relaxed) and the best-fit demographic model (among constant size, exponential growth, skyride and skygrid models), were selected via path sampling (PS) and stepping-stone (SS) methods (100 steps with 5 million iterations each) [37]. In all cases, the uncorrelated lognormal (UCLN) relaxed-clock model was preferred. For the model-based demographic inference, the non-parametric skygrid model was preferred for DENV-1 and DENV-3, and the parametric Exponential growth model for DENV-2 and DENV-4 [S4 Table].

Viral demographic curves were reconstructed for each serotype using either the isolates from Colombia (n = 235) sampled between 1982 and 2015 or the isolates from the region of
Santander (n = 160) collected between 1998 and 2015 in order to have a measure of dengue genetic diversity and its fluctuations over time.

The pattern of DENV spread was identified using a standard continuous-time Markov chain (CTMC) coupled with the Bayesian stochastic search variable selection (BSSVS) procedure [38]. We assumed that the transition rates between locations were reversible (symmetrical model). Different schemes of discrete geographical locations were used [S3 Fig] to account for any sampling bias: in the first scheme, the country of isolation was used. In the second scheme, neighboring countries were grouped in regions except for Colombia and Venezuela: Andes (Bolivia, Ecuador and Peru), Central America, Greater Antilles, Lesser Antilles, North America and Southern Cone (Argentina, Brazil and Paraguay). In the final scheme, the Andean region and Southern cone regions were merged into one region named South American region. A Bayes factor test (BF > 9) was used to recognize well-supported diffusion rates using the SpreaD3 v0.9.6 software [39]; well-supported diffusion routes concerning Colombia were compared among schemes.

Bayesian approach was implemented in BEAST v1.8.3 software package [40]. For every analysis, five independent Markov Chain Monte Carlo (MCMC) were used with at least 100 million generations each. Samples were combined and diagnosed using Tracer v1.6 (http://tree.bio.ed.ac.uk/software/tracer) until the reached effective sample sizes over 200 for all parameters [41]. Maximum clade credibility trees (MCC) were summarized using TreeAnnotator v1.8.1.

# Results

## DENV genotypes circulating in Colombia

A total of 235 full-length E-gene sequences from Colombia viruses isolated in Colombia were used in this study: 193 from the region of Santander (160 from the metropolitan area of Buqueparamanga and 33 from the Province of Ocaña); 31 sequences were from the states of Antioquia (n = 21), Valle del Cauca (n = 6), Bolivar (n = 2), Guaviare (n = 1), and Cesar (n = 1); and 11 sequences for which the exact state was not available.

The nucleotide diversity (π)—the average number of nucleotide differences per site in pairwise comparisons among sequences—for the Colombian strains was similar for DENV-2 and DENV-4 (0.020) and lower for DENV-3 (0.015), and DENV-1 (0.011). The haplotype diversity (h)—the probability that two randomly sampled alleles are different—was high and nearly similar for all serotypes (DENV-1, 0.98; DENV-2, 1.0; DENV-3, 0.99; DENV-4, 0.99).

## Evolutionary history of DENV in Colombia

We documented DENV-1 genotype V, DENV-2 genotype Asian/American, DENV-3 genotype III, and DENV-4 genotype II in Colombia [S2 Fig]. The estimated mean rates of nucleotide substitution per site and per year (expressed as the mean and the 95% highest posterior density interval, HPD) were 6.14 x 10⁻⁴ (5.07–7.11), 7.30 x 10⁻⁴ (6.11–8.58), 8.80 x 10⁻⁴ (6.53–9.66), and 7.48 x 10⁻⁴ (5.96–9.04) for DENV-1 to DENV-4, respectively. These estimates were in line with those previously reported for the same genotypes circulating in the Americas [42–46]. For every genotype, we evidenced multiple introductions into Colombia as detailed below.

## DENV-1

The MCC tree for DENV-1 included 79 E-gene sequences from Colombia and evidenced several introductions into the country [Fig 2]. Two different virus strains isolated in 1985 (AF425616) and 1996 (AF425617) were related to viruses circulating in the Lesser Antilles and
<figure><img src="image_5.png" /><figcaption>**Fig 2. Maximum clade credibility tree of DENV-1.** Branches colored according to the most probable posterior location of their child nodes. Monophyletic clades from Colombia highlighted in grey. Node posterior probability (bold), ancestral state probability and the mean estimated time to the most recent common ancestor (tMRCA) with the respective 95% highest posterior density interval (HPD) are indicated for relevant clades. *: virus strains isolated in 1985 (AF425616) and 1996 (AF425617) available in Genbank.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g002

had an estimated time to the most recent common ancestor (TMRCA) around the early 1970s. Most of the strains isolated during the 1990s and the majority of the strains isolated during the 2000s (all of them sampled in the region of Santander) grouped in a different clade that inter-mixed with strains from Venezuela (1994–2008) and had a TMRCA around 1993 (95% HPD 1992–1995). This clade also included strains from other South America regions (Brazil/Paraguay/Argentina, 2008–2011), Greater Antilles (Dominica/Puerto Rico/Haiti, 2006–2014) and the United States of America (2013). Venezuela appeared as the most probable ancestral location (PP = 1.00) for the latest introductions of DENV-1 into Colombia.

## DENV-2

The majority of Colombian strains (38 out of 42) were from the Santander Region, and the remaining isolates were from the states of Antioquia (2004, n = 1), Guaviare (2005, n = 1) and
<figure><img src="image_6.png" /><figcaption>**Fig 3. Maximum clade credibility tree of DENV-2.** Branches colored according to the most probable posterior location of their child nodes. Monophyletic clades from Colombia highlighted in grey. Node posterior probability (bold), ancestral state probability and the mean estimated time to the most recent common ancestor (tMRCA) with the respective 95% highest posterior density interval (HPD) are indicated for relevant clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g003

Valle del Cauca (2015, n = 1). The MCC tree evidenced at least five different introductions of DENV-2 into the country [Fig 3]. A single Colombian strain isolated in 1993 (location unknown) grouped with viruses from the Lesser Antilles (Trinidad and Tobago from 1997) and together they had a TMRCA around 1988 (95% HPD: 1988–1992). Another two strains isolated Santander (1998) grouped either with a few strains from Venezuela (1990–1996) or with strains from Central American (1999–2012) and, in each case, they had an estimated TMRCA around 1988 (95% HPD: 1986–1989 and 1986–1993, respectively). Viral strains sampled in the 1990s and thenceforth, grouped with the remaining strains from Venezuela (1991–2008) and with one strain from Costa Rica (2003); this group had a TMRCA around 1987 (95% HPD: 1986–1988) and was most likely introduced from the Lesser Antilles (PP = 0.93). Finally, a single strain isolated in 2007 from unknown locality did not fall in the previous group and instead clustered with viruses from Central America with a TMRCA around 2000 (95% HPD: 1999–2001); this case represents the most recent DENV-2 documented introduction into the country and possible a lineage that circulated in low-levels.
## DENV-3

Our analysis included 89 E-gene sequences from Colombia and we identified two distinct lineages that co-circulated [Fig 4]. The first lineage included strains from Santander (2001–2015, n = 38), Norte of Santander (2005–2006, n = 5) Antioquia (2002–2007, n = 8), and one strain from Valle del Cauca (2014) that grouped together with strains from Venezuela (2000–2008) and the Greater Antilles (Puerto Rico, 2006; Cuba, 2001). This lineage was probably introduced from the Greater Antilles (PP = 0.81) and its TMRCA was around 1995 (95% HPD: 1994–1998). The second lineage consisted of the remaining 37 strains from multiple states from Colombia (Santander (2007–2015, n = 18), Antioquia (2003–2009, n = 12); Valle del Cauca (2004–2015, n = 3), Bolivar (2006 and 2015), and two strains for which the precise location was not available (2007, n = 2)) that grouped together with strains from Puerto Rico (2004) and Venezuela (2005). This second lineage was probably introduced from Central America (PP = 0.67) and its TMRCA was around 1999 (95% HPD: 1995–2001).

## DENV-4

The MCC tree for DENV-4 included 25 E-gene sequences from Colombia and evidenced at least two introductions [Fig 5]. The majority (n = 19) of strains were from the Santander

<figure><img src="image_7.png" /><figcaption>**Fig 4. Maximum clade credibility tree of DENV-3.** Branches colored according to the most probable posterior location of their child nodes. Monophyletic clades from Colombia highlighted in grey. Node posterior probability (bold), ancestral state probability and the mean estimated time to the most recent common ancestor (tMRCA) with the respective 95% highest posterior density interval (HPD) are indicated for relevant clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g004
<figure><img src="image_8.png" /><figcaption>**Fig 5. Maximum clade credibility tree of DENV-4.** Branches colored according to the most probable posterior location of their child nodes. Monophyletic clades from Colombia highlighted in grey. Node posterior probability (bold), ancestral state probability and the mean estimated time to the most recent common ancestor (tMRCA) with the respective 95% highest posterior density interval (HPD) are indicated for relevant clades.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g005

Region, two strains (2015) were from the states of Valle del Cauca and Cesar, and there were four strains for which the precise location was not available. The two oldest available isolates from Colombia (KC963424 and GU289913) sampled in 1982 grouped with strains from the Greater Antilles, North, South and Central America, and had a TMRCA around 1977 (95% BCI:1976–1977). Most Colombian strains sampled in 2000–2015 (20 out of 25) grouped with strains from Venezuela (1995–2008), and the majority of strains from other South American countries (Ecuador, Peru and Brazil, 2006–2013). The inferred TMRCA for this group was around 1993 (95% BCI: 1990–1995), and it was most likely introduced from Venezuela (PP = 0.65). The remaining three strains from Colombia, which were isolated in 1996, 1997 and 2006, grouped with one strain from Venezuela (1995), one from Ecuador (2000) and one from Peru (2000). This small group had a TMRCA around 1990 (95% HPD: 1985–1992). The strains from the region of Santander and the neighbor state of Cesar were closely related to strains from Venezuela and Brazil strains, while the single strain from Valle del Cauca (in the southwest part of Colombia) was closely related to strains from Ecuador and Peru.
<figure><img src="image_9.png" /><figcaption>**Fig 6.** Summary of dengue virus (DENV) lineages spread with respect to Colombia. Colors indicate each serotype. Only routes with Bayes factor support > 9 are presented.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0203090.g006

## Well-supported migration rates concerning Colombia

We estimated the most significant migration routes among geographical locations using a discrete Bayesian phylogeographic analysis. To account for any sampling bias, the reconstructions were done using three different geographic schemes—detailed in methods. Well-supported diffusion routes with respect to Colombia were compared among schemes. In all cases, the recovered phylogeographical pattern underscored a strong link between Colombia and Venezuela in terms of viral diffusion. Likewise, the Caribbean region played an important role in the diffusion of DENV2–DENV4. Lastly, Central America was relevant for the diffusion of DENV2, whereas South America was relevant for the diffusion of DENV4.

The most homogenous schema (scheme C) was also analyzed under an asymmetrical model as a proxy for directionality [Fig 6]. This analysis showed that most DENV lineages where introduced into Colombia from the other regions. Colombia appeared to be important for the diffusion of DENV3 into Venezuela and the Greater Antilles and had a significant bidirectional diffusion of DENV-2 lineages with Venezuela and Central America.

## DENV serotypes demographic history

The Bayesian demographic reconstructions for Colombia are shown in Fig 7. Overall, we observed a steady increase in the effective population size (Ne) over time for DENV-2 and
Total of Colombian sequences

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

Total of the Region of Santander's sequences

<figure><img src="image_14.png" /></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

**Fig 7.** DENV Population dynamics in Colombia and the state of Santander (1998–2015). Left panels: relative genetic diversity of each serotype in Colombia estimated using non-parametric Skygrid coalescent analysis. Right panels: relative genetic diversity of each serotype in the state of Santander overlapped with the respective serotype predominance. National dengue epidemics indicated with triangles. Gray colored areas show the periods that were not sampled. *Ne* denotes effective population size in logarithmic scale.

https://doi.org/10.1371/journal.pone.0203090.g007
DENV-4. There were oscillations in the case of DENV-1 and more pronounced ones in the case of DENV-3, but the general trend is of increased diversity over time. In the case of DENV-3, there is an increase in diversity peaking during the period 2001–2005, a sharp drop after that and a steady increase from 2009 onwards. The demographic reconstruction for the region of Santander, which accounted for 82% of the Colombian sequence data, followed the same dynamics as that recovered with the dataset from the Country.

There was an agreement between the increase of *Ne* of DENV-3 and the re-introduction of this serotype into the country in 2001 that led to high DENV-3 prevalence in the following years. The same type of pattern was not observed for DENV-1, which was highly prevalent during the periods of 1998–1999 and 2007–2008. DENV-2 and DENV-4 are respectively the most and less prevalent serotypes isolated in the region and although there have been changes in the prevalence of both of them over time, we did not observe any visible difference in the growth of *Ne*.

Discussion

This study is the first large-scale analysis of the spatial and temporal dynamics of DENV in Colombia during the period 1998–2015 using newly sequenced and curated genetic data retrieved from GenBank. The novel sequences were sampled in distinct locations of Colombia, but most of them came from the region of Santander. We documented the frequent introduc- tion of dengue lineages into Colombia, the ongoing co-circulation of multiple serotypes and lineages within the country and the local viral population dynamics.

Our phylogenetic analysis confirmed the circulation of one genotype per serotype (DENV-1: Genotype V, DENV-2: Asian/American Genotype, DENV-3: Genotype III and DENV-4: Genotype II) and it is consistent with the results of other studies [16–19]. The presence of DENV-4 genotype I was documented in Brazil in 2013 [47], but we did not find evidence of its circulation in Colombia. Despite multiple strain introductions of DENV-1, DENV-2 and DENV-4, our data indicated that only one lineage predominates up to date, suggesting the turnover of viral lineages over time. With the data at hand we could not address if such turn- overs are due to stochastic die-off or other factors; however, the re-emergence of DENV-2 in an dengue-endemic setting from Indonesia was due to the loss or decrease of herd immunity during a 5-year period where DENV-1 predominated [48]. Similarly, changing serotype preva- lence has been associated with lineage extinction and lineage replacements in Thailand and with dengue epidemics in Singapore [49, 50, 51, 52]. Here, we recorded two distinct lineages of DENV-3 that apparently were introduced at the same time, that had the same level of circula- tion and that overlapped in time and space. It may be the case that, the complete susceptibility of the population to this serotype at the time of the introduction allowed two successful inde- pendent transmission chains. Similarly, the co-circulation of two DENV-3 lineages has been reported in the metropolitan area of Medellin, Colombia [18] and the co-circulation of differ- ent lineages has been demonstrated for DENV-1 and DENV-2 in Brazil [53,54].

Our estimated tMRCAs for each serotype preceded (between 5–7 years) the respective offi- cial epidemiologic reports by the National Institute of Health from Colombia [9,12]. This is likely the outcome of a preliminary silent circulation of virus (i.e., individuals experiencing mild to asymptomatic infections) coupled with local passive surveillance for dengue [2]. The DENV surveillance in Colombia relies on the immediate reporting of fatal cases and the weekly routine reporting of symptomatic ones by the health care providers [55]. Consequently, the true burden of dengue might be underestimated once local physicians may fail to report cases quickly and routinely. Such passive surveillance includes neither a thorough serotyping nor the genetic characterization of strains owing to financial costs and time constraints. All
fatal cases are laboratory-confirmed, but only around 10–20% of laboratory confirmations remain mandatory for the non-fatal ones [55]. Thus, health care providers belonging to the public health sentinel surveillance system send only a small fraction of collected serum samples to the National Institute of Health for laboratory confirmation and serotype identification. As a consequence, an epidemic has often reached or passed its peak before it has been recognized or before it can be controlled, resulting in a vague idea of the influence of the serotypes behind it.

The increased regional commerce and travel in dengue hyper-endemic Latin America likely drive frequent exchanges and importation of DENV through human movement. For instance, it has been shown an effect of air travel in increased dengue transmission in the Americas and Asia [43, 53, 56]. Here, We observed frequent viral lineage exchange among Latin American countries and persistent co-circulation of dengue viral lineages. Human flow in Colombia is also particularly significant as the country is the largest sending country of migrants in South America [57,58]. However, the country seems to act as a sink instead of a source for DENV lineages. The strong association with Venezuela probably reflects geographical proximity, connectivity through trade and migratory history (A significant migrant wave occurred in the mid 80’s to Venezuela, primarily motivated by its economic boom and by the economic difficulties of Colombia at that time). However, we also do not rule out the possibility that this association could be a result of intense sampling in the northeast of Colombia. For example, the DENV-4 sequence from Valle del Cauca (southwestern Colombia) was related to those from Peru and Ecuador. Nevertheless, the region of Santander represents a point of transit with the neighboring country of Venezuela and thus substantial morbidity due to cross-border exchange of dengue-infected people is expected to continue.

Given that DENV transmission is similar to other arboviruses, the aforementioned high levels of viral exchange are consistent with the rapid and recent establishment of chikungunya (CHIKV) and Zika (ZIKV) in the region [59,60,61]. However, whereas the dynamics of DENV are complicated by the four serotypes and complex susceptibility profiles of local populations, recurring CHIKV and ZIKV outbreaks will be functionally related to the turnover rate of susceptible humans in addition to reintroduction rates [62]. These observations underscore that arboviral diseases are a recurring public health problem in Colombia and recognizing the importance of their molecular investigation will strengthen diagnosis and epidemiological integrated surveillance.

The population dynamics of DENV can be driven by viral introduction events that could eventually lead to lineage turnovers as have been documented in other settings [43,53,63,64, 65]. For instance, regular stochastic clade replacement led to recurring homotypic dengue outbreaks in Malaysia [66]. In the same vein, constant DENV introductions and *in situ* evolution contribute to viral diversity in Singapore where replacement of a predominant viral clade, even in the absence of a switch in the predominant serotype, hinted a possible increase in transmission [52]. In our study, and in particular in the Santander Region, we observed the consistently co-circulation of the four serotypes over consecutive years. This underscores an intricate pattern of competition that does not result in complete serotype displacement and may explain the increase of diversity over time for the four serotypes. We observed much rapid rise in the *Ne* of DENV-3 followings its introduction that was also reflected in the pre-dominance of this serotype in the period 2001–2004. Thus, in the Santander Region, both *in situ* evolution and the recurrent introduction of lineages drive local dengue dynamics.

This study had limitations: most of the samples were collected from a single region and this limits our knowledge on the country-wide viral genetic diversity. Likewise, samples from other countries and times are quite heterogeneous and this might also hide additional routes and the relevance of specific regions from the establishment and transmission of dengue in Colombia
and the Americas. Furthermore, giving scarcity of resources, limited data on serotype prevalence restricts the conclusion that could be drawn on the dynamics of DENV. Larger series with long-term follow-up are needed to confirm the effect that a particular serotype might have during epidemics. Nonetheless, even with these limitations, we showed that passive laboratory-based disease surveillance studies allowed us to have a general picture of the diversity and dynamics of DENV.

## Conclusion

We characterized the spatial-temporal dynamics of all dengue serotypes in a highly endemic area of Colombia, documented the co-circulation of a single genotype of each serotype, and the unapparent circulation of every serotype prior its first detection in the country. This study also showed that genetic diversity of serotypes circulating in the country continues to grow due to *in-situ* evolution and recurring introductions of viral strains from different countries in the region. Our study advances the countrywide genomic surveillance to lay the groundwork for the introduction of dengue vaccines and future control initiatives and underscores the continued need for a sensitive surveillance system aiming to collect reliable baseline data.

## Supporting information

**S1 Fig. Phylogenetic signal analysis.** Likelihood mapping diagrams were obtained using TREE-PUZZLE. All datasets had between 67% and 91,3% proportion of well resolved quarters. The index of substitution saturation calculated with DAMBE was below 1.0 (0,04–0,1) and smaller than the critical index of substitution saturation (0,4–0,8) indicating that the sequences did not experience severe substitution saturation.
(SVG)

**S2 Fig. Phylogeographical reconstruction.** Three different schemes were analyzed under symmetrical models. The most homogenous schema (scheme C) was also analyzed under an asymmetrical model. Panels next to maps show the number of sequences included. BF denotes Bayes Factor.
(SVG)

**S3 Fig. Phylogenetic trees constructed with Maximum Likelihood methods.** Trees were rooted using Yellow fever virus; every tree included representative taxa for each serotype and representative taxa for the corresponding genotype. Colored tips represent the samples from Colombia: DENV-1 Genotype V in blue, DENV-2 Genotype Asian/American in red, DENV-3 Genotype III in green, and DENV-4 Genotype II in yellow.
(SVG)

**S1 Table. Primers using to amplify and sequence DENV E-genes.**
(DOCX)

**S2 Table. DENV E gene sequences from Colombian viruses generated in this study.**
(DOCX)

**S3 Table. Clock-likeness using TempEst v 1.5.**
(DOCX)

**S4 Table. Log marginal likelihood estimates for different coalescent model combinations.**
(DOCX)
## Acknowledgments

This project was possible thanks to support from coordinators in healthcare institutions and diagnostic laboratories in Santander and Ocaña. We thank Luis Angel Villar from Unidad de Epidemiología Clínica, Escuela de Medicina, Universidad Industrial de Santander, for his supporting virological surveillance activities in 2014–2016. Financial support: This study was supported by the Colombian Government: Science, Technology and Innovation Fund of Sistema General de Regalias, Santander, Casanare, Valle del Cauca (BPIN 2013000100011).

## Author Contributions

**Conceptualization:** Cinthy L. Jiménez-Silva, Christian Julián Villabona-Arenas, Raquel E. Ocazionez.

**Data curation:** Cinthy L. Jiménez-Silva, María Fernanda Carreño, Ayda Susana Ortiz-Baez, Luz Aida Rey.

**Formal analysis:** Cinthy L. Jiménez-Silva.

**Funding acquisition:** Raquel E. Ocazionez.

**Investigation:** Cinthy L. Jiménez-Silva, María Fernanda Carreño, Ayda Susana Ortiz-Baez.

**Methodology:** Cinthy L. Jiménez-Silva, Christian Julián Villabona-Arenas.

**Project administration:** Raquel E. Ocazionez.

**Resources:** Raquel E. Ocazionez.

**Supervision:** Christian Julián Villabona-Arenas, Raquel E. Ocazionez.

**Visualization:** Cinthy L. Jiménez-Silva, María Fernanda Carreño.

**Writing – original draft:** Cinthy L. Jiménez-Silva.

**Writing – review & editing:** Cinthy L. Jiménez-Silva, Christian Julián Villabona-Arenas, Raquel E. Ocazionez.

## References

1. Sharp TM, Tomashek KM, Read JS, Margolis HS, Waterman SH. A New Look at an Old Disease: Recent Insights into the Global Epidemiology of Dengue. Curr Epidemiol Reports. 2017; 4(1):11–21.

2. Rodriguez-Roche R, Gould EA. Understanding the dengue viruses and progress towards their control. Biomed Res Int. 2013; 2013:20.

3. Stanaway JD, Shepard DS, Undurraga EA, Halasa YA, Coffeng LE, Brady OJ, et al. The Global Burden of Dengue: an Analysis from the Global Burden of Disease Study 2013. Lancet Infect Dis. 2016; 16 (6):712–23. https://doi.org/10.1016/S1473-3099(16)00026-8 PMID: 26874619

4. Weaver SC, Vasilakis N. Molecular evolution of dengue viruses: Contributions of phylogenetics to understanding the history and epidemiology of the preeminent arboviral disease. Infect Genet Evol. 2009; 9(4):523–40. https://doi.org/10.1016/j.meegid.2009.02.003 PMID: 19460319

5. Rico-Hesse R. Microevolution and virulence of dengue viruses. Adv Virus Res. 2003; 59:315–41. PMID: 14696333

6. San Martín JL, Brathwaite O, Zambrano B, Solórzano JO, Bouckenooghe A, Dayan GH, et al. The epidemiology of dengue in the Americas over the last three decades: A worrisome reality. Am J Trop Med Hyg. 2010; 82(1):128–35. https://doi.org/10.4269/ajtmh.2010.09-0346 PMID: 20065008

7. Villar LA, Rojas DP, Besada-Lombana S, Sarti E. Epidemiological Trends of Dengue Disease in Colombia (2000–2011): A Systematic Review. PLoS Negl Trop Dis. 2015; 9(3):e0003499. https://doi.org/10.1371/journal.pntd.0003499 PMID: 25790245

8. PAHO / WHO: Pan-American Health Organization / World Health Organization. Dengue: data, maps and statistics. Available from: https://www.paho.org/hq/index.php?option=comTopics&view=readall&id=3274&lte=4073lang=es. Accessed 29 April, 2018.
9. Padilla JC, Rojas-Álvarez DP, Sáenz-Gómez R. Dengue en Colombia, epidemiología de la reemergencia a la hiperendemia [Dengue in Colombia, epidemiology of the reemergency towards hyperendemia]. Bogotá, D.C.; Guías de Impresión Ltda.: 2012. p.142–5.

10. Pinheiro FP, Corber SJ. Global situation of dengue and dengue haemorrhagic fever, and its emergence in the Americas. World Health Stat Q. 1997; 50(3-4):161-9 PMID: 9477544

11. Ocazionez RE, Cortés FM, Villar LA, Gómez SY. Temporal distribution of dengue virus serotypes in Colombian endemic area and dengue incidence. Re-introduction of dengue-3 associated to mild febrile illness and primary infection. Mem Inst Oswaldo Cruz. 2006; 101(7):725-31.9. PMID: 17160279

12. SIVIGILA: Sistema de Vigilancia, Instituto Nacional de Salud. [National Public Health Surveillance System database]. Available from: https://www.minsalud.gov.co/salud/Paginas/SIVIGILA.aspx [Web Page]. Accessed 29 March, 2018

13. Villabona-arenas CJ, Elvira R, Jimenez O, Silva LJ. Dengue Vaccine: Considerations before Rollout in Colombia. PLoS Negl Trop Dis. 2016;5-9.

14. Rangel SYG, Villabona-Arenas CJ, Pimiento F a T, Miranda-Esquivel DR, Jimenez REO. Dengue virus serotype 3 (genotype III) from Colombia: A perspective of its pathogenic potential. Dengue Bull. 2008; 32(genotype III):126-37.

15. Ocazionez RE, Ortiz AS, Gómez SY, Miranda DR. Virus del dengue serotipo 1 (DENV-1) de Colombia: su contribución a la presentación del dengue en el departamento de Santander [Dengue virus serotype 1 (DENV-1) of Colombia: its contribution to dengue occurrence in the Santander State]. Biomédica. 2013; 33(1):22-30.

16. Méndez J a, Usme-Ciro J a, Domingo C, Rey GJ, Sánchez J a, Tenorio A, et al. Phylogenetic reconstruction of dengue virus type 2 in Colombia. Virol J. 2012; 9(1):64.

17. Méndez JA, Usme-Ciro JA, Domingo C, Rey GJ, Sánchez JA, Tenorio A, et al. Phylogenetic history demonstrates two different lineages of dengue type 1 virus in Colombia. Virology J. 2010; 7:226. https://doi.org/10.1186/1743-422X-7-22.

18. Ospina MC, Diaz FJ, Osorio JE. Prolonged co-circulation of two distinct Dengue virus type 3 lineages in the hyperendemic area of Medellin, Colombia. Am J Trop Med Hyg. 2010; 83(3):672-8. https://doi.org/10.4269/ajtmh.2010.09-0766 PMID: 20810837

19. Villabona-Arenas CJ, Miranda-Esquivel DR, Jimenez REO. Phylogeny of dengue virus type 3 circulating in Colombia between 2001 and 2007. Trop Med Int Hea!. 2009; 14(10):1241-50.

20. Pérez G, Aguilera M, Otero A, Sánchez A, Acosta K, Galvis L. Economía de las grandes ciudades en Colombia: seis estudios de caso [Economy of the big cities in Colombia: six case studies]. Colección de Economía Regional, Banco de la República. 2014. 161 p.

21. Departamento Nacional de Planeación (DANE). Anuario de movimientos internacionales de viajeros de Colombia 2013. [Yearbook of international traveler movements of Colombia 2013]. Available from: http://migracioncolombia.gov.co/phocadownload/Anuario%20de%20Movimientos%20Internacionales.pdf.

22. Villar L, Dayan GH, Arredondo-García JL, Rivera DM, Cunha R, Deseda C, et al. Efficacy of a Tetravalent Dengue Vaccine in Children in Latin America. N Engl J Med. 2015; 372(2):113-23. https://doi.org/10.1056/NEJMoa1411037 PMID: 25365753

23. Alcaldía de Ocaña, Norte de Santander [Ocaña City], Hall, Norte de Santander]. Available from: www.ocana-nortedesantander.gov.co/informacion_general.shtm.

24. Osorio L, Ramirez M, Bonelo A, Villar LA, Parra B. Comparison of the diagnostic accuracy of commercial NS1-based diagnostic tests forearly dengue infection Virol J. 2010 Dec 6; 7:361. https://doi.org/10.1186/1743-422X-7-361 PMID: 21134275

25. Villamor E, Villar LA, Lozano A, Herrera VM, Herrán OF. Vitamin D serostatus and dengue fever progression to dengue hemorrhagic fever/dengue shock syndrome. Epidemiol Infect. 2017; 145(14):2961-70. https://doi.org/10.1017/S0950268817002059 PMID: 28903788

26. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010; 26(19):2460-1. https://doi.org/10.1093/bioinformatics/btq461 PMID: 20709691

27. Edgar R. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004; 5(1):113.

28. Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999; 41:95-8.

29. Guindon S, Dufayard J-F, Lefort V, Anisimova M. New Alogrithms and Methods to Estimate Maximum-Likelihoods Phylogenies: Assessing the performance of PhyML 3.0. Syst Biol. 2010; 59(3):307-21. https://doi.org/10.1093/sysbio/syq010 PMID: 20525638

30. Schmidt HA, Strimmer K, Vingron M. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. 2002; 18(3):502-4. PMID: 11934758
31. Xia X. DAMBE6: New tools for microbial genomics, phylogenetics and molecular evolution. Journal of Heredity. 2017; 108(4):431-437. https://doi.org/10.1093/jhered/esx033 PMID: 28379490

32. Bruen TC, Philippe H, Bryant D. A simple and robust statistical test for detecting the presence of recom-bination. Genetics. 2006; 172(4):2665-81. https://doi.org/10.1534/genetics.105.048975 PMID: 16489234

33. Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009; 25(11):1451-2. https://doi.org/10.1093/bioinformatics/btp187 PMID: 19346325

34. Paradis E, Claude J, Strimmer K. APE: Analyses of phylogenetics and evolution in R language. Bioinformatics. 2004; 20(2):289-90. PMID: 14734327

35. R Core team. R Core Team. Vol. 55, R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL http://www.R-project.org/. 2015. p. 275-86.

36. Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol. 2016; 2(1):vew007. https://doi.org/10.1093/ve/vew007 PMID: 27774300

37. Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko A V., et al. Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol. 2013; 29(5):1970-9.

38. Lemey P, Rambaut A, Drummond AJ, Suchard M a. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009; 5(9).

39. Bielejec F, Baele G, Vrancken B, Suchard MA, Rambaut A, Lemey P. SpreaD3: Interactive Visualization of Spatiotemporal History and Trait Evolutionary Processes. Mol Biol Evol. 2016; 33(8):2167-9. https://doi.org/10.1093/molbev/msw082 PMID: 27189542

40. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012; 29(8):1969-73. https://doi.org/10.1093/molbev/mss075 PMID: 22367748

41. Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006; 4(5):699-710.

42. Mir D, Romero H, Max L, Carvalho F De, Bello G. Spatiotemporal dynamics of DENV-2 Asian-American genotype lineages in the Americas. PLoS ONE 2014; 9(6): e98519. https://doi.org/10.1371/journal.pone.0098519 PMID: 24897118

43. Allicock OM, Lemey P, Tatem AJ, Pybus OG, Bennett SN, Mueller B, et al. Phylogeography and population dynamics of dengue viruses in the Americas. Mol Biol Evol. 2012; 29(6):1533-43. https://doi.org/10.1093/molbev/msr320 PMID: 22319149

44. Villabona-Arenas CJ, Zanotto PMD a. Evolutionary history of Dengue virus type 4: Insights into genotype phylodynamics. Infect Genet Evol. 2011; 11(5):878-85. https://doi.org/10.1016/j.meegid.2011.02.007 PMID: 21335103

45. Bennett SN, Drummond a. J, Kapan DD, Suchard M a., Munoz-Jordán JL, Pybus OG, et al. Epidemic dynamics revealed in dengue evolution. Mol Biol Evol. 2010; 27(4):811-8. https://doi.org/10.1093/molbev/msp285 PMID: 19965886

46. Araújo JMG, Nogueira RMR, Schatzmayr HG, Zanotto PMD a, Bello G. Phylogeography and evolutionary history of dengue virus type 3. Infect Genet Evol. 2009; 9(4):716-25. https://doi.org/10.1016/j.meegid.2008.10.005 PMID: 19010450

47. Campos R de M, Veiga CSB, Meneses MDF, de Souza LM, Fernandes CAS, Malirat V, et al. Emergence of Dengue virus 4 genotypes II b and I in the city of Rio de Janeiro. J Clin Virol. 2013; 56(1):86-8. https://doi.org/10.1016/j.jcv.2012.10.006 PMID: 23127562

48. Zhang C, Jr MPM, Klungthong C, Monkongdee P, Nimmannitya S, Kalayanarooj S, et al. Clade Replacements in Dengue Virus Serotypes 1 and 3 Are Associated with Changing Serotype Prevalence Clade Replacements in Dengue Virus Serotypes 1 and 3 Are Associated with Changing Serotype Prevalence†. J Virol. 2005; 79(24):15123-30. https://doi.org/10.1128/JVI.79.24.15123-15130.2005 PMID: 16306584

49. Kotaki T, Yamanaka A, Mulyatno KC, Churrotin S, Sucipto TH, Labiqah A, et al. Divergence of the dengue virus type 2 Cosmopolitan genotype associated with two predominant serotype shifts between 1 and 2 in Surabaya, Indonesia, 2008-2014. Infect Genet Evol. 2016; 37:88-93. https://doi.org/10.1016/j.meegid.2015.11.002 PMID: 26553170

50. Koh BK, Ng LC, Kita Y, Tang CS, Ang LW, Wong KY, et al. The 2005 dengue epidemic in Singapore: epidemiology, prevention and control. Ann Acad Med Singapore. 2008; 37(7):538-45. PMID: 18695764

51. Lee KS, Lai YL, Lo S, Barkham T, Aw P, Ooi PL, et al. Dengue virus surveillance for early warning, Singapore. Emerg Infect Dis. 2010; 16(5):847-9. https://doi.org/10.3201/eid1605.091006 PMID: 20409381
52. Lee KS, Lo S, Tan SS, Chua R, Tan LK, Xu H, et al. Dengue virus surveillance in Singapore reveals high viral diversity through multiple introductions and in situ evolution. *Infect Genet Evol.* 2012; 12 (1):77–85. https://doi.org/10.1016/j.meegid.2011.10.012 PMID: 22036707

53. Nunes MRT, Palacios G, Faria NR, Sousa EC, Pantoja JA, Rodrigues SG, et al. Air travel is associated with intracontinental spread of dengue virus serotypes 1–3 in Brazil. *PLoS Negl Trop Dis.* 2014 Apr 17; 8(4):e2769. https://doi.org/10.1371/journal.pntd.0002769 PMID: 24743730

54. Drumond BP, Mondini A, Schmidt DJ, de Bronzoni RVM, Bosch I, Nogueira ML. Circulation of different lineages of dengue virus 2, genotype American/Asian in Brazil: dynamics and molecular and phyloge- netic characterization. *PLoS One*. 2013; 8(3):e59422. https://doi.org/10.1371/journal.pone.0059422 PMID: 23533624

55. Instituto Nacional de Salud. Instrucciones para la intensificación de las acciones de vigilancia, preven- ción, atención y control del dengue y dengue grave en Colombia [Instructions for the intensification of surveillance, prevention, care, and control of dengue and severe dengue in Colombia]. Available from: https://www.minsalud.gov.co/sites/rid/Lists/.../RIDE/DE/.../Circular-0008-de-2012.pdf.

56. Tian H, Sun Z, Faria NR, Yang J, Cazelles B, Huang S, et al. Increasing airline travel may facilitate co-circulation of multiple dengue virus serotypes in Asia. *Althouse B*, editor. *PLoS Negl Trop Dis.* 2017 Aug 3; 11(8):e0005694. https://doi.org/10.1371/journal.pntd.0005694 PMID: 28771468

57. International Organization for migration Regional Office for South America. Perfil migratorio de Colom- bia 2012 [Migration profile of Colombia 2012]. Ramírez C, Mendoza L, editors. Bogotá; 2013. 141 p. https://www.iom.int/files/live/sites/iom/.../Perfil-Migratorio-de-Colombia-2012.pdf.

58. (DANE) DAN de E. Estadisticas de demografia y poblacion [Internet]. Available from: http://www.dane.gov.co/index.php/estadisticas-por-tema/demografia-y-poblacion

59. Thézé J, Li T, du Plessis L, Bouquet J, Kraemer MUG, Somasekar S, et al. Genomic Epidemiology Reconstructs the Introduction and Spread of Zika Virus in Central America and Mexico. *Cell Host Microbe.* 2018 13; 23(6):855–864.e7. https://doi.org/10.1016/j.chom.2018.04.017 PMID: 29805095

60. Grubaugh ND, Ladner JT, Kraemer MUG, Dudas G, Tan AL, Gangavarapu K, et al. Genomic epidemiol- ogy reveals multiple introductions of Zika virus into the United States. *Nature.* 2017 15; 546(7658):401– 405. https://doi.org/10.1038/nature22400 PMID: 28538723

61. Faria NR, Quick J, Claro IM, Thézé J, de Jesus JG, Giovanetti M5, et al. Establishment and cryptic transmission of Zika virus in Brazil and the Americas. *Nature.* 2017 15; 546(7658):406–410. https://doi.org/10.1038/nature22401 PMID: 28538727

62. Christofferson RC. Zika Virus Emergence and Expansion: Lessons Learned from Dengue and Chikun- gunya May Not Provide All the Answers. *Am J Trop Med Hyg.* 2016 6; 95(1): 15–18. https://doi.org/10. 4269/ajtmh.15-0866 PMID: 26903610

63. Santiago GA, McElroy-Horne K, Lennon NJ, Santiago LM, Birren BW, Henn MR, et al. Reemergence and decline of dengue virus serotype 3 in Puerto Rico. *J. Infect Dis.* 2012; 206(6):893–901. https://doi.org/10.1093/infdis/jis426 PMID: 22740715

64. Afreen N, Naqvi IH, Broor S, Ahmed A, Kazim SN, Dohare R, et al. Evolutionary Analysis of Dengue Serotype 2 Viruses Using Phylogenetic and Bayesian Methods from New Delhi, India. *PLoS Negl Trop Dis.* 2016 15; 10(3):e0004511. https://doi.org/10.1371/journal.pntd.0004511 PMID: 26977703

65. Manakkadan A, Joseph I, Prasanna RR, Kunju RI, Kailas L, Sreekumar E. Lineage shift in Indian strains of Dengue virus serotype-3 (Genotype III), evidenced by detection of lineage IV strains in clinical cases from Kerala. *ViroJ.* 2013 29; 10:37. https://doi.org/10.1186/1743-422X-10-37 PMID: 23360780

66. Teoh BT, Sam SS, Tan KK, Johari J, Shu MH, Danlami MB, et al. Dengue virus type 1 clade replace- ment in recurring homotypic outbreaks. *BMC Evol Biol.* 2013 28; 13:213. https://doi.org/10.1186/1471- 2148-13-213 PMID: 24073945