Automated detection and classification of polioviruses from nanopore sequencing reads using piranha.

Widespread surveillance, rapid detection, and appropriate intervention will be critical for successful eradication of poliovirus. Using deployable next-generation sequencing (NGS) approaches, such as Oxford Nanopore Technologies' MinION, the time from sample to result can be significantly reduced compared to cell culture and Sanger sequencing. We developed piranha (poliovirus investigation resource automating nanopore haplotype analysis), a 'sequencing reads-to-report' solution to aid routine poliovirus testing of both stool and environmental samples and alleviate the bioinformatic bottleneck that often exists for laboratories adopting novel NGS approaches. Piranha can be used for efficient intratypic differentiation of poliovirus serotypes, for classification of Sabin-like polioviruses, and for detection of wild-type and vaccine-derived polioviruses. It produces interactive, distributable reports, as well as summary comma-separated values files and consensus poliovirus FASTA sequences. Piranha optionally provides phylogenetic analysis, with the ability to incorporate a local database, processing from raw sequencing reads to an interactive, annotated phylogeny in a single step. The reports describe each nanopore sequencing run with interpretable plots, enabling researchers to easily detect the presence of poliovirus in samples and quickly disseminate their results. Poliovirus eradication efforts are hindered by the lack of real-time detection and reporting, and piranha can be used to complement direct detection sequencing approaches.

Comparison of Eleven RNA Extraction Methods for Poliovirus Direct Molecular Detection in Stool Samples.

Direct detection by PCR of poliovirus RNA in stool samples provides a rapid diagnostic and surveillance tool that can replace virus isolation by cell culture in global polio surveillance. The sensitivity of direct detection methods is likely to depend on the choice of RNA extraction method and sample volume. We report a comparative analysis of 11 nucleic acid extraction methods (7 manual and 4 semiautomated) for poliovirus molecular detection using stool samples (n = 59) that had been previously identified as poliovirus positive by cell culture. To assess the effect of RNA recovery methods, extracted RNA using each of the 11 methods was tested with a poliovirus-specific reverse transcription-quantitative PCR (RT-qPCR), a pan-poliovirus RT-PCR (near-whole-genome amplification), a pan-enterovirus RT-PCR (entire capsid region), and a nested VP1 PCR that is the basis of a direct detection method based on nanopore sequencing. We also assessed extracted RNA integrity and quantity. The overall effect of extraction method on poliovirus PCR amplification assays tested in this study was found to be statistically significant (P < 0.001), thus indicating that the choice of RNA extraction method is an important component that needs to be carefully considered for any diagnostic based on nucleic acid amplification. Performance of the methods was generally consistent across the different assays used. Of the 11 extraction methods tested, the MagMAX viral RNA isolation kit used manually or automatically was found to be the preferable method for poliovirus molecular direct detection considering performance, cost, and processing time. IMPORTANCE Poliovirus, the causative agent of poliomyelitis, is a target of global eradication led by the World Health Organization since 1988. Direct molecular detection and genomic sequencing without virus propagation in cell culture is arguably a critical tool in the final stages of polio eradication. Efficient recovery of good-quality viral RNA from stool samples is a prerequisite for direct detection by nucleic acid amplification. We tested 11 nucleic acid extraction methods to identify those facilitating sensitive, fast, simple, and cost-effective extraction, with flexibility for manual and automated protocols considered. Several different PCR assays were used to compare the recovered viral RNA to test suitability for poliovirus direct molecular detection. Our findings highlight the importance of choosing a suitable RNA extraction protocol and provide useful information to diagnostic laboratories and researchers facing the choice of RNA extraction method for direct molecular virus detection from stool.

Molecular detection and genetic characterization of human metapneumovirus strains circulating in Islamabad, Pakistan.

Lower respiratory illness is one of the leading causes of death among children in low- and high-income countries. Human metapneumovirus (hMPV) is a key contributor to respiratory illnesses commonly reported among children and causes serious clinical complications ranging from mild respiratory infections to severe lower respiratory tract anomalies mainly in the form of bronchiolitis and pneumonia. However, due to the lack of a national surveillance system, the clinical significance of hMPV remains obscure in the Pakistani population. This study was conducted to screen throat swabs samples collected from 127 children reported with respiratory symptoms at a tertiary care hospital in Islamabad. Out of 127, 21 (16.5%) samples were positive for hMPV with its genotype distribution as A2a (10%), A2b (20%), B1 (10%), and B2 (60%). Phylogenetic analysis showed that the hMPV viruses were closely related to those reported from neighboring countries including India and China. This work will contribute to a better understanding of this virus, its diagnosis, and the handling of patients in clinical setups. Further studies at a large-scale are warranted for a better understanding of the disease burden and epidemiology of hMPV in Pakistan.

Conserved Antigenic Structure of Contemporary Wild Poliovirus Type 1 Strains Endemic in Pakistan.

BACKGROUND: Elimination of poliovirus in Pakistan and Afghanistan is challenged by notions against the role of oral poliovirus vaccine (OPV) in eradicating contemporary wild poliovirus (WPV) strains. METHODS: A total of 1055 WPV type 1 (WPV1) strains isolated between 2013 and 2018 were categorized into 68 antigenic groups and tested for neutralization by OPV-derived antibodies. Molecular docking was conducted to determine neutralization efficiency of antibodies against WPV. The clinical significance of WPV1 variants was assessed to ascertain their role in patient outcomes. RESULTS: We found that 88% of WPV1 strains isolated from paralytic children belonged to a single antigenic lineage identical to the WPV1 strain detected in 1993. WPV1 antigenic variants were effectively neutralized by OPV-derived antibodies, with geometric mean titers comparable to the neutralization titers found for 3 strains in OPV (OPV1-3, 7.96-9.149 [95% confidence interval, 6.864-10.171]; WPV1 strains, 7.542-8.786 [6.493-9.869]). Docking examination underscored a strong antigen-antibody interaction despite variations within the viral protein 1 epitopes. There was no significant association (P = .78) with clinical prognosis among patients infected with antigenically diverse WPV1 strains and patient outcomes, including death. CONCLUSIONS: Our findings substantiate the robustness of OPV for neutralizing the contemporary WPV1 strains endemic in Pakistan and Afghanistan. Vaccination coverage must be augmented to achieve early eradication.

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.

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.

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.

Rapid and Sensitive Direct Detection and Identification of Poliovirus from Stool and Environmental Surveillance Samples by Use of Nanopore Sequencing.

Global poliovirus surveillance involves virus isolation from stool and environmental samples, intratypic differential (ITD) by PCR, and sequencing of the VP1 region to distinguish vaccine (Sabin), vaccine-derived, and wild-type polioviruses and to ensure an appropriate response. This cell culture algorithm takes 2 to 3 weeks on average between sample receipt and sequencing. Direct detection of viral RNA using PCR allows faster detection but has traditionally faced challenges related to poor sensitivity and difficulties in sequencing common samples containing poliovirus and enterovirus mixtures. We present a nested PCR and nanopore sequencing protocol that allows rapid (<3 days) and sensitive direct detection and sequencing of polioviruses in stool and environmental samples. We developed barcoded primers and a real-time analysis platform that generate accurate VP1 consensus sequences from multiplexed samples. The sensitivity and specificity of our protocol compared with those of cell culture were 90.9% (95% confidence interval, 75.7% to 98.1%) and 99.2% (95.5% to 100.0%) for wild-type 1 poliovirus, 92.5% (79.6% to 98.4%) and 98.7% (95.4% to 99.8%) for vaccine and vaccine-derived serotype 2 poliovirus, and 88.3% (81.2% to 93.5%) and 93.2% (88.6% to 96.3%) for Sabin 1 and 3 poliovirus alone or in mixtures when tested on 155 stool samples in Pakistan. Variant analysis of sequencing reads also allowed the identification of polioviruses and enteroviruses in artificial mixtures and was able to distinguish complex mixtures of polioviruses in environmental samples. The median identity of consensus nanopore sequences with Sanger or Illumina sequences from the same samples was >99.9%. This novel method shows promise as a faster and safer alternative to cell culture for the detection and real-time sequencing of polioviruses in stool and environmental samples.

Emergence of Chikungunya Virus, Pakistan, 2016-2017.

During December 2016-May 2017, an outbreak of chikungunya virus infection occurred across Pakistan. The East/Central/South African genotype was predominant. This study provides baseline data on the virus strain and emphasizes the need for active surveillance and implementation of preventive interventions to contain future outbreaks.

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.

High prevalence of G3 rotavirus in hospitalized children in Rawalpindi, Pakistan during 2014.

Rotavirus A species (RVA) is the leading cause of severe diarrhea among children in both developed and developing countries. Among different RVA G types, humans are most commonly infected with G1, G2, G3, G4 and G9. During 2003-2004, G3 rotavirus termed as "new variant G3" emerged in Japan that later disseminated to multiple countries across the world. Although G3 rotaviruses are now commonly detected globally, they have been rarely reported from Pakistan. We investigated the genetic diversity of G3 strains responsible RVA gastroenteritis in children hospitalized in Rawalpindi, Pakistan during 2014. G3P[8] (18.3%; n = 24) was detected as the most common genotype causing majority of infections in children less than 06 months. Phylogenetic analysis of Pakistani G3 strains showed high amino acid similarity to "new variant G3" and G3 strains reported from China, Russia, USA, Japan, Belgium and Hungary during 2007-2012. Pakistani G3 strains belonged to lineage 3 within sub-lineage 3d, containing an extra N-linked glycosylation site compared to the G3 strain of RotaTeqTM. To our knowledge, this is the first report on the molecular epidemiology of G3 rotavirus strains from Pakistan and calls for immediate response measures to introduce RV vaccine in the routine immunization program of the country on priority.

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.

Outbreak of dengue virus type-3 in Malakand, Pakistan 2015; A laboratory perspective.

An outbreak of dengue fever was reported in Malakand district, Khyber Pakhtunkhwa (KP) province of Pakistan during 2015. Detection of viral RNA by real-time PCR confirmed dengue virus serotype-3 (DENV-3) to be the causative agent causing the outbreak. Phylogenetic analysis based on partial E-NS1 gene sequences showed that the DENV-3 viruses belonged to genotype III with maximum homology with the dengue-3 strains previously reported from Pakistan and India. Our current report provides updated information on molecular epidemiology and phylogenetic analysis of dengue virus serotypes responsible for 2015 outbreak in KP.

Respiratory syncytial virus and influenza are the key viral pathogens in children <2 years hospitalized with bronchiolitis and pneumonia in Islamabad Pakistan.

Pneumonia remains a leading cause of morbidity and mortality in developing countries. Comprehensive surveillance data are needed to review the prevention and control strategies. We conducted active surveillance of acute lower respiratory infections among children aged <2 years hospitalized at two hospitals of Islamabad, Pakistan. Viral etiology was determined using real-time PCR on respiratory specimens collected during March 2011-April 2012. The overall mean age was 7.83 ± 5.25 months while no statistical difference between age or sex distribution of patients with positive and negative viral etiology (p > 0.05). The average weight of the study group was 6.1 ± 2.25 kg. ≥1 viral pathogens were detected in 75% cases. Major respiratory viruses included RSV-A: 44%, RSV-B: 23%, Influenza-A: 24.5%, Influenza-B: 7%, Adenovirus: 8.4% and HmPV: 5.2%. A single, dual or multiple viral pathogens were detected in 43%, 27% and 5.2% patients respectively. Common symptoms were cough (95%), apnoea (84%), fever (78%), wheeze (64.5%), nasal congestion (55%) and rhinorrhea (48%). Among the RSV positive cases, 2-6 months age group had highest detection rate for RSV-A (30%, n = 21/69) and RSV-B (20%, n = 14/69) while patients infected with Influenza-A were in 2.1-6 months age group (61%, 23/38). Statistically significant difference was observed between RSV-positive and negative cases for nutrition status (p = 0.001), cigarette/wood smoke exposure (p = 0.001) and concomitant clinical findings. Most patients had successful outcome on combination therapy with bronchodilators, inhaled steroids and antibiotics. Our findings underscore high burden of ALRI in Pakistan. Interventions targeting viral pathogens coupled with improved diagnostic approaches are critical for better prevention and control.

Dengue Virus Serotypes Circulating in Khyber Pakhtunkhwa Province, Pakistan, 2013-2015.

From 2013 to 2015, the National Institute of Health, Pakistan, received 1,270 blood samples of suspected dengue cases reported from inpatient and outpatient departments of various hospitals in Khyber Pakhtunkhwa (KPK) province. In this study, we determined the circulating dengue virus (DENV) serotypes using real-time reverse transcriptase (RT)-PCR to understand the serotype-based epidemiology of DENV. All four serotypes (DENV-1 [6%], DENV-2 [33%], DENV-3 [47%], and DENV-4 [0.1%]) were found circulating during the study period. Our findings suggest the need for an active surveillance system coupled with the laboratory diagnosis, especially in the chronic endemic areas of the country. Public awareness programs are needed for effective control and prevention of outbreaks in the future.

NS1 antigen: A new beam of light in the early diagnosis of dengue infection.

OBJECTIVE: To evaluate NS1 antigen detection ELISA for the early laboratory diagnosis of dengue virus infection. METHODS: The present study was conducted to evaluate the overall positivity of NS1 antigen detection ELISA and its comparison with viral RNA detection via real time PCR and IgM antibodies detection by ELISA. RESULTS: A total of 1270 serum samples were tested 86% (1097/1270) were detected positive by one or more than one diagnostic test. Out of 1 270, 64% (807/1270) were positive by NS1 ELISA and 52% (662/1270), 51% (646/1270) were positive by real-time RT-PCR and IgM ELISA respectively. CONCLUSIONS: NS1 antigen detection ELISA is highly suitable diagnostic tools and it also has great value for use in outbreak and epidemic situation.

Dengue outbreak in Swat and Mansehra, Pakistan 2013: An epidemiological and diagnostic perspective.

OBJECTIVE: To high light some epidemiological, clinical and diagnostic features of dengue fever during an outbreak and the role of different diagnostic techniques to achieve the highest level of accuracy in results. METHODS: Blood samples (n = 323) were collected along with epidemiological and clinical data from suspected dengue patients who visited different hospitals in Swat and Mansehra district of Pakistan between May-November 2013 during a dengue outbreak. Samples were tested for the detection of viral nucleic acid by real-time PCR, non structural protein-1 (NS1) antigen and IgM antibodies by ELISA. RESULTS: Out of 323 cases with clinical dengue infection, 304 were positive by one or more diagnostic parameter; 201 samples were positive by real-time PCR, 209 were positive by NS1 ELISA and 190 were positive by IgM antibodies. Sensitivities of real-time PCR and NS1 ELISA were comparable for early diagnosis of dengue virus infection, IgM antibody detection assay was found useful for the diagnosis in the samples collected later than day 5 of onset. CONCLUSIONS: The use of real-time PCR or detection of non structural protein NS1 by ELISA followed by IgM antibodies detection can be recommended for early diagnosis of dengue virus infection with a high level of accuracy.