Molecular Epidemiology of Wild Poliovirus Type 1 and Shift in the Historically Reservoir Areas of Pakistan During 2019-2022.

BACKGROUND: Pakistan is one of the two countries endemic for wild poliovirus type 1 (WPV1). Active clinical and environmental wastewater surveillance along with laboratory investigation is an integral and primary component of the polio eradication strategies. The current study is mainly focused on the virological data to understand the current epidemiology of WPV1 in Pakistan during 2019-2022. METHODS: 141,037 stool specimens of patients reported with acute flaccid paralysis (AFP) and 3,171 wastewater samples were tested for poliovirus detection using cell culture and PCR. Phylogenetic analysis of WPV1 was performed using MEGA and Nextstrain. RESULTS: Poliovirus isolates were classified into 15 distinct genetic clusters with multiple transmission lineages. Spatio-temporal trends indicated a significant decline in the incidence of poliomyelitis reported in 58 districts in 2019 to just 3 in 2022. The historical reservoirs in Peshawar, Quetta, and Karachi successfully eliminated the indigenous transmission chains of wild poliovirus active there for years. CONCLUSIONS: Our findings reinforce the evolving epidemiology of poliovirus in Pakistan which is now confined to South KP. All historically known reservoirs in Peshawar, Karachi and Quetta blocs are now free of poliovirus. Intensified clinical and environmental surveillance should be maintained to eliminate the very few remaining transmission lineages and certify the poliovirus eradication by 2026.

Unravelling the complex interplay of age, comorbidities, and multimorbidities in COVID-19 disease progression: Clinical implications and future perspectives.

Introduction: The COVID-19 infection as an inflammatory disease has posed significant challenges to global public health due to multi-factor risks associated with it leading to disease severity and mortality. Understanding the effect of age and comorbidities on overall disease progression is crucial to identify highly susceptible individuals and to develop effective disease management strategies in a resource limited country like Pakistan. Methodology: A retrospective study was conducted on hospitalized COVID-19 patients to assess the prevalence of various comorbidities among different age groups and their effect on disease severity and mortality rate. Results: In this retrospective study, a cohort of 618 hospitalized COVID-19 patients was analyzed, consisting of 387 males (62.6 %) and 231 females (37.4 %). Notably, the young age group (15-24 years), had the lowest frequency of hospitalized COVID-19 patients, while no case was observed in children (≤14 years) showing a significant association (p < 0.001) of age and disease prevalence. Comorbidities were observed in 63.9 % of COVID-19 patients including hypertension (HTN), diabetes mellitus (DM), ischemic heart diseases (IHD), asthma, chronic kidney disease (CKD) and tuberculosis (TB). The most common comorbidities were HTN (42.1 %) followed by DM (33.8 %), IHD (16.5 %), asthma (11.2 %), CKD (7.9 %) and TB (1.9 %).Furthermore, the study revealed a significant association between comorbidities, age groups, and the need for non-invasive ventilation (NIV) (p < 0.001), mechanical ventilation (MV) (p < 0.001), and intensive care unit (ICU) admission (p < 0.001). Patients with specific comorbidities and those in the older age group (≥65 years) demonstrated a higher need for these interventions. However, patients without any comorbidity consistently exhibited the highest cumulative proportion of survival at each time point, indicating better overall survival outcomes. In contrast, patients with multimorbidities of DM/HTN/IHD, HTN/IHD, and DM/HTN/CKD had comparatively lower survival rates and higher mortality rates (p < 0.001). Conclusion: This research highlights the significant impact of age, comorbidities and multimorbidities on the severity and mortality of COVID-19 patients. It highlights the importance of considering these factors in tailoring effective management strategies for patients with COVID-19 or other infectious respiratory diseases.

Serotype and genomic diversity of dengue virus during the 2023 outbreak in Pakistan reveals the circulation of genotype III of DENV-1 and cosmopolitan genotype of DENV-2.

Dengue, a mosquito-borne viral disease, poses a significant public health challenge in Pakistan, with a significant outbreak in 2023, prompting our investigation into the serotype and genomic diversity of the dengue virus (DENV). NS-1 positive blood samples from 153 patients were referred to the National Institute of Health, Pakistan, between July and October 2023. Among these, 98 (64.1%) tested positive using multiplex real-time PCR, with higher prevalence among males (65.8%) and individuals aged 31-40. Serotyping revealed DENV-1 as the predominant serotype (84.7%), followed by DENV-2 (15.3%). Whole-genome sequencing of 18 samples (DENV-1 = 17, DENV-2 = 01) showed that DENV-1 (genotype III) samples were closely related (>99%) to Pakistan outbreak samples (2022), and approx. > 98% with USA (2022), Singapore and China (2016), Bangladesh (2017), and Pakistan (2019). The DENV-2 sequence (cosmopolitan genotype; clade IVA) shared genetic similarity with Pakistan outbreak sequences (2022), approx. > 99% with China and Singapore (2018-2019) and showed divergence from Pakistan sequences (2008-2013). No coinfection with dengue serotypes or other viruses were observed. Comparisons with previous DENV-1 sequences highlighted genetic variations affecting viral replication efficiency (NS2B:K55R) and infectivity (E:M272T). These findings contribute to dengue epidemiology understanding and underscore the importance of ongoing genomic surveillance for future outbreak responses in Pakistan.

Molecular characterization of rotavirus indicates predominance of G9P[4] genotype among children with acute gastroenteritis: First report after vaccine introduction in Pakistan.

Globally, Group A rotavirus (RVA) is the leading cause of acute gastroenteritis in children under 5 years old, with Pakistan having the highest rates of RVA-related morbidity and mortality. The current study aims to determine the genetic diversity of rotavirus and evaluate the impact of Rotarix-vaccine introduction on disease epidemiology in Pakistan. A total of 4749 children, hospitalized with acute gastroenteritis between 2018 and 2020, were tested at four hospitals in Lahore and Karachi. Of the total, 19.3% (918/4749) cases were tested positive for RVA antigen, with the positivity rate varying annually (2018 = 22.7%, 2019 = 14.4%, 2020 = 20.9%). Among RVA-positive children, 66.3% were under 1 year of age. Genotyping of 662 enzyme-linked immuno sorbent assay-positive samples revealed the predominant genotype as G9P[4] (21.4%), followed by G1P[8] (18.9%), G3P[8] (11.4%), G12P[6] (8.7%), G2P[4] (5.7%), G2P[6] (4.8%), and 10.8% had mixed genotypes. Among vaccinated children, genotypes G9P[4] and G12P[6] were more frequently detected, whereas a decline in G2P[4] was observed. Phylogenetic analysis confirmed the continued circulation of indigenous genotypes detected earlier in the country except G9 and P[6] strains. Our findings highlight the predominance of G9P[4] genotype after the vaccine introduction thus emphasizing continual surveillance to monitor the disease burden, viral diversity, and their impact on control of rotavirus gastroenteritis in children.

Genomic characterization of human respiratory syncytial virus circulating in Islamabad, Pakistan, during an outbreak in 2022-2023.

In this study, conducted at the National Institute of Health, Islamabad, during an outbreak of human respiratory syncytial virus (hRSV) from December 2022 to January 2023, the first whole-genome sequences of hRSV isolates from Islamabad, Pakistan, were determined. Out of 10 positive samples, five were sequenced, revealing the presence of two genotypes: RSV-A (GA2.3.5, ON1 strain) and RSV-B (GB5.0.5.a, BA-10 strain). A rare non-synonymous substitution (E232G) in G the protein and N276S in the F protein were found in RSV-A. In RSV-B, the unique mutations K191R, Q209R, and I206M were found in the F protein. These mutations could potentially influence vaccine efficacy and viral pathogenicity. This research underscores the importance of genomic surveillance for understanding RSV diversity and guiding public health responses in Pakistan.

Insights into measles virus: Serological surveillance and molecular characterization.

BACKGROUND: Measles has been a significant public health concern in Pakistan, especially in the Khyber Pakhtunkhwa (KPK) province, where sporadic and silent epidemics continue to challenge existing control measures. This study aimed to estimate the prevalence and investigate the molecular epidemiology of the measles virus (MeV) in KPK and explore the vaccination status among the suspected individuals. METHODS: A cross-sectional study was conducted between February and October 2021. A total of 336 suspected measles cases from the study population were analyzed for IgM antibodies using Enzyme-Linked Immunosorbent Assay (ELISA). Throat swabs were randomly collected from a subset of positive cases for molecular analysis. Phylogenetic analysis of MeV isolates was performed using the neighbor-joining method. The vaccination status of individuals was also recorded. RESULTS: Among the suspected participants, 61.0% (205/336) were ELISA positive for IgM antibodies, with a higher prevalence in males (64.17%) compared to females (57.04%). The majority of cases (36.0%) were observed in infants and toddlers, consistent with previous reports. The majority of IgM-positive cases (71.7%) had not received any dose of measles vaccine, highlighting gaps in vaccine coverage and the need for improved immunization programs. Genetic analysis revealed that all MeV isolates belonged to the B3 genotype, with minor genetic variations from previously reported variants in the region. CONCLUSION: This study provides valuable insights into the genetic epidemiology of the MeV in KPK, Pakistan. The high incidence of measles infection among unvaccinated individuals highlights the urgency of raising awareness about vaccine importance and strengthening routine immunization programs.

Crimean-Congo Hemorrhagic Fever Virus Diversity and Reassortment, Pakistan, 2017-2020.

Sporadic cases and outbreaks of Crimean-Congo hemorrhagic fever (CCHF) have been documented across Pakistan since 1976; however, data regarding the diversity of CCHF virus (CCHFV) in Pakistan is sparse. We whole-genome sequenced 36 CCHFV samples collected from persons infected in Pakistan during 2017-2020. Most CCHF cases were from Rawalpindi (n = 10), followed by Peshawar (n = 7) and Islamabad (n = 4). Phylogenetic analysis revealed the Asia-1 genotype was dominant, but 4 reassorted strains were identified. Strains with reassorted medium gene segments clustered with Asia-2 (n = 2) and Africa-2 (n = 1) genotypes; small segment reassortments clustered with the Asia-2 genotype (n = 2). Reassorted viruses showed close identity with isolates from India, Iran, and Tajikistan, suggesting potential crossborder movement of CCHFV. Improved and continuous human, tick, and animal surveillance is needed to define the diversity of circulating CCHFV strains in Pakistan and prevent transmission.

Integrated approach for detection of SARS-CoV-2 and its variant by utilizing LAMP and ARMS-PCR.

Global impact of COVID-19 pandemic has heightened the urgency for efficient virus detection and identification of variants such as the Q57H mutation. Early and efficient detection of SARS-CoV-2 among densely populated developing countries is paramount objective. Although RT-PCR assays offer accuracy, however, dependence on expansive kits and availability of allied health resources pose an immense challenge for developing countries. In the current study, RT-LAMP based detection of SARS-Cov-2 with subsequent confirmation of Q57H variant through ARMS-PCR was performed. Among the 212 collected samples, 134 yielded positive results, while 78 tested negative using RT-LAMP. Oropharyngeal swabs of suspected individuals were collected and processed for viral RNA isolation. Isolated viral RNA was processed further by using either commercially available WarmStart Master Mix or our in house developed LAMP master mix separately. Subsequently, the end results of each specimen were evaluated by colorimetry. For LAMP assays, primers targeting three genes (ORF1ab, N and S) were designed using PrimerExplorer software. Interestingly, pooling of these three genes in single reaction tube increased sensitivity (95.5%) and specificity (93.5%) of LAMP assay. SARS-CoV-2 positive specimens were screened further for Q57H mutation using ARMS-PCR. Based on amplicon size variation, later confirmed by sequencing, our data showed 18.5% samples positive for Q57H mutation. Hence, these findings strongly advocate use of RT-LAMP-based assay for SARS-CoV-2 screening within suspected general population. Furthermore, ARMS-PCR also provides an efficient mean to detect prevalent mutations against SARS-Cov-2.

Evolutionary analysis of seasonal influenza A viruses in Pakistan 2020-2023.

INTRODUCTION: Influenza A viruses cause global health concerns due to their high amino acid substitution rates. They are linked to yearly seasonal epidemics and occasional pandemics. This study focused on sequencing influenza A virus strains in Pakistan. MATERIALS AND METHODS: We analyzed the genetic characteristics of influenza A(H1N1)pdm09 and A(H3N2) viruses circulating in Pakistan from January 2020 to January 2023. Whole genome sequences from influenza A (n = 126) virus isolates were amplified and sequenced by the Oxford Nanopore (MinION) platform. RESULTS: The HA genes of influenza A(H1N1)pdm09 underwent amino acid substitutions at positions K54Q, A186T, Q189E, E224A, R259K, and K308R in sequenced samples. The HA genes of influenza A(H3N2) had amino acid substitutions at G53D, E83K, D104G, I140M, S205F, A212T, and K276R in the sequenced samples. Furthermore, the HA gene sequences of influenza A(H1N1)pdm09 in this study belonged to subclade 6B.1A.5a.2a. Similarly, the HA gene sequences of influenza A(H3N2) were classified under six subclades (3C.3a.1 and 3C.2a1b.2a [2, 2a.1, 2b, 2c, and 2a.3b]). Notably, amino acid substitutions in other gene segments of influenza A(H1N1)pdm09 and A(H3N2) were also found. CONCLUSION: These findings indicate influenza A(H1N1)pdm09 and A(H3N2) viruses co-circulated during the 2020-2023 influenza season in Pakistan. Continued surveillance is crucial for real-time monitoring of possible high-virulence variation and their relevance to existing vaccine strains.

Genomic characterization of the first mpox and varicella-zoster co-infection in Pakistan (2023) through next-generation sequencing.

The global mpox outbreak spanning 2022-2023 has affected numerous countries worldwide. In this study, we present the first report on the detection, whole-genome sequence, and coinfection of the mpox virus and varicella zoster virus (VZV) from Pakistan. During April-May 2023, samples from 20 suspected cases of mpox were tested at the National Institutes of Health, Islamabad among which 4 tested positive. All four cases had a travel history of Saudi Arabia. All the suspected samples were processed by using a Zymo research kit for DNA extraction, followed by qRT-PCR amplification by using a DaAn Gene detection kit for the mpox virus. Further, two of the positive samples with a low Ct value (<20) were subjected to whole-genome sequencing using a metagenomic approach on the iSeq (Illumina) platform. The sequencing results revealed Clade IIb and genotype A.2.1 of MPXV, which clustered with viruses from Slovenia and the UK in July and June 2022, respectively. Our analysis identified two novel nonsynonymous substitutions in mpox virus, namely V98I in OPG046 and P600S in OPG109. Furthermore, we successfully retrieved the complete genome of VZV from the same sample, belonging to Clade 5. This study represents the first positive case of MPXV in Pakistan and the coinfection of mpox and VZV by using a metagenome approach providing insights into their complete genomes. Our results highlight the importance of surveillance at the point of entries, strengthening lab capacities including next-generation sequencing, and using differential diagnosis for timely and accurate detection of mpox cases.

Whole-genome sequencing of Crimean-Congo hemorrhagic fever virus circulating in Pakistan during 2022.

Pakistan is an endemic country for Crimean-Congo hemorrhagic fever (CCHF) and its Balochistan province is considered a hotspot region for circulation of the virus whereas sporadic cases have been reported from other parts of the country. Our study aims to investigate the genomic diversity of the CCHF virus circulating in Punjab and Khyber Pakhtunkhwa provinces of Pakistan. Between April to September 2022, 46 samples from suspected CCHF patients were tested, with 6 (13%) showing positive RT-PCR results. Among the positive cases, all were male, aged 14-48 years among which 4 were butchers. Three CCHF patients succumbed to the disease. The complete S-M-L-gene fragments of 4 positive samples were sequenced. The S and L segments belonged to the Asia-1 clade and clustered with regional strains from Iran, India, and Afghanistan. One M segment sequence grouped into Africa-2 along with those reported from India during 2016-2019. We report the detection of a reassorted virus (NIH-PAK-CCHF-03|2022) having Asia-1-Africa-2-Asia-1 (S-M-L) segment pattern for the first time from Pakistan. Mutational analysis showed M segments harboring eight mutations (T55A, S80P, T110I, T185A, T189A, A212T, and N239I/T) in the mucin-like domain, five mutations (D250N, T333S, I375V, M401I, A433T), four mutations (N545D, Y657F, K688R, and I824V) in GP38-domain, and three mutations (T1418N, A1431V, and G1449S) in Gc-domain. These findings highlight the significance of whole-genome sequencing of indigenous strains for a better understanding of the CCHFV evolution in Pakistan.

First report of coinfection and whole-genome sequencing of norovirus and sapovirus in an acute gastroenteritis patient from Pakistan.

Acute gastroenteritis is one of the most common diseases in infants and children in developing countries including Pakistan. In Pakistan, rotavirus (RVA) is known to contribute significantly to pediatric diarrheal illness, but the contribution of other viruses is still unclear. In the current study we have identified a case of mixed infection of norovirus (NoV) and sapovirus (SaV) in a 2-year-old child with acute gastroenteritis. The sample was initially processed for the detection of group A RVA through ELISA followed by NoV using RT-PCR assay. The sample tested positive for NoV RNA and was later subjected to whole-genome sequencing using meta-genome approach on Miseq (Illumina) platform. Sequencing results revealed GII.15 genotype of NoV that clustered with viruses from China and USA from 2017 to 2021. We also retrieved the complete genome of SaV (GI.1 genotype) from the same sample and phylogenetic analysis showed clustering with strains reported from Japan, South Korea, US, and Taiwan during 2012-2016. This is the first report from Pakistan that confirms coinfection of NoV and SaV and elucidates their whole genomes. We recommend initiation of NoV and SaV surveillance program to ascertain disease burden and explore genetic diversity, especially as RVA vaccines have been included in national immunization program.

Genomic Characterization of Dengue Virus Outbreak in 2022 from Pakistan.

Pakistan, a dengue-endemic country, has encountered several outbreaks during the past decade. The current study aimed to explore the serotype and genomic diversity of dengue virus responsible for the 2022 outbreak in Pakistan. From August to October 2022, NS-1 positive blood samples (n = 343) were collected from dengue patients, among which, (85%; n = 293) were positive based on RT-PCR. In terms of gender and age, dengue infection was more prevalent in male patients (63%; n = 184), with more adults (21-30 years; n = 94) being infected. The serotyping results revealed DENV-2 to be the most predominant serotype (62%; n = 183), followed by DENV-1 (37%; n = 109) and DENV-3 (0.32%; n = 1). Moreover, a total of 10 samples (DENV-2; n = 8, DENV-1; n = 2) were subjected to whole-genome sequencing. Among these, four were collected in early 2022, and six were collected between August and October 2022. Phylogenetic analysis of DENV-2 sequenced samples (n = 8) revealed a monophyletic clade of cosmopolitan genotype IVA, which is closely related to sequences from China and Singapore 2018, and DENV-1 samples (n = 2) show genotype III, which is closely related to Pakistan isolates from 2019. We also reported the first whole genome sequence of a coinfection case (DENV1-DENV2) in Pakistan detected through a meta-genome approach. Thus, dengue virus dynamics reported in the current study warrant large-scale genomic surveillance to better respond to future outbreaks.

Breakthrough cases of Omicron and Delta variants of SARS-CoV-2 during the fifth wave in Pakistan.

COVID-19 pandemic has severely affected Pakistan with 1,557,134 cases as of August 4, 2022. However, the data regarding breakthrough infections in Pakistan is scant. Hence, the objective was to analyze SARS-CoV-2 breakthrough infections with respect to vaccines and variants during the fifth wave in Pakistan. Therefore, the Department of Virology (NIH, Pakistan) genotyped 2,467 randomly selected individuals between November 2021 and February 2022 using the SNPsig® SARS-CoV-2 (EscapePLEX) kit (PrimerDesign, UK). P681R and K417N mutations were used to distinguish delta and omicron. Data on the patient's age, gender, date of collection, variant, and vaccination status were analyzed using Statistical Package for Social Sciences (SPSS) software. Among 2,467 genotyped samples, Omicron was detected in 58.6% (n = 1445), Delta in 40.4% (n = 998) and undetermined/wildtype variant in 24 samples. The vaccination status of omicron-positive patients showed (49.7%; n = 718/1445) and Delta-positive patients (39.67%; n = 396/998) to be fully vaccinated. Of note, a high percentage 85% of breakthrough cases (n = 947) were identified among fully vaccinated individuals (n = 1114). Among them, 85.9% (n = 617/718) belonged to omicron and 83.3% (n = 330/396) to delta. Moreover, 76.7% (n = 855) of vaccinated individuals (n = 1114) received Sinopharm (n = 432) and Sinovac (n = 423) vaccines. The majority of breakthrough subjects who contracted Omicron were vaccinated with Sinopharm (93.0%, n = 256) and delta with Cansino (100%, n = 44). Individuals vaccinated with Sinovac showed the most frequent breakthrough cases for both Omicron and Delta variant between the 4th and 6th months (n = 278) after primary vaccination as compared to the 7th to 9th months (n = 24) category. While in case of Sinopharm, maximum breakthrough cases occurred between 7th to 9th months (n = 234) as compared to the 4th to 6th months (n = 120) after primary vaccination. Omicron and Delta breakthrough cases in men (n = 364 and 193) are more frequently seen than women (n = 253 and 138) respectively and breakthrough majority cases (n = 392) occurred in individuals aged 18-33 years. Breakthrough cases limiting monitoring in Pakistan impose a substantial constraint on policymakers' ability to take timely effective decisions. Since the current study consists of only a 2,467-genotyped sample, comprehensive data should be analyzed.