Next Article in Journal
Tendon-Derived Mesenchymal Stem Cells (TDSCs) as an In Vitro Model for Virological Studies in Wild Birds
Previous Article in Journal
Clinical and Biomarker Profile Responses to Rehabilitation Treatment in Patients with Long COVID Characterized by Chronic Fatigue
Previous Article in Special Issue
Modeling Infectious Bursal Disease Virus (IBDV) Antigenic Drift In Vitro
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China

1
Institute of Animal Husbandry and Veterinary Medicine, Fujian Academy of Agricultural Sciences, Fujian Animal Disease Control Technology Development Center, Fuzhou 350013, China
2
Institute of Veterinary Medicine, Jiangsu Academy of Agricultural Sciences, National Research Center for Engineering and Technology of Veterinary Bio-Products, Nanjing 210014, China
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(7), 1453; https://doi.org/10.3390/v15071453
Submission received: 23 May 2023 / Revised: 25 June 2023 / Accepted: 26 June 2023 / Published: 27 June 2023
(This article belongs to the Special Issue Avian Viral Immunosuppressive Disease)

Abstract

:
The novel variant IBDV (nVarIBDV, especially genotype A2dB1) mainly affects broilers in China. It causes an infection characterized by the atrophy of the bursa, a decrease in the level of lymphocytes, proliferation of fibrous tissue around the follicle, and severe atrophy of the follicle in the bursa. Poultry vaccinated with live IBDV vaccines do not have the challenge present with bursa atrophy, which is misdiagnosed for nVarIBDV because of the lack of other gross clinical symptoms. The present study sought to explore the potential and reliability of the real-time TaqMan analysis method for the detection and discrimination of the nVarIBDV genotype from that of the non-nVarIBDV, especially in live vaccine strains. This method will help monitor vaccinated poultry to control and manage infection with the nVarIBDV IBDVs. The nucleotide polymorphism in the 5′-UTR region and the vp5/vp2 overlapping region of the segment A sequences of IBDV were used to establish a one-step real-time TaqMan reverse transcription polymerase chain reaction (RT-PCR) method in this study. The results showed that the method accurately distinguished the nVarIBDV and non-nVarIBDV strains (especially live vaccine strains), and there were no cross-reactions with the infectious bronchitis virus (IBV), Newcastle disease virus (NDV), avian influenza virus (AIV), infectious laryngotracheitis virus (ILTV), fowlpox virus (FPV), Mycoplasma gallisepticum (M. gallisepticum), Mycoplasma synoviae (M. synoviae), and IBDV-negative field samples. The method showed a linear dynamic range between 102 and 107 DNA copies/reaction, with an average R2 of 0.99 and an efficiency of 93% for nVarIBDV and an average R2 of 1.00 and an efficiency of 94% for non-nVarIBDV. The method was also used for the detection of 84 clinical bursae of chickens vaccinated with the live vaccine. The results showed that this method accurately distinguished the nVarIBDV and non-nVarIBDV strains (vaccine strains), compared with a strategy based on the sequence analysis of HVRs at the vp2 gene or the reverse transcription PCR (RT-PCR) for the vp5 gene. These findings showed that this one-step real-time TaqMan RT-PCR method provides a rapid, sensitive, specific, and simple approach for detection of infections caused by nVarIBDV and is a useful clinical diagnostic tool for identifying and distinguishing nVarIBDV from non-nVarIBDV, especially live vaccine strains.

1. Introduction

Infectious bursal disease (IBD) is a common viral disease that affects poultry worldwide. Infectious bursal disease virus (IBDV) has a bi-segmented RNA genome (segments A and B) located in a nonenveloped icosahedral capsid. The virus belongs to the Avibirnavirus genus in the Birnaviridae family and was initially reported in Gumboro, USA in 1957 [1]. Studies on neuralization and cross-protection experiments report two serotypes of IBDV. Serotype 1 strains cause infections in chicken. These strains are classified as attenuated strains, classical strains (also known as standard strains), variant strains of USA, very virulent (vv) strains [2,3], and novel variant IBDV (nVarIBDV) [4,5,6], according to their pathogenicity. Serotype 2 strains are nonpathogenic to poultry [6]. Segment A consists of two overlapping open reading frames (ORF). The ORF1 encodes viral non-structural proteins (VP5), while ORF2 encodes a polyprotein (VP2-VP4-VP3) that can be cleaved through autoproteolysis to produce VP2, VP3, and VP4 [7]. VP2 forms the outer surface of the virion and serves as a major immune activation protein that can induce host neutralization antibodies. The vp2 gene has a hypervariable region (HVR, nt 616–1050), which is responsible for the viral antigenic variation and virulence [8,9]. Genomic segment B (2.9 kbp) encodes the RNA-dependent RNA polymerase VP1 [10].
The classic and vvIBDV strains are associated with different levels of mortality, with vvIBDV strains characterized by higher morbidity and mortality than classic strains. vvIBDV strains cause a mortality rate of 40–100% in specific pathogen-free (SPF) chickens, up to 60% in layers, and up to 30% in broilers [11,12,13]. Additionally, the vvIBDV was found in the natural infection of turkey poults for the first time in Egypt [14]. Major symptoms of the disease include muscular hemorrhage, inflammatory exudation, hemorrhage, and yellow staining of the bursae of poultry infected with classic and vvIBDV strains [5,13,15]. The nVarIBDV strains (such as FJ2019-01 and SHG19) are novel pathogenic viruses and are significantly different from the American IBDV variants. These viruses are characterized by the atrophy of the bursa, a decrease in the levels of lymphocytes, macrophage infiltration in the follicle, proliferation of fibrous tissue around the follicle, and severe atrophy of the follicle in SPF chickens [4,5]. Notably, no gross clinical symptoms or mortality were observed in the chicken infected with SHG19 or FJ2019-01 IBDV variants [4,5]. Currently used commercial IBDV vaccines are not effective against infection by these nVarIBDVs [5,16]. The nVarIBDVs cause subclinical infections that increase the susceptibility to infection by other pathogens and induce a poor immune response to vaccines [4].
Vaccination is the most effective control of IBD. However, studies report that poultry vaccinated with immune complex vaccines or live vaccines before infection with the virus presented with atrophy of the bursae, compared with that of non-vaccinated birds [17], and had an increased susceptibility to other pathogens [18]. The nVarIBDVs are not effective for differentiation between infected and vaccinated animals (DIVA). Therefore, it is imperative to explore a strategy for the differentiation of infected animals from vaccinated animals for IBDV for effective monitoring of vaccinated flocks to control infections with the nVarIBDVs. The aim of this study was to evaluate the potential and reliability of one-step real-time TaqMan analysis of the 5′-UTR region and the vp5/vp2 overlapping region of the segment A sequences of IBDV to detect and discriminate the genotype of the nVarIBDV from that of the non-nVarIBDV using allelic discrimination probes, especially live vaccine strains. The findings from the study will provide a basis for improving the current diagnostic capability, and the approach is a rapid, sensitive, and specific method for the effective screening of a large number of samples and for distinguishing the vaccine strains that cause bursae atrophy.

2. Materials and Methods

2.1. Virus Strains

The IBDV FJ2019-01 strain (GenBank: MZ736578 or MZ044944) was stored in our laboratory [5]. The IBDV BC6/85 strain (GenBank: ON286951) was maintained in our laboratory and used to assess the effectiveness of licensed IBDV vaccines currently in use in China [5]. The licensed IBDV live vaccine strains, B87, D78, W2512, M.B., K85, NF8, and CF, and BC6/85, YM (vvIBDV), WH (vvIBDV), FJ2019-01, FJ2019-02 (GenBank: MZ044945), FJ2019-03 (GenBank: MZ044946), FJ2019-04 (GenBank: MZ044947), FJ2019-05 (GenBank: MZ044948), and FJ2021(GenBank: MZ593902) from China were used to evaluate the specificity of the assay. RNA or DNA were extracted from the Newcastle disease virus (NDV), avian influenza virus (AIV), infectious laryngotracheitis virus (ILTV), fowlpox virus (FPV), Mycoplasma gallisepticum (M. gallisepticum), Mycoplasma synoviae (M. synoviae), and infectious bronchitis virus (IBV) to carry out specificity tests.

2.2. Primer and Probe Designs

Complete and partial sequences of IBDV genomes (approximately 80 sequences) were retrieved from the GenBank database. Multiple sequence alignments were carried out using MEGA7 software. A nucleotide polymorphism that distinguished the nVarIBDV genogroup from other IBDV strains (including the classic, vvIBDV, variants of USA, and live vaccine strains) was explored. Primer sets and TaqMan minor groove-binding (MGB) probes were designed, based on the 5′-UTR region and the vp5/vp2 overlapping region of the segment A sequences, according to highly conserved regions observed in the multiple sequence alignment (Figure 1), for subsequent use in real-time qPCR. The primer sequences F: 5′-CCT CCT TCT AYA RYG CTR TCA T-3′ and R: 5′-CGT ATG AAC GGA ACA ATC TG-3′ were used to target the SNP-containing region and to amplify a 105 bp sequence (Figure 1). The A/G change maximized differences in Tm between allele-specific probes. Therefore, the probe sequences used in the study were 5′-TAG AGA TCA GAC GAA CG-3′, corresponding to nVarIBDV labeled with the VIC and MGB quencher groups in the 5′ and 3′ termini (named PnV), and 5′-AGT AGA GAT CAG ACA AA-3′, corresponding to non-nVarIBDV labeled with the FAM and MGB quencher groups in the 5′ and 3′ termini (named PnnV), respectively. The primers and probes were synthesized by Sangon Biotech Co., Ltd. (Shanghai, China).

2.3. Viral RNA Extraction and Reverse Transcription

Viral RNA extraction from samples was performed using the QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. RNA from vaccines was extracted using PBS resuspension of lyophilized vaccine powder. Reverse transcription (RT) was performed with the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA).

2.4. Construction of a Positive Plasmid Standard

Partial segment A of the FJ2019-01 strain and the BC6/85 strain was amplified using the RevertAid RT kit (Thermo Fisher Scientific) with specific primers (10 μM; AU [5′-GGA TAC GAT CGG TCT GAC CCC GGG GGA GTC-3′] and A1542L [5′-GTA GTC TAC ACC TTC CCC AAT TGC AT-3′]) [19]. The purified PCR product was cloned into the pCRTMⅡ-Blunt-TOPO® vector using the Zero Blunt® TOPO® PCR cloning kit (Invitrogen, Carlsbad, CA, USA), and the ligation product was used for the transformation of E. coli, DH5α (Tiangen, Beijing, China). The plasmid DNA was extracted using the Qiagen plasmid Mini kit (Qiagen, Hilden, Germany), and the concentrations of plasmids were determined using the NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). The target copy number of positive plasmids was calculated using the formula: copies/μL = (concentration in ng × 6.023 × 1023)/(genome length × 1 × 109 × 660 dalton/bp). The pCR-BC6/85-A and pCR-FJ2019-01-A were constructed and used as positive controls for the detection of IBDV and for performing detection limit assay.

2.5. Establishment of the One-Step Real-Time TaqMan RT-PCR Method

The one-step real-time TaqMan RT-PCR method was performed on a LightCycler96 (Roche, Basel, Switzerland) to distinguish the nVarIBDV genogroup from the non-nVarIBDV strains (including the classic, vvIBDV, variants of USA, and live vaccine strains). The 20 μL reaction mix contained 5 μL Fast one-step Master Mix (4×) (Thermo Fisher Scientific), 0.50 μL forward or reverse primer (10 μM), 0.25 μL PnV probe (10 μM) or the PnnV probe (10 μM), 5 μL template RNA or plasmid DNA, and 8.5 μL deionized distilled water. The reaction conditions were as follows: 50 °C for 10 min; 95 °C for 20 s; 40 cycles of 95 °C for 3 s; and 60 °C for 30 s. Fluorescence intensities for specific reporter fluorophores were determined at the 60 °C step of each cycle and at the end of the run. All reactions and amplifications were performed and analyzed using the LightCycler® 96 SW 1.1 software.

2.6. Determination of Specificity of Primers and Probe Sets

The specificities of the primer pairs and probes were performed using NCBI BLAST https://www.ncbi.nlm.nih.gov/tools/primer-blast/ (accessed on 18 September 2021) (Bethesda, MD, USA). A BLAST search was performed to predict the in silico primer and probe sequence specificities and to evaluate the occurrence of non-specific homology between the sequences and the IBDV genome or the chicken genome. Previously confirmed IBDV-negative samples were also analyzed. The specificity of the established method was verified using RNA samples isolated from B87, D78, W2512, M.B, K85, NF8, CF, BC6/85, FJ2019-01, FJ2019-02, FJ2019-03, FJ2019-04, FJ2019-05, FJ2021, WH, and YS and RNA or DNA samples from NDV, IBV, AIV, ILTV, FPV, M. gallisepticum, and M. synoviae.

2.7. Assay Detection Limit

Serial dilutions (2 × 106 copies/μL to 2 × 101 copies/μL) of two positive plasmids prepared using deionized water were added in the respective reaction mix to evaluate the sensitivities of primers and probes. A standard curve was generated by plotting the Cq values vs. log10 of 10-fold serial dilutions (107 to 102) of plasmid DNA. Analyses for relative and absolute sensitivities were conducted in triplicate. The standard curve and assay efficiency were established using the LightCycler® 96 SW 1.1 software.
Interference resulting from differences in the proportion of the mixed templates was evaluated by mixing two positive plasmids in various ratios (1:107, 1:106 1:105, 1:104, 1:103, 1:102, 1:0, 107:1, 106:1, 105:1, 104:1, 103:1, 102:1, and 0:1; 1 = 105 copies/reaction). Comparisons of different mixed templates and single templates were performed using the obtained Cq values.

2.8. Clinical Sample Testing

A total of 84 clinical bursal samples were collected from 20- to 40-day-old chickens vaccinated with the live vaccine of W2512, B87, D78, M.B., K85, and CF strains obtained from an infection in 6 chicken flocks reported in 2021 in China. The samples were used to verify the accuracy of the established one-step real-time TaqMan RT-PCR assay. Virus RNA was extracted from bursal homogenates and analyzed using the one-step real-time TaqMan RT-PCR assay. The reverse transcription PCR assays were carried out with an IBDV fluorescent reverse transcription PCR (RT-PCR) detection kit (BIOTECHSY, Beijing, China) to detect the IBDV of the vp5 gene, which is present in all IBDVs, regardless of the genotype, to verify the results. In addition, PCR of HVRs at the vp2 gene was performed using the VP2-F (5′-CCT CAG CTT ACC CAC ATC-3′) and VP2-R (5′-CCT TCC CCA ATT GCA TGG-3′) primers, as described previously [4]. PCR products were sequenced and analyzed using previously reported methods [4,20].

3. Results

3.1. Feasibility of the One-Step Real-Time TaqMan RT-PCR Method

A BLAST analysis was performed for the segment A sequences of several IBDV strains to identify nucleotide polymorphisms that could be used to distinguish between nVarIBDV and non-nVarIBDV strains. Two single nucleotide polymorphisms (SNPs) were detected in this study that could be used to distinguish nCv and non-strains, based on the vp5 region. The C103A and G117A SNPs were selected for their high consistency to differentiate the genogroup from non-strains (including cIBDV, vvIBDV, variants of USA, and live vaccine strains) (Figure 1). This was important because the A/G change maximized Tm differences between allele-specific probes. Two TaqMan-MGB probes were designed to discriminate and separately quantify the nVarIBDV and non-nVarIBDV genotypes. The PnV probe was labeled with the VIC fluorescent dye and was specific for nVarIBDV strains. The PnnV probe was labeled with the 6-FAM fluorescent dye and was specific for non-nVarIBDV strains.
A one-step real-time TaqMan RT-PCR platform was established to simultaneously discriminate and quantify nVarIBDV and non-nVarIBDV genotypes. The results showed that only the positive plasmid of nVarIBDV was recognized on the VIC channel, and the positive plasmid of non-nVarIBDV was recognized on the FAM channel (Figure 2). Notably, no significant cross-reaction was observed between positive plasmids of nVarIBDV and non-nVarIBDV using the non-specific TaqMan-MGB probe (Figure 2). The amplification from nVarIBDV (FJ2019-01) and non-nVarIBDV (BC6/85) was detected only with matching primer and probe combinations, indicating typical “S” amplification curves. The Cq values were all less than 35. No specific amplification was detected from mismatched combinations of probes. Negative controls did not exhibit amplification curves (Figure 2). These findings indicate that the established one-step real-time TaqMan RT-PCR method is highly feasible and suitable for the simultaneous identification of the viral RNA of nVarIBDV and non-nVarIBDV.

3.2. Specificity of the One-Step Real-Time TaqMan RT-PCR Method

A nucleotide BLAST search of each primer and probe only showed homology with the genome regions of the expected IBDV genogroup strain. The IBDV-negative field samples and non-IBDV pathogenies (IBV, NDV, AIV, ILTV, FPV, M. gallisepticum, and M. synoviae) were not amplified, and the negative or blank controls in the assay were also not amplified (Figure 3a). Six nVarIBDV strains, two vvIBDV strains, and seven vaccine strains were analyzed to evaluate the assay performance and diagnostic effectiveness (Figure 3b,c). All strains were accurately diagnosed using this method, indicating 100% sensitivity and effectiveness of the method. This result indicates that the two probes were highly specific, feasible, and suitable for targeting the IBDV genogroup, and no cross-reactivity was observed between the nVarIBDV and non-nVarIBDV strains.

3.3. Sensitivity of the One-Step Real-Time TaqMan RT-PCR Method

The one-step real-time TaqMan RT-PCR method was performed using a series of known concentrations of plasmid standards to explore its sensitivity. The assay using the pCR-FJ2019-01-A plasmid as a template showed a linear dynamic range between 102 and 107 DNA copies/reaction, with an average R2 of 0.99 and an efficiency of 93% (Figure 4a,c). The assay using the pCR-BC6/85-A plasmid as a template exhibited a linear dynamic range between 102 and 107 DNA copies/reaction, with an average R2 of 1.00 and an efficiency of 94% (Figure 4b,c). Additionally, the assay was performed using different ratios of the nVarIBDV template to non-nVarIBDV template. The Cq values obtained from the assay are presented in Table 1. Compared with single templates, low-copy (less 103 copies/reaction) templates did not generate the amplified FAM signal or VIC signal when the mixed templates had another template with high copies. This finding implied a presence of competition for low copies of templates.

3.4. Analysis of the Clinical Samples Using the One-Step Real-Time TaqMan RT-PCR Method

A total of 84 clinical bursae samples were collected from six chicken flocks vaccinated with W2512, B87, D78, M.B., K85, and CF strains. The results showed that 84 samples were positive for the one-step real-time TaqMan RT-PCR method (84/84, 100%) (Table 2). This implies that the established one-step real-time TaqMan RT-PCR method can be used for nVarIBDV and non-nVarIBDV strain detection. A total of 19 samples out of the 84 clinical samples were positive for nVarIBDV (22.62%), 67 samples were positive for non-nVarIBDV (79.76%), and 2 samples were positive for both IBDV (2.38%) (Table 2). Then, RT-PCR was conducted for the vp5 gene to detect IBDV, and the results were positive for all 84 clinical bursae samples. Sequence analysis of HVRs at the vp2 gene was conducted to explore the genotyping capacity of the method. Some IBDV-positive samples (15 of 84) failed the sequence analysis of the HVRs process, due to a low number of copies of the templates. The results showed that IBDV-positive samples that were successfully sequenced were correctly detected by the corresponding specific TaqMan-MGB probe to the genotype (Table 2). This finding demonstrated that the established one-step real-time TaqMan RT-PCR method has a potential genotyping capacity and high effectiveness. The specificity and sensitivity of the one-step real-time TaqMan RT-PCR method were 100% and 100%, respectively, compared with the traditional sequence analysis of HVRs at the vp2 gene (Table 2). The genotype of IBDV is important in distinguishing between the vaccine strains with atrophied bursae and nVarIBDV and can be quickly explored using the one-step real-time TaqMan RT-PCR method, compared with traditional RT-PCR and sequence analysis.

4. Discussion

IBD induced by serotype 1 IBDV is a common infectious disease reported in the chicken industry globally and was also found in turkey poults for the first time in Egypt [14]. Serotype 1 strains are classified as avirulent strains, classical strains (also known as standard strains), variant strains of USA, and very virulent (vv) strains [2,3], according to their pathogenicity. The nVarIBDV strains belonging to the A2dB1 genotype were reported in 2019 in China and have rapidly spread to several provinces in China [4,5,6]. The new reassortment strains (genotype A2dB3) have also been observed in China, with a segment A from nVarIBDV strains and a segment B from HLJ0504-like strains (Genotype A3B3), which shows a similar pathogenicity to nVarIBDV in specific-pathogen-free (SPF) chickens [21]. Currently, various methods are available for the diagnosis of IBDV, such as virus isolation, RT-PCR, real-time RT-PCR, and ELISA. Real-time RT-PCR methods are highly sensitive and specific and mainly require “SYBR Green” or “TaqMan” probes to generate high-resolution melting curves or specific fluorescent signals. The real-time RT-PCR method is currently one of the most promising and effective methods used in control and epidemiological surveillance programs [22,23,24,25,26]. In addition, genetic relatedness is currently widely used to characterize IBDV strains [22,23,25,27,28]. Conventional reverse transcription polymerase chain reactions (RT-PCR) and diverse genotyping techniques, such as nucleotide sequencing and restriction fragment length polymorphism (RFLP) [20,22,29,30], have been replaced with the RT-PCR method using two distinct and specific probes or high-resolution melting curve analysis. This novel RT-PCR method can simultaneously and accurately discriminate differentiations between infectious bursal disease virus strains [22,25,26]; thus, it is more time-saving and economical, compared with conventional methods. Therefore, simpler and more rapid assays for detection and discrimination of all nVarIBDV from other IBDV genotype strains, especially vaccine strains, should be explored to monitor strain spread and improve the control of this disease. Differentiation of vaccines and wild strains is crucial for developing effective control strategies.
Monitoring of nVarIBDV is mainly conducted through sequencing of the full viral genome or sequencing partially conserved regions of vp1 and vp2 genes [4,5]. However, this approach is costly, requires technical expertise, and is characterized by long processing times, which delay control and treatment actions. In this study, an investigation of distinct SNPs at the vp5/vp2 overlapping region of the segment A sequences was conducted to evaluate them as potential targets to rapidly screen and distinguish nVarIBDV from non-nVarIBDV strains, especially live vaccine strains. Specificity of primers and two probes for the one-step real-time TaqMan RT-PCR method was evaluated experimentally, as well as through BLAST analysis. The results demonstrated that these primers and probes had no cross-reactions with virus DNA or RNA and the chicken genome. Moreover, the probes did not exhibit cross-reactions between nVarIBDV and non-nVarIBDV strains, indicating that the method had high specificity. The results using this novel assay showed that the set of primers and probes accurately distinguished the nVarIBDV and non-nVarIBDV strains, and there were no cross-reactions with IBV, NDV, AIV, ILTV, FPV, M. gallisepticum, M. synoviae, and IBDV-negative field samples.
The real-time PCR method has a high detection limit and can be used for IBDV detection or for distinguishing vvIBDV from other non-vvIBDV genotype strains [22,25,26,31]. The method developed in the present study had a lower limit of detection for IBDV identification (102 copies/reaction), compared with the traditional real-time PCR method. The findings showed that the method could detect target templates at a ratio of 105:104 copies/reaction of the two templates to accurately distinguish strains in the interference assay. The high Cq values of the target and no Cq (no amplification) to lower number of target copies in the mixed templates could be attributed to the interference of exponential amplification, due to a higher number of copies of target templates. However, the interference was not evaluated for mixed templates in the TaqMan-MGB real-time RT-PCR assay, which can accurately distinguish all vvIBDV strains from non-vvIBDV strains [22]. The results showed that the method simultaneously determined the identities of nVarIBDV and non-nVarIBDV strains (especially live vaccine strains) in bursa samples, even at a low target template concentration of 102 copies/reaction.
Classic IBDV and vvIBDV infections cause hemorrhage, waxy yellow jelly in the bursa, and/or mortality [5,13,15]. Therefore, it is important to explore the differences in clinical symptoms between classic or vvIBDV strains and live vaccine strains in infected animals using different PCR methods to determine all strains of IBDVs. The nVarIBDV infection causes atrophy of the bursa, induces a decrease in the level of lymphocytes, promotes macrophage infiltration in the follicle, causes proliferation of fibrous tissue around the follicle, and leads to severe atrophy of the follicle in SPF chickens [4,5]. Notably, unchallenged chickens vaccinated with immune complex vaccines present with bursa atrophy with a bursa:body weight index (BBIX = [bursa:body weight ratios]/[bursa:body weight ratios in the negative group]) of 0.59–0.26, compared with those of unvaccinated birds aged 21 to 35 days [17]. Atrophy of the bursa and non-specific gross clinical symptoms minimized the differentiation of the nVarIBDV and live vaccine strains when reverse transcription PCRs (RT-PCR), nanoparticle-assisted PCR, SYBR green, and TaqMan-based real-time RT-PCRs (RT-qPCR) were used to detect all IBDVs [22,25,26,31,32]. The one-step real-time TaqMan RT-PCR method reported in this study effectively and simultaneously distinguished the nVarIBDV from commercially live IBDV vaccine strains and the clinical samples obtained from chickens vaccinated using the live vaccine. The results showed that some bursa samples obtained from chicken were effectively distinguished from the nVarIBDV and vaccine strains. These results imply that the one-step real-time TaqMan RT-PCR method is a highly sensitive, specific, and accurate method for simultaneous and accurate detection of the nVarIBDV in bursa samples obtained from chickens vaccinated with live vaccines when infection with the nVarIBDV is suspected. However, reassortment, gene mutations, and homologous recombination are the main mechanisms in the evolution of double-stranded RNA; the one-step real-time TaqMan RT-PCR method may have deficiencies when detecting the new IBDV variants.

5. Conclusions

This one-step real-time TaqMan RT-PCR method provides a rapid, sensitive, specific, and simple strategy for the detection of infections caused by nVarIBDV and for epidemiological investigation. Additionally, this method is a useful clinical diagnostic tool for identifying and distinguishing nVarIBDV from non-nVarIBDV, especially live vaccine strains.

Author Contributions

Conceptualization, B.H. and C.W. (Chunhe Wan); methodology, B.H. and C.W. (Chenyan Wang); bioinformatics, B.H. and C.W. (Chenyan Wang); resources, G.S.; data curation, C.W. (Chenyan Wang); writing—original draft preparation, B.H.; writing—review and editing, B.H. and C.W. (Chenyan Wang); funding acquisition, B.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Sci-Tech Innovation Team of Fujian Academy of Agricultural Sciences, grant number CXTD2021014-3; and CAAS-FuJian Government 5511 Collaborative Innova-tion Project, Grant No. XTCXGC2021008.

Institutional Review Board Statement

The clinical bursae collected from six chicken flocks were conducted under the guidance of the FAAS’s Institutional Animal Care and Use Committee and were performed in accordance with the animal ethics guidelines and approved protocols. All chickens were euthanized with tiletamine hydrochloride and zolazepam hydrochloride for injection (Virbac, Carros, France), and every effort was made to minimize suffering.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are reported in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eterradossi, N.; Saif, Y.M. Infectious bursal disease. In Diseases of Poultry, 14th ed.; Saif, Y.M., Fadly, A.M., Glison, J.R., McDouglad, L.R., Nolan, L.K., Eds.; Ames IA Blackwell Publishing: Ames, IA, USA, 2017; pp. 185–208. [Google Scholar]
  2. Jackwood, D.J.; Saif, Y.M.; Hughes, J.H. Characteristics and serologic studies of two serotypes of infectious bursal disease virus in turkeys. Avian Dis. 1982, 26, 871–882. [Google Scholar] [CrossRef] [PubMed]
  3. Müller, H.; Islam, M.R.; Raue, R. Research on infectious bursal disease—The past, the present and the future. Vet. Microbiol. 2003, 97, 153–165. [Google Scholar] [CrossRef] [PubMed]
  4. Fan, L.; Wu, T.; Hussain, A.; Gao, Y.; Zeng, X.; Wang, Y.; Gao, L.; Li, K.; Wang, Y.; Liu, C.; et al. Novel variant strains of infectious bursal disease virus isolated in China. Vet. Microbiol. 2019, 230, 212–220. [Google Scholar] [CrossRef] [PubMed]
  5. Hou, B.; Wang, C.; Luo, Z.; Shao, G. Commercial vaccines used in China do not protect against a novel infectious bursal disease virus variant isolated in Fujian. Vet. Rec. 2022, 191, e1840. [Google Scholar] [CrossRef]
  6. Zhang, W.; Wang, X.; Gao, Y.; Qi, X. The Over-40-Years-Epidemic of Infectious Bursal Disease Virus in China. Viruses 2022, 14, 2253. [Google Scholar] [CrossRef]
  7. Petit, S.; Lejal, N.; Huet, J.; Delmas, B. Active Residues and Viral Substrate Cleavage Sites of the Protease of the Birnavirus Infectious Pancreatic Necrosis Virus. J. Virol. 2000, 74, 2057–2066. [Google Scholar] [CrossRef] [Green Version]
  8. Letzel, T.; Coulibaly, F.; Rey, F.A.; Delmas, B.; Jagt, E.; van Loon, A.A.; Mundt, E. Molecular and structural bases for the antigenicity of VP2 of infectious bursal disease virus. J. Virol. 2007, 81, 12827–12835. [Google Scholar] [CrossRef] [Green Version]
  9. Mato, T.; Tatar-Kis, T.; Felfoldi, B.; Jansson, D.S.; Homonnay, Z.; Banyai, K.; Palya, V. Occurrence and spread of a reassortant very virulent genotype of infectious bursal disease virus with altered VP2 amino acid profile and pathogenicity in some European countries. Vet. Microbiol. 2020, 245, 108663. [Google Scholar] [CrossRef]
  10. Spies, U.; Muller, H.; Becht, H. Properties of RNA polymerase activity associated with infectious bursal disease virus and characterization of its reaction products. Virus Res. 1987, 8, 127–140. [Google Scholar] [CrossRef]
  11. Jackwood, D.J.; Sommer-Wagner, S.E.; Stoute, A.S.; Woolcock, P.R.; Crossley, B.M.; Hietala, S.K.; Charlton, B.R. Characteristics of a very virulent infectious bursal disease virus from California. Avian Dis. 2009, 53, 592–600. [Google Scholar] [CrossRef]
  12. van den Berg, T.P.; Morales, D.; Eterradossi, N.; Rivallan, G.; Toquin, D.; Raue, R.; Zierenberg, K.; Zhang, M.F.; Zhu, Y.P.; Wang, C.Q.; et al. Assessment of genetic, antigenic and pathotypic criteria for the characterization of IBDV strains. Avian Pathol. 2004, 33, 470–476. [Google Scholar] [CrossRef] [PubMed]
  13. Van den Berg, T.P.; Gonze, M.; Meulemans, G. Acute infectious bursal disease in poultry: Isolation and characterisation of a highly virulent strain. Avian Pathol. 1991, 20, 133–143. [Google Scholar] [CrossRef]
  14. Mosad, S.M.; Eladl, A.H.; El-Tholoth, M.; Ali, H.S.; Hamed, M.F. Molecular characterization and pathogenicity of very virulent infectious bursal disease virus isolated from naturally infected turkey poults in Egypt. Trop. Anim. Health Prod. 2020, 52, 3819–3831. [Google Scholar] [CrossRef]
  15. Stoute, S.T.; Jackwood, D.J.; Sommer-Wagner, S.E.; Crossley, B.M.; Woolcock, P.R.; Charlton, B.R. Pathogenicity associated with coinfection with very virulent infectious bursal disease and Infectious bursal disease virus strains endemic in the United States. J. Vet. Diagn. Investig. 2013, 25, 352–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Fan, L.; Wu, T.; Wang, Y.; Hussain, A.; Jiang, N.; Gao, L.; Li, K.; Gao, Y.; Liu, C.; Cui, H.; et al. Novel variants of infectious bursal disease virus can severely damage the bursa of fabricius of immunized chickens. Vet. Microbiol. 2020, 240, 108507. [Google Scholar] [CrossRef] [PubMed]
  17. Kelemen, M.; Forgách, K.; Iván, J.; Palya, V.; Süveges, T.; Tóth, B.; Mészáros, J. Pathological and immunological study of an in ovo complex vaccine against infectious bursal disease. Acta Vet. Hung. 2000, 48, 443–454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Arafat, N.; Eladl, A.H.; Mahgoub, H.; El-Shafei, R.A. Effect of infectious bursal disease (IBD) vaccine on Salmonella Enteritidis infected chickens. Vaccine 2017, 35, 3682–3689. [Google Scholar] [CrossRef]
  19. Hussain, A.; Wu, T.; Li, H.; Fan, L.; Li, K.; Gao, L.; Wang, Y.; Gao, Y.; Liu, C.; Cui, H.; et al. Pathogenic Characterization and Full Length Genome Sequence of a Reassortant Infectious Bursal Disease Virus Newly Isolated in Pakistan. Virol. Sin. 2019, 34, 102–105. [Google Scholar] [CrossRef]
  20. Islam, M.R.; Nooruzzaman, M.; Rahman, T.; Mumu, T.T.; Rahman, M.M.; Chowdhury, E.H.; Eterradossi, N.; Muller, H. A unified genotypic classification of infectious bursal disease virus based on both genome segments. Avian Pathol. 2021, 50, 190–206. [Google Scholar] [CrossRef]
  21. Wang, Y.; Jiang, N.; Fan, L.; Niu, X.; Zhang, W.; Huang, M.; Gao, L.; Li, K.; Gao, Y.; Liu, C.; et al. Identification and Pathogenicity Evaluation of a Novel Reassortant Infectious Bursal Disease Virus (Genotype A2dB3). Viruses 2021, 13, 1682. [Google Scholar] [CrossRef]
  22. Tomás, G.; Hernández, M.; Marandino, A.; Panzera, Y.; Maya, L.; Hernández, D.; Pereda, A.; Banda, A.; Villegas, P.; Aguirre, S.; et al. Development and validation of a TaqMan-MGB real-time RT-PCR assay for simultaneous detection and characterization of infectious bursal disease virus. J. Virol. Methods 2012, 185, 101–107. [Google Scholar] [CrossRef]
  23. Wu, C.C.; Rubinelli, P.; Lin, T.L. Molecular detection and differentiation of infectious bursal disease virus. Avian Dis. 2007, 51, 515–526. [Google Scholar] [CrossRef] [PubMed]
  24. Li, Y.P.; Handberg, K.J.; Kabell, S.; Kusk, M.; Zhang, M.F.; Jørgensen, P.H. Relative quantification and detection of different types of infectious bursal disease virus in bursa of Fabricius and cloacal swabs using real time RT-PCR SYBR green technology. Res. Vet. Sci. 2007, 82, 126–133. [Google Scholar] [CrossRef] [PubMed]
  25. Kong, L.L.; Omar, A.R.; Hair Bejo, M.; Ideris, A.; Tan, S.W. Development of SYBR green I based one-step real-time RT-PCR assay for the detection and differentiation of very virulent and classical strains of infectious bursal disease virus. J. Virol. Methods 2009, 161, 271–279. [Google Scholar] [CrossRef] [PubMed]
  26. Ghorashi, S.A.; O’Rourke, D.; Ignjatovic, J.; Noormohammadi, A.H. Differentiation of infectious bursal disease virus strains using real-time RT-PCR and high resolution melt curve analysis. J. Virol. Methods 2011, 171, 264–271. [Google Scholar] [CrossRef]
  27. Hernandez, M.; Banda, A.; Hernandez, D.; Panzera, F.; Perez, R. Detection of very virulent strains of infectious bursal disease virus (vvIBDV) in commercial broilers from Uruguay. Avian Dis. 2006, 50, 624–631. [Google Scholar] [CrossRef]
  28. Kusk, M.; Kabell, S.; Jørgensen, P.H.; Handberg, K.J. Differentiation of five strains of infectious bursal disease virus: Development of a strain-specific multiplex PCR. Vet. Microbiol. 2005, 109, 159–167. [Google Scholar] [CrossRef]
  29. Hernandez, M.; Tomas, G.; Hernandez, D.; Villegas, P.; Banda, A.; Maya, L.; Panzera, Y.; Perez, R. Novel multiplex RT-PCR/RFLP diagnostic test to differentiate low- from high-pathogenic strains and to detect reassortant infectious bursal disease virus. Avian Dis. 2011, 55, 368–374. [Google Scholar] [CrossRef]
  30. Drissi, T.C.; Fellahi, S.; Fassi, F.O.; Gaboun, F.; Khayi, S.; Mentag, R.; Lico, C.; Baschieri, S.; El, H.M.; Ducatez, M. Complete genome analysis and time scale evolution of very virulent infectious bursal disease viruses isolated from recent outbreaks in Morocco. Infect. Genet. Evol. 2020, 77, 104097. [Google Scholar] [CrossRef]
  31. Techera, C.; Tomás, G.; Panzera, Y.; Banda, A.; Perbolianachis, P.; Pérez, R.; Marandino, A. Development of real-time PCR assays for single and simultaneous detection of infectious bursal disease virus and chicken anemia virus. Mol. Cell. Probes 2019, 43, 58–63. [Google Scholar] [CrossRef]
  32. Luan, Q.; Jiang, Z.; Wang, D.; Wang, S.; Yin, Y.; Wang, J. A sensitive triple nanoparticle-assisted PCR assay for detection of fowl adenovirus, infectious bursal disease virus and chicken anemia virus. J. Virol. Methods 2022, 303, 114499. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Alignment of the TaqMan-minor groove-binding (MGB) real-time reverse transcription polymerase chain reaction (RT-PCR)-amplified region. Various representative nVarIBDV and non-nVarIBDV isolates were included in the alignment. Forward and reverse primers and the probe target site are indicated using shadow position of the reference strain FJ2021 (MZ593902) in the alignment.
Figure 1. Alignment of the TaqMan-minor groove-binding (MGB) real-time reverse transcription polymerase chain reaction (RT-PCR)-amplified region. Various representative nVarIBDV and non-nVarIBDV isolates were included in the alignment. Forward and reverse primers and the probe target site are indicated using shadow position of the reference strain FJ2021 (MZ593902) in the alignment.
Viruses 15 01453 g001
Figure 2. The feasibility of the one-step real-time TaqMan RT-PCR method. IBDV FJ2019-01 strain (a), positive plasmid pCR-FJ2019-01-A (b), BC6/85 strain (c), and positive plasmid pCR-BC6/85-A (d) specifically exhibited the amplified VIC signal (dotted line) and the FAM signal (full line). The IBDV-negative field sample (e), AVE buffer (f), and H2O (g) did not exhibit amplified signals during the assay.
Figure 2. The feasibility of the one-step real-time TaqMan RT-PCR method. IBDV FJ2019-01 strain (a), positive plasmid pCR-FJ2019-01-A (b), BC6/85 strain (c), and positive plasmid pCR-BC6/85-A (d) specifically exhibited the amplified VIC signal (dotted line) and the FAM signal (full line). The IBDV-negative field sample (e), AVE buffer (f), and H2O (g) did not exhibit amplified signals during the assay.
Viruses 15 01453 g002
Figure 3. Specificity of the one-step real-time TaqMan RT-PCR assay. (a) FJ2021 and B87 strains were the positive controls for VIC and FAM signals, respectively; IBV, NDV, AIV, ILTV, FPV, M. gallisepticum, M. synoviae, and IBDV-negative field samples did not show amplification signals for VIC or FAM during the assay. (b) Six nVarIBDV strains, BC6/85, and seven vaccine strains were analyzed to specifically generate the amplified VIC signal and FAM signal, respectively. (c) Two vvIBDV strains (WH and YS) were analyzed to specifically generate the amplified FAM signal.
Figure 3. Specificity of the one-step real-time TaqMan RT-PCR assay. (a) FJ2021 and B87 strains were the positive controls for VIC and FAM signals, respectively; IBV, NDV, AIV, ILTV, FPV, M. gallisepticum, M. synoviae, and IBDV-negative field samples did not show amplification signals for VIC or FAM during the assay. (b) Six nVarIBDV strains, BC6/85, and seven vaccine strains were analyzed to specifically generate the amplified VIC signal and FAM signal, respectively. (c) Two vvIBDV strains (WH and YS) were analyzed to specifically generate the amplified FAM signal.
Viruses 15 01453 g003
Figure 4. Sensitivity and standard curves of the one-step real-time TaqMan RT-PCR assay. (a,b) Amplification curves using 10-fold serially diluted template concentrations of pCR-FJ2019-01-A and pCR-BC6/85-A plasmids. (c) The linear dynamic range was established between 102 and 107 copies/reaction for nVarIBDV (dotted line) and non-nVarIBDV (full line) templates. The coefficient of determination (R2) and efficiency of each linear regression curve are indicated.
Figure 4. Sensitivity and standard curves of the one-step real-time TaqMan RT-PCR assay. (a,b) Amplification curves using 10-fold serially diluted template concentrations of pCR-FJ2019-01-A and pCR-BC6/85-A plasmids. (c) The linear dynamic range was established between 102 and 107 copies/reaction for nVarIBDV (dotted line) and non-nVarIBDV (full line) templates. The coefficient of determination (R2) and efficiency of each linear regression curve are indicated.
Viruses 15 01453 g004
Table 1. Mean Cq values of different concentrations of the ratios of templates for nVarIBDV to non-nVarIBDV, as assessed using the one-step real-time TaqMan RT-PCR method.
Table 1. Mean Cq values of different concentrations of the ratios of templates for nVarIBDV to non-nVarIBDV, as assessed using the one-step real-time TaqMan RT-PCR method.
Ratios of nVarIBDV to Non-nVarIBDV Templates (Copies/Reaction)Cq FAMCq VIC
105:10714.0718.94
105:10618.5321.04
105:10520.5621.01
105:10423.3320.87
105:103Negative20.50
105:102Negative21.07
0:10714.06Negative
0:10616.49Negative
0:10520.44Negative
0:10423.41Negative
0:10327.60Negative
0:10231.04Negative
107:10526.7714.08
106:10521.6518.07
105:10521.1321.47
104:10520.9824.33
103:10521.57Negative
102:10520.81Negative
107:0Negative14.58
106:0Negative16.78
105:0Negative20.32
104:0Negative23.53
103:0Negative27.66
102:0Negative31.95
Table 2. Results of the clinical sample analyzed using the one-step real-time TaqMan RT-PCR method, the analysis of HVRs at the vp2 gene, or reverse transcription PCR (RT-PCR) for the vp5 gene.
Table 2. Results of the clinical sample analyzed using the one-step real-time TaqMan RT-PCR method, the analysis of HVRs at the vp2 gene, or reverse transcription PCR (RT-PCR) for the vp5 gene.
Sample No.One-Step Real-Time TaqMan RT-PCRRT-PCR for vp5 GeneHVRs of vp2 Gene
FAM CqVIC Cq
S1negative20.23positivenVarIBDV
S228.39negativepositivenon-nVarIBDV
S329.95negativepositiveN.D.
S432.25negativepositiveN.D.
S527.83negativepositivenon-nVarIBDV
S625.44negativepositivenon-nVarIBDV
S729.95negativepositiveN.D.
S824.01negativepositivenon-nVarIBDV
S926.63negativepositivenon-nVarIBDV
S1027.22negativepositivenon-nVarIBDV
S1125.54negativepositivenon-nVarIBDV
S1228.67negativepositiveN.D.
S13negative22.28positivenVarIBDV
S1426.62negativepositivenon-nVarIBDV
S1529.11negativepositiveN.D.
S1627.55negativepositivenon-nVarIBDV
S1726.71negativepositivenon-nVarIBDV
S1828.95negativepositiveN.D.
S1930.13negativepositiveN.D.
S2025.46negativepositivenon-nVarIBDV
S2127.02negativepositivenon-nVarIBDV
S2225.18negativepositivenon-nVarIBDV
S2322.93negativepositivenon-nVarIBDV
S2427.36negativepositivenon-nVarIBDV
S25negative23.39positivenVarIBDV
S2626.49negativepositivenon-nVarIBDV
S2728.16negativepositivenon-nVarIBDV
S2825.77negativepositivenon-nVarIBDV
S2926.32negativepositivenon-nVarIBDV
S3029.15negativepositiveN.D.
S3127.0329.80positivenon-nVarIBDV
S3221.43negativepositivenon-nVarIBDV
S3325.75negativepositivenon-nVarIBDV
S3423.13negativepositivenon-nVarIBDV
S35negative22.59positivenVarIBDV
S3629.10negativepositiveN.D.
S3726.70negativepositivenon-nVarIBDV
S3830.75negativepositiveN.D.
S3926.60negativepositivenon-nVarIBDV
S4029.02negativepositiveN.D.
S4123.34negativepositivenon-nVarIBDV
S4229.06negativepositiveN.D.
S4322.77negativepositivenon-nVarIBDV
S4426.68negativepositivenon-nVarIBDV
S45negative26.02positivenVarIBDV
S4626.33negativepositivenon-nVarIBDV
S4727.21negativepositivenon-nVarIBDV
S4825.83negativepositivenon-nVarIBDV
S4925.96negativepositivenon-nVarIBDV
S50negative15.19positivenVarIBDV
S51negative14.33positivenVarIBDV
S5222.24negativepositivenon-nVarIBDV
S53negative11.37positivenVarIBDV
S54negative21.53positivenVarIBDV
S55negative18.45positivenVarIBDV
S5627.71negativepositivenon-nVarIBDV
S5732.1429.23positiveN.D.
S58negative16.56positivenVarIBDV
S5924.82negativepositivenon-nVarIBDV
S60negative21.17positivenVarIBDV
S6123.13negativepositivenon-nVarIBDV
S6225.93negativepositivenon-nVarIBDV
S6326.15negativepositivenon-nVarIBDV
S6428.01negativepositivenon-nVarIBDV
S65negative25.05positivenVarIBDV
S6627.97negativepositivenon-nVarIBDV
S6727.44negativepositivenon-nVarIBDV
S6831.07negativepositiveN.D.
S6924.90negativepositivenon-nVarIBDV
S70negative23.50positivenVarIBDV
S7126.68negativepositivenon-nVarIBDV
S7226.96negativepositivenon-nVarIBDV
S7325.34negativepositivenon-nVarIBDV
S7426.65negativepositivenon-nVarIBDV
S7526.49negativepositivenon-nVarIBDV
S7629.37negativepositiveN.D.
S7725.38negativepositivenon-nVarIBDV
S7828.48negativepositiveN.D.
S79negative18.21positivenVarIBDV
S80negative18.58positivenVarIBDV
S8127.61negativepositivenon-nVarIBDV
S8225.41negativepositivenon-nVarIBDV
S8326.81negativepositivenon-nVarIBDV
S84negative21.48positivenVarIBDV
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, C.; Hou, B.; Shao, G.; Wan, C. Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China. Viruses 2023, 15, 1453. https://doi.org/10.3390/v15071453

AMA Style

Wang C, Hou B, Shao G, Wan C. Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China. Viruses. 2023; 15(7):1453. https://doi.org/10.3390/v15071453

Chicago/Turabian Style

Wang, Chenyan, Bo Hou, Guoqing Shao, and Chunhe Wan. 2023. "Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China" Viruses 15, no. 7: 1453. https://doi.org/10.3390/v15071453

APA Style

Wang, C., Hou, B., Shao, G., & Wan, C. (2023). Development of a One-Step Real-Time TaqMan Reverse Transcription Polymerase Chain Reaction (RT-PCR) Assay for the Detection of the Novel Variant Infectious Bursal Disease Virus (nVarIBDV) Circulating in China. Viruses, 15(7), 1453. https://doi.org/10.3390/v15071453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop