Next Article in Journal
Transmission Dynamics of Monkeypox Virus in Nigeria during the Current COVID-19 Pandemic and Estimation of Effective Reproduction Number
Next Article in Special Issue
Enterovirus 71-Associated Infection in South Vietnam: Vaccination Is a Real Solution
Previous Article in Journal
Vaccine-Acquired SARS-CoV-2 Immunity versus Infection-Acquired Immunity: A Comparison of Three COVID-19 Vaccines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Epidemiology of Type 3 Poliovirus AFP Cases in Israel between 1973 and 1988: Whole Genome Sequencing of RNA Extracted Directly from Archived Stocks to Avoid Re-Culturing Neurovirulent Wild Poliovirus

1
Central Virology Laboratory, Public Health Service Laboratories Israel Ministry of Health at Sheba Medical Center, Tel Hashomer, Ramat Gan 52621, Israel
2
Department of Epidemiology and Preventive Medicine, School of Public Health, Sackler Faculty of Medicine, Tel Aviv University, Tel Aviv 69978, Israel
3
US Food and Drug Administration, Center for Biologics Evaluation and Research, Silver Spring, MD 20993, USA
*
Author to whom correspondence should be addressed.
Vaccines 2022, 10(12), 2154; https://doi.org/10.3390/vaccines10122154
Submission received: 2 September 2022 / Revised: 5 December 2022 / Accepted: 12 December 2022 / Published: 15 December 2022

Abstract

:
Background: Poliovirus post-eradication containment of wild-type 2 poliovirus (PV2) requires the destruction of all materials containing, or potentially containing, PV2. Acute flaccid paralysis (AFP) cases in Israel between 1973 and 1988 were caused by all three serotypes; thus, isolates from cases and case-contacts were either PV2 or potentially contaminated with PV2. Aims: To provide a proof-of-concept that whole genome sequences (WGS) of wild-type 3 poliovirus (PV3s) could be salvaged from the RNA extracted directly from archived poliovirus stocks avoiding re-amplification of neurovirulent viruses, we link WGSs to case histories and determine the phylogenetic relationships among the PV3s. Methods: Data retrieved from 427 poliovirus-positive cases reported between 1973 and 1988 identified 85 PV3-associated cases. A total of 71 archived PV3 isolates were available from PV3-positive cases and contacts. WGSs were obtained by NGS from cDNA libraries constructed from RNA extracted directly from archived viral stocks. Sequences were subjected to phylogenetic analysis and linked to case data. Results: WGSs were successfully constructed for 55 isolates. Phylogenetic analysis revealed the circulation of seven lineages of PV3. One lineage, with 23 isolates, presented as an outbreak of six-year duration. Isolates from six other lineages were consistent with subsequent separate introductions, sporadic cases, and limited transmission. Recombinant vaccine-like PV3 recombinants were isolated from some cases. Conclusions: Whole or near-whole genome sequence information, obtained from RNA extracted directly from the archived material, safely provided detailed genetic information linked to patient data from a time when limited sequence information was previously available and revealed the pattern of transmission of wild PV3 in Israel.

1. Introduction

There are three serotypes of poliovirus, poliovirus type 1 (PV1), type 2 (PV2), and type 3 (PV3). Two types of poliovirus vaccine were developed and used globally, since the 1950s to prevent infections with the three serotypes of poliovirus, as well as the subsequent development of poliomyelitis (acute flaccid paralysis, AFP, resulting from damage to Motor neurons and denervation of muscles in a small minority of individuals infected with poliovirus). One vaccine contained inactivated polio virus strains (IPV) [1], and the other contained live attenuated poliovirus strains (OPV) [2]. Both were trivalent, e.g., they contained a representative strain of each of the three serotypes.
Between 1957 and 1987, various preparations and schedules of OPV were used in Israel (exception: two of fourteen health districts exclusively used IPV between 1982 and 1987) [3,4]. A supplementary immunization of almost all individuals up to the age of 40 in 1988 with a dose of OPV in response to an outbreak of wild-type 1 poliovirus that had begun in 1987 prevented further infections and AFP caused by any of the three wild poliovirus serotypes [3,4,5].
Wild-type 2 poliovirus (WPV2) was declared eradicated in 2015 [6,7]. This, in turn, triggered two events. The first was a globally coordinated replacement of the trivalent oral polio vaccine (OPV), which contained representative attenuated strains of Type 1, Type 2, and Type 3 poliovirus, with bivalent polio vaccine (bOPV, oral polio vaccine that only contains representative attenuated strains of Type 1 and Type 3 poliovirus) [8]. The second event was the initiation of phase 2 of the global action plan (GAP III) [9]. Annexes 2, 3, and 6 of GAP III require the containment of all type 2 polio viruses in polio essential facilities (PEFs) and the destruction of all Type 2 polioviruses present in non-essential polio facilities or their transfer from the non-PEFs to PEFs [9]. Most polio laboratories in the WHO global polio laboratory network (GPLN), including The Israel National Center for Polio and non-Polio Enteroviruses, chose to remain non-PEFs, rather than undergo the difficult and expensive procedure of becoming annually certified as a PEF by the WHO.
The Israeli laboratory’s archived collection of 257 wild polioviruses that were isolated from clinical AFP cases and contacts between 1973 and 1988 that were transferred to a PEF at the FDA in White Oaks, MD, for whole genome sequencing (WGS). Whole genome sequences can be generated using next generation quenching (NGS) techniques. NGS employs a number of different strategies to sequence the equivalent of a whole genome in a single reaction without the need to design specific primers for the poliovirus isolate of interest [10]. The NGS strategy used in this study is based on fragmentation of RNA into short 150–250 nt segments, converting the fragments into cDNA using random primers, ligating short sequences of DNA containing targets for a single set of primers to the ends of all cDNA fragments, amplifying each fragment into separate clusters on a flow cell, and then separately determining, in parallel, for all clusters, which of four fluorescently-labeled nucleotides are added to the growing complimentary strand during successive cycles of elongation. Contigs, the formation of a longer sequence from overlapping shorter sequences, representing whole genome sequences are generated in silico by either aligning short sequences reads to reference sequences or by the de novo assembly of overlapping sequences. However, VP1 sequences remain the standard for molecular intratypic differentiation (ITD) of polioviruses.
For safety, polioviruses RNA for sequencing was extracted directly from these samples without any need to grow high titers of neurovirulent viruses in tissue cultures. Genetic and phylogenetic analysis, based on the WGS of representative isolates of PV3 presented here, together with a summary of the demographic and clinical data associated with AFP cases and/or contacts with PV3-positive stools, shed light on the epidemiology of poliovirus infections in Israel between 1973 and 1988. WGS also significantly increased the existing database of sequences of the other three capsid proteins, non-structural genes, and non-coding and regulatory region PV3s that circulated in the Middle East during these years.

2. Materials and Methods

2.1. AFP Cases and Contacts

The Central Virology Laboratory of the Israel Ministry of Health investigated 1020T2.01 AFP cases between 1973 and 1988. The annual incidence of AFP cases and estimated AFP rate per 100,000 children ≤15 years old is presented in Table 1. The list of isolates and AFP cases described in this manuscript is presented in Table 2. Note: Of the 1020T2.01 AFP cases, 427T2.02 cases reported poliovirus-positive stools. Seventy-threeT2.04 AFP cases had PV3(+) stools. There were 93T2.08 individuals with PV3(+) stools. This includes stools from PV3(+) cases and from PV3(+) contacts of PV3(+) cases and PV3(−) cases with ≥1 PV3(+) contact. The Global Polio Surveillance Action Plan, 2018–2020, of the Polio Global Eradication Initiative defines contacts of AFP cases as follows:
“AFP contact sampling is the collection and testing of stool samples from contacts of AFP cases. A contact of an AFP case is defined as a child (preferably younger than five years of age) who likely had direct contact with the AFP case in the week prior to the onset of paralysis and/or in the two-week period after onset of paralysis. …AFP contact sampling is done to increase the sensitivity of the surveillance system to detect circulating polioviruses (wild and/or vaccine-derived) and, during an outbreak, to gain a better understanding of the geographic extent of the transmission.”
Demographic data (Table 3) associated with AFP cases and/or contacts with PV3-positive stools was obtained retrospectively from summarized case reports maintained by the Central Virology Laboratory.

2.2. RNA Extraction from Archived PV3-Positive Isolates for Whole Genomic Sequencing (WGS)

RNA was extracted directly from 1 mL aliquots of archived PV3-positive isolates for two reasons: (1) for safety reasons to avoid the necessity for growing high titter stocks of neurovirulent wild polioviruses; and (2) the unknown viability of 71 archived stocks of PV3s isolated from stools of cases or contact that had been stored at −20 °C for many years. RNA extraction was performed on 1 mL aliquots of archived cell culture supernatant from the earliest passage that had been archived using a MagNA Pure LC2.0 Automatic extractor with MagNA Pure LC Total Nucleic Acid Isolation Kit-High Performance (Roche Diagnostics, Indianapolis, IN, USA) and eluted into 50 µL of elution buffer, according to manufacturer’s instructions. PV3-positive RNA suitable for whole genome sequencing was extracted from 55T2.10 samples. The remaining 16 samples were not subjected to RNA extraction, did not yield RNA suitable for WGS, or only yielded partial sequences of insufficient length to confirm the poliovirus serotype or use for phylogenetic analysis.

2.3. Whole Genome Sequencing (WGS)

cDNA libraries, ranging from 200 to 500 bp in size, were prepared from the 0.5 µg of the total RNA extracted from 53T2.11 archived test tubes of poliovirus isolates from PV3(+) cases plus 2T2.12 isolates from PV3(−)PV(+) cases yielding WGS for 55T2.10 PV3 isolates. Briefly, NEBNext mRNA Sample Prep Master Mix Set 1 (New England BioLabs, Ipswich, MA, USA) was used to prepare RNA libraries, according to the manufacturer’s protocol (NEB), as described by Laasri et al. [11]. RNA was fragmentated by a focused ultrasonicator (Covaris) to generate the fragments of optimal sizes (250–300 nt). cDNA was synthesized using SuperScript III Reverse Transcriptase (Invitrogen, Waltham, MA, USA) and random primers. The cDNA was converted into double stranded cDNA, followed by an end repair procedure (Klenow fragment, T4 polynucleotide kinase, and T4 polymerase) and was ligated to Illumina (San Diego, CA, USA) paired-end (PE) adaptors. A double AMPure bead selection step (Beckman Coulter, Brea, CA, USA) was used for size selection. The DNA libraries were expanded using 15 cycles of PCR, with multiplex indexed primers and purified by magnetic beads (Agencourt AMPure PCR purification system, BeckmanCoulter). Quality and size distribution of the DNA fragments was confirmed on a BioAnalyzer using a high sensitivity kit (Agilent Technologies, Inc., Santa Clara, CA, USA).
Deep sequencing was performed using MiSeq (Illumina), producing paired-end reads 250 nt long, or HiSeq2000 (Illumina) producing 101 nt long paired-end reads. The sequencing reads were analyzed by the in-house ‘swarm’ or High-Performance Integrated Virtual Environment (HIVE) software the DNA fragments [12]. A curated database of 500 reference enteroviruses was used as the template for aligning the short sequences reads to the reference sequences in the database. Contigs containing complete sequences of the coding sequences (P1 which encodes all four capsid proteins and P2 and P3 which encode the nonstructural genes) were generated. Insufficient depth of coverage at either end of the genome made reliable sequence reconstruction of the ends of the 5′ and 3′ coding regions difficult for some of the contigs. In samples where more than one poliovirus was present, the individual consensus sequences were generated using the algorithm described in Karagiannis et al. [13].
Isolates were named as follows: Their typical differentiation (“WPV3” for wild-type 3 poliovirus or “SL3” for Sabin 3-like); the 3-letter country code “ISR” for Israel; the 4-digit year of isolation and two-digit incremental number assigned to indicate the incremental order of the case from which they were isolated during that year; and an uppercase letter indicating the lineage (example: isolate WPV3_ISR_1973_17-A1 is a wild PV3 of lineage A1 and was isolated from the 17th AFP case in 1973). Sub-lineages of isolates from lineage A have an additional number and lowercase letter following the “A”). Genomic sequences were deposited in the GenBank/EMBL/DDBJ database (accession numbers OP137272-OP137324).

2.4. Phylogenetic Analysis

PV3 whole genome sequences were aligned with the Sabin 3 vaccine strain sequence (GenBank access number AY421821) using the Sequencher program v5.4 (Gene Codes Corp., Ann Arbor, MI, USA). Aligned sequences were trimmed in silico to produce aligned sequences of VP1, P1, and the entire coding region (P1-P2-P3) for phylogenetic analysis using the equivalent sequences of Sabin 3 (AY421821) as the outgroups.
Aligned sequences were exported to Clustal x [14], where the sequences were re-aligned, and a neighbor-joining tree created after bootstrapping for 1000 iterations. The resultant trees were visualized using NjPlot v2.3 (http://pbil.univ-lyon1.fr/software/njplot.html last accessed 14 November 2022.) and annotated using Intaglio 3.95 (intag [email protected]; https://pcmacstore.com/en/app/417852764/intaglio%20last%20accessed%2014 last accessed 14 November 2022). Aligned sequences for specific lineages were analyzed for recombination events using the SimPlot program v3.5 (https://sray.med.som.jhmi.edu/SCRoftware/ last accessed 14 November 2022.).

3. Results

AFP case reports from 1973–1988 (n = 1020T2.01) were examined retrospectively. There were 85T2.03 AFP case reports indicating that PV3 had been detected in stools from the case and/or ≥1 contact. PV3 was isolated exclusively from 73T2.04 cases. The remaining 12 cases had PV3-negative stools, but ≥1 PV3-positive contact (four poliovirus-negative cases with PV3(+) contacts and eight PV1-positive and/or PV2-positive-cases with ≥1 PV3-positive contact). Complete or near-complete whole genome sequences were obtained by NGS for 53T2.11 PV3 isolates had been archived from cases or contacts of the 85T2.03 AFP cases (Table 2 and Table 4). In addition, there was one PV1-positive case report, WPV3_ISR_1974-20_A2a, that yielded a unique WPV3 sequence, instead of a PV1 sequence by NGS.
Of the 55T2.10 sequenced isolates, 43T2.13 isolates from 34T2.13 cases or case-contacts, was classified as wild-type 3 polioviruses, WPV3s, since their VP1 sequence was >15% different from the VP1 sequence of Sabin 3 VP1 (AY421821). These isolates are to the right of the vertical dashed line in Figure 1. The remaining 12T2.14 isolates from 11T2.14 archived samples were classified as Sabin 3-like (SL3), since their VP1 sequence diverged from that of Sabin 3 vaccine strain by less than 1% (sequences to the left of the vertical dashed line in Figure 1. Of the 43T2.13 WPV3 isolates, 35 of43t2.13 isolates were obtained from stools from 34T2.13 AFP cases (two different WPV3 sequences were obtained from case 1983-01), and 8 WPV3 isolates were obtained from stools of contacts (6 from 3 of the sequenced cases and two from poliovirus-negative cases). Of the 12T2.14 SL3s, 1T2.16 (SL3_ISR-1978-33) came from a case reporting PV2 and PV3 isolates, 9T2.15 from PV3 cases, 1T2.17 (SL3-ISR-1988-25) came from a PV3 case that also contained an isolate that yielded an incomplete WPV3 whole genome sequence, and 2 isolates were obtained from one archived sample.
A bootstrapped, neighbor-joining phylogenetic tree with Sabin 3 (AY411821) as outgroup was prepared for P1 sequences (2634 nt encoding 878 amino acids) (Figure 1). WPV3 isolates sequences from cases and contacts was segregated into 7T2.18 lineages. The lineages of the PV3 isolates were arbitrarily designated with the letters A to G in alphabetical order, corresponding to the year of isolation of the first isolate in the lineage (name in red text in the figure). Simplot analysis of complete open reading frames (6618 nts) encoding the structural and non-structural genes revealed that the sequences in lineage A could be subdivided into three sub-lineages, A1, A2, and A3, and sub-lineages A2 and A3 could be sub-divided into two and three sub-sub-lineages, respectively, based on major sequence differences in their P2-P3 nonstructural gene coding regions, resulting from different recombination events (not shown). A1 case isolates (n = 15T2.20) were obtained from stools of non-Jewish Arab patients, 14 from Gaza and 1 from the West Bank/East Jerusalem region. Of the twoT2.21 A2 case isolates, both were obtained from non-Jewish Arab patients, one was from Gaza, and the other from the West Bank. For the 6T2.22, the A3 case isolates were also obtained from stools of non-Jewish Arab patients, eight from Gaza, and one from an Israeli health district.
The presumed recombinant regions within the P2-P3 genomic regions were excised in silico and subjected to an NCBI nucleotide BLAST search [15] of the GenBank nucleotide collection (nr/nt) in an attempt to identify the non-polio enterovirus that supplied the recombinant genome segment. It was not surprising that closest matches were ≤86% homologous with contemporary viruses, due to the paucity of the available sequence data. The exception was that sequences from then P2-P3 region of isolates WPV3_ISR_1986-04_F and WPV3_ISR_1988_01_G both shared up to 90% homology with the WPV3s isolated in Finland in 1984.
WGS of 8T2.23 of the 12T2.14 isolates from 11T2.14 cases with SL3 VP1 sequences revealed that they had recombined with sequences from OPV1, OPV2, and/or a non-polio enterovirus (NPEV) in their P2 and/or P3 regions (Figure 2). Seven isolates were obtained from Jewish children in Israeli health districts, with two from Arab children whose geographic location was not available and three from cases lacking information.
In silico analysis of nucleotides associated with reversion to neurovirulence, (1) U472C in 5′UTR, (2) U2034C (VP3 phe91ser), and (3) U2493C (VP1 thr6iso), are presented in Table 5. Four SL3 isolates had fully reverted to neurovirulence, six had reverted at two of three sites, and one had reverted at only one site. Of note, seven were isolates obtained within 30 days of a dose of OPV.
The annual incidence of the different lineages of WPV3 and SL3 isolates was graphed (Figure 3). Isolates of the A lineage, the lineage with the most isolates, continued to circulate between 1973 and 1978. The remaining lineages appeared sporadically over periods of one to three years from 1974 to 1988. The SL3 isolates appeared throughout the time period investigated.
The clinical and demographic data presented in Table 6, Table 7, Table 8 were calculated from case report summaries for the following five categories: (1) all 85T2.03 AFP cases where stools from the case and/or ≥ 1 contact contained PV3; (2) all 73T2.04 AFP cases where PV3 was isolated from the case; (3) all 43T2.13 AFP cases with complete WPV3 capsid sequences; (4) the 23T2.19 cases with isolates that mapped to lineage A or a sub lineage of A; and (5) the 11T2.14 cases that had SL3-like isolates. Summarized case reports were reviewed for vaccination histories of PV3-positive individuals. The vaccination histories and the interval between the last OPV dose and appearance of symptoms are presented in Table 6. Clinical symptoms from case reports where they were recorded are summarized in Table 7. Finally, kinetic changes in anti-polio IgG titers after onset of symptoms (seroconversion or ≥ 4-fold increases), when recorded in case reports, are summarized in Table 8.

4. Discussion

Our current database of time-dated poliovirus sequences reflects the developments in the field of genome sequencing. The transition from Sanger sequencing, through automation of Sanger sequencing, to next generation sequencing has enabled us to easily generate longer and longer sequences culminating in whole genome sequences. Many of the sequences available for the earliest isolates cover only 250 nucleotides, spanning the junction between capsid protein 1 (VP1) and protease 2a. This segment served as the initial standard for sequence-based intratypic differentiation (ITD), e.g., the characterization of poliovirus isolates as either Sabin strain vaccine-like (SL) or wild (WPV), based on comparison to the same region of the corresponding Sabin vaccine strain [16]. Automation of Sanger sequencing enabled ITD comparisons to shift to the 900–906 nucleotides of the complete VP1 sequence [17]. Isolates were classified as Sabin vaccine strain-like (SL) when VP1 sequences diverged from the Sabin vaccine strain by ≤1% for PV1s and PV3s and ≤0.6% for a PV2, as vaccine-derived poliovirus (VDPV) for divergence between 1–15% or 0.6–15% or wild poliovirus WPV, when differences were in excess of 15% [18]. The concatenation of sequences of overlapping genome segments generated in separate Sanger sequencing reactions using specific primers was a method for generating whole genome sequences. It was labor intensive, and in many cases, required the design of new primers whenever genomic sequences became insufficiently complimentary, due to individual substitutions of individual nucleotides or major differences resulting from genetic re-combinations with other polioviruses and/or non-polio enteroviruses.
Many laboratories kept collections of un-sequenced or partially sequenced wild polioviruses with linked case reports that were isolated before complete genome sequences could be easily generated by NGS. Containment of type 2 poliovirus (phase 2 of GAP III), invoked when type 2 wild poliovirus [ was declared eradicated in 2015 6,7,9], required the destruction of all samples containing or potentially containing PV2 or their transfer to a polio essential facility (PEF). Many laboratories chose to destroy samples from enteric infections, including those containing type 1 or 3 polioviruses, in addition to those known to contain type 2, since they could not rule out the possibility that they might also have contained some type 2 poliovirus. The Israel National Center for Polioviruses and non-Polio Enteroviruses, instead, chose to transfer aliquots of annotated polioviruses from AFP cases isolated between 1973 and 1988 to the PEF laboratory at the FDA in the US. This transfer had four aims. The first was to link wild poliovirus whole genome sequences with case histories to improve our understanding of the epidemiology of the last AFP cases caused by wild poliovirus infections in Israel [3,5,17]. The second was to generate complete or near complete genome sequences of polioviruses isolated from a time period when few sequences were publicly available for any polioviruses isolated anywhere. The third was that it would be possible to destroy potentially infected stocks once sequenced while “preserving” them in silico for further study since viable polioviruses can be generated de novo [19]. The fourth was to provide proof of concept that, as we demonstrated for archived samples of another enterovirus, echovirus 11 [11], we could retrieve WGS data from archived stocks of poliovirus, without regard to viability and without the need to re-amplify them in cell culture.
The study, described here, describes the WGS of PV3 polioviruses isolated from stools of AFP cases and/or contacts. Seven distinct lineages of WPV3 were isolated from stools of AFP cases between 1973 and 1988 in Israel based on the phylogenetic analysis of sequences of their P1 regions (Figure 1) coupled with additional sequence information from their P2-P3 regions. Annual AFP rates ≥2 per 100,000 children <15 years old throughout this period (Table 1) indicate that, most likely, all AFP cases were investigated each year. The number and temporal pattern of cases with WPV3 isolates belonging to lineages B through G are consistent with sporadic cases, due to separate introductions and chains of transmission lasting a few years at most. In contrast, the 23T2.19 AFP cases associated with isolates of lineage A that appeared for at least 6 years, starting in 1973, and clearly fit the pattern of an outbreak of prolonged duration. The duration of transmission may have been even longer if related cases occurred before 1973, the beginning of the period of observation. Moreover, members of all of the lineages may have been missed, since the only 43T2.13 whole genomic sequences were recovered from 58T2.09 cases with ≥1 archived WPV3 isolates (Table 2 and Table 4).
Phylogenetic analysis of P1 sequences (Figure 2) suggested that three sub-lineages of Lineage A co-circulated. This was confirmed by SimPlot analysis
SimPlot analysis of isolates of the complete open reading frames (P1, P2, and P3 genomic regions) confirmed the three lineages and revealed that both the A2 and A3 sub-lineages could each be further subdivided into two sub-sub-lineages based on different the recombinations in non-structural genes in P2 and/or P3. Both of the P2-P3 region recombinant sequences of sub-lineage A3 were different from the P2-P3 regions of isolates in sub-lineages A1 and A2 and likely occurred after A3 diverged from A1 and A2. The years during which cases were associated with A3a (1976–1977) and of A3b (1976 to 1978) isolates imply that both circulated in Israel for some time. The three sub-lineages of lineage A may have diverged during endogenous circulation in Israel, although it is more likely that they represent multiple simultaneous importations from a common external reservoir, where they had diverged and recombined. If sequences from such an external reservoir had been available, we could have inferred multiple introductions, rather than divergence during local circulation, if the sequences of each of the two isolates showed higher similarity to the sequences in the reservoir than to each other (see Manor et. al. [20] for an example where such an analysis was successfully applied to two WPV1s isolated from sewage in Gaza, one in February and one in August in 2000).
Isolates from 11T2.14 AFP cases were SL3s, rather than WPV3s, by sequence-based ITD. OneT2.17 of the SL3s was a minor component of a case where WPV3 was the major component. Such a large number of SL3s might be unexpected, since these were AFP cases. The association of the isolates with AFP may have been coincidental, given the short interval between onset of symptoms and a last dose of OPV in some of the cases (Table 6). However, it is possible that some of these SL3s contributed to or caused the AFP, since in 8T2.23 of the isolates from the 11T2.14 cases, SL3 had recombined with sequences from OPV1, OPV2 and/or a non-polio enterovirus (NPEV) (Figure 2), and in silico analysis indicated reversion to the neurovirulent genotype of some or all of the neurovirulence attenuation sites of OPV3. If they did contribute etiologically to the AFP, then in the absence of sequence confirmation, historical reports on the incidence of poliovirus-associated paralytic poliomyelitis (VAPP) may be an underestimation.
The clinical symptoms of cases with PV3-positive stools are listed in Table 7. Of potential interest is that cases associated with lineage A presented with paralysis and/or paresis, while the following symptoms were underrepresented, compared with the other AFP cases associated with PV3s: fever, rash, sore throat, and/or problem swallowing, diarrhea/nausea, respiratory/bulbar involvement, meningitis/encephalitis, unconsciousness, and general paresis.
Seroconversion after a poliovirus infection is generally defined a change from a seronegative antibody titter of <1:8 before infection to a titter ≥8 after infection or a 4-fold rise in antibody titer for individuals who were seropositive before vaccination [21]. However, the Polio Laboratory Manual 4th Edition states that:
“Testing for polio neutralizing antibodies is not recommended for routine use in the diagnosis of poliomyelitis. It has rarely been useful in clarifying questionable virus diagnoses and requires additional field and laboratory time. Interpreting serum antibody titres is difficult with widespread immunization, and the method does not differentiate between antibodies against wild and vaccine strains.”
Many of the case reports included kinetic studies of neutralizing IgG antibody titers.
These results, summarized in Table 8, confirm that seroconversion may be less useful for determining the etiology of a poliomyelitis case than serotyping the poliovirus isolated from the stools of the patient. However, vaccination with OPV prior to appearance of symptoms may confound the interpretation of results obtained by serotyping the poliovirus isolated from stools and confirmation by RT-PCR and WGS ITD.
In conclusion, (1) we provide proof of concept that we can obtain WGS of archived neurovirulent polioviruses, independent of viability by sequencing RNA extracted directly from the archived material, e.g., without the need to grow high titter cultures of the viruses; (2) we provide sequence-based information on the epidemiology of circulation of wild polioviruses; and (3) we provide whole or near-whole genome sequence information about the wild polioviruses that circulated between 1973 and 1988 and related them to case reports of the AFP cases from which they were isolated.

Author Contributions

Conceptualization: L.M.S. Formal analysis: K.C. and L.M.S. Investigation: M.L., R.H., T.Z., D.S., M.W., K.C. and L.M.S. Methodology: K.C., M.L., R.H., T.Z., M.W. and D.S. Software: K.C. Supervision: L.M.S., R.H., D.S., M.W., E.M. and K.C. Writing—original draft: L.M.S. Writing—review and editing: L.M.S. and K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by The Sheba Medical Center Ethical Review Board, Sheba Medical Center, Tel Hashomer, ISRAEL, PROTOCOL: SMC-2027-15; 18 March 2015.

Informed Consent Statement

The Sheba Medical Center Ethical Review Board approved this study (SMC-2027-15) and patient consent was waived due to removal of all links to personal details that could be used to identify individuals and anonymous analyses of data.

Data Availability Statement

Sequences were deposited in GenBank/EMBL/DDBJ (Access numbers OP137272 to OP137324.

Acknowledgments

We acknowledge the support of members of the Central Virology Laboratory and the Epidemiology Department of the Israel Ministry of Health.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Plotkin, S.A.; Vidor, E. Poliovirus vaccine-inactivated. In Vaccines, 5th ed.; Plotkin, S.A., Orenstein, W.A., Offit, P.A., Eds.; Saunders/Elsevier: Philadelphia, PA, USA, 2008; pp. 605–629. [Google Scholar]
  2. Sutter, R.W.; Kew, O.M.; Cochi, S.L. Poliovirus vaccine-live. In Vaccines, 5th ed.; Plotkin, S.A., Orenstein, W.A., Offit, P.A., Eds.; Saunders/Elsevier: Philadelphia, PA, USA, 2008; pp. 631–685. [Google Scholar]
  3. Swartz, T.A. The Epidemiology of Polio in Israel an Historical Perspective; Dyonon Pub. Ltd.: Tel Aviv, Israel, 2008. [Google Scholar]
  4. Shulman, L.M. Polio and Its Epidemiology. In Encyclopedia of Sustainability Science and Technology Vol. Infectious Diseases; Meyers, R.A., SpringerLink (Online service), Eds.; Springer: New York, NY, USA, 2020; pp. 1–73. [Google Scholar]
  5. Slater, P.E.; Orenstein, W.A.; Morag, A.; Avni, A.; Handsher, R.; Green, M.S.; Costin, C.; Yarrow, A.; Rishpon, S.; Havkin, O.; et al. Poliomyelitis outbreak in Israel in 1988: A report with two commentaries. Lancet 1990, 335, 1192–1195, discussion 1196–1198. [Google Scholar] [CrossRef]
  6. GPEI. Global Eradication of Wild Poliovirus Type 2 Declared. 2015. Available online: http://www.polioeradication.org/mediaroom/newsstories/Global-eradication-of-wild-poliovirus-type-2-declared/tabid/526/news/1289/Default.aspx (accessed on 1 September 2016).
  7. WHO. Two out of Three Wild Poliovirus Strains Eradicated. Global Eradication of Wild Poliovirus Type 3 Declared on World Polio Day. 24 October 2019. Available online: https://www.who.int/news-room/feature-stories/detail/two-out-of-three-wild-poliovirus-strains-eradicated (accessed on 1 March 2021).
  8. Previsani, N.; Tangermann, R.H.; Tallis, G.; Jafari, H.S. World Health Organization Guidelines for Containment of Poliovirus Following Type-Specific Polio Eradication-Worldwide. MMWR Morb. Mortal. Wkly. Rep. 2015, 64, 913–917. [Google Scholar] [CrossRef] [Green Version]
  9. WHO. WHO Global Action Plan to Minimize Poliovirus Facility-Associated Risk after Type-Specific Eradication of Wild Polioviruses and Sequential Cessation of Oral Polio Vaccine Use (GAPIII). 2015. Available online: http://www.polioeradication.org/Posteradication/Certification.aspx or https://apps.who.int/iris/handle/10665/208872 (accessed on 30 September 2021).
  10. Metzker, M.L. Sequencing technologies-the next generation. Nat. Rev. Genet. 2010, 11, 31–46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Laassri, M.; Zagorodnyaya, T.; Hassin-Baer, S.; Handsher, R.; Sofer, D.; Weil, M.; Karagiannis, K.; Simonyan, V.; Chumakov, K.; Shulman, L. Evolution of echovirus 11 in a chronically infected immunodeficient patient. PLoS Pathog. 2018, 14, e1006943. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Simonyan, V.; Mazumder, R. High-Performance Integrated Virtual Environment (HIVE) Tools and Applications for Big Data Analysis. Genes 2014, 5, 957–981. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Karagiannis, K.; Simonyan, V.; Chumakov, K.; Mazumder, R. Separation and assembly of deep sequencing data into discrete sub-population genomes. Nucleic Acids Res. 2017, 45, 10989–11003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sayers, E.W.; Beck, J.; Bolton, E.E.; Bourexis, D.; Brister, J.R.; Canese, K.; Comeau, D.C.; Funk, K.; Kim, S.; Klimke, W.; et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2021, 49, D10–D17. [Google Scholar] [CrossRef] [PubMed]
  16. Kew, O.M.; Mulders, M.N.; Lipskaya, G.Y.; de Silva, E.; Pallansch, M.A. Molecular Epidemiology of Polioviruses. Semin. Virol. 1995, 6, 401–405. [Google Scholar] [CrossRef]
  17. Shulman, L.M.; Handsher, R.; Yang, C.F.; Yang, S.J.; Manor, J.; Vonsover, A.; Grossman, Z.; Pallansch, M.; Mendelson, E.; Kew, O.M. Resolution of the pathways of poliovirus type 1 transmission during an outbreak. J. Clin. Microbiol. 2000, 38, 945–952. [Google Scholar] [CrossRef] [PubMed]
  18. GPEI. Classification and Reporting of Vaccine-Derived Polioviruses (VDPV) GPEI Guidelines. 2016. Available online: http://polioeradication.org/wp-content/uploads/2016/09/Reporting-and-Classification-of-VDPVs_Aug2016_EN.pdf (accessed on 1 January 2022).
  19. Cello, J.; Paul, A.V.; Wimmer, E. Chemical synthesis of poliovirus cDNA: Generation of infectious virus in the absence of natural template. Science 2002, 297, 1016–1018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Manor, Y.; Blomqvist, S.; Sofer, D.; Alfandari, J.; Halmut, T.; Abramovitz, B.; Mendelson, E.; Shulman, L.M. Advanced environmental surveillance and molecular analyses indicate separate importations rather than endemic circulation of wild type 1 poliovirus in Gaza district in 2002. Appl. Environ. Microbiol. 2007, 73, 5954–5958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Giri, S.; Kumar, N.; Dhanapal, P.; Venkatesan, J.; Kasirajan, A.; Iturriza-Gomara, M.; John, J.; Abraham, A.M.; Grassly, N.C.; Kang, G. Quantity of Vaccine Poliovirus Shed Determines the Titer of the Serum Neutralizing Antibody Response in Indian Children Who Received Oral Vaccine. J. Infect. Dis. 2018, 217, 1395–1398. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Neighbor-Joining Phylogenetic Tree of the P1 regions of PV3 isolates. The tree was constructed using unweighted pair group method with arithmetic mean (UPGMA). Sabin 3 (text in purple) was set as the outgroup. The vertical dashed line is set at a divergence of 15% from Sabin 1. Bootstrap values (n = 1000) are indicated near the internal nodes of the tree. Black arrows indicate the AFP case for an isolate from a contact of that case. The lineages are indicated by the letter or letters, following the year of isolation and the order in which the stool from the AFP case was collected. Text in red indicate the earliest isolate in each lineage.
Figure 1. Neighbor-Joining Phylogenetic Tree of the P1 regions of PV3 isolates. The tree was constructed using unweighted pair group method with arithmetic mean (UPGMA). Sabin 3 (text in purple) was set as the outgroup. The vertical dashed line is set at a divergence of 15% from Sabin 1. Bootstrap values (n = 1000) are indicated near the internal nodes of the tree. Black arrows indicate the AFP case for an isolate from a contact of that case. The lineages are indicated by the letter or letters, following the year of isolation and the order in which the stool from the AFP case was collected. Text in red indicate the earliest isolate in each lineage.
Vaccines 10 02154 g001
Figure 2. Recombination events in Sabin 3-like polio isolates. The colored rectangles in the P2, P3, and 3′ noncoding regions of the genomes indicate range and the donor of the recombined fragments. The numbers below the rectangles indicate the nucleotide positions at the beginning and ends of the recombined fragments, relative to Sabin 1.
Figure 2. Recombination events in Sabin 3-like polio isolates. The colored rectangles in the P2, P3, and 3′ noncoding regions of the genomes indicate range and the donor of the recombined fragments. The numbers below the rectangles indicate the nucleotide positions at the beginning and ends of the recombined fragments, relative to Sabin 1.
Vaccines 10 02154 g002
Figure 3. Annual incidence of PV3 lineages. The graph indicates the annual number of sequence-confirmed isolates belonging to the different lineages of type 3 poliovirus isolated from stools of AFP cases in Israel between 1973 and 1988. The letters A-G refer to the lineages in Figure 1. The 23T2.19 isolates of Lineage A fall under the dashed line labeled “A”. A1, A2, and A3 are three sub-lineages within Lineage A.
Figure 3. Annual incidence of PV3 lineages. The graph indicates the annual number of sequence-confirmed isolates belonging to the different lineages of type 3 poliovirus isolated from stools of AFP cases in Israel between 1973 and 1988. The letters A-G refer to the lineages in Figure 1. The 23T2.19 isolates of Lineage A fall under the dashed line labeled “A”. A1, A2, and A3 are three sub-lineages within Lineage A.
Vaccines 10 02154 g003
Table 1. Annual AFP incidence, AFP rates, and PV3-positive cases in Israel between 1973 and 1988.
Table 1. Annual AFP incidence, AFP rates, and PV3-positive cases in Israel between 1973 and 1988.
YearAFP CasesTotal
Pop
(×105)
Pop ≤ 15 Yrs (×105)AFP Rate/(×105) ≤15 YrsPV(+) AFP CasesPV3(+) AFP CasesPV(−) Cases PV3(+) ContactsPV(+) PV3(−) Cases PV3(+) Contacts
19738633389359.2261323
1974164342295817.169800
19757734939787.9381112
1976150357510011584912
197776365310237.425500
197872373810476.936400
197977383610747.239400
198070392210986.440100
198142397811143.810200
198232406411382.88200
19833541191153310400
198422420011761.92000
198523426611941.96100
198627433112132.27300
1987124407123414100
198855447712544.423501
Total1020 4277348
Table 2. Cases and isolates described in this manuscript.
Table 2. Cases and isolates described in this manuscript.
IDNumberDescription
AFP Cases
T2.011020AFP cases investigated between 1973 and 1988
T2.02427AFP cases with PV(+) stools
T2.0385AFP cases with PV3(+) stools from the index case and/or contacts
T2.0473AFP cases with PV3(+) positive stools
T2.058AFP PV(+)/PV3(−) cases with PV3(+) contacts
T2.064AFP PV(−) cases with PV3(+) contacts
Archived PV Isolates
T2.07257Archived PV isolates (Types 1, 2, and 3)
T2.0893PV3(+) isolates reported for the index case and/or contacts (isolates from 93 individuals)
T2.0971 of the 93Archived PV3 isolates (isolates: 70 labeled PV3; 1 labeled PV1) from 58 cases with archived PV3(+) samples from the index case and/or contact
Whole Genome Sequences (WGS)
T2.1055PV3(+) isolates with WGS from 52 WPV3(+) cases and 1 WPV1(+) case
T2.1153 of the 55PV3(+) WGS from 51 WPV3(+) cases and 1 WPV1(+) case
T2.122 of the 55WGS available from PV3(-), PV(+) cases
T2.1343 of 55WPV3(+) WGS from 42 archived isolates from 34 cases
T2.1412 of 55SL3(+) WGS from 11 archived isolates from 11 cases—one sample contained two isolates
T2.159 of 12Isolates from VP3(+) cases (including 2 isolates from one case)
T2.161 of 12Isolate from a case with a mixture of PVs
T2.171 of 12Isolate from a case with a major incomplete WPV3 WGS component
Genetic lineages
T2.187Lineages (arbitrarily labeled A to G)
T2.1923Isolates with WGS from Lineage A
T2.2015Isolates in sub lineage A1
T2.212Isolates in sub lineage A2
T2.226Isolates in sub lineage A3
T2.238SL3 a with recombinant genomes
a SL3 = Sabin-like 3—a type 3 polio isolate with a viral capsid protein gene sequence that has diverged from the sequence of the Sabin 3 vaccine seed strain sequence by <1%.
Table 3. Demographic data for cases and/or contacts with PV3-positive stools.
Table 3. Demographic data for cases and/or contacts with PV3-positive stools.
PV3(+) Cases and Contacts a (n = 93T2.08)PV3(+) Case b
(n = 73T2.04)
Ethnicity
Jewish2316
Non-Jewish (Arab)6956
Non-Jewish (Druse)11
Not listed00
Geographical location
Israel (all health districts)3926
West Bank/E. Jerusalem54
Gaza4541
Not listed32
Other (Jordan) 10
Gender
Male56 (60%)47 (64%)
Female33 (36%)23 (32%)
Not recorded/not clear4 (4%)3 (4%)
a cases and/or contacts with PV3(+) stools, plus PV3(−) cases with ≥1 PV3(+) contact with PV3(+) stools. b cases with PV3(+) stools.
Table 4. Isolates from cases and/or contacts, from which complete P1 PV3-positive sequences were obtained by WGS.
Table 4. Isolates from cases and/or contacts, from which complete P1 PV3-positive sequences were obtained by WGS.
Description of Archived SampleNumber
Archived
Number Sequenced
TotalWPV3SL3
Isolates from Cases
PV3-positive53T2.1143T2.133310 a
PV3-positive and PV1-positive2101
PV3-positive and PV2-positive1110
PV3-positive and PV1- and PV2-positive2000
Isolates from PV3-positive Contacts
PV-negative cases4110
PV1-positive cases3000
PV2-positive cases2000
PV3-positive cases7660
Anomalous result
PV3 sequence from a PV1-positive case 1110
a Ten cases with 11 isolates—two closely related SL3 sequences were obtained from one archived isolate.
Table 5. In silico analysis of reversion to neurovirulence of SL3 isolates.
Table 5. In silico analysis of reversion to neurovirulence of SL3 isolates.
5′UTRVP3VP1OPV DosesInterval
between Last OPV and Symptoms
Recombinant
Nucleotide
Substitution
U472CU2034C U2493C
Amino Acid
Substitution
n/aphe91ser thr6iso
SL3_ISR-1978-04Cpheiso1<30 daysNo
SL3_ISR-1978-08U472Cseriso118 daysYes
SL3_ISR-1978-33U472Cpheiso11 dayYes
SL3_ISR-1980-03aU472Cpheiso0unvaccinatedYes
SL3_ISR-1982-02U472Cseriso3176 daysYes
SL3_ISR-1983-05U472Cseriso144 daysYes
SL3_ISR-1983-06U472Cpheiso28 daysYes
SL3_ISR-1987-00U472Cpheiso1<30 daysNo
SL3_ISR-1988-13U472CpheisoUnkYearsYes
SL3_ISR-1988-15U472Cpheiso27 daysYes
SL3_ISR-1988-25U472Cseriso114 daysNo
Grey highlight indicates reversion site unchanged. <30 days = exact dates not recorded. Unk = 32 years old, vaccinated but not recorded.
Table 6. The OPV vaccination histories of cases with WPV3-positive and SL3-positive stool isolates.
Table 6. The OPV vaccination histories of cases with WPV3-positive and SL3-positive stool isolates.
OPV Vaccination History
PV3(+) (n = 85T2.03)PV3(+) Case a
(n = 73T2.04)
PV3(+) Cases with Sequenced
Isolates
PV3(+)
(n = 43T2.13)
WPV3(+)
A1 (n = 15T2.20)
WPV3(+)
A3 (n = 6T2.22)
SL3(+)
(n = 11T2.14)
Unvaccinated b18 b17 b9311
1 dose OPV18138306 c
2 doses OPV d873002
3 doses OPV e232112911
4 doses of OPV 1297040
Not recorded or incomplete f664001
Days between last OPV dose and appearance of symptoms
PV3(+) (n = 55T2.10/85T2.03)PV3(+) case b
(n = 43T2.13/73T2.04)
PV3(+) Cases with Sequenced
Isolates c
PV3(+)
(n = 30/43T2.13)
WPV3(+)
A1 (n = 12/15T2.20)
WPV3(+)
A3 (n = 5/6T2.22)
SL3(+)
(n = 6/11T2.14)
Average21924417516424345
Median11313013115533116
Minimum37724247
Maximum20062006968355421176
Number < 30 days19138215
a AFP cases with PV3(+) stools. Does not include PV3(−) patients with PV3(+) contacts and PV(+)/PV3(−) patients with PV3(+) contacts. b One child received 1 dose of IPV only. c Includes one case with WPV3 partial sequence. d One child received two doses of IPV plus two doses of OPV. e One child received 3 doses of OPV plus one dose of monovalent type 1OPV. f Children with no vaccination history recorded or who were vaccinated, but details of vaccine type and number of doses were not recorded.
Table 7. Summary of clinical symptom associated with PV3-positive cases.
Table 7. Summary of clinical symptom associated with PV3-positive cases.
SymptomPV3(+) (n = 85T2.03)PV3(+) Cases a
(n = 73T2.04)
Sequence Confirmed Cases
PV3(+)
(n = 43T2.13)
WPV3(+)
A1 (n = 15T2.20)
WPV3(+)
A3 (n = 6T2.22)
SL3 (n = 10T2.14) b
Fever14120000
Rash110000
Sore throat or problem swallowing220000
Diarrhea/nausea220001
Respiratory/bulbar443101
Meningitis/encephalitis320001
Unconscious/coma221100
Death330001
General paresis 22219213
Left leg only665110
Right leg only220000
Left and right legs only333101
Left arm only000000
Right arm only220000
Left and right arm only000000
Arms and legs (any
combination)
761001
One leg (not
designated)
220001
General paralysis5644251245
Left leg only232113641
Right leg only642101
Left and right legs
only
171510502
Left arm only000000
Right arm only000000
Left and right arm only 000000
Arms and legs (any combination)1040000
Misc (palsy, arthritis, GB, otitis, facialis)442002
a AFP cases with PV3(+) stools. Does not include PV3(−) patients with PV3(+) contacts and PV(+)/PV3(−) patients with PV3(+) contacts. b Data available for 10 of the 11 cases.
Table 8. The number of sequence-confirmed, PV3-positive cases with >1:8 anti-PV3 IgG titers, <1:8 titers of anti-PV2 IgG and <1:8 titers of anti-PV3 IgG, and cases that seroconverted or had a ≥4-fold increase in anti-PV3 IgG titters.
Table 8. The number of sequence-confirmed, PV3-positive cases with >1:8 anti-PV3 IgG titers, <1:8 titers of anti-PV2 IgG and <1:8 titers of anti-PV3 IgG, and cases that seroconverted or had a ≥4-fold increase in anti-PV3 IgG titters.
LineagePV3
Total
(n = 43T2.13)
Positive a or ≥4-Fold Increase in Anti-PV3 IgG>Anti-PV3 IgG or ≥4-Fold Increase in IgG Not Done or Not Recorded
Positive But No
Increase a
Sero-Conversion
or ≥4-Fold Increase
Total Anti-PV1 IgGAnti-PV2 IgG
A2359 b11 b216
A11556 b11 b201
A22000002
A36033013
B31231 c00
C3033000
D1000001
E1011000
F2112000
G11011 d00
SL311 e03 3032
SL33011011
SL3rec8 e022031
PV3-positive Not sequenced30211136 f9 g7
a Anti-PV3 >1:8; anti-PV1; and anti-PV2 both <1:8. b In one case, the endpoint titers for anti-PV3 were ≥1:16,384, and the patients had 1:1024 titers against the other two serotypes. c One isolate PV3-positive by NGS was recorded as PV1- and PV3-positive and had >4-fold increases in anti-PV1 (1:256–>1:16,384) and anti-PV3 antibodies (<1:4–>1:4096). d O < PV1 (1:4096); anti-PV2 (>1:16,384); anti-PV3(1:16,384). e one SL3 was the minor component with a WPV3 as the major component. f One isolate had a 1:16,384 titter against PV1, seroconverted to 1:16,384 for anti-PV2, and was only 1:256 for anti-PV3, another seroconverted (1:4096–>1:10,384) for anti-PV1 and was less than 1:8 for anti-PV3. g 3 cases had ≥1:16,384 anti-PV2 titers and lower anti-PV3 titers and in a 4th case both anti-PV2 and PV3 were ≥1:16,384.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shulman, L.M.; Laassri, M.; Handsher, R.; Zagorodnyaya, T.; Sofer, D.; Weil, M.; Mendelson, E.; Chumakov, K. Epidemiology of Type 3 Poliovirus AFP Cases in Israel between 1973 and 1988: Whole Genome Sequencing of RNA Extracted Directly from Archived Stocks to Avoid Re-Culturing Neurovirulent Wild Poliovirus. Vaccines 2022, 10, 2154. https://doi.org/10.3390/vaccines10122154

AMA Style

Shulman LM, Laassri M, Handsher R, Zagorodnyaya T, Sofer D, Weil M, Mendelson E, Chumakov K. Epidemiology of Type 3 Poliovirus AFP Cases in Israel between 1973 and 1988: Whole Genome Sequencing of RNA Extracted Directly from Archived Stocks to Avoid Re-Culturing Neurovirulent Wild Poliovirus. Vaccines. 2022; 10(12):2154. https://doi.org/10.3390/vaccines10122154

Chicago/Turabian Style

Shulman, Lester M., Majid Laassri, Rachel Handsher, Tatiana Zagorodnyaya, Danit Sofer, Merav Weil, Ella Mendelson, and Konstantin Chumakov. 2022. "Epidemiology of Type 3 Poliovirus AFP Cases in Israel between 1973 and 1988: Whole Genome Sequencing of RNA Extracted Directly from Archived Stocks to Avoid Re-Culturing Neurovirulent Wild Poliovirus" Vaccines 10, no. 12: 2154. https://doi.org/10.3390/vaccines10122154

APA Style

Shulman, L. M., Laassri, M., Handsher, R., Zagorodnyaya, T., Sofer, D., Weil, M., Mendelson, E., & Chumakov, K. (2022). Epidemiology of Type 3 Poliovirus AFP Cases in Israel between 1973 and 1988: Whole Genome Sequencing of RNA Extracted Directly from Archived Stocks to Avoid Re-Culturing Neurovirulent Wild Poliovirus. Vaccines, 10(12), 2154. https://doi.org/10.3390/vaccines10122154

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop