The pan-genome of the emerging multidrug-resistant pathogen Corynebacterium striatum.

Corynebacterium striatum, a common constituent of the human skin microbiome, is now considered an emerging multidrug-resistant pathogen of immunocompromised and chronically ill patients. However, little is known about the molecular mechanisms in the transition from colonization to the multidrug-resistant (MDR) invasive phenotype in clinical isolates. This study performed a comprehensive pan-genomic analysis of C. striatum, including isolates from "normal skin microbiome" and from MDR infections, to gain insights into genetic factors contributing to pathogenicity and multidrug resistance in this species. For this, three novel genome sequences were obtained from clinical isolates of C. striatum of patients from Brazil, and other 24 complete or draft C. striatum genomes were retrieved from GenBank, including the ATCC6940 isolate from the Human Microbiome Project. Analysis of C. striatum strains demonstrated the presence of an open pan-genome (α = 0.852803) containing 3816 gene families, including 15 antimicrobial resistance (AMR) genes and 32 putative virulence factors. The core and accessory genomes included 1297 and 1307 genes, respectively. The identified AMR genes are primarily associated with resistance to aminoglycosides and tetracyclines. Of these, 66.6% are present in genomic islands, and four AMR genes, including aac(6')-ib7, are located in a class 1-integron. In conclusion, our data indicated that C. striatum possesses genomic characteristics favorable to the invasive phenotype, with high genomic plasticity, a robust genetic arsenal for iron acquisition, and important virulence determinants and AMR genes present in mobile genetic elements.

Suboptimal bioinformatic predictions of antimicrobial resistance from whole-genome sequences in multidrug-resistant Corynebacterium isolates.

Herein, we combined different bioinformatics tools and databases (BV-BRC, ResFinder, RAST, and KmerResistance) to perform a prediction of antimicrobial resistance (AMR) in the genomic sequences of 107 Corynebacterium striatum isolates for which trustable antimicrobial susceptibility (AST) phenotypes could be retrieved. Then, the reliabilities of the AMR predictions were evaluated by different metrics: area under the ROC curve (AUC); Major Error Rates (MERs) and Very Major Error Rates (VMERs); Matthews Correlation Coefficient (MCC); F1-Score; and Accuracy. Out of 15 genes that were reliably detected in the C. striatum isolates, only tetW yielded predictive values for tetracycline resistance that were acceptable considering Food and Drug Administration (FDA)'s criteria for quality (MER < 3.0% and VMER with a 95% C.I. ≤1.5-≤7.5); this was accompanied by a MCC score higher than 0.9 for this gene. Noteworthy, our results indicate that other commonly used metrics (AUC, F1-score, and Accuracy) may render overoptimistic evaluations of AMR-prediction reliabilities on imbalanced datasets. Accordingly, out of 10 genes tested by PCR on additional multidrug-resistant Corynebacterium spp. isolates (n = 18), the tetW gene rendered the best agreement values with AST profiles (94.11%). Overall, our results indicate that genome-based AMR prediction can still be challenging for MDR clinical isolates of emerging Corynebacterium spp.

Pan-genomic analysis of Corynebacterium amycolatum gives insights into molecular mechanisms underpinning the transition to a pathogenic phenotype.

Corynebacterium amycolatum is a nonlipophilic coryneform which is increasingly being recognized as a relevant human and animal pathogen showing multidrug resistance to commonly used antibiotics. However, little is known about the molecular mechanisms involved in transition from colonization to the MDR invasive phenotype in clinical isolates. In this study, we performed a comprehensive pan-genomic analysis of C. amycolatum, including 26 isolates from different countries. We obtained the novel genome sequences of 8 of them, which are multidrug resistant clinical isolates from Spain and Tunisia. They were analyzed together with other 18 complete or draft C. amycolatum genomes retrieved from GenBank. The species C. amycolatum presented an open pan-genome (α = 0.854905), with 3,280 gene families, being 1,690 (51.52%) in the core genome, 1,121 related to accessory genes (34.17%), and 469 related to unique genes (14.29%). Although some classic corynebacterial virulence factors are absent in the species C. amycolatum, we did identify genes associated with immune evasion, toxin, and antiphagocytosis among the predicted putative virulence factors. Additionally, we found genomic evidence for extensive acquisition of antimicrobial resistance genes through genomic islands.