key: cord-0072558-zk13m01p
authors: Peng, Zhong; Liu, Junyang; Liang, Wan; Wang, Fei; Wang, Li; Wang, Xueying; Hua, Lin; Chen, Huanchun; Wilson, Brenda A.; Wang, Jia; Wu, Bin
title: Development of an Online Tool for Pasteurella multocida Genotyping and Genotypes of Pasteurella multocida From Different Hosts
date: 2021-12-17
journal: Front Vet Sci
DOI: 10.3389/fvets.2021.771157
sha: 5e9682260b0d21fa6ce33b7303b54f186490dc01
doc_id: 72558
cord_uid: zk13m01p

Pasteurella multocida is a versatile zoonotic pathogen. Multiple systems have been applied to type P. multocida from different diseases in different hosts. Recently, we found that assigning P. multocida strains by combining their capsular, lipopolysaccharide, and MLST genotypes (marked as capsular: lipopolysaccharide: MLST genotype) could help address the biological characteristics of P. multocida circulation in different hosts. However, there is still lack of a rapid and efficient tool to diagnose P. multocida according to this system. Here, we developed an intelligent genotyping platform PmGT for P. multocida strains according to their whole genome sequences using the web 2.0 technologies. By using PmGT, we determined capsular genotypes, LPS genotypes, and MLST genotypes as well as the main virulence factor genes (VFGs) of P. multocida isolates from different host species based on their whole genome sequences published on NCBI. The results revealed a closer association between the genotypes and pasteurellosis rather than between genotypes and host species. With the advent of high-quality, inexpensive DNA sequencing, PmGT represents a more efficient tool for P. multocida diagnosis in both epidemiological studies and clinical settings.

Rapid and accurate diagnosis of sources of infections is critical for both medical and veterinary activities, and it is important for improved understanding of disease mechanisms and measures to control the illness (1) . Microbial typing is an important link for the diagnosis of pathogens associated with diseases. The most widely used typing methods consist of serological typing systems and PCR-based molecular typing methods (2, 3) . The establishment of discriminatory typing systems help in the understanding and control of pathogens, especially those with multiple serovars and/or genotypes from different environmental or host sources. Whole genome sequencing combined with the high-end computational technology is such an emerging approach for microbial diagnosis (4) . Using the whole genome sequencing technologies, it is possible to determine the causative agent of infectious diseases rapidly and accurately, including newly emerged ones (5, 6) . However, interpretation of the sequencing results to formulate a definitive diagnosis still requires technical experts with computational and bioinformatics skills. Therefore, a practical, automated platform that combines whole genome sequencing with computational technologies to provide diagnostic outcomes would be beneficial in advancing the field.

Pasteurella multocida is an important zoonotic pathogen and it can colonize and cause infections in a wide range of domestic and wild animals including food producing animals (e.g., poultry, pigs, beef, sheep) and companion animals (e.g., cats and dogs) as well as in humans (7) (8) (9) . Animal diseases associated with P. multocida such as fowl cholera in poultry and other birds, progressive atrophic rhinitis and pneumonic pasteurellosis in pigs, haemorrhagic septicaemia and respiratory diseases in cattle and buffalos, leporine atrophic rhinitis and pneumonic pasteurellosis, are of great economic significance in agriculture (9) . In humans, opportunistic infections of soft tissue, including wound dermonecrosis, respiratory disease with chronic pulmonary, urinary tract infection and bacteremic meningitis have also been reported (9) . Most of these infections are associated with animal biting, scratching, kissing, and/or licking (10) (11) (12) . In this regard, P. multocida represents a risk to public health. P. multocida strains from different hosts are serologically classified into five serogroups (A, B, D, E, F) (13) (14) (15) and/or 16 serovars (serovars 1 to 16) (16), according to their capsular and lipopolysaccharide (LPS) antigens, respectively. However, these two traditional serological typing methods require high-quantity antisera that are challenging to prepare, particularly for clinical use, such those methods are no longer widely used for large-scale epidemiological studies (7, 17) .

In 2001, a multiplex PCR-based method was established to type the five serogroups into five capsular genotypes (A, B, D, E, F) (18) , and in 2015, another multiplex PCR-based method was also developed to classified the 16 serovars into eight LPS genotypes (L1∼L8) (19) . In 2004 and 2010, two multilocus sequencing typing systems were also developed to genotype P. multocida strains (https://pubmlst.org/pmultocida/) from multiple mammalian hosts and birds, respectively (20, 21) . In 2017, a virulence genotyping system based on the detection of different virulence factor gene (VFG) profiles was also reported for distinguishing P. multocida strains from different hosts (22) . Compared to the traditional serological typing methods, these molecular DNA-based typing systems are indeed highly effective and accurate, and they are now widely used to determine the epidemiological and genetic characteristics of clinical isolates (23) (24) (25) (26) (27) .

Despite of more than 135 years of research, differences on the molecular biological characteristics of P. multocida prevalence in different host species remain to be addressed. For example, P. multocida type A strains have been recovered from avian species, pigs, bovine species, and many other host species (8, 9) , but little is known about differences on those type A isolates from different hosts. Recently, we developed a system to assign P. multocida strains from different host species by combining their capsular, LPS, and MLST genotypes (marked as capsular genotype: LPS genotype: MLST genotype), as well as determine the VFG profiles, which contributes to address the molecular biological characteristics of P. multocida prevalence in different host species (7, 23, 27) . However, this strategy requires bioinformatics experts for data analysis and interpretation. Here, we report the development of an automated platform to type P. multocida strains from multiple hosts that combines the use of whole genome sequencing.

P. multocida strains used in this study include one isolate of bovine origin (strain HB01), one isolate of avian origin (strain HB02), and 50 isolates of porcine origin (strains HB03, HN04, HN05, HN06, HN07, HNA01∼HNA22, HND01∼HND21, HNF01 and HNF02) (Supplementary Table S1 ). All of these strains are from our laboratory collection, for which we have previously sequenced their whole genome sequences (27) (28) (29) (30) .

Nucleotide sequences specific for the determination of P. The nucleotide sequences of 23 types of virulence genes commonly detected in P. multocida epidemiological studies, including those encoding fimbriae and other adhesins (ptfA, fimA, hsf-1, hsf-2, pfhA, and tadD), toxin (toxA), iron acquisition proteins (exbB, exbD, tonB, hgbA, hgbB, fur, and tbpA), sialidases (nanB and nanH), hyaluronidase (pmHAS), outer membrane proteins (OMPs) (ompA, ompH, oma87, and plpB), and superoxide dismutase (sodA and sodC), were amplified from the genomic DNA of P. multocida HN06 and HB01 by PCR assays using the protocols documented elsewhere (23, 31) . These nucleotide sequences were deposited in GenBank under accession numbers MT570167∼ MT570189 (Supplementary Text 1) .

The publicly available whole genome sequences of 262 P. multocida strains from bovine species [n = 106; including those recovered bovine haemorrhagic septicaemia cases (32)], avian species (n = 39), porcine species (n = 66), leporine species (n = 20), ovine species (n = 6), humans (n = 13), canines (n = 3), murine species (n = 2), horses (n = 2), cats (n = 2), alpacas (n = 2) and 1 synthetic DNA sequence in NCBI genome database were downloaded for use (Supplementary Table S1 ).

The PmGT platform was integrated on a CentOS server, mainly providing two kinds of online services: genotyping tool, and data query and display. To establish the genotyping online service, we first used Apache (https://www.apache.org) as the web container. Then, we downloaded the BLAST package (ftp://ftp. ncbi.nlm.nih.gov/blast/executables/LATEST/) from NCBI, which was thereafter installed and configured on the web container. PHP was used as the server-side language and the browserside script used jQuery, which is a fast, small, and featurerich JavaScript library. The view pages were constructed with Hypertext Markup Language (HTML) and Cascading Style Sheets (CSS). For the target strain, the format of the sequence was first verified by the web user interface and then the sequence data was uploaded to the server through the PHP program which subsequently called the localized BLAST to align the uploaded sequence with the reference database. The nucleotide sequences specific for the determination of P. multocida strains, capsular genotypes, LPS genotypes, and the 23 types of virulence factor genes (VFGs) were packaged and used as the reference database for sequence alignment. Finally, the result was returned and displayed in the web page. In addition, if the user selected the option of "MLST 

genotyping, " the http request function "curl_setopt" in PHP was used to request PubMLST's RESTful interface (http://rest. pubmlst.org/db/pubmlst_Pmultocida_seqdef/sequence) and the function "curl_exec" was used to catch the response which thereafter was parsed to the result and displayed in the genotyping page.

PCR Detection of Capsular Genotypes, LPS Genotypes, MLST Genotypes, and Virulence Genes of P. multocida Strains From Pigs Capsular genotypes and LPS genotypes of P. multocida strains from our laboratory collection were determined using multiplex PCR-based assays, as documented elsewhere (18, 19) . Profiles of 23 types of virulence genes mentioned above were determined by PCR assays, as described previously (23) . Sequence types (STs) were determined according to the protocols described in Pasteurella multocida MLST database (https://pubmlst.org/ organisms/pasteurella-multocida/multi-host).

Nucleotide sequences specific for P. multocida and its capsular genotypes, LPS genotypes, as well as VFGs were publicly available in GenBank under accession numbers MN938443-MN938455 and MT570167∼MT570189. The typing system developed in the present study is available at: http://vetinfo.hzau.edu.cn/PmGT.

The general process for genotyping is summarized as: when a query sequence is submitted via the web user interface, this sequence will be then submitted to the CentOS server via HTTP protocol. Thereafter, the sequence is evaluated by the PHP program, and the passed sequence will be BLASTed against the genotype database to yield a result, which will be returned to the webpage through the PHP program (Figures 1A,B) .

Through the above procedures, the genotyping module of PmGT (http://vetinfo.hzau.edu.cn/PmGT) was developed (Figure 1) . Currently, PmGT provides the above services includes five menus: (1) the "Home" page gives a brief introduction of P. multocida etiological characteristics to help the users understand the bacterium; (2) the "Isolates" page displays the genotypes of P. multocida strains based on their whole genome sequences that are publicly available in NCBI; this page also provides the link for the users to download the genomes of these P. multocida strains from NCBI; (3) the "Genotyping" page enables the users to determine whether a putative isolate is a P. multocida and genotype P. multocida strains by using the whole genome sequence assembled from the sequencing reads ( Figure 1C) ; (4) the "About" page summarizes the guidelines for the use of this web tool; (5) the "Contact" page provides the contact information of the developers. 

To test the accuracy of PmGT, we used two methods to type 52 P. multocida isolates (HB01, HB02, HB03, HN04, HN05,   HN06 , HN07, HNA01∼HNA22, HND01∼HND21, HNF01, and HNF02) from our laboratory collection (27) . First, we submitted their whole genome sequences to PmGT for genotyping. As a comparison, we also determined the capsular genotypes, LPS genotypes, sequence types, as well as the profile of the abovementioned 23-kinds of virulence genes by using PCR assays. All these 52 strains were genotyped by PmGT and through this online genotyping platform ( Table 1) . Genotyping by PCR assays confirmed these capsular, LPS, and MLST genotypes. PCR results of capsular and LPS genotypes are provided in Supplementary Figures S1, S2 .

Determination of the 23 types of virulence genes for each of the 52 strains by using this online system revealed that several genes (ptfA, fimA, oma87, and sodC) were broadly presented in the genome sequences genotyped (Figure 2) . However, several genes (hsf-1, hsf-2, pfhA, and tadD) were heterogeneously distributed, and in particularly, none of the 52 sequences genotyped carried the toxA or tbpA genes ( Figure 2) . These results were also confirmed by PCR assays (Supplementary Table S1 ).

To understand the genotypes of P. multocida strains circulation in different host species, the 262 whole genome sequences of P. multocida strains were genotyped by PmGT. The results revealed that P. multocida isolates from different hosts displayed a certain preference for "capsular/LPS/MLST genotypes" (Figure 3) . For example, most of the porcine strains were determined as capsular genotypes A (52%) and D (39%), LPS genotypes L3 (36%) and L6 (61%), sequence types ST3 (29%), ST11 (22%), and ST10 (34%), respectively; while most of the genotyped bovine strains were determined as capsular genotypes A (72%) and B (28%), LPS genotypes L3 (67%) and L2 (27%), and sequence types ST1 (59%) and ST44 (25%), respectively (Figure 3) . When combining the capsular genotypes and the LPS genotypes, it revealed that most of the genotyped avian P. multocida were typed as A:L1 and A:L3, while most of the genotyped bovine P. multocida were typed as A:L3 and B:L2; the genotyped porcine P. multocida mainly belonged to D:L6, A:L3, and A:L6; while the genotyped leporine P. multocida mainly belonged to A:L3; most of the genotyped human P. multocida were typed as A:L3 and A:L1 (Figure 4A ). If the capsular genotypes, LPS genotypes, and MLST genotypes were combined, most of the genotyped avian P. multocida were typed as A:L1:ST128, while most of the genotyped bovine P. multocida were typed as A:L3:ST1 and B:L2:ST44; the genotyped porcine P. multocida mainly belonged to D:L6:ST11, A:L3:ST3, and A:L6:ST10; while the genotyped leporine P. multocida mainly belonged to A:L3:ST12 (Figure 4) .

Virulence genotyping using the system developed herein revealed that the presence of multiple VFGs, including ptfA, fimA, hsf-2, exbB, exbD, tonB, hgbA, hgbB, fur, nanB, nanH, ompA, ompH, oma87, plpB, sodA, and sodC, was a broad characteristic of P. multocida strains from multiple host species (Figure 5) . However, several VFGs were only determined in the genome sequences of P. multocida from certain hosts. For example, toxA, a gene encoding a dermonecrotic toxin, was found only in strains from pig, sheep, and alpacas, while tbpA, a transferrin binding protein coding gene, was found only in strains from cattle, sheep, and alpacas ( Figure 5) . 

P. multocida is the causative agent of multiple diseases with a wide spectrum of host species, including humans and other primates (7) (8) (9) . In addition, P. multocida isolates recovered from different hosts with different diseases can be classified in many different serovars/genotypes according to different typing systems (7, 9) . Relying on only one or two typing systems is difficult to address the characteristics of P. multocida isolates from different host species and/or their association with different diseases. For example, P. multocida isolates from different host species might have the same capsular genotypes but possess different LPS genotypes and/or MLST genotypes; even those from different host species that share the same capsular, LPS, and MLST genotypes might carry different VFGs (27, 33) . Therefore, we have proposed a combined "capsular: LPS: MLST" genotyping system that includes virulence genotyping to discriminate P. multocida isolates from different host species and/or those associated with different diseases (7) . However, this combined genotyping system is multiplex PCR-based and is laborious and time-consuming.

Advances in bioinformatics and bioinformatical tools enable the application of whole genome sequence data for inclusion of various demographic information for bacterial characterization, such as capsular and LPS genotyping; the presence of adhesins, toxins, or other virulence factors (34) . In the present study, we reported the development of a genotyping platform for distinguishing P. multocida isolates according to the bacterial whole genome sequences. Validation of the PmGT platform was performed on a collection of P. multocida isolates from our laboratory. Results revealed that this genotyping system provides consistent results of determining the capsular-, LPS-, MLST genotypes, and VFGs, as compared with that obtained using multiplex PCR-based typing systems. Compared to the multiplex PCR-based typing systems (18, 19, 21, 22) and traditional serological typing systems (13, 16) , this genotyping system takes less time to yield results and does not require high-quality antisera, which represents a more efficient and cost-saving tool for characterizing P. multocida isolates in both epidemiological studies and clinical settings.

By using PmGT, the capsular-, LPS-, MLST genotypes, and VFGs of P. multocida strains from different hosts were determined according to the whole genome sequences. These results agree with those of the epidemiological studies (23, 24, 26, 35) . For example, P. multocida serovars B: 2 and A: 3 strains are frequently associated with bovine haemorrhagic septicaemia and respiratory diseases, respectively (36, 37) . It is known that P. multocida serogroups A and B are assigned to capsular genotypes A and B by multiplex PCR, respectively (18); while P. multocida Heddleston serovars 2 and 3 are assigned to LPS genotypes L2 and L3 by multiplex PCR, respectively (19) . That is why the capsular: LPS genotypes of most of the bovine strains were determined as A: L3 and B: L2, respectively. In addition, P. multocida strains isolated from bovine haemorrhagic septicaemia are commonly determined as ST122 (38), this sequence type can be reassigned to ST44 by using the multihost MLST database (27) . These findings could explain why P. multocida strains associated with bovine haemorrhagic septicaemia were typed as capsular: LPS: MLST genotype B: L2: ST44. Similar findings were also observed in P. multocida strains from the other host species. In particularly, most of the P. multocida strains from pigs were determined as capsular: LPS: MLST genotypes D: L6: ST11, A: L3: ST3, and A: L6: ST10. These results are also in agreement with the results of our previously epidemiological study (23) , suggesting that these three genotypes, particularly genotype D/L6/ST11, are likely to be strongly associated with swine respiratory diseases. However, during our test we also found the capsular-, LPS-, and/or MLST-genotypes of several strains could not be determined by PmGT according to the whole genome sequences. After check the data we put forward several reasons to explain this result: (1) most of these nontypeable genomes are sequenced and assembled using the second-generation sequencing technologies and the quality of these genomes are not high, some of the genes used for capsular/LPS/MLST genotyping fell within the gaps between genome contigs in the assemblies (7); (2) the genome sequences might be those of the capsular nontypeable strains reported (23, 39) ; (3) several strains belong to novel sequence types and the current Pasteurella multocida MLST database do not include these sequence types.

In conclusion, we developed an online platform for P. multocida genotyping (PmGT platform), which combines whole genome sequence analysis tools with web 2.0 technologies. By using this system, we determined the genotypes of P. multocida isolates from different host species. Overall, this system represents a more convenient tool for P. multocida diagnosis in both epidemiological studies and clinical settings. More importantly, our study provides an example to develop rapid and efficient tools for bacterial diagnosis by using their whole genome sequences in the coming age of artificial intelligence.

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https:// www.ncbi.nlm.nih.gov/genbank/, MN938443-MN938455 and MT570166∼MT570166.

ZP, BWi, JW, and BWu contributed to conception and design of the study. ZP, JL, WL, FW, LW, XW, and LH performed the experiments. ZP, JL, LW, and JW performed the statistical analysis. ZP wrote the first draft of the manuscript. ZP, HC, BWi, JW, and BWu revised the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version. 

We sincerely acknowledge Huazhong Agricultural University College of Informatics for providing the CentOS server and the PubMLST database for the RESTful port for connection.

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

Why microbial diagnostics need more than money

Advances in the diagnosis and treatment of Clostridium difficile infections

The GenMark ePlex((R)): another weapon in the syndromic arsenal for infection diagnosis

The diagnosis of infectious diseases by whole genome next generation sequencing: a new era is opening

A genomic perspective on the origin and emergence of SARS-CoV-2

Rapid whole-genome sequencing of bacterial pathogens in the clinical microbiology laboratory-pipe dream or reality?

Pasteurella multocida: genotypes and genomics

Pasteurella multocida: diseases and pathogenesis

Pasteurella multocida: from zoonosis to cellular microbiology

Dog licks baby Baby gets Pasteurella multocida meningitis

Pasteurella multocida from a dog causing Ludwig's angina

Beware of dogs licking ears

Studies on Pasteurella multocida. I A hemagglutination test for the identification of serological types

A new serological type of Pasteurella multocida from Central Africa

A new capsule serogroup of Pasteurella multocida

Fowl cholera: gel diffusion precipitin test for serotyping Pasteruella multocida from avian species

Molecular typing methods for Pasteurella multocida-a review

Genetic organization of Pasteurella multocida cap loci and development of a multiplex capsular PCR typing system

Development of a rapid multiplex PCR assay to genotype Pasteurella multocida strains by use of the lipopolysaccharide outer core biosynthesis locus

Characterisation of bovine strains of Pasteurella multocida and comparison with isolates of avian, ovine and porcine origin

Development of a multi-locus sequence typing scheme for avian isolates of Pasteurella multocida

Characterization of Pasteurella multocida associated with ovine pneumonia using multi-locus sequence typing (MLST) and virulenceassociated gene profile analysis and comparison with porcine isolates

A capsule/lipopolysaccharide/MLST genotype D/L6/ST11 of Pasteurella multocida is likely to be strongly associated with swine respiratory disease in China

Investigation of genetic diversity and epidemiological characteristics of Pasteurella multocida isolates from poultry in southwest China by population structure, multi-locus sequence typing and virulence-associated gene profile analysis

Virulence gene profiling of porcine Pasteurella multocida isolates of Assam

Characterization of Pasteurella multocida involved in rabbit infections

Genetic and phylogenetic characteristics of Pasteurella multocida isolates from different host species

Genomic characterization of Pasteurella multocida HB01, a serotype A bovine isolate from China

Complete genome sequence of Pasteurella multocida HN06, a toxigenic strain of serogroup D

Experimental pathogenicity and complete genome characterization of a pig origin Pasteurella multocida serogroup F isolate HN07

Occurrence of virulence factors and antimicrobial resistance in Pasteurella multocida strains isolated from slaughter cattle in Iran

Comparative genomic analysis of Asian Haemorrhagic Septicaemia-associated strains of Pasteurella multocida identifies more than 90 Haemorrhagic Septicaemia-specific genes

Virulence gene profiling and ompA sequence analysis of Pasteurella multocida and their correlation with host species

Evolutionary history of the global emergence of the escherichia coli epidemic clone ST131

Virulence genotype of Pasteurella multocida strains isolated from different hosts with various disease status

A review of hemorrhagic septicemia in cattle and buffalo

Isolation and antimicrobial susceptibilities of bacterial pathogens from bovine pneumonia: 1994-2002

Multilocus sequence typing of a global collection of Pasteurella multocida isolates from cattle and other host species demonstrates niche association

Isolation, antimicrobial resistance, and virulence genes of Pasteurella multocida strains from swine in China

Supplementary Text 1 | Nucleotide sequences and their GenBank accession numbers for the construction of a comparative database for P. multocida genotyping.Supplementary Figure S1 | PCR results of capsular genotypes of tested P. multocida strains in the present study.Supplementary Figure S2 | PCR results of LPS genotypes of tested P. multocida strains in the present study.Supplementary Table S1 | Pasteurella multocida genome sequences used in the present study.

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.