Abstract

Nucleoprotein (N) gene from rabies virus (RABV) is a useful sequence target for variant studies. Several specific RABV variants have been characterized in different mammalian hosts such as skunk, dog, and bats by using anti-nucleocapsid monoclonal antibodies (MAbs) via indirect fluorescent antibody (IFA) test, a technique not available in many laboratories in Mexico. In the present study, a total of 158 sequences of N gene from RABV were used to design eight pairs of primers (four external and four internal primers), for typing four different RABV variants (dog, skunk, vampire bat, and nonhematophagous bat) which are most common in Mexico. The results indicate that the primer and the typing variant from the brain samples, submitted to nested and/or real-time PCR, are in agreement in all four singleplex reactions, and the designed primer pairs are an alternative for use in specific variant RABV typing.

1. Introduction

Despite the significant progress for prevention of the rabies disease and its control in the developing countries, this disease still causes over 60 thousand human deaths every year. Rabies disease is caused by infection with viruses of the family Rhabdoviridae, genus Lyssavirus [1]. Until now, fourteen species of Lyssavirus have been described in the world. Actually, the rabies virus (RABV) is the only one present in the American continent [26].

Although all mammals are susceptible to lyssaviruses, bats and carnivores are the major Lyssavirus reservoirs. In the Americas, distinct RABV variants are associated with different animals, such as foxes, coyotes, raccoons, skunks, and multiple species of nonhematophagous (frugivorous, insectivorous) and hematophagous bats [712]. In Mexico, we have been faced with less than ideal surveillance in animal populations. The reduced resources available are prioritized for diseases with overwhelming human morbidity and mortality. Accurate diagnosis and determination of RABV variants are paramount components of surveillance system and frequently are important from the perspective of veterinary and public health, when the source of exposure needs to be determined and relevant control strategies need to be implemented [13].

The direct fluorescent antibody test (FAT) is the “gold standard” for rabies diagnosis [14]; the modern conjugates used in FAT are able to detect antigens of all lyssaviruses described to date [15, 16]. Virus variants associated with certain host species can be distinguished by application of anti-nucleocapsid monoclonal antibodies (MAbs) via indirect IFA. The MAbs are still commonly used in Latin American countries, particularly in the laboratories lacking established molecular techniques [17]; these have been applied to Mexican rabies virus samples and provide data regarding the most likely reservoir species involved in rabies transmission and dissemination. Even though there has been a decrease in dog rabies, as a result of massive dog vaccination in Mexico, there is a high risk of an increase of human rabies cases transmitted from wild reservoirs as well as the simultaneous presence of more than one reservoir and more than one virus variant [1820].

The reactivity of certain viral isolates does not match the reactivity patterns in some cases [8, 1719]. Other molecular assays like the restriction analysis of RT-PCR amplified fragments of RABV genes were suggested for the differentiation of two major RABV variants but suffered from low specificity. Amplification and sequencing of viral genes followed by their phylogenetic analysis have provided more robust characterization. However, this approach requires expensive equipment and experienced laboratory staff, and it takes a relatively long processing time (typically at least 10–12 hours) [2023]. So, in this study, we designed eight pairs of primers of the RABV associated variant, which were used in a nested endpoint RT-PCR (four external and four internal primers) for the real-time RT-PCR assay, in order to detect and type the major RABV variants present in Mexico.

2. Materials and Methods

2.1. Primer Design

Primer design was based on the alignment constructed with ClustalW using complete RABV N gene sequences available in GeneBank associated variant; these were designed in consensus region. Two pairs of the primers (external and internal) were designed for each of the RABV variants associated, and the maximum average entropy (Hx) and the maximum entropy of each position were calculated using Bio Edit v7.2.5.

Two external primers and two internal primers were designed for dogs variant; 36 N gene sequences were obtained from different Mexican states; for the vampire bats variant, 18 N gene sequences were considered from Mexican states; for the nonhematophagous bat variant, the primer design comprised 50 N gene complete sequences from hosts Eptesicus, Myotis, and Nycticeius genera, distributed close to Mexico [21]; these genera are distributed from North America to Central America and have high diversity; in the case of skunks variant, 4 Mexican RABV sequences were considered; 34 RABV sequences were from USA and 13 CASK RABV sequences were from USA related to Mexican skunk rabies virus; previous studies consider two variants circulating in Mexico, MEXSK-2 and MEXSK-1 [22], located in South Baja California (SBC skunk) and Central Mexico; these are closely related and circulate predominantly in spotted skunks [23] (Table 1).

2.2. Samples

Twenty-three brain samples collected in Mexico were used as follows: nine brain samples tested negative by FAT and fourteen tested positive by FAT and typed by MAbs: these RABV isolated consisted of six samples of dog brain, one sample of skunk, two samples of cow, two samples of vampire bat, and three samples of nonhematophagous bat (Table 2).

2.3. Nucleic Acid Extraction

The brain tissues (approximately 3 mm3) were homogenized in 200 μL of lysis/binding buffer using MagNA Lyser Green Beads (Roche, Germany) and MagNA Lyser (Roche Applied Science, Germany). Total RNA was extracted from the homogenates using MagNA Pure LC Total Nucleic Acid Isolation kit (Roche, Germany) and MagNA Pure LC 2.0 (Roche Applied Science, Germany) following the manufacturer’s instructions. Total RNA was eluted in 200 μL buffer elution and quantified into Nanodrop (Invitrogen); RNA concentration was calculated considering  ng/μL, and all samples were adjusted at 20 ng/L of final concentration with elution buffer.

2.4. Nested RT-PCR Amplification and Sequence Determination

The reverse transcription reaction was performed using four singleplex reactions with external primers and SuperScript®III Platinum® One-Step qRT-PCR kit (Invitrogen) in 50 μL of reaction mixture containing 25 μL of 2x reaction mix, 1 μL of forward sense primer (10 μM), 1 μL of reverse sense primer (10 μM), 1 μL SuperScript III RT/Platinum TaqMix, 17 μL of DEPC-treated water, and 5 μL RNA extracted (20 ng/μL). Amplification was performed in C1000 Thermal Cycler (Bio-Rad, USA) using the following program: one cycle of RT at 50°C for 30 min, followed by denaturation at 92°C for 3 min, 35 cycles with denaturation at 92°C for 30 s, annealing at primer-specific temperature (Table 3) for 30 s, and elongation at 72°C for 1 min, with the final extension at 72°C for 4 min.

For nested PCR, 1 μL of the primary amplification products was added to a new singleplex PCR reaction using internal primers and Taq DNA polymerase kit; in a 50 μL total volume, add 5 μL 10x PCR buffer, 1 μL of 1x dNTP mix (200 μM of each dNTP), 1 μL internal forward primer, 1 μL internal reverse primer, 0.25 μL Taq DNA polymerase (1.25 units/reaction), 1 μL of primary amplification, and 40.75 μL RNase-free water. The thermal program consisted of a first cycle of 2 min at 94°C, followed by 35 repetitive cycles of denaturation of 1 min at 93°C, 1 min of annealing at the primer-specific temperature (Table 3), 1 min of elongation at 72°C, and the final elongation at 72°C for 4 min. The four singleplex RT-PCR and four nested PCR products were analyzed in 1-2% agarose gel. Bands of the expected size were excised, purified, and cloned in TOPO-TA vector (Invitrogen, Carlsbad, USA). The resulting plasmids were purified from E. coli colonies using Pure Link Quick Plasmid Miniprep kit (Invitrogen), sequenced with the universal M13 primers (Macrogen, Korea), and analyzed with MEGA 6.06 [24].

2.5. Real-Time RT-PCR with SYBR Green

The one-step real-time PCR was performed using internal primers and LCFastStart RNA Master SYBR Green I kit (Roche, Germany) in 20 μL of total volume, four singleplex reactions including 100 ng total of total RNA and 0.01 μM of each internal pair of primers for RABV associated variant (Table 3). Amplification was performed in LightCycler 2.0 (Roche, Germany) using the following program: one cycle of RT at 55°C for 30 min, followed by denaturation at 95°C for 30 s, 40 cycles with denaturation at 95°C for 10 s, annealing at 60°C for 15 s, and elongation at 72°C for 25 s. The measurement of the fluorescent signal was carried out during the extension phase at 530 nm. By the end of the amplification test, an analysis of the dissociation curves from the product was made to ensure the absence of hairpin and dimer formation. Hybridization temperatures and primer concentrations were optimized for each reaction based on the preliminary standardization experiments.

3. Results

The set of primers for specific RABV variants was designed aligning the sequence of N gene region. Four regions highly conserved were selected for two external primers and two internal primers designed, with more than 90% of conservation, for each variant and high variability between variants. Two mismatches were permitted for the primers design, and less was possible for the internal primers located at the primers beginning or end.

As the maximum entropy values increased, the number of identified conserved regions, their length, the coverage of conserved regions, and the average length of single conserved regions also increased. Two external primers and two internal primers were designed for dog-associated variant; the alignment presented a high conservation level (Figure 1). The maximum average entropy (Hx) was 0.04 and the maximum entropy of each position was 0.97. In the case of the set of primers for skunk-associated variant, the alignment presented high conservation level (Figure 2). The maximum average entropy (Hx) was 0.19 and the maximum entropy of each position was 0.99. The alignment of the set of primers for vampire bat-associated variant showed a high conservation level (Figure 3). The maximum average entropy (Hx) was 0.04 and the maximum entropy of each position was 0.98. Finally, the alignment of the set of primers for nonhematophagous bat-associated variant presented high conservation level (Figure 4). The maximum average entropy (Hx) was 0.07 and the maximum entropy of each position was 0.97 on average; this means that, at the same position of every base, a few sequences of alignment of the associated variant differed from the others and thus were considered conservative; with respect to the maximum average entropy, the variant associated with more differences was the nonhematophagous bat-associated variant (Figure 4).

The sequences and locations of the two pairs of variant-specific primers are listed in Table 3. Even when the melting temperature is similar between them, the sequence is dependent on the specific host variant.

All brain samples from Mexican host mammals were diagnosed as negative or positive in FAT, with a corresponding signal in the nested RT-PCR assay. A positive control for each variant was performed using a positive example previously MAbs tested. In the first step, the external amplification produced a single band of 608–1187 bp, while the second amplification of the primary PCR products with the internal primer showed products of 200–400 bp. To complement the nested information, one-step RT-PCR as well as the second nested RT-PCR was performed with external primers and the same samples (Figure 5).

The optimal annealing temperature for external RT-PCR was in the range 48–55°C and was 56°C in the case of nested RT-PCR. Optimal concentration of Mg2+ was in the order of 2.5–3 mM for both RT-PCR and nested RT-PCR reaction mixtures.

3.1. RT-PCR SYBR Green

For real-time RT-PCR assay, we used the internal primers (Table 3). This technique had an optimal annealing temperature of 60°C from four pairs of primers, and the dissociation temperature curves were as follows: 85.50°C for skunk specific variant, 80.19°C for dog-associated variant, 83.96°C for bat-associated variant, and 85.23°C for vampire bat specific variant (Figure 6).

3.2. Comparison of Diagnostic Methods

A total of twenty-three samples were assessed as follows: nine negative-control samples performed by nested or real-time RT-PCR assays showed no positive detection with the internal and external primers; previously, fourteen positive-RABV variants samples were tested by FAT; eleven of them were categorized with monoclonal antibodies resulting in six positive to variant 1 (V1), three to V8, one to V5, two atypical variants, and two undetected. The RABV specific variant characterizations to dog, vampire bat, nonhematophagous bat, and skunk were determined by real-time RT-PCR using external primers. However, to prevent cross-reaction and to increase sensitivity in the nested PCR, internal primers were used to confirm thus the abovementioned variants. In addition, real-time PCR detection using internal primers confirmed the reservoir variant with dissociation temperatures of 60°C. The amplified fragments were sequenced with a subsequent analysis by BLAST; this analysis confirmed the reservoir for nested PCR and real-time RT-PCR (Table 4).

3.3. Sensitivity of nRT-PCR and SYBR Green

Twenty-three brain samples were analyzed as follows: nine negative-control samples and fourteen positive samples were confirmed by nucleotide sequencing. Regarding the results of the nested PCR and real-time PCR assays of the brain samples, they showed 100% sensitivity (100% CI: 76.84% to 100.00%) and 100% specificity (100% CI: 66.37% to 100%).

4. Discussion

In some studies of antigenic characterization of rabies virus, a panel of eight anti-N protein monoclonal antibodies (MAbs) has been used, which can differentiate between eleven distinct variants harbored by a variety of terrestrial and chiropteran hosts [25, 26]. Application of this panel to rabies virus collections from many Latin American countries has identified two major variants, associated with dog and vampire bat (Desmodus rotundus), as well as other variants associated with several insectivorous bats, including the free-tailed bat (Tadarida brasiliensis) and the hoary bat (Lasiurus cinereus) [27].

Real-time RT-PCR techniques have been used for diagnosis and genotyping of all the Lyssavirus genus including the RABV [28]. The nested RT-PCR assay, which requires both multiple transfers of material and substantial time, is sufficient to detect virus from each virus-positive brain sample [29] and therefore still offers a useful tool for variants rabies diagnosis where conventional PCR technology exists.

The real-time RT-PCR detection with SYBR Green, where the specificity is being given by the primers, is an easy-to-use assay to detect infected brain material in a single tube test and, consequently, is an attractive option for laboratory use as a screening surveillance tool. In the present study, these latest technologies for typing the RABV variants depending on the host (vampire bat, skunk, dog, and bat) were used.

The real-time RT-PCR detection with SYBR Green, whose specificity is given by the primers, is an easy-to-use assay to detect infected brain material in a single tube test and, consequently, is an attractive option for laboratory use as a screening surveillance tool. In the present study, these latest technologies for typing the RABV associated variants on the host (vampire bat, skunk, dog, and bat) were used.

In the design with highly specific primers from the conserved region from the nucleoprotein of RABV, the maximum average entropy (Hx) was in the order of 0.03–0.19 and the maximum entropy of each position was 0.97–0.99. In addition, the positions of different primers in N gene sequence are close but different for variant host, increasing the specificity.

The current gold standard test has been and is the fluorescent antibody test (FAT), which uses a conjugated monoclonal antibody against the RABV nucleoprotein. Although it is cheaper, some laboratories have no access to MAbs but have PCR and/or real-time technology.

In the characterization of the antigenic variants (AgV) with MAbs in the dog samples, the dog variant-specific primers identified the dog variant (V1). This result matched both the nRT-PCR and SYBR Green primers at 100%. Similarly, the skunk samples matched the same percentage with skunk variant-specific primers. Furthermore, in the bovine samples where the MAbs detection identified the skunk variant (V8), the determined host by nRT-PCR and SYBR Green was diagnosed as positive with the vampire primers, with this last result confirmed by sequences and AC Numbers JQ037818 to JQ037831 (Table 1). The real-time RT-PCR result coincides with some other studies where rabies transmission from vampire bats to bovines has been described.

The MAbs detection in nonhematophagous bat was V5 bat, a result which coincides with both nRT-PCR and SYBR Green with the bat primer. The vampire bat 110 sample was determined as atypical and the vampire bat 3919 was not determined with the MABs; however, both samples were diagnosed as positive with the bat primers for nRT-PCR and SYBR Green.

In some cases, the classification of certain rabies virus isolates by monoclonal panel can obtain nontypical reactivity patterns and is not assigned to any known variant, as found in certain Argentinian rabies viruses [12]. The application of molecular genetic techniques for characterization of viral collections can assist in resolving such typing difficulties.

In the hematophagous bat samples determined as atypical and the one not determined with the MAbs, it was concluded that the host was a vampire bat by nRT-PCR and SYBR Green detection. This may have occurred due to the high sensitivity of the RT-PCR molecular technique, as it has been shown in studies where positive results in brains analysis were demonstrated by nRT-PCR and negative results by FAT [30, 31]. In all results, the host was confirmed by the amplicon sequencing. The access numbers are shown in Table 4.

According to the RABV variant detection, the external primers and internal primers detect a specific variant and do not present cross-reaction between them, and the final result is given for the internal primer reaction in nested and/or RT-PCR real time, as they were obtained in different samples (Table 5).

In addition, this study showed 100% sensitivity and 100% specificity assessed by nRT-PCR and real-time RT-PCR with SYBR Green. These findings are an early estimate by what is required of a greater number of related studies, increasing the number of samples to obtain better sensitivity and specificity evaluation. However, this assay could be useful, for institutions without access to MAbs and those that have PCR and/or real-time technology as an alternative.

The relevance of the present study falls in the rabies virus typing from original host-brains samples and the association with the variant-specific host performed by nested endpoint PCR or real-time RT-PCR assays. Previous studies report the detection in decomposed brains from dogs and humans samples [32]; in humans exhumed between 8 and 30 days after burial [27]; in wolves by nested RT-PCR [33]; in mice previously infected by heminested RT-PCR [29]; and in bats and herbivores by RT-PCR.

The sequence obtained for this study, a splitting between the urban rabies (dog) and the sylvan rabies (bat, vampire bat, and skunk), was shown in Tables 4 and 5; the results were according to the primers designed associated variant for the dog (urban rabies) and bat, skunk, and vampire bat (sylvan rabies).

5. Conclusion

This study describes the development of an alternative tool for RABV typifying in real-time RT-PCR and/or nested RT-PCR, considering dog, skunk, vampire bat, and nonhematophagous bat specific variants.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The authors would like to acknowledge Lizdah Ivette García Rodríguez and José A. Valdes-Zúñiga, Unidad de Enseñanza, Investigación y Calidad, for supporting the administrative project permission and formats, IPN: COFAA SIP.