Comparative phylogenetic analysis and transcriptomic profiling of Dengue (DENV-3 genotype I) outbreak in 2021 in Bangladesh

Background The next-generation sequencing (NGS) technology facilitates in-depth study of host-pathogen metatranscriptome. We, therefore, implicated phylodynamic and transcriptomic approaches through NGS technology to know/understand the dengue virus (DENV) origin and host response with dengue fever. Methods In this study, blood serum RNA was extracted from 21 dengue patients and 3 healthy individuals. Total transcriptomic data were analyzed for phylogenetic, phylodynamic, differential express gene (DEG), and gene ontology (GO) using respective bioinformatics tools. Results The viral genome sequence revealed dengue viral genome size ranges 10647 to 10707 nucleotide. Phylogenetic and phylodynamic analysis showed that the 2021 epidemic isolates were DENV-3 genotype-I and maintained as a new clade in compared to 2019 epidemic. Transcriptome analysis showed a total of 2686 genes were DEG in dengue patients compared to control with a q-value < 0.05. DESeq2 plot counts function of the top 24 genes with the smallest q-values of differential gene expression of RNA-seq data showed that 11 genes were upregulated, whereas 13 genes were downregulated. GO analysis showed a significant upregulation (p = < 0.001) in a process of multicellular organismal, nervous system, sensory perception of chemical stimulus, and G protein-coupled receptor signaling pathways in the dengue patients. However, there were a significant downregulation (p = < 0.001) of intracellular component, cellular anatomical entity, and protein-containing complex in dengue patients. Most importantly, there was a significant increase of a class of immunoregulatory proteins in dengue patients in compared to the controls, with increased GO of immune system process. In addition, upregulation of toll receptor (TLR) signaling pathways were found in dengue patients. These TLR pathways were particularly involved for the activation of innate system coupled with adaptive immune system that probably involved the rapid elimination of dengue virus infected cells. These differentially expressed genes could be further investigated for target based prophylactic interventions for dengue. Conclusion This is a first report describing DENV complete genomic features and differentially expressed genes in patients in Bangladesh. These genes may have diagnostic and therapeutic values for dengue infection. Continual genomic surveillance is required to further investigate the shift in dominant genotypes in relation to viral pathogenesis. Supplementary Information The online version contains supplementary material available at 10.1186/s12985-023-02030-1.


Introduction
The dengue virus is transmitted by mosquitoes and causes a global outbreaks and epidemics in tropical and subtropical areas, with nearly 400 million cases annually, posing an immediate threat to human health in developing countries like Bangladesh [1].Though first dengue report in Bangladesh was at 1964 but the outbreaks began in 2000 [2] caused by DENV-3 [3] and continuing with circulating serotype/genotype.In 2000-2009 DENV-3 genotype II were prevalent but 2013-2016 DENV-1 and DENV-2 were the most prevalent serotypes in Bangladesh [4].The serotype DENV-3 was predominant in outbreaks in 2017-2019 and reported to have clad shift to genotype I [5].The World Health Organization (WHO) reports that 500,000 cases of severe dengue each year require hospitalization, mainly among children [6].Clinically, dengue virus (DENV) diseases can be classified in to two categories, uncomplicated dengue or dengue fever (DF), dengue hemorrhagic fever (DHF), or dengue with warning signs [1].A hallmark of DHF is increased vascular permeability, resulting in plasma leakage, rash, bleeding, and cardiovascular collapse.The leading cause of death and morbidity in DHF is vascular leakage and its secondary complications [7].DENV consists of four genetically associated serotypes (DENV-1 to DENV-4), which can be further subdivided into genotypes [8].Infection with one type of dengue virus confers life-long immunity; subsequent infections with other types of the virus do not confer immunity [9].The shift of a predominantly circulating dengue serotype and/or genotype is responsible for increasing the incidence and severity of dengue outbreaks worldwide [10].Dengue serotyping/genotyping shifting facilitate infectivity via a non-neutralizing cross-reactive antibody that enhance disease by binds to the second infecting DENV virus serotype result of accelerate virus entry into cells through an Fc receptor-mediated endocytosis and by suppressing intracellular innate responses against the virus [11,12].Similarly, DENV vaccine recipients who have previously been exposed to DENV may be at risk for serious illness because of the increased virus production and suppressed antiviral defenses [13].Live attenuated vaccines elicit both the innate immune system and provide antigens to stimulate adaptive immunity [14].DENV infection triggers an innate immune response, which alters gene expression profiles and leads to host-pathogen interactions [15].This interaction is mediated by several pattern-recognition receptors (PRRs), including the endosome-associated Toll-like receptor-3 (TLR-3) and TLR-7 and the cytoplasmic retinoic acid-inducible gene I (RIG-I), NOD-like receptor protein 3 (NLRP-3)-specific inflammasome and melanoma differentiation-associated protein 5 (MDA-5), are responsible for activating downstream signaling pathways via recognition of their ligands [16].Activating downstream signaling pathways can result in the production of type I and type III interferons (IFNs) and proinflammatory cytokines, and hundreds of interferon-stimulated genes (ISGs) that help prevent DENV infection [17,18].IFN production is intrinsically linked to the transcriptional and posttranslational modification and regulate gene expression in the signaling pathways that stimulate DENV-specific adaptive immunity [19].A number of classical innate immunity pathways, such as RNAi, TLR, and JAK/STAT pathways, can inhibit viral infection and immunity in mosquitoes [20].Dengue virus inhibits IFN-α and IFN-β signaling by suppressing the JAK-STAT pathway, resulting in reduced host immunity against the virus [21].However, DENV infections alter gene expression profiles in dengue infected patients.
The combination of several factors such as reduced host defense, increased virus uptake and delayed viral clearance show synergism to produce higher virus titer that causes severe outcome.Early identification is more challenging due to the diverse clinical manifestation of dengue infection, as well as its substantial similarities to other febrile viral infections [7].To control the dengue diseases, it is necessary to identify the dengue serotype, genotype, diseases process, pathogenesis, and host response to dengue infection.The genomic characteristic of DENV outbreak and host response to dengue infection in Bangladesh is largely unknown.So we have therefore chosen to study the genomic characteristic and host response of DENV infection.pathways were found in dengue patients.These TLR pathways were particularly involved for the activation of innate system coupled with adaptive immune system that probably involved the rapid elimination of dengue virus infected cells.These differentially expressed genes could be further investigated for target based prophylactic interventions for dengue.Conclusion This is a first report describing DENV complete genomic features and differentially expressed genes in patients in Bangladesh.These genes may have diagnostic and therapeutic values for dengue infection.1).By aligning a reference genome (GenBank accession no.NC_001475.2),the consensus was obtained using the Burrows-Wheeler Aligner (BWA v0.7.17) [24], SAMtools v1.12 [25], and BEDTools v2.30.0 [26].Snippy was used to assess indel and area-wise mutation coverage, and the genome was repaired accordingly.Freebays (https://github.com/freebayes/freebayes) is the variation caller of snippy, with the minimum number of reads covering a site to be considered (default = 10) and the minimum VCF variant call "quality" (default = 100) [27].GENOME DETECTIVE VIRUS TOOL (https://www.genomedetective.com/app/typingtool/virus/)was also implemented to virus detection and de-novo assembly of the virus [28].Genome Detective begins by categorizing short reads into groupings, or buckets.All readings from a single viral species are placed in the same bucket.DIAMOND [29] is then allocated to each bucket for taxonomic identification.Once all of the readings have been sorted into buckets, each bucket is constructed from scratch using metaSPAdes [30].BLASTx and BLASTn employ contigs to search for possible reference sequences against the National Center for Biotechnology Information (NCBI) RefSeq viral database.Advanced Genome Aligner (AGA) [31] is used to link the contigs for each particular species [31].A report is created that refers to the final contigs and consensus sequences, both of which are accessible as FASTA files [28].DENGUE VIRUS TYPING TOOL (https://www.genomedetective.com/app/typingtool/dengue/) implemented in Genome detective [28] were used to identify virus genotype.d.Phylogenetic analysis: NCBI BLASTn was used to find the most closely related strain, and representative sequences were retrieved for phylogeny.Outgroup DENV-3 genotypes II, III, and V were also collected, and all sequences were genotyped using the DENGUE VIRUS TYPING TOOL (https://www.genomedetective.com/app/typingtool/dengue/).Multiple alignment using fast Fourier transform (MAFFT) online version (v7.0) (https://mafft.cbrc.jp/alignment/server/)were used for multiple sequence alignment with their default parameters [32].The complete viral genome was built using a neighbor joining (NJ) phylogenetic tree [33] with a bootstrap value of 1000 and Jukes-Cantor substitution model [34].Phylo.io(http://phylo.io/) was used to visualize the phylogenetic tree.Another env gene sequence based phylogenetic tree was reconstructed that has more DENV3 isolates reported from Bangladesh [5] which contains 135 isolates.After making ClustalW alignment, a model selection was carried out based on the Akaike information criterion in MEGA11 [35].A neighbor-joining tree was reconstructed using the Tamura-Nei model with Gamma distributed rates among sites.1000 bootstrap replications were used to test the phylogeny.
e. Phylodynamic analysis: Bayesian inference (BI) analyses were conducted by using Timra-Nei with a gamma-distributed rate variation substitution model (TN93 + G).The Markov Chain Monte Carlo (MCMC) algorithm in the BEAST package v.1.10.4 was used [36].The calibration point was the date of isolation of each sample.Runs were performed using the Bayesian skyline as the tree prior under the uncorrelated relaxed molecular clock [37].The evolutionary analysis was run for 100 million steps and the trees were sampled every 1,000 states.Convergence of the MCMC chains were inspected using TRACER v.1.7.2 (http://tree.bio.ed.ac.uk) [38].Posterior trees were summarized by discarding the first 10% of the sampled trees and choosing the Maximum Clade Credibility (MCC) were summarized using TreeAnnotator v.1.10.4.The final tree was then visualized and plotted using FigTree v.1.4.4 (http://tree.bio.ed.ac.uk).A posterior mean for each evolutionary parameter was presented along with a 95% Bayesian credibility interval.f.Recombination analysis: The DENV-3 recombination was detected using seven approaches, including RDP [39], Chimaera [40], BootScan [41], 3Seq [42], NCONV [43], MaxChi [44], and SiScan [45], all of which are accessible in the Recombination Detection Program (RDP version 5) (http://web.cbio.uct.ac.za/~darren/rdp.html)[46], with p < 0.01.If a sequence was discovered by at least three ways using the multiple comparison correction setting option, the recombinant event was identified [47].g.Human gene expression analysis: RNA-Seq Alignment tools from illumine basespase (https:// www.illumina.com/products/by-type/informaticsproducts/basespace-sequence-hub/apps/rna-seqalignment.html) were applied for human gene expression.Here, sample reads were mapped by using the quick universal RNA-seq aligner, STAR [48].The aligner utilized the NCBI GRCh38 human genome as a reference genome.Salmon, a rapid and bias-aware transcript expression quantification tool [49] was used to quantification of expressed genes from transcriptome samples.h.Statistical analysis: DESeq2 (release 3.13) in R (version 4.0) was used for differential gene expression and statistical analysis.DESeq2 [50] was chosen since it is a popular parametric tool with a detailed and constantly updated user guideline.(http:// bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html).Because DESeq2 inherently corrects for library size, it is critical to supply un-normalized raw read counts as input.To adjust for discrepancies in sequencing depth, we utilized a variance stabilizing transformation [51].The Benjamini-Hochberg technique [52] was used to modify P-values for multiple testing.For the selection of DE genes, a false discovery rate (FDR) adjusted p-value < 0.05 was used.Ggplot2 R package was used to visualize the data.EnhancedVolcano R package was also used for the volcano plot.The data was shown using the R tool ggplot2 [53]

Results
Dengue is a one of the most common mosquito-borne viral diseases that affect millions of people each year globally.WHO state that Asia is responsible for 75% of all dengue cases worldwide, especially in countries such as the Philippines, Indonesia, and Thailand.The Dengue seroprevalence in Bangladesh is lower than that of the Southeast Asian countries, but the trend is changing rapidly.Bangladesh experienced its largest dengue outbreak on record during the year 2019 (Lancet).However, there is a lack of genomic surveillance of Dengue in the country.As a result, in this study, we compared RNA-seq data from dengue patients to examine the whole genome sequence, genotype, and genomic characterization of the dengue virus, as well as differential gene expression, overrepresented pathways, and upstream regulators.This study includes a total of 21 dengue patients and 3 healthy controls.Both male (n = 13) and female (n = 8) patients with ages ranging from 3 to 61 years were enrolled in the study (Table 1).The average days of symptoms were 4.5.At the time of patient recruitment, a blood sample was taken.

Metagenomic sequencing and annotation of the dengue virus genomes
The complete viral genomes of 21 samples positive for Dengue virus were assembled by metagenomic sequencing.The average coverage of the assembled genomes ranged from 28.15 × to 82474.8 × and the genome assembly length varied between 10,647 bp (sample: BCSIR-IBS-D19) and 10,707 bp (sample: BCSIR-05) (Table 1).All the positive samples (n = 21) belonged to DENV-3 Genotype I (Table 1; Fig. 1).No reads from the healthy control (n = 3) aligned with the Dengue virus (Table 1).

Phylogenetics and phylodynamics of dengue virus in Bangladesh
After analysing these sequences, we investigated their phylogenetic relationships to existing dengue virus genomes based on both the complete genomes and the envelope protein E. In both phylogenetic trees, the genomes from this study belonged to DENV-3 Genotype I and the genomes of these samples were highly similar and belonged to a cluster of samples of DENV-3 genotype I originated from China and Thailand that also includes some sequences from Bangladesh collected between 2009 and 2017 (Fig. 2.A).In both the trees, isolates from Genotype II, III, and I were clustered separately.The env gene tree has more DENV-3 isolates from the previously reported Dengue epidemic in Bangladesh (Supplementary Fig. S1).It shows that 2021 epidemic samples are closely related to 2017 Dengue epidemic in Bangladesh.
Bayesian inference analysis comparing all the sequences produced a phylogenetic tree that has similar topography to the previously described NJ tree.The Bangladeshi DENV-3 sequences are distributed in two clades (Fig. 1), characterized as genotype I for DENV3 epidemic of 2017 and 2021, and the previous epidemics in genotype II that contains isolates from 2009, 2008, 2006, 2002, and 2000.This confirms a previously reported clad shift in DENV3 epidemic in Bangladesh [5].However, the tree also indicates that the 2021 epidemic isolates, which have strong similarities among themselves, are well separated from 2017 epidemic isolates in Bangladesh but interleaved by two DENV3 isolates sampled in 2019.These two isolates, 19XN13542 and 19XN14065, were sampled from travelers return from Bangladesh based on the NCBI GenBank description.This gives us an observational time point between 2017 and 2021 and indicates that the DENV3 was evolving in Bangladesh in-between times.The estimated median tMRCA from genotype I DENV3 sequences from Bangladesh (epidemic of 2017 and 2021) was 2013 (95% HPD 2010-2014).This was right after the 2006-2009 DENV3 occurrence from genotype II reported in the previous study, where previously reported clad-based genotype shift presumably happened.

Principal component analysis (PCA) and hierarchical clustering analysis
We performed PCA to investigate the clustering of the samples in each group (patients and healthy); whether they clustered within the same group or with other groups.First, we used HTSeq [56] to count reads that uniquely aligned to one gene, and PCA plots were then generated by importing the data into DESeq2 [50].The PCA results demonstrated that most of the positive samples clustered together distinctly separated from the samples of the control group (Fig. 3).However, the control group did not cluster so closely, likely due to lower sample number (n = 3) in the control group.Additionally, PCA scree plots showed that principal component 1 (PC1) and 2 (PC2) accounted for 36% and 12%, respectively, of the analysed data (Fig. 3).

Differentially expressed genes (DEGs)
A total of 19,267 genes were identified by at least one read from all the samples including dengue positive and negative (control) samples.Positive samples aligned to a total of 17,375 genes whereas a total of 12,314 genes aligned in the negative control samples.Among them, 6,953 and 1,892 unique genes found in dengue-positive and healthy control samples respectively (Fig. 4, Supplementary Data 1).Overall, 2,686 differentially expressed genes (DEGs), with a q-value < 0.05, were detected in the analysis of DESeq2 (Supplementary Data 1).

Gene Ontology (GO) and pathway analysis of DEGs
To characterize the GO terms, including biological processes, cellular components, molecular functions, and functional pathways of DEGs, we performed over-representation tests in PANTHER version 11.1 (Fig. 8 and Supplementary Data 2).We used the GO-Slim PAN-THER annotation data set, which represents phylogenetically inferred annotations [57].

Discussion
Dengue virus (DENV) is one of the leading threats to public health causing nearly 390 million infections every year worldwide [58], with an estimated 2.5 billion people at risk of dengue infection.This mosquito-borne viral disease affects more than 100 endemic countries globally with most disease burden limited to tropical and subtropical regions [59].The recent upsurge in the global occurrence of dengue epidemics has prompted research into dengue virus biology and pathology to develop therapeutics for infection prevention and control.Despite Fig. 7 The top 24 DEGs based on p-value adjusted (padj 2.48E-06) were displayed.The Y-axis shows the normalized gene counts, while the X-axis shows the groupings (negative and positive).Details of the analysis can be found in Supplementary Data 1 numerous efforts to unravel the DENV pathogenesis, the underlying mechanisms of this disease remain elusive and must be investigated further.In this study, we depict the first comprehensive picture of complete genomic characteristics as well as differentially expressed genes in dengue patients in Bangladesh.The expression levels may be implicated to predict the disease progression and, importantly, to design prophylactic therapeutics for the disease.
Dengue infections occur with four distinct serotypes namely DENV1, DENV2, DENV3 and DENV4, posing a major health threat globally.Although mostly the primary infection causes mild to asymptomatic infection with activation of immune responses against DENV serotypes, the severity of the disease is enhanced via heterotypic infection by various serotypes due to antibodydependent enhancement (ADE) [60].Dengue poses a significant health and economic burden in the South-East Asia region.Dengue virus has now become hyper endemic in this region for multiple serotypes and genotypes [61].Co-circulation of all four DENV serotypes is causing the coinfection of these countries with multiple serotypes due to shifts in the dominant serotype [61].Notably, shift of dominant serotypes from DENV-2 in 2013 to DENV-3 in 2016 occurred in India.Moreover, in 2017, studies conducted in the country revealed that DENV-2 and DENV-3 were dominant in the southern part of India while DENV-1 and DENV-2 were prevalent in the northern part of India [62][63][64].Recently, the shifting of predominant serotypes also observed in Thailand (from DENV-2 and DENV-4 to DENV-3), [65] Sri Lanka (from DENV-1 to DENV-2) [66], Myanmar (from predominant serotype 1 to co-dominant serotypes 1, 2, and 4) [67], Indonesia, Bhutan (from DENV-3 to DENV-1), and in Nepal [68][69][70] .
In Bangladesh, DENV-3 of genotype II was predominant till around 2000, which is closely related to the strains from Thailand and Myanmar [71,72].From 2013 to 2016, DENV-2 was predominant, with the coexistence of DENV-1.Neither DENV-3 nor DENV-4 were detected during this period [4,73].However, genomic characterization of the outbreak in 2019 having 112,000 cases and 129 deaths revealed that DENV type 3 genotype I was predominant, accounting for 93% of tested samples followed by DENV-3 genotype III were detected in two samples from separate districts, and only one DENV-2 cosmopolitan genotype in the Dhaka, capital city of Bangladesh [74].In this study we extensively analyze the Fig. 8 Pathway enrichment and GO phrase analysis for negative and positive groups found in PANTHER Functional classification viewed in gene list analyses.The X-axis of the stack plot shows the GO term, and the Y-axis depicts the proportions of genes for each GO term in each group.The positive group is represented by orange, while the negative group is represented by green.Details of the analysis can be found in Supplementary Data 2 transcriptome as well as alternatively spliced transcriptome of blood samples collected from participants with and without dengue infection.We also reveal the whole genome sequence (WGS) of dengue viruses from our samples to determine serotype, genotype, phylogenetic relation and phylodynamics.Similarly, the present study revealed the predominance of DENV-3 genotype I in all the dengue positive samples (n = 21) from the outbreak of 2021.Interestingly, the viruses closely resemble the strains from China and Thailand.These findings indicate the shift of predominant serotypes from DENV-2 to DENV-3 in the recent outbreaks.The viruses were likely introduced from China or Thailand as previous reports of introduction of the DENV-2 Cosmopolitan genotype from India [5].These all suggest the continuous crossborder movement of the virus into the South-East Asian region and a requirement for coordinated monitoring and surveillance of DENV [61,75].Additionally, these epidemic cycles are potentially responsible for shifting in predominant serotypes and genotypes and have often been associated with higher dengue incidence and severe cases in the territory.
There is a lack of a comprehensive spectrum of gene expression profiles that may contribute to the setting up of dengue severity globally.Several studies have reported gene expression and molecular responses in dengue patients [76][77][78][79].Most of these data produced through the microarray technology lack of a comprehensive spectrum of gene expression patterns due to limited accuracy of expression measurements, especially for transcripts with low abundance [80].However, recently, the RNA-Seq technology presents a useful tool in disease study with more accuracy at detecting low abundance transcripts and identifying genetic variants [81].RNA-Seq also has a wider dynamic range covering the identification of more differentially expressed genes with a higher fold-change in transcript levels of a host [80].Thus, in the present study, using the RNA-Seq based platform, we provide the first comprehensive overview of the differentially expressed genes in dengue patients compared to healthy controls in Bangladesh.Based on the transcriptomic signatures elicited by dengue patients, we demonstrate differential gene expression patterns.Gene ontology (GO) analysis depicts the function of the differentially expressed gene.
Our results demonstrated that more genes are expressed in dengue positive samples (n = 17,375) compared to the healthy controls (n = 12,314) with a count of 6,953 and 1,892 unique genes, respectively.These indicate the higher physiological activities stimulated by a viral infection in the patients [77,78].Identification of 2,686 differentially expressed genes (DEGs), with a q-value < 0.05, refers significant alterations in gene expression in dengue patients.Moreover, the heatmap shows distinct clusters of gene expression patterns between patients and healthy individuals (Fig. 5).Though the positive samples were subclustered into two groups, no correlation were observed between age, sex and gene expression pattern, thus indicating a similar host response independent of the factors [77].Further, the plot counts function in DESeq2 identified the top 24 genes with the smallest q-values which include KCNQ1OT1, GRIN2B, TMSB4X, TTN, F13A1, TSIX, PPBP, WDR87, IGFN1, B2M, YWHAZ, MBNL1, TUBB1, LINC01206, NCOA4, LOC440300, LUZP6, MTPN, OST4, RSU1, C22orf46, MEG3, NAP1L1 and ZNF793.Among them, compared to the negative samples, the genes KCNQ1OT1, GRIN2B, TTN, TSIX, WDR87, IGFN1, LINC01206, LOC440300, C22orf46, MEG3, and ZNF793 were upregulated in dengue-positive samples, whereas the others downregulated.These genes can be targeted for diagnostic and therapeutic interventions for dengue.
Within the biological process, the most characterized GO category, significant upregulation (p = < 0.001) of multicellular organismal process, nervous system process, sensory perception of the chemical stimulus, and G protein-coupled receptor signalling pathway in the dengue patients correlates with the stimulation of various physiological pathways in the patients [82].However, within the cellular component, cellular anatomical entity, protein-containing complex and intracellular significantly decreased (p = < 0.001) in dengue patients.Cellular metabolic process, cellular process, biological process and cellular protein metabolic process family were significantly downregulated in all the positive samples, which was contrasting with the previous report [79].Importantly, overall, defense/immunity protein class increased in patients compared to the healthy controls, with increased GO of immune system process.Cytokine activity significantly (p = < 0.001) increased (fold change 2.250) in the patients.These indicate post-infection immune stimulation in dengue patients, which corroborated with previous findings [76,78,79].Strikingly, number of genes decreased related to B cell activation (patients: 54, healthy: 58) and T cell activation (patients: 72, healthy: 79).To investigate these phenomena, comprehensive studies are warranted aiming at different time points of the infection.Moreover, upregulation of blood coagulation, toll receptor signalling pathway, and inflammation mediated by chemokine and cytokine signalling pathway represent the innate response to the infection [78].Transforming growth factor-beta 1 (TGF-β1) was highly activated upstream regulator in the positive samples could be associated with apoptosis of the infected cells.These differentially expressed genes could be further investigated for target based prophylactic interventions for dengue.
Collectively, these studies demonstrate specific transcriptional pathways and immune cell types implicated in inducing IFNs and ISGs, emphasizing the importance of gene regulation mechanisms at the splicing and isoform levels for eliciting an immunological response.
Our study shows for the first time that a comprehensive view of transcriptome expression in plasma of dengue patients in Bangladesh.It offers a valuable resource for understanding transcriptome that might be involved in the manifestation and progression of dengue.It might help to discover valuable gene for diagnosis, therapeutic and key transcriptional regulatory factors that will block viral replication and transmission in mosquitoes.We also clearly showed that clade shifting which might be a predictor of dengue severity.

Conclusion
Dengue remained a major public health threat due to consistent increase in the cases and deaths of infections in Bangladesh.DENV type 3 genotype I is associated with the recent dengue outbreaks in the country.The circulating virus strains closely resemble with the strains from China and Thailand, indicating possible introduction of the DENV-3 genotype I from these countries.Moreover, our results highlight the differentially expressed signature genes associated with dengue, which could be implicated for diagnostic and prophylactic interventions for dengue.Our study has a limitation of small sample size, and this limited our statistical power and our ability to refine validate, and formally compare genetic classifiers.However, we depict the first comprehensive picture of complete genomic characteristics as well as differentially expressed genes in dengue patients in Bangladesh.Further investigation and continuous genomic surveillance are warranted to explore the shift in predominant genotypes and viral pathogenesis.

Fig. 2 Fig. 4
Fig.2Neighbor joining (NJ) phylogenetic tree of dengue whole genome and envelope protein E. (A) A neighbor joining (NJ) phylogenetic tree with a bootstrap value of 1000 and the Jukes-Cantor substitution model were used to construct the whole viral genome and (B) the envelope protein E phylogenetic tree of 95 DENV-3 strains.The DENV-3 genotype I strain sequenced in this study is shown in red.To root the tree, the full sequences of DENV-3 genotypes II, III, and V were employed

Fig. 3
Fig. 3 Dengue positive and negative groups were subjected to principal component analysis (PCA).The variance stabilizing transformation in DESeq2 was used to identify principal component 1 (PC1) and principal component 2 (PC2).The percentage of variance represents the amount of variance explained by PC1 and PC2 Continual genomic surveillance is required to further investigate the shift in dominant genotypes in relation to viral pathogenesis.

Table 1
Sample and sequence information of the collected samples in this study