Introduction

Dengue fever (DF) is an acute, febrile, and self-limiting systemic viral infection transmitted by the Aedes aegypti and Aedes albopictus mosquitos worldwide. The disease is endemic in tropical and subtropical regions such as Southeast Asia, the Eastern Mediterranean Region, Sub-Saharan Africa, and Latin America, and affects about 40% of the global population1. It has been a major health concern for decades, causing about 390 million illnesses annually and killing up to 40% of severe patients2,3. The first instance of dengue fever in Pakistan was recorded in November 1994, when a hospitalized woman presented with multiple dengue symptoms. Since then, Pakistan has seen multiple cyclical dengue epidemics in different regions. This can be attributed to the dense population, poor socioeconomic conditions, poor healthcare systems, inadequate water and sanitation facilities, and frequent natural calamities such as floods, earthquakes, and climate changes4,5.

Dengue fever can present as mild fever, dengue hemorrhagic fever (DHF), or dengue shock (DSS) syndrome. DHF and DSS have been recorded in several patients in Pakistan during dengue epidemics in recent decades. The initial symptoms of DF include fever, headache, muscle soreness, discomfort, a transient rash, petechiae, ecchymosis, and venipuncture bleeding. In the critical phase, plasma leakage causes shock, fluid accumulation, and severe bleeding. Some patients get acute liver failure, renal dysfunction, and encephalopathy with convulsions, as well as lower gastrointestinal bleeding6. There is no specific medication for dengue fever, however, mild fever can be treated with antipyretics, but severe infections necessitate hospitalization and enough fluids7. Since the initial outbreak, the disease’s clinical symptoms have evolved over temporal changes and are constantly developing. These various clinical signs are assumed to be caused by the prevalence of a certain serotype (DENV 1–4) at a given time8,9. Pakistan lacks a surveillance program for vector-borne diseases like dengue fever and therefore relies on a passive approach limited to select areas. Hence, the characterization of the dengue virus serotypes and surveillance of associated clinical features of patients in Karachi would serve as a scientific basis for dengue and other flaviviral infection control and prevention, as well as dengue immunizations in Pakistan. In this study, we performed serotyping of dengue virus and clinical surveillance of dengue patients in Karachi, Pakistan, from 2019 to 2023.

Materials and methods

Study design, ethical approval, sample collection and preparation

This cross-sectional study was conducted between January 2019 and March 2023 at the National Institute of Virology, Dr. Panjwani Center for Molecular Medicine and Drug Research, University of Karachi in collaboration with Indus Hospital, Abbasi Shaheed Hospital, and Essa Diagnostic Laboratory in Karachi. The study protocol was approved by the institutional independent ethics committees and collaborating hospitals (ICCBS/IEC-042/HB-2018/ProtocoV1.0) and all methods were performed under the relevant guidelines and regulations.

The social science statistics online tool estimated a sample size of 385, with a 95% confidence level, a margin of error of ± 4, and a 0.5 estimated proportion. We adopted the purposive sampling approach including patients who experienced symptoms during infection and had their immune-chromatographic (ICT) results confirmed positive. Overall, 1500 samples were included in this study. Clinical observations were also recorded, including white blood cell (WBC) and platelet counts, and clinical signs and symptoms of selected 486 patients. All samples were collected using commercially available ethylenediamine tetraacetic acid (EDTA) anticoagulant tubes. Informed consent was obtained from all subjects and their legal guardian(s) about the DENV epidemiological and clinical surveillance before they were included in the study. The plasma was separated by centrifugation at 3000 rpm for five minutes, then collected into a fresh tube and stored at −20°C until needed. The demographic and clinical information of the patient was obtained from guardians and on-duty medical officers, respectively.

Extraction of RNA and synthesis complementary DNA (cDNA)

This study considered genotyping samples that tested positive for the NS1 (nonstructural protein 1) antigen of the dengue virus. QIAamp viral RNA isolation kit (Qiagen, Heldin, Germany) was used to extract RNA from 140 µl plasma sample in 50 µl elution buffer as per the manufacturer’s protocol. The quality and quantity of the extracted RNA were determined using NanoDrop Multiscan Go Spectrophotometer µDrop™ Plate. RevertAid First strand cDNA synthesis kit (Thermo Scientific Inc.) was used to synthesize cDNA from the extracted RNA. Briefly, RNA 1 µg, was incubated with DNase I in a 10 µL reaction mixture for 30 min at 37 ℃ followed by 10 min of incubation with 50 mM EDTA at 65 ℃. This RNA was used to synthesize cDNA in 20 µL reaction mixture at a thermal profile: incubation at 25 °C for five minutes followed by 42°C for 60 min, and termination at 70 °C for 5 min in the T100 Thermal Cycler (Bio-Rad). The resulting cDNA was used to classify all four dengue virus serotypes using multiplex PCR.

Multiplex PCR

The first PCR reaction was carried out with dengue virus universal primers D1 and D210. The positive samples were further analyzed with specialized serotyping primers (TS1, TS2, TS3, and DEN4) based on their highest homology to four serotypes11. All primers were purchased from the macrogen’s commercially available services through a local supplier (Table 1). Briefly, a 25 µl reaction was prepared using absolute mastermix Molequle-On (Cat# AM-M-0011250, New Zealand) and PCR was performed for 35 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 60 s, extension at 72 °C for 120 s12. Amplified PCR products (511 bp, 482 bp, and 119 bp gene fragments) were visualized in 2% TAE (1× Tris-acetate-EDTA) agarose gel. The respective regions of different serotypes in the gel, characterized by Thermo Scientific GeneRuler 100 bp DNA Ladder, were excised and purified using an MQ PCR Product Purification Kit (Cat# PPK-M-002100, New Zealand).

Sanger sequencing

A total of 14 samples were subjected to Sanger sequencing using D1 and D2 primers. BigDye® Direct Sanger Sequencing Kit was used for sequencing PCR of these samples as per the manufacturer’s protocol. Purified PCR sequencing products were subjected to Sanger sequencing. Data were analyzed by the nucleotide basic local alignment search tool (BLAST).

Table 1 Oligonucleotide primers and their sequences used in conventional PCR and multiplex PCR for serotyping of dengue virus.

Results

Demographic information of patients with DENV infection

This study analyzed 1500 dengue patients from Indus Hospital, Abbasi Shaheed Hospital, and Essa Diagnostic Laboratory in Karachi, from 2019 to 2023. The study included 53.87% (808/1500) males and 46.13% (692/1500) females. Patients were divided into three groups based on age differences, i.e., 1–17 years, 18–60 years, and > 60 years with contributions of 18.93% (284/1500), 69.14% (1037/1500), and 11.93% (179/1500) respectively. The distribution of samples among the three groups of the two genders was not significantly different (P = 0.357, χ2 = 2.059, 2).

Trends in seropositivity of patients with DENV infection

An analysis of seropositivity and co-presence of NS1 with IgM and IgG was made using results from ICT. The results showed that 57.8% (531/919) of patients were positive for NS1, 39.8% (488/1226) for IgM, 89.77% (1018/1135) for IgG, 39.2% (138/352) for NS1 + IgM, and 89.9% (96/108) for NS1 + IgG, and 85.7% (359/419) for IgM + IgG. The chi-square test revealed a significant difference (P = 0.001, χ2 = 54.83, 2) in positivity across immune-chromatographic detection (Fig. 2A; Table 2).

Table 2 Immune-chromatographic detection of dengue cases in Karachi, Pakistan, from 2019 to 2023.

Phylogenetic analysis and serotyping

Since only the serotype region (D1/D2 amplified region, Fig. 1A, B) was targeted for sequencing, mutational analysis was not in the scope rather confirmation of serotypes and serotyping primers were aimed. The set criteria were 10 samples of each serotype but no sample of serotypes 3 and 4 could be detected using serotyping primers. Moreover, many samples were rejected due to poor quality of DNA determined by Qubit Fluorometric quantification. A total of 14 samples were successfully sequenced and phylogenetic analysis confirmed four samples of DENV 1 and 10 samples of DENV 2. This result was in 100% agreement with PCR results representing 20% and 80% of the dengue virus circulating serotypes, 1 and 2 respectively (Fig. 2B).

Fig. 1
figure 1

shows the ethidium bromide-stained gel electrophoresis image of the amplified products of D1 and TS1, 482 bp (A) and D1 and TS2, 119 bp (B) of DEN 1 and DEN 2 serotypes.

Laboratory findings for patients with the DENV infection

The study investigated WBCs and platelet counts, and the clinical manifestations of 486 patients enrolled at Indus Hospital. Males and females contributed to 54.53% (265/486) and 45.47% (221/486) of the samples respectively. Similarly, the age groups 1–17, 18–60, and > 60 contributed 19.54% (95/486), 68.51% (333/486), and 11.93% (58/486) respectively. Chi-square test revealed no significant difference in WBCs (P = 0.734, χ2 = 0.618, 2, P = 0.267, χ2 = 2.641, 2, and P = 0.356, χ2 = 2.064, 2) as well as platelets count (P value 0.272, χ2 = 2.603, 2, P value 0.638, χ2 = 0.898, 2, and P value 0.158, χ2 = 3.686, 2) among the three age groups of two genders (Table 3).

Fig. 2
figure 2

shows trends in the seropositivity of dengue patients (A) and clinical serotypes of dengue virus (B) in Karachi from 2019 to 2023.

Clinical surveillance of patients with DENV infection

Data on clinical presentations among the three groups showed that shortness of breath (SOB) (P value < 0.001, χ2 = 44.30, 2), fever (P value < 0.001, χ2 = 14.09, 2), body ache (P value < 0.002, χ2 = 12. 19, 2), headache (P value < 0.049, χ2 = 6.003, 2), nausea (P value < 0.003, χ2 = 11.13, 2), vomiting (P value < 0.001, χ2 = 13.03, 2), diarrhea (P value < 0.025, χ2 = 7.36, 2), and epistaxis (P value < 0.025, χ2 = 7.32, 2) revealed significant difference but weakness, chest pain, headache, Skin rashes, retro-orbital pain, abdominal pain, gum bleeding, epistaxis, ear canal bleeding, hematemesis, and melena were not significant (Table 4).

We also analyzed the effect of gender on significant clinical features. We found that SOB (P value < 0.918, χ2 = 0.010, 1), fever (P value < 0.565, χ2 = 330, 1), body ache (P value < 0.359, χ2 = 0.841, 1), head ache (P value < 0.057, χ2 = 3.611, 2), nausea (P value < 0.150, χ2 = 2.069, 1), diarrhea (P value < 0.202, χ2 = 1.621, 1), and epistaxis (P value < 0.410, χ2 = 0.677, 1) were not significant, however, vomiting (P value < 0.036, χ2 = 4.366, 1) was significantly more common in females than males (Table 5).

We also examined data on the duration of the significant clinical features. A significantly high number of patients had a fever (P value < 0.011, χ2 = 16.55, 3) for ≥ 1 – ≤ 7 days. However, the duration of SOB, body aches, nausea, vomiting, and diarrhea was not significant across the three age groups (Table 6).

Table 3 WBCs and platelet counts of DENV patients in Karachi, Pakistan, from 2019 to 2023.
Table 4 Clinical features observed in DENV patients in Karachi, Pakistan, from 2019 to 2023.
Table 5 Gender effects on clinical features in dengue patients in Karachi, Pakistan, from 2019 to 2023.
Table 6 Persistence of clinical features in DENV patients in Karachi, Pakistan, from 2019 to 2023.

Discussion

Pakistan, an Eastern Mediterranean country, has seen repeated dengue fever outbreaks in recent decades13,14,15. It lacks a surveillance program for vector-borne diseases like malaria and dengue, instead relying on a passive strategy restricted to specific areas. Karachi, one of the most populous cities of Pakistan, has a significant burden of Aedes aegypti and Aedes albopictus due to suitable breeding sites, weather, and rapid and unplanned urbanization of squatter areas1,16,17. This study examined the infectivity rate and clinical characteristics of dengue patients in Karachi from 2019 to 2023. We found that the seropositivity of the dengue patients was 57.8%, 39.8%, and 89.77% for NS1, IgM, and IgG respectively (Fig. 2A; Table 2), with 20 and 80% circulation of DENV 1 and DENV 2 serotypes (Fig. 2B). It showed that males are at a greater risk of dengue infections (53.87%) than females (46.13%), though not significantly. DENV infections varied with age, being highest in the 18–60 years (69.14%) and the lowest in the > 60 years (11.93%) age group. Clinical features such as SOB, fever, body pains, headache, nausea, vomiting, diarrhea, and epistaxis differed significantly amongst the three groups (Table 4), whereas vomiting differed significantly between the two genders only (Table 5). The duration of fever varied significantly between ≥ 1–≤ 7 days among the three age groups (Table 6).

Studies conducted around the world revealed a variety of demographic trends in dengue infections. In North America either an equal or higher prevalence of dengue was seen in females than males. In contrast, in the Asian population, men are more affected than females. We also found a higher infection rate in males than females, but the difference was not significant. The difference between these observations can be related to prevailing social gender norms in different regions. For example, in Pakistan males have traditionally been responsible for outdoor activities, whereas women mostly work as housewives and are dressed to conceal their entire bodies, influencing their risk of contracting vector-borne disease18,19,20,21. Age group 18–60 has a high incidence of dengue infection followed by 1–17 similar to previously reported studies19,22,23,24.

Different studies use a variety of approaches to detect dengue virus infection. NS1, a glycoprotein of DENV, is implicated in host defense evasion, pathogenesis, and early stages of viral replication25. Sera from infected patients can carry it during the acute phase of the disease, and some studies have found its presence for up to 9–10 days after dengue symptoms manifest26,27,28. On the other hand, the host defense can be harnessed by detecting IgM and IgG antibodies for the presence of a primary or secondary infection. The sensitivity of each diagnostic method depends on whether the infection is primary or secondary, as well as the dengue serotype29.

We found that the seroprevalence of NS1 was 57.8%, IgM 39.8%, and IgG 89.77% (Fig. 2A; Table 2). In one study in the Khyber Pakhtunkhwa province, IgM was 32% and IgG was 20%30, whilst another reported that NS1 was 25.26, IgM 41.05%, and IgG 8.42%31. A study conducted in Multan showed IgM and IgG at 54.17% and 27.6%32. A study conducted in Karachi showed NS1 at 73.5% while a few at 2.6%, and 5.3% for IgM and IgG respectively, however, the sample size was small4. The disparities between our study and other reported studies could be attributed to changes in sample size, study duration, and temporal variation. The higher IgG values observed in our study could be attributed to unnoticed mild dengue infections that occurred over time in Karachi.

Typically, dengue manifests itself as fever, pain, or rash however bleeding and gastrointestinal complications might occur at varying degrees33. Studies showed fever in all patients of dengue with the most common associated features being nausea and/or vomiting, severe headache, anorexia, and retro-orbital pain34,35,36 as we found in the study Similarly, consistent with the previous studies37,38, the current study also shows a decreasing proportion of mucosal bleeding (epistaxis, gum bleeding, and otorrhagia), and GI bleeding (hematemesis and/or melena). Furthermore, unlike other studies that identified skin rash as a critical feature, we did not observe it significantly (Table 4). The decreased prevalence of skin rash could be attributed to prevalent secondary infection rather than primary infection24,39, as evidenced by 89.7% IgG and DENV 1 and DENV 2 serotype circulation40.

Vomiting and diarrhea were the most common digestive symptoms in this study, however abdominal pain and weakness were not frequent (Table 4). Studies also showed a significant increase in diarrhea and abdominal pain over the past decade41,42,43. Khan et al. and Guzman and Kouri 2003 observed that abdominal pain was more likely to occur in DEN-3 serotypes than in DEN-1 and DEN-240,44. We found DENV 1 and 2 serotypes were more prevalent in 2019–2023 in Karachi, which could explain the decrease in abdominal pain45, nevertheless, further evidence from substantial serotype-specific studies is needed.

Limitations

Diagnosis of dengue infections based on NS1 antigen, IgM, and IgG antibody assay may produce false positives from cross-reactivity with other flaviviral infections46. We performed RT-PCR confirmation on a small number of samples intended to be sequenced. The number sequenced was also small, due to cost constraints, which are some of the drawbacks of our study. The strength of this study is the inclusion of serologically confirmed positive cases, a large number of female samples, patients of all ages, a long period, and a detailed clinical profile of the enrolled cases, which can assist local health authorities in allocating resources for capacity development.