Comparative genomics of quinolone-resistant Escherichia coli from broilers and humans in Norway.

BACKGROUND: The usage of fluoroquinolones in Norwegian livestock production is very low, including in broiler production. Historically, quinolone-resistant Escherichia coli (QREC) isolated from Norwegian production animals rarely occur. However, with the introduction of a selective screening method for QREC in the Norwegian monitoring programme for antimicrobial resistance in the veterinary sector in 2014; 89.5% of broiler caecal samples and 70.7% of broiler meat samples were positive. This triggered the concern if there could be possible links between broiler and human reservoirs of QREC. We are addressing this by characterizing genomes of QREC from humans (healthy carriers and patients) and broiler isolates (meat and caecum). RESULTS: The most frequent mechanism for quinolone resistance in both broiler and human E. coli isolates were mutations in the chromosomally located gyrA and parC genes, although plasmid mediated quinolone resistance (PMQR) was also identified. There was some relatedness of the isolates within human and broiler groups, but little between these two groups. Further, some overlap was seen for isolates with the same sequence type isolated from broiler and humans, but overall, the SNP distance was high. CONCLUSION: Based on data from this study, QREC from broiler makes a limited contribution to the incidence of QREC in humans in Norway.

Tenacibaculosis in Norwegian Atlantic salmon (Salmo salar) cage-farmed in cold sea water is primarily associated with Tenacibaculum finnmarkense genomovar finnmarkense.

Skin conditions associated with Tenacibaculum spp. constitute a significant threat to the health and welfare of sea-farmed Atlantic salmon (Salmo salar L.) in Norway. Fifteen presumptive tenacibaculosis outbreaks distributed along the Norwegian coast during the late winter and spring of 2018 were investigated. Bacteriological culture confirmed the presence of Tenacibaculum spp. Seventy-six isolates cultured from individual fish were selected and subjected to whole-genome sequencing and MALDI-TOF MS analysis. Average nucleotide identity and MALDI-TOF analyses confirmed the presence of T. finnmarkense and T. dicentrarchi, with further division of T. finnmarkense into genomovars (gv.) finnmarkense and ulcerans. Core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses identified the presence of a genetically conserved cluster of gv. finnmarkense isolates against a background of relatively genetically diverse gv. finnmarkense and gv. ulcerans isolates in 13 of the 15 studied cases. This clustering strongly suggests a link between T. finnmarkense gv. finnmarkense and development of clinical tenacibaculosis in sea-farmed Norwegian salmon in the late winter and spring. Analysis of 25 Tenacibaculum isolates collected during the spring of 2019 from similar cases identified a similar distribution of genotypes. Low water temperatures were common to all cases, and most incidences involved relatively small fish shortly after sea transfer, suggesting that these fish are particularly predisposed to Tenacibaculum infection.

Dissemination of Quinolone-Resistant Escherichia coli in the Norwegian Broiler and Pig Production Chains and Possible Persistence in the Broiler Production Environment.

In Norway, the use of quinolones in livestock populations is very low, and prophylactic use is prohibited. Despite this, quinolone-resistant Escherichia coli (QREC) isolates are present at low levels in several animal species. The source of these QREC isolates is unknown. The aim of this study was to characterize and compare QREC isolates from different animal species to identify putative factors that may promote the occurrence of QREC. A total of 280 QREC isolates, from broilers, pigs, red foxes, and wild birds, were whole-genome sequenced and analyzed. Well-known chromosomal and plasmid-mediated resistance mechanisms were identified. In addition, mutations in marR, marA, and rpoB causing novel amino acid substitutions in their respective proteins were detected. Phylogenetic analyses were used to determine the relationships between the isolates. Quinolone resistance mechanism patterns appeared to follow sequence type groups. Similar QREC isolates with similar resistance mechanism patterns were detected from the samples, and further phylogenetic analysis indicated close evolutionary relationships between specific isolates from different sources. This suggests the dissemination of highly similar QREC isolates between animal species and also the persistence of QREC strains within the broiler production chain. This highlights the importance of both control measures at the top of the production chain as well as biosecurity measures to avoid the further dissemination and persistence of QREC in these environments.IMPORTANCE Since antimicrobial usage is low in Norwegian animal husbandry, Norway is an ideal country to study antimicrobial resistance in the absence of selective pressure from antimicrobial usage. In particular, the usage of quinolones is very low, which makes it possible to investigate the spread and development of quinolone resistance in natural environments. Comparison of quinolone-resistant E. coli (QREC) isolates from livestock and wild animals in light of this low quinolone usage provides new insights into the development and dissemination of QREC in both natural and production environments. With this information, preventive measures may be taken to prevent further dissemination within Norwegian livestock and between other animals, thus maintaining the favorable situation in Norway.

Tenacibaculum piscium sp. nov., isolated from skin ulcers of sea-farmed fish, and description of Tenacibaculum finnmarkense sp. nov. with subdivision into genomovars finnmarkense and ulcerans.

Results of previous multilocus sequence and whole-genome-based analyses have suggested that a homogeneous group of isolates belonging to the genus Tenacibaculum, represented by strain TNO020T and associated with skin ulcer development in sea-farmed fish, represents an as-yet-undescribed species. Comparative whole-genome analysis performed in the present study clustered five isolates, including TNO020T, in a distinct lineage within the genus Tenacibaculum. Phenotypic differences, high intra-cluster average nucleotide identity (ANI) values and low ANI values with other Tenacibaculum species support the proposal of a novel species, for which we propose the name Tenacibaculum piscium sp. nov. with strain TNO020T (=CCUG 73833T=NCIMB 15240T) as the type strain. Further, large-scale genome analyses confirmed the existence of two different phylogenetic lineages within 'T. finnmarkense', a species effectively but not validly published previously. ANI values just above the species delineation threshold of 95-96 % confirmed that both lineages belong to the same species. This result was also supported by DNA-DNA hybridization values. Phenotypically, the two conspecific lineages are distinguishable by differences in growth temperature range and ability to degrade l-proline. For the group of isolates already commonly known as 'T. finnmarkense', we propose the name Tenacibaculum finnmarkense sp. nov., with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain. We further propose the subdivision of T. finnmarkense sp. nov. into two genomovars, T. finnmarkense genomovar finnmarkense with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain and T. finnmarkense genomovar ulcerans with strain TNO010T (=CCUG 73832T=NCIMB 15239T) as the type strain.

Multilocus Variable-Number Tandem-Repeat Analysis of Yersinia ruckeri Confirms the Existence of Host Specificity, Geographic Endemism, and Anthropogenic Dissemination of Virulent Clones.

A multilocus variable-number tandem-repeat analysis (MLVA) assay was developed for epizootiological study of the internationally significant fish pathogen Yersinia ruckeri, which causes yersiniosis in salmonids. The assay involves amplification of 10 variable-number tandem-repeat (VNTR) loci in two five-plex PCRs, followed by capillary electrophoresis. A collection of 484 Y. ruckeri isolates, originating from various biological sources and collected from four continents over 7 decades, was analyzed. Minimum-spanning-tree cluster analysis of MLVA profiles separated the studied population into nine major clonal complexes and a number of minor clusters and singletons. The major clonal complexes could be associated with host species, geographic origin, and serotype. A single large clonal complex of serotype O1 isolates dominating the yersiniosis situation in international rainbow trout farming suggests anthropogenic spread of this clone, possibly related to transport of fish. Moreover, subclustering within this clonal complex indicates putative transmission routes and multiple biotype shift events. In contrast to the situation in rainbow trout, Y. ruckeri strains associated with disease in Atlantic salmon appear as more or less geographically isolated clonal complexes. A single complex of serotype O1 exclusive to Norway was found to be responsible for almost all major yersiniosis outbreaks in modern Norwegian salmon farming, and site-specific subclustering further indicates persistent colonization of freshwater farms in Norway. Identification of genetically diverse Y. ruckeri isolates from clinically healthy fish and environmental sources also suggests the widespread existence of less-virulent or avirulent strains.IMPORTANCE This comprehensive population study substantially improves our understanding of the epizootiological history and nature of an internationally important fish-pathogenic bacterium. The MLVA assay developed and presented represents a high-resolution typing tool particularly well suited for Yersinia ruckeri infection tracing, selection of strains for vaccine inclusion, and risk assessment. The ability of the assay to separate isolates into geographically linked and/or possibly host-specific clusters reflects its potential utility for maintenance of national biosecurity. The MLVA is internationally applicable and robust, and it provides clear, unambiguous, and easily interpreted results. Typing is reasonably inexpensive, with a moderate technological requirement, and may be completed from a harvested colony within a single working day. As the resulting MLVA profiles are readily portable, any Y. ruckeri strain may rapidly be placed in a global epizootiological context.

NucDiff: in-depth characterization and annotation of differences between two sets of DNA sequences.

BACKGROUND: Comparing sets of sequences is a situation frequently encountered in bioinformatics, examples being comparing an assembly to a reference genome, or two genomes to each other. The purpose of the comparison is usually to find where the two sets differ, e.g. to find where a subsequence is repeated or deleted, or where insertions have been introduced. Such comparisons can be done using whole-genome alignments. Several tools for making such alignments exist, but none of them 1) provides detailed information about the types and locations of all differences between the two sets of sequences, 2) enables visualisation of alignment results at different levels of detail, and 3) carefully takes genomic repeats into consideration. RESULTS: We here present NucDiff, a tool aimed at locating and categorizing differences between two sets of closely related DNA sequences. NucDiff is able to deal with very fragmented genomes, repeated sequences, and various local differences and structural rearrangements. NucDiff determines differences by a rigorous analysis of alignment results obtained by the NUCmer, delta-filter and show-snps programs in the MUMmer sequence alignment package. All differences found are categorized according to a carefully defined classification scheme covering all possible differences between two sequences. Information about the differences is made available as GFF3 files, thus enabling visualisation using genome browsers as well as usage of the results as a component in an analysis pipeline. NucDiff was tested with varying parameters for the alignment step and compared with existing alternatives, called QUAST and dnadiff. CONCLUSIONS: We have developed a whole genome alignment difference classification scheme together with the program NucDiff for finding such differences. The proposed classification scheme is comprehensive and can be used by other tools. NucDiff performs comparably to QUAST and dnadiff but gives much more detailed results that can easily be visualized. NucDiff is freely available on https://github.com/uio-cels/NucDiff under the MPL license.

The genome sequence of Atlantic cod reveals a unique immune system.

Atlantic cod (Gadus morhua) is a large, cold-adapted teleost that sustains long-standing commercial fisheries and incipient aquaculture. Here we present the genome sequence of Atlantic cod, showing evidence for complex thermal adaptations in its haemoglobin gene cluster and an unusual immune architecture compared to other sequenced vertebrates. The genome assembly was obtained exclusively by 454 sequencing of shotgun and paired-end libraries, and automated annotation identified 22,154 genes. The major histocompatibility complex (MHC) II is a conserved feature of the adaptive immune system of jawed vertebrates, but we show that Atlantic cod has lost the genes for MHC II, CD4 and invariant chain (Ii) that are essential for the function of this pathway. Nevertheless, Atlantic cod is not exceptionally susceptible to disease under natural conditions. We find a highly expanded number of MHC I genes and a unique composition of its Toll-like receptor (TLR) families. This indicates how the Atlantic cod immune system has evolved compensatory mechanisms in both adaptive and innate immunity in the absence of MHC II. These observations affect fundamental assumptions about the evolution of the adaptive immune system and its components in vertebrates.

Draft genome sequences of 13 Tenacibaculum maritimum isolates from farmed Norwegian Cyclopterus lumpus (lumpfish) and Scophthalmus maximus (turbot).

Thirteen bacterial isolates of Tenacibaculum maritimum were sequenced and assembled. The strains were isolated from four disease outbreaks in farmed marine fish in Norway. Eight isolates were from Cyclopterus lumpus (lumpfish), and five were from Scophthalmus maximus (turbot). Overall, sequence similarity did not correlate with host species or geographic location.

Polyether ionophore resistance in a one health perspective.

Antimicrobial resistance is a major threat to human health and must be approached from a One Health perspective. Use of antimicrobials in animal husbandry can lead to dissemination and persistence of resistance in human pathogens. Polyether ionophores (PIs) have antimicrobial activities and are among the most extensively used feed additives for major production animals. Recent discoveries of genetically encoded PI resistance mechanisms and co-localization of resistance mechanisms against PIs and antimicrobials used in human medicine on transferrable plasmids, have raised concerns that use of PIs as feed additives bear potential risks for human health. This review summarizes the current knowledge on PI resistance and discusses the potential consequences of PI-usage as feed additives in a One Health perspective.

Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV).

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

Comparisons of infant Escherichia coli isolates link genomic profiles with adaptation to the ecological niche.

BACKGROUND: Despite being one of the most intensely studied model organisms, many questions still remain about the evolutionary biology and ecology of Escherichia coli. An important step toward achieving a more complete understanding of E.coli biology entails elucidating relationships between gene content and adaptation to the ecological niche. RESULTS: Here, we present genome comparisons of 16 E.coli strains that represent commensals and pathogens isolated from infants during a specific time period in Trondheim, Norway. Using differential gene content, we characterized enrichment profiles of the collection of strains relating to phylogeny, early vs. late colonization, pathogenicity and growth rate. We found clear gene content distinctions relating to the various grouping criteria. We also found that different categories of strains use different genetic elements for similar biological processes. The sequenced genomes included two pairs of strains where each pair was isolated from the same infant at different time points. One pair, in which the strains were isolated four months apart, showed maintenance of an early colonizer genome profile but also gene content and codon usage changes toward the late colonizer profile. Lastly, we placed our sequenced isolates into a broader genomic context by comparing them with 25 published E.coli genomes that represent a variety of pathotypes and commensal strains. This analysis demonstrated the importance of geography in shaping strain level gene content profiles. CONCLUSIONS: Our results indicate a general pattern where alternative genetic pathways lead toward a consistent ecological role for E.coli as a species. Within this framework however, we saw selection shaping the coding repertoire of E.coli strains toward distinct ecotypes with different phenotypic properties.

From green to red: horizontal gene transfer of the phycoerythrin gene cluster between Planktothrix strains.

Horizontal gene transfer is common in cyanobacteria, and transfer of large gene clusters may lead to acquisition of new functions and conceivably niche adaption. In the present study, we demonstrate that horizontal gene transfer between closely related Planktothrix strains can explain the production of the same oligopeptide isoforms by strains of different colors. Comparison of the genomes of eight Planktothrix strains revealed that strains producing the same oligopeptide isoforms are closely related, regardless of color. We have investigated genes involved in the synthesis of the photosynthetic pigments phycocyanin and phycoerythrin, which are responsible for green and red appearance, respectively. Sequence comparisons suggest the transfer of a functional phycoerythrin gene cluster generating a red phenotype in a strain that is otherwise more closely related to green strains. Our data show that the insertion of a DNA fragment containing the 19.7-kb phycoerythrin gene cluster has been facilitated by homologous recombination, also replacing a region of the phycocyanin operon. These findings demonstrate that large DNA fragments spanning entire functional gene clusters can be effectively transferred between closely related cyanobacterial strains and result in a changed phenotype. Further, the results shed new light on the discussion of the role of horizontal gene transfer in the sporadic distribution of large gene clusters in cyanobacteria, as well as the appearance of red and green strains.

Escherichia coli multilocus sequence type 38 from humans and broiler production represent distinct monophyletic groups.

Escherichia coli belonging to multilocus sequence type 38 (ST38) is a well-known cause of extra-intestinal infections in humans, and are frequently associated with resistance to extended-spectrum cephalosporins (ESCs). Resistance to carbapenems, mediated by blaOXA-genes has also been reported in this ST. Recently, the European Centre for Disease Prevention and Control (ECDC) released a rapid risk assessment on the increased detection of OXA-244 producing E. coli ST38 in humans, requesting further knowledge to determine the source. ST38 is also one of the most common STs among ESC-resistant E. coli from broiler production. Our aim was to investigate the genetic characteristics and relationship between E. coli ST38 from broiler production and humans, and to investigate if there has been a potential spillover between these sources. A total of 288 E. coli ST38 genomes isolated from humans in Europe (collected 2009-2019) and from Nordic broiler production (collected 2011-2014) were analyzed. The results showed distinct monophyletic clades associated to humans and broiler production. Furthermore, there were differences in the ESC resistance genes present in E. coli ST38 from the two sources. The blaOXA-244 gene was not present in E. coli from broiler production. Our results show that ST38 from humans and broiler production belong to well-separated clades, and suggest that the increased detection of OXA-244-producing E. coli ST38 in humans is not associated with spillover from broiler production.

Phylogeography and host specificity of Pasteurellaceae pathogenic to sea-farmed fish in the north-east Atlantic.

The present study was undertaken to address the recent spate of pasteurellosis outbreaks among sea-farmed Atlantic salmon (Salmo salar) in Norway and Scotland, coinciding with sporadic disease episodes in lumpfish (Cyclopterus lumpus) used for delousing purposes in salmon farms. Genome assemblies from 86 bacterial isolates cultured from diseased salmon or lumpfish confirmed them all as bona fide members of the Pasteurellaceae family, with phylogenetic reconstruction dividing them into two distinct branches sharing <88% average nucleotide identity. These branches therefore constitute two separate species, namely Pasteurella skyensis and the as-yet invalidly named "Pasteurella atlantica". Both species further stratify into multiple discrete genomovars (gv.) and/or lineages, each being nearly or fully exclusive to a particular host, geographic region, and/or time period. Pasteurellosis in lumpfish is, irrespective of spatiotemporal origin, linked almost exclusively to the highly conserved "P. atlantica gv. cyclopteri" (Pac). In contrast, pasteurellosis in Norwegian sea-farmed salmon, dominated since the late-1980s by "P. atlantica gv. salmonicida" (Pas), first saw three specific lineages (Pas-1, -2, and -3) causing separate, geographically restricted, and short-lived outbreaks, before a fourth (Pas-4) emerged recently and became more widely disseminated. A similar situation involving P. skyensis (Ps) has apparently been unfolding in Scottish salmon farming since the mid-1990s, where two historic (Ps-1 and -2) and one contemporary (Ps-3) lineages have been recorded. While the epidemiology underlying all these outbreaks/epizootics remains unclear, repeated detection of 16S rRNA gene amplicons very closely related to P. skyensis and "P. atlantica" from at least five cetacean species worldwide raises the question as to whether marine mammals may play a part, possibly as reservoirs. In fact, the close relationship between the studied isolates and Phocoenobacter uteri associated with harbor porpoise (Phocoena phocoena), and their relatively distant relationship with other members of the genus Pasteurella, suggests that both P. skyensis and "P. atlantica" should be moved to the genus Phocoenobacter.

Alterations of microbiota in urine from women with interstitial cystitis.

BACKGROUND: Interstitial Cystitis (IC) is a chronic inflammatory condition of the bladder with unknown etiology. The aim of this study was to characterize the microbial community present in the urine from IC female patients by 454 high throughput sequencing of the 16S variable regions V1V2 and V6. The taxonomical composition, richness and diversity of the IC microbiota were determined and compared to the microbial profile of asymptomatic healthy female (HF) urine. RESULTS: The composition and distribution of bacterial sequences differed between the urine microbiota of IC patients and HFs. Reduced sequence richness and diversity were found in IC patient urine, and a significant difference in the community structure of IC urine in relation to HF urine was observed. More than 90% of the IC sequence reads were identified as belonging to the bacterial genus Lactobacillus, a marked increase compared to 60% in HF urine. CONCLUSION: The 16S rDNA sequence data demonstrates a shift in the composition of the bacterial community in IC urine. The reduced microbial diversity and richness is accompanied by a higher abundance of the bacterial genus Lactobacillus, compared to HF urine. This study demonstrates that high throughput sequencing analysis of urine microbiota in IC patients is a powerful tool towards a better understanding of this enigmatic disease.

Draft Genome Sequences of 23 Tenacibaculum Isolates from Farmed Norwegian Lumpfish, Cyclopterus lumpus.

Draft genome sequences of 23 Tenacibaculum sp. strains that were isolated from Cyclopterus lumpus (lumpfish) were investigated to elucidate possible routes of transmission between Salmo salar (Atlantic salmon) and lumpfish.

Closed Genome Sequences of Providencia alcalifaciens Isolates from Dogs.

Eight Providencia alcalifaciens isolates from eight different dogs in Norway with acute hemorrhagic diarrhea were sequenced. Based on Illumina and Oxford Nanopore Technologies sequencing, all of the genomes were complete and closed after hybrid assembly.

Assessing diversity of the female urine microbiota by high throughput sequencing of 16S rDNA amplicons.

BACKGROUND: Urine within the urinary tract is commonly regarded as "sterile" in cultivation terms. Here, we present a comprehensive in-depth study of bacterial 16S rDNA sequences associated with urine from healthy females by means of culture-independent high-throughput sequencing techniques. RESULTS: Sequencing of the V1V2 and V6 regions of the 16S ribosomal RNA gene using the 454 GS FLX system was performed to characterize the possible bacterial composition in 8 culture-negative (<100,000 CFU/ml) healthy female urine specimens. Sequences were compared to 16S rRNA databases and showed significant diversity, with the predominant genera detected being Lactobacillus, Prevotella and Gardnerella. The bacterial profiles in the female urine samples studied were complex; considerable variation between individuals was observed and a common microbial signature was not evident. Notably, a significant amount of sequences belonging to bacteria with a known pathogenic potential was observed. The number of operational taxonomic units (OTUs) for individual samples varied substantially and was in the range of 20-500. CONCLUSIONS: Normal female urine displays a noticeable and variable bacterial 16S rDNA sequence richness, which includes fastidious and anaerobic bacteria previously shown to be associated with female urogenital pathology.

One Health Surveillance Codex: promoting the adoption of One Health solutions within and across European countries.

Cross-sector communication, collaboration and knowledge exchange are still significant challenges for practical adoption of the One Health paradigm. To address these needs the "One Health Surveillance Codex" (OHS Codex) was established to provide a framework for the One Health community to continuously share practical solutions (e.g. tools, technical resources, guidance documents and experiences) applicable for national and international stakeholders from different One Health Surveillance sectors. Currently, the OHS Codex provides a number of resources that support the adoption of the OH paradigm in areas linked to the harmonization and interpretation of surveillance data. The OHS Codex framework comprises four high-level "action" principles, which respectively support collaboration, knowledge exchange, data interoperability, and dissemination. These principles match well with priority areas identified in the "Tripartite Guide to Addressing Zoonotic Diseases in Countries" published by WHO, FAO and OIE. Within each of the four principles, the OHS Codex provides a collection of useful resources as well as pointers to success stories for the application of these resources. As the OHS Codex is designed as an open community framework, it will continuously evolve and adapt to the needs of the OH community in the future.

An Official Outbreak Investigation of Acute Haemorrhagic Diarrhoea in Dogs in Norway Points to Providencia alcalifaciens as a Likely Cause.

An outbreak investigation was initiated in September 2019, following a notification to the Norwegian Food Safety Authority (NFSA) of an unusually high number of dogs with acute haemorrhagic diarrhoea (AHD) in Oslo. Diagnostic testing by reporting veterinarians had not detected a cause. The official investigation sought to identify a possible common cause, the extent of the outbreak and prevent spread. Epidemiological data were collected through a survey to veterinarians and interviews with dog owners. Diagnostic investigations included necropsies and microbiological, parasitological and toxicological analysis of faecal samples and food. In total, 511 dogs with acute haemorrhagic diarrhoea were registered between 1 August and 1 October. Results indicated a common point source for affected dogs, but were inconclusive with regard to common exposures. A notable finding was that 134 of 325 faecal samples (41%) cultured positive for Providencia alcalifaciens. Whole genome sequencing (WGS) of 75 P. alcalifaciens isolates from 73 dogs revealed that strains from 51 dogs belonged to the same WGS clone. Findings point to P. alcalifaciens as implicated in the outbreak, but investigations are needed to reveal the pathogenic potential of P. alcalifaciens in dogs and its epidemiology.