Abstract
Corynebacterium jeikeium is an opportunistic pathogen which has been noted for significant genomic diversity. The population structure within this species remains poorly understood. Here, we explore the relationships among 15 clinical isolates of C. jeikeium (reference strains K411 and ATCC 43734, and 13 primary isolates collected over a period of 7 years) through genetic, genomic, and phenotypic studies. We report a high degree of divergence among strains based on 16S ribosomal RNA (rRNA) gene and rpoB gene sequence analysis, supporting the presence of genetically distinct subgroups. Whole genome sequencing indicates genomic-level dissimilarity among subgroups, which qualify as four separate and distinct Corynebacterium species based on an average nucleotide identity (ANIb) threshold of <95%. Functional distinctions in antibiotic susceptibilities and metabolic profiles characterize two of these genomospecies, allowing their differentiation from others through routine laboratory testing. The remaining genomospecies can be classified through a biphasic approach integrating phenotypic testing and rpoB gene sequencing. The genomospecies predominantly recovered from patient specimens does not include either of the existing C. jeikeium reference strains, implying that studies of this pathogen would benefit from examination of representatives from the primary disease-causing group. The clinically dominant genomospecies also has the smallest genome size and gene repertoire, suggesting the possibility of increased virulence relative to the other genomospecies. The ability to classify isolates to one of the four C. jeikeium genomospecies in a clinical context provides diagnostic information for tailoring antimicrobial therapy and may aid in identification of species-specific disease associations.