Genetic and Phylogenetic Characteristics of Pasteurella multocida Isolates From Different Host Species

Pasteurella multocida is a leading cause of respiratory diseases in many host species. To understand the genetic characteristics of P. multocida strains isolated from different host species, we sequenced the genomic DNA of P. multocida isolated from pigs and analyzed the genetic characteristics of strains from avian species, bovine species, pigs, and rabbits using whole genome sequence (WGS) data. Our results found that a capsular: lipopolysaccharide (LPS): multilocus sequence typing (MLST) genotype A: L1: ST129 (43.75%) was predominant in avian P. multocida; while genotypes B: L2: ST122 (60.00%) and A: L3: ST79 (30.00%) were predominate in bovine P. multocida; genotype D: L6: ST50 (37.50%) in porcine P. multocida; and genotype A: L3: ST9 (76.47%) in rabbit P. multocida. Comparative genomic analysis of P. multocida from different host species found that there are no genes in the P. multocida genome that are specific to any type of host. Phylogenetic analysis using either whole-genome single nucleotide polymorphisms (SNPs) or the set of SNPs present in all single-copy core genes across genomes showed that P. multocida strains with the same LPS genotype and MLST genotype were clustered together, suggesting the combining both the LPS and MLST typing schemes better explained the topology seen in the P. multocida phylogeny.

Introduction

Pasteurella multocida infects a wide spectrum of domestic and wild animals such as poultry and wild birds, pigs, cattle and buffaloes, rabbits, small ruminants, cats (including house cats and large wild cats), dogs and other mammals (Wilkie et al., 2012; Wilson and Ho, 2013). P. multocida isolates are associated with clinical manifestations ranging from asymptomatic or mild chronic upper respiratory inflammation to acute, pneumonic and/or disseminated disease (Wilson and Ho, 2013). While less common, P. multocida can also cause infections in humans through animal bites and/or scratches (Wilson and Ho, 2013).

There are many virulence factors contributing to the pathogenesis of P. multocida. Currently known virulence factors include genes involved in the formation of the capsule, lipopolysaccharide (LPS), fimbriae and adhesins, toxins, iron regulated and iron acquisition proteins, sialic acid metabolism, hyaluronidase, and outer membrane proteins (OMPs) (Harper et al., 2006). The two key surface components, capsule and LPS, form the main typing basis of P. multocida. Serologically, the bacterium is classified into five capsular serogroups (A, B, D, E, F) and/or 16 somatic serotypes according to its capsule and/or LPS antigens, respectively (Carter, 1955; Heddleston et al., 1972). These traditional serological typing methods are too complicated to be used, because the preparation of the high-sensitive anti-serum required for these methods is very difficult (Peng et al., 2016). Therefore, molecular typing methods have been developed to help assign P. multocida into five capsular genotypes (A, B, D, E, F) (Townsend et al., 2001) and eight LPS genotypes (L1-L8) (Harper et al., 2015). In addition, the multilocus sequence typing method has been also designed for P. multocida (Subaaharan et al., 2010), and it has been widely used in epidemiology and surveillance (Hotchkiss et al., 2011; Moustafa et al., 2013; Fernández-Rojas et al., 2014; Peng et al., 2018).

It is documented that the infection of P. multocida isolates displays host predilection (Wilkie et al., 2012; Wilson and Ho, 2013). For example, a P. multocida B: 2 strain can kill cattle and buffaloes at a low dose, but it has no effect on chickens, even at very high doses (Aktories et al., 2012). Isolates from non-avian hosts generally do not cause symptoms of fowl cholera in birds (Ahir et al., 2011; Peng et al., 2017). However, little is known about the genetic characteristics of P. multocida strains circulating in different hosts. In this study, the genomic DNA of P. multocida strains isolated from pigs was sequenced, and the genetic characteristics of P. multocida isolates from avian hosts, bovine hosts, pigs and rabbits were determined using the previously published whole genome sequences (WGSs) data. The aim of this study is to understand the genetic characteristics of P. multocida strains circulating in different host species.

Materials and Methods

P. multocida Strains and Whole Genome Sequencing

A total of 47 P. multocida strains were selected for whole genome sequencing in this study (Table 1). Most of these strains were isolated from pigs with respiratory disorders in China (Peng et al., 2018), with the exception of strain HN04, which is a capsular type B isolate from a swine haemorrhagic septicaemia case. The capsular types and LPS genotypes of these 47 P. multocida isolates were determined through PCR assays, as described previously (Townsend et al., 2001; Harper et al., 2015).

TABLE 1

Table 1. Pasteurella multocida strains discussed in this study.

The genomic DNA (gDNA) of the 47 P. multocida strains were isolated using a QIAGEN Blood and Cell Culture DNA Midi Kit (Lot NO. 157037192, QIAGEN, Germany) with QIAGEN Genomic-tip 100/G (Lot NO. 157024372, QIAGEN, Germany). DNA quantity and quality were evaluated by electrophoresis on a 1% agarose gel and evaluation using a NanoDrop2000 (Thermo Scientific, Waltham, USA), respectively.

DNA libraries were constructed using NEBNext^®Ultra™ II DNA Library Prep Kit. Sequencing was carried out on an Illumina Hiseq Xten platform (Illumina Inc., San Diego, USA) at Guangdong Magigene Biotechnology Co. LTD (Guangzhou, China), using the pair-end 150-bp sequencing protocol. The strategy yielded 5,308,336~12,159,104 raw reads. The raw reads were then filtered to eliminate reads with low quality according to the following criteria: low quality base pairs at each terminal of the reads (Quality-Value < 20) were removed; reads with over short length (parameter setting at 50 bp), > 10% Ns (parameter setting at 10 bp) or > 15 bp overlap overlap with Illumina TruSeq adapter sequences (parameter setting at 15 bp) were removed. After that, 4,732,078~10,715,272 clean reads (Q20% = 100, Q30% ≥ 95.64) were produced. High quality reads were de novo assembled via SPAdes v3.9.0 (Bankevich et al., 2012) to generate contigs.

In addition to the 47 porcine P. multocida sequenced herein, we also retrieved genome sequences of P. multocida from different hosts from the National Center for Biotechnology Information (NCBI) genome database (https://www.ncbi.nlm.nih.gov/genome/genomes/912). By the time of this study design (Assessed at December 31, 2017), there are 132 P. multocida genome sequences publicly available. To ensure the reliability and quality of our analysis, only those high-quality genomes assemblies with contigs numbers ≤ 50 plus coverage ≥50 × were selected. While not meeting the criterion of selection, the two avian P. multocida isolates P1059 (GenBank accession number: CM001581) and X73 (GenBank accession number: CM001580) were still included, because their phenotypical and genetic characteristics have been published (Johnson et al., 2013). The strategy resulted in 62 P. multocida genome sequences, which include 16 avian, 20 bovine, 9 porcine, and 17 rabbit isolates (Table 1). A total of 109 P. multocida genome sequences were finally used in this study, including 47 sequenced herein plus 62 retrieved from NCBI database (Table 1).

Draft genome sequences of the 47 porcine P. multocida were annotated using prokka v1.12 (Seemann, 2014) following the NCBI Prokaryotic Genome Annotation Pipeline (Angiuoli et al., 2008). The 47 annotated genomes were submitted to GenBank accession numbers are displayed in Table 1, which are highlighted in blue text.

Genotyping

The capsular genotypes and LPS genotypes of the 47 P. multocida strains sequenced were determined using PCR methods (Townsend et al., 2001; Harper et al., 2015) and are documented on our previous publication (Peng et al., 2018). Capsular genotypes and LPS genotypes of the 62 P. multocida genomes from NCBI genome database were determined by blasting the oligonucleotide sequences designed for capsular typing (Townsend et al., 2001) and LPS genotyping (Harper et al., 2015) against their WGSs, respectively. MLST genotypes of the 109 P. multocida strains were assigned by performing blast of their WGSs against the P. multocida MLST Databases (https://pubmlst.org/pmultocida). Currently, there are two separate P. multocida MLST schemes on the database, i.e., Multi-host MLST and RIRDC MLST. As for some bovine isolates (BUKK, Islm, Pesh, PVAcc, THA, THD, THF, TX1, V1) included in this study, an existing study has reported the use of RIRDC MLST in determining their MLST genotypes (Moustafa et al., 2013), we use the same MLST scheme to determine the P. multocida MLST genotypes in this study.

Phylogenetic Trees

The phylogenetic relationship between P. multocida strains from different host species was predicted by analyzing the whole-genome single nucleotide polymorphisms (WG-SNPs) as well as the set of SNPs present in all single-copy core genes across genomes (CG-SNPs). The WG-SNPs were identified by comparing each of the WGSs against the reference P. multocida ATCC 43137 genome sequence (GenBank accession NO. CP008918) using MAFFT (version 7.222) software (Katoh and Standley, 2013). A phylogenetic tree based on these WG-SNPs was exported using Parsnp (version 1.2) software (Treangen et al., 2014). Another phylogenetic tree based on CG-SNPs was also constructed. In this strategy, gene contents of each of the 109 P. multocida genomes were clustered using Roary v3.11.0 software (Page et al., 2015) with a criterion of 90% DNA identity plus 90% gene coverage as the minimum criteria for a match to define core genes at first, then SNPs within all the single-copy core genes were compared to generate a neighbor-joining tree using MEGA6 with Kimura 2-parameter model and 1000 bootstrap (Tamura et al., 2013).

Core and Pan-Genome Analysis

Gene contents of each of the 109 P. multocida genomes were clustered using Roary v3.11.0 (Page et al., 2015). Homologous genes were identified by using 90% DNA identity and 90% gene coverage as the minimum criteria for a match. Because potential genes might be disrupted by the gaps between the contigs, we defined those genes present in 99~100% of the strains as core genes. Dispensable genes were defined as those genes shared by at least two strains but shared by <99% strains, while unique genes were defined those present only in a single strain. Host-specific genes were defined as those genes shared by all strains isolated from the same host species but absent in the strains from the other host species. Core genes, dispensable genes, and strain-specific were annotated using the Clusters of Orthologous Groups of proteins (COGs) database (http://www.ncbi.nlm.nih.gov/COG/) for functional analysis (Galperin et al., 2015).

Distribution of Main Virulence Factors Associated Genes

Genes identified as virulence factors (VFGs) were predicted by performing BLAST analysis of the genome sequence against the bacterial Virulence Factor Database (VFDB) (Chen et al., 2005, 2012). After that, a BLAST score ratio (BSR) analysis was performed on these predicted VFGs, and the presence of a VFG in a genome was determined based upon the BSR analysis with a normalized ratio ≥0.8, as described previously (Rasko et al., 2008). A heat-map showing the BSR ratio values of each gene across the genomes of all bacterial strains was generated.

Results

General Features of the 47 Porcine P. multocida Genomes

Whole genome sequencing yielded approximately 796.25~1823.87 Mbp raw data for the 47 porcine P. multocida isolates. The data filter strategy produced approximately 701.86~1589.23 Mb clean data for assembly. Sequences assembled using SPAdes v3.9.0 (Bankevich et al., 2012) generated approximately 18~66 contigs for the 47 porcine P. multocida isolates, with an average of 25 contigs for each strain. These contigs represented a 161~631-kb N50 fragment with an average contig size of approximately 35~130 kb. The predicted genome sizes ranged from 2.22 to 2.49 Mbp, with average G+C content ranging from 40.04 to 40.67% for the 47 P. multocida isolates. The average genome size and G+C content for each of the P. multocida sequenced were 2.35 Mbp and 40.3%, respectively. The number of coding DNA sequences among the genome sequences ranged from 2,027 to 2,330. The general features of the 47 sequenced P. multocida isolates from pigs and their accession numbers as well as those public available genomes of P. multocida from different hosts are shown in Table 1.

Genotypical Characteristics of the P. multocida Isolates From Different Hosts

Capsular Genotypes

Capsular typing strategy identified two categories of capsular genotypes (A and F) for the 16 avian P. multocida isolates, two categories of capsular genotypes (A and B) for the 20 bovine P. multocida isolates, four categories of capsular genotypes (A, B, D, and F) for the 56 porcine P. multocida isolates, and two categories of capsular genotypes (A and D) for the 17 rabbit P. multocida isolates. Type A was the most common capsular genotype for the avian P. multocida isolates (93.75% of the isolates) and rabbit P. multocida isolates (94.12% of the isolates); types A and B were the common capsular genotypes for the bovine P. multocida isolates (40.00% of the isolates and 60.00% of the isolates, respectively); types A and D were the common capsular genotypes for the porcine P. multocida isolates (48.21% of the isolates and 42.86% of the isolates, respectively; see Table 1 and Image 1 in Supplementary Materials).

LPS Genotypes

LPS typing strategy identified five categories of LPS genotypes (L1, L2, L3, L5, and L6) for the 16 avian P. multocida isolates. L1 (56.25%) was the dominant genotype with the highest percentage, followed by L3 (25.00%). For the 20 bovine P. multocida isolates, three categories of LPS genotypes (L1, L2, L3) were identified, and L2 (60.00%) was the mostly common LPS genotypes, followed by L3 (35.00%). There were also three categories of LPS genotypes (L2, L3, L6) defined for the 56 porcine P. multocida isolates. Among them, L6 (66.07%) ranked the highest position, followed by L3 (32.14%). For the 17 rabbit P. multocida isolates, two categories of LPS genotypes (L3 and L6) were identified, and L3 (76.47%) was the most cited (see Table 1 and Image 2 in Supplementary Materials).

When combining capsular types with LPS genotypes, the mostly commonly identified capsule: LPS genotype for the 16 avian isolates was A: L1 (56.25%), followed by A: L3 (18.75%). For the 20 bovine isolates, type B: L2 (60.00%) was the most common, followed by A: L3 (35.00%). Types D: L6 (42.86%), A: L3 (25.00%), and A: L6 (23.21%) were the dominant capsule: LPS genotypes for the 56 porcine isolates, and A: L3 (76.47%) was the most dominant capsule: LPS genotypes for the 17 rabbit isolates (Figure 1).

FIGURE 1

Figure 1. Distribution of different capsule: LPS genotypes in P. multocida isolated from different hosts.

MLST Genotypes

RIRDC MLST assigned eight categories of MLST genotypes (ST8, ST9, ST27, ST53, ST60, ST129, ST156, ST159) for the 16 avian P. multocida, and ST129 (43.75%) was the predominate MLST genotype. The 20 bovine P. multocida were assigned into four categories of MLST genotypes (ST65, ST79, ST80, ST122). For these genotypes, ST122 (60.00%) was the most dominant. For the 56 porcine P. multocida, a total of eight categories of MLST genotypes (ST7, ST9, ST13, ST27, ST50, ST74, ST122, ST287) were present, and ST50 (37.50%), ST13 (23.21%), and ST74 (21.43%) were the common MLST genotypes. The 17 rabbit P. multocida were assigned into five categories of MLST genotypes (ST9, ST27, ST50, ST204, and ST298), and ST9 (52.94%) was the most prevalent (see Table 1 and Image 3 in Supplementary Materials).

When combining MLST genotypes with capsular types and LPS genotypes, the results showed that a capsular: LPS: MLST genotype A: L1: ST129 (43.75%) was predominant among the 16 avian isolates; and B: L2: ST122 (60.00%) as well as A: L3: ST79 (30.00%) were predominate among the 20 bovine isolates. For the 56 porcine isolates, D: L6: ST50 (37.50%) was the predominate genotype; while for the 17 rabbit isolates, A: L3: ST9 (76.47%) was the predominate (Figure 2).

FIGURE 2

Figure 2. Distribution of different capsule: LPS: MLST genotypes in P. multocida isolated from different hosts.

Phylogeny

We generated two phylogenetic trees among the P. multocida isolates from different isolates collected in this study, and they were based on two criteria, respectively. The first tree was generated using 21,765 SNPs across the WGSs of the P. multocida comparing against the reference P. multocida ATCC 43137 genome. Using this strategy, it is clearly showed that P. multocida strains with the same LPS plus MLST genotypes were phylogenetically clustered together (Figure 3). The second tree was constructed based on the SNPs within all single-copy core genes among the comparison genomes. Based on this tree, P. multocida strains with the same LPS plus MLST genotypes were also phylogenetically clustered together (Figure 4).

FIGURE 3

Figure 3. Phylogenetic tree generated based on the SNPs across the whole genome sequences of the P. multocida comparing against the reference P. multocida ATCC 43137 genome. P. multocida capsular type A strains were marked in red, type B strains were marked in green, type D strains were marked in blue, and type F strains were marked in purple. Avian isolates are marked with “”; Bovine isolates are marked with “”; Porcine isolates are marked with “”; Rabbit isolates are marked with “”.

FIGURE 4

Figure 4. Phylogenetic tree generated based on the SNPs within the single copy core genes among the comparison P. multocida genomes. P. multocida capsular type A strains were marked in red, type B strains were marked in green, type D strains were marked in blue, and type F strains were marked in purple. Avian isolates are marked with “”; Bovine isolates are marked with “”; Porcine isolates are marked with “”; Rabbit isolates are marked with “”.

Core and Pan-Genome Analysis

All genes encoded by the 109 P. multocida genomes were used for the determination of homologous and unique genes. These analyses identified a shared set of 1,806 genes (core genes), 1,841 dispensable genes, and 609 strain-specific genes (Figure 5). According to COG functional prediction, the core genes mainly participated in translation, ribosomal structure and biogenesis (J), amino acid transport and metabolism (E), carbohydrate transport and metabolism (G), cell wall/membrane/envelope biogenesis (M), energy production and conversion (C), coenzyme transport and metabolism (H), inorganic ion transport and metabolism (P), and posttranslational modification, protein turnover, chaperones (O); while the dispensable genes mainly involved in mobilome (prophages, transposons) (X), replication, recombination and repair (L), and transcription (K); the strain-specific genes mainly participated in mobilome (X), replication, recombination and repair (L), defense mechanisms (V), and transcription (K) (Figure 5).

FIGURE 5

Figure 5. Distribution of core and dispensable genes among the 109 P. multocida and their COG functions.

To figure out genes carried by P. multocida genome that are associated with host predilection, we first defined a host predilection related gene as one present among the strains from one host species but absent in strains from all the other host species. This criterion found that there were no genes specific to any hosts. Therefore, we next defined a host predilection related gene was determined as one gene present among 90% of the strains from one host species but absent in 90% of the strains from all the other host species. However, this criterion also found that there were no genes specific to any hosts.

Distribution of Virulence Factor-Associated Genes

A total of 432 putative VFGs were identified by performing BLAST analysis of the 109 P. multocida genome sequences against the VFDB database, and these VFGs were assigned into 366 VFG categories (see Table S1 in Supplementary Materials). According to BSR analysis, the distribution of VFGs between different P. multocida strains was mainly reflected on the genes responsible for the synthesis of capsule, LPS and fimbrial low-molecular-weight proteins; however, none of these genes were specific to any hosts (see Image 4 in Supplementary Materials). Interestingly, some VFGs displayed a certain level of LPS: MLST genotype-preference. For example, the hyaluronidase encoding gene pmHAS was absent in the genotypes L6: ST50, L6: ST287, and L2: ST122 strains; the surface fibrils encoded by the gene hsf-2 was missing in genotypes L6: ST50 and L3: ST13 strains (see Image 5 in Supplementary Materials).

Discussion

Birds, pigs, cattle and buffaloes, as well as rabbits are the natural hosts of P. multocida (Wilkie et al., 2012; Wilson and Ho, 2013). In this study, we genotyped P. multocida strains from these host species using their WGS data. Types A: L1 (56.25%) and A: L3 (18.75%) were the common capsule: LPS genotypes for avian P. multocida. It is widely documented that P. multocida strains causing avian pasteurellosis (fowl cholera) are most frequently designated Carter: Heddleston serotypes A:1, A:3, or A:4 (Adler et al., 1999; Wilson and Ho, 2013), and P. multocida Heddleston serotype 1 and 3/4 strains produced LPS genotypes L1 and L3, respectively (Harper et al., 2011, 2013). Therefore, our results are in agreement with those previously documents (Adler et al., 1999; Wilson and Ho, 2013). It has been noted that P. multocida Carter: Heddleston serotype B: 2 strains are commonly responsible for bovine haemorrhagic septicaemia while serovar A: 3 strains are associated with bovine respiratory diseases (Welsh et al., 2004; Dey et al., 2007; Shivachandra et al., 2011; Lainson et al., 2013; Wilson and Ho, 2013). Both haemorrhagic septicaemia and respiratory diseases including pneumonia are the common infections associated with P. multocida in bovine species (bovine pasteurellosis) (Wilson and Ho, 2013). The P. multocida Heddleston serotype 2 strains produced LPS genotypes L2 (St Michael et al., 2009), and Heddleston serotype 3 strains produced LPS genotype L3 (Harper et al., 2013). This is why genotypes B: L2 (60.00%) and A: L3 (35.00%) possess the first two ranks among the P. multocida strains of bovine source. The predominate capsule: LPS genotypes defined among the porcine isolates were types D: L6 (42.86%), A: L3 (25.00%), and A: L6 (23.21%). Those are because P. multocida strains with porcine pasteurellosis are frequently designated capsular types A and D (Tang et al., 2009; Liu et al., 2017; Peng et al., 2018), and LPS serovars 3, 4, and 12 (Pijoan et al., 1983; Lainson et al., 2002; Jamaludin et al., 2005). In particular, these results are in agreement with the latest molecular epidemical data in China (Peng et al., 2018). Interestingly, the capsular: LPS genotype of an isolate (strain HN04) from a swine haemorrhagic septicaemia case was B: L2. This finding suggests this type of P. multocida might also be able to cause haemorrhagic septicaemia in pigs. The mostly common capsule: LPS genotype identified in rabbit P. multocida isolates was A: L3 (76.47%). This result is similar to the previous publication (Townsend et al., 2001), as genotype A: L3 corresponds to P. multocida Carter: Heddleston serotype A: 3 or A: 4 (Harper et al., 2013).

The predominate MLST genotype of P. multocida from avian species was ST129 (43.75%), and this result is in agreement with Wang et al. (2013). The difference between the two studies is that Wang et al. only defined ST129 from the avian P. multocida strains but we identified multiple MLST genotypes in addition to ST129; other MLST genotypes such as ST8, ST9, ST27, ST53, ST60 are also common in avian P. multocida (Subaaharan et al., 2010). These findings suggest the prevalence of P. multocida genotypes from different geographic regions might be different. ST122 was the common MLST genotype in P. multocida isolates from bovine species (Table 1). This genotype is widely documented to be associated with bovine haemorrhagic septicaemia, a common pasteurellosis in bovine species (Hotchkiss et al., 2011). Interestingly, the MLST genotype of the porcine isolate HN04 from swine haemorrhagic septicaemia case was also ST122, suggesting that ST122 might be particularly associated with the disease other than the host species. The mostly common MLST genotypes of P. multocida from pigs defined in this study were ST50 (37.50%), ST13 (23.21%), and ST74 (19.64%). These results are quite different from the date of a recent molecular epidemiological study (Peng et al., 2018), but this could be explained by the two studies using two different MLST databases to determine the MLST genotypes. When we used the Multiple host MLST database to re-assigned these strains, we found that all ST50, ST13, and ST74 strains were re-assigned as ST11, ST3, and ST10, respectively. These results are in agreement with the latest epidemical data from China (Peng et al., 2018), and also in agreement with the data from Spain (García-Alvarez et al., 2017).

The predominate MLST genotype in P. multocida isolates from rabbit species was ST9 (76.47%), which was re-assigned as ST12 using the Multiple host MLST database. Previous studies have documented that ST12 was the common MLST genotype in rabbit P. multocida (García-Alvarez et al., 2015; Massacci et al., 2018). Interestingly, a previous study documented that the ST11 (assigned using Multiple host MLST) clone of P. multocida is widely spread among farmed rabbits in the Iberian Peninsula and demonstrates respiratory niche association (García-Alvarez et al., 2015). However, the percentage of the ST11 clone (defined as ST50 using RIRDC MLST) defined in this study was very low, only 5.88% (one isolate). According to NCBI genome database (https://www.ncbi.nlm.nih.gov/genome/genomes/912), the WGSs data from the 17 rabbit P. multocida using in this study are from France, these findings also suggest that the prevalence of P. multocida genotypes in different geographic regions might be different.

When combining the capsular genotypes, LPS genotypes and the MLST genotypes, it is interesting to see that although some P. multocida strains from different hosts share the same capsule: LPS genotypes, their MLST genotypes are quite different. The MLST genotypes of P. multocida capsule: LPS genotype A: L1 strains from avian species were ST129 and ST69, but the MLST genotype of those type A: L1 strains from bovine species was ST65 (Table 1). However, P. multocida strains with the same capsule: LPS: MLST genotype are able to spread in different host species but cause similar disease symptoms. For example, the capsule: LPS: MLST genotype B: L2: ST122 strain is associated with both bovine haemorrhagic septicaemia and swine haemorrhagic septicaemia; the F: L3: ST9 clone can lead to both avian respiratory disease and porcine respiratory disease; the D: L6: ST50 clone is associated with both the porcine respiratory disease and the rabbit respiratory disease (Table 1).

Although many publications have documented that the infection of P. multocida isolates displays host predilection (Kubatzky, 2012; Wilkie et al., 2012; Wilson and Ho, 2013), these analyses using comparative genomics found that there were no genes or VFGs specific to any isolates for any specific type of host (see Image 4 in Supplementary Materials). Though the dermonecrotic toxin encoding gene toxA was only present in one isolate from pigs (see Image 5 in Supplementary Materials), this gene has been also identified in P. multocida isolates from bovine species, avian species, small ruminants, sheep, and goats (Ewers et al., 2006; Sarangi et al., 2016; Shirzad Aski and Tabatabaei, 2016). These findings suggest that toxA is not associated with host predilection. Interestingly, it is documented that the presence of toxA displays a significant association to the disease status in swine, even though it is detected in P. multocida isolates from small ruminants, cattle, and poultry as well (Ewers et al., 2006). This suggests that the disease symptoms caused by P. multocida with specific VFGs might be also associated with the host itself and/or the locations that the bacteria are found on the host. However, more studies are necessary for further confirmation. Similar to toxA, the transferrin binding protein A encoding gene tbpA was also present in isolates from bovine species (see Image 5 in Supplementary Materials). However, this gene has been also reported to be identified in P. multocida isolates from ovine species and rabbits (Ewers et al., 2006; Ferreira et al., 2012). These findings suggest that tbpA is also not host specific.

It has been proposed that there is little or no correlation between the phylogenetic groups and phenotypic characteristics (Boyce et al., 2012; Moustafa et al., 2015; Cao et al., 2017). However, in this study, using both the SNPs across the WGSs and the SNPs within all single-copy core genes among the comparison genomes, we found that P. multocida isolates with the same LPS genotype and MLST genotype could be phylogenetically grouped, though these P. multocida strains were isolated from different geographic regions or host species, and/or associated with different disease symptoms (Figures 3, 4). While more data are needed for the further confirmation, these findings suggest the combining both the LPS and MLST typing schemes might better explain the topology seen in the P. multocida phylogeny.

In conclusion, we performed whole genome sequencing on P. multocida strains isolated from pigs and determined the genetic characteristics of P. multocida isolates from avian species, bovine species, pigs and rabbits in this study. Genotyping using WGSs data suggest that P. multocida isolates from different host species displayed a certain preference for “capsule/LPS/MLST genotypes.” Although host predilection in P. multocida infections are widely documented, our analyses found there were no genes specific to any hosts. In addition, our phylogenetic analyses found that P. multocida isolates with the same LPS genotype and MLST genotype could be phylogenetically grouped, which suggests that the combining both the LPS and MLST typing schemes might better explain the topology seen in the P. multocida phylogeny.

Author Contributions

ZP, WL, RZ, HC, and BW participated in the conception and design of the work; ZP, WL, ZfX, ZhX, and ZL performed whole genome sequencing and data analysis; FW and LH contributed to bacterial genomic DNA isolation and purification; ZP wrote the manuscript; RZ, HC, and BW helped to revised the manuscript. All authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We sincerely acknowledge Dr. Yi-Wei Tang at Memorial Sloan Kettering Cancer Center (New York, United States) for his comments on this manuscript. We also thank Ms. Xiaoyun Yao and Mr. Taoyu Zhou at Guangdong Magigene Biotechnology Co. Ltd. (Guangzhou, China) for their help to the whole genome sequencing project. This study was supported in part by the National Key Technology Research and Development Program of the Ministry of Science and Technology of China (grants 2015BAD12B04 and 2014BAD20B01).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01408/full#supplementary-material

Table S1. Virulence factor-associated genes of the 109 P. multocida isolates predicted by the VFDB database.

Image 1. Distribution of different capsule genotypes in P. multocida isolated from different host species.

Image 2. Distribution of different LPS genotypes in P. multocida isolated from different species.

Image 3. Distribution of different MLST genotypes in P. multocida isolated from different species.

Image 4. Heat-map showing the distribution of virulence factor-associated genes among P. multocida isolates from different host species.

Image 5. Heat-map showing the distribution of the 27 main kinds of virulence factor-associated genes among P. multocida isolates from different host species.

References

Adler, B., Bulach, D., Chung, J., Doughty, S., Hunt, M., Rajakumar, K., et al. (1999). Candidate vaccine antigens and genes in Pasteurella multocida. J. Biotechnol. 73, 83–90. doi: 10.1016/S0168-1656(99)00111-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahir, V. B., Roy, A., Jhala, M. K., Bhanderi, B. B., Mathakiya, R. A., Bhatt, V. D., et al. (2011). Genome sequence of Pasteurella multocida subsp. gallicida Anand1_poultry. J. Bacteriol. 193:5604. doi: 10.1128/JB.05706-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Aktories, K., Orth, J. H., and Adler, B. (2012). Pasteurella Multocida: Molecular Biology, Toxins and Infection. Heidelberg: Springer-Verlag.

Google Scholar

Angiuoli, S. V., Gussman, A., Klimke, W., Cochrane, G., Field, D., Garrity, G., et al. (2008). Toward an online repository of Standard Operating Procedures (SOPs) for (meta)genomic annotation. OMICS 12, 137–141. doi: 10.1089/omi.2008.0017

PubMed Abstract | CrossRef Full Text | Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyce, J. D., Seemann, T., Adler, B., and Harper, M. (2012). Pathogenomics of Pasteurella multocida. Curr. Top. Microbiol. Immunol. 361, 23–38. doi: 10.1007/82_2012_203

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, P., Guo, D., Liu, J., Jiang, Q., Xu, Z., and Qu, L. (2017). Genome-wide analyses reveal genes subject to positive selection in Pasteurella multocida. Front. Microbiol. 8:961. doi: 10.3389/fmicb.2017.00961

PubMed Abstract | CrossRef Full Text | Google Scholar

Carter, G. R. (1955). Studies on Pasteurella multocida. I. A hemagglutination test for the identification of serological types. Am. J. Vet. Res. 16, 481–484.

PubMed Abstract | Google Scholar

Chen, L., Xiong, Z., Sun, L., Yang, J., and Jin, Q. (2012). VFDB 2012 update: toward the genetic diversity and molecular evolution of bacterial virulence factors. Nucleic Acids Res. 40, D641–D645. doi: 10.1093/nar/gkr989

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Yang, J., Yu, J., Yao, Z., Sun, L., Shen, Y., et al. (2005). VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 33, D325–D328. doi: 10.1093/nar/gki008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dey, S., Singh, V. P., Kumar, A. A., Sharma, B., Srivastava, S. K., and Singh, N. (2007). Comparative sequence analysis of 16S rRNA gene of Pasteurella multocida serogroup B isolates from different animal species. Res. Vet. Sci. 83, 1–4. doi: 10.1016/j.rvsc.2006.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ewers, C., Lübke-Becker, A., Bethe, A., Kiebling, S., Filter, M., and Wieler, L. H. (2006). Virulence genotype of Pasteurella multocida strains isolated from different hosts with various disease status. Vet. Microbiol. 114, 304–317. doi: 10.1016/j.vetmic.2005.12.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernández-Rojas, M. A., Vaca, S., Reyes-Lopez, M., de la Garza, M., Aguilar-Romero, F., Zenteno, E., et al. (2014). Outer membrane vesicles of Pasteurella multocida contain virulence factors. Microbiology 3, 711–717. doi: 10.1002/mbo3.201

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferreira, T. S., Felizardo, M. R., Sena de Gobbi, D. D., Gomes, C. R., Nogueira Filsner, P. H., Moreno, M., et al. (2012). Virulence genes and antimicrobial resistance profiles of Pasteurella multocida strains isolated from rabbits in Brazil. ScientificWorldJournal 2012:685028. doi: 10.1100/2012/685028

PubMed Abstract | CrossRef Full Text | Google Scholar

Galperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2015). Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 43, D261–D269. doi: 10.1093/nar/gku1223

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Alvarez, A., Chaves, F., Fernandez, A., Sanz, C., Borobia, M., and Cid, D. (2015). An ST11 clone of Pasteurella multocida, widely spread among farmed rabbits in the Iberian Peninsula, demonstrates respiratory niche association. Infect. Genet. Evol. 34, 81–87. doi: 10.1016/j.meegid.2015.07.018

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Alvarez, A., Vela, A. I., San Martin, E., Chaves, F., Fernandez-Garayzabal, J. F., Lucas, D., et al. (2017). Characterization of Pasteurella multocida associated with ovine pneumonia using multi-locus sequence typing (MLST) and virulence-associated gene profile analysis and comparison with porcine isolates. Vet. Microbiol. 204, 180–187. doi: 10.1016/j.vetmic.2017.04.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Harper, M., Boyce, J. D., and Adler, B. (2006). Pasteurella multocida pathogenesis: 125 years after Pasteur. FEMS Microbiol. Lett. 265, 1–10. doi: 10.1111/j.1574-6968.2006.00442.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Harper, M., John, M., Turni, C., Edmunds, M., St Michael, F., Adler, B., et al. (2015). Development of a rapid multiplex PCR assay to genotype Pasteurella multocida strains by use of the lipopolysaccharide outer core biosynthesis locus. J. Clin. Microbiol. 53, 477–485. doi: 10.1128/JCM.02824-14

PubMed Abstract | CrossRef Full Text

Harper, M., St Michael, F., John, M., Vinogradov, E., Adler, B., Boyce, J. D., et al. (2011). Pasteurella multocida Heddleston serovars 1 and 14 express different lipopolysaccharide structures but share the same lipopolysaccharide biosynthesis outer core locus. Vet. Microbiol. 150, 289–296. doi: 10.1016/j.vetmic.2011.01.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Harper, M., St Michael, F., John, M., Vinogradov, E., Steen, J. A., van Dorsten, L., et al. (2013). Pasteurella multocida Heddleston serovar 3 and 4 strains share a common lipopolysaccharide biosynthesis locus but display both inter- and intrastrain lipopolysaccharide heterogeneity. J. Bacteriol. 195, 4854–4864. doi: 10.1128/JB.00779-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Heddleston, K. L., Gallagher, J. E., and Rebers, P. A. (1972). Fowl cholera: gel diffusion precipitin test for serotyping Pasteruella multocida from avian species. Avian Dis. 16, 925–936. doi: 10.2307/1588773

PubMed Abstract | CrossRef Full Text | Google Scholar

Hotchkiss, E. J., Hodgson, J. C., Lainson, F. A., and Zadoks, R. N. (2011). Multilocus sequence typing of a global collection of Pasteurella multocida isolates from cattle and other host species demonstrates niche association. BMC Microbiol. 11:115. doi: 10.1186/1471-2180-11-115

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamaludin, R., Blackall, P. J., Hansen, M. F., Humphrey, S., and Styles, M. (2005). Phenotypic and genotypic characterisation of Pasteurella multocida isolated from pigs at slaughter in New Zealand. N. Z. Vet. J. 53, 203–207. doi: 10.1080/00480169.2005.36505

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, T. J., Abrahante, J. E., Hunter, S. S., Hauglund, M., Tatum, F. M., Maheswaran, S. K., et al. (2013). 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. doi: 10.1186/1471-2180-13-106

PubMed Abstract | CrossRef Full Text | Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubatzky, K. F. (2012). Pasteurella multocida and immune cells. Curr. Top. Microbiol. Immunol. 361, 53–72. doi: 10.1007/82_2012_204

PubMed Abstract | CrossRef Full Text | Google Scholar

Lainson, F. A., Aitchison, K. D., Donachie, W., and Thomson, J. R. (2002). Typing of Pasteurella multocida isolated from pigs with and without porcine dermatitis and nephropathy syndrome. J. Clin. Microbiol. 40, 588–593. doi: 10.1128/JCM.40.2.588-593.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Lainson, F. A., Dagleish, M. P., Fontaine, M. C., Bayne, C., and Hodgson, J. C. (2013). Draft genome sequence of Pasteurella multocida A:3 strain 671/90. Genome Announc. 1:e00803–13. doi: 10.1128/genomeA.00803-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Zhao, Z., Xi, X., Xue, Q., Long, T., and Xue, Y. (2017). Occurrence of Pasteurella multocida among pigs with respiratory disease in China between 2011 and 2015. Ir. Vet. J. 70:2. doi: 10.1186/s13620-016-0080-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Massacci, F. R., Magistrali, C. F., Cucco, L., Curcio, L., Bano, L., Mangili, P., et al. (2018). Characterization of Pasteurella multocida involved in rabbit infections. Vet. Microbiol. 213, 66–72. doi: 10.1016/j.vetmic.2017.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Moustafa, A. M., Bennett, M. D., Edwards, J., Azim, K., Mesaik, M. A., Choudhary, M. I., et al. (2013). Molecular typing of haemorrhagic septicaemia-associated Pasteurella multocida isolates from Pakistan and Thailand using multilocus sequence typing and pulsed-field gel electrophoresis. Res. Vet. Sci. 95, 986–990. doi: 10.1016/j.rvsc.2013.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Moustafa, A. M., Seemann, T., Gladman, S., Adler, B., Harper, M., Boyce, J. D., et al. (2015). Comparative genomic analysis of asian haemorrhagic septicaemia-associated strains of Pasteurella multocida identifies more than 90 haemorrhagic septicaemia-specific genes. PLoS ONE 10:e0130296. doi: 10.1371/journal.pone.0130296

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Z., Liang, W., Wang, Y., Liu, W., Zhang, H., Yu, T., et al. (2017). Experimental pathogenicity and complete genome characterization of a pig origin Pasteurella multocida serogroup F isolate HN07. Vet. Microbiol. 198, 23–33. doi: 10.1016/j.vetmic.2016.11.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Z., Liang, W., and Wu, B. (2016). Molecular typing methods for Pasteurella multocida-A review. Acta Microbiol. Sin. 56, 1521–1529. doi: 10.13343/j.cnki.wsxb.20160002

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Z., Wang, H., Liang, W., Chen, Y., Tang, X., Chen, H., et al. (2018). A capsule/lipopolysaccharide/MLST genotype D/L6/ST11 of Pasteurella multocida is likely to be strongly associated with swine respiratory disease in China. Arch. Microbiol. 200, 107–118. doi: 10.1007/s00203-017-1421-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pijoan, C., Morrison, R. B., and Hilley, H. D. (1983). Serotyping of Pasteurella multocida isolated from swine lungs collected at slaughter. J. Clin. Microbiol. 17, 1074–1076.

PubMed Abstract | Google Scholar

Rasko, D. A., Rosovitz, M. J., Myers, G. S., Mongodin, E. F., Fricke, W. F., Gajer, P., et al. (2008). The pangenome structure of Escherichia coli: comparative genomic analysis of E. coli commensal and pathogenic isolates. J Bacteriol. 190, 6881–6893. doi: 10.1128/JB.00619-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarangi, L. N., Thomas, P., Gupta, S. K., Kumar, S., Viswas, K. N., and Singh, V. P. (2016). Molecular epidemiology of Pasteurella multocida circulating in India by multilocus sequence typing. Transbound. Emerg. Dis. 63:e286–e292. doi: 10.1111/tbed.12270

PubMed Abstract | CrossRef Full Text | Google Scholar

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | CrossRef Full Text | Google Scholar

Shirzad Aski, H., and Tabatabaei, M. (2016). Occurrence of virulence-associated genes in Pasteurella multocida isolates obtained from different hosts. Microb. Pathog. 96, 52–57. doi: 10.1016/j.micpath.2016.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Shivachandra, S. B., Viswas, K. N., and Kumar, A. A. (2011). A review of hemorrhagic septicemia in cattle and buffalo. Anim. Health Res. Rev. 12, 67–82. doi: 10.1017/S146625231100003X

PubMed Abstract | CrossRef Full Text | Google Scholar

St Michael, F., Harper, M., Parnas, H., John, M., Stupak, J., Vinogradov, E., et al. (2009). Structural and genetic basis for the serological differentiation of Pasteurella multocida Heddleston serotypes 2 and 5. J. Bacteriol. 191, 6950–6959. doi: 10.1128/JB.00787-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Subaaharan, S., Blackall, L. L., and Blackall, P. J. (2010). Development of a multi-locus sequence typing scheme for avian isolates of Pasteurella multocida. Vet. Microbiol. 141, 354–361. doi: 10.1016/j.vetmic.2010.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729. doi: 10.1093/molbev/mst197

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, X., Zhao, Z., Hu, J., Wu, B., Cai, X., He, Q., et al. (2009). Isolation, antimicrobial resistance, and virulence genes of Pasteurella multocida strains from swine in China. J. Clin. Microbiol. 47, 951–958. doi: 10.1128/JCM.02029-08

PubMed Abstract | CrossRef Full Text | Google Scholar

Townsend, K. M., Boyce, J. D., Chung, J. Y., Frost, A. J., and Adler, B. (2001). Genetic organization of Pasteurella multocida cap Loci and development of a multiplex capsular PCR typing system. J. Clin. Microbiol. 39, 924–929. doi: 10.1128/JCM.39.3.924-929.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Treangen, T. J., Ondov, B. D., Koren, S., and Phillippy, A. M. (2014). The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 15:524. doi: 10.1186/s13059-014-0524-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhu, J., Lu, C., Wu, B., Liu, D., Hang, W., et al. (2013). Evidence of circulation of an epidemic strain of Pasteurella multocida in Jiangsu, China by multi-locus sequence typing (MLST). Infect. Genet. Evol. 20, 34–38. doi: 10.1016/j.meegid.2013.07.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Welsh, R. D., Dye, L. B., Payton, M. E., and Confer, A. W. (2004). Isolation and antimicrobial susceptibilities of bacterial pathogens from bovine pneumonia: 1994–2002. J. Vet. Diagn. Invest. 16, 426–431. doi: 10.1177/104063870401600510

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkie, I., Harper, M., Boyce, J., and Adler, B. (2012). Pasteurella multocida: diseases and pathogenesis. Curr. Top. Microbiol. Immunol. 361, 1–22. doi: 10.1007/82_2012_216

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, B. A., and Ho, M. (2013). Pasteurella multocida: from zoonosis to cellular microbiology. Clin. Microbiol. Rev. 26, 631–655. doi: 10.1128/CMR.00024-13

PubMed Abstract | CrossRef Full Text | Google Scholar