One-year Surveillance of Human Enteric Viruses in Raw and Treated Wastewaters, Downstream River Waters, and Drinking Waters.

Human enteric viruses are a major cause of waterborne diseases, and can be transmitted by contaminated water of all kinds, including drinking and recreational water. The objectives of the present study were to assess the occurrence of enteric viruses (enterovirus, norovirus, adenovirus, hepatitis A and E virus) in raw and treated wastewaters, in rivers receiving wastewater discharges, and in drinking waters. Wastewater treatment plants' (WWTP) pathogen removal efficiencies by adenovirus quantitative real-time PCR and the presence of infectious enterovirus, by cell culture assays, in treated wastewaters and in surface waters were also evaluated. A total of 90 water samples were collected: raw and treated wastewaters (treated effluents and ultrafiltered water reused for industrial purposes), water from two rivers receiving treated discharges, and drinking water. Nested PCR assays were used for the identification of viral DNA/RNA, followed by direct amplicon sequencing. All raw sewage samples (21/21), 61.9 % of treated wastewater samples (13/21), and 25 % of ultrafiltered water samples (3/12) were contaminated with at least one viral family. Multiple virus families and genera were frequently detected. Mean positive PCRs per sample decreased significantly from raw to treated sewage and to ultrafiltered waters. Moreover, quantitative adenovirus data showed a reduction in excess of 99 % in viral genome copies following wastewater treatment. In surface waters, 78.6 % (22/28) of samples tested positive for one or more viruses by molecular methods, but enterovirus-specific infectivity assays did not reveal infectious particles in these samples. All drinking water samples tested negative for all viruses, demonstrating the effectiveness of treatment in removing viral pathogens from drinking water. Integrated strategies to manage water from all sources are crucial to ensure water quality.

Hepatitis A and E Viruses in Wastewaters, in River Waters, and in Bivalve Molluscs in Italy.

Several studies have reported the detection of hepatitis A (HAV) and E (HEV) virus in sewage waters, indicating a possibility of contamination of aquatic environments. The objective of the present study was to assess the occurrence of HAV and HEV in different water environments, following the route of contamination from raw sewage through treated effluent to the surface waters receiving wastewater discharges . Bivalve molluscan shellfish samples were also analyzed, as sentinel of marine pollution. Samples were tested by RT-PCR nested type in the VP1/2A junction for HAV, and in the ORF1 and ORF2 regions for HEV. Hepatitis A RNA was detected in 12 water samples: 7/21 (33.3%) raw sewage samples, 3/21 (14.3%) treated sewage samples, and 2/27 (7.4%) river water samples. Five sequences were classified as genotype IA, while the remaining 7 sequences belonged to genotype IB. In bivalves, HAV was detected in 13/56 samples (23.2%), 12 genotype IB and one genotype IA. Whether the presence of HAV in the matrices tested indicates the potential for waterborne and foodborne transmission is unknown, since infectivity of the virus was not demonstrated. HEV was detected in one raw sewage sample and in one river sample, both belonging to genotype 3. Sequences were similar to sequences detected previously in Italy in patients with autochthonous HEV (no travel history) and in animals (swine). To our knowledge, this is the first detection of HEV in river waters in Italy, suggesting that surface water can be a potential source for exposure .

Hepatitis E virus genotypes 1 and 3 in wastewater samples in Tunisia.

Hepatitis E represents an important public-health concern throughout the world. It is one of the leading causes of hepatitis in North Africa, Asia and the Middle East. In Tunisia, the true burden of HEV infection is still unknown. The objectives of the present study were to assess the occurrence of hepatitis E virus in Tunisia through the monitoring of urban sewage and to characterize the strains identified using molecular assays. A total of 150 sewage samples (raw and treated) were collected from three wastewater treatment plants located in the regions of Monastir and Mahdia and analyzed by nested RT-PCR using a qualitative assay targeting the methyltransferase gene in ORF1. Of these, only three samples (2 %) were found to be positive for HEV, one belonging to genotype 1 and two to genotype 3. The results of the present study indicate a low level of virus excretion among the Tunisian population. Both genotypes 1 and 3 are circulating in this country, however, possibly causing sporadic infections. The presence of the zoonotic genotype 3, known to be transmitted to humans mainly by swine and demonstrated in Tunisia for the first time in this work, raises the question of possible reservoir species, since pork products are not consumed in this country, pigs are not bred, and wild boar is not endemic. Further studies will be needed to gather information on the occurrence and diversity of HEV strains circulating among humans and animals in Tunisia, and on possible animal reservoirs.

Frequent and abundant Merkel cell polyomavirus detection in urban wastewaters in Italy.

Viruses strongly associated with human cancer have recently been detected in urban sewages and other water environments worldwide. The aim of the present study was to assess the presence of Merkel cell polyomavirus (MCPyV), a newly discovered, potentially oncogenic human virus, in urban sewage samples collected at wastewater treatment plants (WTPs) in Italy. A total of 131 raw sewage samples were collected from 21 WTPs in nine Italian regions and analyzed by both qualitative (PCR/nested) and quantitative (Real-Time qRT-PCR) methods. Of these, 66 samples (50.3 %) were positive for MCPyV by the qualitative assay. Quantitative data showed high viral loads in wastewaters (mean, 1.5E + 05 genome copies/liter). High concentrations of MCPyV were found in all WTPs under study, suggesting a wide circulation of the virus and thus the need for further studies to assess possible waterborne MCPyV transmission.

Qualitative and Quantitative Assessment of Hepatitis A Virus in Wastewaters in Tunisia.

Hepatitis A causes substantial morbidity in both industrialized and non-industrialized countries and represents an important health problem in several southern Mediterranean countries. The objectives of the study were as follows: (a) to assess the occurrence of hepatitis A virus (HAV) in Tunisia through the monitoring of urban wastewaters collected at wastewater treatment plants (WTPs); (b) to characterize environmental strains; and (c) to estimate the viral load in raw and treated sewages, in order to evaluate the potential impact on superficial waters receiving discharges. A total of 150 raw and treated wastewaters were collected from three WTPs and analyzed by both qualitative (RT-PCR/nested) and quantitative (qRT-PCR) methods. Of these, 100 (66%) were found to be positive for HAV by the qualitative assay: 68.3% in influents and 64.7% in effluents. The vast majority of HAV sequences belonged to sub-genotype IA, with 11 different strains detected found to be identical to clinical strains isolated from Tunisian patients with acute hepatitis. Five unique variants were also detected, not previously reported in clinical cases. Only two IB strains were found, confirming the rarity of this sub-genotype in this country. The results of the present study indicate a wide circulation of the pathogen in the population, most probably in the form of asymptomatic infections, a finding consistent with the classification of the country as having intermediate/high endemicity. Quantitative data showed high viral loads in influents (3.5E+05 genome copies/liter, mean value) as well as effluents (2.5E+05 genome copies/liter, mean value), suggesting that contaminated water could be a critical element in transmission.

GIV Noroviruses in Wastewaters and in Stool Specimens from Hospitalized Patients.

Noroviruses (NoVs) are important human pathogens associated with foodborne and waterborne gastroenteritis. These viruses are genetically highly heterogeneous, with more than forty genotypes within three genogroups (GI, GII, and GIV) identified in humans. However, the vast majority of human infections are associated with variants of a unique genotype, GII.4. Aside from these NoV strains of epidemiological relevance, NoV strains of genogroup GIV (Alphatron-like) are reported in a sporadic fashion and their overall prevalence in the community is unknown and this likely reflects the lack of specific diagnostic tools. We analyzed raw sewages collected from 32 wastewater treatment plants distributed throughout Italy (307 samples) and stool specimens collected from hospitalized patients with clinical signs of diarrhea of unknown etiology (285 samples). By using specific qualitative and quantitative RT-PCR assays, 21.8 % of the sewage samples and 3.2 % of the stool specimens tested positive for GIV NoVs. The number of genome copies in fecal samples ranged from 5.08 × 104 to 1.73× 106/g of feces. Sequence analysis showed limited genetic variability in human GIV viruses. The presence of GIV NoV both in sewage and in clinical samples confirms that not only GI and GII NoVs but also GIV strains are circulating in humans. Monitoring of GIV NoV is recommended in order to understand the dynamics of circulation in human populations, environmental contamination, and potential health risks.

The impact of temperature on the inactivation of enteric viruses in food and water: a review.

Temperature is considered as the major factor determining virus inactivation in the environment. Food industries, therefore, widely apply temperature as virus inactivating parameter. This review encompasses an overview of viral inactivation and virus genome degradation data from published literature as well as a statistical analysis and the development of empirical formulae to predict virus inactivation. A total of 658 data (time to obtain a first log(10) reduction) were collected from 76 published studies with 563 data on virus infectivity and 95 data on genome degradation. Linear model fitting was applied to analyse the effects of temperature, virus species, detection method (cell culture or molecular methods), matrix (simple or complex) and temperature category (<50 and ≥50°C). As expected, virus inactivation was found to be faster at temperatures ≥50°C than at temperatures <50°C, but there was also a significant temperature-matrix effect. Virus inactivation appeared to occur faster in complex than in simple matrices. In general, bacteriophages PRD1 and PhiX174 appeared to be highly persistent whatever the matrix or the temperature, which makes them useful indicators for virus inactivation studies. The virus genome was shown to be more resistant than infectious virus. Simple empirical formulas were developed that can be used to predict virus inactivation and genome degradation for untested temperatures, time points or even virus strains.

Molecular characterization of adenovirus from clinical samples through analysis of the hexon and fiber genes.

Human adenoviruses (HAdVs) are common pathogens associated with a variety of clinical manifestations. Although most infections are self-limiting, HAdVs can cause severe or lethal infections in immunocompromised as well as in healthy individuals. Several HAdVs have recently been characterized as emerging pathogens. In Italy, epidemiological, and especially molecular epidemiological, information on this pathogen is scarce. This study describes the characterization by cell culture, PCR and phylogenetic analysis of HAdV strains originating from a small collection of clinical samples gathered between 2008 and 2009. The distribution of different HAdV species was studied and the possible presence of newly emerging types was ascertained. A broad-range primer pair was used, targeting a portion of the hexon gene, in combination with species-specific primer pairs targeting a portion of the fiber gene. Human and animal reference AdV strains were included in the study. The broad-range assay identified all HAdV strains (study and reference samples), as well as three out of four animal AdV reference strains. Seven different types belonging to three HAdV species (B, C and F) were identified in the study samples. Species C was by far the most frequent. Two co-infections were detected, each with two serotypes within species C (types 1/2 and 2/6). The combined use of these two PCR assays--allowing not only the identification of known types but also, potentially, the discovery of newly emerging ones--can provide valuable epidemiological information on the spread of HAdVs.

Detection and molecular characterization of noroviruses from five sewage treatment plants in central Italy.

Noroviruses (NoVs) are the most frequent etiological agents of non-bacterial gastroenteritis. These viruses are transmitted through the fecal-oral route, leading to high viral loads in sewages. The objective of this paper was to study the environmental occurrence of the most prevalent NoV strains in different wastewater treatment plants. In addition, molecular characterization of the isolated strains was performed. Two different PCR-based methods were carried out and a novel strategy was used to verify the level of RT-PCR inhibition. From May to September 2007, a total of 97 inflow and outflow samples were collected from five wastewater treatment plants in central Italy. We detected NoV by nested PCR in 96.9% of influent samples: 89.1% contained both genogroups; 4.7% contained only GI and 3.1% only GII. In effluents, we detected NoV in 78.8% of samples: 30.3% contained both genogroups, and 48.5% contained only GI. The major genotypes detected by sequencing analyses were GI/2, GI/5, GII/b, GII/4 and GII/6. This work confirms the wide circulation of NoVs in Italy with a predominance of GI strains, and the widespread distribution of NoV variants in both raw and treated wastewater.

Detection of genogroup IV noroviruses in environmental and clinical samples and partial sequencing through rapid amplification of cDNA ends.

Noroviruses (NoVs) give rise to clinically relevant gastroenteritis in all age groups and are widely distributed in both clinical and environmental settings. NoVs are classified into five genogroups (GI to GV), of which GI, GII and GIV infect humans. While data on the epidemiology of human NoVs GI and GII have been steadily increasing, very little information has been published on the spread of GIV in either the health care system or the environment, resulting in a lack of information about its clinical significance and pathogenesis. In order to investigate the distribution of GIV strains in the environment, we analyzed sewage samples collected from five treatment plants, by using newly designed nested RT-PCR assays. A collection of clinical stool samples, originating from pediatric patients with symptoms of acute gastroenteritis, previously analyzed in our laboratory for the presence of NoV GI or GII, was also analyzed for the presence of GIV norovirus. Results of this work attest to the presence of GIV in both clinical and environmental contexts and underline the importance of routinely screening for this genogroup, along with GI and GII, in order to better understand its distribution, prevalence and role during epidemics, which is probably underestimated.

Molecular identification and genetic analysis of Norovirus genogroups I and II in water environments: comparative analysis of different reverse transcription-PCR assays.

Noroviruses have received increased attention in recent years because their role as etiologic agents in acute gastroenteritis outbreaks is now clearly established. Our inability to grow them in cell culture and the lack of an animal model hinder the characterization of these viruses. More recently, molecular approaches have been used to study the genetic relationships that exist among them. In the present study, environmental samples from seawater, estuarine water, and effluents of sewage treatment plants were analyzed in order to evaluate the role of environmental surface contamination as a possible vehicle for transmission of norovirus genogroups I and II. Novel broad-range reverse transcription-PCR/nested assays targeting the region coding for the RNA-dependent RNA polymerase were developed, amplifying fragments of 516 bp and 687 bp in the nested reactions for genogroups II and I, respectively. The assays were evaluated and compared against widely used published assays. The newly designed assays provide long regions for high-confidence BLAST searches in public databases and therefore are useful diagnostic tools for molecular diagnosis and typing of human noroviruses in clinical and environmental samples, as well as for the study of molecular epidemiology and the evolution of these viruses.

Validation of rt-PCR assays for molecular characterization of porcine teschoviruses and enteroviruses.

Porcine enteroviruses (PEVs) and teschoviruses (PTVs) are described as causative agents of neurological disorders, fertility disorders and dermal lesions of swine. Difficulties in the serological detection of these viruses may lead to a significant underestimation of infections with clinical symptoms. With the recent availability of genome sequence data for all the serotypes, molecular diagnosis is a possibility. The present study describes a new approach to molecular 'serotyping' of PTVs and PEV-B viruses, involving the amplification and sequencing of a genomic fragment of the VP1 coding region. A molecular characterization of Italian entero-teschovirus isolates was performed using a set of previously published and newly designed polymerase chain reaction primers. A total of 33 porcine isolates and 10 reference strains were analysed. Porcine enterovirus-B samples were first diagnosed as positive for enterovirus by amplification of the 5'-non-translated region. Samples were then typed by amplification and sequencing of a portion of the VP1 coding region. Porcine enterovirus-A and PTVs were detected by a published assay in the 5'-NC region that allows them to be differentiated according to the size of amplification product, using the same set of primers. For serotype characterization of PTV, we evaluated four different regions: the N terminus of the capsid protein VP2, the region encoding for RNA-dependent RNA polymerase, and the capsid VP1 and VP4 regions. The newly designed primers in the VP1 region was proved to be broad in range and suitable for serotype assessment and therefore constitute a useful diagnostic tool for molecular diagnosis of porcine teschovirus/enterovirus strains and for the study of molecular epidemiology and evolution of these viruses.

Genetic diversity of bacterial strains isolated from soils, contaminated with polycyclic aromatic hydrocarbons, by 16S rRNA gene sequencing and amplified fragment length polymorphism fingerprinting.

In order to study microbial diversity in a polycyclic aromatic hydrocarbon-impacted soil, 14 bacterial strains were analyzed by 16S rRNA gene sequencing and amplified fragment length polymorphism (AFLP) analysis. Bacterial strains isolated from two different hydrocarbon-polluted sites were identified to the species level by 16S rRNA full-gene sequencing using MicroSeq 16S rRNA gene sequencing. Their genome was subsequently analyzed by high-resolution genotyping with AFLP analysis, in order to monitor species variability and to differentiate closely related strains. Cluster analysis based on AFLP fingerprinting showed intra-specific polymorphism, even among strains with 100% 16S rRNA gene sequence identity. The results show that AFLP is a powerful, highly reproducible and discriminatory tool for revealing genetic relationships in bacterial populations. The ability to differentiate and track related closely microbes is fundamental for studying structure and dynamics of microbial communities in contaminated ecosystems.

Validation of a pXO2-A PCR assay to explore diversity among Italian isolates of Bacillus anthracis strains closely related to the live, attenuated Carbosap vaccine.

Several circulating Bacillus anthracis strains isolated in Italy and belonging to the A1.a cluster, genotype 3 (A1.a-3) are genotypically indistinguishable from Carbosap, a live attenuated vaccine strain, containing both pXO1 and pXO2 plasmids. The genotype was assessed by using eight-locus multilocus variable-number tandem repeat analysis. We describe here the use of a ninth locus able to explore variability among strains that have the same genotype. It is important to be able to genotype the wild isolate of B. anthracis strains from outbreaks of anthrax in areas where Carbosap vaccination of cattle and sheep is common practice. A total of 27 representative field strains isolated in Italy and four vaccinal strains, namely, Carbosap, Sterne, Pasteur I, and Pasteur II, were characterized by a ninth marker, called pXO2-A. Twenty-three field strains were genotype 3 and therefore identical to Carbosap. The marker was in the pXO2 plasmid and is based on the polymorphism of the already-known VX2-3 locus. Detection was obtained by PCR with fluorescence-labeled forward primers in order to produce appropriate fragments for capillary electrophoresis with an ABI 310 genetic analyzer. Genetic relationships showed heterogeneity in all of the examined samples. Interestingly, with respect to genotype 3, samples grouped into eight different subtypes, A to H, and the subtype G, had only two samples indistinguishable from Carbosap. The results of the present study confirm the validity of a hierarchical progressive protocol for discrimination among closely related isolates.

Molecular characterization of low pathogenicity H7N3 avian influenza viruses isolated in Italy.

The complete coding regions of the surface glycoproteins, nucleoprotein (NP), polymerase 2 (PB2), and matrix (M) of A/turkey/214845/02 and A/turkey/220158/99 (H7N3) low pathogenicity avian influenza (LPAI) viruses isolated in October 2002 in Italy were amplified and sequenced to determine the epidemiologic relationships with an A/turkey/Italy/4603/99 (H7N1/4603/99) LPAI virus isolated during the 1999-2001 epizootic in Italy. The hemagglutinin (HA) of H7N3 viruses showed 97.8% nucleotide similarity with A/turkey/Italy/4603/99 (H7N1), and NP, M, and PB2 gene similarities were 93.6%, 98.2%, and 96.2%, respectively. Phylogenetic analyses of HA, PB2, and M genes showed that H7N3 and H7N1 viruses were closely related. Sequence analysis revealed a 23 amino acid deletion in the stalk of the neuraminidase of H7N3 viruses and a unique deletion of amino acid glycine in position 17 in the NP gene of H7N1 virus.

Experimental infection of calves with bovine viral diarrhoea virus type-2 (BVDV-2) isolated from a contaminated vaccine.

A non-cytopathic strain of BVDV-2 was isolated from a batch of live infectious bovine rhinotracheitis (IBR) vaccine, and inoculated intranasally into four 3-month-old calves. Severe signs of disease developed by days 4 and 6 in three of the calves, free of BVDV and antibodies to BVDV, that had been exposed to the virus. These calves survived the acute phase of the infection and progressively recovered. BVDV was consistently isolated, or the respective viral RNA was detected, in the buffy coats from blood samples collected starting from days 2 or 4 up to days 11 or 14 after the experimental infection. Viral RNA was also detected in sera from these infected calves until the presence in the serum of virus neutralizing antibodies was demonstrated. By contrast, the only calf having pre-existing neutralizing antibodies to BVDV at the start of the study was protected from the disease. No virus was detected at any time after experimental inoculation of this calf. Genomic characterization of the BVDV-2 isolated in cell cultures, or detected in sera from the experimentally infected animals, revealed 100%, homology in the nucleotide sequence with the BVDV-2 detected as a contaminant of the live IBR virus vaccine. These findings provided evidence of the infective nature of the contaminant BVDV-2 and of its potential to generate disease outbreaks when inoculated into susceptible animals.

Genetic heterogeneity of bovine viral diarrhoea virus in Italy.

The genetic characteristics, of 38 field isolates of bovine viral diarrhoea virus (BVDV) collected in 1999 from sick or healthy and persistently infected cattle of dairy farms situated in northern Italy, were investigated. A partial 5'-untranslated region (5'-UTR) sequence of each isolate was determined and a phylogenetic analysis was performed. All the isolates were classified as belonging to the BVDV-1 genotype and could be assigned to different BVDV-1 groups, namely BVDV-1b (n = 20), BVDV-1d (n = 6) and BVDV-1e (n = 10). Two remaining isolates could be classified as BVDV-1f and BVDV-1h, respectively. These results provided evidence for genetic heterogeneity of BVDV in Italy, and contribute to a better knowledge of the circulation of BVDV strains, and to their classification.

Genetic diversity at alkB locus in Brucella abortus.

DNA polymorphism of the alkB gene, a DNA repair gene, was assessed by PCR on Brucella abortus biovars 1 (strains 99, S19, 45/20, RB51 and 2308), 3 (Tulya strain), 5 (B3196 strain) and 6 (870 strain). A DNA repetitive element, named IS711, was detected in all studied biovars 1 and its complete nucleotide sequence was determined. We found that the element in alkB gene, bounded by 14 bp imperfect inverted repeats (IRs), is 840 bp long and appears to duplicate a consensus target site, CTAG. Analysing its nucleotide sequence of both forward and reverse strands, more than 10 open reading frames (ORFs) were found. Two potential transposase coding regions were chosen comparing all possible ORFs with the database. Comparing IS711 elements isolated from Brucella species, including both those characterized in our work and the published ones, differences in length and in nucleotide composition were observed among Brucella species, members of the same species and within the same strain. Our results confirm the heterogeneity of IS711 elements in Brucella genus and suggest the possibility to use this element to assess gene and genome diversity and to identify new molecular markers for Brucella species.

Sequence analysis of the genes encoding for the major virulence factors of Bacillus anthracis vaccine strain 'Carbosap'.

AIMS: This study was performed to analyse the molecular characteristics of genes encoding for the major virulence factors in Bacillus anthracis vaccine strain 'Carbosap' compared with the wild B. anthracis strain, to evaluate the basis of attenuation. METHODS AND RESULTS: The molecular characteristics of the B. anthracis 'Carbosap' vaccine strain, used as vaccine in Italy, were analysed in comparison with a B. anthracis virulent strain. Despite the presence of the two virulence plasmids pXO1 and pXO2, the 'Carbosap' strain proved to be protective for cattle. The presence of the regulatory genes atxA and pagR and the gerX operon, known to be involved in the virulence, was verified. In addition, all genes were sequenced. The results showed that no molecular differences between 'Carbosap' and the virulent strain were evident. CONCLUSIONS: The results of this study indicate that the attenuation of the 'Carbosap' vaccine strain is not due to the lack of virulence genes or to modifications occurring on the sequence of these genes. Therefore, other virulence factors, still unknown, could be involved in the pathogenic mechanisms. SIGNIFICANCE AND IMPACT OF THE STUDY: This paper adds new information regarding the molecular characteristics of the vaccine strain 'Carbosap' and highlights the need to better understand the virulence factors involved in the pathogenicity of B. anthracis strains.

Use of polymerase chain reaction to identify Brucella abortus strain RB51 among Brucella field isolates from cattle in Italy.

Brucella abortus strain RB51, a rough mutant of the B. abortus 2308 virulent strain, was recently approved in the United States as the official vaccine for brucellosis in cattle. Following recent evidence of unauthorized use of RB51 vaccine in Italy, where the use of vaccines for brucellosis is no longer allowed, the suitability of an RB51-specific polymerase chain reaction assay for identifying the RB51 strain among Brucella field isolates from cattle in Italy was investigated. The oligonucleotide primers used in this study, belonging to a six-primer cocktail for Brucella species previously described by other authors, allowed the amplification of a 364-base pair (bp) fragment specific for RB51 and its parent strain 2308, and a 498-bp product specific for B. abortus. In addition, unresolved bands ranging from 600 to 700 bp were observed from RB51 strain. Brucella abortus biovars 1, 2 and 4 have only one specific sensitive 498-bp band. The B. abortus biovars 3, 5 and 6 did not give any signal. The 498-bp product from a reference Brucella strain was sequenced and submitted to EMBL with the accession number AJ271969 while the 364-bp fragment from RB51 strain was submitted to EMBL database with accession number AJ271968. The sequence studies confirmed the specificity of the detected fragments. No amplification was obtained by testing DNA from strains antigenically related to Brucella, such as Yersinia enterocolitica O:9, Escherichia coli O:157, Salmonella urbana and Pasteurella multocida. The results of this study indicate that this technique, in combination with specific serological tests, could be a useful diagnostic method to verify the use of RB51 vaccine and can contribute to the creation of a databank of circulating strains.