Comparative genome analysis of an avirulent and two virulent strains of avian Pasteurella multocida reveals candidate genes involved in fitness and pathogenicity
BMC Microbiology volume 13, Article number: 106 (2013)
Pasteurella multocida is the etiologic agent of fowl cholera, a highly contagious and severe disease of poultry causing significant mortality and morbidity throughout the world. All types of poultry are susceptible to fowl cholera. Turkeys are most susceptible to the peracute/acute forms of the disease while chickens are most susceptible to the acute and chronic forms of the disease. The whole genome of the Pm70 strain of P. multocida was sequenced and annotated in 2001. The Pm70 strain is not virulent to chickens and turkeys. In contrast, strains X73 and P1059 are highly virulent to turkeys, chickens, and other poultry species. In this study, we sequenced the genomes of P. multocida strains X73 and P1059 and undertook a detailed comparative genome analysis with the avirulent Pm70 strain. The goal of this study was to identify candidate genes in the virulent strains that may be involved in pathogenicity of fowl cholera disease.
Comparison of virulent versus avirulent avian P. multocida genomes revealed 336 unique genes among the P1059 and/or X73 genomes compared to strain Pm70. Genes of interest within this subset included those encoding an L-fucose transport and utilization system, several novel sugar transport systems, and several novel hemagglutinins including one designated PfhB4. Additionally, substantial amino acid variation was observed in many core outer membrane proteins and single nucleotide polymorphism analysis confirmed a higher dN/dS ratio within proteins localized to the outer membrane.
Comparative analyses of highly virulent versus avirulent avian P. multocida identified a number of genomic differences that may shed light on the ability of highly virulent strains to cause disease in the avian host, including those that could be associated with enhanced virulence or fitness.
Avian pasteurellosis, also known as fowl cholera is a highly contagious, systemic, and severe disease affecting wild and domestic birds frequently resulting in high mortality and morbidity. The disease is of major economic importance throughout the world in areas of domestic poultry production [1–3]. The causative agent of fowl cholera is Pasteurella multocida, a Gram-negative bacterium. Carter [4, 5] identified five capsular types of P. multocida based on differences in capsular antigens and designated them as A, B, D, E, and F serogroups. Heddleston and co-workers classified the bacterium into 16 somatic types based on differences in the lipopolysaccharide antigens . In 1981, a standard system for identifying serotypes of P. multocida was recommended that combined both the Carter capsular typing and Heddleston somatic typing systems to designate serotypes  and a serotype is designated by its capsular type followed by the somatic type. Using this system, the most common serotypes causing fowl cholera in the United States are A:1, A:3, and A:3.4 . While there are no indications that any particular serotype is more or less virulent than others the virulence of avian isolates of most common serotypes appears to vary considerably .
Fowl cholera disease can occur in peracute/acute and subacute/chronic forms . All types of poultry are susceptible to the disease, although among them turkeys, pheasants and partridges are highly susceptible to peracute/acute forms of disease whereas chickens are relatively more resistant . In chickens, the most common forms of the disease are acute and chronic. In peracute/acute disease there is sudden death due to terminal – stage bacteremia and endotoxic shock [1, 3]. Signs of acute cholera have been reproduced by injection of endotoxin from P. multocida[12–14]. Post-mortem findings are dominated by general septicemic lesions. [1, 2]. In chronic disease, signs are principally due to localized infections of leg or wing joints, comb, wattles and subcutaneous tissue of the head [2, 10]. The completed genome of P. multocida strain Pm70 has been available for over eleven years  and has greatly facilitated subsequent genomic-based approaches towards better understanding the underlying genetic mechanisms related to virulence and fitness. This complete genome sequence has been used in the study of specific enzymes , microarray analyses of differentially expressed genes [17–20], proteomic analyses [21, 22], study of virulence factors [16, 23–25], reverse vaccinology approaches , and as a reference for assembly and comparison to other genomes. While the Pm70 genome sequence has been a great asset in our studies, progress has been modest in the identification and understanding of P. multocida virulence . Even today, very little is known about the totality of the mechanisms behind P. multocida’s ability to cause disease. The Pm70 strain was isolated from the oviduct of a layer chicken in 1976 from Texas (personal communication- RE. Briggs). This strain belongs to serotype F:3  and not A:3 as reported earlier , is avirulent and does not cause experimental fowl cholera disease in chickens . In contrast, other strains of P. multocida have been isolated, such as strains X73 and the P1059, that are highly virulent to chickens, turkeys, and other poultry species [29, 30].
Additional P. multocida strains of bovine, avian, and porcine origin have recently been sequenced, which was the subject of a recent comparative review . The authors noted, based on the nine genomes sequenced to date, there was “no clear correlation between phylogenetic relatedness and host predilection or disease”. Information is sparse on the location and characterization of the genes responsible for differences in virulence of avian and other P. multocida. In another recent review, Boyce et al.  speculated that the combination of additional P. multocida genome sequences and advances in our ability to genetically manipulate the organism will facilitate major advances in our understanding of disease pathogenesis. To that end, we undertook a detailed comparative genome analysis of two virulent strains (X73 and P1059) and avirulent strain Pm70 of P. multocida. The goal of this study is to enable narrowed identification of a repertoire of unique genes present in the highly pathogenic avian strains that may play a role in virulence. This information will also facilitate the design of improved modified live vaccine candidates with defined mutations that can be evaluated as immunoprophylactic agent(s) to control P. multocida-caused disease in avian and other host species.
The strains sequenced in this study included P. multocida strains P1059 (ATCC# 15742) and X73 (ATCC# 11039). Strain P1059 is a well characterized pathogenic strain isolated from the liver of a turkey that died of fowl cholera . Strain X73 is also a well characterized pathogenic strain isolated from the liver of a chicken that died of fowl cholera . For comparative purposes, the avirulent Pm70 strain was used . There are several reasons why we selected strains P1059 and X73 in this study. First, both strains are highly virulent to chickens, turkeys and other poultry species. Second, they are of different serotypes (P1059 = A:3; and X73 = A:1) and different immunologic types . Thirdly, they are reference serotype strains that are readily available to investigators and there is abundant literature on the biology of these two strains [1, 11, 30, 33–35].
Genome sequencing and annotation
Sequencing of strains P1059 and X73 was performed using 454 Life Sciences pyrosequencing at the National Animal Disease Center in Ames, Iowa. The following data sets were generated for each strain: GS- FLX, with 270,010 shotgun reads of average length 240 bp yielding 64,827,159 bp for P1059; and GS-FLX, with 227,030 shotgun reads of average length of 240 bp, yielding 54,398,540 bp for X73. Reads were de novo assembled into scaffolds using Newbler 2.3 . The draft sequences of these genomes are deposited under the following accession numbers: P1059 [Genbank:AMBQ01000000] and X73 [Genbank:AMBP01000000].
Annotation of P1059 and X73 was performed using publicly available tools. Putative coding regions were predicted using GeneMarkS . Gene function was assigned using HMMER3 against Pfam-A 24.0, RPS-BLASTp against CDD, and BLASTp against all microbial proteins [38, 39]. tRNA genes were identified using tRNAscan-SE . rRNA genes were identified using RNAmmer . For analysis of the shared and unique proteins in the P. multocida genomes sequenced, BlastP was used with a similarity cutoff of 90% identity over 90% of the protein as an arbitrary designation for similar versus dissimilar proteins. For genomic island analysis, whole genome alignments were performed using MAUVE to identify regions present in strains P1059 and X73 but absent from strain Pm70 . Linear and circular genomic maps were generated using XPlasMap and Circos . Single nucleotide polymorphism (SNP) analysis was performed using SNPeff .
Results and discussion
Overview of the P. multocidaP1059 and X73 genomes
A total of 270,010 reads were used to draft assemble strain P1059, resulting in a single scaffold of 27 large contigs (> 500 bp) of approximately 27-fold coverage and an estimated genome size of 2.4 Mb. A total of 227,030 reads were used to draft assemble strain X73, resulting in 17 large contigs (> 500 bp) of approximately 23-fold coverage and an estimated genome size of approximately 2.4 Mb. No plasmids were identified in either strain sequenced. The contigs generated were then aligned to strain Pm70 to generate collinear draft sequences and subsequently compare the three avian source genomes.
Unique regions of virulent avian P. multocida
The draft genomes of strains P1059 and X73 were found to contain 2,144 and 2,085 predicted proteins, respectively. Along with strain Pm70, the genomes all contained 51 tRNA-carrying genes and 4 rRNA-carrying operons. The genomes of the three avian P. multocida strains contained a remarkably high number of shared proteins (1,848), which comprised 86.2-90.7% of the predicted proteins of the three avian strains using a BlastP similarity cut-off of 90% (Figure 1). Compared to strain Pm70, a total of 336 unique proteins were identified in either strains P1059 or X73, of which 61 were contained within both genomes (Table 1). Most of the 61 shared proteins were small predicted proteins of unknown function and located individually throughout the P. multocida genome that could be attributed to differences in annotation approaches (Figure 2). Also, most of the predicted proteins identified were present in one or more sequenced P. multocida from the NCBI database that were not from avian hosts. However, one noteworthy region of difference shared by P1059 and X73, but absent from Pm70 and other strains of non-avian source, was located between the core genes deoC and rfaD in both P1059 and X73 (P1059 – 01496 to 01503; X73 – 01400 to 01407). This region contained ten predicted proteins with similarity to systems involved in the transport and utilization of L-fucose. L-fucose is an important component of host mucin and has shown to be a chemoattractant for certain bacterial species, such as Campylobacter jejuni. Moreover, the ability to utilize L-fucose by C. jejuni has been shown to confer a fitness advantage for avian strains in low nutrient environments such as the respiratory tract [45, 46]. Comparison of available P. multocida sequences suggests that the presence of this region may be a defining feature of pathogenic avian-source P. multocida strains, as it was present in P1059, X73, and P. multocida subsp. gallicida strain Anand1 isolated from a chicken in India  but absent from strains Pm70, pathogenic bovine-source strain 36950  and pathogenic swine source strains 3480 and HN06 . Other studies have demonstrated an ability of avian-source P. multocida to ferment L-fucose, further suggesting that the majority of avian-source P. multocida strains harbor this system [9, 33, 50]. Other bacteria inhabiting the respiratory tracts of poultry have been identified to utilize L-fucose, such as Gallibacterium anatis, suggesting that such capabilities may be advantageous for respiratory bacterial pathogens of poultry . Such systems could play a role in increased fitness and/or virulence capability of strains P1059 and X73 in the avian host.
Twelve proteins were also identified that were present in both strains P1059 and X73 at greater than 90% amino acid similarity, but at less than 90% similarity in strain Pm70 (Table 2). Among the twelve proteins identified were several membrane-associated proteins, including LspB, PfhB3, Opa, and SprT. The presence of divergent protein sequences that are membrane-associated is suggestive of adaptation of P. multocida strains towards particular hosts.
There were also predicted proteins identified as unique to strains P1059 (148 total) and X73 (127 total) compared to strain Pm70. Many of these proteins were again of unknown function and/or associated with prophage-like elements (Additional file 1: Table S1 and Additional file 2: Table S2). However, some systems unique to each strain were noteworthy. In strain P1059, one unique region contained six genes predicted as involved in the transport and modification of citrate, and the conversion of citrate to oxaloacetate via citrate lyase (00080 to 00085). This system was absent in all other sequenced P. multocida genomes. The conversion of citrate to oxaloacetate is linked to citrate fermentation. Also unique to strain P1059, but present in strains 36950, 3480, and HN06, are four genes involved in xylose ABC transport system with a transcriptional repressor (01538 to 01541). Present in strains X73 and 36950 was a putative toxin-antitoxin system similar to the HipAB systems (genes 02005 and 02006). Finally, genes for several novel proteins with similarity to the previously described Pfh-type filamentous hemagglutinins were identified in strains P1059 and X73. Strain P1059 contained a novel predicted filamentous hemagglutinin (designated PfhB4 – gene # 00523) that shares similarity with PfhB1 and PfhB2 from P. multocida. PfhB4 has conserved domains related to hemagglutination activity, two-partner secretion, hemagglutinin repeats, and toxicity. PfhB4 is present only in strains P1059, HN06, and 3480 (Figure 3). It is adjacent to a putative hemolysin secretion/activation protein (gene # 00522) and thus appears to represent a novel two-partner system involved in hemagglutination. PfhB2 of strain P1059 has been shown to play an important role in either colonization or invasion in the turkey model . Also, vaccination with recombinant P1059 PfhB2 peptides cross protected turkeys against an X73 challenge . PfhB2 was present in strain Pm70, P1059, and X73, but was only 90% similar in the latter two as compared to Pm70. Overall, the presence of unique genes/systems related to metabolism and adhesion could provide strains such as P1059 with additional tools for increased fitness leading to higher virulence.
Of the 127 unique proteins identified in strain X73 were five genes for a galactitol-specific phosphotransferase and utilization system (00310 to 00316), only present in strain X73; three genes for a TRAP dicarboxylate transporter system (01441 to 01443), also present in strain 36950; and six genes for a novel simple sugar D-allose transport and utilization systems (00951 to 00956), only present in strain X73. Such systems could again provide additional means of energy production in a resource-limited environment.
Known virulence factors and antigens
Comparisons were performed for several known virulence factors and outer membrane proteins that are important for P. multocida pathogenesis, functionality, and vaccine development . These comparisons revealed some noteworthy aspects relative to their presence and evolution in P. multocida. For example, the hemoglobin receptors hgbA and hgbB were present in all sequenced P. multocida genomes, but are significantly different in their amino acid similarities (Table 3). HgbA and HgbB have been shown to exhibit hemoglobin binding properties [53, 54]. Their incomplete distribution reported in previous studies could be attributed to genetic variation rather than complete absence of these genes . The outer membrane porins ompH1 and ompH2 were also present in all sequenced strains, with ompH2 more highly conserved than ompH1 with respect to amino acid similarity. Furthermore, a third outer membrane porin ompH3 was present in all sequenced strains except strain X73, but was highly conserved within these strains. The ptfA gene, encoding a type 4 fimbrial subunit, was highly conserved in all sequenced strains, as was comE encoding a fibronectin-binding protein. The pfhB1 gene, encoding a filamentous hemagglutinin protein, was present in strains Pm70, P1059, X73, and 3480. PfhB1 was highly conserved among these strains. PfhB2, a second filamentous hemagglutinin, was present in strains Pm70, P1059, X73, and 36950. This protein was more variable in amino acid sequence among these strains (Figure 3). Two other genes encoding filamentous hemaggultinins, pfhB3 and pfhB4, were absent in strain Pm70, with pfhB3 present in strains P1059, X73, and 36950, and pfhB4 present in strains P1059, HN06, and 3480. Finally, lipoproteins plpP, plpB, and plpD were present in all sequenced strains, and all were highly conserved except plpP, whose product shared only 82-98% amino acid similarity between strains.
Single nucleotide polymorphisms
The three avian source P. multocida genomes were also compared for SNPs within the conserved regions of their genomes using MAUVE , and the SNPs were analyzed for their coding effects using SNPeff  (Table 4). A total of 31,021 SNPs were identified between strains Pm70 and P1059, and 26,705 SNPs were identified between strains Pm70 and X73. The density of SNPs varied considerably across the P. multocida genome, with some regions containing a much higher density of SNPs than the rest of the core genome (Figure 4). This suggests that some regions of the genome are under diversifying selection, while the majority of the genome is under neutral or purifying selection. The ratio between non-synonymous to synonymous substitutions (dN/dS) is commonly used as a measure of purifying versus diversifying selection . The overall dN/dS ratios of all coding regions of strains P1059 and X73 compared to strain Pm70 were 0.40 and 0.38, respectively. Proteins were then divided into groups based upon predicted subcellular localization of each protein using PSORT-B version 3.0. Using this approach, the dN/dS ratios varied considerably, with higher ratios (0.76-0.93) found within proteins predicted as extracellular or outer membrane . Amongst specific outer membrane proteins, the highest dN/dS ratios were observed within PfhB2, HgbA, HemR, pm0591 (a secreted effector protein), pm0803 (an iron-regulated outer membrane protein), TadD-F (pilus assembly proteins), RcpB-C (pilus assembly proteins), and PlpP. The higher dN/dS ratios observed among this subset of extracellular and outer membrane proteins is suggestive that they are under diversifying selection due to interactions with the host immune system, although further analyses would be required to confirm this observation.
The Heddleston somatic typing system classifies P. multocida into 16 somatic types based on antigenic differences in the lipopolysaccharide (LPS) . Good progress has been made in understanding the structural basis for the LPS typing scheme. The genes and the transferases required for the biosynthesis of the somatic type-specific outer core region of the LPS has been identified in strains of P. multocida strains representing various somatic types [24, 58–63]. Since endotoxin (LPS) is a key virulence factor in P. multocida, we examined each gene involved in LPS biosynthesis in the X73 and P1059 strains and compared with the Pm70 strain. All three strains produced two glycoforms simultaneously, termed glycoforms A and B. Both X73 and P1059 contained the inner core biosynthetic complement of genes, including kdtA (P1059-01455; X73- 01363), hptA (opsX; P1059-02017; X73- 01921), kdkA (P1059-01451; X73-01359), hptC (rfaF; P1059-02018 ; X73-01922), hptD (P1059-01443; X73-01351 ) and gctA (P1059-01456; X73-01364). The gene that encodes for the enzyme which catalyzes the attachment of phosphoethanolamine to L-α-D Heptose −11 (Pm70-pm0223) was present only in strains P1059 and Pm70. There appeared to be some variation in the hptD gene between Pm70 and the X73 and P1059 strains although it was generally conserved between strains. Linking the inner core to the outer core is the hptE gene, present in both X73 and P1059 (X73-01185; P1059-01293). The outer core structure expressed by X73, P1059 and Pm70 strains are structurally distinct and distal part of the molecule because in all three strains a polymeric O antigen was absent. The X73 strain but not P1059 and Pm70 express an outer core oligosaccharide that contains two terminal galactose residues, with phosphocholine (PCho). Present in X73 but absent from Pm70 and P1059 were the outer core biosynthetic genes involved in phosphocholine (PCho) biosynthesis genes for somatic type 1. As reported previously , these genes include pcgA (X73-01180), pcgB (X73-01182), pcgC (X73-01181), and pcgD (X73-01183) as well as gatA (X73-01184). X73 attaches a phosphoethanolamine (PEtn) residue to the terminal galactose. Studies have shown  that PCho on the LPS is important for virulence of X73 strain to chickens. However, a clear role for PEtn has not been defined. Present in the outer core of Pm70 and P1059, but absent in X73, were the biosynthetic genes for somatic type 3. These genes include losA (Pm70-Pm1143; P1059-01292); (Pm70-Pm1138; P1059-01287); (Pm70-Pm1139; P1059-01288); (Pm70-Pm1140; P1059-01289); and (Pm70- Pm1141; P1059-01290).
In summary, comparative analyses of highly virulent versus avirulent P. multocida identified a number of genomic differences that may shed light on the ability of highly virulent strains to cause disease in the avian host. Most of the differences observed involved the presence of additional systems in virulent avian-source strains P1059 and/or X73 that appear to play metabolic roles. Such systems might enhance the fitness of these strains in the avian extraintestinal compartment, but without experimental evidence this is purely a speculative observation. This work does, however, underscore the need to utilize such genomic data towards targeted molecular approaches to better understand the role of horizontal gene transfer in the pathogenesis of this organism. Also, it is evident given the high degree of large sequence and single nucleotide polymorphisms in P. multocida that focused studies need to be conducted to appreciate adaptation of these strains to their respective hosts.
Christensen JP, Bisgaard M: Avian pasteurellosis: taxonomy of the organisms involved and aspects of pathogenesis. Avian Path. 1997, 26: 461-483. 10.1080/03079459708419228.
Christenson JP, Bisgaard M: Fowl Cholera. Rev Sci Tech. 2000, 19: 626-637.
Wilkie IW, Harper M, Boyce JD, Adler B: Pasteurella multocida: Diseases and Pathogenesis. Curr Top Microbiol Immunol. 2012, 361: 1-22. 10.1007/82_2012_216.
Carter GR: Studies on Pasteurella multocida. A hemagglutination test for the identification of serological types. Amer J Vet Res. 1955, 16: 481-484.
Carter GR: A new serological type of Pasteurella multocida from Central Africa. Vet Rec. 1961, 73: 1052-
Heddleston KL, Gallagher JE, Rebers PA: Fowl cholera: Gel diffusion precipitin test for serotyping Pasteurella multocida from avian species. Avian Dis. 1972, 16: 925-936. 10.2307/1588773.
Carter GR, Chengappa MM: Recommendations for a standard system of designating serotypes of Pasteurella multocida. Proceedings of the 24th Amer. Assoc. Veterinary Laboratory Diagnosticians. 1981, 24: 37-42.
Rhodes KR, Rimler RB: Somatic serotypes of Pasteurella multocida strains isolated from avian hosts (1976–1988). Avian Dis. 1990, 34: 193-195. 10.2307/1591351.
Lee MD, Wooley RE, Glisson JR, Brown J: Comparison of Pasteurella multocida serotype 3,4 isolates from turkeys with fowl cholera. Avian Dis. 1988, 32: 501-508. 10.2307/1590918.
Webster LT: The epidemiology of fowl cholera. J Exp Med. 1930, 51: 219-223. 10.1084/jem.51.2.219.
Petersen KD, Christensen JP, Permin A, Bisgaard M: Virulence of Pasteurella multocida subsp. multocida isolated from outbreaks of fowl cholera in wild birds for domestic poultry and game birds. Avian Pathol. 2001, 30: 27-31. 10.1080/03079450020023168.
Heddleston KL, Rebers PA: Properties of free endotoxin from Pasteurella multocida. Am J Vet Res. 1975, 36: 573-574.
Rhodes KR, Rimler RB: Effect of Pasteurella multocida endotoxins on turkey poults. Avian Dis. 1987, 31: 523-526. 10.2307/1590734.
Heddleston KL, Rebers PA, Ritchie AE: Immunizing and toxic properties of particulate antigens from two immunogenic types of Pasteurella multocida of avian origin. J Immunol. 1966, 96: 124-133.
May BJ, Zhang Q, Li LL, Paustian ML, Whitman TS, Kapur V: Complete genome sequence of Pasteurella multocida Pm70. Proc Natl Acad Sci USA. 2001, 98: 3460-3465. 10.1073/pnas.051634598.
Steenbergen SM, Lichtensteiger CA, Caughlan R, Garfinkle J, Fuller TE, Vimr ER: Sialic acid metabolism and systemic pasteurellosis. Infect Immun. 2005, 73: 1284-1294. 10.1128/IAI.73.3.1284-1294.2005.
Steen JA, Steen JA, Harrison P, Seemann T, Wilkie I, Harper M, Adler B, Boyce JD: Fis is essential for capsule production in Pasteurella multocida and regulates expression of other important virulence factors. PLoS Pathog. 2010, 6: e1000750-10.1371/journal.ppat.1000750.
Nanduri B, Shack LA, Burgess SC, Lawrence ML: The transcriptional response of Pasteurella multocida to three classes of antibiotics. BMC Genomics. 2009, 14 (10 Suppl 2): S4-
Boyce JD, Wilkie L, Harper M, Paustian ML, Kapur V, Adler B: Genomic scale analysis of Pasteurella multocida gene expression during growth within liver tissue of chickens with fowl cholera. Microbes Infect. 2004, 6: 290-298. 10.1016/j.micinf.2003.12.004.
Paustian ML, May BJ, Kapur V: Transcriptional response of Pasteurella multocida to nutrient limitation. J Bacteriol. 2002, 184: 3734-3739. 10.1128/JB.184.13.3734-3739.2002.
Nanduri B, Lawrence ML, Peddinti DS, Burgess SC: Effects of subminimum inhibitory concentrations of antibiotics on the Pasteurella multocida proteome: a systems approach. Comp Funct Genomics. 2008, 10.1155/2008/254836.
E-Komon T, Burchmore R, Herzyk P, Davies R: Predicting the outer membrane proteome of Pasteurella multocida based on consensus prediction enhanced by results integration and manual confirmation. BMC Bioinformatics. 2012, 13: 63-80. 10.1186/1471-2105-13-63.
Harper M, Cox A, St Michael F, Parnas H, Wilkie I, Blackall PJ, Adler B, Boyce JD: Decoration of Pasteurella multocida lipopolysaccharide with phosphocholine is important for virulence. J Bacteriol. 2007, 189: 7384-7391. 10.1128/JB.00948-07.
St Michael F, Vinogradov E, Li J, Cox AD: Structural analysis of the lipopolysaccharide from Pasteurella multocida genome strain Pm70 and identification of the putative lipopolysaccharide glycosyltransferases. Glycobiology. 2005, 15: 323-333.
Bosch M, Garrido ME, de Rozas AM P, Badiiola I, Barbe J, Llagostera M: Pasteurella multocida contains multiple immunogenic haemin- and haemoglobin-binding proteins. Vet Microbiol. 2004, 99: 102-112.
Hatfaludi T, Al-Hasani K, Gong L, Boyce JD, Ford M, Wilkie IW, Quinsey N, Dunstone MA, Hoke DE, Adler B: Screening of 71 P. multocida proteins for protective efficacy in a fowl cholera infection model and characterization of the protective antigen PlpE. PLoS One. 2012, 7: e39973-10.1371/journal.pone.0039973.
Ewers C, Becker AL, Bethe A, Kiebling S, Filter M, Wieler LH: Virulence genotype of Pasteurella multocida strains isolated from different hosts with various disease status. Vet Microbiol. 2006, 114: 304-317. 10.1016/j.vetmic.2005.12.012.
Al-Hasani K, Boyce J, McCarl VP, Bottornley S, Wilkie I, Adler B: Identification of novel immunogens in Pasteurella multocida. Microb Cell Fact. 2007, 6: 1-5. 10.1186/1475-2859-6-1.
Rhoades KR: The microscopic lesions of acute fowl cholera in mature chickens. Avian Dis. 1964, 8: 658-665. 10.2307/1587954.
Heddleston KL: Studies on pasteurellosis. V. Two immunogenic types of Pasteurella multocida associated with fowl cholera. Avian Dis. 1962, 6: 315-321. 10.2307/1587901.
Boyce JD, Seemann T, Adler B, Harper M: Pathogenomics of Pasteurella multocida. Curr Top Microbiol Immunol. 2012, 361: 23-38. 10.1007/82_2012_203.
Boyce JD, Harper M, Wilkie IW, Adler B: Pasteurella. Pathogenesis of Bacterial Infections in Animals. Chapter 17. Edited by: Gyles CL, Prescott JF, Songer JG, Thoen CO. 2010, Ames, Iowa: Wiley-Blackwell, 4,
Lee MD, Wooley RE, Brown J, Glisson JR: A survey of potential virulence markers from avian strains of Pasteurella multocida. Vet Microbiol. 1991, 26: 213-225. 10.1016/0378-1135(91)90015-8.
Tatum FM, Yersin AG, Briggs RE: Construction and virulence of a Pasteurella multocida fhaB2 mutant in turkeys. Microb Pathog. 2005, 39: 9-17. 10.1016/j.micpath.2005.05.003.
Tatum FM, Tabatabai LB, Briggs RE: Cross-protection against fowl cholera disease with the use of recombinant Pasteurella multocida FHAB2 peptides. Avian Dis. 2012, 56: 589-591. 10.1637/9991-111611-ResNote.1.
Kumar S, Blaxter ML: Comparing de novo assemblers for 454 transcriptome data. BMC Genomics. 2010, 11: 571-10.1186/1471-2164-11-571.
Besemer J, Lomsadze A, Borodovsky M: GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001, 29: 2607-2618. 10.1093/nar/29.12.2607.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
McClure MA, Smith C, Elton P: Parameterization studies for the SAM and HMMER methods of hidden Markov model generation. Proc Int Conf Intell Syst Mol Biol. 1996, 4: 155-64.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.
Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW: RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35: 3100-3108. 10.1093/nar/gkm160.
Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14: 1394-1403. 10.1101/gr.2289704.
Darzentas N: Circoletto: visualizing sequence similarity with Circos. Bioinformatics. 2010, 26: 2620-2621. 10.1093/bioinformatics/btq484.
Cingolani P, Platts A, le Wang L, Coon M, Nguyen T, Wang L, Land SJ, Lu X, Ruden DM: A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012, 6 (2): 80-92. 10.4161/fly.19695.
Muraoka WT, Zhang Q: Phenotypic and genotypic evidence for L-fucose utilization by Campylobacter jejuni. J Bacteriol. 2011, 193: 1065-1075. 10.1128/JB.01252-10.
Stahl M, Friis LM, Nothaft H, Liu X, Li J, Szymanski CM, Stintzi A: L-fucose utilization provides Campylobacter jejuni with a competitive advantage. Proc Natl Acad Sci USA. 2011, 108: 7194-7199. 10.1073/pnas.1014125108.
Ahir VB, Roy A, Jhala MK, Bhanderi BB, Mathakiya RA, Bhatt VD, Padiya KB, Jakhesara SJ, Koringa PG, Joshi CG: Genome sequence of Pasteurella multocida subsp. gallicida Anand1_poultry. J Bacteriol. 2011, 193: 5604-10.1128/JB.05706-11.
Michael GB, Kadlec K, Sweeney MT, Brzuszkiewicz E, Liesegang H, Daniel R, Murray RW, Watts JL, Schwarz S: ICEPmu1, an integrative conjugative element (ICE) of Pasteurella multocida: structure and transfer. J Antimicrob Chemother. 2012, 67: 91-100. 10.1093/jac/dkr411.
Liu W, Yang M, Xu Z, Zheng H, Liang W, Zhou R, Wu B, Chen H: Complete genome sequence of Pasteurella multocida HN06, a toxigenic strain of serogroup D. J Bacteriol. 2012, 194: 3292-3293. 10.1128/JB.00215-12.
Muhairwa AP, Christensen JP, Bisgaard M: Investigations on the carrier rate of Pasteurella multocida in healthy commercial poultry flocks and flocks affected by fowl cholera. Avian Pathol. 2000, 29: 133-142. 10.1080/03079450094162.
Christensen H, Bisgaard M, Bojesen AM, Mutters R, Olsen JE: Genetic relationships among avian isolates classified as Pasteurella haemolytica, Actinobacillus salpingitidis' or Pasteurella anatis with proposal of Gallibacterium anatis gen. nov., comb. nov. and description of additional genomospecies within Gallibacterium gen. nov. Int J Syst Evol Microbiol. 2003, 53 (Pt 1): 275-87.
Hatfaludi T, Al-Hasani K, Boyce JD, Adler B: Outer membrane proteins of Pasteurella multocida. Vet Microbiol. 2010, 14: 1-17.
Bosch M, Garrido ME, Llagostera M, Perez De Rozas AM, Badiola I, Barbe J: Characterization of the Pasteurella multocida hgbA gene encoding a hemoglobin-binding protein. Infect Immun. 2002, 70: 5955-64. 10.1128/IAI.70.11.5955-5964.2002.
Cox AJ, Hunt ML, Boyce JD, Adler B: Functional characterization of HgbB, a new hemoglobin binding protein of Pasteurella multocida. Microb Pathog. 2003, 34: 287-96. 10.1016/S0882-4010(03)00062-7.
Garcia N, Fernandez-Garayzabal JF, Goyache J, Dominguez L, Vela AI: Associations between biovar and virulence factor genes in Pasteurella multocida isolates from pigs in Spain. Vet Rec. 2011, 169: 362-10.1136/vr.d4869.
Rocha EP, Smith JM, Hurst LD, Holden MT, Cooper JE, Smith NH, Feil EJ: Comparisons of dN/dS are time dependent for closely related bacterial genomes. J Theor Biol. 2006, 239: 226-35. 10.1016/j.jtbi.2005.08.037.
Gardy JL, Spencer C, Wang K, Ester M, Tusnady GE, Simon I, Hua S, DeFays K, Lambert C, Nakai K: PSORT-B: Improving protein subcellular localization prediction for Gram-negative bacteria. Nucleic Acids Res. 2003, 31: 3613-7. 10.1093/nar/gkg602.
Harper M, Cox AD, Adler B, Boyce JD: Pasteurella multocida lipopolysaccharide: The long and the short of it. Vet Microbiol. 2011, 153: 109-15. 10.1016/j.vetmic.2011.05.022.
Harper M, St Michael F, John M, Vinogradov E, Adler B, Boyce JD, Cox AD: Pasteurella multocida Heddleston serovars 1 and 14 express different lipopolysaccharide structures but share the the same lipopolysaccharide biosynthesis outer core locus. Vet Microbiol. 2011, 150: 289-96. 10.1016/j.vetmic.2011.01.028.
Harper M, St Michael F, Vinogradov E, John M, Boyce JD, Adler B, Cox AD: Characterization of the lipopolysaccharide from Pasteurella multocida Heddleston serovar 9; identification of a proposed bi-functional dTDP-3-acetamido-3,6-dideoxy-a-D-glucose biosynthesis enzyme. Glycobiology. 2012, 22: 332-44. 10.1093/glycob/cwr147.
St Michael F, Harper M, Parnas H, John M, Stupak J, Vinogradov E, Adler B, Boyce JD, Cox AD: Structural and genetic basis for the serological differentiation of Pasteurella multocida Heddleston serotypes 2 and 5. J Bacteriol. 2009, 191: 6950-59. 10.1128/JB.00787-09.
St Michael F, Li J, Cox AD: Structural analysis of the core oligosaccharide from Pasteurella multocida strain X73. Carbohydr Res. 2005, 340: 1253-57. 10.1016/j.carres.2005.02.014.
Harper M, St Michael F, Vinogradov E, John M, Steen JA, Van Dorsten L, Boyce JD, Adler B, Cox AD: Structure and biosynthetic locus of the lipopolysaccharide outer core produced by Pasteurella multocida serovars 8 and 13 and the identification of a novel phosphoglycero moity. Glycobiology. 2013, 23: 286-294. 10.1093/glycob/cws154.
This work was supported by the Biotechnology Research and Development Corporation, Peoria, Illinois, USA. Tools for comparative genome analysis were provided through support of the Minnesota Supercomputing Institute.
The authors declare that they have no competing interests.
TJJ performed the genomic analysis, and was the primary author of this study. JEA participated in bioinformatics analyses, including sequence annotation, alignments and pathway reconstruction. SSH formatted and prepared assemblies and annotations for submission to GenBank. MH was involved in analyzing the genome sequences. FMT participated in the editorial review of the manuscript. SKM coordinated this study and helped to draft the manuscript. REB conceived this study, performed the genome sequences data and participated in writing of the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Coding regions present in Pasteurella multocida strain P1059 but absent from strains Pm70 and X73, excluding prophage-associated regions. (PDF 98 KB)
Additional file 2: Table S2: Coding regions present in Pasteurella multocida strain X73 but absent from strains Pm70 and P1059, excluding prophage-associated regions. (PDF 106 KB)
About this article
Cite this article
Johnson, T.J., Abrahante, J.E., Hunter, S.S. et al. Comparative genome analysis of an avirulent and two virulent strains of avian Pasteurella multocida reveals candidate genes involved in fitness and pathogenicity. BMC Microbiol 13, 106 (2013). https://doi.org/10.1186/1471-2180-13-106