2023-24 dengue outbreak in Valle del Cauca, Colombia caused by multiple virus serotypes and lineages ==================================================================================================== * Nathan D. Grubaugh * Daniela Torres-Hernández * Mónica A. Murillo-Ortiz * Diana M. Dávalos * Pio Lopez * Isabel C. Hurtado * Mallery I. Breban * Ellie Bourgikos * Verity Hill * Eduardo López-Medina ## Abstract Global dengue cases rapidly rose to record levels in 2023-24. We investigated this trend in Valle del Cauca, Colombia to determine if specific dengue virus serotypes or lineages were responsible for the large outbreak. We detected all four serotypes and multiple lineages, suggesting that other factors, such as climatic conditions, are likely responsible. ## Introduction Reported cases caused by dengue virus (DENV: *Flaviviridae*; *Orthoflavivirus*), composed of four genetically distinct serotypes (DENV-1-4), are on the rise. In 2023, 4.6 million dengue cases were reported, a record at the time and a 64% increase over 2022 (1). Those numbers were quickly surpassed in 2024 with almost 10 million dengue cases through June. Of those, ∼8.4 were reported from Brazil alone; however, large outbreaks have been reported by many countries, including Colombia. The Valle del Cauca State Health Department in Colombia reported ∼56,000 dengue cases through May, 2024, compared to ∼23,000 for all of 2023 and less than 5,000 in 2022 (**Figure 1A**). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.16.24310413/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2024/07/18/2024.07.16.24310413/F1) Figure 1. Dengue cases by virus serotype and lineage in Valle del Cauca, Colombia (**A**) Monthly dengue cases reported by Valle del Cauca State Health Department in Colombia. Samples from confirmed dengue cases (*n*=266) sorted by DENV (**A**) serotype per month by RT-qPCR (unknown = below limits of PCR detection) and (**B**) lineage by amplicon-based sequencing (listed as serotype, genotype, and lineage). The cause of the substantial increase in reported dengue cases is likely multifaceted. Warming temperatures caused by climate change are known to increase the transmission potential and expand the geographical range of the primary mosquito vector, *Aedes aegypti* (2). Moreover, Indian Ocean sea surface temperature anomalies, especially El Niño events, are associated with dengue epidemics (3). There was a powerful El Niño–Southern Oscillation event in 2023-24, the first since 2015-16. Moreover, new DENV introductions, perhaps due to the resumption of travel post the COVID-19 pandemic (4), could be reaching large susceptible populations. For example, DENV-3 was rarely detected in the Americas for the last 10 years before a recent introduction into the Caribbean from Asia (5,6). To investigate if there was a specific DENV serotype or lineage contributing to the recent surge in cases, we sequenced DENV infections diagnosed in Valle del Cauca, Colombia from April 2023 to May 2024. ## Methods and Results The Valle del Cauca State Health Department in Colombia reported 966 dengue cases in 2019, 655 in 2020, 8,940 in 2021, 4,630 in 2022, 22,988 in 2023, and 56,355 cases through May, 2024 (**Figure 1A**). To determine what DENV serotypes and lineages caused some of the recent cases, we collected serum (150-500 µL) from 266 confirmed dengue cases diagnosed at Hospital Universitario del Valle (HUV) from April 2023 to May 2024, aged <1 to 77 years (average = 16 years). The samples were shipped to Yale University in New Haven, Connecticut for molecular processing. Our study was approved by the Research Ethics Committee at HUV and Yale University Human Research Protection Program (Protocol ID: 2000033281) and all participants signed an informed consent. We extracted RNA from 140 µL of serum per sample using the QIAamp Viral RNA Mini Kit. We initially determined the DENV serotype using a multiplexed RT-qPCR assay (7) before attempting pan-serotype whole-genome amplicon sequencing with DengueSeq (8). Bioinformatic analysis, including primer trimming and consensus sequence generation, was conducted using an iVar pipeline as described in Vogels et al. (8). We assigned DENV lineages to samples with at least 5% genome completeness using the new nomenclature system (9) implemented in Genome Detective ([dengue-lineages.org](http://dengue-lineages.org)). We assigned serotypes to 185 of the 266 samples (70%; **Figure 1B**) and assigned lineages to 171 samples (64%; **Figure 1C**; **Appendix Table**). We detected all four DENV serotypes from the 2023-24 samples (DENV-1 = 35, DENV-2 = 85, DENV-3 = 63, DENV-4 = 2, unknown = 81; **Figure 1B**). DENV-1, -2, and -3 were of relatively equal proportions for part of 2023, but then DENV-1 decreased as DENV-2 increased in frequency from late 2023 to early 2024. We also detected multiple lineages per serotype (except for DENV-4), with DENV-3 genotype III lineage C.1 (3III_C.1, according to (9)), DENV-2III_D.2, DENV-2II_F.1.2, and DENV-1V_D.1 being the most common (**Figure 1C**). To further investigate the DENV lineages, we performed phylogenetic analysis using 79 sequenced samples for which we achieved >70% genome coverage (DENV-1 = 10 sequences, DENV-2 = 38, DENV-3 = 31, DENV-4 = 0). We combined our data with a background dataset downloaded from GenBank, and then downsampled the data per serotype so that we kept all sequences from Colombia, 5 per year from the other countries in the Americas, and 1 per year from each country outside of the Americas (DENV-1 = 1007 sequences, DENV-2 = 965, DENV-3 = 542). We analyzed the sequences using the Nextstrain bioinformatic and phylogenetic framework (10), including constructing maximum likelihood trees using IQTree (11) with a GTR substitution model. Our sequencing data can be found at NCBI BioProject PRJNA1132139 and the GenBank accession numbers are listed in **Appendix Table**. Our alignments, trees, and Nextstrain outputs are available at: github.com/grubaughlab/DENV-genomics/tree/master/paper_2024-CO. Our DENV-1 phylogenetic analysis reveals the co-circulation of two distinct lineages, DENV-1V_D.1 and D.2 (**Figure 2**). Both have been previously detected in Colombia and elsewhere in South America (5,12), representing ongoing local and regional persistence and diversification of these lineages for the past ∼15-20 years. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.16.24310413/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2024/07/18/2024.07.16.24310413/F2) Figure 2. Time-resolved maximum-likelihood phylogeny of DENV-1. (**A**) Full reconstruction of 1007 DENV-1 sequences colored by sampling location. (**B**) Zoom on the DENV-1V_D clade highlighting sequences from Cali, Colombia (black). Our phylogenetic analysis of DENV-2 presents a more complicated picture of three genetic clusters and three individual sequences dispersed among two defined lineages, DENV-2III\_D.2 and DENV-2II\_F.1.2 (**Figure 3**). Lineage 2III_D.2 is a descendent of the original DENV-2 genotype III (“Asian-American”) introduced in the Americas during the late 1970s that subsequently became established throughout the region, including Colombia (13). DENV-2 genotype II (“Cosmopolitan”) was recently introduced into the Americas from Asia, first detected during a dengue outbreak in Peru in 2019 (14). Detection of DENV-2II_F.1.2 in Valle del Cauca demonstrates that the emerging “Cosmopolitan” genotype can become established alongside the existing “Asian-American” genotype. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.16.24310413/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2024/07/18/2024.07.16.24310413/F3) Figure 3. Time-resolved maximum-likelihood phylogeny of DENV-2. (**A**) Full reconstruction of 965 DENV-2 sequences colored by sampling location. Zoom on the (**B**) DENV-2III_D.2 and (**C**) DENV-2II_F.1.2 clades highlighting sequences from Cali, Colombia (black). One hypothesis for the sudden increase in dengue cases was the introduction and rapid spread of a new DENV-3 lineage from Asia (5). DENV-3 in the Americas can often go decade or more between detections (4,6,15). Therefore, detecting an emerging lineage (DENV-3III_B.3) in the Caribbean (5), Brazil (6), Nicaragua (4), and elsewhere in the Americas was of concern. While we detected one dengue case from Valle del Cauca in January 2024 with a likely 3III_B.3 infection (18% genome coverage), 61 of 63 (97%) of DENV-3 infections were lineage 3III_C.1 (**Figure 1C**). DENV-3III_C was likely first introduced into the Americas in the early 1990s (13) and we show that it has persisted through long periods of low detection (**Figure 4**). This includes sporadic detections of 3III\_C.1 in Colombia since the early 2000s. Therefore, our results suggest that populations in the Americas may be susceptible to DENV-3 in general and not just the emerging 3III\_B.3 lineage. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/18/2024.07.16.24310413/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2024/07/18/2024.07.16.24310413/F4) Figure 4. Time-resolved maximum-likelihood phylogeny of DENV-3. (**A**) Full reconstruction of 542 DENV-3 sequences colored by sampling location. (**B**) Zoom on the DENV-3III_C.1 clade highlighting sequences from Cali, Colombia (black). ## Conclusions We investigated the DENV infections from Valle del Cauca, Colombia to determine if a specific virus serotype or lineage was in part driving the record number of dengue cases (1). We detected all four serotypes, with DENV-1, -2, and -3 sharing dominance, and at least eight separate defined lineages (9). These included multiple lineages of DENV-1 genotype V and DENV-2 genotype III (“Asian-America”) which have circulated in the Americas for ∼40 years (13), as well as the emerging DENV-2II_F.1.2 (“Cosmopolitan”) lineage. Moreover, despite the rapid spread of a new DENV-3III_B.3 lineage in the Americas (4–6), we found that the dominant DENV-3 lineage was 3III_C.1 that has been sporadically detected in Colombia for ∼20 years. Our results show that the conditions for transmission in Valle del Cauca are highly suitable for a diversity of DENV serotypes and lineages, and thus that the specific viruses are not the primary driver of the large outbreak. ## Data Availability Our sequencing data can be found at NCBI BioProject PRJNA1132139 and the GenBank accession numbers are listed in Appendix Table. Our alignments, trees, and Nextstrain outputs are available at: github.com/grubaughlab/DENV-genomics/tree/master/paper_2024-CO. [https://www.github.com/grubaughlab/DENV-genomics/tree/master/paper\_2024-CO](https://www.github.com/grubaughlab/DENV-genomics/tree/master/paper_2024-CO) ## Acknowledgement This publication was made possible by the National Institute of Allergy and Infectious Diseases of the National Institutes of Health (NIH) under Award Number DP2AI176740 (NDG) and an NIH Shared Equipment grant 1S10OD028669-01 award to the Yale Center for Genome Analysis. The findings and conclusions in this report are those of the author(s) and do not necessarily represent the official position of the NIH. ## Footnotes * We fixed oversights in the authorship: added Pio Lopez from CEIP and corrected Eduardo Lopez-Medina's name * Received July 16, 2024. * Revision received July 18, 2024. * Accepted July 18, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.PAHO/WHO Data - Dengue [Internet]. [cited 2024 Jul 1]. Available from: [https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en.html](https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en.html) 2. 2.Messina JP, Brady OJ, Golding N, Kraemer MUG, Wint GRW, Ray SE, et al. The current and future global distribution and population at risk of dengue. Nature microbiology [Internet]. 2019 Sep [cited 2024 Jun 28];4(9). Available from: [https://pubmed.ncbi.nlm.nih.gov/31182801/](https://pubmed.ncbi.nlm.nih.gov/31182801/) 3. 3.Chen Y, Xu Y, Wang L, Liang Y, Li N, Lourenço J, et al. Indian Ocean temperature anomalies predict long-term global dengue trends. Science [Internet]. 2024 May 10 [cited 2024 Jul 1]; Available from: [https://www.science.org/doi/10.1126/science.adj4427](https://www.science.org/doi/10.1126/science.adj4427) 4. 4.Cerpas C, Vásquez G, Moreira H, Juarez JG, Coloma J, Harris E, et al. Introduction of New Dengue Virus Lineages of Multiple Serotypes after COVID-19 Pandemic, Nicaragua, 2022 - Volume 30, Number 6— June 2024 - Emerging Infectious Diseases journal - CDC. [cited 2024 Jul 1]; Available from: [https://www.nc.cdc.gov/eid/article/30/6/pdfs/23-1553.pdf](https://www.nc.cdc.gov/eid/article/30/6/pdfs/23-1553.pdf) 5. 5.Taylor-Salmon E, Hill V, Paul LM, Koch RT, Breban MI, Chaguza C, et al. Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean. Nat Commun. 2024 Apr 25;15(1):1–14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-024-45766-2&link_type=DOI) 6. 6.Naveca FG, Santiago GA, Maito RM, Ribeiro Meneses CA, do Nascimento VA, de Souza VC, et al. Reemergence of Dengue Virus Serotype 3, Brazil, 2023. Emerg Infect Dis. 2023 Jul;29(7):1482–4. 7. 7.Santiago GA, Vergne E, Quiles Y, Cosme J, Vazquez J, Medina JF, et al. Analytical and Clinical Performance of the CDC Real Time RT-PCR Assay for Detection and Typing of Dengue Virus. PLoS Negl Trop Dis. 2013 Jul 11;7(7):e2311. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0002311&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23875046&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.16.24310413.atom) 8. 8.Vogels CBF, Hill V, Breban MI, Chaguza C, Paul LM, Sodeinde A, et al. DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus. BMC Genomics. 2024 May 1;25(1):1–16. 9. 9.Hill V, Cleemput S, Fonseca V, Tegally H, Brito AF, Gifford R, et al. A new lineage nomenclature to aid genomic surveillance of dengue virus [Internet]. medRxiv. 2024 [cited 2024 Jul 1]. p. 2024.05.16.24307504. Available from: [https://www.medrxiv.org/content/10.1101/2024.05.16.24307504v1.abstract](https://www.medrxiv.org/content/10.1101/2024.05.16.24307504v1.abstract) 10. 10.Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 2018 Dec 1;34(23):4121–3. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10/gdkbqx&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.16.24310413.atom) 11. 11.Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol Biol Evol. 2020 May 1;37(5):1530–4. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msaa015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32556291&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.16.24310413.atom) 12. 12.Carrillo-Hernandez MY, Ruiz-Saenz J, Jaimes-Villamizar L, Robledo-Restrepo SM, Martinez-Gutierrez M. Phylogenetic and evolutionary analysis of dengue virus serotypes circulating at the Colombian-Venezuelan border during 2015-2016 and 2018-2019. PLoS One. 2021 May 28;16(5):e0252379. 13. 13.Allicock OM, Lemey P, Tatem AJ, Pybus OG, Bennett SN, Mueller BA, et al. Phylogeography and population dynamics of dengue viruses in the Americas. Mol Biol Evol. 2012 Jun;29(6):1533–43. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msr320&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22319149&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.16.24310413.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000304200600006&link_type=ISI) 14. 14.García MP, Padilla C, Figueroa D, Manrique C, Cabezas C. Emergence of the Cosmopolitan genotype of dengue virus serotype 2 (DENV2) in Madre de Dios, Peru, 2019. Rev Peru Med Exp Salud Publica. 2022 Jun 24;39(1):126–8. 15. 15.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 Sep 15;206(6):893–901. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jis426&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22740715&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F18%2F2024.07.16.24310413.atom)