Comparative Genomic Analysis of Pathogenic Factors of Pectobacterium Species Isolated in South Korea Using Whole-Genome Sequencing

Abstract

In this study, we conducted whole-genome sequencing with six species of Pectobacterium composed of seven strains, JR1.1, BP201601.1, JK2.1, HNP201719, MYP201603, PZ1, and HC, for the analysis of pathogenic factors associated with the genome of Pectobacterium. The genome sizes ranged from 4,724,337 bp to 5,208,618 bp, with the GC content ranging from 50.4% to 52.3%. The average nucleotide identity was 98% among the two Pectobacterium species and ranged from 88% to 96% among the remaining six species. A similar distribution was observed in the carbohydrate-active enzymes (CAZymes) class and extracellular plant cell wall degrading enzymes (PCWDEs). HC showed the highest number of enzymes in CAZymes and the lowest number in the extracellular PCWDEs. Six strains showed four subsets, and HC demonstrated three subsets, except hasDEF, in type I secretion system, while the type II secretion system of the seven strains was conserved. Components of human pathogens, such as Salmonella pathogenicity island 1 type type III secretion system (T3SS) and effectors, were identified in PZ1; T3SSa was not identified in HC. Two putative effectors, including hrpK, were identified in seven strains along with dspEF. We also identified 13 structural genes, six regulator genes, and five accessory genes in the type VI secretion system (T6SS) gene cluster of six Pectobacterium species, along with the loss of T6SS in PZ1. HC had two subsets, and JK2.1 had three subsets of T6SS. With the GxSxG motif, the phospholipase A gene did locate among tssID and duf4123 genes in the T6SSa cluster of all strains. Important domains were identified in the vgrG/paar islands, including duf4123, duf2235, vrr-nuc, and duf3396.