1. Introduction
According to the World Health Organization (WHO), Dengue virus (DENV) stands as the most relevant arbovirus in the world, posing a threat to 3.9 billion people in 120 countries [1]. Annually, over 390 million infections are estimated to occur, with tropical and subtropical areas being the most affected ones, particularly in urban and developing regions such as countries in South America and Asia [1]. The International Committee of Taxonomy of Viruses (ICTV) classifies DENV within the Flaviviridae family, Orthoflavivirus genus, and it is currently known as the Orthoflavivirus dengue species [2]. Like other members of the Orthoflavivirus genus, the DENV genome is a single-strand, positive-sense, ~11 kb long RNA, which encodes three structural proteins (which are known as Pre-membrane/membrane, Capsid and Envelope) and seven non-structural (NS) proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5). DENV is classified into four serotypes (DENV-1 to DENV-4) due to antigenic and genetic differences, with each serotype further exhibiting various genotypes [3]. DENV serotypes can be defined as viruses that present 65% similarity in the amino acid sequence, while genotypes refer to groups of DENV that diverge no more than 6% in the nucleotide sequence [3]. The four DENV serotypes entered the human population at distinct moments, with DENV-2 being the oldest [4]. Specifically, DENV-1 presents five distinct genotypes, which coalesced at distinct moments. For example, genotype V is the one that predominates in Brazil, and it was believed that this genotype originated in 1935 [5].
Infection with any serotype can lead to dengue fever, a disease whose clinical manifestation includes flu-like symptoms such as fever, nausea, vomiting, rashes and pain (more specifically headaches and retro-orbital pain). It can progress to dengue with warning signs, which usually manifests 2 days after the fever ceases, and it is characterized by the presence of blood in vomit or stool, abdominal pain and bleeding from mucous membranes, such as the nose and gums, and ultimately progresses to a life-threatening condition known as severe dengue, which is manifested by the presence of shock, internal bleeding and respiratory distress, which can lead to death [1]. It is also estimated that one out of four infections are symptomatic [1]. In Brazil, from 2013 to 2022, a substantial number of possible dengue cases amounting to 9,938,730, with 5.836 deaths, were reported [6,7]. During this period, there was a clear pattern of alternating circulating serotypes between DENV-1 and DENV-2 in the country, except for the year 2013, when DENV-4 was the main circulating serotype [6,7]. In 2023, DENV-3 re-emerged in the country, with three autochthonous cases detected in the State of Roraima and one imported case detected in Paraná [8]. It is also well established that the clade replacement phenomenon, which is defined as the alternance among DENV serotypes, genotypes and lineages in a given location, plays an important role in dengue fever severity and viral maintenance [6]. The impact of DENV clade replacement on the number of dengue cases and severity was described in many different regions, such as Brazil, Panamá, Vietnam and Thailand [9,10].
In this scenario, the municipality of Araraquara emerged as a significant contributor to the escalating number of dengue cases [11]. Up until 2005, only sporadic cases were reported in the municipality. However, the panorama of dengue fever underwent a significant change in 2006, when serotype 3 was introduced into the region, leading to the first reported outbreak [11,12]. Subsequently, in 2010, the introduction of serotypes 1 and 2 resulted in a significant rise in the number of cases and deaths [11]. A similar situation was reported in 2023 by Gularte JS et al., 2023. In this study, it was reported that more than 67,000 cases of dengue fever in Rio Grande do Sul State were registered in 2022, and genomic analyses performed by Gularte and colleagues confirmed the predominance of DENV-1 genotype V in the region [13].
Particularly noteworthy were the years 2015, 2019 and 2022 in Araraquara, which reported 8209, 23,538 and 21,070 cases, respectively [14]. During the period of 2014–2018, DENV-1 was the main serotype circulating in the municipality, followed by clade replacement in 2019, when DENV-2 triggered a massive outbreak [15]. In this study, we performed genome sequencing analysis of DENV using samples collected from the municipality of Araraquara in 2015, 2016, 2019, 2021 and 2022 to assess the molecular characteristics of the circulating viral strains during this period.
2. Materials and Methods
2.1. Ethical Statement
This study was approved by the ethical committee of the Medicine School of the University of São Paulo, Approval under No. 4.519.364, encompassing authorization for sample collection, viral detection and sequencing of viral genomes.
2.2. Study Population
The study population is part of a cohort of children and youth (2–16 years old) that has been under observation since 2014 from Araraquara. Briefly, this cohort was built to assess the immunological status of the youth population in the municipality, as well as to evaluate, annually, the number of dengue cases in the studied group. The cohort consisted of 3514 individuals, who had serum collected at the time of entry in the study. After initial sample collection for serological status assessment, blood samples were collected again and sent to the Institute of Tropical Medicina (IMT-USP) if symptoms befitting dengue fever were reported by the family of the participants (named “fever case” samples). The complete cohort study is better described elsewhere [15]. These samples were tested by qPCR (described in following subtopic) and/or NS1 antigen and serology and from 2014 to 2018, and a total of 257 qPCR DENV-1 cases were confirmed. In the present study, we randomly included 50 samples representative of the positivity detected by year.
2.3. DENV RNA Detection and Sample Selection for Next-Generation Sequencing (NGS)
Viral RNA was extracted from 140 µL of plasma and/or serum using the QIAamp Viral RNA MiniKit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions, obtaining a total of 60 µL of RNA. For the initial detection of viral RNA, a probe-based qPCR capable of detecting the 4 DENV serotypes was performed [16,17]. Following the confirmation of DENV infection, samples were typed through specific qPCR [17,18]. DENV-1 samples with a cycle threshold (CT) < 30 were eligible for the Next-Generation Sequencing (NGS).
2.4. DENV-1 Complete Genome Amplification
The primers to amplify the complete genome of DENV-1 were designed using the online tool PrimalScheme (available at primalscheme.com) and a polifasta file containing reference DENV-1 complete genomes as input. After defining the desired amplicon size of 400 base pairs and an overlap of 75 nucleotides, the tool generated two sets of primers suitable for implementation in a multiplex scheme, arranged into two pools. These primer sets enabled the amplification of the complete viral genome [19].
Following the cDNA synthesis utilizing First Strand cDNA ProtoScript II (New England Biolabs, Oxford, UK), the two multiplex PCR reactions were conducted under the conditions described elsewhere [20]. Amplicons were purified using 1× AMPure beads (Beckman Coulter, High Wycombe, UK) and quantified using the Qubit 4 fluorometer (Invitrogen, Carlsbad, CA, USA). A DNA concentration >2 ng DNA/µL per sample was considered suitable for whole-genome sequencing.
The libraries were generated using the genomic DNA sequencing kit SQK-LSK108, where 250 ng of barcoded PCR products were pooled in equimolar proportions. The libraries were loaded onto an Oxford Nanopore Flow Cell R9.4 and the data were collected after ~30 h [20].
2.5. Bioinformatics Workflow and Consensus Generation
The base-calling of the FAST5 reads containing the raw signal data was conducted using Guppy base-caller integrated within MinKNOW (
2.6. Envelope (E) Gene Sequencing
The DENV genomes that were not suitable for Next-Generation Sequencing (NGS) underwent partial envelope gene sequencing using the Sanger method. We opted for this approach to incorporate the highest possible number of samples from Araraquara into the analysis. To amplify ~800 pb of the envelope, we performed PCR using the primers and protocol described by Bruycker-Nogueira F et al., 2015, with some modifications [22]. Briefly, the PCR was performed using 14.75 µL of Ultrapure Nuclease-Free Water (NFW), 1× PCR Buffer, MgCl2 1.5 mM, 0.25 uM of each primer, 0.4 µM of dNTP and 1U of PlatinumTaq Polymerase in a final volume of 25 µL.
Amplicons were sequenced using BigDye Terminator v3.1 on the ABI 3730 DNA analyzer (Thermo Fisher Scientific, Waltham, MA, USA). Electropherograms were analyzed and the consensus sequences were assembled using the Geneious Prime software v.2023.2.1 (Dotmatics, Hertfordshire, UK).
2.7. Multiple Sequence Alignment and ML Phylogenetic Reconstruction
The dataset comprised 284 partial E gene sequences (~800 nucleotides long) retrieved from GenBank (sequence list available in the Supplementary Materials, Table S1) containing sampling date, host and country information and 37 sequences generated in this study (either by NGS or Sanger). Sequences were multiple aligned using MAFFT v7.407 and visually inspected using the SeaView v5.0.4 software [23,24]. To first classify DENV-1 sequences into genotypes, we reconstructed a Maximum Likelihood phylogenetic tree using the IQ-TREE 2 under the GTR + I + G nucleotide substitution model, as it was the best-fitting one selected by ModelFinder implemented in IQTREE [25]. The phylogenetic tree was visualized with FigTree v.1.4.4 [26].
2.8. Coalescent Analysis and Molecular Clock
The time to the most recent common ancestor (tMRCA) of the DENV-1 genotypes, and particularly of the Araraquara ones, was estimated using the Bayesian Markov Chain Monte Carlo implemented in BEAST 1.10.4 [27,28]. The Bayesian Skyline (BSL) coalescent method was performed with a relaxed uncorrelated exponential molecular clock using time-stamped data scaled in years using the GTR + I + G nucleotide substitution model. The convergence of parameters was inspected with Tracer v.1.7.2 with uncertainties addressed as 95% highest probability density (HPD) intervals.
After 100 million runs, the trees sampled at every 10,000 steps were summarized in a maximum clade credibility (MCC) tree with 10% burn-in using TreeAnnotator and visualized with FigTree v.1.4.4. In addition to the BSL method, the exponential growth tree prior was run for comparison.
3. Results
A total of 37 samples were positive on the qPCR for DENV-1 with CT < 30, median 24.5 (18.4–30.7), and were submitted to the whole-genome amplification. The number of samples per year was as follows: 2015 (n = 26, collected from March to May), 2016 (n = 1, collected in April), 2019 (n = 2, from February), 2021 (n = 2, collected in June) and 2022 (n = 6, collected in May). Twenty-eight out of the thirty-seven samples submitted to the whole-genome amplification had a sufficient DNA concentration and were sequenced using the Nanopore method. The remaining samples (n = 9) were submitted to Sanger sequencing (envelope gene). Maximum likelihood analysis of DENV-1 showed a clear separation of the distinct genotypes and the genotype V division into lineage I and lineage II, previously classified by Bruycker-Nogueira and colleagues [9] (results not shown). The genomes sequenced in this study can be found in the Genbank public database under the accession numbers OR553381 → OR553391, OR545518 → OR545534 and OR520604.
We then ran the same dataset in BEAST to perform the temporal estimates. According to the BLS estimates, DENV-1 coalesced approximately 138 years ago (95% HPD 102–188 ya) in 1884, while the origin of the genotype V was estimated to be around 57 years ago (95% HPD 50–67 ya) in 1965. Similar data were found under the exponential growth.
The MCC tree also showed the two different lineages of genotype V [9]. The genomes from 2015 and 2016 were classified as lineage II (LII), while genomes collected from 2019 onwards were classified as lineage I (LI). According to Bayesian estimates, lineage II coalesced in 1982 (around 41 years ago), lineage I coalesced in early 1977 (around 46 years ago) and lineages Ia and Ib both coalesced in 2001 (Figure 1).
All but five viruses from Araraquara 2015/2016 clustered together in a branch that probably emerged in 2009 (95% HPD 2006–2011). Two out of the remaining samples from 2015 clustered with samples from Central (2013) and Southeast (2013) regions. This second cluster consists primarily of Brazilian DENV strains, comprising viruses from Goiás and Rio de Janeiro (2013), Ceará and Sergipe (2014) and São Paulo (2016). Notably, these two clusters constitute the entire lineage II. Viruses from Araraquara, 2019, belong to lineage Ib and are grouped with other viruses from Southeast Brazil (2016 and 2018).
The cluster also includes viruses from the North (2010) and Northeast (2014) regions of the country. The tMRCA for this group was estimated at 16.5 years ago, in 2006 (95% HPD 2004–2008). It is also worth noting that, in contrast to lineage II, some of these genomes were closer to other South American countries, such as Venezuela and Colombia. The DENV-1 sequences obtained in 2021 and 2022 clustered together and were close to other genomes from the Southeast (2016, 2021, 2022 and 2023) and Central regions (2016) of the country. The node comprising these sequences is dated to 9.5 years ago (95% HPD 8–10 ya), in 2013.
Finally, the dates of key DENV-1 events were re-estimated. According to our results, DENV-1 had its origin 138 years ago, estimated in 1884 (95% HPD 1833–1920), and genotype V arose 58 years ago, in 1965 (95% HPD 1955–1972). The approximated origin of genotype V lineages and sublineages as well as the tMRCA of the other four genotypes is described in Table 1.
4. Discussion
The phylogenetic structure of DENV in endemic regions is marked by lineage turnovers, which cause significant impacts on the epidemiology of dengue due to the susceptibilities to cross-reactive immune responses against the currently circulating virus. Brazil annually reports the highest number of dengue fever cases and deaths, and clade replacement is believed to play a great part in these numbers; therefore, constant vigilance for the introduction of new clades in endemic regions is of utmost importance [6,9,10,29,30,31,32,33].
Araraquara, located in the central region of São Paulo state, is among the most industrialized cities in the state and reported the first autochthonous dengue case in 1996. However, information regarding the serotype was not available until 2006, when DENV-3 was described in the region. Since then, the number of cases progressively rose, leading to severe outbreaks and epidemics every year [9,11]. After the emergence of DENV-1 in 2010, Araraquara experienced an extended predominance of this serotype for several years with a remarkable number of reported cases over the years [11]. Particularly, from 2014 to 2018, DENV-1 emerged as the predominant serotype in Araraquara, while DENV-2 dominated from late 2018 to 2019, leading to more than 20 thousand reported cases [11,14,15].
Nonetheless, in 2020, DENV-1 reclaimed its status as the prevailing serotype and has continued to do so up to the present time [34]. The phylogenetic analysis using longitudinal DENV-1 samples collected in an 8-year time span performed in this study represents the first description of DENV-1 lineages circulating in the municipality of Araraquara. During this period, we detected the presence of the two lineages of DENV-1 genotype V in Araraquara, L I and L II [29]. Initially, DENV-1 was represented exclusively by lineage II and almost all viruses clustered together. Only three samples (two from 2015 and one from 2016) did not cluster with the remaining viruses collected in 2015, though we cannot definitively determine whether this is due to a distinct introduction of a very similar strain or a result of locally accumulated diversity. From 2019 onwards, viruses from L II were no longer sampled, and all circulating viruses until 2022 belonged to lineage I. Particularly, viruses from 2019 clustered within sub-clade L Ib, while those from 2021 and 2022 belonged to the sub-clade L Ia, clearly revealing a second clade replacement. Indeed, the most recent spread of lineage I across the country is also supported by a study by Gularte et al., 2023, in which the phylogenetic analyses of 70 DENV-1 genomes from 2022 clustered within lineage I [13].
The presence of distinct lineages and sublineages of DENV-1 in Araraquara might have played a critical role in sustaining the prevalence of this serotype in the region for such an extended period, as it has been seen elsewhere [32]. Unfortunately, we do not have samples prior to 2015 and also no samples were collected in 2017 and 2018, so it is impossible to precisely determine how many clade/lineage replacements occurred since DENV-1’s emergence in Araraquara. Nevertheless, the continual emergence of distinct lineages might be the key to the observed time frame of DENV-1 predominance. The isolates detected in 2019 (lineage Ib) were not able to spread, maybe due to the dominance of DENV-2 at that time. But two years later, the number of DENV-1-susceptible hosts was sufficient to allow the DENV-1 lineage Ia to spread and replace DENV-2, causing again a major outbreak.
The phylogeny also evidenced two distinct clusters within lineage I, sufficiently divergent from the previously described L Ia and L Ib [9], which could be named L Ic and Id here, but the low posterior probability did not allow us to properly propose these classifications. Although these clusters also comprise Brazilian viruses, no sample from Araraquara was classified as these new sub-lineages.
While several studies point out that clade replacement is a stochastic event and does not depend on selection, there is also evidence that strains of DENV can exhibit different fitnesses either in mosquitoes or in vertebrate hosts [10,35,36,37]. For instance, during a clade replacement in Indonesia within the Cosmopolitan genotype of DENV-2, the increased number of cases was associated with a higher fitness of the new clade in mammalian cells [38]. Another relevant topic of clade replacement was demonstrated in a study by Katzelnick et al., 2021. It was shown that within one serotype, genotypes can also present significant divergence at the antigenic level, with one genotype being more antigenically related to sequences from other serotypes than to sequences from the same serotype. For instance, DENV-3 genotype II is much more similar to other serotypes than to DENV-3 genotype III. These findings support the hypothesis that within one serotype, genetic and antigenic divergences play important roles in serotype maintenance and number of cases/severe cases [39].
Our group recently reported a clade replacement within the DENV-2 American/Asian genotype (III) in São Paulo, where the introduction of the BR-4 lineage was associated with major outbreaks in 2019 and rapid expansion of the lineage in the country [31]. Although we have not evaluated the number of deaths and severe dengue cases in the 2015 and 2022 DENV-1 epidemics, the higher number of cases caused by lineage I compared to lineage II was noticeable (i.e., 8209 in 2015 and 21,070 reported cases in 2022).
Our Bayesian analysis also enabled us to re-estimate the dates of key DENV-1 events. According to the estimates, the origins of serotype 1 as well as the genotypes I to IV were similar to those reported elsewhere [5]. However, we noticed that the clustering pattern of two sequences (AF425620, from Cote d’Ivoire and AF425625, from Nigeria), previously classified as genotype V [40], were not consistent and grouped near genotype III in the ML (not shown) and MCC trees. The previous reports that considered these sequences to be genotype V dated the origin of this genotype back to 1935 [34,35]. However, according to our estimation, genotype V might have emerged much more recently, around 1968. It has to be noted that both previous studies included a limited number of genotype III reference sequences [5,40]. Consequently, the lower genetic diversity of genotypes III and V from the previous phylogenetic analysis might have hindered accurate classification during the reconstructions. In an effort to explore the potential factors underlying the discrepancies observed in previous phylogenetic classifications, we tested these samples for recombination and estimated the genetic distances between them and others from genotypes III and V. RDP [41] found no evidence of recombination in both sequences. Genetic distance analysis using the Kimura 2 parameter revealed >6% distance between these genomes and those from genotype V, but a 5.2% distance from viruses from genotype III.
Based on previous estimates that classify DENV as distinct genotypes when the genomes differ by around 6–8% at the nucleotide level [3,42], one might suggest that these viruses might be more appropriately classified as genotype III. Finally, the Dengue genotyping tool (available at
We did not go further on the analysis of these two genomes but it is likely they represent either atypical genotype III viruses or sylvatic strains, so we opted to include them as genotype III in our tMRCA estimates.
5. Conclusions
In summary, this investigation has unveiled the presence of the two major lineages of DENV-1 genotype V in Araraquara, SP, and described two clade replacements within this genotype during 2015–2022. These lineages exhibit a history of propagation across various regions of the country before emerging in the municipality. The implications are suggestive that these lineages might have influenced the divergent case counts observed during each epidemic occurrence and played a critical role in the maintenance of DENV-1 circulation. These results emphasize the importance of genomic surveillance as an additional tool in understanding the relationship between viral lineages and the magnitude of outbreaks. Furthermore, our findings have yielded novel insights into the origins of DENV-1 genotype V—a strain of utmost significance within the Brazilian context.
Conceptualization: C.S.d.S. and C.M.R.; Methodology: G.S.C., I.M.C., J.G.d.J., T.M.C., C.A.M.d.S., A.C.F., A.V.d.P., E.J.d.A.L., E.C.S. and C.M.R.; Formal analysis: C.S.d.S. and C.M.R.; Investigation: G.S.C., I.M.C., J.G.d.J., T.M.C., C.A.M.d.S., Â.A.C., M.I., A.C.R., A.C.F., A.V.d.P., W.M.F., E.J.d.A.L., E.C.S. and C.M.R.; Resources: E.C.S. and C.M.R.; Writing—original draft, C.S.d.S. and C.M.R.; Writing—review & editing, C.S.d.S. and C.M.R.; Supervision: C.M.R. All authors have read and agreed to the published version of the manuscript.
This study was conducted in accordance with the Declaration of Helsinki, and approved by approved by the ethical committee of the Medicine School of the University of São Paulo, Approval # 4.519.364.
Informed consent was obtained from all subjects involved in the study.
Data are contained within the article and
The authors declare no conflict of interest.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Maximum clade credibility tree of DENV-1. The genotypes I to V are identified by colorful internal blocks. Genotype V, which includes most of the reference sequences and all Brazilian DENV-1, is also classified into lineages and sub-lineages. Small black dots denote the most important nodes with posterior probability >0.7. The tree is scaled by years, represented by vertical line denoting the present (time zero) to 150 years before the present. Sequences from Araraquara were colored according to the year of sampling, 2015–2022 (colorful circles).
DENV-1 and genotypes’ tMRCAs estimated under the Bayesian Skyline Coalescent and Exponential growth tree priors. The tMRCA and 95% High-Probability Densities are shown for both models, genotypes and lineages. The tMRCA is presented in years.
| Group | tMRCA Bayesian Skyline (HPD 95%) | tMRCA Exponential Growth (HPD 95%) |
|---|---|---|
| DENV-1 | 1884 (1833–1920) | 1886 (1842–1916) |
| Genotype I | 1966 (1953–1974) | 1967 (1954–1976) |
| Genotype II | 1931 (1907–1948) | 1930 (1909–1947) |
| Genotype III | 1945 (1927–1958) | 1943 (1926–1955) |
| Genotype IV | 1944 (1921–1959) | 1944 (1923–1959) |
| Genotype V | 1965 (1955–1972) | 1966 (1957–1973) |
| Lineage I | 1977 (1971–1978) | 1976 (1973–1978) |
| Lineage Ia | 2001 (1998–2005) | 2002 (1999–2005) |
| Lineage Ib | 2001 (1998–2005) | 2002 (2000–2005) |
| Lineage II | 1982 (1977–1985) | 1976 (1968–1982) |
| Araraquara 2015 | 2009 (2006–2011) | 2008 (2007–2011) |
Supplementary Materials
The following supporting information can be downloaded at
References
1. World Health Organization. Dengue and Severe Dengue. Available online: https://www.who.int/news-room/fact-sheets/detail/dengue-and-severe-dengue (accessed on 17 March 2023).
2. Lefkowitz, E.J.; Dempsey, D.M.; Hendrickson, R.C.; Orton, R.J.; Siddell, S.G.; Smith, D.B. Virus taxonomy: The database of the International Committee on Taxonomy of Viruses (ICTV). Nucleic Acids Res.; 2018; 46, pp. D708-D717. [DOI: https://dx.doi.org/10.1093/nar/gkx932] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29040670]
3. Holmes, E.; Twiddy, S. The origin, emergence and evolutionary genetics of dengue virus. Infect. Genet. Evol.; 2003; 3, pp. 19-28. [DOI: https://dx.doi.org/10.1016/S1567-1348(03)00004-2] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12797969]
4. Zanotto, P.M.; Gould, E.A.; Gao, G.F.; Harvey, P.H.; Holmes, E.C. Population dynamics of flaviviruses revealed by molecular phylogenies. Proc. Natl. Acad. Sci. USA; 1996; 93, pp. 548-553. [DOI: https://dx.doi.org/10.1073/pnas.93.2.548] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8570593]
5. Villabona-Arenas, C.J.; de Andrade Zanotto, P.M. Worldwide Spread of Dengue Virus Type 1. PLoS ONE; 2013; 8, e62649. [DOI: https://dx.doi.org/10.1371/journal.pone.0062649] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23675416]
6. Ministério da Saúde. Óbito por Arboviroses no Brasil, 2008 a 2019. Available online: http://plataforma.saude.gov.br/anomalias-congenitas/boletim-epidemiologico-SVS-33-2020.pdf (accessed on 25 August 2020).
7. Souza, C.S.; Romano, C.M. Dengue in the cooling off period of the COVID-19 epidemic in Brazil: From the shadows to the spotlight. Rev. Inst. Med. Trop. São Paulo; 2022; 64, e44. [DOI: https://dx.doi.org/10.1590/s1678-9946202264044] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35730870]
8. Naveca, F.G.; Santiago, G.A.; Maito, R.M.; Meneses, C.A.R.; Nascimento, V.A.; Souza, V.C.; do Nascimento, F.O.; Silva, D.; Mejía, M.; Gonçalves, L. et al. Reemergence of Dengue Virus Serotype 3, Brazil, 2023. medRxiv; 2023; [DOI: https://dx.doi.org/10.3201/eid2907.230595] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/37347884]
9. de Bruycker-Nogueira, F.; Souza, T.M.A.; Chouin-Carneiro, T.; da Costa Faria, N.R.; Santos, J.B.; Torres, M.C.; Ramalho, I.L.C.; de Aguiar, S.F.; Nogueira, R.M.R.; de Filippis, A.M.B. et al. DENV-1 Genotype V in Brazil: Spatiotemporal dispersion pattern reveals continuous co-circulation of distinct lineages until 2016. Sci. Rep.; 2018; 8, 17160. [DOI: https://dx.doi.org/10.1038/s41598-018-35622-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30464188]
10. Díaz, Y.; Chen-Germán, M.; Quiroz, E.; Carrera, J.-P.; Cisneros, J.; Moreno, B.; Cerezo, L.; Martinez-Torres, A.O.; Moreno, L.; de Mosca, I.B. et al. Molecular Epidemiology of Dengue in Panama: 25 Years of Circulation. Viruses; 2019; 11, 764. [DOI: https://dx.doi.org/10.3390/v11080764]
11. Ferreira, A.C.; Chiaravalloti Neto, F.; Mondini, A. Dengue in Araraquara, state of São Paulo: Epidemiology, climate and Aedes aegypti infestation. Rev. Saúde Pública; 2018; 52, 18. [DOI: https://dx.doi.org/10.11606/S1518-8787.2018052000414]
12. Centro de Vigilância Epidemiológica “Prof. Alexandre Vranjac”. Dados Estatísticos. São Paulo: CVE. 2015; Available online: http://portal.saude.sp.gov.br/cve-centro-devigilancia-epidemiologica-prof.-alexandre-vranjac/areas-de-vigilancia/doencas-de-transmissaopor-vetores-e-zoonoses/agravos/dengue/dados-estatisticos (accessed on 15 September 2016).
13. Gularte, J.S.; Sacchetto, L.; Demoliner, M.; Girardi, V.; da Silva, M.S.; Filippi, M.; Pereira, V.M.d.A.G.; Hansen, A.W.; da Silva, L.L.; Fleck, J.D. et al. DENV-1 genotype V linked to the 2022 dengue epidemic in Southern Brazil. J. Clin. Virol.; 2023; 168, 105599. [DOI: https://dx.doi.org/10.1016/j.jcv.2023.105599]
14. BRASILMinistério da SaúdeSecretaria de Vigilância em Saúde. Sistema de Informação de Agravos de Notificação (SINAN). Available online: http://www.portalsinan.saude.gov.br/sinan-dengue-chikungunya (accessed on 25 January 2023).
15. Luna, E.J.; Figueiredo, G.M.; Levi, J.E.; Campos, S.R.; Felix, A.C.; e Souza, N.S.; Figueiredo, W.M.; Costa, A.A.; Cardoso, M.R.; Pannuti, C.S. A cohort study to assess the incidence of dengue, Brazil, 2014–2018. Acta Trop.; 2020; 204, 105313. [DOI: https://dx.doi.org/10.1016/j.actatropica.2019.105313] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31863736]
16. Huhtamo, E.; Hasu, E.; Uzcátegui, N.Y.; Erra, E.; Nikkari, S.; Kantele, A.; Vapalahti, O.; Piiparinen, H. Early diagnosis of dengue in travelers: Comparison of a novel real-time RT-PCR, NS1 antigen detection and serology. J. Clin. Virol.; 2010; 47, pp. 49-53. [DOI: https://dx.doi.org/10.1016/j.jcv.2009.11.001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19963435]
17. Luna, E.; Figueiredo, G.; Levi, J.; Campos, S.; Felix, A.; Souza, N.; Figueiredo, W.; Costa, A.; Cardoso, M.; Pannuti, C. Data on dengue incidence in South-eastern Brazil, 2014–2018. Data Brief; 2020; 29, 105266. [DOI: https://dx.doi.org/10.1016/j.dib.2020.105266] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32095494]
18. Callahan, J.D.; Wu, S.-J.L.; Dion-Schultz, A.; Mangold, B.E.; Peruski, L.F.; Watts, D.M.; Porter, K.R.; Murphy, G.R.; Suharyono, W.; King, C.-C. et al. Development and evaluation of serotype- and group-specific fluorogenic reverse transcriptase PCR (TaqMan) assays for dengue virus. J. Clin. Microbiol.; 2001; 39, pp. 4119-4124. [DOI: https://dx.doi.org/10.1128/JCM.39.11.4119-4124.2001] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11682539]
19. Claro, I.M. Em tempo Real, Rápida Detecção e Sequenciamento de Arbovírus No Brasil. Ph.D. Thesis; Universidade de São Paulo: São Paulo, Brazil, 2021; [DOI: https://dx.doi.org/10.11606/T.5.2021.tde-07012022-140752]
20. Quick, J.; Grubaugh, N.D.; Pullan, S.T.; Claro, I.M.; Smith, A.D.; Gangavarapu, K.; Oliveira, G.; Robles-Sikisaka, R.; Rogers, T.F.; A Beutler, N. et al. Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples. Nat. Protoc.; 2017; 12, pp. 1261-1276. [DOI: https://dx.doi.org/10.1038/nprot.2017.066] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28538739]
21. Liu, C.H.; Di, Y.P. Analysis of RNA Sequencing Data Using CLC Genomics Workbench. Methods Mol. Biol.; 2020; 2102, pp. 61-113. [DOI: https://dx.doi.org/10.1007/978-1-0716-0223-2_4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31989550]
22. de Bruycker-Nogueira, F.; Nogueira, R.M.R.; Faria, N.R.d.C.; Simões, J.B.S.; Nunes, P.C.G.; de Filippis, A.M.B.; dos Santos, F.B. Insights of the genetic diversity of DENV-1 detected in Brazil in 25years: Analysis of the envelope domain III allows lineages characterization. Infect. Genet. Evol.; 2015; 34, pp. 126-136. [DOI: https://dx.doi.org/10.1016/j.meegid.2015.07.007]
23. Katoh, K.; Standley, D.M. MAFFT Multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol.; 2013; 30, pp. 772-780. [DOI: https://dx.doi.org/10.1093/molbev/mst010]
24. Gouy, M.; Guindon, S.; Gascuel, O. Seaview version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol.; 2009; 27, pp. 221-224. [DOI: https://dx.doi.org/10.1093/molbev/msp259]
25. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; von Haeseler, A.; Lanfear, R. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol.; 2020; 37, pp. 1530-1534. [DOI: https://dx.doi.org/10.1093/molbev/msaa015]
26. Rambaut, A. FigTree v1.3.1. Institute of Evolutionary Biology, University of Edinburgh, Edinburgh. 2010; Available online: http://tree.bio.ed.ac.uk/software/figtree/ (accessed on 20 July 2020).
27. Ferreira, M.A.R.; Suchard, M.A. Bayesian analysis of elapsed times in continuous-time Markov chains. Can. J. Stat.; 2008; 36, pp. 355-368. [DOI: https://dx.doi.org/10.1002/cjs.5550360302]
28. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol.; 2012; 29, pp. 1969-1973. [DOI: https://dx.doi.org/10.1093/molbev/mss075] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22367748]
29. Teoh, B.-T.; Sam, S.-S.; Tan, K.-K.; Johari, J.; Shu, M.-H.; Danlami, M.B.; Abd-Jamil, J.; MatRahim, N.; Mahadi, N.M.; AbuBakar, S. Dengue virus type 1 clade replacement in recurring homotypic outbreaks. BMC Evol. Biol.; 2013; 13, 213. [DOI: https://dx.doi.org/10.1186/1471-2148-13-213] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24073945]
30. European Center for Disease Control. Dengue Worldwide Overview. Available online: https://www.ecdc.europa.eu/en/dengue-monthly (accessed on 23 August 2023).
31. de Jesus, J.G.; Dutra, K.R.; da Silva Sales, F.C.; Claro, I.M.; Terzian, A.C.; da Silva Candido, D.; Hill, S.C.; Thézé, J.; Torres, C.; D’agostini, T.L. et al. Genomic detection of a virus lineage replacement event of dengue virus serotype 2 in Brazil, 2019. Mem. Inst. Oswaldo Cruz; 2020; 115, e190423. [DOI: https://dx.doi.org/10.1590/0074-02760190423] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32428189]
32. Nunes, P.C.G.; Sampaio, S.A.F.; da Costa, N.R.; de Mendonça, M.C.L.; da Rocha Queiroz Lima, M.; Araujo, S.E.M.; dos Santos, F.B.; Simões, J.B.S.; de Santis Gonçalves, B.; Nogueira, R.M.R. et al. Dengue severity associated with age and a new lineage of dengue virus-type 2 during an outbreak in Rio De Janeiro, Brazil. J. Med. Virol.; 2016; 88, pp. 1130-1136. [DOI: https://dx.doi.org/10.1002/jmv.24464] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27061403]
33. de Souza, U.J.B.; Macedo, Y.d.S.M.; dos Santos, R.N.; Cardoso, F.D.P.; Galvão, J.D.; Gabev, E.E.; Franco, A.C.; Roehe, P.M.; Spilki, F.R.; Campos, F.S. Circulation of Dengue Virus Serotype 1 Genotype V and Dengue Virus Serotype 2 Genotype III in Tocantins State, Northern Brazil, 2021–2022. Viruses; 2023; 15, 2136. [DOI: https://dx.doi.org/10.3390/v15112136] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/38005815]
34. Ministério da Saúde. Monitoramento dos Casos de Arboviroses Até a Semana Epidemiológica 51 de 2022. Available online: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/epidemiologicos/edicoes/2022/boletim-epidemiologico-vol-53-no48/view (accessed on 29 December 2022).
35. Thu, H.M.; Lowry, K.; Jiang, L.; Hlaing, T.; Holmes, E.C.; Aaskov, J. Lineage extinction and replacement in dengue type 1 virus populations are due to stochastic events rather than to natural selection. Virology; 2005; 336, pp. 163-172. [DOI: https://dx.doi.org/10.1016/j.virol.2005.03.018]
36. Armstrong, P.M.; Rico-Hesse, R. Efficiency of dengue serotype 2 virus strains to infect and disseminate in Aedes aegypti. Am. J. Trop. Med. Hyg.; 2003; 68, pp. 539-544. [DOI: https://dx.doi.org/10.4269/ajtmh.2003.68.539]
37. Cologna, R.; Rico-Hesse, R. American genotype structures decrease dengue virus output from human monocytes and dendritic cells. J. Virol.; 2003; 77, pp. 3929-3938. [DOI: https://dx.doi.org/10.1128/JVI.77.7.3929-3938.2003]
38. Tan, C.H.; Hapuarachchi, H.C.; Tan, L.K.; Wong, P.S.J.; Li, M.Z.I.; Wong, W.Y.; Ng, L.C. Lineage Replacement Associated with Fitness Gain in Mammalian Cells and Aedes aegypti: A Catalyst for Dengue Virus Type 2 Transmission. Microorganisms; 2022; 10, 1100. [DOI: https://dx.doi.org/10.3390/microorganisms10061100]
39. Katzelnick, L.C.; Escoto, A.C.; Huang, A.T.; Garcia-Carreras, B.; Chowdhury, N.; Berry, I.M.; Chavez, C.; Buchy, P.; Duong, V.; Dussart, P. et al. Antigenic evolution of dengue viruses over 20 years. Science; 2021; 374, pp. 999-1004. [DOI: https://dx.doi.org/10.1126/science.abk0058] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34793238]
40. Goncalvez, A.P.; Escalante, A.A.; Pujol, F.H.; Ludert, J.E.; Tovar, D.; Salas, R.A.; Liprandi, F. Diversity and evolution of the envelope gene of dengue virus type 1. Virology; 2002; 303, pp. 110-119. [DOI: https://dx.doi.org/10.1006/viro.2002.1686] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12482662]
41. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol.; 2015; 1, vev003. [DOI: https://dx.doi.org/10.1093/ve/vev003] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27774277]
42. Rico-Hesse, R. Molecular evolution and distribution of dengue viruses type 1 and 2 in nature. Virology; 1990; 174, pp. 479-493. [DOI: https://dx.doi.org/10.1016/0042-6822(90)90102-W]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Thaís Moura Coletti 5 ; Camila Alves Maia da Silva 5
; Ângela Aparecida Costa 6 ; Inenami, Marta 6 ; Ribeiro, Andreia C 6 ; Felix, Alvina Clara 1 ; Anderson Vicente de Paula 1 ; Figueiredo, Walter M 6 ; Expedito José de Albuquerque Luna 7 ; Sabino, Ester C 5 ; Romano, Camila M 8 1 Laboratório de Virologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
2 Laboratório de Virologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
3 Laboratório de Parasitologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
4 Laboratório de Parasitologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
5 Laboratório de Parasitologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
6 Serviço Especial de Saúde de Araraquara-SESA, Faculdade de Saúde Pública da USP, São Paulo 01246-904, SP, Brazil;
7 Departamento de Medicina Preventiva/Instituto de Medicina Tropical, Faculdade de Medicina, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;
8 Laboratório de Virologia, Instituto de Medicina Tropical, Universidade de São Paulo, São Paulo 05403-000, SP, Brazil;