key: cord-0741016-hzbh7s2t
authors: Wang, Yin; Feng, Yuan; Zheng, Wanglong; Noll, Lance; Porter, Elizabeth; Potter, Megan; Cino, Giselle; Peddireddi, Lalitha; Liu, Xuming; Anderson, Gary; Bai, Jianfa
title: 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
date: 2019-03-20
journal: J Virol Methods
DOI: 10.1016/j.jviromet.2019.03.011
sha: e2aeb76496f7eb1cca212bd737dd6db43f81d5f0
doc_id: 741016
cord_uid: hzbh7s2t

A multiplex quantitative real-time polymerase chain reaction (mqPCR) assay was developed and validated for the detection and differentiation of porcine circovirus type 3 (PCV3) and type 2 (PCV2) strains. The assay coverage was 97.9% (184/188) for PCV3 and 99.1% (1889/1907) for PCV2 sequences that were available from the current GenBank database. The PCR amplification efficiencies were 98–99% for plasmids, and 92–96% for diagnostic samples, with correlation coefficients all greater than 0.99. The limit of detection (LOD) determined as plasmid copies per reaction was 17 for PCV3 and 14 for PCV2. The assay specifically detected the targeted viruses without cross reacting to each other or to other common porcine viruses. Among 336 swine clinical samples collected in 2018, 101 (30.1%) were PCV3 positive, 56 (16.7%) were PCV2 positive and 18 (5.4%) were co-positives. Sixty selected PCV3 positives were confirmed by Sanger sequencing, and 53 of the 56 PCV2 positive samples were tested positive by another validated PCR assay.

Different from the non-pathogenic PCV1 strains (Tischer et al., 1974) , PCV2 is considered a major swine pathogen causing porcine circovirus-associated diseases (PCVAD) including post-weaning multisystemic wasting syndrome (PMWS), porcine dermatitis and nephropathy syndrome (PDNS), porcine respiratory disease complex (PRDC), enteritis and reproductive failure (Opriessnig et al., 2007; Tischer et al., 1986) . Since its first identification in 1998, PCV2 has been reported worldwide and has caused significant economic losses to the swine industry (Allan et al., 1998) .

In 2015, another porcine circovirus, named PCV3, was found in the US that can cause PCVAD-like clinical symptoms similar to those caused by PCV2. The symptoms include PDNS and reproductive failure, and cardiac and multi-organ inflammation (Palinski et al., 2017; Phan et al., 2017) . PCV3 has also been detected in several other countries including Brazil, China, Thailand, Sweden, Denmark, Italy, Poland and Spain (Chen et al., 2017; Franzo et al., 2018; Kedkovid et al., 2018; Stadejek et al., 2017; Tochetto et al., 2018; Xu et al., 2018; Ye et al., 2018) .

Porcine circoviruses are small non-enveloped viruses with a circular single-stranded DNA genome, belonging to genus Circovirus in the family of Circoviridae. In the viral genome, ORF1 codes for the replicase (Rep) protein and ORF2 for the capsid (Cap) protein. The Rep protein is a non-structural protein and functions for the viral replication, while the structural Cap protein dominates immunogenicity (Cheung, 2003; Mankertz et al., 1998a, b; Nawagitgul et al., 2002) . As a DNA virus, PCV2 has high evolutionary mutation rate of 1.2 × 10 −3 substitutions/ site/year, and has been constantly evolving. With the emerging viral strains, PCV2 has been divided into 5 genotypes, namely PCV2a-2e strains, according to the diversity level of the ORF2 nucleotide sequences (Franzo et al., 2016; Segales et al., 2008; Xiao et al., 2016) . The continued mutation in the PCV2 genome (Eddicks et al., 2015; Yang et al., 2018) to rapidly detect and differentiate the two important circoviruses, with significantly improved diagnostic coverage to current field strains.

Cell culture isolates of PCV2a, PCV2b and PCV2d, porcine reproductive and respiratory syndrome virus type 2 (PRRSV-2) and Seneca Valley virus 1 (SVV-1) that were previously confirmed by culture and sequencing were used in this study.

According to the 188 PCV3 whole genome sequences that were available from the GenBank database, one set each of primers and probes were designed from ORF1 and ORF2. Both PCV3 probes were labeled with 5′-FAM and 3′-BHQ1. Based on 1907 PCV2 whole genome sequences from the GenBank database, two sets of primers and probes were chosen from ORF1 and ORF3 respectively. The PCV2 probes were labeled with 5′-VIC and 3′-BHQ1. A conserved swine gene, serum beta-2-microglobulin (SB2M) was applied as internal control in the assay, coding for a small membrane protein that may be involved in immune system regulation (Xie et al., 2003) . The SB2M probe was labeled with 5′-Cy5 and 3′-BHQ2. Information of all primers and probes are given in Table 1 .

The viral DNA was extracted from 140 μl of clinical samples or cell culture by QIAamp Viral RNA Mini Kit (Qiagen, MD) or ZR Viral DNA/ RNA Kit (Zymo Research, CA) according to the manufacturer's recommendations, and stored at −80°C until use. To ensure that the obtained PCV3 genome was from a single genome, and not amplified from more than one genome, a full length PCV3 genome was amplified by a single pair of tail-to-tail primers that overlapped at a unique PstI site in the PCV3 genome (illustrated in Fig. 1 ). The full-genome amplicon of PCV3 was then digested with the PstI restriction enzyme and cloned into pACYC177 cloning vector. The PCV2 fragments containing the assay targets were also amplified, and the PCR product was cloned into the pCR ™ 2.1 vector using the original TA Cloning kit (Invitrogen, CA). The presence of cloned inserts was confirmed by gel electrophoresis and subsequent DNA sequencing. The primers used for cloning are also listed in Table 1 .

All PCR reactions were performed in a 20 μL reaction composed of 4 μL of DNA samples prepared as described (Shi et al., 2016) , 0.4 μM each of forward and reverse PCR primers, 0.2 μM each of probes, and 10 μl of 2X iQ™ Multiplex Powermix (Bio-Rad, CA). The parameters for thermocycling start with an initial denaturation at 94°C for 10 min, followed by 45 cycles of 94°C for 15 s and 60°C for 45 s. The cycle threshold (Ct) values were generated with CFX96 TouchTM Real-Time PCR Detection System and standard curve results were analyzed with Bio-Rad CFX Manager 3.0.

The analytical sensitivity was determined by generating standard curves with triplicates of 10-fold serial dilutions of control constructs and positive clinical samples, and a cultured PCV2d isolate. To obtain accurate LODs, 2-fold serial dilutions were prepared from the last 10fold dilution to fine-tune the least concentration that still amplifies. The assay specificity was evaluated with cell cultures (PCV2a, PCV2b, PCV2d, PRRSV-2 and SVV-1) and clinical samples positive to specific pathogens (PCV3, PRRSV-2, groups A, B and C swine rotaviruses, swine Sketch of PCV3 genome amplification using tail-to-tail overlapping primers that were designed to have an unique type II restriction endonuclease, PstI recognition site in the overlapping region that was later used for fullgenome cloning. This is to ensure that the genome amplified is from the same genome, and not from separate genomes.

influenza virus (SIV), porcine epidemic diarrhea virus (PEDV), porcine delta coronavirus (PDCoV), porcine parvovirus (PPV) and porcine parainfluenza virus type 1 (PPIV-1) (summarized in Table3).

To compare the performances of the singular PCV2 and PCV3 assays with the multiplex assay, standard curves were also generated for both singular assays and at the multiplexed condition that composed PCV2, PCV3 and the internal control. The performance of multiplex and singular assays was compared using PCV3 and PCV2 standard control constructs and a PCV3 positive sample and a cultured PCV2 isolate. The individual viral DNA was subjected to singular assays while a mixture of equal molarity of PCV2 and PCV3 DNA was used for the multiplex assay. Standard curves were generated and PCR efficiency analysis between the singular and the multiplex assays was performed.

A total of 336 clinical porcine samples collected in 2018 with different sample types including serum, oral fluid, feces and stomached organ tissue homogenates of heart, lung, thymus, tonsil, liver, spleen, kidney, fetus, placenta, lymph nodes and intestine, were collected from Kansas State Veterinary Diagnostic Laboratory (KSVDL). The mqPCR assays were performed on these samples for PCV3 and PCV2 detections, and for potential coinfection identifications.

Based on 188 PCV3 full genome sequences in the GenBank, two sets of primers and probes, targeting the respective cap gene and rep gene, were designed for PCV3 detection. The cap gene set matched 92.6% (174/188) and the rep gene set matches 89.4% (168/188) of the strains with an overall coverage of 97.9%. The four mismatching strains all had a single nucleotide variation in the primers or probes (data not shown). The PCV2 sets were designed based on 1907 PCV2 whole genome sequences from the GenBank, including five genotypes, PCV2a-2e. To ensure high coverage of divergent PCV2 genomes, two sets of primers and probes were designed in ORF1 and ORF3 genes with coverages of 94.8% (1808/1907) and 90.5% (1726/1907), respectively. The combined coverage of the two PCV2 sets was 99.1% (1889/1907) based on an in silico analysis. 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 template, and analyzed by plotting their Cts against dilution factors. The results indicated that the PCR amplification efficiencies were 98.9% for PCV3 and 98.5% for PCV2 with correlation coefficient (R 2 ) both greater than 0.995 ( Fig. 2A) . For more accurate determination of limit of detection (LOD) of the assay, the cloned standard controls from the last 10-fold dilution were further diluted by 2-fold serial dilutions. The results indicated that LODs were 17 copies per PCR reaction for PCV3 and 14 copies per reaction for PCV2.

To make sure that the assay sensitivity was not compromised when tested on virus and diagnostic samples, a PCV2d cell culture isolate and PCV3 positive clinical samples (we have not been able to culture PCV3) were also used. The standard curves generated by 10-fold dilutions were plotted with Cts versus dilution factors and showed that the amplification efficiencies were 92.9% for PCV3 and 95.3% for PCV2 with R 2 greater than 0.995 for both viruses. The LOD of the PCV2d cell culture isolate was around 1.4 TCID50 per reaction. A more accurate Ct cutoff for detection limit was also determined by 2-fold serial dilutions starting with the last 10-fold dilution that still generated signals. The results indicated that the cutoff for both PCV3 and PCV2 positive for diagnostic samples was Ct 37. The standard curve of the internal control SB2M showed that there is no inhibition to PCR amplification (Fig. 2B) .

From mqPCR data, 60 selected PCV3 positive samples with Ct values ranging from 18 to 35 were verified by Sanger sequencing, indicating a 100% diagnostic sensitivity. For PCV2, 53 of 56 positive samples were verified as positives by a previously validated PCV2 qPCR assay that is currently used at KSVDL. Thus for the 53 positives identified by the current KSVDL assay, the diagnostic sensitivity of the newly developed assay is also 100%. The new assay may have higher diagnostic sensitivity as it detected three more positive samples.

As shown in Fig. 2 , the mqPCR generated similar correlation coefficients and PCR amplification efficiencies to those generated by singular qPCR reactions. Both multiplex and singular assays had R 2 greater than 0.994 and PCR amplification efficiencies between 92.9% and 96.9%. The Ct values for serial dilutions of clinical positive samples in singular and multiplex reactions were nearly identical (Table 2) indicating that multiplexing is not reducing the assay's sensitivity.

The specificity of primers and probes was evaluated by an in silico analysis with NCBI primer design tool, which presented a unique viral target for each set of assay. Assay specificity was further tested experimentally with a PCV3 dominated clinical sample as confirmed by next generation sequencing, viral isolates (PCV2a, PCV2b, PCV2d, PRRSV-2, and SVV-1), and clinical samples that were previously tested positive to non-target, specific pathogens. The results demonstrated that assay specifically detected PCV2 and PCV3 positive samples without cross detection and no positive targets were detected from samples that were positive to other common porcine viruses that included 16 PEDV positive samples, 16 PRRSV positives, 7 SIV positives, 1 PPIV positive, 1 PDCoV positive, 1 PPV positive, and one positive sample each to group A, group B and group C rotaviruses. Results showed that there was no signal generated on those non-target positive samples indicating a good specificity of the assay (Table3).

The 336 clinical porcine samples collected in 2018 with different sample types were tested for PCV2 and PCV3 to determine assay performance and viral prevalence in the field. Of the 336 samples, 56 (16.7%) were PCV2 positives, 101 (30.1%) were PCV3 positives. Among the 139 PCV2 and PCV3 positive samples, 18 (5.4% of total number of samples, and 12.9% of positive samples) were positive for both PCV3 and PCV2 viruses ( Table 4 ). The house keeping gene, SB2M, that was included as an internal control, generated Cts on these diagnostic samples mostly in between 20-30, indicating an efficient nucleic acid extraction and no PCR inhibition observed.

PCVAD is a common disease in swine production systems and has Y. Wang, et al. Journal of Virological Methods 269 (2019) 7-12 caused significant economic losses. PCV2 is a major pathogen involved in PCVAD (Segales et al., 2008) . A new porcine circovirus, PCV3, was identified in 2015 and is causing PCVAD-like clinical signs, including PDNS and reproductive failure, and cardiac and multi-organ inflammation, and thus has attracted researchers' and producers' attention (Chen et al., 2017; Palinski et al., 2017; Phan et al., 2017; Stadejek et al., 2017) . Considering the similar clinical symptoms caused by PCV2 and PCV3, and that PCV2 assays built many years ago are no longer covering the majority of field strains, there was a need to develop a mqPCR assay that was based on the most current sequencing data and capable of detecting the two viruses with high diagnostic coverage to field strains. Although several molecular detection methods were developed for rapid detection of PCV3 and PCV2 (Kim et al., 2017; Li et al., 2018; Zhang et al., 2018) , our mqPCR assay showed apparent advantages. With the help of bioinformatics tools, 1907 PCV2 full genome sequences that were currently available in the GenBank were downloaded and analyzed to achieve high coverage in the assay design stage. We believe that the coverage of the primers and probes over available sequences in the designing stage is the best estimation of assay's future diagnostic sensitivity against field strains (Bai et al., 2018) . Compared to the three published duplex assays, our assay has much higher coverages: 97.9% (184/188) for PCV3 and 99.1% (1889/1907) for PCV2. All 4 PCV3 strains that our primers or probes mismatch to are all single nucleotide variations, and very often these oligoes that have single nucleotide mismatches will not affect their ability of binding to the templates. Therefore, the actual diagnostic sensitivity for PCV3 can be higher than this estimation. Using two non-overlapping targets for a given virus may not increase the assay's analytical sensitivity, but it will increase field strain coverage, and will help to detect strains with additional mutations, as the chance for a strain to have mutations on both target sites is small. Also, compared to the LOD of 50 copies per reaction in Kim et al. (2017) and 90 copies per reaction in Zhang et al. (2018) , our assay appeared to be more sensitive: 17 copies per reaction of PCV3 and 14 copies per reaction of PCV2. The LOD of PCV2 was also evaluated with cell culture, with an LOD of around 1.4 TCID50 per reaction. In addition, the internal control, swine SB2M gene, is included in our assay to monitor nucleic acid extraction efficiencies and potential PCR inhibitions in order to reduce the false-negative rate, which were not used in the previous studies.

In our study, different types of porcine samples were collected and subjected to the mqPCR testing. The viral targets were detected from porcine serum, oral fluid, feces, intestines and tonsil, indicating that our assay can be used for a wide range of sample types that are encountered in routine diagnostic operations. It is interesting to see, with our limited data, that the PCV3 positive rate was 30.1%, which was much higher than the PCV2 positive rate of 16.7%. We currently do not have the data to indicate whether the application of PCV2 vaccines in recent years play a role in reducing PCV2 prevalence. We will keep monitoring the field samples to see if this prevalence data holds true in the future.

In conclusion, the newly developed and validated mqPCR assay enables us to perform rapid, sensitive and specific detection and differentiation of PCV3 and PCV2 strains in clinical samples. In silico analysis and clinical sample testing indicated that the assay has high strain coverage for both viruses. Our limited prevalence data indicated that PCV3 strains have been widely distributed, and may be more prevalent than PCV2 strains currently in the US.

The authors declare that they have no competing interests. 

Isolation of porcine circovirus-like viruses from pigs with a wasting disease in the USA and Europe

A multiplex real-time PCR assay, based on invA and pagC genes, for the detection and quantification of Salmonella enterica from cattle lymph nodes

Detection and genome sequencing of porcine circovirus 3 in neonatal pigs with congenital tremors in South China

The essential and nonessential transcription units for viral protein synthesis and DNA replication of porcine circovirus type 2

Detection of a new cluster of porcine circovirus type 2b strains in domestic pigs in Germany

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

Full-genome sequencing of porcine circovirus 3 field strains from Denmark, Italy and Spain demonstrates a high within-Europe genetic heterogeneity

Porcine circovirus type 3 (PCV3) infection in grower pigs from a Thai farm suffering from porcine respiratory disease complex (PRDC)

Multiplex real-time polymerase chain reaction for the differential detection of porcine circovirus 2 and 3

A duplex real-time PCR assay for the simultaneous detection of porcine circovirus 2 and circovirus 3

Identification of a protein essential for replication of porcine circovirus

Transcription analysis of porcine circovirus (PCV)

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

Porcine circovirus type 2 associated disease: update on current terminology, clinical manifestations, pathogenesis, diagnosis, and intervention strategies

A novel porcine circovirus distantly related to known circoviruses is associated with porcine dermatitis and nephropathy syndrome and reproductive failure

Erratum to: detection of a novel circovirus PCV3 in pigs with cardiac and multi-systemic inflammation

PCV-2 genotype definition and nomenclature

A multiplex real-time PCR panel assay for simultaneous detection and differentiation of 12 common swine viruses

First detection of porcine circovirus type 3 on commercial pig farms in Poland

Characterization of papovavirus-and picornavirus-like particles in permanent pig kidney cell lines

Studies on epidemiology and pathogenicity of porcine circovirus

Full-Genome Sequence of Porcine Circovirus type 3 recovered from serum of sows with stillbirths in Brazil

Beta 2-microglobulin as a negative regulator of the immune system: high concentrations of the protein inhibit in vitro generation of functional dendritic cells

Detection and phylogenetic analysis of porcine circovirus type 3 in central China

Phylogenetic and genetic variation analyses of porcine circovirus type 2 isolated from China

Detection and genetic characterisation of porcine circovirus 3 from pigs in Sweden

Phylogenetic analysis of porcine circovirus type 3 and porcine circovirus type 2 in China detected by duplex nanoparticle-assisted PCR

This research is supported by Swine Health Information Center and Kansas State Veterinary Diagnostic Laboratory. We thank Dr. Ying Fang and Dr. Richard Hesse of our department for sharing the viral strains.