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.

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.

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.

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.

Dyakonov waves generation at uniaxial chiral-plasma interface

Numerical analysis of Dyakonov waves generation has been carried out at uniaxial chiral-plasma planar interface. The extended electromagnetic wave theory is utilized, and an impedance boundary conditions approach is employed to obtain characteristics equation. Effective mode index and attenuation under the different values of collisional frequency, plasma frequency and chirality in the THz frequency range for three cases for uniaxial chiral media are discussed. These results can be used in the field of photonics and integrated optics to fabricate nanophotonic devices in the THz frequency range.

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-September 2022, 46 samples from suspected CCHF patients were tested, with six (13%) showing positive RT-PCR results. Among the positive cases, all were male, aged 14-48 years among which four were butchers. Three CCHF patients succumbed to the disease. The complete S-M-L-gene fragments of four 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. This article is protected by copyright. All rights reserved.

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.

Prevalence of Gram positive bacteria in the affected individuals of Otitis media with effusion from the indigenous population of Southern Punjab, Pakistan: first report

Otitis media with effusion (OME) is a type of otitis media (OM) characterized by the presence of fluid behind intact tympanic membrane and is one of the most common diseases of early childhood. It is an infectious disease associated with the presence of many pathogenic bacteria in the middle ear of affected individuals. This study was aimed to determine the prevalence of Gram-positive bacteria from the middle ear of OME patients in the population of Southern Punjab, Pakistan. The incidence of OME under comprehensive healthcare setting was investigated in patients who consulted at the department of ear, throat and nose, Bahawal Victoria Hospital (BVH), Bahawalpur, from December, 2019 to May, 2021. Ear swabs were taken from affected and normal individuals. After culturing bacteria from the ear swabs, microscopic analysis and biochemical tests were performed to characterize the cultured Gram-positive bacteria. Out of 352 patients examined, 109 (30.9%) patients had OME. Age of the participants ranged from 14 to 50 years; individuals between the ages of 14 and 22 years had the highest infection rates, while individuals between 40 and 50 years had the lowest rate of infection. Tympanic membrane perforation, fever, cough, sore throat, ear pain and hearing problem showed association with symptoms of OME. Microscopic analysis and biochemical characterization showed the presence of streptococci and staphylococci in all the studied samples. The most frequently isolated bacteria were Streptococcus pneumoniae, Streptococcus pyogenes, and Staphylococcus aureus with percentage of 53.3%, 20% and 13.3% respectively. Enterococcus faecalis (6.6%) and Staphylococcus epidermidis (6.6%) were also identified in the studied samples. This study will help in the better medical administration of OME affected individuals.

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.

Genomic diversity of SARS-CoV-2 in Pakistan during the fourth wave of pandemic

The emergence of different variants of concern of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) has resulted in upsurges of coronavirus disease 2019 (COVID‐19) cases around the globe. Pakistan faced the fourth wave of COVID‐19 from July to August 2021 with 314,786 cases. To understand the genomic diversity of circulating SARS‐CoV‐2 strains during the fourth wave of the pandemic in Pakistan, this study was conducted. The samples from 140 COVID‐19‐positive patients were subjected to whole‐genome sequencing using the iSeq Sequencer by Illumina. The results showed that 97% (n = 136) of isolates belonged to the delta variant while three isolates belonged to alpha and only one isolate belonged to the beta variant. Among delta variant cases, 20.5% (n = 28) isolates were showing B.1.617.2 while 23.5% (n = 25), 17.59% (n = 19), 14.81% (n = 16), and 13.89% (n = 15) of isolates were showing AY.108, AY.43 AY.127, and AY.125 lineages, respectively. Islamabad was found to be the most affected city with 65% (n = 89) of delta variant cases, followed by Karachi (17%, n = 23), and Rawalpindi (10%, n = 14). Apart from the characteristic spike mutations (T19R, L452R, T478K, P681R, and D950N) of the delta variant, the sublineages exhibited other spike mutations as E156del, G142D, T95I, A222V, G446V, K529N, N532S, Q613H, and V483A. The phylogenetic analysis revealed the introductions from Singapore, the United Kingdom, and Germany. This study highlights the circulation of delta variants (B.1.617.2 and sublineages) during the fourth wave of pandemic in Pakistan.

Molecular identification and characterization of Measles virus genotypes circulating in Afghanistan, 2008-2018

The measles virus remains a leading cause of morbidity and mortality in children under five years of age. Molecular identification of circulating wild-type measles virus (MV) strains is a vital component of the measles elimination program. We received 159 oral swab samples from Afghanistan during 2008-2018. Viral RNA was extracted, followed by one-step RT-PCR and positive amplicons were subject to sequencing for genotype identification. Out of 159 total samples, 52% (83/159) were detected positive by RT-PCR. Genotype D4 was identified from 2.4% (2/83), genotype H1, 4.8% (4/83) and genotype B3, 92.7% (77/83) cases respectively. This article is protected by copyright. All rights reserved.

Genomic surveillance reveals the detection of SARS-CoV-2 delta, beta, and gamma VOCs during the third wave in Pakistan

SARS-CoV-2 variants of concern (VOCs) have emerged worldwide and gained significant importance due to their high transmissibility and global spread, thus meriting close monitoring. In Pakistan, limited information is available on circulation of these variants as the alpha variant has been reported the main circulating lineage. The current study was designed to detect and explore the genomic diversity of SARS-CoV-2 lineages circulating during the third wave of the pandemic in the indigenous population. From May 01 to June 09, 2021, a total of 16 689 samples were tested using TaqPath™ COVID-19 kit for the presence of SARS-CoV-2. Overall, 2562 samples (15.4%) were COVID-19 positive. Out of these positive samples, 2124 (12.7%) did not show the spike gene amplification (spike gene target failure ([SGTF]), whereas 438 (2.6%) showed spike gene amplification (non-SGTF). A subset (n = 58/438) of non-SGTF samples were randomly selected for whole-genome sequencing. Among VOCs, 45% (n = 26/58) were delta, 46% (n = 27/58) were beta, and one was gamma variant. The delta variant cases were reported mainly from Islamabad (n = 15; 58%) followed by Rawalpindi and Azad Kashmir (n = 1; 4% each). Beta variant cases originated mainly from Karachi (n = 8; 30%) and Islamabad (n = 11; 41%) and the gamma variant case was reported in a traveler from Italy. The delta, beta, and gamma variants possessed lineage-specific spike mutations. Notably, two rare mutations (E484Q and L5F) were found in the delta variant. Furthermore, in the beta variant, two significant rare non-synonymous spike mutations (A879S and K444R) were also reported. High prevalence of beta and delta variants in local population may increase the number of cases in the near future and provides an early warning to national health authorities to take timely decisions and devise suitable interventions to contain a possible fourth wave.

Chikungunya virus: Molecular epidemiology of nonstructural proteins in Pakistan

Chikungunya virus (CHIKV) is considered a public health problem due to its rapid spread and high morbidity. In 2016-2017 an outbreak of CHIKV was occurred in Pakistan but the data regarding the genomic diversity of CHIKV was not reported. Hence, the current study aimed to determine the genetic diversity of CHIKVs in Pakistan. A cross sectional study was carried out using sera of infected CHIKV patients (n = 1549) during the outbreak in Pakistan (2016-2018). Nucleotide sequencing of non-structural genes of CHIKV from eight isolates were performed followed by phylogenetic analysis using Bayesian method. Phylogenetic analysis suggested that the Pakistani CHIKV strains belonged to Indian Ocean Lineage (IOL) of genotype ECSA and C1.3a clade. Furthermore, the Pakistani isolates showed several key mutations (nsP2-H130Y, nsP2-E145D, nsP4-S55N and nsP4- R85G) corresponding to mutations reported in 2016 Indian strains of CHIKV. The molecular analysis revealed high evolutionary potential of CHIKV strains as well as better understanding of enhanced virulence and pathogenesis of this outbreak. The study highlights the need to continue surveillance in order to understand viral diversity over time and to devise preventive measures to limit diseases transmission in the region.

Genomic Epidemiology of SARS-CoV-2 in Pakistan

COVID-19 has swept globally and Pakistan is no exception. To investigate the initial introductions and transmissions of the SARS-CoV-2 in Pakistan, we performed the largest genomic epidemiology study of COVID-19 in Pakistan and generated 150 complete SARS-CoV-2 genome sequences from samples collected from March 16 to June 1, 2020. We identified a total of 347 mutated positions, 31 of which were over-represented in Pakistan. Meanwhile, we found over 1000 intra-host single-nucleotide variants (iSNVs). Several of them occurred concurrently, indicating possible interactions among them or coevolution. Some of the high-frequency iSNVs in Pakistan were not observed in the global population, suggesting strong purifying selections. The genomic epidemiology revealed five distinctive spreading clusters. The largest cluster consisted of 74 viruses which were derived from different geographic locations of Pakistan and formed a deep hierarchical structure, indicating an extensive and persistent nation-wide transmission of the virus that was probably attributed to a signature mutation (G8371T in ORF1ab) of this cluster. Furthermore, 28 putative international introductions were identified, several of which are consistent with the epidemiological investigations. In all, this study has inferred the possible pathways of introductions and transmissions of SARS-CoV-2 in Pakistan, which could aid ongoing and future viral surveillance and COVID-19 control.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Negative impact of the COVID-19 pandemic on routine childhood immunization: experience from Pakistan

The COVID-19 pandemic has disrupted routine childhood immunization coverage for many vaccine-preventable diseases at the global, regional and national levels. In Pakistan, more than 40 million children were unable to receive measles vaccination and around 50 million children missed polio vaccination in 2020. Moreover, more than 36 million children did not receive vitamin A supplementation in 2020. Implementing catch-up vaccination services to prevent outbreaks of vaccine-preventable diseases, particularly polio and measles, in malnourished children is urgently required to prevent a global resurgence of these deadly diseases.

Muhammad Suleman Rana and colleagues from the National Institute of Health in Pakistan discuss the urgent need to implement catch-up vaccination programmes for measles and polio to prevent resurgence of these deadly diseases.

Genetic and evolutionary analysis of SARS-CoV-2 circulating in the region surrounding Islamabad, Pakistan

Genomic epidemiology of Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) has provided global epidemiological insight into the COVID-19 pandemic since it began. Sequencing of the virus has been performed at scale, with many countries depositing data into open access repositories to enable in-depth global phylogenetic analysis. To contribute to these efforts, we established an Oxford Nanopore Technologies (ONT) sequencing capability at the National Institutes of Health (NIH), Pakistan. This study highlights multiple SARS-CoV-2 lineages co-circulating during the peak of a second COVID-19 wave in Pakistan (Nov 2020-Feb 2021), with virus origins traced to the United States of America and Saudi Arabia. Ten SARS-CoV-2 positive samples were used for ONT library preparation. Sequence and phylogenetic analysis determined that the patients were infected with lineage B.1.1.250, originally identified in the United Kingdom and Bangladesh during March and April of 2020, and in circulation until the time of this study in Europe, USA and Australia. Lineage B.1.261 was originally identified in Saudi Arabia with widespread local dissemination in Pakistan. One sample clustered with the parental B.1 lineage and the other with lineage B.6 originally from Singapore. In the future, monitoring the evolutionary dynamics of circulating lineages in Pakistan will enable improved tracing of the viral spread, changing trends of their expansion trajectories, persistence, changes in their demographic dynamics, and provide guidance for better implementation of control measures.

Detection of SARs-CoV-2 in wastewater using the existing environmental surveillance network: A potential supplementary system for monitoring COVID-19 transmission

The ongoing COVID-19 pandemic is caused by SARs-CoV-2. The virus is transmitted from person to person through droplet infections i.e. when infected person is in close contact with another person. In January 2020, first report of detection of SARS-CoV-2 in faeces, has made it clear that human wastewater might contain this virus. This may illustrate the probability of environmentally facilitated transmission, mainly the sewage, however, environmental conditions that could facilitate faecal oral transmission is not yet clear. We used existing Pakistan polio environment surveillance network to investigate presence of SARs-CoV-2 using three commercially available kits and E-Gene detection published assay for surety and confirmatory of positivity. A Two-phase separation method is used for sample clarification and concentration. An additional high-speed centrifugation (14000Xg for 30 min) step was introduced, prior RNA extraction, to increase viral RNA yield resulting a decrease in Cq value. A total of 78 wastewater samples collected from 38 districts across Pakistan, 74 wastewater samples from existing polio environment surveillance sites, 3 from drains of COVID-19 infected areas and 1 from COVID 19 quarantine center drainage, were tested for presence of SARs-CoV-2. 21 wastewater samples (27%) from 13 districts turned to be positive on RT-qPCR. SARs-COV-2 RNA positive samples from areas with COVID 19 patients and quarantine center strengthen the findings and use of wastewater surveillance in future. Furthermore, sequence data of partial ORF 1a generated from COVID 19 patient quarantine center drainage sample also reinforce our findings that SARs-CoV-2 can be detected in wastewater. This study finding indicates that SARs-CoV-2 detection through wastewater surveillance has an epidemiologic potential that can be used as supplementary system to monitor viral tracking and circulation in cities with lower COVID-19 testing capacity or heavily populated areas where door-to-door tracing may not be possible. However, attention is needed on virus concentration and detection assay to increase the sensitivity. Development of highly sensitive assay will be an indicator for virus monitoring and to provide early warning signs.

Proliferation of SARS-CoV-2 B.1.1.7 Variant in Pakistan-A Short Surveillance Account

The emergence of a more transmissible variant of SARS-CoV-2 (B1. 1.7) in the United Kingdom (UK) during late 2020 has raised major public health concerns. Several mutations have been reported in the genome of the B.1.1.7 variant including the N501Y and 69-70deletion in the Spike region that has implications on virus transmissibility and diagnostics. Although the B.1.1.7 variant has been reported by several countries, only three cases have been reported in Pakistan through whole-genome sequencing. Therefore, the objective of the study was to investigate the circulation of B.1.1.7 variant of concern (VOC) in Pakistani population. We used a two-step strategy for the detection of B.1.1.7 with initial screening through TaqPath^TM COVID-19 CE-IVD RT-PCR kit (ThermoFisher Scientific, Waltham, US) followed by partial spike (S) gene sequencing of a subset of samples having the spike gene target failure (SGTF). From January 01, 2021, to February 21, 2021, a total of 2,650 samples were tested for SARS-CoV-2 and 70.4% (n = 1,867) showed amplification of all the 3 genes (ORF, N, and S). Notably, 29.6% (n=783) samples have been SGTF that represented numbers from all the four provinces and suggest a rather low frequency during the first 3 weeks of January (n = 10, n = 13, and n = 1, respectively). However, the numbers have started to increase in the last week of January, 2021. During February, 726 (93%) cases of SGTF were reported with a peak (n = 345) found during the 3rd week. Based on the partial sequencing of SGTF samples 93.5% (n = 29/31) showed the characteristic N501Y, A570D, P681H, and T716I mutations found in the B.1.1.7 variant. In conclusion, our findings showed an upsurge of B.1.1.7 cases in Pakistan during February, 2021 affecting 15 districts and warranting large scale genomic surveillance, strengthening of laboratory network and implementation of appropriate control measures in the country.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Genetic Epidemiology Reveals 3 Chronic Reservoir Areas With Recurrent Population Mobility Challenging Poliovirus Eradication in Pakistan

BACKGROUND

Pakistan is among 3 countries endemic for wild poliovirus type 1 (WPV1) circulation that are still struggling for eradication of poliomyelitis. Active clinical and environmental surveillance with meticulous laboratory investigations provide insights into poliovirus transmission patterns and genomic diversity to inform decisions for strategic operations required to achieve eradication.

METHODS

We analyzed epidemiological and virological data to comprehend the current epidemiological status of WPV1 in Pakistan during 2015-2017. Stool specimens of patients with acute flaccid paralysis (AFP) and sewage samples collected from 60 environmental sites were tested. Viral culturing, intratypic differentiation by real-time polymerase chain reaction, and nucleic acid sequencing of the VP1 region of the poliovirus genome to determine genetic relatedness among WPV1 strains were applied.

RESULTS

Poliovirus isolates were grouped into 11 distinct clusters, which had ≥95% nucleotide homology in the VP1 coding region. Most of the poliovirus burden was shared by 3 major reservoirs: Karachi, Peshawar, and Quetta block (64.2% in 2015, 75.4% in 2016, and 76.7% in 2017).

CONCLUSIONS

Environmental surveillance reveals importations and pockets of unimmunized children that dictate intensive target mop-up campaigns to contain poliovirus transmission. A decrease in the number of orphan isolates reflects effective combination of AFP and environmental surveillance in Pakistan. The genetic data reflect sustained transmission within reservoir areas, further expanded by periodic importations to areas of high immunity reflected by immediate termination of imported viruses. Improved immunization coverage with high-quality surveillance is vital for global certification of polio eradication.

Whole-genome sequencing of SARS-CoV-2 reveals the detection of G614 variant in Pakistan

Since its emergence in China, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide including Pakistan. During the pandemic, whole genome sequencing has played an important role in understanding the evolution and genomic diversity of SARS-CoV-2. Although an unprecedented number of SARS-CoV-2 full genomes have been submitted in GISAID and NCBI, data from Pakistan is scarce. We report the sequencing, genomic characterization, and phylogenetic analysis of five SARS-CoV-2 strains isolated from patients in Pakistan. The oropharyngeal swabs of patients that were confirmed positive for SARS-CoV-2 through real-time RT-PCR at National Institute of Health, Pakistan, were selected for whole-genome sequencing. Sequencing was performed using NEBNext Ultra II Directional RNA Library Prep kit for Illumina (NEW ENGLAND BioLabs Inc., MA, US) and Illumina iSeq 100 instrument (Illumina, San Diego, US). Based on whole-genome analysis, three Pakistani SARS-CoV-2 strains clustered into the 20A (GH) clade along with the strains from Oman, Slovakia, United States, and Pakistani strain EPI_ISL_513925. The two 19B (S)-clade strains were closely related to viruses from India and Oman. Overall, twenty-nine amino acid mutations were detected in the current study genome sequences, including fifteen missense and four novel mutations. Notably, we have found a D614G (aspartic acid to glycine) mutation in spike protein of the sequences from the GH clade. The G614 variant carrying the characteristic D614G mutation has been shown to be more infectious that lead to its rapid spread worldwide. This report highlights the detection of GH and S clade strains and G614 variant from Pakistan warranting large-scale whole-genome sequencing of strains prevalent in different regions to understand virus evolution and to explore their genetic diversity.

Epidemiology of Chikungunya virus isolates 2016-2018 in Pakistan

The chikungunya virus (CHIKV) is a mosquito-transmitted alphavirus, which has infected millions of people in Africa, Asia, Americas, and Europe since it remerged in India and Indian Ocean regions in 2005-2006. The purpose of this study was to evaluate the genetic diversity and evolutionary changes in CHIKV from 2016 to 2018 in Pakistan. Blood specimens were collected and processed following the Centers for Disease Control and Prevention Trioplex Protocol. Sequencing and phylogenetic analysis of complete coding sequence of representative isolates from the CHIKV outbreak was carried out during December 2016 to July 2018, a total of 1549 samples were received, out of which 50% (n = 774) were found positive for CHIKV RNA. Mean age of chikungunya positive patients was 31.8 ± 15.7 years and most affected were between 21 and 40 years of age. The Pakistan CHIKV strains clustered with the Indian Ocean sublineage of East/Central/South African with cocirculation of some variants In the structural proteins region, two noteworthy changes (A226V and D284E) were observed in the membrane fusion glycoprotein E1. Key substitutions in the neutralizing epitopes site and a few changes indicative of adaptive to other insect cells were also detected in Pakistani strains. This study provides the emerging trend of CHIKV in the country for early identification of potential variants of high virulence and preventive measures for vector borne disease especially in the endemic areas.

Epidemiological and molecular investigation of a measles outbreak in Punjab, Pakistan, 2013-2015

Despite the availability of an effective vaccine, the measles virus continues to cause significant morbidity and mortality in children worldwide. Molecular characterization of wild-type measles strains is an invaluable component of epidemiological studies or surveillance systems that provides important information pertinent to outbreak linkages and transmission pathways. Serum samples and throat swabs were collected from suspected measles cases from the Punjab province of Pakistan (2013-2015) and further tested for measles immunoglobulin M (IgM) through enzyme-linked immunosorbent assay and reverse-transcriptase polymerase chain reaction for molecular characterization. Among the total of 5415 blood samples, 59% tested positive for measles IgM. Males had a higher infection rate (55%) than females (45%), and the highest frequency of positive cases (63%) was found in the age group of 0 to 5 years. Partial sequencing of the nucleoprotein gene showed that 27 strains belonged to the B3 genotype, whereas 2 viruses were identified as D4. On phylogenetic analysis, Pakistani B3 strains were found to be closely related to previously reported indigenous strains and those from neighboring countries of Iran and Qatar. This is the first report on the detection of the measles B3 genotype from Punjab, Pakistan. The current study shows a high burden of measles infections in Punjab province owing to poor routine immunization coverage in major cities. It is imperative that national health authorities adopt strategic steps on an urgent basis for improvement of routine immunization coverage. Molecular epidemiology of the measles viruses circulating in different parts of the country can provide useful data to manage future outbreaks.

Rotavirus surveillance in Pakistan during 2015-2016 reveals high prevalence of G12P[6]

The G12 rotavirus genotype has emerged globally since their first detection in 1987 from the Philippines; however it remains a rare cause of gastroenteritis in Pakistan. Rotavirus surveillance conducted during 2015-2016, assessed 3446 children <5 years hospitalized for gastroenteritis and found 802 (23.2%) positive on ELISA. Genotyping of a subset of positive samples (n = 319) revealed G12P[6] (11.28%) as the third most common G/P combination following G3P[8] (28.5%) and G1P[8] (12.5%); G2P[4] (10.65%) and G3P[6] (8.15%) were other frequently detected strains. Phylogenetic analysis of G12 strains from Pakistan revealed high genetic similarity to G12 strains from Italy, Thailand, Korea, and Great Britain as well as local strains within G12 lineage III. In conclusion, G12P[6] was a major contributor of RVA gastroenteritis in Pakistani children. Robust surveillance after the introduction of rotavirus vaccines will help determine the evolution of G12 and other circulating genotypes in the country.

Homozygous XYLT2 variants as a cause of spondyloocular syndrome

Spondyloocular syndrome (SOS) is a rare autosomal recessive, skeletal disorder. Two recent studies have shown that it is the result of biallelic sequence variants in the XYLT2 gene with pleiotropic effects in multiple organs, including retina, heart muscle, inner ear, cartilage, and bone. The XYLT2 gene encodes xylosyltransferase 2, which catalyzes the transfer of xylose (monosaccharide) to the core protein of proteoglycans (PGs) leading to initiating the process of PG assembly. SOS was originally characterized in 2 families A and B of Iraqi and Turkish origin, respectively. Using DNA from affected members of the same 2 families, we performed whole exome sequencing, which revealed 2 novel homozygous missense variants (c.1159C > T, p.Arg387Trp) and (c.2548G > C, p.Asp850His). Our findings extend the body of evidence that SOS is caused by homozygous variants in the XYLT2 gene. In addition, this report has extended the phenotypic description of SOS by adding follow-up data from 5 affected individuals in one of the two families, presented here.