Background {#Sec1}
==========

Poliomyelitis is a highly contagious neurological disease caused by three distinct serotypes of poliovirus. There are two excellent vaccines that protect against poliomyelitis: inactivated poliovirus vaccine (IPV) and live oral poliovirus vaccine (OPV). Both vaccines played key roles in eliminating paralytic poliomyelitis from most countries worldwide.

The World Health Organization's (WHO) Global Polio Eradication Initiative (GPEI) has made significant progress toward eradicating the disease. It has stopped transmission of wild polioviruses types 2 and 3 \[[@CR1]\] and eliminated wild type-1 poliovirus except for endemic regions in Afghanistan and Pakistan. This remarkable progress was achieved using multiple doses of trivalent live attenuated oral poliovirus vaccine (OPV), which stimulates robust systemic and mucosal immunity \[[@CR2]\] and serves to immunize close contacts as well as vaccinees themselves. All Recipients of OPV excrete mutant variants of Sabin strains in stools; those vaccine-derived polioviruses can evolve into pathogenic variants that circulate in the population (cVDPV) and infect unvaccinated people \[[@CR3]\]. Unfortunately, about one in every million OPV vaccinees and their contacts have developed vaccine-associated paralytic poliomyelitis (VAPP), an unacceptable complication in countries from which the wild-type viruses have been eradicated. In addition, persons with certain kinds of B-cell immunodeficiency can become chronically infected and excrete vaccine-derived virus for a long time \[[@CR4]--[@CR7]\]. Consequently, many countries stopped vaccinating with OPV. To address this issue, in April 2016 the WHO recommended gradually phasing out use of trivalent OPV and replacing it with a bivalent OPV (bOPV) containing only serotypes 1 and 3, plus administering at least one dose of IPV \[[@CR8], [@CR9]\]. An important limitation of IPV is that it elicits poor intestinal immunity. Therefore, children receiving bOPV with only limited number doses of IPV may not acquire sufficient mucosal immunity to prevent infections with type-2 poliovirus, putting them at risk for paralytic disease and amplifying transmission of the type 2 of poliovirus, should it be reintroduced into the environment \[[@CR10]--[@CR12]\]. Efforts are now underway to improve IPV's ability to induce mucosal immunity \[[@CR13]\] but those preparations of IPV are not yet available. Better methods to evaluate mucosal immunity will be increasingly important to evaluate improved vaccines. A straightforward way to evaluate intestinal immunity might be to challenge IPV recipients with OPV and then identify and titrate the viruses they excrete in stool. Surveillance of polioviruses in patients with acute flaccid paralysis and in sewage samples are also important parts of any polio eradication program. In addition, industry needs better tests to identify, titrate and estimate dynamic inactivation of polioviruses during vaccine production. Any new test for polioviruses should be rapid, robust, have high throughput, and offer multiplex simultaneous detection and accurate titration for each of the three serotypes of poliovirus.

Currently many viruses are quantified by traditional culture in susceptible cell monolayers, either by plaque assay or by quantal terminal-dilution in multi-well plates. Plaque assays require tedious visual counting of plaques that vary in size, complicating automated reading. Terminal-dilution quantal methods require monitoring of cytopathic effects (CPE), often over a long period \[[@CR14]--[@CR16]\]. In addition, when several virus types are present in a mixture, all other viruses except the one of interest must be neutralized, further complicating the task. Therefore, conventional virus titration techniques are not well suited to screen the many samples generated during clinical vaccine trials, environmental surveys, and quality control of vaccine production \[[@CR17]\].

We previously compared a real-time PCR-based assay (osRT-PCR) with a conventional quantal terminal-dilution method using CPE-based endpoint assay (CCID~50~). Cells were infected with poliovirus followed by PCR quantitation of viral RNA in cell lysates before CPE appeared \[[@CR18]\]; both methods yielded similar results. Based on that observation, we developed a multiplex PCR-based titration (MPBT) assay that facilitates rapid titrations of OPV strains. Hep-2C cells in 96-well plates are exposed to serial dilutions of virus followed by 2 days of incubation. Cell lysates are then subjected to multiplex quantification by qmosRT-PCR \[[@CR19]\]. Endpoints obtained using this PCR-based assay were compared to those from conventional CCID~50~ assays and found to be similar. The MPBT method is simple, rapid, robust, reproducible, sensitive, and suitable for multiplex titration of viruses.

Results {#Sec2}
=======

Comparative titrations of different lots of Sabin OPV strains with MPBT and CCID~50~ assays {#Sec3}
-------------------------------------------------------------------------------------------

Four different lots of Sabin strains were titrated in simple format (only one virus titrated in each run) and in multiplex format (all three Sabin OPV strains were mixed and titrated simultaneously as one sample) using both MPBT and CCID~50~ assays for simple format and MPBT assay alone for the multiplex format. Results are summarized in Table [1](#Tab1){ref-type="table"}. The titers calculated for each virus were similar for both assays and for both monotype and multiplex formats. The four lots of Sabin OPV strains analyzed in this experiment had high titers, around 8 log~10~ CCID~50~/ml, and similar titers were obtained by MPBT and CCID~50~ assay for samples with titers less than 100 CCID~50~ /ml (see Table [2](#Tab2){ref-type="table"}). These results show that the MPBT assay generated results similar to those of the CCID~50~ assay and worked well for multiplex titration of OPV viruses. Table 1Titration of different lots of Sabin strains with both MPBT and CCID50 assaysSabin type (lot)Single titration^a^Multiplex titration^a^MeanSDRSD (%)MPBT assayCCID50 assayMPBT assay1 (1)8.758.759.138.880.2253.922 (1)8.638.638.638.630.000.003 (1)9.369.2599.200.1844.471 (2)8.68.638.628.620.023.522 (2)8.198.048.188.140.0819.493 (2)8.578.428.498.490.0817.411 (3)8.528.438.568.500.0715.422 (3)8.317.978.268.180.1844.233 (3)8.438.398.348.390.0510.411 (4)7.777.778.077.870.1741.522 (4)7.737.658.147.840.2666.523 (4)8.17.947.847.960.1330.90^a^The titers are in Log10CCID50/ml, *RSD* Relative standard deviation, *SD* Standard deviation of log10 titers, Single titration; Each Sabin Strain was titrated separately, Multiplex titration; All 3 Sabin strains were titrated in the same reaction Table 2Determination of the low limit of titration (LOT) of OPV viruses by MPBT assay and its comparison with LOT of CCID50 assayAssaySabin strain combinationsTiter of the analyzed viruses (CCID~50~/ml)100502510510.1MPBT1125.8981.2835.4831.62NT19.95UDCCID50158.4987.1050.1222.3915.8517.78UDMPBT2123.0351.2928.1822.3917.78UDUDCCID5095.5036.3135.4819.9517.78UDUDMPBT3104.7056.2330.9015.8519.95UDUDCCID5095.5051.3033.8816.9815.85UDUDMPBT1 2 3104.7139.8136.3117.7815.85UDUD74.1374.1321.8819.9517.78UDUD56.2328.1821.8819.9515.85UDUD*NT* Not tested, *UD* Undetermined

MPBT specificity, sensitivity, and ability to determine amounts of each Sabin OPV strain in a mixture {#Sec4}
-----------------------------------------------------------------------------------------------------

Previously we characterized qmosRT-PCR and generated standard calibration curves by testing Sabin viruses of known titers (expressed as CCID~50~/ml) \[[@CR19]\]. RNA was extracted from the three OPV viruses and serial ten-fold dilutions prepared from individual virus RNAs, combined RNA from two viruses, and combined RNA from all three viruses; samples were subjected to quantitative simplex one-step RT-PCR, duplex one-step RT-PCR, or triplex one-step RT-PCR, depending on the combinations of RNAs tested in the same reaction to generate standard curves. All curves showed good linearity with R-squared values exceeding 0.95. The linear ranges were 9 log~10~ for single-type PCR, 8--9 log~10~ for duplex PCR and 7--8 log~10~ for triplex PCR. These results showed that both monospecific and multiplex PCRs were very specific and sensitive. The limit of quantification (Based on viral RNA quantification) of qmosRT-PCR for three Sabin OPV strains mixed together fell between 0.29--2.86, 0.13--1.26 and 0.36--3.60 CCID~50~/ml for types 1, 2, and 3 respectively \[[@CR19]\].

In this work, virus dilutions containing 0.1 to 100 CCID~50~/ml were used to determine the sensitivity of the MPBT assay. For single-virus titrations, we compared MPBT and CCID~50~ assays. Results, summarized in Table [2](#Tab2){ref-type="table"}, showed that the limit of titration (LOT) of single-virus titrations were 0.1 to 1 CCID~50~/ml for Sabin 1 and 1 to 5 CCID~50~/ml for Sabin 2 and 3 for both MPBT and conventional CCID~50~ assays. When all three Sabin strains were titrated together in the same reaction, LOTs of the MPBT assay were 1--5 CCID~50~/ml for Sabin 1, 2, and 3. Both assays had similar sensitivity for titrations of a single virus. While CCID~50~ assays cannot titrate more than one virus per reaction, the MPBT assays were able to titrate each Sabin strain mixed in the same sample with high sensitivity and specificity.

Correlations between MPBT and CCID~50~ assays {#Sec5}
---------------------------------------------

We assessed correlations between MPBT and CCID~50~ assays using three samples (one for each Sabin strain) with titers previously determined by CCID~50~ assay. The viruses with known titers were mixed, serially diluted ten-fold, and each dilution subjected to MPBT assay titration as described above. The results of MPBT assays plotted against the known CCID~50~ titers of the corresponding dilutions for each Sabin strain are presented in Fig. [1](#Fig1){ref-type="fig"}. Excellent correlations were observed between the two assays for each of the three Sabin strains with R-squared value of 0.99 for Sabin 1 virus and 1.0 for Sabin 2 and 3 viruses. The estimated slopes and associated 95% confidence intervals are 0.95 (0.85, 1.06), 0.95 (0.88, 1.02), and 1.02 (0.99, 1.05) for Sabin 1, 2, and 3 viruses, respectively, suggesting no proportional bias between the two methods. Fig. 1Evaluation of correlations between MPBT and CCID~50~ assays using three Sabin strains with known titers (determined by CCID~50~ assay)

Consistency of MPBT assay {#Sec6}
-------------------------

To assess consistency of the MPBT assay, ten-fold dilutions, we titrated one lot of Sabin 1, 2 and 3 mixture simultaneously five times on the same day. On another day, the same operator titrated the same lot of the three Sabin strains simultaneously five times using two-fold dilutions. Result of these experiments, summarized in Table [3](#Tab3){ref-type="table"}, showed that the MPBT assay generated consistent results for all three Sabin strains: the repeats with ten-fold dilutions and with two-fold dilutions. We observed no significant differences in titers repeatedly obtained with ten-fold and two-fold dilutions. The maximum percent of relative standard of variation observed was 65 for Sabin 2 virus repeatedly titrated with ten-fold dilutions. in short, the MBPT assay yielded consistent and reproducible titers. Table 3Results of consistency evaluation of MPBT assay (titers are expressed on log10CCID50/ml)Sabin strainRepeat codesMeanSDRSD %12345Sabin 1^a^8.758.639.108.758.758.800.1842.69Sabin 1^b^8.668.778.849.078.518.770.2151.03Sabin 2^a^8.638.638.759.258.758.800.2664.91Sabin 2^b^8.288.398.518.518.288.390.1226.96Sabin 3^a^8.888.758.889.009.258.950.1945.56Sabin ^b^8.328.138.438.248.358.290.1126.74*SD* Standard deviation of log10 titers, *RSD* Relative standard deviation^a^Virus was titrated with 10-fold dilution^b^virus was titrated with 2-fold dilution

Robustness of virus titrations by the MPBT assay {#Sec7}
------------------------------------------------

We studied the robustness of the MPBT assay to ensure that measurements of virus titer remained unaffected by small variations deliberately introduced into the procedure. To do this, we performed MPBT assays for the three Sabin strains simultaneously but seeded different numbers of cells into the wells (1--4 × 10^4^). Results of these titrations, summarized in Table [4](#Tab4){ref-type="table"}, showed that, despite four-fold variation in cell numbers per well, MPBT assays generated similar results, with relative standard deviation percentages for each virus of 88 for Sabin 1, 16 for Sabin 2 and 36 for Sabin 3. Table 4Evaluation of the effect of cell numbers per well variation on the MPBT assay results (log10CCID50/ml)Sabin StrainCell numbers per wellMeanSDRSD %1 10^4^2 10^4^4 10^4^18.59.128.628.750.3387.9828.58.388.58.460.0716.0538.628.888.888.790.1535.62*SD* Standard deviation of log10 titers, *RSD* Relative standard deviation

We also investigated effects of delay between diluting virus and adding cells in MPBT assays. Delays of 15 min to 1 h showed no effect on results. Titers, summarized in Table [5](#Tab5){ref-type="table"}, were similar after delays of 15, 30, 45 and 60 min, demonstrating that the MPBT assay was robust when all three Sabin strains were titrated simultaneously in the same reaction. Table 5Evaluation of the effect of delay between diluting virus and adding cells on MPBT assay results (log10CCID50/ml)Sabin StrainDelay in minutesMeanSDRSD %1530456018.758.889.138.888.910.1637.8628.758.638.638.508.630.1023.8439.258.889.139.009.070.1638.15*SD* Standard deviation of log10 titers, *RSD* Relative standard deviation

MPBT assays with blind samples {#Sec8}
------------------------------

To validate MPBT assay, polioviruses of all three serotypes with known CCID~50~ titers were mixed in different combinations and concentrations and tested by a single operator using the MPBT assay in blind format. The results are presented in Table [6](#Tab6){ref-type="table"}. The assay specifically identified each type of Sabin OPV virus in all combinations, and the titers for each Sabin strain were those expected based on the dilution tested. Table 6Virus titer determination in samples prepared with different combination of Sabin strains and analyzed in blind format with MPBT assaySamples codesSpiked CCID50 titer^a^ (Sabin type)Results of MPBT assaySabin typeTiter^a^A7.63 (1)17.5B7.04 (2)27C7.42 (3)37.75D7.63 (1)17.387.04 (2)26.75E7.63 (1)17.387.42 (3)37.5F7.04 (2)26.887.42 (3)37.25G7.63 (1)17.257.04 (2)27.257.42 (3)37.75^a^Titers of Sabin strains are expressed on log10CCID50/ml

MPBT assays of stool samples collected from a clinical trial of OPV2 {#Sec9}
--------------------------------------------------------------------

To determine the ability of the MPBT assay to titrate OPV virus strains in clinical samples, we used supernatants of stool samples collected during a clinical trial of monovalent OPV2, part of the FIDEC study \[[@CR20]\]. Samples were diluted two-fold with DMEM and filtered through 0.22-μm spin-tube filters (COSTAR). The filtrates were analyzed with both MPBT and CCID~50~ assays. Results, summarized in Table [7](#Tab7){ref-type="table"}, show that titers for each sample were similar with both assays. This finding demonstrates that the MPBT assay can be used to detect and titrate OPV virus shed in the stool during clinical trials and in the environment. Table 7Analysis of stool samples collected from clinical trial of OPV2 with MPBT and CCID50 assaysStool sample codesMPBT assay (log~10~CCID~50~/ml)CCID50 assay (log~10~CCID~50~/ml)13.433.5123.773.2933.863.6744.053.952.12.2362.492.8173.553.43

Discussion {#Sec10}
==========

Sensitive, specific, reproducible and high-throughput assays of infectious viruses are needed to detect polioviruses in clinical and environmental samples and essential to monitor the manufacturing of live and inactivated poliovirus vaccines \[[@CR21]\]. The viral plaque assay was developed in 1952 by Renato Dulbecco and became widely used to quantify viruses \[[@CR22]--[@CR25]\]. It is one of the most widely used method in virology to purify a clonal population of virus or to determine viral titer as plaque formation unite (PFU) per milliliter \[[@CR26]--[@CR30]\]. Since most plaques contain the progeny of a single virus particle, clonal populations can be purified by isolating virus from individual plaques. However, it is usually necessary to stain monolayers to enhance the contrast between adjacent living cells and the plaques. Therefore, viral plaque assays are typically time-consuming, generally requiring four to 10 days, depending on the virus. Furthermore, plaque assays work only for viruses capable of infecting cells in adherent monolayers, multiplying and lysing them; some viruses do not form plaques in culture. CCID~50~ assays are also used to titrate viruses that cause CPE in tissue cultures---usually over five to 20 days---and may be effective even when many cells in the culture remain viable. CCID~50~ titrations are also relatively cheap, economical of experimental materials and reagents and other consumable supplies, and easy to analyze \[[@CR27]\]. However, not all viruses cause CPE in cells they infect. Results of both plaque and quantal CCID~50~ assays are somewhat variable; relative errors with the plaque assay can be more than 10%, and CCID~50~ titers often have 35% error rate \[[@CR31]\]. In general, CCID~50~ assays are used more often than plaque assays to titrate viruses for research and diagnosis \[[@CR27], [@CR31]--[@CR33]\], but both are difficult to automate, and neither is especially suitable for large-scale analyses of clinical and environmental samples. Without further modification, neither assay deals easily with samples containing more than one virus, such as trivalent OPV viruses.

Quality control of poliovirus vaccine production, including in-process confirmation of virus identity and study of viral inactivation dynamics, also requires rapid testing of many samples. Manufacturers use seroneutralization identity testing of two serially diluted samples and a reference virus in cell culture medium tested with and without addition of specific monoclonal antibodies. The identity of the serotype is determined from the difference between infectivity titers of samples with and without specific antibodies \[[@CR34], [@CR35]\]. This method is time-consuming, labor-intensive, and error-prone, sometimes yielding false negative results. To overcome these limitations, encouraged by our previous finding that real-time PCR quantitation of virus from lysates of infected cell cultures gave results similar to those with a conventional CCID~50~ assay \[[@CR18]\], we developed a quantitative multiplex one-step RT-PCR (qmosRT-PCR) \[[@CR19]\] and a multiplex PCR-based titration (MPBT) assay for multiplex simultaneous titration of all three infectious Sabin OPV poliovirus strains.

Similar to CCID~50~ assays, the MPBT assay tests serially diluted samples of virus mixed with susceptible cells in replicate wells of 96-well plates incubated at 36 °C, except that instead of waiting for the appearance of CPE, medium is discarded after 36--42 h of incubation, cells lysed with 0.9% triton, and then viral nucleic acids in cell lysates are quantified by qmosRT-PCR assay; results are used to determine virus titers expressed as CCID~50~ units. Therefore, virus titer is not determined based on proportionality to Ct values obtained by qmosRT-PCR, but rather qmosRT-PCR is used as a readout method to determine the presence or absence of viral replication. Therefore qmosRT-PCR results are expressed in yes / no format, allowing a simple Karber formula to be used similar to conventional CCID~50~ assay. As a result this assay detects only live viruses and not nucleic acids that may be present in dead particles.

MPBT proved to be very sensitive, detecting an equivalent of 1--5 CCID~50~/ml of OPV types 2 and 3 alone and OPV types 1, 2 and 3 in mixtures and 0.1--1 CCID~50~/ml of OPV type 1 virus tested alone. The MPBA assay detected OPV types 1, 2 and 3 alone with sensitivities similar to those of the CCID~50~ assay (Table [1](#Tab1){ref-type="table"}). In addition, the MPBT method generated titers similar to those of CCID~50~ assays for all combinations of Sabin OPV poliovirus strains s (Tables [2](#Tab2){ref-type="table"} and [6](#Tab6){ref-type="table"}).

We found excellent correlations between the two assays for all the three Sabin OPV poliovirus strains with an R-squared value of 0.99 for Sabin 1 and R-squared values of 1.0 for each of Sabin 2 and 3 OPV polioviruses (Fig. [1](#Fig1){ref-type="fig"}). The MPBT assay generated consistent results for all three Sabin OPV poliovirus strains (Table [3](#Tab3){ref-type="table"}). The MPBT assay proved to be very robust; even a one-hour delay between diluting the viruses and adding cells and a four-fold variation in numbers of cells added had no effect on titers of viruses (Tables [4](#Tab4){ref-type="table"} and [5](#Tab5){ref-type="table"}). This method reduced the time needed to titrate OPV polioviruses from the seven to 10 days for conventional CCID~50~ to only two to 3 days for MPBT. In contrast to CCID~50~ assay, the MPBT assay accurately titrated all three Sabin OPV poliovirus strains mixed in the same sample.

Conclusions {#Sec11}
===========

The MPBT assay described in this communication offers a simple and rapid alternative to traditional CCID~50~ assays to detect, identify, and titrate either individual or combined serotypes of Sabin OPV polioviruses. The MPBT method is well suited to detect and quantify polioviruses in the many fecal samples collected during clinical trials of new poliovirus vaccines (IPV/OPV), and during routine clinical and environmental poliovirus surveillance programs. The MPBT assay can also be applied during manufacture of poliovirus vaccines, for in-process quality control by identifying serotypes of vaccine polioviruses (potentially replacing the current seroneutralization identity test) and monitoring inactivation dynamics. This assay is suitable to automate, facilitating high-throughput applications, improving consistency of titrations, and saving time, cost, and labor. The MPBT method can also be applied to detect and titrate other viruses, including those that produce no CPE.

Methods {#Sec12}
=======

Viruses, Hep-2C cells and clinical samples {#Sec13}
------------------------------------------

Lots of US neurovirulence poliovirus reference vaccines (containing Sabin type 1, 2, and 3 OPV strains having GenBank accession numbers: AY184219, AY184220, and AY184221 respectively) were used to compare results of CCID~50~ and MPBT assays.

HEp-2C cells (ATCC® CCL-23™), derived from a human carcinoma, were cultured in 175-cm^2^ flasks at + 37 °C ± 2 °C in Dulbecco's modified eagle medium (DMEM; Gibco) supplemented with 5% fetal calf serum (FCS; Gibco) and penicillin-streptomycin (100 U/mL and 100 μg/ mL, respectively; Gibco). Viable cells were sampled, stained with Trypan Blue (Invitrogen) and counted with a hemocytometer (Countess™; ThermoFischer). Hep-2C cells were used for both the CCID~50~ and MPBT assays.

Anonymized stool samples used in this work were collected from a clinical trial of OPV2 ("Fighting Infectious Diseases in Emerging Countries" \[FIDEC\]) \[[@CR20]\] that was approved by ethics committees and as appropriate by National Regulatory Authorities, the Colorado Multiple Institutional Review Board, and the Western Institutional Review Board.

Virus titration by CCID~50~ assay {#Sec14}
---------------------------------

Virus strains used in OPV were quantified by endpoint dilution titration on HEp-2C cells and the results expressed as log~10~ CCID~50~/mL. (CCID~50~ is cell culture 50% infectious dose, defined as that dilution of virus required to infect 50% of the cell monolayers.) Serial ten-fold or two-fold dilutions of viral samples were prepared in 100 μl of DMEM supplemented with 2% FCS in 96-well plates. One-hundred-μl aliquots of cell suspension containing 2 × 10^4^ HEp-2C cells in DMEM with 2% FCS were added to each well of diluted virus in replicate wells of 96-well plates. Virus-infected plates were incubated for 10 days at 36 °C in a 5% CO~2~ humid atmosphere; wells were periodically inspected for CPE. Wells showing CPE were counted on day ten after infection and virus titers calculated using the Spearman-Karber formula \[[@CR36]\].

Multiplex PCR-based titration (MPBT) assay {#Sec15}
------------------------------------------

Viruses were diluted and mixed with cells in 96-well plates as for the CCID~50~ assay described above, with the following exceptions: plates were incubated for only 42 h, after which the medium was discarded and 50-μl aliquots of lysis solution (DMEM with 0.9% Triton X-100) were added to each well. The plates were then sealed with foil and stored at -80 °C prior to quantitative multiplex one-step RT-PCR (qmosRT-PCR) analysis \[[@CR19]\].

The plates with the triton-lysed cells were thawed for 30 min at room temperature and briefly centrifuged in a 5810R centrifuge (Eppendorf, San Diego, CA) at 1000 rpm for 1 min to collect the lysate at the bottom of the wells. Dilutions of 1:10 were then prepared by adding 90 μl of molecular biology grade water (5-PRIME, Gaithersburg, MD) to 10 μl of the cell lysates.

All three poliovirus serotypes were quantified in the same qmosRT-PCR reaction as described previously \[[@CR19]\]. Briefly, TaqMan probes with different dyes (FAM, VIC, and NED fluorescent reporters and a non-fluorescent quencher; Applied Biosystems) were used in the qmosRT-PCR to discriminate between different poliovirus serotypes. Two microliters of diluted lysate were added to each well of a PCR plate, and the multiplex quantitative RT-PCR assay was performed as described \[[@CR19]\]. The samples that have Cts less or equal than 40 are considered positive and samples with Cts higher than 40 are considered negative. The PCR reaction considered non-valid if positive control was negative and/or negative control was positive. Each well of the plate was scored as positive or negative for the presence of the virus activity and used for titers calculation of each OPV serotype according to Spearman-Karber formula \[[@CR36]\]. The Spearman-Karber formula used is as follows: log CCID50 = −(L-d(S-0.5)) where, L is the log of lowest dilution in test, d is the difference between log dilutions and S is the sum of proportions of positive tests (i.e. cell cultures positive for the presence of the virus by PCR).

bOPV

:   Bivalent oral poliovirus vaccine

CCID50

:   Cell culture infectious dose 50%

CPE

:   Cytopathic effect

GPEI

:   Global Polio Eradication Initiative

IPV

:   Inactivated poliovirus vaccine

MPBT

:   Multiplex PCR-based titration

OPV

:   Trivalent oral poliovirus vaccine

PFU

:   Plaque formation unite

qmosRT-PCR

:   Quantitative multiplex one-step reverse-transcriptase polymerase chain reaction

RSD

:   Relative standard deviation

RT-PCR

:   Reverse-transcriptase polymerase chain reaction

WHO

:   World Health Organization

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

We thank Dr. David Asher for critical and careful review of the manuscript, and Dr. Tsai-Lien Lin for statistical advice and data analysis.

HM, ER and TZ performed the experiments. HM, TL and ML: Analysis of results, HM and ML designed the study. ML and KC supervised the work and wrote the article. All the authors have read and approved the final version of this manuscript.

This work was not funded by external institutions or grants.

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Anonymized stool samples used in this work were collected from a clinical trial of OPV2 ("Fighting Infectious Diseases in Emerging Countries" \[FIDEC\]) \[[@CR20]\] that was approved by ethics committees and as appropriate by National Regulatory Authorities, the Colorado Multiple Institutional Review Board, and the Western Institutional Review Board.

Not applicable.

The authors declare that they have no competing interests.
