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

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.


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.

AFP Cases and Contacts
The Central Virology Laboratory of the Israel Ministry of Health investigated 1020 T2 . 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 1020 T2 . 01 AFP cases, 427 T2 . 02 cases reported poliovirus-positive stools. Seventy-three T2 . 04 AFP cases had PV3(+) stools. There were 93 T2 . 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." https://polioeradication.org/tools-andlibrary/resources-for-polio-eradicators/gpei-tools-protocols-and-guidelines/, last accessed 1 September 2022.
Demographic data (Table 3) associated with AFP cases and/or contacts with PV3positive stools was obtained retrospectively from summarized case reports maintained by the Central Virology Laboratory. 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 55 T2 . 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.

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 53 T2 . 11 archived test tubes of poliovirus isolates from PV3(+) cases plus 2 T2 . 12 isolates from PV3(−)PV(+) cases yielding WGS for 55 T2 . 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).

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.

Results
AFP case reports from 1973-1988 (n = 1020 T2 . 01 ) were examined retrospectively. There were 85 T2 . 03 AFP case reports indicating that PV3 had been detected in stools from the case and/or ≥1 contact. PV3 was isolated exclusively from 73 T2 . 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 53 T2 . 11 PV3 isolates had been archived from cases or contacts of the 85 T2 . 03 AFP cases (Tables 2 and 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. Table 4. Isolates from cases and/or contacts, from which complete P1 PV3-positive sequences were obtained by WGS. Of the 55 T2 . 10 sequenced isolates, 43 T2 . 13 isolates from 34 T2 . 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 12 T2 . 14 isolates from 11 T2 . 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 43 T2 . 13 WPV3 isolates, 35 of43 t2 . 13 isolates were obtained from stools from 34 T2 . 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 12 T2 . 14 SL3s, 1 T2 . 16 (SL3_ISR-1978-33) came from a case reporting PV2 and PV3 isolates, 9 T2 . 15 from PV3 cases, 1 T2 . 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 7 T2 . 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 = 15 T2 . 20 ) were obtained from stools of non-Jewish Arab patients, 14 from Gaza and 1 from the West Bank/East Jerusalem region. Of the two T2 . 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 6 T2 . 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 8 T2 . 23 of the 12 T2 . 14 isolates from 11 T2 . 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.
OR PEER REVIEW 10 of 16 In silico analysis of nucleotides associated with reversion to neurovirulence, (1) U 472 C in 5 UTR, (2) U 2034 C (VP3 phe 91 ser), and (3) U 2493 C (VP1 thr 6 iso), 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.   1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1987 1988 1986 0 A1 Ye ar Number of isolates  Figure 1. The 23 T2 . 19 isolates of Lineage A fall under the dashed line labeled "A". A1, A2, and A3 are three sub-lineages within Lineage A.
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 85 T2 . 03 AFP cases where stools from the case and/or ≥ 1 contact contained PV3; (2) all 73 T2 . 04 AFP cases where PV3 was isolated from the case; (3) all 43 T2 . 13 AFP cases with complete WPV3 capsid sequences; (4) the 23 T2 . 19 cases with isolates that mapped to lineage A or a sub lineage of A; and (5) the 11 T2 . 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. Table 6. The OPV vaccination histories of cases with WPV3-positive and SL3-positive stool isolates. 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.

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 po-liovirus 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 23 T2 . 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 43 T2 . 13 whole genomic sequences were recovered from 58 T2 . 09 cases with ≥1 archived WPV3 isolates (Tables 2 and 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 sublineages 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 11 T2 . 14 AFP cases were SL3s, rather than WPV3s, by sequence-based ITD. One T2 . 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 8 T2 . 23 of the isolates from the 11 T2 . 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. 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.