- Research article
- Open Access
Duplex fluorescence melting curve analysis as a new tool for rapid detection and differentiation of genotype I, II and Bartha-K61 vaccine strains of pseudorabies virus
BMC Veterinary Researchvolume 14, Article number: 372 (2018)
Recently, pseudorabies (PR) outbreaks have been reported in a large number of swine herds vaccinated with the Bartha-K61 vaccine in China, the current pseudorabies virus (PRV) belonging to Genotype II is differential genetically from Bartha-K61 vaccine belonging to Genotype I. Furthermore, it has been proved that the Bartha-K61 vaccine cannot provide sufficient protection against the current PRVs in China. Therefore, the accurate and rapid identification of PRVs is essential. The objective of this study is to develop a duplex fluorescence melting curve analysis (FMCA) capable of rapid, simple, high-throughput differentiation of Chinese, European/American and Bartha-K61 vaccine strains of PRV.
Primers 6F/6R and probes P1/P2, combined with three recombinant plasmids p-B (Bartha-K61), p-N (Genotype I), and p-H (Genotype II), were used to establish the Bicolor FMCA. FAM Tm values (probe P1) and HEX (probe P2) channels of p-B were used as reference values. Tm differences (ΔTm) between detected samples and reference plasmid p-B were calculated in each channel. Bartha-K61 vaccine samples had ΔTm values of ±1 °C in both FAM and HEX channels, Genotype I samples had ΔTm values of ±1 °C in the FAM channel and 4.38 ± 1 °C in the HEX channel, and Genotype II samples had ΔTm values of 6.52 ± 1 °C in the FAM channel and 4.38 ± 1 °C in the HEX channel. The minimum detection limit of the duplex FMCA was approximately 1 × 100 copies per reaction for p-B, p-N, and p-H. The duplex FMCA technique was used to detect and different 198 suspected clinical samples, of which 18 (9%) were positive for Genotype II strains and eight (4%) were positive for Bartha-K61 vaccine strains, and the results were compared with sequencing and phylogenetic analyses, which confirmed that the Bicolor FMCA worked correctly for all samples.
In this study, we developed a duplex FMCA of dual-labeled, self-quenched probes that was performed for rapid detection and differentiation of Genotype I, II and Bartha-K61 vaccine strains of PRV. The duplex FMCA was rapid, simple, and high-throughput, and will likely prove useful for molecular epidemiological investigations and pathogen surveillance of PRV.
Pseudorabies (PR), also known as Aujeszky’s disease (AD), is caused by pseudorabies virus (PRV), which is also known as suid herpesvirus 1 (SuHV-1) or Aujeszky’s disease virus (ADV). PR is an economically important disease in the pig industry worldwide. Swine are the only known natural host of PRV, but nearly all mammals are susceptible to infection . PRV infection is characterised by neurological symptoms and death in newborn piglets, respiratory disorders in fattening pigs, and reproductive failure in pregnant pigs . PR can be controlled using attenuated Bartha-based vaccines and in combination with a worldwide eradication program. However, since the latter part of 2011, PR outbreaks have been reported in a large number of swine herds vaccinated with the Bartha-K61 vaccine in China [3,4,5,6,7,8]. The Bartha-K61 vaccine did not provide full protection against the current PRVs [3, 9]. Current evidence suggests PRV can be divided into two distinct clusters, with Chinese strains belonging to genotype II and European/American strains in genotype I, which contains Bartha-K61 vaccine [10,11,12,13]. Therefore, the accurate differentiation of the PRVs is essential.
Traditionally, etiological diagnostic methods for PR include virus isolation , DNA hybridisation [15,16,17], immunofluorescence , loop-mediated isothermal amplification (LAMP) assay , real-time polymerase chain reaction (PCR) assay , and recombinase polymerase amplification assays (RPA) , none of which are able to distinguish between PRV strains. Several assays, including the LAMP assay , nanoparticle-assisted PCR assay , and two multiplex real-time PCR (qPCR) assays , have been developed to differentiate between wild-type (WT) PRV and gene-deleted PRV vaccines. However, these assays are unable to discriminate between WT PRV strains. Triplex real-time PCR  is applicable for differential detection of classical, variant and Bartha-K61 vaccine strains of PRV. However, classical PRV strains in the present study include both Chinese strains and European/American PRVs, which are evolutionarily distinct clades [10, 24].
Probe-based fluorescence melting curve analysis (FMCA) is a powerful tool for single-nucleotide polymorphism (SNP) genotyping of target sequences based on melting temperature generated by thermal denaturation of the probe-target hybrid. Recently, FMCA was used in the detection of bacterial resistance  and virus genotyping . The aim of the present study was to develop a duplex FMCA capable of rapid, simple, high-throughput differentiation of Chinese, European/American and Bartha-K61 vaccine strains of PRV.
The developed Bicolor FMCA of dual-labeled, self-quenched probes was able to discriminate 92 samples within 2 h when the template DNA was ready. This assay could identify Chinese strains (Genotype II), European/American strains (Genotype I) and Bartha-K61 vaccine strains of PRV based on Tm value and fluorescence type, as shown in Fig. 1b and Fig. 1c. For each sample, two Tm values were obtained from the FAM and HEX channels, corresponding to probes P1 and P2, respectively. P1 was designed to match Genotype II but mismatch Genotype I (including Bartha-K61 vaccine strains), while P2 was designed to match Genotype I and Genotype II but mismatch Bartha-K61 vaccine strains. As shown in Table 1, in the FAM channel, Genotype I (including Bartha-K61 vaccine strains) yielded Tm values of 70.41 ± 0.49 °C (0.49 °C = 3SD, n = 63), and Genotype II Tm values were 76.93 ± 0.48 °C (n = 32). In the HEX channel, Bartha-K61 vaccine strains gave Tm values of 71.62 ± 0.18 °C (n = 32), while Genotype I and Genotype II Tm values were 76.00 ± 0.48 °C (n = 63). In each assay, recombinant plasmids p-B (Bartha-K61 Vaccine), p-N (Genotype I), and p-H (Genotype II) were used as positive controls, and the NTC sample served as the negative control. In order to simplify data analysis, Tm values in FAM and HEX channels of recombinant plasmid p-B was used as reference values, and Tm differences (ΔTm) between detected samples and reference plasmid p-B were calculated in each channel. As shown in Table 1, Tm values were highly reproducible, with a maximum 3SD value of 0.49 (ΔTm = 6.52 °C, FAM channel) and 0.48 (ΔTm = 4.38 °C, HEX channel). Bartha-K61 vaccine samples had ΔTm values of ±1 °C in both the FAM and HEX channels, Genotype I samples had ΔTm values of ±1 °C in the FAM channel and 4.38 ± 1 °C in the HEX channel, and Genotype II samples had ΔTm values of 6.52 ± 1 °C in the FAM channel and 4.38 ± 1 °C in the HEX channel.
Sensitivity, specificity and reproducibility
Serial dilutions of the three recombinant plasmids p-B, p-N, and p-H ranged from 1 × 109 to 1 × 100 copies/μL when tested with the Bicolor FMCA, and fluorescence signals corresponding to melting curves were obtained between 1 × 100 and 1 × 109 copies per reaction for each plasmid. No specific signals were detected and corresponding melting curves were not obtained with the negative control (Additional file 1). Therefore, the detection limit of the Bicolor FMCA was approximately 1 × 100 copies per reaction for recombinant plasmids p-B, p-N, and p-H.
The specificity of the Bicolor FMCA was validated using DNA or cDNA samples from three PRV isolates, five other viruses known to cause similar clinical symptoms in pig, and three positive control plasmids (Fig. 2). Based on the data analysis method described above, detected samples included two Bartha-K61 vaccine samples (including plasmid p-B and Bartha-K61), two Genotype I samples (including plasmid p-N and NIA3), and two Genotype II samples (including plasmid p-H and HDDJ). However, no specific melting peak was detected with other non-targeted pig viruses such as CSFV, PPV, PCV2, PRRSV, and JEV in the FAM and HEX channels. These results suggest that the designed primers and probes were highly specific and selective for their target viruses and exhibited no cross-reactivity with other viruses.
To evaluate the run-to-run reproducibility of the Bicolor FMCA, given the narrow Tm windows, the assay was run on different days using PRV Bartha-K61, NIA3, and HDDJ samples. Melting curve analysis for samples run on different days yielded a maximum 3SD value for Tm of no more than 0.5 for all melting peaks (Table 1), demonstrating the high reproducibility of this assay.
Agreement among the bicolor FMCA, sequencing and phylogenetic analysis
Of 198 suspected clinical samples were tested by the Bicolor FMCA technique, 18 (9%) were positive for Genotype II strains and eight (4%) were positive for Bartha-K61 vaccine strains (Table 2), but Genotype I strains were not detected.
To confirm that the clustering pattern resulting from the difference-melt curves of the different genotypes, the gC gene of each sample was sequenced and compared with consensus sequences (Additional file 2). Ten control sequences were included from the GenBank databases when sequences were aligned. Sequences of the 26 positive samples aligned well with the GenBank controls and with each other, confirming that the Bicolor FMCA worked correctly for all samples (Additional files 2 and 3). The relationship between individuals is shown in the cladogram.
To additionally confirm the eight Bartha-K61 strains, the fragment encompassing the gE gene was amplified by the primers gE-F and gE-R. The result showed that no fragment was amplified from the eight Bartha-K61 strains (Additional file 4). Furthermore, the eight samples were obtained from Bartha-K61–vaccinated pigs.
Recently, PR outbreaks have been reported in a large number of swine herds vaccinated with the Bartha-K61 vaccine in China [3,4,5,6,7,8]. And the phylogenetic analysis based on full-length PRV sequences and based on only the gC gene and other genes revealed that PRV could be divided into two major genotypes, with Genotype I containing isolates from Europe, America, Australia, and Asia, and Genotype II consisting mainly of isolates from China, including classical strains such as Fa and Ea, as well as recently emerged variant strains [10,11,12,13]. Furthermore, it has been proved that the Bartha-K61 vaccine belonging to Genotype I cannot provide full protection against the current PRVs in China belonging to Genotype II [3, 9]. Therefore, it is very important to develop a method to identify pigs infected with PRV genotypes or immunized with the PRV Bartha-K61 vaccine strain.
In this study, we developed a simple and rapid method for detection and differentiation of PRV strains using a duplex FMCA approach with dual-labeled, self-quenched probes. This method proved to be easier, faster and more straightforward than conventional differentiation methods such as virus isolation combined with sequence analysis, which is often used to distinguish different strains .
The TaqMan probe has a random coil conformation that leads to fluorescence quenching or only weak fluorescence unless the probe is either hybridised with its target or digested. In classic TaqMan real-time PCR, generation of the fluorescence signal is based on probe degradation by the 5′-3′ exonuclease activity of DNA Taq polymerase in each round of amplification . In FMCA, after PCR, a larger amount of fluorescence is emitted upon probe hybridisation with target sequences because the fluorophore is no longer in close proximity to the quencher. In this study, TaKaRa LA Taq DNA Polymerase combined Taq DNA Polymerase with 5′-3′ exonuclease activity and a DNA proofreading polymerase with 3′-5′ exonuclease activity were employed during PCR analyses. Hence, as described previously , the effects of hydrolysis of the probe during PCR on the melting curve profile require assessment. We compared the melting curves and melting peaks obtained from experiments in which probes were added either before or after the PCR run (Additional file 5). The results showed that the background fluorescence intensities in FAM and HEX channels increased when probes were added before the PCR run at the beginning of the melting curve analysis (Additional file 5.A and C) because some probes were digested during PCR by the 5′-3′ exonuclease activity of TaKaRa LA Taq DNA Polymerase. However, this did not affect the Tm value of the melting curve or the melting peak of FAM and HEX channels (Additional file 5). Additionally, this did not affect the −dF/dT value of the FAM channel when temperature equal to the Tm value (Additional file 5B), or slightly below the −dF/dT value of the HEX channel (Additional file 5.D). In summary, TaKaRa LA Taq was shown to be suitable for this assay.
In this study, the limit of detection of the Bicolor FMCA was determined to be 1 × 100 copies per reaction for recombinant plasmids p-B, p-N and p-H. Some factors contributed greatly to this. Firstly, only one pair of primers was used during the PCR process, and post-multiplex FMCA could reduce some of the nonspecific interactions if two or more pairs of primers are used . Secondly, the asymmetric PCR post-multiplex FMCA was optimised in terms of the ratio of primer 6F to 6R. In order to maximise the amplified single strand used to hybridise with the probes while retaining the basic PCR reaction conditions, the ratio of primers is important, as described previously [26, 28, 30,31,32,33,34,35,36,37,38]. The results showed that detection by the Bicolor FMCA was highest when the ratio of primer 6F to 6R was 1:10 (0.04 μM: 0.4 μM) and probe P1 and P2 was 0.2 μM. Thirdly, the PCR amplification system containing TaKaRa LA Taq and GC Buffer I improved the PCR amplification efficiency significantly over the conventional PCR amplification system for the high GC content of the PRV genome.
Since the 1990s, most pig herds in China have been immunised with the Bartha-K61 vaccine, which is a live attenuated strain (Bartha, 1961) . However, co-infection with more than one PRV genotype has not been widely reported, possibly due to limited differential detection methods in the past. Using the dual presence of PRV recombinant plasmids p-B and p-H as an example, we tested the ability to distinguish the two genotypes when they are simultaneously present in the same reaction. By varying their relative ratios, the results showed that levels as low as 10% (1:9) of one genotype in the presence of 90% (9:1) of the other genotype could be detected when the overall copy number was 108 copies per reaction (Fig. 3).
It should be pointed out that this method is only applicable for genotypes from Chinese, European/American and Bartha-K61 vaccine strains of PRV because it was developed based on SNP differences. Thus, the method is not suitable for the genotyping of PRV genetically engineered vaccine strains such as HB-98 (MF434040), HB-2000 (MF434041), and SA215 (EU719635.1) strains. Also, the method is unable to subdivide either Chinese strain into ancient and recently emerged strains, or European/American PRVs into six subgroups . A more inclusive method needs to be developed in the future for this purpose.
In summary, the newly developed Duplex FMCA is a specific and sensitive method for the differential detection of Chinese, European/American and Bartha-K61 vaccine strains of PRV, and will likely prove useful for molecular epidemiological investigations and pathogen surveillance of PR.
Strain collection and DNA extraction
A total of 28 PRV samples (Table 2) and five other viruses were tested in this study. Three positive strains were included to establish the Bicolor FMCA: PRV Bartha-K61 (Bartha-K61 Vaccine, Boehringer CO. LTD, Germany), NIA3 (Genotype I, the NIA3 strain), and HDDJ (Genotype II, the HDDJ strain) and other 25 field samples were isolated in our laboratory; the other five viruses are known to cause similar clinical symptoms in pigs and consisted of classical swine fever virus (CSFV) C strain (Winsun Bio, Guangzhou, China), porcine parvovirus (PPV) GD3 strain (isolated in our laboratory), porcine circovirus type 2 (PCV2) BL strain (isolated in our laboratory), porcine reproductive and porcine respiratory syndrome virus (PRRSV) GDr180 strain (Winsun Bio, Guangzhou, China), and Japanese encephalitis virus (JEV) SA14–14-2 strain (Wuhan Keqian Biology., Ltd., Wuhan, China). Viral RNA/DNA was extracted using the MiniBEST Viral RNA/DNA Extraction Kit version 4.0 (Takara, Dalian, China) according to the manufacturer’s instructions.
Primer and probe design
Main primers and probes were designed to target the gC gene of PRV (Fig. 1a). Primers 6F (5’-GCCAACGCCTCCCTCGCCCAC-3′) and 6R (5’-CGACGCACACCGCCCGGAAGG-3′), and probes P1 (BHQ1–5’-CGCCAACGGCACCGAGGTCCGCAGCGC-3’-FAM) and P2 (HEX-5’-TCGGCCTGAGCGCGCCGCCCGTC-3’-BHQ1; differential nucleotides are underlined) were used to establish the Bicolor FMCA. Probe P1 was designed to differentiate between Genotype II and Genotype I strains (including Bartha-K61 vaccine strains), and probe P2 was designed to distinguish Bartha-K61 vaccine strains from other strains. The single nucleotide polymorphisms (SNPs, Fig. 1a), C699G (probe P1) and A756G (probe P2), in the gC gene, were chose based on a sequence alignment of all available PRV gC gene by a BLAST search in GenBank databases from the National Center for Biotechnology Information (NCBI). The primers wai3F (5’-CCGTTTCCTGATTCACGCCCACG-3′) and wai4R (5’-CGATGGCTCGGTTCAACGC-3′) and wai4BR (5’-CGGGTCGGACTCGCTGTCGTTTA-3′) were used to generate recombinant plasmids, which were used to determine the limit of detection (LOD) of the assay. The primers were also used for sequencing of the gC gene. Two reverse primers wai4R and wai4BR bind to similar regions of the PRV genome, but primer wai4R is specific for Chinese strains (Genotype II), combined with primer wai3F, used to amplify a 1594-bp amplicon; while primer wai4BR is specific for European/American strains (Genotype I), combined with primer wai3F, used to amplify a 1577-bp amplicon. This experimental design could allow both PRV genotypes (I and II) to be amplified.
The primers gE-F (5’-TTGAGACCATGCGGCCCTTTCTGCT-3′) and gE-R (5′- GACCGGTTCTCCCGGTATTTAAGCG-3′) were used to confirm the Bartha-K61 strains detected in field samples. The primers amplified no fragment from Bartha-K61 strain, or a 1763-bp fragment encompassing the gE gene of PRV wild strains, respectively.
All primer and probe sequences were analysed using BLAST to confirm specificity.
Duplex real-time PCR and melting curve analysis
For each sample, 1 μL of template DNA was added to a 10 μL PCR mixture containing 1 × GC Buffer I, 400 μM of each deoxynucleoside triphosphate (dNTP), 0.04/0.4 μM primers 6F/6R, 0.2 μM probe P1/P2, and 0.5 U of TaKaRa LA Taq (Takara, Dalian, China). A No Template Control (NTC) was used as the negative control. Amplification and melting curve analysis was performed in a Roche LightCycler 96 System (Roche, Switzerland) with initial denaturation at 95 °C for 3 min, followed by amplification using 55 cycles of 95 °C for 20 s, 60 °C for 20 s, and 72 °C for 20 s. Melting curve analysis was initiated with denaturation at 95 °C for 10 s, followed by hybridisation at 40 °C for 60 s and a stepwise temperature increment from 40 °C to 97 °C at an Ramp rate of 0.13 °C/s with a duration of 1 s and five readings per °C. Fluorescence was measured at each step in both FAM and HEX channels. Melting curves were obtained by plotting the negative derivative of the fluorescence intensity with respect to temperature (−dF/dT) versus temperature (T). The Tm value from each probe was automatically obtained by identifying the peak of the corresponding melting curve. For best results, 11 ratios of forward primers to reverse primers (0.5:10, 1:10, 2:10, 3:10, 4:10, 5:10, 6:10, 7:10, 8:10, 9:10, 10:10) and four concentrations of two probes (0.1, 0.2, 0.3, 0.4 μM) were set. The optimization of Tm (s) was at a best ratio of 0.04/0.4 μM primers 6F/6R and a concentration of 0.2 μM probe P1/P2.
The analytical sensitivity of the Bicolor FMCA was estimated by serial dilution experiments using three recombinant plasmids p-B, p-N, and p-H. The three plasmids were generated using the pMD18-T Vector Cloning Kit (Takara, Dalian, China), and the insert fragments of the plasmids were amplified from Bartha-K61, NIA3, and HDDJ using primers wai3F/wai4R/wai4BR with the TaKaRa LA PCR Kit version 2.1, respectively. Each 25 μL PCR contained 1 × GC Buffer I, 400 μM of each deoxynucleoside triphosphate (dNTP), 0.4 μM primers wai3F/wai4R/wai4BR, 1.25 U of TaKaRa LA Taq (Takara, Dalian, China), and 2 μL of template DNA. Thermal conditions involved initial denaturation at 95 °C for 3 min, followed by amplification using 45 cycles of 95 °C for 40 s, 62 °C for 40 s, and 72 °C for 90 s. Each plasmid was 10-fold serially diluted with water to concentrations of 1 × 109 copies/μL to 1 × 100 copies/μL. The limit of detection for each genotype was determined from the lowest concentration measured by the Bicolor FMCA.
Specificity and reproducibility
In the Bicolor FMCA, primers 6F/6R and probes P1/P2 were tested for specificity using three PRV positive samples including Bartha-K61, NIA3, HDDJ, and five other viruses known to cause similar clinical symptoms in pig including CSFV C strain, PPV GD3 strain, PCV2 BL strain, PRRSV GDr180 strain, JEV SA14–14-2 strain. DNA samples from three recombinant plasmids p-B, p-N, and p-H were used as positive controls. To evaluate the reproducibility of the assay, Bicolor FMCA was performed on different days using different DNA templates of each positive control sample PRV Bartha-K61, NIA3, and HDDJ. The Tm value of each genotype was measured in each channel.
Screening and genotyping of clinical samples
The Bicolor FMCA was used to screen 198 suspected clinical samples received for laboratory diagnosis of PRV infection from different geographical locations in China. For confirmation of Bicolor FMCA results, each sample was rechecked by virus isolation, followed by DNA sequencing for positive samples.
The genotype of each sample was determined by sequencing using primers wai3F/wai4R/wai4BR with the PCR protocol described as above. PCR products were subjected to bi-directional sequencing (Sangon, Shanghai, China) and the results were analysed using Clustal W with BioEdit software (version 7.0.0, Tom Hall; Ibis Biosciences, Carlsbad, CA, USA). Phylogenetic analysis was performed by the neighbour-joining method with 1000 bootstrap replicates, based on the Kimura 2-parameter model, with joining comparison of the positive samples and control gC gene sequences using MEGA 5.2 software .
The Bartha-K61 strains detected in field samples was additionally confirmed based on the deleted gE gene by the primers gE-F and gE-R. The PCR protocol was following the amplified gC gene described above.
Fluorescence melting curve analysis
Mettenleiter TC. Pseudorabies virus A2 - Mahy, Brian W.J. In: Regenmortel MHVV, editor. Encyclopedia of Virology (Third Edition). Edn. Oxford: Academic Press; 2008. p. 341–51.
Beran GW, Davies EB, Arambulo PV 3rd, Will LA, Hill HT, Rock DL. Persistence of pseudorabies virus in infected swine. J Am Vet Med Assoc. 1980;176(10 Pt 1):998–1000.
An TQ, Peng JM, Tian ZJ, Zhao HY, Li N, Liu YM, Chen JZ, Leng CL, Sun Y, Chang D, et al. Pseudorabies virus variant in Bartha-K61-vaccinated pigs, China, 2012. Emerg Infect Dis. 2013;19(11):1749–55.
Wu R, Bai C, Sun J, Chang S, Zhang X. Emergence of virulent pseudorabies virus infection in northern China. J Vet Sci. 2013;14(3):363–5.
Yu X, Zhou Z, Hu D, Zhang Q, Han T, Li X, Gu X, Yuan L, Zhang S, Wang B, et al. Pathogenic pseudorabies virus, China, 2012. Emerg Infect Dis. 2014;20(1):102–4.
Wang Y, Qiao S, Li X, Xie W, Guo J, Li Q, Liu X, Hou J, Xu Y, Wang L, et al. Molecular epidemiology of outbreak-associated pseudorabies virus (PRV) strains in Central China. Virus Genes. 2015;50(3):401–9.
Yang QY, Sun Z, Tan FF, Guo LH, Wang YZ, Wang J, Wang ZY, Wang LL, Li XD, Xiao Y, et al. Pathogenicity of a currently circulating Chinese variant pseudorabies virus in pigs. World journal of virology. 2016;5(1):23–30.
Fan J, Zeng X, Zhang G, Wu Q, Niu J, Sun B, Xie Q, Ma J. Molecular characterization and phylogenetic analysis of pseudorabies virus variants isolated from Guangdong province of southern China during 2013-2014. J Vet Sci. 2016;17(3):369–75.
Wang CH, Yuan J, Qin HY, Luo Y, Cong X, Li Y, Chen J, Li S, Sun Y, Qiu HJ. A novel gE-deleted pseudorabies virus (PRV) provides rapid and complete protection from lethal challenge with the PRV variant emerging in Bartha-K61-vaccinated swine population in China. Vaccine. 2014;32(27):3379–85.
Ye C, Zhang QZ, Tian ZJ, Zheng H, Zhao K, Liu F, Guo JC, Tong W, Jiang CG, Wang SJ, et al. Genomic characterization of emergent pseudorabies virus in China reveals marked sequence divergence: evidence for the existence of two major genotypes. Virology. 2015;483:32–43.
Wu XM, Chen QY, Chen RJ, Che YL, Wang LB, Wang CY, Yan S, Liu YT, Xiu JS, Zhou LJ. Pathogenicity and whole genome sequence analysis of a pseudorabies virus strain FJ-2012 isolated from Fujian, southern China. The Canadian journal of infectious diseases & medical microbiology = Journal canadien des maladies infectieuses et de la microbiologie medicale. 2017:2017, 9073172.
Song C, Gao L, Bai W, Zha X, Yin G, Shu X. Molecular epidemiology of pseudorabies virus in Yunnan and the sequence analysis of its gD gene. Virus Genes. 2017;53(3):392–9.
Liu H, Li XT, Hu B, Deng XY, Zhang L, Lian SZ, Zhang HL, Lv S, Xue XH, Lu RG, et al. Outbreak of severe pseudorabies virus infection in pig-offal-fed farmed mink in Liaoning Province, China. Arch Virol. 2017;162(3):863–6.
Balasch M, Pujols J, Segales J, Pumarola M. Aujeszky's disease (pseudorabies) virus detection in cerebrospinal fluid in experimentally infected pigs. Vet Microbiol. 1998;60(2–4):99–106.
McFarlane RG, Thawley DG. DNA hybridization procedure to detect pseudorabies virus DNA in swine tissue. Am J Vet Res. 1985;46(5):1133–6.
Belak S, Rockborn G, Wierup M, Belak K, Berg M, Linne T. Aujeszky's disease in pigs diagnosed by a simple method of nucleic acid hybridization. Zentralbl Veterinarmed B. 1987;34(7):519–29.
Wongwatcharadumrong R, Kunavongkrit A, Neumuller M, Linne T. Comparison of virus isolation, immunofluorescence and DNA probe hybridization for detection of pseudorabies virus in experimentally infected pigs. Zentralbl Veterinarmed B. 1992;39(2):91–6.
En FX, Wei X, Jian L, Qin C. Loop-mediated isothermal amplification establishment for detection of pseudorabies virus. J Virol Methods. 2008;151(1):35–9.
Zanella EL, Miller LC, Lager KM, Bigelow TT. Evaluation of a real-time polymerase chain reaction assay for pseudorabies virus surveillance purposes. J Vet Diagn Investig. 2012;24(4):739–45.
Yang Y, Qin X, Zhang W, Li Z, Zhang S, Li Y, Zhang Z. Development of an isothermal recombinase polymerase amplification assay for rapid detection of pseudorabies virus. Mol Cell Probes. 2017;33:32–5.
Zhang CF, Cui SJ, Zhu C. Loop-mediated isothermal amplification for rapid detection and differentiation of wild-type pseudorabies and gene-deleted virus vaccines. J Virol Methods. 2010;169(1):239–43.
Ma X, Cui Y, Qiu Z, Zhang B, Cui S. A nanoparticle-assisted PCR assay to improve the sensitivity for rapid detection and differentiation of wild-type pseudorabies virus and gene-deleted vaccine strains. J Virol Methods. 2013;193(2):374–8.
Wernike K, Beer M, Freuling CM, Klupp B, Mettenleiter TC, Muller T, Hoffmann B. Molecular double-check strategy for the identification and characterization of Suid herpesvirus 1. J Virol Methods. 2014;209:110–5.
Meng XY, Luo Y, Liu Y, Shao L, Sun Y, Li Y, Li S, Ji S, Qiu HJ. A triplex real-time PCR for differential detection of classical, variant and Bartha-K61 vaccine strains of pseudorabies virus. Arch Virol. 2016;161(9):2425–30.
Liu Z, Zhang J, Rao S, Sun L, Zhang J, Liu R, Zheng G, Ma X, Hou S, Zhuang X, et al. Heptaplex PCR melting curve analysis for rapid detection of plasmid-mediated AmpC beta-lactamase genes. J Microbiol Methods. 2015;110:1–6.
Athar MA, Xu Y, Xie X, Xu Z, Ahmad V, Hayder Z, Hussain SS, Liao Y, Li Q. Rapid detection of HCV genotyping 1a, 1b, 2a, 3a, 3b and 6a in a single reaction using two-melting temperature codes by a real-time PCR-based assay. J Virol Methods. 2015;222:85–90.
McGuigan FE, Ralston SH. Single nucleotide polymorphism detection: allelic discrimination using TaqMan. Psychiatr Genet. 2002;12(3):133–6.
El Housni H, Heimann P, Parma J, Vassart G. Single-nucleotide polymorphism genotyping by melting analysis of dual-labeled probes: examples using factor V Leiden and prothrombin 20210A mutations. Clin Chem. 2003;49(10):1669–72.
Elnifro EM, Ashshi AM, Cooper RJ, Klapper PE, Multiplex PCR. Optimization and application in diagnostic virology. Clin Microbiol Rev. 2000;13(4):559–70.
Crockett AO, Wittwer CT. Fluorescein-labeled oligonucleotides for real-time pcr: using the inherent quenching of deoxyguanosine nucleotides. Anal Biochem. 2001;290(1):89–97.
Wittwer CT, Herrmann MG, Moss AA, Rasmussen RP. Continuous fluorescence monitoring of rapid cycle DNA amplification. Biotechniques. 1997;22(1):130–1 134-138.
Logan JM, Edwards KJ, Saunders NA, Stanley J. Rapid identification of campylobacter spp. by melting peak analysis of biprobes in real-time PCR. J Clin Microbiol. 2001;39(6):2227–32.
Howell WM, Jobs M, Brookes AJ. iFRET: an improved fluorescence system for DNA-melting analysis. Genome Res. 2002;12(9):1401–7.
French DJ, Archard CL, Brown T, McDowell DG. HyBeacon probes: a new tool for DNA sequence detection and allele discrimination. Mol Cell Probes. 2001;15(6):363–74.
Tyagi S, Kramer FR. Molecular beacons: probes that fluoresce upon hybridization. Nat Biotechnol. 1996;14(3):303–8.
Belousov YS, Welch RA, Sanders S, Mills A, Kulchenko A, Dempcy R, Afonina IA, Walburger DK, Glaser CL, Yadavalli S, et al. Single nucleotide polymorphism genotyping by two colour melting curve analysis using the MGB eclipse probe system in challenging sequence environment. Human genomics. 2004;1(3):209–17.
Ahn JJ, Kim Y, Lee SY, Hong JY, Kim GW, Hwang SY. Fluorescence melting curve analysis using self-quenching dual-labeled peptide nucleic acid probes for simultaneously identifying multiple DNA sequences. Anal Biochem. 2015;484:143–7.
Lukhtanov EA, Lokhov SG, Gorn VV, Podyminogin MA, Mahoney W. Novel DNA probes with low background and high hybridization-triggered fluorescence. Nucleic Acids Res. 2007;35(5):e30.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and Maximum Parsimony Methods. Molecular Biology and Evolution. 2011;28(10):2731–9.
We thank Dr. Gali Bingga, Vocational and Technical College, Inner Mongolia Agricultural University, for assistance with discussions of the manuscript, and Dr. Jingyun Ma, College of Animal Science, South China Agricultural University, for assistance collecting samples.
This work was supported by the Science and Technology Planning Project of Guangdong Province, China [grant numbers: 2016A040403084, 2015A030401066, 2017A020208030, 2016B020234006, 2016A040403085]. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article. All datasets are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Sensitivity of the duplex FMCA method. Melting curves from duplex FMCA with recombinant control plasmids (A) p-B, (C) p-N, and (E) p-H in the FAM channel, and (B) p-B, (D) p-N, and (F) p-H in the HEX channel, ranging from 1 × 100 to 1 × 109 copies per reaction (from bottom to top). (TIF 1335 kb)
Alignment of gC genes from the detected PRV samples used for Bicolor FMCA. Strains indicated with black discs were used for Bicolor FMCA. (TIF 1008 kb)
Phylogenetic analysis (B) of gC genes from the detected PRV samples used for Bicolor FMCA. Strains indicated with black discs were used for Bicolor FMCA. (TIF 203 kb)
The Bartha-K61 strains detected in field samples were identified by PCR. The PCR results showed no fragments were amplified from the eight Bartha-K61 strains detected in field samples using the primers gE-F and gE-R. Lane 1, DL 2000 marker; Lane 2, Bartha-K61; Lane 3, NIA3; Lane 4, HDDJ; Lane 5, NTC; Lane 6, GDHY5; Lane 7, GDZQ4; Lane 8, JXGZ4; Lane 9, GDSG; Lane 10, GDMM; Lane 11 GDFS; Lane 12, GDZS; Lane 13, GDQY3. (TIF 366 kb)
Influence of probe hydrolysis on the melting curve profile. (A) Melting curves and (B) Melting peaks in the FAM channel, and (C) Melting curves and (D) Melting peaks in the HEX channel for recombinant plasmid p-H. Red and green lines indicate the addition of probes P1 and P2 before and after PCR for plasmid p-H, and blue and pink lines indicate the addition of probes P1 and P2 before and after PCR for NTC. (TIF 1001 kb)