
==== Front
Bioinform Biol Insights
Bioinform Biol Insights
BBI
spbbi
Bioinformatics and Biology Insights
1177-9322
SAGE Publications Sage UK: London, England

10.1177/11779322241272399
10.1177_11779322241272399
Research Article
Emergence of SARS-CoV-2 Variants Are Induced by Coinfections With Dengue
https://orcid.org/0000-0003-1185-6720
Al-Emran Hassan M 1
Rahman Fazlur 2
Sarkar Laxmi 3
Das Prosanto Kumar 3
Mondol Provakar 3
Yesmin Suriya 3
Sultana Pipasha 3
Ahammed Toukir 3
Parvez Rasel 3
Hasan Md Shazid 3
Sarkar Shovon Lal 3
Rahman M Shaminur 3
Hossain Anamica 4
Rahman Mahmudur 5
Islam Ovinu Kibria 3
Islam Md Tanvir 3
Nigar Shireen 6
Akter Selina 3
Rubayet Ul Alam A S M 3
Rahman Mohammad Mahfuzur 7
Jahid Iqbal Kabir 38
Hossain M Anwar 8
1 Department of Biomedical Engineering, Jashore University of Science and Technology, Jashore, Bangladesh
2 Department of Accounting and Information Systems, Jashore University of Science and Technology, Jashore, Bangladesh
3 Department of Microbiology, Jashore University of Science and Technology, Jashore, Bangladesh
4 Department of Microbiology, Dhaka University, Dhaka, Bangladesh
5 Abhaynagar Upazila Health Complex, Jashore, Bangladesh
6 Department of Nutrition and Food Technology, Jashore University of Science and Technology, Jashore, Bangladesh
7 Department of Climate and Disaster Management, Jashore University of Science and Technology, Jashore, Bangladesh
8 Genome Centre, Jashore University of Science and Technology, Jashore, Bangladesh
M Anwar Hossain, Genome Centre, Jashore University of Science and Technology, Jashore 7408, Bangladesh. Email: hossaina@du.ac.bd
Hassan M Al-Emran, Department of Biomedical Engineering, Jashore University of Science and Technology, Jashore 7408, Bangladesh. Email: al-emran@just.edu.bd
11 9 2024
2024
18 117793222412723998 4 2024
29 6 2024
© The Author(s) 2024
2024
SAGE Publications Ltd unless otherwise noted. Manuscript content on this site is licensed under Creative Commons Licenses
https://creativecommons.org/licenses/by-nc/4.0/ This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License (https://creativecommons.org/licenses/by-nc/4.0/) which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/open-access-at-sage).
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) that emerged in late 2019 has accumulated a series of point mutations and evolved into several variants of concern (VOCs), some of which are more transmissible and potentially more severe than the original strain. The most notable VOCs are Alpha, Beta, Gamma, Delta, and Omicron, which have spread to various parts of the world. This study conducted surveillance in Jashore, Bangladesh to identify the prevalence of SARS-CoV-2 coinfected with dengue virus and their genomic effect on the emergence of VOCs. A hospital-based COVID-19 surveillance from June to August, 2021 identified 9 453 positive patients in the surveillance area. The study enrolled 572 randomly selected COVID-19-positive patients, of which 11 (2%) had dengue viral coinfection. Whole genome sequences of SARS-CoV-2 were analyzed and compared between coinfection positive and negative group. In addition, we extracted 185 genome sequences from GISAID to investigate the cross-correlation function between SARS-CoV-2 mutations and VOC; multiple ARIMAX(p,d,q) models were developed to estimate the average number of amino acid (aa) substitution among different SARS-CoV-2 VOCs. The results of the study showed that the coinfection group had an average of 30.6 (±1.7) aa substitutions in SARS-CoV-2, whereas the dengue-negative COVID-19 group had that average of 25.6 (±1.8; P < .01). The coinfection group showed a significant difference of aa substitutions in open reading frame (ORF) and N-protein when compared to dengue-negative group (P = .03). Our ARIMAX models estimated that the emergence of SARS-CoV-2 variants Delta required additional 9 to 12 aa substitutions than Alpha, Beta, or Gamma variant. The emergence of Omicron accumulated additional 19 (95% confidence interval [CI]: 15.74, 21.95) aa substitution than Delta. Increased number of point mutations in SARS-CoV-2 genome identified from coinfected cases could be due to the compromised immune function of host and induced adaptability of pathogens during coinfections. As a result, new variants might be emerged when series of coinfection events occur during concurrent two epidemics.

SARS-CoV-2
dengue
coinfection
genome sequence
variants of concern
ARIMAX model
Jashore University of Science and Technology https://doi.org/10.13039/501100016172 JUST/Research Cell/Research Project/2020-21/FOET11 cover-dateJanuary-December 2024
typesetterts1
==== Body
pmcIntroduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus has been under constant observance due to its continuously emerging variants since December 2019. 1 The WHO-defined variants of concern (VOCs) until 2022 were Alpha (September 2020), Beta (September 2020), Gamma (December 2020), Delta (December 2020), Eta (December 2020), and Omicron (November 2021). 2 One of the most convincing hypotheses prevailing that the variants are originated from chronic infections in immunocompromised individuals where the virus is able to establish a persistent infection.3,4 Those VOCs showed better potential in terms of pathogenicity, virulence, and transmission rate and recent variants exhibit lower antibody neutralization sensitivity.5,6 In addition to SARS-CoV-2, other seasonal viruses such as influenza, dengue, and respiratory syncytial virus (RSV) continue to circulate in the population.7-9 Those coinfections or secondary infections are often associated with disease severity and mortality 10 although their frequency is not well known.

One study reported a high case fatality rate among COVID-19 and flu viral coinfections (r = 0.2). 7 The coinfection of dengue and COVID-19 is a matter of concern, especially in dengue-endemic regions, as it can pose serious challenge to the health system. 8 ,11-15

Although the clinical manifestations of dengue,16,17 and COVID-19 are similar in some patients, the diagnosis may not be affected by dengue virus. However, the generation of cross-reactive antibodies during dengue virus infection18,19 influence the antibody-dependent enhancement and make the cases extremely complicated.14,20,21 Analyzing the whole-genome sequences of SARS-CoV-2 among the Dengue virus coinfection cases can help in understanding the potential effects of coinfections on the evolution of SARS-CoV-2 virus. However, such studies are still in their early stages, and more research is needed to establish any link between coinfection and genomic mutations. 22 The objective of this study is to explore the whole-genome sequences of SARS-CoV-2 among the dengue virus coinfection cases, SARS-CoV-2 mutation analysis, and predict their correlation with coinfection and emerging variants.

Methods

Sample size calculation

Proportion was used to calculate the minimum sample size using the adjusted formula. 23 We have considered a proportion of the indicators as 0.5 (P = 0.5) to give a conservative estimate of the sample size, 5% two-sided statistical level of significance (Z0.05/2 = 1.96, 5% margin of error (e = 0.05) and corrected size of sample for the nonresponse as well as large population size using design effect (DEFF = 1.50) and Cochran’s corrected sample size determination formula (equation 2), respectively, which results in a sample size of 544

(1) n0=Zα22×p^(1−p^)e2×DEFF

The adjusted sample size formula is given as

(2) n=n01+(n0−1N)

Study population

The study was conducted during the COVID-19 national surveillance at Genome Centre, Jashore University of Science and Technology (JUST), Bangladesh with the left-over samples. 10 Genome Center JUST tested a total of 27 993 samples from June 1, 2021 to August 31, 2021 and reported 9 453 (34%) COVID-19-positive cases in Jashore, Bangladesh. 24 Randomly selected 750 cases (8%, N = 9 453) were included for the verbal interview to reach the estimated sample size. The inclusion criteria in this coinfection study were SARS-CoV-2-positive cases diagnosed by Novel Coronavirus Nucleic Acid Diagnostic Kit (Sansure Biotech Inc., China) following the manufacturer’s instruction. The cases were excluded from the study if the patients are non-responsive over phone call or unable to provide the informed consent.

Spatial distribution mapping

Spatial distributions of COVID-19 cases were plotted using the ArcGIS (version 10.6) software. The approximate geolocations of the COVID-19-positive cases were pointed in Google Earth Pro platform.

Whole-genome sequencing

The SARS-CoV-2 genome sequencing was performed on samples from dengue coinfection–positive and age-, sex-, and time-matched dengue-negative cases. The sequencing procedure were performed as described elsewhere. 2 > In brief, 280 µL of nasopharyngeal samples were processed for viral RNA extraction using QIAamp Viral RNA Mini Kit according to manufacturer’s instruction. First-strand cDNAs were from the extracted viral RNA. The Ion AmpliSeq™ SARSCoV-2 Research Panel (Thermo Fisher Scientific, USA) was used to prepare libraries and target amplification. Each amplified sample was partially digested with FuPa reagents and ligated with Ion P1 adapter and Ion Xpress™ Barcode Adapters 1–16 kit (Ion Torrent™, Thermo Fisher Scientific, USA). All equimolar libraries were pooled for the preparation of template-positive Ion Sphere™ Particles (ISPs) using Ion 530™ Kit—OT2 (Thermo Fisher Scientific, USA) on the Ion One Touch™ 2 System and finally sequenced in Ion S5™ System.

Statistical analysis

Descriptive analysis was conducted from the study patients using frequency distribution to explore the sociodemographic profile, COVID-19-related experience, and information about hospitalization. The coinfection cases were observed as rare events and provided a smaller size of a sample (11 out of 572 COVID-19-positive cases). In this circumstance, this study adopted non-parametric test procedure to obtain trustworthy inferences. 25 This study employed non-parametric Chi-square test with Yates continuity correction to compare respondents’ characteristics and clinical features between COVID-19-positive and coinfected groups. In addition, the study employed non-parametric Mann-Whitney U-test to compare the mean amino acid (aa) substitution between COVID-19-positive and coinfected groups. The non-parametric Kruskal-Wallis H test were performed to identify the equality of mean mutations including pairwise variant comparison by protein change.

Time-series analysis of aa substitution

A total of 185 SARS-COV-2 whole-genome data were extracted from the GISAID database from January 2021 to October 2022. Each month we collected five to seven high-quality genome sequences mostly from Bangladeshi COVID-19 patients. We organized the genomic information based on their WHO-defined VOC and nonsynonymous mutations. We developed an ARIMA model of nonsynonymous mutations in S, open reading frame (ORF), N, M, E protein, and overall protein change applying usual procedure of identification, estimation, and diagnostics. 26 It is to be noted 27 that the number of mutations depends on the previous number of mutations as well as the previous VOC, which exhibits a nonstationary time series. This method involves successive differentiation and/or transformation to make the series stationary in order to determine the order of integration (d), as well as the development of autocorrelation and partial autocorrelation functions to determine the autoregressive (p) and moving average (q) order of the original series of SARS-CoV-2 mutation. In addition, the cross-correlation function was incorporated to identify the significant lag (l) of WHO-defined VOCs to explain the aa substitution. Let Mt and Ct be the integrated (if different) or non-integrated series of mutation and VOC respectively; thus, the ARIMAX (p, d, q) model can be written as

(3) Mt=α1Mt−1+α2Mt−2+⋯+αpMt−p+εt+θ1εt−1+θ2εt−2+⋯+θqεt−q+βlCt−l

where ε t is the white noise and α, θ, and β are the parameters of AR process, MA process, and exogenous variables. 27 The Bayesian information criteria were employed to identify the optimal model with the correct combination of p, d, and q order; accordingly, the coefficient of exogenous variable estimates the average number of aa substitutions for a shift to VOC groups. It should be noted that the residual assumptions of autocorrelation, normality, and stationarity were tested after the estimation of model parameters.

Results

The prevalence of COVID-19-positive cases were 34% (9 453 out of 27 993) from June 2021 to August 2021 in Jashore, Bangladesh. This study included 750 patients out of 9 453 positive cases, of which 153 cases (20.4%) were nonresponsive in phone call and 25 personnel (3.3%) did not provide verbal consent. The remaining 572 cases provided the verbal consent for the interview and were included in the study. The study found 11 cases of dengue-COVID-19 coinfection (1.9%) out of 572 COVID-19-positive cases in the study area.

Demographic features of the patients

The male-female ratio of the COVID-19-positive patients was 271 (47.4%): 301 (52.6%). The coinfection was more common among the male patients (P = .089). Among the dengue-COVID-19 coinfection-positive cases, clinical manifestations appeared in 100% of cases compared to 87% of negative cases. Coinfected patients were more likely to experience dry cough, loss of smell and taste, and headache (P < .05). Fever appeared in 100% coinfected patients compared to 78% of the dengue-negative cases, and muscle pain was more common among the coinfected cases (45% vs 27%). There was no reported death among the coinfected cases, whereas 4 COVID-19-positive cases died out of 561 (0.7%; see Table 1).

Table 1. Descriptive statistical analysis of the study-enrolled patients.

Variables	Dengue results	χ2-Statistics (P-value)	
Positive	Negative	
Sex	% (n)	% (n)		
 Male	72.7 (8)	46.9 (263)	1.947 (.163)	
 Female	27.3 (3)	53.1 (298)		
Symptoms appearance				
 Appeared	100.0 (11)	87.0 (488)	0.680 (.409)	
 Not appeared/cannot remember	—	13.0 (73)		
Appeared symptoms (multiple response)		
 Fever	100.0 (11)	78.1 (438)	1.911 (.167)	
 Dry cough	72.7 (8)	28.2 (158)	8.350 (.004)	
 Running nose	9.1 (1)	29.6 (166)	—	
 Muscle pain	45.5 (5)	27.1 (152)	1.021 (.312)	
 Difficulty/short breathing	45.5 (5)	39.0 (219)	0.014 (.905)	
 Loss of smell	90.9 (10)	23.2 (130)	23.239 (< .01)	
 Loss of taste	90.9 (10)	28.5 (160)	17.228 (< .01)	
 Loss of appetite	18.2 (2)	19.4 (109)	—	
 Sore throat	36.4 (4)	32.8 (184)	< 0.01 (.999)	
 Headache	63.6 (7)	6.8 (38)	40.601 (< .01)	
 Itching	—	3.0 (17)	—	
 Diarrhea	9.1 (1)	1.2 (7)	—	
 Blurred vision	—	2.3 (13)	—	
 Pain or swelling of legs and hands	9.1 (1)	5.0 (28)	—	
Infection outcome				
 Completely recovered	81.8 (9)	89.1 (500)	0.079 (.779)	
 Still weak	18.2 (2)	10.2 (57)		
 Death	—	0.7 (4)		

Spatial distribution of COVID-19 and dengue

Coronavirus disease 2019 was most prevalent in Avoynagar (N = 305) and Jashore Sadar (N = 85) during the study period, where the Dengue-COVID-19 coinfections rates were 2.9% (n = 9) and 2.3% (n = 2), respectively (see Figure 1).

Figure 1. Geospatial location of Dengue-COVID-19-coinfected cases in Jashore, Bangladesh from June-August 2021.

In Vivo whole-genome sequence analysis

Seven out of 11 Dengue-COVID-19-coinfected patients and five out of nine dengue-negative COVID-19-positive patients were selected for SARS-CoV-2 genomic analysis, which provided greater than 29 500 nucleotide sequence data (see Table 2). Remaining eight samples provided only 20 000 to 26 000 genomic data and were not considered for the mutation analysis. The average number of total nucleotide substitutions among the coinfected cases were 40.9 (±2.8) compared to 34.2 (±1.7) among COVID-19 cases (P < .001). The average differences of aa substitutions were 30.6 (±1.7) among the coinfected group and 25.6 (±1.8) among the latter group. Non-parametric Mann-Whitney U-test found significant differences of aa substitution in coinfected group (P < .01). The ORF (P = .003) and N-protein (P = .03) also found significant differences between the groups (see Figure 2).

Table 2. Genomic mutation analysis of Dengue-COVID-19 coinfection cases compared to COVID-19 cases.

ID	COVID-19 positive with	clade	Nextclade _pango	S-protein mutations	ORF-protein mutations	N-protein mutations	All non-syn. mutations	Unique mutations	
DC02	Dengue negative	21J (Delta)	B.1.617.2	8	14	3	26	0	
DC05	Dengue negative	21A (Delta)	B.1.617.2	7	14	4	26	0	
DC06	Dengue negative	21J (Delta)	B.1.617.2	8	14	5	28	2(Spike_F92S, M_K205R)	
DC13	Dengue negative	21A (Delta)	B.1.617.2	7	11	4	23	0	
DC14	Dengue negative	21J (Delta)	AY.122	7	14	3	25	0	
DP01	Dengue positive	21J (Delta)	B.1.617.2	9	17	4	31	0	
DP03	Dengue positive	21A (Delta)	B.1.617.2	7	19	5	32	0	
DP04	Dengue positive	21A (Delta)	B.1.617.2	7	15	5	28	0	
DP05	Dengue positive	21A (Delta)	B.1.617.2	7	16	5	29	0	
DP06	Dengue positive	21A (Delta)	B.1.617.2	7	17	5	30	0	
DP07	Dengue positive	21A (Delta)	B.1.617.2	8	17	5	31	0	
DP08	Dengue positive	21A (Delta)	B.1.617.2	9	17	6	33	1 (Spike_F329C)	
COVID-19, coronavirus disease 2019; ORF, open reading frame.

Figure 2. Non-parametric Mann-Whitney U-test results indicate significant differences of amino acid substitution in SARS CoV-2 proteins between Dengue-COVID-19-coinfected and dengue-negative COVID-19 patients. All samples were collected in Jashore, Bangladesh.

Protein heat map

The heat map analysis presented the changes of aa in the respective protein between coinfection-positive and -negative cases. Membrane proteins had one aa substitution, while nucleocapsid proteins had seven, non-structural proteins (ORF) had 45, and spike proteins had 15. Among the non-structural proteins, NS7a_L116 F, NSP12_P323L, NSP2_A318V, and NSP3_H1274Y were more commonly found in coinfected cases. On the contrary, NSP14_A394V, NSP3_A488S, NSP3_P1228L, NSP4_T492I, and NSP6_T77A were commonly found in dengue-negative cases. In spike proteins, few aa substitutions were observed in five strains (see Figure 3).

Figure 3. WGS analysis of SARS CoV-2 generated heat map of amino acid substitutions among dengue-positive and -negative groups.

In silico time-series analysis of mean aa substitution

From 185 SARS-CoV-2 genome sequence extracted from GISAID, 21 were Alpha, 20 were Beta, 21 were Gamma, 59 were Delta, 20 were Eta, and 44 were Omicron (see Figure 4). Mutation analysis revealed that aa substitutions in spike and ORF protein was increasing continuously (see Figure 5). The non-parametric Kruskal-Wallis H test results identify different pairs of VOCs in which the mean aa substitution differs significantly. The aa substitution in Alpha, Beta, Gamma, and Eta differs significantly than Delta and Omicron. Omicron variants also had significant aa substitutions (see Table 3).

Figure 4. Time-series plot with frequency of amino acid substitution in WHO-defined SARS CoV-2 variants of concern (VOC) detected in Bangladesh. All SARS CoV-2 genomic data were extracted from GISAID.

Figure 5. Frequency of mutations in SARS-CoV-2 protein types from January 2021 to October 2022. All SARS CoV-2 genomic data were extracted from GISAID.

Table 3. Mean ± standard deviation of aa substitution in SARS-CoV-2 strains and Kruskal-Wallis H test for the equality of means including pairwise comparison by protein change.

WHO VOC	S-protein	ORF-protein	N-protein	M-protein	E-protein	Overall	
Alpha	7.56a ± 0.62	9.57a ± 1.43	4.10a ± 0.30	1.00a ± < 0.01	—	20.24a ± 2.59	
Beta	8.65b ± 0.99	10.45a ± 1.28	1.15b ± 0.37	—	1.00 ± < 0.00	21.25ab ± 1.89	
Gamma	12.05c ± 1.02	7.95b ± 1.32	2.86c ± 0.85	—	—	22.86b ± 2.43	
Delta	9.44b ± 2.96	17.00c ± 3.25	3.76d ± 0.70	1.05a ± 0.29	1.00 ± < 0.01	31.31c ± 4.34	
Eta	6.05d ± 1.00	5.85d ± 2.32	2.70ce ± 0.57	1.00a ± < 0.01	1.00 ± < 0.01	15.70d ± 2.75	
Omicron	26.02e ± 5.99	16.84c ± 3.20	4.09af ± 0.88	2.14b ± 0.67	1.20 ± 0.40	50.20e ± 8.50	
H-stat (P-value)	168.06 (< .01)	99.78 (< .01)	63.89 (< .01)	55.58 (< .01)	2.13 (.105)	204.62 (< .01)	
Overall	13.10 ± 8.22	13.18 ± 5.10	3.38 ± 1.15	1.43 ± 0.69	1.12 ± 0.32	30.81 ± 13.04	
aa, amino acid; ORF, open reading frame; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2; VOC, variant of concern; WHO, World Health Organization.

Means within a column followed by different superscripted letters are statistically significantly (P < .05) different based on non-parametric pairwise comparison.

Our time-series analysis found the parameters estimates of the exogenous variables with 95% confidence interval (CI) for each aa substitution (see Table 4). According to the estimated results, emergence of SARS-CoV-2 variants from Alpha to Delta or Omicron required an average of 11.1 or 29.5 additional mutations, respectively. An average of 10.1 and 28.3 additional mutations were required from Beta variant to Delta and Omicron, respectively. From Gamma variant, Delta and Omicron accumulated an average of 8.5 and 27.0 mutations, respectively. From Delta or Eta to Omicron, 18.9 and 33.9 additional mutations took place (see Figure 6). Average mutations with 95% CI in protein category (S, ORF, N, M, and E proteins) were also mentioned in Table 4.

Table 4. Estimated parameters, 95% confidence interval of the exogenous variables in the estimated ARIMAX model.

Reference	Explanatory variable	S protein	ORF protein	N protein	M protein	E protein	Overall	
Alpha	Alpha to Beta	1.13*** (0.75, 1.50)	0.97 (–1.15, 3.08)	–3.00*** (–3.98, –2.02)			1.01 (–0.35, 2.37)	
	Alpha to Gamma	4.00*** (1.62, 6.38)	–1.62*** (–2.44, –0.80)	–1.22*** (–1.70, –0.75)			2.62*** (1.12, 4.12)	
	Alpha to Delta	0.72 (–1.37, 2.81)	6.53*** (3.34, 9.71)	–0.33 (–0.64, –0.03)	0.05, (–0.35, 0.45)		11.07*** (9.10, 13.03)	
	Alpha to Eta	–1.10 (–2.19, –0.01)	–2.15 (–4.63, 0.34)	–1.40*** (–1.67, –1.12)			–4.54*** (–6.15, –2.92)	
	Alpha to Omicron	10.05*** (2.98, 17.11)	8.38*** (4.51, 12.24)	0.01 (–0.48, 0.49)	0.10 (–0.83, 1.03)		29.46*** (24.82, 34.09)	
Beta	Beta to Gamma	3.40*** (2.79, 4.01)	–2.50*** (–3.28, –1.71)	2.00*** (0.42, 3.58)			1.66*** (0.50, 2.83)	
	Beta to Delta	–0.11 (–2.31, 2.10)	5.56*** (2.53, 8.77)	2.61*** (2.29, 2.93)			10.06*** (8.10, 12.01)	
	Beta to Eta	–2.59*** (–2.98, –2.21)	–4.60*** (–5.70, –3.50)	1.55*** (1.26, 1.84)			–5.55*** (–6.99, –4.11)	
	Beta to Omicron	10.12** (1.46, 18.79)	7.52*** (3.71, 11.34)	3.00*** (1.42, 4.58)		–0.02 (–0.35, 0.30)	28.30*** (23.31, 33.30)	
Gamma	Gamma to Delta	–2.62*** (–3.94, –1.30)	8.31*** (5.38, 11.24)	0.91*** (0.54, 1.27)			8.45*** (6.50, 10.40)	
	Gamma to Eta	–6.00*** (–6.61, –5.39)	–2.10*** (–3.24, –0.97)	–0.16 (–0.59, 0.28)			–7.19*** (–8.34, –6.05)	
	Gamma to Omicron	9.07*** (3.32, 14.82)	9.80*** (6.07, 13.53)	1.00 (–0.82, 2.82)			26.98*** (22.16, 31.80)	
Delta	Delta to Eta	–3.37*** (–4.71, –2.03)	–11.47***(–13.65, –9.29)	–1.06*** (–1.40, –0.72)	–0.05 (–0.19, 0.08)		–15.61*** (–17.6, –13.4)	
	Delta to Omicron	10.29*** (5.31, 15.26)	–0.05 (–2.03, 1.92)	0.33 (–0.07, 0.73)	0.10 (–0.53, 0.74)	–0.04 (–0.26, 0.18)	18.85*** (15.74, 21.95)	
Eta	Eta to Omicron	14.89*** (9.05, 20.73)	11.08*** (7.28, 14.88)	1.07** (0.24, 1.90)	0.10 (–0.66, 0.86)	–0.04 (–0.26, 0.18)	33.92*** (28.97, 38.88)	
ORF, open reading frame.

All SARS CoV-2 genomic data were extracted from GISAID.

*** and ** Indicates that the relationship between WHO clade category and SARS-CoV-2 genomic mutation is statistically significant at 1% and 5% level, respectively.

Figure 6. Coefficients of variants with 95% confidence interval among the different variants of concern. The differences were plotted based on non-synonymous mutations in SARS-CoV-2. Significant (P < .01) changes were highlighted in red marks.

Discussion

Dengue infection can cause a wide range of clinical manifestations, from asymptomatic or mild fever to potentially fatal diseases, such as dengue hemorrhagic fever or dengue shock syndrome. The immune response to dengue infection involves B cells, T cells, and interferons, but the virus uses various strategies to evade this response. 28 In response to dengue infection, cytokines such as interferon alpha beta (IFN-αβ), interleukin 6 (IL-6), IL-8, migration inhibitory factor (MIF), tumor necrosis factor alpha (TNF-α), IFN-γ, IL-10, and vascular endothelial growth factor (VEGF) can regulate viral replication and direct immune cells to the infection site. If this response fails, viremia develops and an amplified cytokine response leads to disease severity and affects other organs. 29 Coinfection with SARS-CoV-2 manifested severe symptoms in our study patients. We observed significant changes in clinical manifestations, specifically headache, dry cough, loss of taste and smell. The clinical outcome of the patient remained similar to some symptoms except for weakness, fever, and muscle pain (see Table 1). However, we did not perform the non-parametric chi-square test in symptoms that had low respondents. Furthermore, we identified low number of reported death cases (0% vs 0.7%) like other studies regarding the fatality in coinfected group.8,30

The presence of dengue virus coinfection with SARS-CoV-2 may have significant implications on the severity and outcome of the disease. The hypercoagulable condition resulting from high cytokine levels and ACE2 downregulation, which are typical of SARS-CoV-2, can exacerbate the interruption of the fibrinolysis system caused by dengue, potentially leading to hemoptysis. 31 The ACE2 receptors are also present in liver. 32 The impact of dengue results in fatty liver and acute liver failure. The whole spectrum of liver illness caused by dengue can progress from an asymptomatic increase in hepatic transaminases to the manifestation of abrupt liver failure. Between 4% and 52% of adult dengue victims get liver enlargement. Raised liver function results in persistent abnormality among 45% to 96% patients. 33 Furthermore, the use of ACE2 as a receptor by SARS-CoV- 2 may cause severe damage to the liver, which is already weakened by dengue. A weaker liver, elevated cytokine levels, and downregulated ACE2 all contribute to a suppressed immune system, which can prolong the duration of the virus in the body. Through its interactions with the renin-angiotensin system (RAS), ACE2 controls the immunological response.34-36 This suppressed immune system may also facilitate the longer survival and replication of other viruses, including SARS-CoV-2, leading to the emergence of new variants. This study found that the mean aa substitution in SARS-CoV-2 viruses was significantly higher among patients infected with both dengue and SARS-CoV-2 viruses compared to those who were infected with SARS-CoV-2 virus only (P < .01; see Figure 2). This finding suggests that coinfection may be paving the way for the emergence of new variants of the virus, which could have implications for the effectiveness of the current treatments. 37

There are significant differences in the aa sequences of various proteins between coinfection-positive and -negative cases. The changes were observed in the membrane, nucleocapsid, nonstructural (ORF), and spike proteins (see Table 2). The aa substitutions found between the groups exhibited some differences in several proteins such as NS7a (ORF7a), NSP1, NSP2, NSP3, NSP4, NSP6, NSP12, NSP14, and spike proteins (see Figure 3). Most of them are for replication and adaptation of the virus into a new host. 38 ,39 The SARS-CoV-2 ORF7a has immunomodulatory function and it can counteract the antiviral effect the host cell has on the budding virions. The NS7a_L116 F mutation resides in the end-loop of the C-terminal domain. The “indel” mutation in ORF7a usually causes destabilization of protein structure and thus can enhance hindrance from human immune response. 40 This particular mutation seen in our study sequences might be able to affect viral replication and thus can lead to antigenic drift. 41 The NSP2_A318V mutation was reported in Malaysia and other countries, but no specific effect of the mutations has been reported yet.42-44 The NSP3_H1274Y mutation was reported in Pakistan 44 and the Australian Capital Territory (ACT) region at high levels. 45 This mutation resides in another important deletion mutation (S1265-(Serine→deletion) and L1266I (Leucine→Isoleucine) in the βSM region) reported in a recent study. 46 This disordered βSM region has, however, no definite functional role. The NSP12_P323L, coevolving D614G and P323L mutations in SARS-CoV-2 are associated with the severity of COVID-19, 47 infectivity and survival ability of the viral strain. 48

This study has a potential limitation where the prevalence (1.9%) of coinfection may have been underestimated due to only clinically suspected cases being recommended for the dengue virus testing with NS1 antigen and IgG antibodies. Moreover, we did not consider any case as coinfection positive if any patient was unable to show both positive test records. It is also important to note that a large number of selected cases were nonresponsive in phone calls and did not provide the consent, which may introduce some bias into the sample (see Figure 1).

The study used an ARIMAX model from the GISAID repository to estimate effect of point mutations in SARS-CoV-2 variants over time. We designed the model in each protein change by keeping each VOC category as a reference to avoid multi-collinearity. The order of ARIMAX model with minimum Bayesian information criteria including smaller residual variance from the indication of autocorrelation and partial autocorrelation function of the stationary series of mutation (see Table 5). ARIMAX (2,1,0) model in S-Protein for Alpha to Omicron indicates that the later variant might be evolved in the time of lag 2; more specifically, previous two time points. ARIMAX (0,1,1) indicates that the mutation is related with the errors of lag 1 and ARIMAX (0,0,0) indicates the white noise process (see Table 5). However, we were interested to estimate the coefficient of exogenous variables using proper ARIMAX model instead of forecasting out-of-sample mutations. The negative coefficient with significant effect of corresponding exogenous variables indicated that Eta did not emerge from Alpha, Beta, Gamma, or Delta variants. It could be generated from the wild-type variant.

Table 5. The autoregressive (p) and moving average (q) order of the integrated (d) or non-integrated (d = 0) series of mutation for each protein change with reference of WHO-defined variants of concern category.

Reference	Exogenous variable	S protein	ORF protein	N protein	M protein	E protein	Overall	
Alpha	Alpha to Beta	ARIMAX (0,1,2)	ARIMAX (2,1,0)	ARIMAX (0,1,0)			ARIMAX (0,0,0)	
	Alpha to Gamma	ARIMAX (0,1,0)	ARIMAX (0,1,1)	ARIMAX (1,0,0)			ARIMAX (0,1,1)	
	Alpha to Delta	ARIMAX (4,1,1)	ARIMAX (0,1,1)	ARIMAX (0,0,0)	ARIMAX (0,0,0)		ARIMAX (0,1,1)	
	Alpha to Eta	ARIMAX (3,1,0)	ARIMAX (2,1,0)	ARIMAX (0,1,1)			ARIMAX (0,1,1)	
	Alpha to Omicron	ARIMAX (2,1,0)	ARIMAX (0,1,1)	ARIMAX (0,0,1)	ARIMAX (0,1,1)		ARIMAX (0,1,3)	
Beta	Beta to Gamma	ARIMAX (0,1,1)	ARIMAX (0,1,1)	ARIMAX (0,1,0)			ARIMAX (1,1,2)	
	Beta to Delta	ARIMAX (0,1,1)	ARIMAX (0,1,1)	ARIMAX (0,1,1)			ARIMAX (0,1,1)	
	Beta to Eta	ARIMAX (0,1,2)	ARIMAX (0,1,2)	ARIMAX (0,1,1)			ARIMAX (0,1,1)	
	Beta to Omicron	ARIMAX (1,1,1)	ARIMAX (0,1,1)	ARIMAX (0,1,0)		ARIMAX (2,1,1)	ARIMAX (0,1,2)	
Gamma	Gamma to Delta	ARIMAX (1,0,1)	ARIMAX (0,1,1)	ARIMAX (0,1,1)			ARIMAX (0,1,1)	
	Gamma to Eta	ARIMAX (0,1,1)	ARIMAX (0,1,1)	ARIMAX (0,0,0)			ARIMAX (0,1,2)	
	Gamma to Omicron	ARIMAX (0,1,1)	ARIMAX (0,1,1)	ARIMAX (0,1,0)			ARIMAX (0,1,2)	
Delta	Delta to Eta	ARIMAX (1,0,1)	ARIMAX (0,1,1)	ARIMAX (0,1,1)	ARIMAX (0,0,0)		ARIMAX (0,1,1)	
	Delta to Omicron	ARIMAX (0,1,1)	ARIMAX (1,0,1)	ARIMAX (1,1,1)	ARIMAX (0,1,1)	ARIMAX (2,1,2)	ARIMAX (0,1,2)	
Eta	Eta to Omicron	ARIMAX (2,1,0)	ARIMAX (0,1,1)	ARIMAX (0,1,2)	ARIMAX (0,1,1)	ARIMAX (2,1,2)	ARIMAX (0,1,2)	
ORF, open reading frame.

The in silico analysis of SARS-CoV-2 sequences found that the frequency of aa substitution in Spike and ORF proteins increased over time (see Figure 4). Our time-series plot indicated that the emergence of new VOC was mostly due to the mutations in spike protein and ORF protein (see Figure 5). The study also found that the Delta variants mostly emerged from changes in ORF proteins, while the Omicron variants emerged with changes in both Spike and ORF proteins. Our dengue and COVID-19 coinfection study was placed from June to August of 2021 when the Delta variant was predominant in Bangladesh. The study estimated that the SARS-CoV-2 Omicron variant accumulated an average of 18.9 (95% CI: 15.74, 21.95) non-synonymous mutations than Delta. Our hospital-based surveillance observed that one coinfection cycle caused five non-synonymous mutations on top (mean 25.6 vs 30.6) and suggested that existing COVID-19 pandemics might come up with new variants of concern soon. 49 As for the potential emergence of new variants at the end of the dengue season, it is indeed possible that the seasonal spread of dengue could infect along with the SARS-CoV-2 and evolve as a new variant. One of our studies is reporting that the 2023 dengue outbreak had a higher morbidity rate compared to previous years (2019-2022) in Southwest of Bangladesh. 50 That outbreak was caused by DENV-2 Cosmopolitan genotype that was evolved as a new variant among C clade. The outbreaks of 2017 and 2018 in Bangladesh were also caused by DENV-2 Cosmopolitan genotype within B and C clade, 51 although the 2023 variant did not resemble any previous Bangladeshi C clade variants (2017, 2018, and 2019), suggesting a shift of variants.

Conclusions

The SARS-CoV-2, despite having proof-reading activity, continues to produce new variants since 2020. The RdRp of SARS-CoV-2 can correct errors during replication, which contribute to the stability of the virus genome. Nonetheless, the emergence of new variants took place due to its high rate of replication and mutation of the virus, combined with selective pressure from the host immune responses and coinfection with other viruses. Although there was no recombination event evident in our study, our mutation analysis found a significant difference in the frequency of point mutations among the coinfected cases. The findings of our study noted that coinfection with two different viruses can potentially lead to the emergence of new variants. The emergence of new variants and their potential impact on the ongoing pandemic further highlights the need for research on coinfection and their potential effects on viral evolution and pathogenicity. Continuous monitoring and study of the evolution of the virus and its variants is required to predict and prepare for potential future outbreaks. This can help public health officials and policymakers develop effective measures to control the spread of the virus, including developing new vaccines and treatments.

We would like to acknowledge Genome center, Jashore University of Science and Technology and the Ministry of DG Health, Bangladesh Government.

Author Contributions: HMA contributed toward conceptualization, methodology, validation, formal analysis, investigation, visualization, writing original draft, funding acquisition, and project administration. FR contributed to software, formal analysis, and data curation. LS, PM, SY, PS, TA, RP, and MR were involved in the investigation. PKD was involved in the investigation, resources, validation, and formal analysis. MSH was involved in investigation, resources, and validation. MSR contributed to formal analysis and data curation. SLS contributed to investigation, resources, and validation. ASMRUA contributed to writing review and editing, resources. AH, OKI, MTI, SN, and SA contributed to writing review and editing. MMR contributed to writing review and editing, resources, validation, and formal analysis. IKJ contributed to formal analysis, validation, writing review and editing, resources, supervision, and project administration. MAH contributed to writing review and editing, supervision.

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding: The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The study was funded (JUST/Research Cell/Research Project/2020-21/FOET11) by Jashore University of Science and Technology, Jashore-7408, Bangladesh. This funding does not include any personal cost or salary.

Ethical Approval: The study has gone through the Ethical Review Board of Jashore University of Science and Technology (ERC no: ERC/FBST/JUST/2021-62) and approved the study. The participants have been informed about the purpose and the objectives of this study. Verbal consent was obtained from optimistic respondents before they were enrolled in this study, and written informed consent has been taken before the processing of the samples. We ensured that all ethical considerations have been taken into account while conducting the research involving human participants. This includes obtaining the informed consent, protecting participants’ privacy and confidentiality, and ensuring that the study was conducted in an ethical and responsible manner. The research was conducted in accordance with the applicable laws and regulation, and that any potential risks or ethical concerns were addressed and managed appropriately.

ORCID iD: Hassan M Al-Emran https://orcid.org/0000-0003-1185-6720

Data Availability Statement: The SARS-CoV-2 whole-genome sequences data used in this coinfection study were deposited in the GISAID database (accession no. EPI_ISL_17226262, EPI_ISL_17226263, EPI_ISL_17226264, EPI_ISL_17226265, EPI_ISL_17226266, EPI_ISL_17226267, EPI_ISL_17226268, EPI_ISL_17226269, EPI_ISL_17226270, EPI_ISL_17226271, EPI_ISL_17226272, and EPI_ISL_17229387).
==== Refs
References

1 Islam OK Al-Emran HM Hasan MS Anwar A Jahid MIK Hossain MA. Emergence of European and North American mutant variants of SARS-CoV-2 in South-East Asia. Transbound Emerg Dis. 2021;68 :824-832. doi:10.1111/tbed.13748 32701194
2 World Health Organization. WHO Announces Simple, Easy-to-Say Labels for SARS-CoV-2 Variants of Interest and Concern. World Health Organization; 2021. Accessed March 13, 2023. https://www.who.int/news/item/31-05-2021-who-announces-simple-easy-to-say-labels-for-sars-cov-2-variants-of-interest-and-concern
3 Alam ASMRU Islam OK Hasan MS , et al . Dominant clade-featured SARS-CoV-2 co-occurring mutations reveal plausible epistasis: an in silico based hypothetical model. J Med Virol. 2022;94 :1035-1049. doi:10.1002/jmv.27416 34676891
4 Carabelli AM Peacock TP Thorne LG , et al . SARS-CoV-2 variant biology: immune escape, transmission and fitness. Nat Rev Microbiol. 2023;21 :162-177. doi:10.1038/s41579-022-00841-7 36653446
5 Davies NG Kucharski AJ Eggo RM Gimma A Edmunds WJ ; Centre for the Mathematical Modelling of Infectious Diseases COVID-19 Working Group. Effects of non-pharmaceutical interventions on COVID-19 cases, deaths, and demand for hospital services in the UK: a modelling study. Lancet Public Health. 2020;5 :e375-e385. doi:10.1016/S2468-2667(20)30133-X
6 Wang W Xu Y Gao R , et al . Detection of SARS-CoV-2 in different types of clinical specimens. JAMA. 2020;323 :1843-1844. doi:10.1001/jama.2020.3786 32159775
7 Alosaimi B Naeem A Hamed ME , et al . Influenza co-infection associated with severity and mortality in COVID-19 patients. Virol J. 2021;18 :127. doi:10.1186/s12985-021-01594-0 34127006
8 Tsheten T Clements ACA Gray DJ Adhikary RK Wangdi K. Clinical features and outcomes of COVID-19 and dengue co-infection: a systematic review. BMC Infect Dis. 2021;21 :729. doi:10.1186/s12879-021-06409-9 34340682
9 Zhu X Ge Y Wu T , et al . Co-infection with respiratory pathogens among COVID-2019 cases. Virus Res. 2020;285 :198005. doi:10.1016/j.virusres.2020.198005 32408156
10 Al-Emran HM Rahman S Hasan MS , et al . Microbiome analysis revealing microbial interactions and secondary bacterial infections in COVID-19 patients comorbidly affected by Type 2 diabetes. J Med Virol. 2023;95 :e28234. doi:10.1002/jmv.28234
11 Senior Correspondent. Bangladesh reports five deaths from dengue fever, the highest for a day in 2022. bdnews24.com. September 6, 2022. Accessed March 13, 2023. https://bdnews24.com/health/wmgf7thteq
12 Hariadi P Lokida D Menur Naysilla A , et al . Coinfection with SARS-CoV-2 and dengue virus: a case report highlighting diagnostic challenges. Front Trop Dis. 2022;3 :801276. doi:10.3389/fitd.2022.801276
13 Khan S Akbar SMF Nishizono A. Co-existence of a pandemic (SARS-CoV-2) and an epidemic (Dengue virus) at some focal points in Southeast Asia: pathogenic importance, preparedness, and strategy of tackling. Lancet Reg Health Southeast Asia. 2022;4 :100046. doi:10.1016/j.lansea.2022.100046 35873345
14 Lorenz C Azevedo TS Chiaravalloti-Neto F. COVID-19 and dengue fever: a dangerous combination for the health system in Brazil. Travel Med Infect Dis. 2020;35 :101659. doi:10.1016/j.tmaid.2020.101659 32278756
15 Nasomsong W Luvira V Phiboonbanakit D. Case report: dengue and COVID-19 coinfection in Thailand. Am J Trop Med Hyg. 2021;104 :487-489. doi:10.4269/ajtmh.20-1340
16 Munoz-Jordan J Cardona J Beltrán M , et al . Evaluation of serologic cross-reactivity between dengue virus and SARS-CoV-2 in patients with acute febrile illness—United States and Puerto Rico, April 2020-March 2021. MMWR Morb Mortal Wkly Rep. 2022;71 :375-377. doi:10.15585/mmwr.mm7110a3 35271558
17 Yan G Lee CK Lam LTM , et al . Covert COVID-19 and false-positive dengue serology in Singapore. Lancet Infect Dis. 2020;20 :536. doi:10.1016/S1473-3099(20)30158-4 32145189
18 Harapan H Ryan M Yohan B , et al . Covid-19 and dengue: double punches for dengue-endemic countries in Asia. Rev Med Virol. 2021;31 :e2161. doi:10.1002/rmv.2161
19 Nath H Mallick A Roy S , et al . COVID-19 serum can be cross-reactive and neutralizing against the dengue virus, as observed by the dengue virus neutralization test. Int J Infect Dis. 2022;122 :576-584. doi:10.1016/j.ijid.2022.07.013 35811081
20 Haqqi A Awan UA Ali M Saqib MAN Ahmed H Afzal MS. COVID-19 and dengue virus coepidemics in Pakistan: a dangerous combination for an overburdened healthcare system. J Med Virol. 2021;93 :80-82. doi:10.1002/jmv.26144 32510175
21 Joob B Wiwanitkit V. Hemorrhagic problem among the patients with COVID-19: clinical summary of 41 Thai infected patients. Clin Appl Thromb Hemost. 2020;26 :1076029620918308. doi:10.1177/1076029620918308
22 Sharma D Notarte KI Fernandez RA Lippi G Gromiha MM Henry BM . In silico evaluation of the impact of Omicron variant of concern sublineage BA.4 and BA.5 on the sensitivity of RT-qPCR assays for SARS-CoV-2 detection using whole genome sequencing. J Med Virol. 2023;95 :e28241. doi:10.1002/jmv.28241
23 Heinisch O. Cochran W.G. : sampling techniques, 2. Aufl. John Wiley and Sons, New York, London 1963. Preis s. Biom Z. 1965;7 :203. doi:10.1002/BIMJ.19650070312
24 Al-Emran HM Hasan MS Ahasan Setu MA , et al . Genomic analysis of SARS-CoV-2 variants of concern identified from the ChAdOx1 nCoV-19 immunized patients from Southwest part of Bangladesh. J Infect Public Health. 2022;15 :156-163. doi:10.1016/j.jiph.2021.12.002 34952247
25 Brown L . Statistics and evidence-based healthcare. In: Bennett P Williamson C Sykes L MacIntyre DA Dixon PH , eds. Basic Science in Obstetrics and Gynaecology: A Textbook for MRCOG—Part 1. 4th ed. Elsevier; 2010 :289-304. doi:10.1016/B978-0-443-10281-3.00018-X
26 Zhang J Zhao Z Xue Y Chen Z Ma X Zhou Q. Time Series Analysis. 2017. In: ﻿Fang J, eds. Handbook of Medical Statistics; 2017:269-299. doi:10.1142/9789813148963_0009
27 Schroeder L de Souza EM Rosset C , et al . Fire association with respiratory disease and COVID-19 complications in the State of Pará, Brazil. Lancet Reg Health Am. 2022;6 :100102-100114. doi:10.1016/j.lana.2021.100102 34870262
28 Kao YT Lai MMC Yu CY . How dengue virus circumvents innate immunity. Front Immunol. 2018;9 :2860-2868. doi:10.3389/fimmu.2018.02860 30564245
29 Srikiatkhachorn A Mathew A Rothman AL. Immune-mediated cytokine storm and its role in severe dengue. Semin Immunopathol. 2017;39 :563-574. doi:10.1007/s00281-017-0625-1 28401256
30 Estofolete CF Machado LF Zini N , et al . Presentation of fatal stroke due to SARS-CoV-2 and dengue virus coinfection. J Med Virol. 2021;93 :1770-1775. doi:10.1002/JMV.26476 32881018
31 Amin MA Islam MT Khan II Khan ZH Qadri F Khan AI. COVID-19 and dengue infection in Bangladesh: a case of coinfection where hemoptysis as first presentation. Clin Case Rep. 2022;10 :e05252-05255. doi:10.1002/ccr3.5252
32 Leowattana W. Angiotensin-converting enzyme 2 receptors, chronic liver diseases, common medications, and clinical outcomes in coronavirus disease 2019 patients. World J Virol. 2021;10 :86-96. doi:10.5501/WJV.V10.I3.86 34079691
33 Venugopalan DRV. When dengue stings the liver. Published 2023. Accessed June 13, 2024. https://www.deccanherald.com/features/when-dengue-stings-the-liver-1182512.html
34 Marcinkiewicz J Witkowski JM Olszanecki R. The dual role of the immune system in the course of COVID-19. The fatal impact of the aging immune system. Cent Eur J Immunol. 2021;46 :1-9. doi:10.5114/CEJI.2021.105240 33897278
35 Qu L Chen C Yin T , et al . ACE2 and innate immunity in the regulation of SARS-CoV-2-induced acute lung injury: a review. Int J Mol Sci. 2021;22 :11483. doi:10.3390/IJMS222111483 34768911
36 Wang J Kaplan N Wysocki J , et al . The ACE2-deficient mouse: a model for a cytokine storm-driven inflammation. FASEB J. 2020;34 :10505-10515. doi:10.1096/FJ.202001020R 32725927
37 Al-Emran HM Rahman F Sarkar L , et al . Emergence of SARS-CoV-2 variants of concern induced by dengue viral coinfections. Accessed August 17, 2024. https://papers.ssrn.com/sol3/papers.cfm?abstract_id=4401701
38 Chen Y Liu Q Guo D. Emerging coronaviruses: genome structure, replication, and pathogenesis. J Med Virol. 2020;92 :418-423.31967327
39 Khailany RA Safdar M Ozaslan M. Genomic characterization of a novel SARS-CoV-2. Gene Rep. 2020;19 :100682. doi:10.1016/j.genrep.2020.100682 32300673
40 Nemudryi A Nemudraia A Wiegand T , et al . SARS-CoV-2 genomic surveillance identifies naturally occurring truncation of ORF7a that limits immune suppression. Cell Rep. 2021;35 :109197. doi:10.1016/j.celrep.2021.109197 34043946
41 Sahni C Basu Roy Chowdhury P Devadas D , et al . SARS-CoV-2 mutations responsible for immune evasion leading to breakthrough infection. Cureus. 2022;14 :e29544. doi:10.7759/cureus.29544
42 Tan KK Tan JY Wong JE , et al . Emergence of B.1.524(G) SARS-CoV-2 in Malaysia during the third COVID-19 epidemic wave. Sci Rep. 2021;11 :22105. doi:10.1038/s41598-021-01223-4 34764315
43 Patro LPP Sathyaseelan C Uttamrao PP Rathinavelan T . The evolving proteome of SARS-CoV-2 predominantly uses mutation combination strategy for survival. Comput Struct Biotechnol J. 2021;19 :3864-3875. doi:10.1016/j.csbj.2021.05.054 34109017
44 Umair M Ikram A Salman M , et al . Genomic surveillance reveals the detection of SARS-CoV-2 delta, beta, and gamma VOCs during the third wave in Pakistan. J Med Virol. 2022;94 :1115-1129.34726786
45 Tomaszewski T Gurtler V Caetano-Anollés K Caetano-Anollés G . Chapter 8: the emergence of SARS-CoV-2 variants of concern in Australia by haplotype coalescence reveals a continental link to COVID-19 seasonality. In: Pavia CS Gurtler V , eds. COVID-19: Biomedical Perspectives (Vol. 50 ). Academic Press; 2022 :233-268. doi:10.1016/bs.mim.2022.03.003
46 Hossain A Akter S Rashid AA Khair S Alam ASMRU . Unique mutations in SARS-CoV-2 Omicron subvariants’ non-spike proteins: potential impacts on viral pathogenesis and host immune evasion. Microb Pathog. 2022;170 :105699. doi:10.1016/j.micpath.2022.105699 35944840
47 Biswas SK Mudi SR. Spike protein D614G and RdRp P323L: the SARS-CoV-2 mutations associated with severity of COVID-19. Genomics Inform. 2020;18 :e44. doi:10.5808/GI.2020.18.4.e44
48 Ilmjärv S Abdul F Acosta-Gutiérrez S , et al . Concurrent mutations in RNA-dependent RNA polymerase and spike protein emerged as the epidemiologically most successful SARS-CoV-2 variant. Sci Rep. 2021;11 :13705. doi:10.1038/s41598-021-91662-w 34210996
49 Pilapil JD Notarte KI Yeung KL. The dominance of co-circulating SARS-CoV-2 variants in wastewater. Int J Hyg Environ Health. 2023;253 :114224. doi:10.1016/J.IJHEH.2023.114224 37523818
50 Setu MAA Das PK Ahammed T , et al . Dengue virus serotype 2 cosmopolitan C genotype reemerges with a new subclade in southwest region of Bangladesh. medRxiv. 2023.
51 Hossain MS Noman A Al Mamun SMAA Mosabbir AA. Twenty-two years of dengue outbreaks in Bangladesh: epidemiology, clinical spectrum, serotypes, and future disease risks. Trop Med Health. 2023;51 :37. doi:10.1186/s41182-023-00528-6 37434247
