key: cord-0959190-3azfz0w9
authors: Wang, Yin; Noll, Lance; Porter, Elizabeth; Stoy, Colin; Dong, Junsheng; Anderson, Joe; Fu, Jinping; Pogranichniy, Roman; Woodworth, Jason; Peddireddi, Lalitha; Bai, Jianfa
title: Development of a differential multiplex real-time PCR assay for porcine circovirus type 2 (PCV2) genotypes PCV2a, PCV2b and PCV2d
date: 2020-09-12
journal: J Virol Methods
DOI: 10.1016/j.jviromet.2020.113971
sha: 86eaebbc3df12dd534114279cb4361768f9ea0a4
doc_id: 959190
cord_uid: 3azfz0w9

A multiplex quantitative real-time polymerase chain reaction (mqPCR) assay was developed and validated for detection and differentiation of porcine circovirus type 2 (PCV2) genotypes, PCV2a, PCV2b and PCV2d. Single nucleotide polymorphism in primers or probes was deployed for different genotype detections, while conserved sequence in the 3 end of a primer and in the middle of a probe was used for the targeted genotype. In silico analysis of 2601 PCV2 ORF2 sequences showed that the predicted strain coverage of the assay was 93.4 % (409/438) for PCV2a, 95.1 % (1161/1221) for PCV2b and 93.6 % (882/942) for PCV2d strains. The PCR amplification efficiencies were 94.5 %, 100.2 %, and 99.2 % for PCV2a, PCV2b and PCV2d, respectively, with correlation coefficients >0.995 for all genotypes. The limits of detection (LOD) were 1.58 × 10(−4) TCID50/mL for PCV2a, 5.62 × 10(−4) TCID50/mL for PCV2b, and 3.16 × 10(−3) TCID50/mL for PCV2d. Sanger sequencing of 74 randomly selected PCV2 positive clinical samples confirmed the genotypes of strains identified by the mqPCR. Validation with clinical samples co-positive for target and non-target pathogens demonstrated that the mqPCR assay specifically detected targeted viruses without cross reacting to each other or to other common porcine viruses.

Porcine circovirus type 2 (PCV2) plays a significant role in porcine circovirus associated diseases (PCVAD) and is one of the most economically important porcine viral pathogens (Gillespie et al., 2009) . Since it was first isolated as the causative agent of postweaning multisystemic wasting syndrome in 1998 (Ellis et al., 1998) , PCV2 has been a highly prevalent disease worldwide, including in Asia, America and Europe (Saporiti et al., 2020; Wang et al., 2019a; Zheng et al., 2020) .

With a high substitution rate, 5 major genotypes of PCV2 have been identified based on ORF2 sequences: PCV2a, PCV2b, PCV2c, PCV2d and PCV2e (Davies et al., 2016) . Among the five genotypes, PCV2a, PCV2b and PCV2d are the most commonly circulating genotypes in the USA (Wang et al., 2019b; Xiao et al., 2016) , in China (Hou et al., 2019; Lv et al., 2020) and in European countries (Saporiti et al., 2020) . There are two major genotype shifting events that have occurred globally: PCV2b replaced PCV2a as the dominant genotype around 2003, and PCV2d became the most prevalent genotype beginning in 2012 (Xiao et al., 2016) . The first shift is related to increased severity of clinical PCVAD (Beach and Meng, 2012) , while the second shift may be related to inappropriate vaccine applications (Karuppannan and Opriessnig, 2017) .

ORF2-targeted gene sequencing is widely used as a gold standard for genotyping of PCV2 strains (Wang et al., 2019b) . Other methods are developed as well for the differentiation of PCV2a J o u r n a l P r e -p r o o f and PCV2b strains, such as real time PCR, loop-mediated isothermal amplification method (LAMP) and restriction fragment length polymorphism (RFLP) (Guo et al., 2010; Qiu et al., 2012; Wozniak et al., 2019; Xiao et al., 2016) . However, there's no assay that can simultaneously identify genotypes PCV2a, PCV2b and PCV2d. Therefore, based on the current PCV2 ORF2 sequences available from GenBank database, we have developed a multiplex quantitative real time PCR for rapid detection and differentiation of the three genotypes with high sensitivity, specificity and strain coverage.

The sequence dataset was established from our previous study (Wang et al., 2019b) .

Briefly, all available PCV2 ORF2 sequences were downloaded from the GenBank database (https://www.ncbi.nlm.nih.gov/ genbank/). The sequence alignment was performed with MAFFT (Katoh, Misawa, Kuma, & Miyata, 2002) . Low quality sequences, as indicated by the presence of premature stop codons, and short sequences were removed from the analysis. The sequence identity matrix calculation was performed using BioEdit 7.2.5 (https://bioedit.software.

informer.com/). The phylogenetic analysis was conducted with MEGA 7.0.26 (Kumar, Stecher, and Tamura, 2016) . The best fit model was selected based on Bayesian information criterion (BIC).

The maximum likelihood phylogenetic trees were constructed using the best substitution pattern with the lowest BIC scores. The reliability of clusters formed in the tree was evaluated by performing 500 bootstrapping reiterations.

J o u r n a l P r e -p r o o f Based on the sequences from the dataset, one set of primers and probe was designed for each genotype (PCV2a, PCV2b and PCV2d). The PCV2a, PCV2b and PCV2d probes were labeled with 5′-Texas Red and 3′-BHQ1, 5′-VIC and 3′-BHQ1, and 5′-FAM and 3′-BHQ1, respectively (Table 1) .

Cell culture isolates of PCV2a, PCV2b and PCV2d from the Virology Section of Kansas State Veterinary Diagnostic Laboratory were propagated as previously described (Pogranichniy et al., 2002) , which were genotyped by sequencing. The viral DNA was extracted from 140 μl of clinical samples or cell culture by QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's recommendations, and stored at −80 °C until use. Table 1 .

All PCR reactions were performed in a 20 μL total reaction volume composed of 5 μL of DNA samples prepared as described (Wang, Das, et al., 

Standard curves were generated with triplicates of 10-fold serial dilutions of the cell culture isolates of PCV2a, PCV2b and PCV2d to evaluate the analytical sensitivity of the assay. The median tissue culture infectious dose (TCID50) of the highest dilution that still generated positive Ct values was considered the limit of detection (LOD) for cultured viruses.

The specificity of the assay was first evaluated in silico with the online NCBI primer designing tool, Primer-Blast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/), followed by confirmation with Sanger sequencing of PCV2 ORF2. The conserved primers flanking the ORF2 were designed and listed in Table 1 . Furthermore, the cell culture isolates and clinical samples that were positive to specific swine pathogens were tested to evaluate the specificity of the assay (Table   2) 

J o u r n a l P r e -p r o o f Based on our previous study (Wang et al., 2019b) , PCV2 strains were separated into 12 clusters, namely Cluster 1-12, by phylogenetic analysis; Cluster 8, Cluster 3 and Cluster 1 corresponded to genotypes, PCV2a, PCV2b and PCV2d, respectively. The three genotypes were identified as common strains currently circulating in the USA (Wang et al., 2019b; Xiao et al., 2016 ). An unrooted phylogenetic tree generated with PCV2 ORF2 sequences indicated that the three genotypes grouped into three separate clusters (Figure 1 ). The triangle markers in the tree indicated the ORF2 sequences generated in our lab that were used to validate this genotyping real time PCR assay; they share 99.4-100 % identity with the sequences collected from the GenBank.

More complete analysis with 438 PCV2a strains, 942 PCV2b ORF2 strains and 1221

PCV2d strains has identified one set of primers and probe for each genotype (Illustrated in Figure   2 ). In silico analysis showed the PCR set can perfectly match 75.3% of PCV2a strains, 86.3% of PCV2b strains, and 78.1% of PCV2d strains. Assuming a single nucleotide mismatch to a primer (excluding the more critical 5 bp in the 3' end) or a probe (excluding the more critical 5 bp in the middle) in the assay can still amplify and generate signal (Wang et al., 2020b) , the strain coverages would be increased to 93.4%, 95.1% and 93.6% for the three genotypes, respectively. Information of all primers and probes are shown in Table 1 .

Analytical sensitivity of the mqPCR assay was analyzed using standard curves generated by three replications of 10-fold serial dilutions of the cell culture isolates. The data was presented by plotting the Ct values against log dilution factors. The PCR amplification efficiencies were 94.5% for PCV2a, 100.2% for PCV2b and 99.2% for PCV2d, with correlation coefficients (R 2 ) all greater than 0.995 (Figure 3 ). The LODs were 1.58 x 10 -4 TCID50/ml for PCV2a, 5.62 x 10 -4 TCID50/ml for PCV2b, and 3.16 x 10 -3 TCID50/ml for PCV2d.

The specificity of primers and probes was first tested by in silico analysis using Primer-Blast, which determined sequences were unique to their respective assay targets. Then, assay specificity was evaluated by comparison with Sanger sequencing results. Seventy-four PCV2 positive samples were tested by the genotyping real-time PCR; results showed 8 PCV2a, 6 PCV2b and 60 PCV2d strains. Genotypes of all 74 samples were confirmed by Sanger sequencing (Table   2 ). Specificity was also tested with clinical samples that previously tested positive to non-target pathogens. The results demonstrated that the assay specifically detected positive samples and identified PCV2a, PCV2b or PCV2d genotypes without cross-detecting each other. Furthermore, no positive signals were generated from clinical samples that were positive to PCV3 (n=3), PRRSV-2 (n=6), SIV (n=2), PPIV (n=5), RVA (n=2), RVB (n=1), RVC (n=1), PEDV (n=16), and TGEV (n=2), indicating a good specificity of the assay (Table 3) .

PCV2 is a small DNA virus with a circular genome of 1767-1768 nucleotides. The ORF2 gene (702 nucleotides) encodes the Capsid protein, which dominates immunogenicity (Nawagitgul et al., 2002) .Extensive investigations have been done on genetic variations caused by pointmutations and recombinations (Firth et al., 2009; Franzo et al., 2016a) , association of genotypes and disease severity Opriessnig et al., 2006; Opriessnig et al., 2008) , and viral evolution under the selection pressure of vaccination (Franzo et al., 2016b; Karuppannan and Opriessnig, 2017; Xiao et al., 2016) . Apparently genotyping of PCV2 is an important diagnostic tool for the study of viral pathogenesis and epidemiology of the disease, which will provide key information towards formulating strategies for disease management and vaccine development and applications.

By analyzing whole genome sequences, sequence variations among different genotypes were found mostly in the capsid gene, ORF2, which is consistent with previous studies (Cheung et al., 2007; Olvera, Cortey, and Segales, 2007; Wang et al., 2019b) . In this study, all PCV2 ORF2 sequences were downloaded from the current GenBank database and included in the genetic analysis to ensure high strain coverage. Due to the high rate of nucleotide substitution in the genome, PCV2 displays substantial genetic variations (Firth et al., 2009 ). In our previous study, the inter-cluster and intra-cluster identities were as low as 83.6% and 87.9%, respectively. (Wang et al., 2019b) . Because of this high sequence variation within each genotype and relatively high genome homology between genotypes, the development of a molecular genotyping assay with high strain coverage can be challenging. In our design, the three sets of primers and probes can match 75.3-86.3% strains of each genotype. In general, an assay with single nucleotide mismatches occurring in the middle, especially in the 5′ end of a primer, or in the two ends of a probe may still amplify and generate a signal (Bru, Martin-Laurent, and Philippot, 2008; Forney, Zhou, and Brown, 2004; Wang et al., 2020a) . Therefore, when single nucleotide mismatch in such situations is considered in the design of the primers and probes, strain coverage can potentially increase to 93.4%-95.1%. With that said, this genotyping assay should not replace general detection assays that may have higher strain coverage (Wang et al., 2019a) .

Although some PCV2 genotyping assays have been published, these assays only differentiate between PCV2a and PCV2b strains (Wozniak et al., 2019; Xiao et al., 2016) . For that reason, Sanger sequencing of the ORF2 gene, which is a labor intensive and time-consuming process (Lv et al., 2020; Zheng et al., 2020) , is still widely used for PCV2 genotyping. We have developed this genotyping mqPCR assay that can use a single PCR reaction to differentiate the three major PCV genotypes, (PCV2a, PCV2b and PCV2d) circulating in the US. This simplified J o u r n a l P r e -p r o o f method is a low-cost and faster turnaround protocol, and should facilitate the genotyping process for field PCV2 strains. We have already developed a PCV2 general assay using the conserved ORF1 gene as the target (Wang et al., 2019a) , which has enabled us to have a much higher strain coverage (99.1%) for PCV2 detections. By utilizing that generic test for general detection and surveillance, followed by the genotyping PCR reported here, we should be able to detect and differentiate the majority of field PCV2 strains.

In conclusion, the newly developed and validated mqPCR assay allows for rapid, sensitive and specific detection and differentiation of the most common PCV2 genotypes, (PCV2a, PCV2b and PCV2d) in clinical samples. In silico analysis and clinical sample testing indicated that the assay has relatively high strain coverage. It can be used as an alternative method for PCV2 genotyping.

Author Statement: All authors have read the revision and response to reviewers' comment document, and are agreed with the submission of this revised manuscript.

The authors declare no competing interests. J o u r n a l P r e -p r o o f b: The lower boundary of the range of percentage match is generated from perfect matches in a primer pair and its corresponding probe against total sequences analyzed; the higher boundary of the range is generated by allowing presence of single nucleotide variation in each primer (but keep the 5bp conserved in the 3 end), or in the probe (but keep the 5pb conserved in the middle). 

Phylogenetic characterization of porcine circovirus type 2 in PMWS and PDNS Korean pigs between

Efficacy and future prospects of commercially available and experimental vaccines against porcine circovirus type 2 (PCV2)

Quantification of the Detrimental Effect of a Single Primer-Template Mismatch by Real-Time PCR Using the 16S rRNA Gene as an Example

Detection of two porcine circovirus type 2 genotypic groups in United States swine herds

Diagnostic phylogenetics reveals a new Porcine circovirus 2 cluster

Isolation of circovirus from lesions of pigs with postweaning multisystemic wasting syndrome

Insights into the Evolutionary History of an Emerging Livestock Pathogen: Porcine Circovirus 2

Molecular microbial ecology: land of the one-eyed king

Phylodynamic analysis of porcine circovirus type 2 reveals global waves of emerging genotypes and the circulation of recombinant forms

Porcine circovirus type 2 (PCV2) evolution before and after the vaccination introduction: A large scale epidemiological study

Porcine circovirus type 2 and porcine circovirus-associated disease

Porcine circovirus type 2 (PCV2): genetic variation and newly emerging genotypes in China

Genetic variation and phylogenetic analysis of Porcine circovirus type 2 in China from

Porcine Circovirus Type 2 (PCV2) Vaccines in the Context of Current Molecular Epidemiology

MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets

Molecular epidemiology and genetic variation analyses of porcine circovirus type 2 isolated from Yunnan Province in China from 2016-2019

Modified indirect porcine circovirus (PCV) type 2-based and recombinant capsid protein (ORF2)-based enzyme-linked immunosorbent assays for detection of antibodies to PCV

Molecular evolution of porcine circovirus type 2 genomes: phylogeny and clonality

Genetic and experimental comparison of porcine circovirus type 2 (PCV2) isolates from cases with and without PCV2-associated lesions provides evidence for differences in virulence

Differences in virulence among porcine circovirus type 2 isolates are unrelated to cluster type 2a or 2b and prior infection provides heterologous protection

Case-control study on the association of porcine circovirus type 2 and other swine viral pathogens with postweaning multisystemic wasting syndrome

Development of a loopmediated isothermal amplification method to rapidly detect porcine circovirus genotypes 2a and 2b

Detection and genotyping of Porcine circovirus 2 (PCV-2) and detection of Porcine circovirus 3 (PCV-3) in sera from fattening pigs of different European countries

Development and evaluation of multiplex real-time RT-PCR assays for the detection and differentiation of foot-and-mouth disease virus and Seneca Valley virus 1

A multiplex real-time PCR assay for the detection and differentiation of the newly emerged porcine circovirus type 3 and continuously evolving type 2 strains in the United States

Genetic diversity and prevalence of porcine circovirus type 3 (PCV3) and type 2 (PCV2) in the Midwest of the USA during

Development of a real-time PCR assay for detection of African swine fever virus with an endogenous internal control

Real-Time PCR Detection Patterns of Porcine Circovirus Type 2 (PCV2) in Polish Farms with Different Statuses of Vaccination against PCV2

PCV2d-2 is the predominant type of PCV2 DNA in pig samples collected in the US during

Phylogenetic analysis of porcine circovirus type 2 (PCV2) between

J o u r n a l P r e -p r o o f