This protocol presents a qRT-PCR-based approach for determining the rabies virus nucleoprotein (N) gene copy number within various bovine brain anatomical structures using in vitro transcription.
Bovine paralytic rabies (BPR) is a form of viral encephalitis that is of substantial economic importance throughout Latin America, where it poses a major zoonotic risk. Here, our objective was to utilize a laboratory protocol to determine the relative copy number of the rabies virus (RABV) genome in different bovine brain anatomical structures using quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). qRT-PCR quantifies the specific number of gene copies present in a sample based on fluorescence emitted after amplification that is directly proportional to the amount of target nucleic acid present in the sample. This method is advantageous owing to its short duration, reduced risk of contamination, and potential to detect viral nucleic acids in different samples more easily compared to other techniques. The brains of six rabid animals were divided into six anatomical structures, namely the Ammon’s horn, cerebellum, cortex, medulla, pons, and thalamus. All brains were identified as positive for RABV antigens based on a direct immunofluorescence test. The same anatomical structures from the brains of four RABV-negative bovines were also assessed. RNA was extracted from each structure and used for qRT-PCR. An assay was performed to determine the copy numbers of RABV genes using an in vitro transcribed nucleoprotein gene. The standard curve used to quantify viral RNA exhibited an efficiency of 100% and linearity of 0.99. Analysis revealed that the cortex, medulla, and thalamus were the ideal CNS portions for use in RABV detection, based on the observation that these structures possessed the highest levels of RABV. The test specificity was 100%. All samples were positive, no false positives were detected. This method can be used to detect RABV in samples that contain low levels of RABV during diagnosis of BPR.
Rabies can be confirmed ante-mortem and post-mortem by various techniques that enable the detection of viral nucleic acids in the brain, skin, urine, or saliva1. Detection of the rabies virus nucleoprotein (N) gene is primarily used for rabies diagnosis by molecular tests. This gene is also used for viral genotyping. Rabies can be diagnosed in animals using any portion of the affected brain; however, to exclude the possibility of rabies, tissue from at least two regions in the brain must be tested2. Several diagnostic methods exist for rabies detection in animals; however, the direct immunofluorescence test remains the standard reference technique3. Other tests include biological tests that incorporate mouse inoculation, infection in tissue culture, and polymerase chain reaction (PCR)4. All these techniques are recommended by the World Health Organization (WHO) and World Organization of Animal Health (OIE) for the diagnosis of rabies in humans and animals, respectively5.
Nucleic acid detection and amplification techniques have revolutionized the diagnosis of rabies in recent years6 and these techniques play an important role in the ante mortem diagnosis of human rabies. Several PCR-based tests have been evaluated to complement conventional tests for ante-mortem and post-mortem rabies diagnoses7,8,9,10. Most assays target rabies viral nucleoprotein gene for amplification which is the most highly conserved region in the viral genome1,11. In the last 20 years, various molecular assays have been developed to diagnose RABV, and some of these assays have been used for virus characterization. Most trials have aimed to detect conserved genes within the viral genome, most commonly by using conventional or quantitative real-time polymerase chain reaction (qRT-PCR) assays12,13.
PCR is a highly sensitive diagnostic technique that can detect the genome of a given pathogen within tissues, even when these tissues are decomposed. Using PCR-based approaches, minimal quantities of an infectious agent can be detected in a clinical sample through the selective and repetitive amplification of a DNA nucleotide sequence14. qRT-PCR that incorporates fluorescent probes (e.g., TaqMan) or DNA binding dyes (e.g., SYBR Green) has been used in trials to diagnose RABV both ante– and post– mortem with high sensitivity; however, such an approach requires specialized equipment. To overcome this limitation, reverse transcription loop-mediated isothermal amplification (RT-LAMP) has been suggested, based on its low cost, simplicity, and desirable characteristics for the detection of RABV. This assay is particularly important as it can be used in developing countries15.
qRT-PCR is based on the detection and quantification of a molecule, where fluorescent signal increases are in direct proportion to the amount of PCR product in a single reaction. As the number of copies of the nucleic acid target increases, so does the fluorescence. Non-specific intercalating dyes such as the SYBR Green DNA-binding dye or sequence-specific oligonucleotide probes carrying a fluorophore and a quencher are commonly used to provide the fluorescent readout in qRT-PCR. This assay offers advantages over conventional RT-PCR that include a shorter test time (2–4 h), reduced risk of contamination due to the closed tube system (lack of post-PCR manipulation of amplified products), and the ability to detect different targets simultaneously16. qRT-PCR can be used to diagnose rabies ante-mortem from saliva and other samples. This assay can also be used as a universal real-time test for the detection of different Lyssavirus species or lineages of RABV15. In this combo RT-PCR approach, two reactions are used. The first detects different genetic linages of RABV, and the second detects the Lyssavirus species. Both steps involve qRT-PCR assays, where the first uses hybridization probes and the second uses dyes15. Owing to the large number of tests that have demonstrated successful molecular detection of RABV using this technique, the current OIE Terrestrial Manual (2018) recommends the use of PCR for the molecular detection of RABV17.
Mexico is a country with considerable livestock potential. The states with the highest livestock production contain both humid tropical regions and dry regions that are at risk for rabies outbreaks due to the presence of the vampire bat Desmodus rotundus, the main transmitter of rabies. Therefore, it is essential to develop more tools for the prevention and control of bovine paralytic rabies (BPR) in México. Based on this, the aim of this study was to use quantitative qRT-PCR to determine the number of viral particles in different anatomical structures of bovine brains following death due to rabies infection.
Six brains obtained from RABV-positive bovines were donated by an external laboratory for the use in the development of qRT-PCR protocol described below. The bovine brain structure samples were cold-chain transported to the INIFAP CENID-MA laboratory BSL-2 facility and stored at -80 °C until use. Brains were obtained from animals from the states of Campeche, Yucatán, and Querétaro. Prior to the receipt, various structures were dissected from the brains. These structures included the Ammon’s horn, cerebellum, cortex, medulla, pons, and thalamus18,19. Genetic material was extracted as described below. RABV diagnosis was confirmed using direct fluorescent antibodies (DFA)20. As a positive control, mouse brains that were inoculated with RABV21 were used. Additionally, four RABV-negative cattle brains (as determined by DFA) were used as negative controls.
This study was approved by and conducted in strict accordance with the recommendations for the use of animals provided by the Institutional Animal Care and Use Committee (IACUC) of the Centro Nacional de Investigación Disciplinaria en Microbiología Animal (CENID-MA).
1. Samples
2. Direct fluorescent antibodies (DFAs) to confirm rabies
NOTE: Detection of the rabies antigen by DFA is a qualitative method to determine the presence of the rabies nucleoprotein using fluorescein-labeled antibodies. The test was performed using the protocol provided by Dean et al. (1966)20, as described below.
3. Generating positive control for RT-PCR
4. RNA extraction
NOTE: Total RNA was extracted directly from bovine brains using an organic extraction method according to following protocol.
5. cDNA synthesis and PCR
6. In vitro transcription
NOTE: In vitro transcription generates mRNA of a target gene. Use a primer pair that amplifies the complete RABV N gene and one that is used as a positive control for the qRT-PCR assay. These primers are designed to amplify the complete RABV N gene and to add a promoter that recognizes the T7 polymerase (TAATACGACTCACTATAG).
7. Real-time reverse transcription polymerase chain reaction (qRT-PCR)
8. Calibration curve or standard curve
9. qRT-PCR of biological samples
DFA results showed 100%, 100%, 83.3%, 66%, and 50% positivity for RABV in the cortex, thalamus, medulla, pons, and horn, respectively. These results confirmed the previous results, and at least three of the structures dissected from each brain were positive for RABV. A representative positive DAF staining is shown in Figure 1.
Figure 2 shows the amplification of a fragment of the RABV N gene (step 5.5) with the primers first ones reported by Loza-Rubio et al.11. This showed that the material obtained to be used in amplification was of good quality.
The size of the PCR-amplified fragment is shown in Figure 3. The equation of the standard curve shows the relationship between the amount of RABV genetic content in a sample and the size of the PCR-amplified fragment. The amount of RABV genetic content (in ng) was deduced by interpolation of the Ct values in the calibration curve using the equation presented in Figure 4. Using this equation, the detection limit of 1 x 105 μg/μL of viral RNA was found to correspond to 36 copies of the RABV genome. These results demonstrate that the assay is sensitive and can be used to detect the virus in samples that contain a low number of copies RABV N gene. The efficiency of the assay was calculated using the slope of the standard curve, which was -3.28. The samples used for qRT-PCR showed amplification of the RABV N gene.
To determine the number of viral copy present, the Ct values for each sample were interpolated using the calibration curve equation. The result obtained (in nanograms) was substituted into the formula described below from the following ref.24:
Number of copies RABV N gene = (ng sample * 6.022 x 1023) / (104 bp (length of PCR product) * 1 x 109 * 650).
Comparatively, the results of qRT-PCR assays were consistent with those obtained using DFA. According to the results obtained in the qRT-PCR test in cases C and D (Table 1), the thalamus is the structure that possesses the highest number of copies of the RABV N gene. This suggests that early infection of hypothalamic and thalamic neurons is important in the development of the disease, given that these neurons control the vegetative functions of an animal. The copy numbers of the RABV N gene that were detected in each structure of each sample are shown in Table 1. The sensitivity and specificity of the qRT-PCR test were both 100%, as all samples were positive. This result is consistent with the initial diagnoses of the samples.
Figure 1: Bovine brain tissues showing positive staining according to the direct immunofluorescence test that detects rabies virus protein N in infected tissues. An apple green coloration is observed upon staining, where the antibody binds to the rabies antigen. Scale bar: 50 µm Please click here to view a larger version of this figure.
Figure 2: 1.5% agarose gel showing the amplification products. Amplification of a 761 bp fragment of the RABV N gene to demonstrate the presence and quality of the cDNA to be used in in vitro transcription. Lane 1 contains a molecular weight marker, line 2: amplification of positive control; lane 3: negative control (non-infected sample); line 4: amplification of cDNA used for in vitro transcription. Please click here to view a larger version of this figure.
Figure 3: 1.5% agarose gel showing the amplification products. Amplification of the complete N gene is shown in lane 2. Lane 1 contains a molecular weight marker, and lane 3 contains the negative amplification control. Please click here to view a larger version of this figure.
Figure 4: Standard curve of the qRT-PCR using in vitro transcribed N gene mRNA to quantify the number of copies of the RABV N gene. Please click here to view a larger version of this figure.
ID brain | Structure | DAF | Number copies (x109) |
A | Ammon’s horn | + | 0.261 |
Cerebellum | + | 0.0663 | |
Cortex | + | 0.02 | |
Medulla | + | 0.146 | |
Pons | + | 30.7 | |
Thalamus | + | 108 | |
B | Ammon’s horn | – | 0.0251 |
Cerebellum | + | 4.64 | |
Cortex | + | 12.4 | |
Medulla | + | 0.175 | |
Pons | ‒ | 0.0721 | |
Thalamus | + | 121 | |
C | Ammon’s horn | ‒ | 0.168 |
Cerebellum | + | 0.0239 | |
Cortex | + | 21.5 | |
Medulla | + | 17.2 | |
Pons | + | 1.32 | |
Thalamus | + | 102 | |
D | Ammon’s horn | + | 11.8 |
Cerebellum | + | 0.0239 | |
Cortex | + | 8.72 | |
Medulla | + | 1.25 | |
Pons | ‒ | 0.0165 | |
Thalamus | + | 33.4 | |
E | Ammon’s horn | + | 4.15 |
Cerebellum | + | 6.78 | |
Cortex | + | 1.53 | |
Medulla | + | 1.41 | |
Pons | + | 0.025 | |
Thalamus | + | 95 | |
F | Ammon’s horn | + | 0.0221 |
Cerebellum | ‒ | 0.00023 | |
Cortex | + | 4.95 | |
Medulla | ‒ | 0.000556 | |
Pons | + | 1.02 | |
Thalamus | + | 89 |
Table 1: Determination of the number of copies of the RABV N gene within each brain structure. Results were obtained from six anatomical structures of RABV-positive bovine brains diagnosed using DFA. The words highlighted in bold indicate the structures with the most copies. The results were expressed as number of copies of the N gene per milligram of tissue.
Previous studies have shown that DFA can only detect RABV within seven days of the sample being stored at room temperature (21 °C)14. In contrast, this work demonstrated that the sensitivity of RT-PCR begins to decrease after the samples have been exposed to room temperature for 12 days. Therefore, the RABV genome can be detected by qRT-PCR in samples exposed to room temperature for up to 23 days. This demonstrates that the sensitivity of qRT-PCR is relatively higher for more highly decomposed samples. However, it remains necessary to develop simpler methods that do not compromise specificity25. The present research demonstrated a molecular assay that possesses higher sensitivity than DFA, based on the observation that it was able to detect the viral genome regardless of the specific brain structure.
The advantages of the RT-PCR technique for the use in detecting infectious agents have been highlighted in several scientific studies. However, this technique can exhibit some disadvantages, such as cost of execution, as it requires specialized equipment. Trained personnel are required to handle molecular biology assays. Additionally, it should be mentioned that as a highly sensitive technique, some of the steps are critical to obtain the best results. One of these is the proper handling of samples for RNA extraction. This step is crucial, as RNA is easily degraded, and the results could, therefore, be affected. Consider maintaining the sample in cold conditions (approximately 4 °C) during the process. In addition, consider that several rounds of PCR application are carried out as mentioned in step 6.1 to generate the transcript in vitro. This is because a large amount of genetic material is necessary for the reaction to yield substantial (more than milligram) quantities of DNA. Another crucial aspect of this assay is the requirement for purification of the genetic material in the columns, as DNA may also be lost during this step.
qRT-PCR is known to be as specific as the DFA test, which is the standard test approved by the WHO and the OIE. However, molecular assays employing RT-qPCR have been shown to have a higher sensitivity and to thus allow for analysis of samples with a higher degree of decomposition to facilitate the detection of a relatively small number of copies of the viral genome. It is also much faster, reduces the risk of contamination, and can be used ante-mortem26,27.
Bingham and van Der Merwe (2002)28 conducted a study using DFA to determine which brain structures possessed higher concentrations of the RABV N gene. A total of 252 brain structures from several species were evaluated. The thalamus, pons, and medulla were the most relevant brain structures, as they were positive in all evaluated brain samples. These results agree with the results presented here, which were consistent with DFA results in the following structures: cortex and thalamus, 100%; medulla, 83.3%; pons, 66%; and horn, 50%. In contrast, previous studies reported no false negatives. In this work, some negative results were detected in complete brain structures that had previously been diagnosed as positive. This may be because the portion of the sample used for the test did not contain viral particles identified by the antibody or the more sensitive molecular assay. Another possible explanation is that the samples were not properly transported or stored prior to reception at the laboratory.
The qRT-PCR assay conducted in this study could detect as few as 36.3 copies of the RABV N gene per milligram of tissue. The most suitable brain structures for qRT-PCR included the cortex and thalamus. This is based on the observations that higher concentrations of the RABV N gene were detected in these structures and that they were also found to be positive using the RABV reference test. The test was able to detect very low concentrations (0.00023 x 109 copies) of the virus in various samples. In addition to quantifying the viral particles in the sample, the test appears to provide a good alternative diagnostic and research tool for the detection of RABV.
With the results obtained using the rabies virus viral genome detection protocol presented here, it is anticipated that this technique can be applied for the molecular diagnosis of the virus in different types of samples, as it is capable of detecting minimal amounts of the viral genome. Its biggest advantage over other methods such as DFA is that it is more sensitive and can be used ante-mortem, as has been suggested by other authors.
The authors have nothing to disclose.
This work was supported by the National Institute of Agricultural, Forestry, and Livestock Research (INIFAP). We thank Jerzayn Fraustro Esquivel for his collaboration in the development of the video associated with this document.
Chloroform | SIGMA | C7559 | Facilitates recovery of the aqueous phase of PCRs which have been overlaid with mineral oil. |
DNA Clean & Concentrator-500 | Zimo | D4031 | The DNA Clean & Concentrator-500 (DCC-500) is designed for the rapid, large format purification and concentration of up to 500 µg of high quality DNA from samples including large-scale restriction endonuclease digestions and impure DNA preparations. |
Ethanol | Amresco | 193-500 | Purification of nucleic acids |
FastStart High Fidelity PCR System kit, dNTPack | Roche | 3553400001 | High fidelity enzyme for the amplification of PCR products avoiding random base changes |
GelDoc XR | BioRad | XR+ | Analyzes larger protein and DNA gels |
GelRed | Biotium | 41003 | A new generation of nucleic acid gel stains, they possess novel chemical features designed to minimize the chance for the dyes to interact with nucleic acids in living cells. |
iCycler Thermal Cycler Gradient | BioRad | 582BR | Temperature can be monitored and controlled by instrument algorithm, in-sample probe, or sample block modes |
iTaq Universal Probes One-Step Kit | BioRad | 1725141 | Reaction mixture to carry out PCR reactions in real time using TaqMan type hybridization probes |
Isopropyl alcohol | Amresco | 0918-500 | Precipitation of nucleic acids |
QIAquick gel extraction kit | Qiagen | 28706 | QIAquick Kits contain a silica membrane assembly for binding of DNA in high-salt buffer and elution with low-salt buffer or water. The purification procedure removes primers, nucleotides, enzymes, mineral oil, salts, agarose, ethidium bromide, and other impurities from DNA samples |
M-MLV Reverse Transcriptase | Invitrogen | 28025-021 | Moloney Murine Leukemia Virus Reverse Transcriptase (M-MLV RT) uses singlestranded RNA or DNA in the presence of a primer to synthesize a complementary DNA strand. |
NanoDrop 2000 | Thermo-Scientific | ND2000 | Microvolume Spectrophotometer |
Oligo(dT)18 primer | Invitrogen | SO132 | The oligo (dT)18 primer is a synthetic single-stranded 18-mer oligonucleotide with 5'- and 3'-hydroxyl ends. |
RiboMAX Large Scale RNA Production Systems kit SP6 and T7 | Promega | P1300 | In vitro transcription reactions are used to synthesize microgram amounts of RNA probes from recombinant DNA templates. Most transcription reactions designed to generate RNA probes are optimized to maximize incorporation of radiolabeled ribonucleotides rather than to produce large amounts of RNA |
Taq DNA Polymerase | Invitrogen | #EP0402 | Taq DNA Polymerase is a highly thermostable DNA polymerase of the thermophilic bacterium Thermus aquaticus. The enzyme catalyzes 5’ and 3’ synthesis of DNA |
TRIzol reagent | Invitrogen | 15596026 | The TRIzol reagent is a complete, ready-to-use reagent, designed for the isolation of high quality total RNA or the simultaneous isolation of RNA, DNA and proteins from a variety of biological samples. |