Ethics statement
This study was conducted according to the animal welfare guidelines of the World Organization for Animal Health and was approved by the Animal Welfare Committee of the China Animal Health and Epidemiology Center (CAHEC). CAHEC has permission to study the activities of HPAIV. The swab samples were collected after being granted permission by several relevant parties, including the Ministry of Agriculture and Rural Affairs of China, CAHEC, the relevant veterinary sections of the provincial and national governments, and the relevant farm owners.
Samples and extraction of viral nucleic acids
All the nucleic acid of AIVs used in this study were identified in our laboratory (the National Avian Influenza Professional Laboratory in China Animal Health and Epidemiology Center) and were stored at -80°C. Newcastle disease virus (NDV) and infectious bronchitis virus (IBV) were maintained in our laboratory. The samples were centrifuged at 12,000 × g for 10 min and the supernatant from each sample was used for RNA extraction on the QIAxtractor platform using a cador Pathogen 96 QIAcube HT kit (Qiagen). The extracted RNA was stored at -80°C for subsequent tests.
Preparation of viral RNA standard
An H5 subtype AIV viral RNA standard was developed using the reference strain A/duck/Yunnan/5310/2006(H5N1) (GenBank accession number CY030889). A 1776-bp fragment of the whole HA gene of H5 subtype AIV was constructed and digested with enzyme. The digested plasmid was purified with DNA recovery kit (Thermo Fisher Scientific) and transcribed in vitro with RNA Production System-T7 kit. The products were extracted and purified by Rneasy MiNi Kit to remove the heteroproteins and various ions in the system. The concentration of cRNA standard transcribed in vitro was accurately measured by Nano Drop nucleic acid quantitative analyzer, and the copy number was calculated. The cRNA standard was continuously diluted with RNase free H2O in 10-fold gradient. The sensitivity of the RT-RAA was evaluated by real-time fluorescence detection, using the diluted cRNA standard ranging from 106–101 copies/μL.
Design of H5 RT-RAA primers and exo-probes
To detect H5 subtype AIV, a total of 4636 available HA gene sequences of H5 subtype AIV obtained from GenBank database were aligned, which contained the HA gene sequence of all currently circulating branches of H5 subtype AIV, and highly conserved regions were subsequently identified with Molecular Evolutionary Genetics Analysis (MEGA) software 6.0 [11] for the design the gene-specific primers and probes. Primers were designed using OLIGO 7 software [12] and showed no major nonspecific sequence similarities by BLAST analysis. Three H5 forward primers and five reverse primers were designed to select the best primers and probes in combination. The appropriate primers and probe combination were selected by sequence analysis, which were shown in Table 1. The 30th base at the 5′ end of the probe was labeled with the FAM fluorophore. The 30th base was connected to the abasic site tetrahydrofuran (THF). The 31st base was labeled with the BHQ1 quencher, and the 3′ end was modified by 3ʹ block. All the primers and probes were synthesized by Sangon Biotech.
Table 1
Primer and probe sequences used for RT-RAA and RT-qPCR assays
Primer
|
Squence (5'-3')
|
Size (bp)
|
Gene
|
Source
|
H5- F
|
CAGTTTGAGGCYGTTGGAAGGGAATTTAAYAA
|
32
|
HA
|
This study
|
H5- R
|
CTTGTCRTAAAGGTTCTTGACATTTGAGTCAT
|
32
|
HA
|
This study
|
H5- P
|
CTAGATGTCTGGACTTATAATGCTGAACT/i6FAMdT/ /THF/ /iBHQ1dT/GGTTCTCATGGAAAAT[C3-spacer]
|
47
|
HA
|
This study
|
H5+1456
|
ACGTATGACTATCCACAATACTCAG
|
25
|
HA
|
[13]
|
H5-1685
|
AGACCAGCTACCATGATTGC
|
20
|
HA
|
[13]
|
H5+1637
|
FAM-TCAACAGTGGCGAGTTCCCTAGCA-TAMRA
|
24
|
HA
|
[13]
|
RT-RAA for detection of H5 subtype AIV
The appropriate primers and exo-probes were screened, and the RT-RAA reaction was performed with an RT exo kit in 50 μL reaction mixture including the necessary enzymes and reagents for RT and DNA amplification in lyophilized pellets (Jiangsu Qitian Bio-Tech Co. Ltd.). The reaction mixture contained the following: 2 μL RNA template, 25 μL rehydration buffer, 15.7 μL deionized distilled water, 2.5 μL magnesium acetate, 2.1 μL each primer (10 μM) and 0.6 μL target-specific RT-RAA exo-probe (the probe concentration is 40 ng/μL). For amplification, the tubes were then transferred to a tube holder in an RT-RAA fluorescence detection device (QT-RAA-F7200; Jiangsu Qitian Bio-Tech Co. Ltd.) set at 39°C for 20 min. Each run included nuclease-free water as a negative control. The results of the RT-RAA interpretation criteria were determined by the slope of amplification curve of the RT-RAA fluorescence detection device (QT-RAA-F7200), the test result is positive when the slope k value was greater than or equal to 20, in addition, the test result is negative.
Analytical sensitivity and analytical specificity of RT-RAA
A dilution range of H5 subtype HA plasmid standard was used to select the appropriate RAA primers and exo-probe combination. The selected appropriate primers and exo-probes were verified by H5 RT-RAA analysis. To determine the analytical sensitivity of the H5 RT-RAA assay, we detected the H5 subtype AIV molecular RNA standard over a dilution range of 107–101 copies/μL under the optimal RT-RAA conditions, with eight replicates for each dilution.
The analytical specificity assay of the RT-RAA for H5 subtype AIV was evaluated using four H5-positive AIVs (one strain in clade 2.3.2.1, one in clade 7, and 55 in clade 2.3.4.4), 10 other subtype AIVs (H1N2, H3N2, H4N2, H6N2, H7N3, H7N9, H9N2, H10N7, H11N9), two NDVs and two IBVs. These viruses are the main respiratory viruses affecting birds and were previously identified by our laboratory. The details of all the viruses tested are listed in Table 2.
Table 2
Some information about the samples used for analytical sensitivity and analytical specificity assay in the study
Sample
|
Virus
|
HA subtype(clade)
|
H5 RT-RAA assay
|
H5 RT-qPCR assay
|
K144
|
AIV
|
H5N1(2.3.2.1)
|
+
|
+
|
QD1
|
AIV
|
H5N2(7)
|
+
|
+
|
G2324
|
AIV
|
H5N6(2.3.4.4)
|
+
|
+
|
G2084
|
AIV
|
H5N6(2.3.4.4)
|
+
|
+
|
Q221
|
AIV
|
H1N2
|
-
|
-
|
X1330
|
AIV
|
H3N2
|
-
|
-
|
P174
|
AIV
|
H4N2
|
-
|
-
|
A1267
|
AIV
|
H6N2
|
-
|
-
|
H7N3
|
AIV
|
H7N3
|
-
|
-
|
1605
|
AIV
|
H7N9
|
-
|
-
|
X169
|
AIV
|
H9N2
|
-
|
-
|
H9
|
AIV
|
H9N2
|
-
|
-
|
T55
|
AIV
|
H10N2
|
-
|
-
|
S82
|
AIV
|
H11N2
|
-
|
-
|
ND
|
NDV
|
/
|
-
|
-
|
JS1816
|
NDV
|
/
|
-
|
-
|
M41
|
IBV
|
/
|
-
|
-
|
H52
|
IBV
|
/
|
-
|
-
|
Detection and evaluation of clinical samples by H5 RT-RAA
We collected about 420 oropharyngeal and cloacal swab samples from 29 sites in live poultry markets in six provinces, which were immediately placed in 1 mL antibiotic-containing PBS as described above and then stored at -80°C until total nucleic acids were extracted with the above viral RNA extraction kit [21]. We evaluated the performance of the RT-RAA assay in 420 avian clinical samples and compared it with the published RT-qPCR method for H5 subtype AIV. The primers and probe of the RT-qPCR assay are listed in Table 1 [13]. The RT-qPCR assay was conducted with the One Step PrimeScript RT-PCR Kit (Takara).
Statistical analysis
To determine the RT-RAA detection limit, a probit analysis was performed at a 95% confidence interval (CI), and the k and p values of RT-qPCR and RT-RAA were calculated [14]. In addition, we calculated the sensitivity and specificity of RT-qPCR and RT-RAA for detection in clinical samples of poultry. All statistical analyses were performed in SPSS 21.0 (IBM).