Development of a multiplex PCR assay for the detection of key genes associated with Pasteurella multocida subspecies

Abstract

The ability to distinguish among the subspecies of Pasteurella multocida isolates is important epidemiologically; however, classification at the subspecies level based on the results of conventional biochemical tests (fermentation of sorbitol and dulcitol) is reportedly not accurate in all cases. Therefore, we developed a rapid, multiplex PCR assay to differentiate among the 3 subspecies of P. multocida. The PCR assay includes the P. multocida species–specific primers KMT1SP6 and KMT1T7 as an internal amplification control, with a newly designed gatD (galactitol-1-phosphate-5-dehydrogenase)-specific primer pair (unique for subsp. gallicida), and primers targeting a 16S rRNA gene region specific for subsp. septica. The subspecies specificity of the PCR was demonstrated by applying the test to a collection of 70 P. multocida isolates, including the Heddleston serovar reference strains; all isolates and strains were assigned correctly. The PCR assay is a sensitive, specific, and highly effective method for the identification of P. multocida subspecies, and an alternative to biochemical test–based differentiation. A possible relationship was noticed between P. multocida subspecies and lipopolysaccharide (LPS) genotype; all but one of the subsp. gallicida strains were isolated only from avian hosts and represented L1 LPS genotype. Subsp. multocida and subsp. septica isolates were classified into 5 and 4 different LPS genotypes, respectively, of which L3 was the only LPS genotype shared between these 2 subspecies.

Pasteurella multocida is a gram-negative bacterium that can infect many avian and mammalian species causing a variety of diseases. ¹⁸ The most common P. multocida–associated diseases are fowl cholera, atrophic rhinitis in pigs, and hemorrhagic septicemia of cattle and buffalo. The organism also causes pneumonia in ruminants, horses, and pigs, and is responsible for the development of pasteurellosis in rabbits. ¹⁸ Humans may also be infected, with wound infection injuries (bites or scratches) from pets ⁹ being the most frequent clinical manifestation of the disease. ¹⁸

Within P. multocida, 3 subspecies—P. multocida subsp. multocida, P. multocida subsp. septica, and P. multocida subsp. gallicida—have been distinguished by using DNA-DNA hybridization. ¹¹ The subspecies also differ in their fermentation properties. Isolates of subsp. multocida can utilize sorbitol, but not dulcitol. In contrast, isolates of subsp. septica differ from subsp. multocida by not fermenting sorbitol; subsp. gallicida isolates differ by fermenting dulcitol. ¹¹ Subspeciation of P. multocida strains is important for epidemiologic investigations. In previous studies, subsp. multocida was the subspecies found most frequently in pigs and poultry^1,4; isolates of feline and canine origin generally belonged to subsp. septica. ¹⁰

Comparative sequence analysis of the 16S rRNA gene and key genes of basic cellular functions support the existence of only 2 taxonomic units within the species, a complex group comprising subsp. multocida/gallicida and subsp. septica.^2,10,13,14 This phylogenetic divergence in the 16S rRNA gene sequence allows the identification of subsp. septica with molecular methods. ¹⁰ Subspecies classification based on the results of conventional biochemical tests (sorbitol and dulcitol) is reportedly not accurate in all cases,^3,10,14,17 and 16S ribosomal RNA gene–based molecular methods are recommended for precise differentiation in addition to biochemical tests. ¹⁰ A PCR fingerprinting method has also been used to separate P. multocida subspecies. ⁵ However, inter-laboratory standardization and reproducibility of PCR-based fingerprinting methods are difficult.

Fermentation of dulcitol occurs only in subsp. gallicida.^1,4,11 The genetic background of dulcitol (galactitol) utilization was studied in detail in Salmonella and revealed that the genes responsible for this metabolism are part of a gene cluster that encodes a phosphotransferase system. ¹² This galactitol-specific phosphotransferase system was also identified in P. multocida reference strain X73 (subsp. gallicida), but not in P1059 (subsp. multocida). ⁸ A study found that genes responsible for galactitol metabolism are characteristic only of strains isolated from fowl cholera, and are not carried by strains isolated from other host species and diseases. ⁷ Another study on P. multocida pangenomic diversity also identified an association between subspecies and genes responsible for sorbitol and dulcitol metabolism. ¹⁵ We aimed to design a multiplex PCR scheme based on the detection of key genes associated with the P. multocida subspecies.

We included reference strains representing Heddleston serovars 1–16 (except 15), and a total of 55 field isolates of P. multocida from various host species (Table 1). The isolates were cultured on Columbia agar plates (Becton Dickinson) supplemented with 5% sheep blood under aerobic conditions at 37°C for 24 h. Their identity was confirmed by a species-specific PCR assay. ¹⁶ Dulcitol and sorbitol tests were performed as described before.^1,4 The lipopolysaccharide (LPS) genotype was determined using the LPS multiplex PCR. ⁶

Table 1.

Origin of Pasteurella multocida strains and isolates used in our study.

Subspecies	Strain ID*/No. of field isolates	Origin		Dulcitol	Sorbitol	PCR result			LPS genotype
		Country	Host			KMT1T7–KMT1SP6	gatD F–gatD R	S152 F–S991 R
gallicida	X73 (1)	Ref.	Avian	+	+	+	+	–	L1
	P2225 (14)	Ref.	Bovine	+	+	+	+	–	L1
	2	F	Avian	±†	+	+	+	–	L1
	14	F	Avian	+	+	+	+	–	L1
multocida	M1404 (2)	Ref.	Bison	–	+	+	–	–	L2
	P1059 (3)	Ref.	Avian	–	+	+	–	–	L3
	P1662 (4)	Ref.	Avian	–	+	+	–	–	L3
	P1702 (5)	Ref.	Avian	–	+	+	–	–	L2
	P2095 (9)	Ref.	Avian	–	+	+	–	–	L5
	P2100 (10)	Ref.	Avian	–	+	+	–	–	L6
	P2723 (16)	Ref.	Avian	–	+	+	–	–	L8
	P1573 (12)	Ref.	Human	–	+	+	–	–	L6
	P903 (11)	Ref.	Swine	–	+	+	–	–	L6
	1	F	Avian	–	+	+	–	–	L2
	1	F	Avian	–	+	+	–	–	L3
	2	F	Avian	–	+	+	–	–	L3
	1	F	Avian	–	+	+	–	–	L3
	1	H	Bovine	–	+	+	–	–	L2
	2	H	Bovine	–	+	+	–	–	L3
	2	H	Bovine	–	+	+	–	–	L6
	2	H	Caprine	–	+	+	–	–	L3
	1	H	Sheep	–	+	+	–	–	L6
	3	H	Sheep	–	+	+	–	–	L3
	1	H	Swine	–	+	+	–	–	L2
	7	H	Swine	–	+	+	–	–	L3
septica	P2192 (6)	Ref.	Avian	–	+	+	–	+	L4
	P1997 (7)	Ref.	Avian	–	+	+	–	+	L4
	P1581 (8)	Ref.	Avian	–	–	+	–	+	L7
	P1591 (13)	Ref.	Human	–	–	+	–	+	L7
	2	F	Avian	–	+	+	–	+	L3
	1	F	Avian	–	+	+	–	+	L1
	9	F	Avian	–	+	+	–	+	L1
	3	H	feline	–	+	+	–	+	L1

F = France; H = Hungary; LPS = lipopolysaccharide; ref. = Heddleston serovar reference strain.

*

 The serovar of the Heddleston reference strains is indicated in parentheses.

†

 Weak positive color reaction in the biochemical test.

Primers specific for the galactitol-1-phosphate-5-dehydrogenase gene (gatD) were designed (Table 2) from the complete genome sequence of X73 reference strain (GenBank NZ_CM001580). Primers specific for the 16S rRNA gene of P. multocida subsp. septica were designed from the sequences of reference strains available in GenBank (AF294410–AF294423). The primers were synthesized by Sigma-Aldrich. Primers gatD F and gatD R were designed to produce a 527-bp PCR fragment from the gatD gene, whereas primers S152 F and S991 R amplified an 840-bp fragment from the 16S rRNA gene of subsp. septica isolates. Amplicon sizes were calculated from the DNA sequences of subsp. gallicida and subsp. septica reference strains.

Table 2.

Sequences of the oligonucleotides used in the Pasteurella multocida subspecies multiplex PCR assay.

Subspecies	Gene	Name	Sequence (5′–3′)	Amplimer size (bp)	Position	Reference
gallicida	gatD	gatD F	TTGGTGGTATGTGCACCGTT	527	238–257*	Our study
	gatD	gatD R	TTAGCACGTGGACCAGCAAT		764–745*	Our study
septica	16S rRNA	S152 F	GCGTAATCTCTGAGGAGTAAAGGG	840	152–175†	Our study
	16S rRNA	S991 R	CATCTCTGCAAGCTTCTCTGG		991–971†	Our study
All subspecies	KMT1	KMT1T7	ATCCGCTATTTACCCAGTGG	460		16
All subspecies	KMT1	KMT1SP6	GCTGTAAACGAACTCGCCAC			16

*

 These sequences match the coding strand of the gatD gene of P. multocida X73 strain with GenBank NZ_CM001580.

†

 These sequences match the coding strand of the 16S rRNA gene of P. multocida subsp. septica strain CCUG 17977 with GenBank NR_041810.

Lysates of pure cultures were prepared for all strains used in our study. A loopful of cultured bacterial growth was suspended in 100 µL of PCR-grade water (VWR) and heated in a thermocycler for 20 min at 99°C. Cellular debris was pelleted by centrifugation, and the supernatant was used as DNA template for PCR amplification. The lysates were stored at –20°C until use.

PCR was performed in a total volume of 25 µL containing 100 ng of template DNA, 2.5 µL of 10× DreamTaq buffer, 1.2 µL of 25 mM MgCl₂, 0.5 µL of 10 mM dNTPs, 1 µL of forward and reverse primers (10 pmol/µL), and 1 U of DreamTaq polymerase (Thermo Scientific). Reactions were performed in a thermocycler (C1000; Bio-Rad). DNA was amplified for 25 cycles of denaturation at 95°C for 30 s, annealing at 58°C for 30 s, and extension at 72°C for 1 min; an initial denaturation step (95°C, 10 min) and a final extension step (72°C, 7 min) were also performed. Each reaction mixture was analyzed by electrophoresis in 1.5% agarose gel (SeaKem; FMC) stained with GelRed (Biotium) and visualized under UV light. The results were confirmed by sequencing of PCR products. Sequencing was performed by Macrogen Europe. Nucleotide sequences were analyzed using Geneious Prime software (v.2020.2.4, https://www.geneious.com) and confirmed by BLAST (https://blast.ncbi.nlm.nih.gov/Blast.cgi). The GenBank accessions for sequences obtained in our study are MZ476034–MZ476037.

Following sequence determination and analysis, a multiplex PCR assay was developed that contained the P. multocida species–specific primers KMT1SP6 and KMT1T7 as an internal amplification control, ¹⁶ with the newly designed gatD gene–specific primer pair (unique for subsp. gallicida), and S152 F and S991 R primers (targeting 16S rRNA gene regions specific for subsp. septica). Amplification and electrophoresis conditions were the same as described above. Amplimer lengths allowed clear size discrimination (Fig. 1).

Figure 1.

Multiplex subspecies PCR typing system for Pasteurella multocida.

To test the sensitivity of our assay, 10-fold dilutions of genomic DNA (10⁶–10⁰) of reference strains X73 (subsp. gallicida) and P1581 (subsp. septica) were used. DNA was quantified (NanoDrop spectrophotometer; Thermo Scientific), and template copy numbers per reaction were calculated using an online calculator (http://cels.uri.edu/gsc/cndna.html). The subspecies specificity of the PCR assay was determined by applying the test to a collection of 70 P. multocida isolates, including 15 of the 16 Heddleston serovar reference strains. For all strains or isolates belonging to subsp. gallicida or subsp. septica, a subspecies-specific fragment of the expected size and a 460-bp species-specific fragment were amplified. For all strains or isolates of subsp. multocida, only the species-specific fragment was produced (Fig. 1). Sensitivities of the assays were 10⁵ genomic equivalents per reaction (data not shown).

All P. multocida field isolates typed as subsp. gallicida in the biochemical tests (dulcitol and sorbitol positive) were allocated to the corresponding subspecies by the PCR test. In accordance with previous studies,^2,10,13,14 the ability to utilize sorbitol was found among strains belonging to subsp. septica. Thus, we confirmed the proposition that the 16S rRNA gene could be used to identify subsp. septica instead of sorbitol and dulcitol tests. Isolates belonging to subsp. multocida uniformly gave their characteristic phenotypic reactions (fermentation of sorbitol, but not dulcitol). Almost all subsp. gallicida strains had been isolated from avian hosts (with the only exception of P2225, the serovar 14 reference strain of bovine origin) and classified exclusively into L1 LPS outer core biosynthesis loci. On the other hand, subsp. multocida and subsp. septica isolates were classified into 5 and 4 different LPS genotypes, respectively. Interestingly, L3 was the only LPS genotype shared between the latter 2 subspecies.

The identification of P. multocida subspecies has been the subject of research for decades. Our results indicate that our developed PCR test is a sensitive, specific, and highly effective method for the differentiation of P. multocida subspecies. Problems related to the traditional biochemical tests (i.e., the existence of sorbitol-positive subsp. septica strains reported by others^10,17) are avoided. This method could be an alternative to subspecies differentiation by biochemical tests in epidemiologic investigations.