{'fraction_non_alphanumeric': 0.03314636911510233, 'fraction_numerical': 0.025359508323885675, 'mean_word_length': 4.730723200547102, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Quantitative RT-PCR Assays for Quantification of Undesirable Mutants in the Novel Type 2 Oral Poliovirus Vaccine

Emergence of mutations is an inherent property of RNA viruses with several implications for their replication, pathogenesis, and evolutionary adaptation. Oral poliovirus vaccine (OPV), developed by Albert Sabin, is composed of live attenuated polioviruses of three serotypes that can revert to neurovirulence during replication in cell culture and in vaccine recipients. Recently, a new modified variant of Sabin 2 virus was developed by introducing changes in its genome, making it more genetically stable to prevent the reversion. The new strain was used to manufacture novel OPV2 (nOPV2), which was approved by the World Health Organization for emergency use to stop outbreaks caused by circulating vaccine-derived poliovirus (cVDPV2). Manufacture of this improved vaccine requires close attention to the genetic heterogenicity to ensure that the levels of the undesirable mutations are limited. Preliminary studies using whole-genome Illumina sequencing (NGS) identified several genomic sites where mutations tend to occur with regularity. They include VP1-I143T amino acid change at the secondary attenuation site; VP1-N171D, a substitution that modestly increases neurovirulence in mice; and VP1-E295K, which may reduce the immunogenicity of the nOPV2. Therefore, to ensure the molecular consistency of vaccine batches, the content of these mutants must be quantified and kept within specifications. To do this, we have developed quantitative, multiplex, one-step reverse-transcriptase polymerase chain reactions (qmosRT-PCRs) as simple methods for quantification of these mutations. Each method uses specific short TaqMan probes with different dyes for the analysis of both mutants and non-mutants in the same sample. The quantification is done using calibration curves developed using validated reference materials. To evaluate the sensitivity and the linearity of the qmosRT-PCR method, the mutant viruses were spiked in non-mutant viruses, and nOPV2 batches were used to validate the method. The spiked samples and the nOPV2 batches were analyzed by qmosRT-PCR and NGS assays. The results showed that qmosRT-PCR is sensitive enough to detect around 1% of mutants. The percentages of mutants determined by qmosRT-PCR correlate well with the results of the NGS. Further, the analysis of the nOPV2 batches showed that the results of qmosRT-PCR correlated well with the results of NGS. In conclusion, the qmosRT-PCR is a specific, sensitive, and linear method. It could be used for quality control of the nOPV2 batches.


Introduction
OPV is a live attenuated vaccine that protects well against poliomyelitis by stimulating both humoral and mucosal immunity. However, OPV is composed of three poliovirus Sabin strains that are genetically unstable and can revert to neurovirulence during replication in both cell cultures and vaccine recipients, leading to vaccine-associated paralytic of California, San Francisco, CA, USA. These plasmids were used for the recovery the nOPV2c1 mutant viruses as described below.

Recovery of nOPV2-c1 Mutant Variants
nOPV2-c1 mutant viruses were recovered from the plasmids containing T7 promotor and genomes of nOPV2 VP1-I 143 T, VP1-N 171 D, and VP1-E 295 K mutant variants as described previously [5,22] with some changes. Briefly, plasmid clones of nOPV2-c1 mutants were digested with Hind III, purified with SPRI beads (1.8×), and eluted in 45 µL H 2 O. RNA transcripts were generated from the linearized plasmids using the MEGAscript™ T7 Transcription Kit (Thermo Fisher Scientific, South San Francisco, CA, USA). To remove unincorporated nucleotides and most proteins, the nOPV2-c1 -mutant transcripts were precipitated with LiCl and the pellet was washed once with 1 mL of 70% ethanol and re-centrifuged to maximize removal of unincorporated nucleotides. Then, RNA was dissolved in sterile diethylpyrocarbonate-treated water and stored at −80 • C for downstream use. The purified RNA transcripts were analyzed by electrophoresis on 1% agarose. Then, 2-3 µg of RNA transcript was added to 10 6 HEp-2C cells and subjected to electroporation. The transfected cells were planted in 6-well plates and incubated at 34 • C and 5% CO 2 until the appearance of a complete cytopathic effect (CPE). After the appearance of the CPE, the plate was freeze-thawed three times, and the supernatants were used to transfect a monolayer of HEp-2C cells into 25 cm 2 flasks. The cells were incubated until the appearance of the CPE, and the supernatants were prepared as mentioned above, aliquoted, and stored at −80 • C.
The recovered nOPV2-c1 mutant viruses were sequenced with the Illumina sequencing platform as described below to confirm their nucleotide sequences.

Primers and TaqMan Oligo Probes Used for nOPV2-c1 Mutant Variant Detection and Quantification
The forward and reverse primers for each of the nOPV2 VP1-I 143 T (corresponding to U 2970 C nucleotide change), VP1-N 171 D (corresponding to A 3053 G nucleotide change), and VP1-E 295 K (corresponding to G 3425 A nucleotide change) mutant variants were selected to flank the targeted mutations (Table 1). PCR amplification with these primers resulted in DNA fragments of 74, 60, and 61 nucleotides long for nOPV2 VP1-143T, VP1-171D, and VP1-295K mutants, respectively. The TaqMan oligoprobes were designed to be short (less than or equal to 13 nucleotides), with a conjugated minor groove binder (MGB) and low melting temperature, to allow discrimination of a single point mutation ( Table 1). The primers were synthesized by Integrated DNA Technologies IDT (Durham, NC, USA) and TaqMan probes were synthesized by Thermo Fisher Scientific (South San Francisco, CA, USA). Table 1. Primers and probes used in qmosRT-PCR assay for quantification of the nOPV2-c1 VP1-143T, VP1-171D, and VP1-295K mutants.

qmosRT-PCR Amplification
Viral RNA was extracted from nOPV2 viruses using a QIAamp viral RNA mini kit (QIAGEN, Chatsworth, CA, USA) according to the manufacturer's protocol. The qmosRT-PCR reactions were prepared with the QuantiFast Multiplex RT-PCR Kit (QIAGEN, Valencia, CA, USA) in a final volume of 25 µL using 2 µL of viral RNA. Mutant and nonmutant TaqMan probes were used at final concentrations of 25 nM each in a mixture with the forward and reverse primers ( Table 1) at concentrations of 0.8 µM each. The RNA sample of the nOPV2-c1 virus was used as positive control, and water was used as negative control. RNAs and/or plasmids containing the genomes of nOPV2-c1 mutant and non-mutant viruses were used as reference standards: reference standards, mutant, and non-mutant plasmids/RNAs with known genome copy (GC) numbers per mL were used for extrapolation of GC numbers for mutants and non-mutants of the nOPV2 test samples.
The positive control and reference standard samples were run in duplicate, and the negative control and test samples were run in three repeats. The qmosRT-PCR procedure was performed using a real-time PCR System ViiA7 (Thermo Fisher Scientific, South San Francisco, CA, USA) at the following thermocycler conditions: one cycle of incubation for 20 min at 50 • C and 5 min at 95 • C, followed by 40 cycles each consisting of 15 s at 95 • C, 15 s at 50 • C, and 30 s at 60 • C.

Quantitative RT-PCR for Virus Genome Copy Number Determination
To quantify the genome copy (GC) number in each sample, a qmosRT-PCR was used. Briefly, the qmosRT-PCR reactions were prepared with a QuantiFast Multiplex RT-PCR Kit (QIAGEN, Valencia, CA, USA) in 96-well optical plates in a final volume of 25 µL using 2 µL of RNA of test and control samples and 2 µL of DNA-plasmid for standard references. The RNAs of nOPV2-c1 virus were used as positive controls, and water was used as negative control. Plasmid containing the genome of nOPV2-c1 mutant virus with a known GC number was used as a standard reference for extrapolation of the GC for nOPV2 mutants as test samples, and plasmid containing the genome of nOPV2-c1 non-mutant virus with a known GC number was used as standard reference for extrapolation of the GC number for nOPV2-c1 non-mutants as test samples.
All control and standard reference samples were run in duplicate, and the test samples were run in three repeats. The specific primer pairs and probes used for each virus are presented Table 1. The qmosRT-PCR procedure was performed as described above.

Deep Sequencing Analysis
RNA was extracted from nOPV2 viruses as described above and used for RNA library preparation. The RNA library was prepared using the TruSeq NEBNext Ultra II RNA Library Prep Kit for Illumina (New England BioLabs, Ipswich, MA, USA). Briefly, 100 ng of RNA was fragmented to generate a mean fragment distribution of 500 nt and priming was Vaccines 2022, 10, 1394 5 of 15 performed in one reaction using the buffer provided in the kit. The first and second strands of DNA were synthesized according to the manufacturer's protocol. The resulting DNA fragments were subjected to end repair and ligated to Illumina paired end adaptors. The ligated products were size selected using AMPure XP beads, then amplified using eight cycles of PCR with multiplex indexed primers and purified by magnetic beads (Agencourt AMPure PCR purification system, Beckman Coulter, Brea, CA, USA). Then, the DNA libraries were analyzed for the size and quality with a 4200 Tape Station system (Agilent Technologies, Santa Clara, CA, USA) using a high-sensitivity D1000 Screen Tape. After obtaining the Qubit concentration for the libraries and the mean peak sizes from the Agilent Tape Station profile, paired-end sequencing was performed using a MiSeq System (Illumina, San Diego, CA, USA), producing 250 nt reads. The raw sequencing reads were analyzed with the specialized High-Performance Integrated Virtual Environment (HIVE) platform developed in-house [23]. The RNA sequences of nOPV2-c1 and nOPV2-c2 with GenBank accession numbers MZ245455 and MN654096, respectively, were used as genome references for HIVE alignment and mutations profiling.

Assay Design
The primers were selected based on the sequence of nOPV2-c1. They flanked the targeted mutations and produced 60-74 bp long amplicons depending on the mutation position (Table 1).
Short TaqMan probes were designed with the targeted mutation in the middle to discriminate between mutants and non-mutants with a single point mutation, as we showed previously that mutation in the middle of the oligonucleotide destabilizes the complex between the oligonucleotide and the template DNA [15,24]. The TaqMan probes include a minor groove binder (MGB) moiety at the 3 end that increases the melting temperature (Tm) of the probe and stabilizes probe/target complexes, as well as a nonfluorescent quencher (NFQ) to absorb energy from the fluorescent dye label at the other end of the probe. FAM and VIC dyes are banded to the 5 end of the probes to discriminate between mutants and non-mutants by dye, and there is no interference between these dyes.
To assess the content of the mutants, a 100% mutant virus and 100% non-mutant virus were recovered from plasmids as described above. The plasmids and/or extracted RNAs from these viruses were used in each PCR run to generate standard curves for mutants and non-mutants ( Figure 1A). The quantities of mutants and non-mutants extrapolated from the standard curves were used to calculate the percentages of mutants in each sample.
The baseline and the Ct-threshold could be set separately for mutants and non-mutants; the baseline was adjusted as described by the Thermo Fisher guidelines (http://surf.ed. ac.uk/wp-content/uploads/2014/02/Setting-baselines-and-thresholds-.pdf accessed on 21 August 2022), and the Ct-threshold was set wisely above the noise and higher than any cross-amplification ( Figure 1B). With these settings, the method could preserve its sensitivity and was inherently specific (as the Ct-threshold was set higher than any crossamplification between mutant and non-mutant reference standards).

Evaluation of the Sensitivity and Linearity and Comparison between Mutant Variant Quantifications Using qmosRT-PCR and NGS
To evaluate the linearity and sensitivity of the quantification of the VP1-N 171 D mutants in nOPV2 batches with the qmosRT-PCR assay using plasmids as reference standards, we spiked plasmid that contained the mutation with non-mutant plasmid based on their genome copy (GC) numbers. The resulting percentages of mutants were 50, 25, 12.5, 6.3, 3.1, 1.6, 0.8, 0.4, 0.2, and 0.1. The spiked samples were subjected to the qmosRT-PCR assay and Illumina sequencing (NGS) as described above, and the results for the mutant percentages generated by both methods are presented in Table 2. Both methods were able to detect less than 1% of mutants, although the PCR was more accurate for spike fractions of 3% or higher ( Table 2). An excellent correlation was noted between the expected results and the results of the NGS and qmosRT-PCR methods, with R 2 = 1.00 (Figure 2A,B).
Vaccines 2022, 10, x FOR PEER REVIEW 6 of 15 preserve its sensitivity and was inherently specific (as the Ct-threshold was set higher than any cross-amplification between mutant and non-mutant reference standards).

Evaluation of the Sensitivity and Linearity and Comparison between Mutant Variant Quantifications Using qmosRT-PCR and NGS
To evaluate the linearity and sensitivity of the quantification of the VP1-N171D mutants in nOPV2 batches with the qmosRT-PCR assay using plasmids as reference standards, we spiked plasmid that contained the mutation with non-mutant plasmid based on their genome copy (GC) numbers. The resulting percentages of mutants were 50, 25, 12.5, 6.3, 3.1, 1.6, 0.8, 0.4, 0.2, and 0.1. The spiked samples were subjected to the qmosRT-PCR assay and Illumina sequencing (NGS) as described above, and the results for the mutant percentages generated by both methods are presented in Table 2. Both methods were able to detect less than 1% of mutants, although the PCR was more accurate for spike fractions of 3% or higher ( Table 2). An excellent correlation was noted between the expected results and the results of the NGS and qmosRT-PCR methods, with R 2 = 1.00 (Figure 2A,B).   Similarly, we spiked the mutant virus in the non-mutant virus; first, both viruses were titrated with a conventional CCID 50 assay and the MPBT assay [25] and, as both methods generate similar results, the obtained titers were averaged and used to spike mutant virus in non-mutant virus. As the genome copy (GC) and CCID 50 can differ, the GCs of the spiked viruses and the nOPV2 RNA references were determined by subjecting their stocks to quantitative RT-PCR using plasmids with known GCs as reference standards, as described in Section 2.5 above. The percentages of the spiked samples were recalculated based on the GCs of the virus stocks and they are presented in the first columns in Tables 3-5 (and  in Supplement Tables S1-S3). Two sets of the same spiked samples were prepared and subjected to qmosRT-PCR assays using both plasmids and RNAs of mutants and nonmutants as reference standards. The results of the qmosRT-PCR assays using plasmids as reference standards are presented in Tables 3-6  qmosRT-PCR using RNAs as reference standards are presented as supplements in Tables  S1-S4 and Figures S1-  Similarly, we spiked the mutant virus in the non-mutant virus; first, both viruses were titrated with a conventional CCID50 assay and the MPBT assay [25] and, as both methods generate similar results, the obtained titers were averaged and used to spike mutant virus in non-mutant virus. As the genome copy (GC) and CCID50 can differ, the GCs of the spiked viruses and the nOPV2 RNA references were determined by subjecting their stocks to quantitative RT-PCR using plasmids with known GCs as reference standards, as described in Section 2.5 above. The percentages of the spiked samples were recalculated based on the GCs of the virus stocks and they are presented in the first columns in Tables 3-5 (and in Supplement Tables S3-S5). Two sets of the same spiked samples were prepared and subjected to qmosRT-PCR assays using both plasmids and RNAs of mutants and non-mutants as reference standards. The results of the qmosRT-PCR assays using plasmids as reference standards are presented in Tables 3-6 and in  Figures 3-5, and the results of the qmosRT-PCR using RNAs as reference standards are presented as supplements in Tables S3-S6 and Figures S3-S5. To quantify the levels of the VP1-N171D mutants, the RNA was extracted from the spiked samples and subjected to NGS and qmosRT-PCR assays. The generated results are presented in Table 3 (and in Supplement Table S1). These results showed that NGS generated similar results as the expected percentages prepared by spiking based on GCs. The qmosRT-PCR method was able to consistently detect mutations at levels above 1% (GC). The expected percentages of mutants based on GCs and those determined by NGS correlated very well with the percentages determined with the qmosRT-PCR method using plasmid standard curves, with R 2 = 0.99 for both approaches (as shown in Figure  3A  To quantify the levels of the VP1-N 171 D mutants, the RNA was extracted from the spiked samples and subjected to NGS and qmosRT-PCR assays. The generated results are presented in Table 3 (and in Supplement Table S1). These results showed that NGS generated similar results as the expected percentages prepared by spiking based on GCs. The qmosRT-PCR method was able to consistently detect mutations at levels above 1% (GC). The expected percentages of mutants based on GCs and those determined by NGS correlated very well with the percentages determined with the qmosRT-PCR method using plasmid standard curves, with R 2 = 0.99 for both approaches (as shown in Figure 3A,B; see also the corresponding Supplementary Figure S1A,B, where R 2 = 0.99).     (Table 3) are presented with different colors.
For the quantification of the VP1-E295K mutants in nOPV2-c1, we spiked the mutant virus in non-mutant virus using CCID50 values; both viruses were titrated as mentioned above and used to spike mutants in non-mutants. Two sets of the same spiked samples were prepared on different days. As for the N171D mutants, fractions were also calculated based on the GC concentrations of the two stocks. Each spiked set of samples was Note: GC #, genome copy number; NGS, next-generation sequencing (Illumina sequencing); SD, standard deviation; ND, not detected.  (Table 4) are presented with different colors.
To evaluate the quantification of VP1-I143T mutants in nOPV2-c1, we spiked the mutant virus into the non-mutant virus in a similar way as described above for the VP1-N171D mutants. The titrated viruses were used to spike mutants into non-mutants. As mentioned above, GC concentration fractions were calculated for subsequent use. Two sets of spiked samples with the same percentages of mutants were prepared on different days by the same operator and subjected to three runs of the qmosRT-PCR assay on different days; in total, six runs of qmosRT-PCR were performed for the two sets. Each sample was analyzed in three replicates in each run. The two sets of the spiked samples were also subjected to Illumina sequencing to confirm the percentages of mutants in the spiked samples and to compare them to the result of the qmosRT-PCR. The results of the qmosRT-PCR and Illumina sequencing are summarized in Table 5 (Table 4) are presented with different colors.
Vaccines 2022, 10, x FOR PEER REVIEW 10 of 15 Table S3). The qmosRT-PCR results include the mutants' percentages calculated based on GC units. The qmosRT-PCR assay generated consistent results and was able to detect less than 1% of mutants in the spiked samples. Further, the expected percentages of mutants (GCs) and those determined by NGS correlated well with those determined with the qmosRT-PCR method with R 2 = 0.97, as shown in Figure 5A,B (and in Figure S3 A,B, where R 2 = 0.95).   (Table 5) are presented with different colors.

Evaluation of the Level of Mutant Variants in nOPV2 Vaccine Lots
To evaluate the assays for the quantification of VP1-I143T, VP1-N171D, and VP1-E295K mutants in production stocks of nOPV2, different batches of the nOPV2 candidate 1 (nOPV2-c1) and candidate 2 (nOPV2-c2) vaccines with varying mutant percentages (manufactured by Bio Farma, Indonesia), as mentioned above in the Materials and Methods section, were tested by qmosRT-PCR and Illumina sequencing (NGS). Even though these assays were designed to specifically quantify the mutants in nOPV2-c1, nOPV2-c2 lots were tested in this experiment. The alignment of the primer and probe sequences to the nOPV2-c2 genome revealed the presence of one mismatch in the reverse For the quantification of the VP1-E 295 K mutants in nOPV2-c1, we spiked the mutant virus in non-mutant virus using CCID 50 values; both viruses were titrated as mentioned above and used to spike mutants in non-mutants. Two sets of the same spiked samples were prepared on different days. As for the N171D mutants, fractions were also calculated based on the GC concentrations of the two stocks. Each spiked set of samples was subjected to three runs of the qmosRT-PCR assay on different days by the same operator. Each sample was analyzed in three replicates in each run. The two spiked sets were also subjected to Illumina sequencing as described above to confirm the percentages of the spiked samples and the results of the qmosRT-PCR. The results of the qmosRT-PCR and Illumina sequencing are summarized in Table 4 (and in the correspondent Supplementary  Table S2). The qmosRT-PCR assay generated consistent results and was able to detect less than 1% of mutants in the spiked samples. Quantitatively similar results were obtained for spikes above 2% (GC). The expected percentages of mutants (GC) and those determined by NGS correlated well with percentages determined with the qmosRT-PCR method using GC units, with R 2 = 1.00 and 0.99, respectively, as shown in Figure 4A,B (and in Supplementary  Figure S2A,B, where R 2 = 1:00).
To evaluate the quantification of VP1-I 143 T mutants in nOPV2-c1, we spiked the mutant virus into the non-mutant virus in a similar way as described above for the VP1-N 171 D mutants. The titrated viruses were used to spike mutants into non-mutants. As mentioned above, GC concentration fractions were calculated for subsequent use. Two sets of spiked samples with the same percentages of mutants were prepared on different days by the same operator and subjected to three runs of the qmosRT-PCR assay on different days; in total, six runs of qmosRT-PCR were performed for the two sets. Each sample was analyzed in three replicates in each run. The two sets of the spiked samples were also subjected to Illumina sequencing to confirm the percentages of mutants in the spiked samples and to compare them to the result of the qmosRT-PCR. The results of the qmosRT-PCR and Illumina sequencing are summarized in Table 5 (and in Supplementary Table S3). The qmosRT-PCR results include the mutants' percentages calculated based on GC units. The qmosRT-PCR assay generated consistent results and was able to detect less than 1% of mutants in the spiked samples. Further, the expected percentages of mutants (GCs) and those determined by NGS correlated well with those determined with the qmosRT-PCR method with R 2 = 0.97, as shown in Figure 5A,B (and in Figure S3A,B, where R 2 = 0.95).

Evaluation of the Level of Mutant Variants in nOPV2 Vaccine Lots
To evaluate the assays for the quantification of VP1-I 143 T, VP1-N 171 D, and VP1-E 295 K mutants in production stocks of nOPV2, different batches of the nOPV2 candidate 1 (nOPV2-c1) and candidate 2 (nOPV2-c2) vaccines with varying mutant percentages (manufactured by Bio Farma, Indonesia), as mentioned above in the Materials and Methods section, were tested by qmosRT-PCR and Illumina sequencing (NGS). Even though these assays were designed to specifically quantify the mutants in nOPV2-c1, nOPV2-c2 lots were tested in this experiment. The alignment of the primer and probe sequences to the nOPV2-c2 genome revealed the presence of one mismatch in the reverse primer (Pr143s-3011R) used for quantification of VP1-I143T mutants and two mismatches in the forward primer (nOPV2c1_3021F) used for quantification of VP1-N171D mutants.
The NGS and qmosRT-PCR results using plasmids as reference standards are summarized in Table 6, and the results of the NGS and qmosRT-PCR using RNA as reference standards are summarized in Supplementary Table S4. The results in Table 6 show that the levels of VP1-I 143 T mutants were 0.00-1.93% and 0.00-1.60%, as detected by NGS and qmosRT-PCR, respectively. The levels of VP1-N 171 D mutants were 0.00-28.83% and 0.00-29.17%, as detected by NGS and qmosRT-PCR, respectively. The levels of VP1-E 295 K mutants were 2.07-45.83% and 0.08-62.49, as detected by NGS and qmosRT-PCR, respectively. In the nOPV2 lots tested, the nOPV2_295 mutants were detected at higher levels in comparison to the other mutants, followed by the nOPV2-171 mutants and then the VP1-I 143 T mutants, which were detected at a very low level by both PCR and NGS.
Furthermore, the results (Table 6) show that there was a good correlation between the levels of VP1-N 171 D and VP1-E 295 K mutants determined by NGS and by the qmosRT-PCR assays, with R 2 ≥ 0.97. The coefficient of correlation R 2 observed between the levels of VP1-I 143 T mutants detected by NGS and by qmosRT-PCR assay was 0.85, and this was mostly due to the low frequency of VP1-I 143 T mutants. Similar results were generated by the qmosRT-PCR assay using viral RNAs as reference standards (Supplementary Table S4). The presence of the mismatches in the nOPV2 genome primers did not impact the quantification of the mutants in the nOPV2-c2 lots.
These results show that the qmosRT-PCR assay can be used to quantify levels of nOPV2 mutants as a part of quality control of nOPV2 lots.

Discussion
Virus strains used in the manufacture of viral vaccines undergo mutational changes, which may cause reversion of attenuated strains to virulence, as well as changes in the antigenic properties that may affect the immunogenicity of vaccines [26][27][28]. Therefore, monitoring of the genetic consistency of vaccines is an important task to ensure their quality, safety, and potency.
Sabin strains used in the manufacture of OPV vaccine represent a case study of reversion of attenuated vaccine strains to virulence. Molecular studies performed in the 1980s showed that U 472 C mutation in Sabin type 3 virus [29], A 481 G mutation in Sabin type 2 virus [30], and G 480 A and U 525 C mutations in Sabin type 1 virus [31] accumulate during virus growth in vitro and in vivo and are responsible for reversion to virulence. Similar mutations are found in vaccine-derived strains isolated from cases of vaccine-associated paralytic polio and, therefore, the presence of these revertants in batches of polio vaccine is considered a safety risk. A method based on PCR and restriction enzyme cleavage (MAPREC) was developed to identify batches of oral polio vaccine (OPV) with unacceptable levels of neurovirulent revertants [32,33] and recommended by the WHO as an in vitro method of choice for lot release of OPV. However, MAPREC suffers from some shortcomings, including various technical challenges at multiple steps of its complex protocol.
Previously, we developed a method for quantification of mutants in Sabin 3 virus based on allele-specific primers that consist of three segments: the mutant-discriminating 3 -end, linked by a flexible polyinosine stretch to a specific segment that serves as an anchor to increase stability of binding and prevent non-specific priming [34]. Such a method needs the design of tethered allele-specific primers, extensive screening for the working primers, and optimization before it starts to work properly.
nOPV2 was developed as a genetically stable vaccine to control outbreaks caused by cVDPV2. A recent study [41] comparing Sabin 2 virus and nOPV2 demonstrated significantly lower neurovirulence, as determined in a transgenic mouse assay, in the virus shed by recipients of nOPV2 compared to recipients of conventional OPV2. Furthermore, while shed Sabin 2 virus had the anticipated A481G reversion in the primary attenuation site in domain V in the 5 UTR, associated with increased mouse neurovirulence, the stabilized domain V in the nOPV2 viruses did not show polymorphisms consistent with reversion to neurovirulence [41].
We previously reported on the use of NGS for quality control of nOPV2, with quantitative specifications imposed for the three mutations of interest studied in this work [21]. While NGS may be a viable approach for lot release, maintaining NGS assays and equipment in a validated state for indefinite routine release is, to our knowledge, unprecedented. As such, we sought to explore alternatives.
Quantitative PCR and NGS are broadly used to detect mutant variants in viruses and viral vaccines. Quantitative PCR has the advantage of being easy, faster, and cheaper to use to evaluate up to about four SNPs. NGS, in contrast, allows simultaneous analysis of many genomic loci; it is, however, more technically demanding and expensive.
The method proposed in this communication is based on the use of two short minor groove binder probes; one specific for the mutants and the second specific for the nonmutants. Both are run in the same quantitative RT-PCR reaction with forward and reverse primers (Table 1), as described above. The results of the qmosRT-PCR, used to determine quantities of VP1 I 143 T, N 171 D, and E 295 K mutant variants that may affect the safety and immunogenicity of nOPV2, demonstrated a very good correlation with the results of the Illumina sequencing. This suggests that the qmosRT-PCR assay can be used for quantification of mutant variants in nOPV2 batches. Even though the method was designed to specifically quantify the mutants in the nOPV2-c1, the results for the quantification of the mutants in the nOPV2-c2 were comparable to those generated by NGS.
The method enables the quantification of both mutants and non-mutants in the same sample and in the same reaction, reducing time and labor and simplifying the method.
The proposed method is simple, sensitive, and could be used to test about 18 samples simultaneously, along with appropriate controls and reference samples. In addition, the use of different dyes to discriminate between mutants and non-mutants makes this method much easier to perform. The qmosRT-PCR method produces consistent results (with less than twofold difference from run to run) and is able to detect about 1% of mutant variants.
The quantification of the mutant variants in nOPV2 lots correlated well with the results of the NGS, indicating that the qmosRT-PCR method could be used for quality control of the nOPV2 vaccine.
In conclusion, the quantitative RT-PCR for quantification of undesirable mutants in the nOPV2 vaccine described in this communication offers a simple and rapid method for quantification of mutants in nOPV2 virus. The assays for quantitation of nOPV2 VP1-I 143 T, VP1-N 171 D, and VP1-E 295 K mutants were designed specifically to be applied for quality control during the manufacture of the nOPV2 vaccine.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/vaccines10091394/s1, Figure S1: Results of qmosRT-PCR assay using viral RNA reference standards for quantification of nOPV2-c1_171 mutants generated from the virus-spiked samples: (A) the mutant percentages generated by qmosRT-PCR are plotted against the expected percentages of mutants and (B) the mutant percentages generated by qmosRT-PCR are plotted against the results of Illumina sequencing (NGS). The correlation line is drawn for the averages of the qmosRT-PCR results. Figure S2: Results of qmosRT-PCR assay using viral RNA reference standards for quantification of nOPV2-c1_295 mutants generated from the virus-spiked samples: (A) the mutant percentages generated by qmosRT-PCR are plotted against the expected percentages of mutants and (B) the mutant percentages generated by qmosRT-PCR are plotted against the results of Illumina sequencing (NGS). The correlation line is drawn for the averages of the qmosRT-PCR results. Figure S3: Results of qmosRT-PCR assay using viral RNA reference standards for quantification of nOPV2-c1_143 mutants generated from the virus-spiked samples: (A) the mutant percentages generated by qmosRT-PCR are plotted against the expected percentages of mutants and (B) the mutant percentages generated by qmosRT-PCR are plotted against the results of Illumina sequencing (NGS). The correlation line is drawn for the averages of the qmosRT-PCR results. Table  S1: Quantification of nOPV2-c1_171 mutants in the virus-spiked samples by qmosRT-PCR assay using viral RNA as reference standards. Table S2: Quantification of nOPV2-c1_295 mutants in the virus-spiked samples by qmosRT-PCR assay using viral RNAs as reference standards. Table S3: Quantification of nOPV2-c1_143 mutants in the virus-spiked samples by QmosRT-PCR assay using viral RNAs as reference standards. TableS4: Quantification of nOPV2 mutants in nOPV2 batches by qmosRT-PCR assay using viral RNA as reference standards. [CrossRef]