Genome-Wide Analyses Reveal Genes Subject to Positive Selection in Pasteurella multocida

Comparative genome analysis of Pasteurella multocida strains of porcine origin

Pasteurella multocida causes acute/chronic pasteurellosis in porcine, resulting in considerable economic losses globally. The draft genomes of two Indian strains NIVEDIPm17 (serogroup D) and NIVEDIPm36 (serogroup A) were sequenced. A total of 2182-2284 coding sequences (CDSs) were predicted along with 5-6 rRNA and 45-46 tRNA genes in the genomes. Multilocus sequence analysis and LPS genotyping showed the presence of ST50: genotype 07 and ST74: genotype 06 in NIVEDIPm17 and NIVEDIPm36, respectively. Pangenome analysis of 61 strains showed the presence of 1653 core genes, 167 soft core genes, 750 shell genes, and 1820 cloud genes. Analysis of virulence-associated genes in 61 genomes indicated the presence of nanB, exbB, exbD, ptfA, ompA, ompH, fur, plpB, fimA, sodA, sodC, tonB, and omp87 in all strains. The 61 genomes contained genes encoding tetracycline (54%), streptomycin (48%), sulphonamide (28%), tigecycline (25%), chloramphenicol (21%), amikacin (7%), cephalosporin (5%), and trimethoprim (5%) resistance. Multilocus sequence type revealed that ST50 was the most common (34%), followed by ST74 (26%), ST13 (24%), ST287 (5%), ST09 (5%), ST122 (3%), and ST07 (2%). Single-nucleotide polymorphism and core genome-based phylogenetic analysis clustered the strains into three major clusters. In conclusion, we described the various virulence factors, mobile genetic elements, and antimicrobial resistance genes in the pangenome of P. multocida of porcine origin, besides the rare presence of LPS genotype 7 in serogroup D.

Pmorf0222, a Virulence Factor in Pasteurella multocida, Activates Nuclear Factor Kappa B and Mitogen-Activated Protein Kinase via Toll-Like Receptor 1/2

Pasteurella multocida primarily causes hemorrhagic septicemia and pneumonia in poultry and livestock. Identification of the relevant virulence factors is therefore essential for understanding its pathogenicity. Pmorf0222, encoding the PM0222 protein, is located on a specific prophage island of the pathogenic strain C48-1 of P. multocida. Its role in the pathogenesis of P. multocida infection is still unknown. The proinflammatory cytokine plays an important role in P. multocida infection; therefore, murine peritoneal exudate macrophages were treated with the purified recombinant PM0222, which induced the secretion of tumor necrosis factor alpha (TNF-α) and interleukin-1β (IL-1β) via the Toll-like receptor 1/2 (TLR1/2)-nuclear factor kappa B (NF-κB)/mitogen-activated protein kinase (MAPK) signaling and inflammasome activation. Additionally, the mutant strain and complemented strain were evaluated in the mouse model with P. multocida infection, and PM0222 was identified as a virulence factor, which was secreted by outer membrane vesicles of P. multocida. Further results revealed that Pmorf0222 affected the synthesis of the capsule, adhesion, serum sensitivity, and biofilm formation. Thus, we identified Pmorf0222 as a novel virulence factor in the C48-1 strain of P. multocida, explaining the high pathogenicity of this pathogenic strain.

Experimental pathogenicity and comparative genome analysis of high- and low-virulence strains of rabbit-origin Pasteurella multocida

Pasteurella multocida, the causative pathogen of rabbit pasteurellosis, causes significant economic losses in the commercial rabbit industry. However, the associated pathogenic mechanism of P. multocida remains unclear. The aim of this study is to compare the genomes and pathogenicity of high- and low-virulence strains of P. multocida to advance the current understanding of rabbit pasteurellosis. The high-virulence strain rapidly proliferates in the lung and spleen of infected mice within approximately 9 h, maintaining a high bacterial load until host death. Meanwhile, the low-virulence strain only proliferates in mouse organs for a short time, with the bacterial load beginning to decrease 13 h post-infection. Moreover, the expressions of inflammatory cytokines MCP-1, TNF-α, and IL-1β are upregulated in all infected mouse lung and spleen tissue, however, the high-virulence strain induced significantly higher expression than the low-virulence strain. Histopathological analysis revealed greater inflammation and tissue lesions in the lung and spleen of mice infected with the high-virulence strain. Two pathogenicity-associated regions unique to the genome of the high-virulence strain harbor approximately 199 genes, including functional genes related to virulence factors, such as lipopolysaccharide biosynthesis, iron acquisition, biosynthesis of outer membrane proteins, and adhesion. These two genomic regions are shared by three previously sequenced, highly virulent P. multocida strains in rabbits. In conclusion, the increased pathogenicity of high-virulence P. multocida may be due to the presence of virulence-associated genes in two unique genomic regions, resulting in strong proliferative activity, significant inflammation, and pathological lesions in the mouse model. These findings provide important insights regarding the pathogenic mechanism underlying rabbit pasteurellosis.

Comparative genome analysis of Pasteurella multocida serogroup B:2 strains causing haemorrhagic septicaemia (HS) in bovines

Pasteurella multocida, a Gram-negative bacterium with ubiquitous nature, is known to affect wide range of host species worldwide with varied clinical manifestations including haemorrhagic septicaemia (HS) in bovines. Although, HS causing P. multocida strains were identified and characterized by conventional tools and PCR assays, diverse strains are indistinguishable by these tools in the face of disease outbreaks. In this study, draft genomes of three virulent P. multocida serotype B:2 strains (NIVEDIPm32, NIVEDIPm34 and NIVEDIPm35) were analyzed following whole genome sequencing, assembly, annotation and compared them with existing global genomes (n = 43) of bovine origin in the database. Three draft genomes of NIVEDIPm strains consisted of 40-52 contigs with GC content of ∼40.4%. The genome size and predicted genes content was ∼2.3 Mb and 2181-2189, respectively. Besides, the presence of various mobile genetic elements, antimicrobial resistance genes and biofilm related genes suggested their vital roles in virulence; further, adaptation to the host immune system as well as host pathogen interaction. Multi locus sequence analysis based on RIRDC scheme showed the presence of ST122 in all the three strains. wgMLST based phylogenic analysis suggested that HS causing Indian virulent field strains differed geographically and showed diversity from existing HS vaccine strain P52. The phylogenetic tree revealed that North Indian strains share high similarity with strains of Pakistan than South Indian Strain. Notably, a high divergence of SNPs between the HS causing circulating virulent strains of India and current HS vaccine strain P52 suggested an imminent need for relook in to HS vaccination strategy for livestock in India.

A genome-wide scan of wastewater E. coli for genes under positive selection: focusing on mechanisms of antibiotic resistance

Antibiotic resistance is a global health threat and consequently, there is a need to understand the mechanisms driving its emergence. Here, we hypothesize that genes and mutations under positive selection may contribute to antibiotic resistance. We explored wastewater E. coli, whose genomes are highly diverse. We subjected 92 genomes to a statistical analysis for positively selected genes. We obtained 75 genes under positive selection and explored their potential for antibiotic resistance. We found that eight genes have functions relating to antibiotic resistance, such as biofilm formation, membrane permeability, and bacterial persistence. Finally, we correlated the presence/absence of non-synonymous mutations in positively selected sites of the genes with a function in resistance against 20 most prescribed antibiotics. We identified mutations associated with antibiotic resistance in two genes: the porin ompC and the bacterial persistence gene hipA. These mutations are located at the surface of the proteins and may hence have a direct effect on structure and function. For hipA, we hypothesize that the mutations influence its interaction with hipB and that they enhance the capacity for dormancy as a strategy to evade antibiotics. Overall, genomic data and positive selection analyses uncover novel insights into mechanisms driving antibiotic resistance.

Global Genome Diversity and Recombination in Mycoplasma pneumoniae

Genomic changes in Mycoplasma pneumoniae caused by adaptation to environmental or ecologic pressures are poorly understood. We collected M. pneumoniae from children who had confirmed pneumonia in Taiwan during 2017-2020. We used whole-genome sequencing to compare these isolates with a worldwide collection of current and historical clinical strains for characterizing population structures. A phylogenetic tree for 284 strains showed that all sequenced strains consisted of 5 clades: T1-1 (sequence type [ST]1), T1-2 (mainly ST3), T1-3 (ST17), T2-1 (mainly ST2), and T2-2 (mainly ST14). We identified a putative recombination block containing 6 genes (MPN366‒371). Macrolide resistance involving 23S rRNA mutations was detected for each clade. Clonal expansion of macrolide resistance occurred mostly within subtype 1 strains, of which clade T1-2 showed the highest recombination rate and genome diversity. Functional characterization of recombined regions provided clarification of the biologic role of these recombination events in the evolution of M. pneumoniae.

First report from Bangladesh on genetic diversity of multidrug-resistant Pasteurella multocida type B:2 in fowl cholera

Background and Aim:

Fowl cholera (FC) caused by Pasteurella multocida is a highly contagious bacterial disease of global importance for poultry production. The severity and incidence of FC caused by P. multocida may vary considerably depending on several factors associated with the host (including species and age of infected birds), the environment, and the bacterial strain. This study aimed to investigate the genetic diversity of multidrug-resistant P. multocida strains isolated from FC outbreaks in laying hens from commercial farms of Bangladesh.

Materials and Methods:

We collected 57 samples of suspected FC, including 36 live and 21 dead laying hens. P. multocida isolates were characterized by biochemical and molecular-biological methods.

Results:

Twenty-two strains of P. multocida were isolated from these samples through phenotypic and genotypic characterization. The strains were grouped into two distinct random amplification of polymorphic DNA (RAPD) biotypes harboring a range of pathogenic genes; exbB, ompH, ptfA, nanB, sodC, and hgbA. In this study, 90.90% and 81.82% P. multocida strains were multidrug-resistant and biofilm formers, respectively. Whole-genome sequencing of the two representative RAPD phylotypes confirmed as P. multocida type B: L2:ST122, harboring a number of virulence factors-associated genes (VFGs), and antimicrobial resistance (AMR) genes (ARGs). In addition, pan-genome analysis revealed 90 unique genes in the genomes of P. multocida predicted to be associated with versatile metabolic functions, pathogenicity, virulence, and AMR.

Conclusion:

This is first-ever report on the association of P. multocida genotype B: L2:ST122 and related VFGs and ARGs in the pathogenesis of FC in laying hens. This study also provides a genetic context for future researches on the evolutionary diversity of P. multocida strains and their host adaptation.

Anti-Foc RT4 Activity of a Newly Isolated Streptomyces sp. 5-10 From a Medicinal Plant (Curculigo capitulata)

Fusarium wilt of banana caused by Fusarium oxysporum f. sp. cubense (Foc) is a disastrous soil-borne fungal disease. Foc tropical race 4 (Foc TR4) can infect almost all banana cultivars. Until now, there is a shortage of safety and effective control methods and commercial banana cultivars with a resistance against Foc TR4. Biocontrol using environmentally friendly microbes is a promising strategy for the management of Foc TR4. Here, a strain 5-10, newly isolated from a medicinal plant (Curculigo capitulata), exhibited a high antifungal activity against Foc TR4. Combing the morphological characteristics and molecular identification, strain 5-10 was classified as a Streptomyces genus. The sequenced genome revealed that more than 39 gene clusters were involved in the biosynthesis of secondary metabolites. Some multidrug resistance gene clusters were also identified such as mdtD, vatB, and vgaE. To improve the anti-Foc TR4 activity of the strain 5-10 extracts, an optimization method of fermentation broth was established. Antifungal activity increased by 72.13% under the fermentation system containing 2.86 g/L of NaCl and 11.57% of inoculation amount. After being treated with the strain 5-10 extracts, the Foc TR4 hyphae shrinked, deformed, and ruptured. The membrane integrity and cell ultrastructure incurred irreversible damage. Streptomyces sp. 5-10 extracts play a fungicidal role in Foc TR4. Hence, Streptomyces sp. 5-10 will be a potential biocontrol agent to manage fungal diseases by exploring the microbial fertilizer.

Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection

Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.

Pasteurella multocida specific bacteriophage suppresses P. multocida-induced inflammation: identification of genes related to bacteriophage signaling by Pasteurella multocida-infected swine nasal turbinate cells

BACKGROUND

Although Pasteurella multocida is highly prevalent pathogen in animals and plays an important role in swine respiratory diseases, only a few studies on the use of bacteriophages specific to Pasteurella multocida disease have been reported.

OBJECTIVE

The object of this study was to investigate the therapeutic effect of specific P. multocida bacteriophages and to identify genes related to bacteriophage signaling utilizing RNA microarrays in swine nasal turbinate cells.

METHODS

Pas-MUP-1 phages were applied 24 h prior to P. multocida infection (1 × 10⁷ cfu/ml) at several concentrations of bacterial infection. Cells were incubated to detect cytokines and 24 h to detect mucin production. And real-time quantitative PCR was performed to examine related genes expression. To determine the change of total gene expression based on P. multocida and Pas-MUP-1 treatment, we performed RNA sequencing experiments.

RESULTS

We found that P. multocida-infected PT-K75 cells show increased gene expression of IL-1β, IL-6, and Muc1 in a dose-dependent manner. Interestingly, these genes resulted in decreased expression in P. multocida pretreated with the P. multocida-specific Pas-MUP-1 bacteriophage. RNA sequencing analysis revealed that bacteriophage administration regulated genes associated with immune and inflammatory responses, and the regulated genes were dramatically concentrated in the cytokine/chemokine-based signaling pathways. Pas-MUP-1 treatment was shown to regulate P. multocida induced gene expression in the bacteria.

CONCLUSION

These results suggest the specific bacteriophage has therapeutic potential as an alternative to antibiotic treatment to defend against P. multocida infection by altering inflammatory gene expression profiles.

Pasteurella multocida: Genotypes and Genomics

Pasteurella multocida is a highly versatile pathogen capable of causing infections in a wide range of domestic and wild animals as well as in humans and nonhuman primates. Despite over 135 years of research, the molecular basis for the myriad manifestations of P. multocida pathogenesis and the determinants of P. multocida phylogeny remain poorly defined. The current availability of multiple P. multocida genome sequences now makes it possible to delve into the underlying genetic mechanisms of P. multocida fitness and virulence. Using whole-genome sequences, the genotypes, including the capsular genotypes, lipopolysaccharide (LPS) genotypes, and multilocus sequence types, as well as virulence factor-encoding genes of P. multocida isolates from different clinical presentations can be characterized rapidly and accurately. Putative genetic factors that contribute to virulence, fitness, host specificity, and disease predilection can also be identified through comparative genome analysis of different P. multocida isolates. However, although some knowledge about genotypes, fitness, and pathogenesis has been gained from the recent whole-genome sequencing and comparative analysis studies of P. multocida, there is still a long way to go before we fully understand the pathogenic mechanisms of this important zoonotic pathogen. The quality of several available genome sequences is low, as they are assemblies with relatively low coverage, and genomes of P. multocida isolates from some uncommon host species are still limited or lacking. Here, we review recent advances, as well as continuing knowledge gaps, in our understanding of determinants contributing to virulence, fitness, host specificity, disease predilection, and phylogeny of P. multocida.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Genomic analysis of Leptospira interrogans serovar Paidjan and Dadas isolates from carrier dogs and comparative genomic analysis to detect genes under positive selection

BACKGROUND

Leptospirosis is an emerging infectious disease worldwide that can cause high morbidity and mortality rates in humans and animals. The causative spirochetes have reservoirs in mammalian hosts, but there has been limited analysis of the genomes of isolates recovered from animals. The aims of this study were to characterize genomic features of two Leptospira interrogans strains recently isolated from asymptomatic dogs in Thailand (strains CUDO5 and CDUO8), and to perform comparative genome analyses with other strains. Molecular adaptive evolution in L. interrogans as signaled by positive selection also was analyzed.

RESULTS

Whole genome sequence analysis revealed that strains CUDO5 and CUDO8 had genome sizes of approximately 4.9 Mbp with 35.1% GC contents. Using monoclonal antibodies, strains CUDO5 and CUDO8 were identified as serovars Paidjan and Dadas, respectively. These strains harbored genes known to be associated with acute and chronic infections. Using Single Nucleotide Polymorphisms phylogeny (SNPs) with 97 L. interrogans strains, CUDO5 and CUDO8 had closest genetic relatedness with each other. Nevertheless, the serovar determinant region (rfb locus) showed variations in the genes encoding sugar biosynthesis. Amongst 13 representative L. interrogans strains examined for molecular adaptive evolution through positive selection under the site-model of Phylogenetic Analysis of Maximum Likelihood, genes responsible for iron acquisition (tlyA and hbpA), motility (fliN2, flgK, and flhB) and thermal adaptation (lpxD1) were under increased selective pressure.

CONCLUSIONS

L. interrogans serovar Paidjan strain CUDO5 and serovar Dadas strain CUDO8 had close genetic relatedness as analyzed by SNPs phylogeny. They contained genes with established roles in acute and chronic leptospirosis. The rfb locus in both serovars showed gene variation associated with sugar biosynthesis. Positive selection analysis indicated that genes encoding factors involved in motility, temperature adaptation, and iron acquisition were under strong positive selection in L. interrogans. These may be associated with adaptation in the early stages of infection.

Comparative analysis reveals the Genomic Islands in Pasteurella multocida population genetics: on Symbiosis and adaptability

BACKGROUND

Pasteurella multocida (P. multocida) is a widespread opportunistic pathogen that infects human and various animals. Genomic Islands (GIs) are one of the most important mobile components that quickly help bacteria acquire large fragments of foreign genes. However, the effects of GIs on P. multocida are unknown in the evolution of bacterial populations.

RESULTS

Ten avian-sourced P. multocida obtained through high-throughput sequencing together with 104 publicly available P. multocida genomes were used to analyse their population genetics, thus constructed a pan-genome containing 3948 protein-coding genes. Through the pan-genome, the open evolutionary pattern of P. multocida was revealed, and the functional components of 944 core genes, 2439 accessory genes and 565 unique genes were analysed. In addition, a total of 280 GIs were predicted in all strains. Combined with the pan-genome of P. multocida, the GIs accounted for 5.8% of the core genes in the pan-genome, mainly related to functional metabolic activities; the accessory genes accounted for 42.3%, mainly for the enrichment of adaptive genes; and the unique genes accounted for 35.4%, containing some defence mechanism-related genes.

CONCLUSIONS

The effects of GIs on the population genetics of P. multocida evolution and adaptation to the environment are reflected by the proportion and function of the pan-genome acquired from GIs, and the large quantities of GI data will aid in additional population genetics studies.

Genome rearrangements and selection in multi-chromosome bacteria Burkholderia spp

BACKGROUND

The genus Burkholderia consists of species that occupy remarkably diverse ecological niches. Its best known members are important pathogens, B. mallei and B. pseudomallei, which cause glanders and melioidosis, respectively. Burkholderia genomes are unusual due to their multichromosomal organization, generally comprised of 2-3 chromosomes.

RESULTS

We performed integrated genomic analysis of 127 Burkholderia strains. The pan-genome is open with the saturation to be reached between 86,000 and 88,000 genes. The reconstructed rearrangements indicate a strong avoidance of intra-replichore inversions that is likely caused by selection against the transfer of large groups of genes between the leading and the lagging strands. Translocated genes also tend to retain their position in the leading or the lagging strand, and this selection is stronger for large syntenies. Integrated reconstruction of chromosome rearrangements in the context of strains phylogeny reveals parallel rearrangements that may indicate inversion-based phase variation and integration of new genomic islands. In particular, we detected parallel inversions in the second chromosomes of B. pseudomallei with breakpoints formed by genes encoding membrane components of multidrug resistance complex, that may be linked to a phase variation mechanism. Two genomic islands, spreading horizontally between chromosomes, were detected in the B. cepacia group.

CONCLUSIONS

This study demonstrates the power of integrated analysis of pan-genomes, chromosome rearrangements, and selection regimes. Non-random inversion patterns indicate selective pressure, inversions are particularly frequent in a recent pathogen B. mallei, and, together with periods of positive selection at other branches, may indicate adaptation to new niches. One such adaptation could be a possible phase variation mechanism in B. pseudomallei.

Iron-associated protein interaction networks reveal the key functional modules related to survival and virulence of Pasteurella multocida

Pasteurella multocida causes respiratory infectious diseases in a multitude of birds and mammals. A number of virulence-associated genes were reported across different strains of P. multocida, including those involved in the iron transport and metabolism. Comparative iron-associated genes of P. multocida among different animal hosts towards their interaction networks have not been fully revealed. Therefore, this study aimed to identify the iron-associated genes from core- and pan-genomes of fourteen P. multocida strains and to construct iron-associated protein interaction networks using genome-scale network analysis which might be associated with the virulence. Results showed that these fourteen strains had 1587 genes in the core-genome and 3400 genes constituting their pan-genome. Out of these, 2651 genes associated with iron transport and metabolism were selected to construct the protein interaction networks and 361 genes were incorporated into the iron-associated protein interaction network (iPIN) consisting of nine different iron-associated functional modules. After comparing with the virulence factor database (VFDB), 21 virulence-associated proteins were determined and 11 of these belonged to the heme biosynthesis module. From this study, the core heme biosynthesis module and the core outer membrane hemoglobin receptor HgbA were proposed as candidate targets to design novel antibiotics and vaccines for preventing pasteurellosis across the serotypes or animal hosts for enhanced precision agriculture to ensure sustainability in food security.

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.

Pan-genomic approach shows insight of genetic divergence and pathogenic-adaptation of Pasteurella multocida

Pasteurella multocida is a gram-negative, non-motile bacterial pathogen, which is associated with chronic and acute infections as snuffles, pneumonia, atrophic rhinitis, fowl cholera and hemorrhagic septicemia. These diseases affect a wide range of domestic animals, leading to significant morbidity and mortality and causing significant economic losses worldwide. Due to the interest in deciphering the genetic diversity and process adaptive between P. multocida strains, this work aimed was to perform a pan-genome analysis to evidence horizontal gene transfer and positive selection among 23 P. multocida strains isolated from distinct diseases and hosts. The results revealed an open pan-genome containing 3585 genes and an accessory genome presenting 1200 genes. The phylogenomic analysis based on the presence/absence of genes and islands exhibit high levels of plasticity, which reflects a high intraspecific diversity and a possible adaptive mechanism responsible for the specific disease manifestation between the established groups (pneumonia, fowl cholera, hemorrhagic septicemia and snuffles). Additionally, we identified differences in accessory genes among groups, which are involved in sugar metabolism and transport systems, virulence-related genes and a high concentration of hypothetical proteins. However, there was no specific indispensable functional mechanism to decisively correlate the presence of genes and their adaptation to a specific host/disease. Also, positive selection was found only for two genes from sub-group hemorrhagic septicemia, serotype B. This comprehensive comparative genome analysis will provide new insights of horizontal gene transfers that play an essential role in the diversification and adaptation mechanism into P. multocida species to a specific disease.