Epidemiological and Genetic Evolutionary Characteristics of Dengue Virus in Hainan Province,China,Based on Whole-Genome Sequencing and Serological Data

Abstract

Background:

Dengue fever is one of the most widely distributed vector-borne infectious diseases globally, prevalent in tropical and subtropical regions. An estimated 3.9 billion people in 128 countries are at risk of infection globally, and the 2024 outbreak was the most severe. As a tropical region, Hainan Province in China serves as a significant endemic area for dengue fever, where epidemic prevention and control remain critical.

Objective:

To analyze the transmission dynamics and pathogen variation characteristics of the 2024 dengue fever outbreak in Hainan Province.

Methods:

Acute-phase serum samples from 18 confirmed dengue cases during the 2024 outbreak in Hainan Province were collected and serotyped using TaqMan real-time PCR. The viral envelope (E) gene and representative whole genomes were sequenced and compared with domestic reference strains to assess epidemic trends. Phylogenetic and molecular clock analyses were performed based on E gene and full-genome sequences to infer spatiotemporal transmission dynamics. Sequence alignment and homology modeling were used to identify E gene mutations and evaluate their potential effects on E protein functional domains and viral pathogenicity.

Results:

Among the 18 samples, 17 were dengue virus (DENV)-1, and 1 was DENV-3. Phylogenetic analysis revealed that the primary strains in this outbreak were highly homologous to those circulating in Guangdong Province; however, two strains exhibited significant genetic differences, suggesting possible independent introduction. Additionally, two unique DENV-1 E gene mutations (D147N and S338L) were identified, which we hypothesize may be associated with viral characteristics, subject to further functional verification.

Conclusion:

This study systematically reveals the epidemiological characteristics of the 2024 dengue fever outbreak and provides a scientific basis for formulating local control strategies and a research foundation for further exploring the molecular-level pathogenicity and transmissibility of DENV.

Keywords

dengue fever 2024 dengue outbreak E gene Hainan Province (China)phylogenetic analysis

Introduction

Dengue fever is an acute, vector-borne infectious disease caused by the dengue virus (DENV), primarily transmitted by Aedes aegypti and Aedes albopictus mosquitoes. It is endemic in tropical and subtropical regions, affecting approximately 3.9 billion people in 128 countries (Cattarino et al., 2020). Common symptoms include fever, headache, rash, muscle and joint pain, nausea, and vomiting (Katzelnick et al., 2017). Severe or secondary infections with different serotypes can result in dengue hemorrhagic fever or shock syndrome (Zompi and Harris, 2013). Currently, no antiviral drugs are approved for treatment, making prevention the primary approach to reducing disease incidence (Akshatha et al., 2021).

DENV is a spherical, single-stranded positive-sense RNA virus belonging to the family Flaviviridae and genus Flavivirus (Diamond and Pierson, 2015). Its ∼10.7 kb genome contains a single open reading frame encoding three structural proteins—capsid, envelope (E), and precursor membrane—and seven nonstructural proteins (Iglesias and Gamarnik, 2011). The E protein, comprising approximately 495 amino acids, plays a central role in the viral life cycle by mediating receptor recognition and binding (Kuhn et al., 2002) and facilitating fusion between the viral E and host cell membrane (Anasir et al., 2020; Zhang et al., 2004). As the primary target of neutralizing antibodies, mutations in the E protein can influence viral infectivity, immune escape, and transmission efficiency (Sessions et al., 2015). Notably, amino acid substitutions at positions 324, 351, and 380 in domain III of the DENV-1 E protein have been shown to enhance viral infectivity, stability, and immune evasion in mammals, thereby promoting the global spread of this genotype (Jiang et al., 2024). Consequently, in-depth investigation of E protein mutations is critical for elucidating viral evolution, assessing epidemic risk, and guiding vaccine and antibody therapeutic development (Ali et al., 2025).

Based on antigenic differences in the E gene, DENV can be classified into four serotypes (DENV-1 to DENV-4), with amino acid sequence differences ranging from 25% to 40% (World Health Organization, 2009). Each serotype can be further subdivided into multiple genotypes based on genetic variations. For example, DENV-1 includes genotypes I–V, with genotype I further subdivided into lineages I–IV (Du et al., 2021). However, long-lasting cross-immune protection between different serotypes is lacking, resulting in individuals acquiring long-term immunity to one serotype after infection but remaining susceptible to others, which may even lead to more severe disease (Katzelnick et al., 2017). Therefore, accurately distinguishing DENV serotypes and genotypes is critically important for epidemic surveillance, viral tracing, vaccine development, and clinical warnings for severe cases (Katzelnick et al., 2015).

In 2024, the global dengue fever epidemic exhibited unprecedented growth, with over 12 million cumulative reported cases and more than 8000 deaths, representing a 230% increase compared to the same period in 2023 (Venkatesan, 2024), surpassing all historical records. Against this backdrop, Hainan Province, due to its frequent international exchanges, has become a typical high-incidence area for dengue fever in China’s tropical regions and a key focus for control strategies. In 2024, local dengue fever cases emerged successively in various hospitals in the region, yet understanding of the viral genomic information remains highly limited. Analyzing the viral genomic information from local cases in Hainan Province will enable the precise determination of the sources, types, and variations of circulating strains, providing a scientific basis for formulating regional control strategies and effectively addressing the dual challenges of international importation and local transmission. This study conducted viral amplification and sequence analysis on 18 serum samples from the Haikou City Center for Disease Control and Prevention, constructing phylogenetic trees and performing molecular clock analyses to examine the epidemiological characteristics and transmission dynamics of dengue fever in Hainan Province, thereby offering scientific evidence and recommendations for subsequent control strategies.

Materials and Methods

Sample collection and processing

Between October 11 and November 13, 2024, 18 serum samples were collected from patients with dengue-like symptoms at sentinel hospitals and community health service centers in 4 administrative districts of Haikou City (Longhua, Xiuying, Meilan, and Qiongshan). Patients were included if they tested positive for DENV RNA by RT-PCR with a cycle threshold (Ct) value < 25. The cohort comprised 12 males and 6 females aged 21–55 years. Of the sequenced strains, 17 were DENV-1 and 1 was DENV-3, with detailed demographic and epidemiological information provided in Supplementary Table S1. Whole blood samples were collected in tubes without anticoagulants, centrifuged, and the separated serum was transferred to sterile 2 mL tubes. Initial screening was performed using immunochromatography and double-antibody sandwich assays, with final confirmation by molecular detection. Serum samples were transported under cold-chain conditions, delivered to the laboratory within 24 h, and stored at −80°C until further analysis.

Detection and serotyping of DENV

Viral RNA was extracted from serum samples using the Viral RNA Mini Kit (Qiagen). DENV diagnosis and serotyping were performed using an RT-Quantitative Polymerase Chain Reaction (qPCR) kit (Takara), with primers and probes designed according to the national health standard (WS 216–2018). RT-qPCR was conducted in a total reaction volume of 20 µL, containing premix, primers, probes, and RNase-free water. Primers specific for the four DENV serotypes were combined with probes labeled with different fluorophores (Supplementary Table S2). The PCR cycling conditions were as follows: reverse transcription at 52°C for 5 min, 95°C for 10 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s.

Isolation, confirmation, and amplification of the DENV E gene

A total of 500 µL of C6/36 cells at a density of 1 × 10⁵/mL were seeded into a 6-well plate and cultured in a 28°C, 5% CO₂ incubator for 1 day until a monolayer was formed. The medium was aspirated, and the cells were washed twice with phosphate-buffered saline (Gibco). Subsequently, a 20-fold diluted patient serum was added, with continued incubation at 28°C, 5% CO₂, and daily observation of cytopathic effects (CPE). After 7–10 days of culture, significant CPEs were observed, and the supernatant from the affected cells was harvested and stored in a −80°C freezer for subsequent experiments. The viral supernatant was then diluted 10-fold to infect C6/36 cells, and after 24 h, immunofluorescence detection was performed. Viral nucleic acids were extracted, and one-step reverse transcription-PCR was used to specifically amplify the E genes of DENV-1 and DENV-3, with primers as shown in Supplementary Table S3. Nucleic acid electrophoresis was conducted using 1% agarose gel to confirm the correct amplification bands (DENV-1: 1740 bp, DENV-3: 2335 bp) (Fraenkel et al., 2024).

Sequencing of the E gene and full-length virus

After amplifying the E gene from the isolated viral samples (derived from the first passage [P1] following serum inoculation onto C6/36 cells), sequencing of the full-length E gene sequence was commissioned to Sangon Biotech (Shanghai) Co., Ltd., using the Sanger sequencing method. The resulting sequences were uploaded to the National Center for Biotechnology Information (NCBI) database. For whole-genome sequencing we preselected samples with RT-qPCR Ct <25 and adequate RNA quality; from this high-quality subset we prioritized serotype representation, temporal and geographic coverage, and available clinical severity (see Supplementary Table S1), employing the Illumina method for full-length viral sequencing. Genome assembly was performed using MEGAHIT v1.2.9 (completed by ProbeGene Medical Technology [Zhejiang] Co., Ltd.) with default parameters (minimum k-mer size of 21). Detailed sequencing metrics, including genome coverage and average sequencing depth for the five whole-genome samples, are provided in Supplementary Table S5. The assembled sequences were submitted to NCBI.

Phylogenetic and molecular characteristic analysis

The isolated E gene sequences were combined with representative serotype sequences from GenBank, and Basic Local Alignment Search Tool for Nucleotides (BLASTn) (v2.16.0+) was used to verify the correctness of the protein-coding sequences (CDS) individually. Subsequently, multiple sequence alignment of the reference sequences and the E gene sequences obtained in this study was performed using Multiple Alignment using Fast Fourier (MAFFT) (v7.526). Phylogenetic trees were inferred by the Maximum Likelihood method using IQ-TREE 2 (v2.1.3). The optimal substitution model was automatically selected based on the Bayesian information criterion. Branch support was assessed via ultrafast bootstrap analysis (UFBoot2) with 1000 replicates.

Four generated DENV-1 whole-genome sequences from this study were combined with 30 representative sequences retrieved from the NCBI database. These reference sequences were predominantly sampled from Southeast Asia across various years. The earliest available sequence, EU848545.1_USA_Hawaii_1944, was included to root the tree. CDS were aligned using the ClustalW algorithm implemented in MEGA 12.

The temporal signal was evaluated via root-to-tip regression in TempEst, demonstrating a strong temporal correlation ( $R^{2} = 0.904$ Supplementary Fig. S2). Spatiotemporal evolutionary dynamics were reconstructed using BEAST X v10.5.0 (Baele et al., 2025). The evolutionary timescale was estimated using a strict molecular clock model, assuming a constant substitution rate across all branches. The initial clock rate parameter was set to 0.001 and constrained by a Continuous-Time Markov Chain (CTMC) scale reference prior. The analysis employed a coalescent exponential growth model as the tree prior and the GTR + R + I nucleotide substitution model.

Markov chain Monte Carlo chains were run for a total length of 50,000,000 generations. Two independent runs were conducted to assess reproducibility, and convergence was evaluated in Tracer, confirming highly effective sample size values for all parameters. The runs were combined using LogCombiner, and the maximum clade credibility (MCC) tree was generated with TreeAnnotator for visualization.

Structural modeling and mutation analysis of the E protein

Homology modeling was performed on the E protein variations identified in this study using the SWISS-MODEL platform and the DENV-1 amino acid sequences. For the prefusion dimer, the template Protein Data Bank Identifier (PDB ID) 7a3r was used, achieving a sequence identity of 95.95%, Global Model Quality Estimate Transform (GMQE) of 0.68, and QMEANDisCo of 0.74 ± 0.05. For the postfusion trimer, PDB ID 4gsx was selected, with a sequence identity of 96.35%, GMQE of 0.67, and QMEANDisCo of 0.72 ± 0.05. These values confirm the models’ reliability. The mutation sites (S338L and D147N) were annotated, and the structures were visualized using Python-based Molecular visualization system (PyMOL) (Delano, 2002; Nayak et al., 2009; Renner et al., 2021; Waterhouse et al., 2018).

Results

Virus isolation and amplification

During the 2024 dengue outbreak in Hainan Province, 18 acute-phase serum samples (HK-001–HK-018) confirmed by the Haikou Centers for Disease Control and Prevention (CDC) were collected. Sixteen samples were successfully isolated, with RT-qPCR serotyping identifying 15 DENV-1 and 1 DENV-3 infections. All isolates were propagated in A. albopictus C6/36 cells, where CPEs were observed from day 5 postinoculation, intensified by day 7, and became extensive by day 10 (Fig. 1). Viral supernatants were subsequently harvested and confirmed by immunofluorescence assay (Fig. 2).

FIG. 1.

Cytopathic effects of virus isolation in C6/36 cells. Serum samples were inoculated onto monolayer C6/36 cells. On day 5 postinoculation, cytopathic effects (CPE) were initially observed, manifesting as cell fusion, aggregation, and detachment (10× objective). This figure displays microscopic images from day 5 (Panel (A)), day 7 (Panel (B)), and day 10 (Panel (C)), illustrating the progressive development of CPE. Subpanels (A-1), (B-1), and (C-1) represent the negative control group; (A-2), (B-2), and (C-2) represent the sample-infected group.

FIG. 2.

Immunofluorescence assay of dengue virus infection. C6/36 cells postvirus inoculation (using a 100-fold dilution of the virus supernatant cultured for 10 days, 24 h postinfection) were fixed and subjected to immunofluorescence analysis. Panel (A): Uninfected control group. Panel (B): Results of immunofluorescence staining using anti-DENV NS3 monoclonal antibody (purchased from GeneTex). Images were acquired using a fluorescence microscope (10× objective). DENV, dengue fever.

Phylogenetic analysis and molecular clock model

To elucidate the molecular epidemiology of DENV in Hainan Province, we utilized the E gene sequences of 16 DENV strains determined in this study (GenBank accession numbers: PV478122–PV478135, PX317703, and PX317704; Supplementary Table S4). A comprehensive phylogenetic tree including all identified serotypes is presented in Supplementary Figure S1. Given the predominance of DENV-1 in our samples, we further constructed a serotype-specific Maximum Likelihood tree for DENV-1 to resolve the topological relationships with higher resolution (Fig. 3). The genetic evolutionary analysis of the E gene elucidated information on DENV isolates from the Haikou region in Hainan Province. Based on the E gene sequences, 15 strains (HK-001, HK-002 to HK-010, and HK-013 to HK-018) were classified as DENV-1 genotype I, with 14 local Hainan isolates located in the DENV-1 genotype I (genotype I) evolutionary branch; this suggests that the primary circulating strain in this outbreak was DENV-1. This branch formed a tight cluster with the 2023 circulating strain from Guangdong Province (PP563845.1), indicating that this outbreak may have been transmitted from Guangdong to the local area, subsequently establishing cryptic circulation locally. Although HK-010 also belongs to the DENV-1 genotype I, it did not cluster with the other DENV-1 strains in this study and shares 99.23% homology with sequences from Southeast Asian countries such as Thailand and India. HK-002, located in the same branch, is a DENV-3 strain, closely related to recent strains from India and the Philippines.

FIG. 3.

Phylogenetic analysis of DENV-1 strains based on E gene sequences. The Maximum Likelihood (ML) tree was inferred using IQ-TREE 2 (v2.1.3) based on E gene sequences of 16 DENV isolates from the 2024 Hainan outbreak (GenBank accession nos. PV478122–PV478135, PX317703, and PX317704) and representative reference strains. Multiple sequence alignment was performed using MAFFT (v7.526). The optimal substitution model was selected based on the Bayesian information criterion (BIC). Ultrafast bootstrap support values (UFBoot2) calculated from 1000 replicates are shown at nodes (values ≥70% are displayed). The scale bar represents nucleotide substitutions per site. The strains from the Hainan outbreak are highlighted in red, while all other reference strains are shown in blue.

To further resolve the evolutionary dynamics indicated by the E gene phylogeny, an MCC tree was reconstructed based on the whole genomes of four representative isolates: PX352011 (HK-DENV017), PV478125 (HK-DENV006), PX317703 (HK-DENV001), and PX352012 (HK-DENV010). The whole-genome topology (Fig. 4) was consistent with the E gene analysis. HK-DENV017, HK-DENV006, and HK-DENV001 formed a monophyletic clade that clustered tightly with recent strains from Guangdong, including the late-2023 Foshan isolate (PP563845.1). This genomic proximity supports the epidemiological link to Guangdong. In contrast, HK-DENV010 exhibited distinct genetic divergence and segregated into a separate lineage. It clustered alongside a 2018 Indian isolate (MZ312924.1), corroborating the independent origin previously suggested by the E gene data.

FIG. 4.

Time-scaled phylogenomic analysis of DENV-1. A maximum clade credibility (MCC) tree was reconstructed using BEAST X v10.5.0 based on the coding sequences (CDS) of four representative isolates from the current outbreak (HK-DENV001, HK-DENV006, HK-DENV010, and HK-DENV017) and 30 global reference genomes. The reference dataset primarily includes representative strains from Southeast Asia and China, along with historical sequences (such as a 1944 Hawaii isolate) to root the tree. Sequence alignment was performed using the ClustalW algorithm in MEGA 12. Following temporal signal assessment via TempEst (R² = 0.904), the evolutionary timescale was estimated under a strict molecular clock model. The horizontal axis represents evolutionary time in years, spanning from approximately 1899 to 2024. Branch colors correspond to the estimated mean node height, reflecting the temporal progression from ancestral lineages to recent divergence events.

E gene mutation analysis

To compare the genetic characteristics of the isolated DENV-1 strains, amino acid-based multiple sequence alignment was performed using reference DENV-1 sequences from the same phylogenetic branch. High sequence similarity was observed in the E protein region; however, the 14 isolates harbored 2 unique amino acid substitutions, D147N and S338L. D147 is located in domain I, which is associated with conformational changes and membrane fusion, whereas S338 resides in domain III, an immunoglobulin-like domain containing key receptor-binding sites and neutralizing antibody epitopes involved in tissue tropism and immune escape (Hu et al., 2019) (Figs. 5 and 6).

FIG. 5.

Schematic diagram of DENV-1 E protein conformation. Domain I (light blue) consists of three discontinuous segments; Domain II (pink) includes two discontinuous segments; Domain III (dark gray) is located at the N-terminus of the black stem region. The C-terminal light gray segment represents the transmembrane anchor region. Two unique mutation sites are marked in red.

FIG. 6.

Structural prediction of DENV-1 E protein prefusion and postfusion conformations and localization of key mutation sites. (A) Prefusion conformation (dimer) modeled using the template PDB ID 7a3r (Seq ID: 95.95%, GMQE = 0.68, QMEANDisCo = 0.74 ± 0.05). (B) Postfusion conformation (trimer) modeled using the template PDB ID 4gsx (Seq ID: 96.35%, GMQE = 0.67, QMEANDisCo = 0.72 ± 0.05). In both panels, the left side presents the protein backbone in cartoon representation, while the right side shows the molecular surface representation. Key amino acid substitutions identified in this study (S338L and D147N) are highlighted in red and labeled. The fusion loop region is indicated to define the structural orientation. Models were constructed using SWISS-MODEL and visualized via PyMOL. Scale bar = 10 Å.

Discussion

China has experienced multiple dengue outbreaks since the disease reemerged in Foshan, Guangdong, in 1978 (Xiong and Chen, 2014). Following its inclusion in the list of notifiable diseases in the 1990s (Lai et al., 2015), outbreaks were mainly concentrated in southern provinces, though the geographical range and frequency of outbreaks have gradually increased in recent years (Ni et al., 2024). Hainan Province has faced three major dengue outbreaks: 1978 (type 3), 1985–1988 (type 2), and 1991 (type 1) (Fan et al., 1989; Liu et al., 2021; Qiu et al., 1991). In July 2024, the WHO classified DENV as a high-risk pathogen, urging global prevention and control efforts (World Health Organization, 2024). Despite repeated cases in Hainan, there is still limited understanding of the genetic characteristics and epidemic patterns of local strains.

Phylogenetic analysis indicated that the 2024 dengue outbreak in Hainan was dominated by DENV-1 genotype I, forming a tight cluster with the Guangdong strain (PP563845.1). While this genetic similarity suggests a potential introduction from Guangdong, the exact transmission route remains unclear without detailed mobility and vector surveillance data. Among the 18 specimens, one was a DENV-3 infection (HK-002), and another (HK-010) did not match the local circulating strains, suggesting independent transmission chains possibly introduced from Southeast Asia. These strains may be in the early stages of transmission and require close monitoring.

The E protein is the primary antigenic component of DENV and the main target of host immune pressure (Katzelnick et al., 2015; Sessions et al., 2015), comprising three domains (I–III; Table 1). Comparative analysis with reference DENV-1 E proteins identified two mutations: D147N in domain I and S338L in domain III. Structural analysis of the prefusion E protein (HK-008) showed that residues D147 and S338 are surface-exposed, suggesting potential involvement in host receptor interactions and altered binding affinity.

Table 1.

Structural Domain Architecture of Dengue Virus-1 Envelope Glycoprotein

Structural domain	Amino acid range (Meng et al., 2023)	Functional characteristics
Domain I (DI)	1–51, 132–187, 278–293 (discontinuous segments)	Located in the central region of the E protein, composed of β-sheets, stabilized by disulfide bonds, involved in conformational changes and membrane fusion functions (S et al., 2021).
Domain II (DII)	52–131, 188–277 (discontinuous segments)	Folded into a finger-like structure, contains the fusion loop, responsible for the fusion process between the virus and host cell membrane (Park et al., 2020; S et al., 2021).
Domain III (DIII)	294–393 (continuous region)	An independent and complete immunoglobulin-like structure, contains receptor-binding sites and neutralizing antibody epitopes, determines viral tissue tropism and immune escape (Ali et al., 2025; Pierson et al., 2008).

Domain I consists of discontinuous segments in the E protein and participates in viral conformational changes and membrane fusion processes. Studies have shown that T226K/G228E enhances viral infectivity in mosquito or mammalian hosts (Chen et al., 2022), while mutations at position 160 in the E protein increase antibody neutralization activity, leading to replacement of circulating genotypes (Wang et al., 2016). These findings indicate that changes in domain I can affect viral infectivity; therefore, we speculate that the change at D147 might influence the virus’s infectivity (Table 2). Domain III contains type- and subtype-specific structural antigenic epitopes that can induce neutralizing antibodies (Park et al., 2020). Additionally, domain III features an independent and complete immunoglobulin-like structure, containing receptor-binding regions and neutralizing antibody epitopes, playing important roles in the viral life cycle.

Table 2.

Predicted Functional Impact of Key Mutations in Dengue Virus-1 E Protein Domains

Mutation site	Structural domain	Predicted primary impact	Supporting studies
D147N	Domain I (DI)	May enhance viral stability and membrane fusion efficiency	DI-T226K/G228E studies (Chen et al., 2022)
S338L	Domain III (DIII)	May enhance receptor-binding capability, promoting immune escape	DIII-324/351/380 site studies (Jiang et al., 2024)

Thus, the S338 mutation may contribute to immune escape by enhancing receptor-binding capability; however, this inference is based on structural modeling and requires experimental validation. This hypothesis is supported by recent studies on the 2023 autochthonous outbreak in Italy (Carletti et al., 2024) and broader analyses of host–virus interactions (Ghosh et al., 2025), which highlight the critical role of E protein surface mutations in viral adaptation and antibody evasion. The co-occurrence of D147N and S338L is reported here for the first time, and their effects on viral transmission and pathogenicity, including potential synergistic interactions, warrant further investigation using reverse genetics and in vivo models.

From a public health perspective, these findings underscore the need to strengthen local control strategies in Hainan by integrating sustained mosquito surveillance with real-time genomic monitoring for early detection of viral introductions. In addition to conventional vector control, targeted vaccination of high-risk populations may represent an important complementary approach to reducing outbreak risk in this tropical region.

This study has several limitations. First, the limited sample size and the lack of detailed transmission dynamics restrict attribution analyses and may not fully capture the genetic diversity circulating in the population, although the samples were sufficient for initial pathogen identification and phylogenetic placement. Second, although amino acid substitutions (D147N and S338L) were identified in functionally important regions, their biological effects were inferred solely from structural modeling, without experimental validation, such as neutralization or pseudovirus infectivity assays. Therefore, any inferred effects on viral transmissibility or immune escape should be considered hypothesis-generating and require further in vitro and in vivo investigation.

The 2024 epidemic underscores the growing challenge of dengue in China. As the southernmost province, Hainan’s persistently high temperatures, abundant rainfall, and frequent trade exchanges increase its susceptibility to viral introduction and establishment. Surveillance data indicate that A. albopictus has largely replaced A. aegypti as the predominant vector in the region, a situation further complicated by documented insecticide resistance in local mosquito populations (Chen et al., 2020; Zhao et al., 2023). In addition, the tropical climate accelerates transmission by shortening the viral extrinsic incubation period and sustaining high vector densities, resulting in a substantially prolonged transmission season compared with temperate regions of China (Yue et al., 2021).

Against this high-risk ecological background, our study characterizes the molecular and epidemiological features of the 2024 Haikou outbreak. By clarifying viral origins and transmission chains, these findings provide a scientific basis for optimizing local control strategies, including targeting insecticide-resistant vectors, and offer a framework for future research into DENV pathogenicity and transmissibility.

Ethics Review

The experimental protocol was approved by the Ethics Review Committee of the Naval Medical University, having passed all legal and ethical reviews. All experiments were conducted in accordance with relevant guidelines and regulations.

Authors’ Contributions

Z.-M.C.: Investigation (lead), data curation (lead), and writing—original draft (lead). B.-H.X.: Methodology (lead) and data curation (lead). P.L.: Resources (lead) and data validation (lead). X.-H.N.: Conceptualization (lead). Z.-H.Z.: Formal analysis (lead). Q.Z.: Investigation (supporting). H.R. and Y.-Z.Z.: Supervision (lead) and writing—review and editing (equal for both). All authors contributed critically to the drafts and gave final approval for publication.

Footnotes

Author Disclosure Statement

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article.

Funding Information

This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.

Supplemental Material

References

Ali

, Saleem

, Abdaal

, et al. Dengue virus proteins: Dual drivers of pathogenesis and vaccine development. J Clin Virol Plus, 2025; 5(3):100227; doi: 10.1016/j.jcvp.2025.100227

Anasir

, Ramanathan

, Poh

. Structure-based design of antivirals against envelope glycoprotein of dengue virus. Viruses, 2020; 12(4):367; doi: 10.3390/v12040367

Akshatha

, Pujar

, Sethu

, et al. Dengue structural proteins as antiviral drug targets: Current status in the drug discovery & development. Eur J Med Chem, 2021; 221:113527; doi: 10.1016/j.ejmech.2021.113527

Baele

, Ji

, Hassler

, et al. BEAST X for Bayesian phylogenetic, phylogeographic and phylodynamic inference. Nat Methods, 2025; 22(8):1653–1656; doi: 10.1038/s41592-025-02751-x

Carletti

, Carli

, Spezia

, et al. Genetic and structural characterization of dengue virus involved in the 2023 autochthonous outbreaks in central Italy. Emerg Microbes Infect, 2024; 13(1):2420734; doi: 10.1080/22221751.2024.2420734

Cattarino

, Rodriguez-Barraquer

, Imai

, et al. Mapping global variation in dengue transmission intensity. Sci Transl Med, 2020; 12(528):eaax4144; doi: 10.1126/scitranslmed.aax4144

Chen

, Li

, Wang

, et al. Widespread multiple insecticide resistance in the major dengue vector Aedes albopictus in Hainan Province. China. Pest Manag Sci, 2020; 76(8):2723–2731; doi: 10.1002/ps.5851

Chen

, Zhang

, Guo

, et al. Neighboring mutation-mediated enhancement of dengue virus infectivity and spread. EMBO Rep, 2022; 23(11):e55671; doi: 10.15252/embr.202255671

Delano

. PyMOL: An open-source molecular graphics tool. ScienceOpen; 2002. Available from: https://www.scienceopen.com/book?vid=4362f9a2-0b29-433f-aa65-51db01f4962f

10.

Diamond

, Pierson

. Molecular insight into dengue virus pathogenesis and its implications for disease control. Cell, 2015; 162(3):488–492; doi: 10.1016/j.cell.2015.07.005

11.

, Zhang

, Hu

, et al. Phylogenetic analysis of the dengue virus strains causing the 2019 dengue fever outbreak in Hainan, China. Virol Sin, 2021; 36(4):636–643; doi: 10.1007/s12250-020-00335-x

12.

Fan

, Yu

, Cosgriff

. The reemergence of dengue in China. Rev Infect Dis, 1989; 11(Suppl 4):S847–S853; doi: 10.1093/clinids/11.supplement_4.s847

13.

Fraenkel

, Nabeshima

, Xayavong

, et al. The development of new primer sets for the amplification and sequencing of the envelope gene of all dengue virus serotypes. Microorganisms, 2024; 12(6):1092; doi: 10.3390/microorganisms12061092

14.

Ghosh

, Mondal

, Sadhukhan

, et al. Dengue virus and the host immune system: A battle of immune modulation, response and evasion. Pathogens, 2025; 14(11):1132; doi: 10.3390/pathogens14111132

15.

, Zhu

, Li

, et al. A broadly neutralizing germline-like human monoclonal antibody against dengue virus envelope domain III. PLoS Pathog, 2019; 15(6):e1007836; doi: 10.1371/journal.ppat.1007836

16.

Iglesias

, Gamarnik

. Dynamic RNA structures in the dengue virus genome. RNA Biol, 2011; 8(2):249–257; doi: 10.4161/rna.8.2.14992

17.

Jiang

, Huang

, Ruan

, et al. Envelope domain III E324, E351, and E380 mutations lever adaptive evolution of DENV-1 genotype I. J Virol, 2024; 98(10):e01183–24; doi: 10.1128/jvi.01183-24

18.

Katzelnick

, Fonville

, Gromowski

, et al. Dengue viruses cluster antigenically but not as discrete serotypes. Science, 2015; 349(6254):1338–1343; doi: 10.1126/science.aac5017

19.

Katzelnick

, Gresh

, Halloran

, et al. Antibody-dependent enhancement of severe dengue disease in humans. Science, 2017; 358(6365):929–932; doi: 10.1126/science.aan6836

20.

Kuhn

, Zhang

, Rossmann

, et al. Structure of dengue virus: Implications for flavivirus organization, maturation, and fusion. Cell, 2002; 108(5):717–725; doi: 10.1016/S0092-8674(02)00660-8

21.

Lai

, Huang

, Zhou

, et al. The changing epidemiology of dengue in China, 1990–2014: A descriptive analysis of 25 years of nationwide surveillance data. BMC Med, 2015; 13:100; doi: 10.1186/s12916-015-0336-1

22.

Liu

, Wu

, Liu

, et al. The origin and molecular epidemiology of dengue fever in Hainan Province, China, 2019. Front Microbiol, 2021; 12:657966; doi: 10.3389/fmicb.2021.657966

23.

Meng

, Hu

, Zhang

, et al. Isolation and genetic evolution of dengue virus from the 2019 outbreak in Xishuangbanna, Yunnan Province, China. Vector Borne Zoonotic Dis, 2023; 23(6):331–340; doi: 10.1089/vbz.2022.0091

24.

Nayak

, Dessau

, Kucera

, et al. Crystal structure of dengue virus type 1 envelope protein in the postfusion conformation and its implications for membrane fusion. J Virol, 2009; 83(9):4338–4344; doi: 10.1128/JVI.02574-08

25.

, Cai

, Ren

, et al. Epidemiological characteristics and transmission dynamics of dengue fever in China. Nat Commun, 2024; 15(1):8060; doi: 10.1038/s41467-024-52460-w

26.

Park

, Lee

, Khai

, et al. Addition of partial envelope domain II into envelope domain III of dengue virus antigen potentiates the induction of virus-neutralizing antibodies and induces protective immunity. Vaccines (Basel), 2020; 8(1):88; doi: 10.3390/vaccines8010088

27.

Pierson

, Fremont

, Kuhn

, et al. Structural insights into the mechanisms of antibody-mediated neutralization of flavivirus infection: Implications for vaccine development. Cell Host Microbe, 2008; 4(3):229–238; doi: 10.1016/j.chom.2008.08.004

28.

Qiu

, Chen

, Ho

, et al. The first epidemic of dengue hemorrhagic fever in the People’s Republic of China. Am J Trop Med Hyg, 1991; 44(4):364–370; doi: 10.4269/ajtmh.1991.44.364

29.

Renner

, Dejnirattisai

, Carrique

, et al. Flavivirus maturation leads to the formation of an occupied lipid pocket in the surface glycoproteins. Nat Commun, 2021; 12(1):1238; doi: 10.1038/s41467-021-21505-9

30.

Sessions

, Wilm

, Kamaraj

, et al. Analysis of dengue virus genetic diversity during human and mosquito infection reveals genetic constraints. PLoS Negl Trop Dis, 2015; 9(9):e0004044; doi: 10.1371/journal.pntd.0004044

31.

Venkatesan

. Global upsurge in dengue in 2024. Lancet Infect Dis, 2024; 24(10):e620; doi: 10.1016/S1473-3099(24)00609-1

32.

Wang

, Katzelnick

, Montoya

, et al. Evolutionarily successful Asian 1 dengue virus 2 lineages contain one substitution in envelope that increases sensitivity to polyclonal antibody neutralization. J Infect Dis, 2016; 213(6):975–984; doi: 10.1093/infdis/jiv536

33.

Waterhouse

, Bertoni

, Bienert

, et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res, 2018; 46(W1):W296–W303; doi: 10.1093/nar/gky427

34.

World Health Organization. Dengue: Guidelines for diagnosis, treatment, prevention and control, New ed. World Health Organization: Geneva; 2009. http://www.ncbi.nlm.nih.gov/books/NBK143157/

35.

World Health Organization. Dengue – Global situation. In: Disease Outbreak News. WHO; 2024. WHO stated that the overall risk at the global level is assessed as high.

36.

Xiong

, Chen

. Epidemiology of dengue fever in China since 1978. Nan Fang Yi Ke Da Xue Xue Bao, 2014; 34(12):1822–1825.

37.

Yue

, Liu

, Ren

, et al. Assessing the suitability for Aedes albopictus and dengue transmission risk in China with a delay differential equation model. PLoS Negl Trop Dis, 2021; 15(6):e0009444; doi: 10.1371/journal.pntd.0009444

38.

Zhang

, Zhang

, Ogata

, et al. Conformational changes of the flavivirus E glycoprotein. Structure, 2004; 12(9):1607–1618; doi: 10.1016/j.str.2004.06.019

39.

Zhao

, Ran

, Xing

, et al. Evolution of knockdown resistance (kdr) mutations of Aedes aegypti and Aedes albopictus in Hainan Island and Leizhou Peninsula, China. Front Cell Infect Microbiol, 2023; 13:1265873; doi: 10.3389/fcimb.2023.1265873

40.

Zompi

, Harris

. Original antigenic sin in dengue revisited. Proc Natl Acad Sci U S A, 2013; 110(22):8761–8762; doi: 10.1073/pnas.1306333110

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.42 MB

0.11 MB

0.36 MB