Whole Genome Sequencing and Comparative Genomics Analysis of Goat-Derived Klebsiella oxytoca

Background: This research aims to enhance the genomic database of Klebsiella oxytoca by identifying virulence genes through the whole genome sequencing and comparative analysis of a goat-derived K. oxytoca (KOHN1) strain, while clarifying the relationship between its genetic evolution and virulence, ultimately providing a theoretical foundation for clinical prevention and diagnosis. Methods: Third-generation Oxford Nanopore Technologies (ONT) sequencing and second-generation Illumina sequencing were used to sequence the strain and analyze the database annotations. Screening for 10 virulence genes was conducted using PCR. Comparative genomic analyses of the strain KOHN1 with four human-derived K. oxytoca model strains were performed using collinearity analysis, taxonomy classification through ANI analysis, and gene function family analysis. Results: The genome size of the KOHN1 strain was 5,817,806 bp, and the GC content was 55.14%. It contained 5227 predicted coding genes, including 25 rRNA genes, 85 tRNA genes, and 53 sRNA genes. A total of 14 type VI secretion system effector proteins and 146 virulence factor-related genes were annotated. Additionally, eight virulence genes-fimA, fimH, entB, mrkD, clpV, rmpA, vgrG, and hcp-were detected through PCR identification. The strain has 448 drug resistance genes, mainly against beta-lactams and fosfomycins. Comparative genomic analysis indicated that its closest relation is the human isolate ASM338647. Conclusions: In this study, the whole genome sequence of a goat-derived K. oxytoca (KOHN1) strain was obtained, revealing its evolutionary relationship with domestic and foreign isolates and providing a reference for future studies on the mechanisms of antimicrobial resistance and the pathogenicity of K. oxytoca.