- Research
- Open access
- Published:
Genomic and phenotypic analysis of a novel clinical isolate of Corynebacterium pyruviciproducens
BMC Microbiology volume 23, Article number: 385 (2023)
Abstract
Background
Corynebacterium pyruviciproducens is a recently described species of Corynebacterium. There are few reports on the microbiological characteristics of the new species, and there is a lack of reports on the genomic analysis of the species.
Results
This study involved a clinical isolate from the pus of a hospital patient with sebaceous gland abscesses. The clinically isolated strain was identified as C. pyruviciproducens strain WYJY-01. In this study, referring to Koch’s postulates, we observed the pathological changes of animal models infected by intraperitoneal injection and subcutaneous injection of pure culture of the strain WYJY-01. Furthermore, the strain WYJY-01 was isolated and cultured again from animal models' subcutaneous abscess drainage fluid. Subsequently, the genomics of the strain WYJY-01 was analyzed. By comparing various gene databases, this study predicted the core secondary metabolite gene cluster of the strain WYJY-01, virulence factor genes carried by prophage, pathogenicity islands, and resistance islands. In addition, the genomes of C. pyruviciproducens strain WYJY-01, ATCC BAA-1742 T, and UMB0763 were analyzed by comparative genomics, and the differential genes of strain WYJY-01 were compared, and their functions were analyzed.
Conclusion
The findings showed that the strain WYJY-01 had pathogenicity, supplementing the phenotype characteristics of C. pyruviciproducens. Meanwhile, this research revealed the possible molecular mechanism of the pathogenicity of the strain WYJY-01 at the gene level through whole genome sequence analysis, providing a molecular basis for further research.
Background
Corynebacterium is a group of Gram-positive bacilli characterized by rod-like expansion at one or both ends, belonging to actinomycetes. Corynebacterium is mainly isolated from human, animal, soil, water, and other environmental resources [1, 2], most of which belong to normal flora, and some belong to conditional pathogens. Only a few bacteria in Corynebacterium are pathogenic to human beings [3], among which Corynebacterium diphtheria is the most specific pathogenic bacterium, and the primary pathogenic substance is diphtheria toxin, which is the expression product of lysogenic infection of β-bacteriophage [4]. Therefore, the prophage in the genetic material of Corynebacterium has a particular influence on the pathogenicity of Corynebacterium. In addition, Corynebacterium grows slowly in various media and needs some particular nutrients for its metabolic growth [5]. Among them, lipophilic Corynebacterium requires special culture conditions, which are difficult to isolate and culture in the routine clinical laboratory, and easy to miss diagnosis and misdiagnosis [6]. It is worth noting that many studies have shown that the main pathogenic bacteria of granulomatous mastitis are Corynebacterium [7]. However, the pathogenesis of the disease has not been studied. The researchers found preliminarily that most of the pathogenic Corynebacterium in this disease were lipophilic bacteria, and there were lipid bubbles around the bacteria in the abscess [8]. Meanwhile, C. pyruviciproducens has been found in the gut microbiomes of humans [9]. Therefore, the pathogenicity of lipophilic Corynebacterium has been widely discussed.
Corynebacterium pyruviciproducens is a recently described species of Corynebacterium. It was isolated from the inguinal abscess specimen collected by the American Medical Research Center in 2006 by Tong J, and was officially named in 2010 [10]. C. pyruviciproducens does not contain diphtheria toxin, tuberculous stearic acid, and other known virulence factors, and can be used as an immunomodulator to promote the host's humoral immune response to pathogenic microorganisms by regulating the function of macrophages [11, 12]. Then, D. Goldenberger et al. cultured urine, blood culture, skin infection, pus, joint cavity effusion, and other specimens from Switzerland and Canada and isolated C. pyruviciproducens. These clinical isolates of C. pyruviciproducens were found to differ from the type strain ATCC BAA-1742Â T in biochemical reaction and antimicrobial susceptibility [13]. As C. pyruviciproducens was discovered recently, there are relatively few literature reports related to this species, and there is a lack of research on its genomic analysis. Futhermore, up to present, there is an absence of studies describing the putative pathogenicity of this species. Therefore, the purpose of this research project is to supplement the microbiological characteristics of the pathogenicity of C. pyruviciproducens, and predict the molecular biological basis related to its microbiological characteristics through whole-genome sequence analysis.
Results
General microbiological characterization of C. pyruviciproducens strain WYJY-01
C. pyruviciproducens strain WYJY-01 isolated and cultured from abscess drainage fluid was cultured on Columbia blood plates at 37 °C for 48 h, and its colony morphology was round and smooth, and the colony center was convex and transparent (Fig. 1A). The strain WYJY-01 was Gram-positive under the light microscope after Gram staining, and its shape was rod-shaped (Fig. 1B). However, the strain WYJY-01 was identified as Kocuria rosea by VITEK 2 Compact full-automatic microbial analyzer. Because Kocuria rosea is the Gram-positive cocci, it was evident that the results were not consistent with each other [14].
The comparison results based on the GenBank database showed that the 16S rRNA gene sequence of the strain WYJY-01 was the closest to Corynebacterium sp. 718260/2011 (GenBank JX501770.1), with a sequence similarity of 100%. The rpoB gene sequence of the strain WYJY-01 was the closest to that of Corynebacterium sp. NML00-0179 (GenBank GU304660.1), with a sequence similarity of 99.75%. These two strains were clinical isolates of C. pyruviciproducens. Subsequently, a phylogenetic analysis of the 16S rRNA gene with a tree constructed based on the Neighbor-Joining method showed that the strain WYJY-01 was a member of the genus Corynebacterium (Fig. 2).
The difference in biochemical characteristics between the strain WYJY-01 and the strain ATCC BAA-1742Â T was (Additional file 1: Table S1): CAMP reaction and reduction of nitrates of strain WYJY-01 are positive; fermentation of D-ribose, D-xylose, D-glucose, maltose and sucrose of strain WYJY-01 were negative.
According to the Clinical and Laboratory Standards Institute (CLSI) guidelines [15], the drug sensitivity of the strain WYJY-01 was measured by the Etest method (Table 1). The strain WYJY-01 was resistant to ceftriaxone, gentamicin, erythromycin, ciprofloxacin, and clindamycin. However, the strain ATCC BAA-1742Â T was sensitive to gentamicin and ciprofloxacin.
C. pyruviciproducens strain WYJY-01 was cultured at 37 °C in nutrient broth (NB) medium and NB medium with 1% Tween 80 added. The OD600 of bacterial liquid in different media was measured at intervals to obtain the growth curve of the strain WYJY-01 (Fig. 1C). According to the growth curve, the strain WYJY-01 belongs to the lipophilic microorganism. The strain WYJY-01 was incubated at 37 °C for about 40 h in NB medium supplemented with 1% Tween 80 at a plateau stage.
Pathogenic phenotype of C. pyruviciproducens strain WYJY-01
The macroscopic changes in the organs of BALB/c mice in each group were observed after being dissected. Compared with BALB/c mice injected intraperitoneally with normal saline (Fig. 3A), BALB/c mice injected intraperitoneally with the strain WYJY-01 showed scattered small abscesses on the surface of the liver (Fig. 3B). Compared with BALB/c mice injected with normal saline subcutaneously (Fig. 3C), BALB/c mice injected with the strain WYJY-01 subcutaneously had abscess formation at the injection site (Fig. 3D).
The liver, spleen, and skin of normal and diseased BALB/c mice were taken and paraffin sectioned. After hematoxylin–eosin (HE) staining, the pathological changes in the tissue structure were observed under the light microscope (Fig. 4). The liver portal area of BALB/c mice injected intraperitoneally with the strain WYJY-01 showed cysts of different sizes and irregular shapes, suggesting cyst formation in the portal area (Fig. 4D) and partial central vein dilatation (Fig. 4G); Lymphocyte and neutrophil infiltration can be seen in the liver (Fig. 4H). The spleen of BALB/c mice injected intraperitoneally with the strain WYJY-01 showed reactive proliferation of splenic corpuscles (Fig. 4E). The subcutaneous tissue of BALB/c mice injected subcutaneously with the strain WYJY-01 showed inflammatory infiltration (Fig. 4F), inflammatory necrosis, cellulose-like exudation, and granulation tissue hyperplasia (Fig. 4I).
After the subcutaneous abscess drainage fluid of BALB/c mice that were subcutaneously injected with the strain WYJY-01 was smeared directly, heated, and fixed. Gram-positive and rod-shaped bacteria were seen after Gram staining. The abscess drainage fluid was cultured on Columbia blood plate at 37 °C for 48 h. The colonies isolated and cultured are round, smooth, central bulge, and transparent (Additional file 2: Fig. S1A). Some of the isolated colonies in the blood plate were Gram-stained. Gram-positive and rod-shaped bacteria (Additional file 2: Fig. S1B) were observed under the light microscope, and their morphology was the same as that of the strain WYJY-01. Moreover, the 16S rRNA gene sequence of the isolated colonies from the subcutaneous abscess drainage fluid is consistent with the 16S rRNA gene sequence of the strain WYJY-01. This proves that the colonies recovered as C. pyruviciproducens strain WYJY-01.
Genomic characteristics of C. pyruviciproducens strain WYJY-01
The genome of the strain WYJY-01 was extracted and electrophoresed through agarose gel (Additional file 3: Fig. S2A). Based on the Oxford Nanopore technology platform, the raw sequencing results of the strain WYJY-01 genome showed that the Reads N50 length was 18,025Â bp (Additional file 3: Fig. S2B), and the MeanQual was 10.7. Therefore, the sequencing result of the strain WYJY-01 genome was high quality.
The assembled contig sequence was aligned with the Nucleotide database to determine the chromosome type. The assembly generated a single circular chromosome with total size of 2,881,050Â bp, and the GC content was 60.76%. The GC content of the strain ATCC BAA-1742Â T genome was 61.1% (BioProject ID PRJNA78965). The GC content of the strain UMB0763 genome was 60.8% (BioProject ID PRJNA316969). It can be seen that the difference in GC content among the three C. pyruviciproducens species was less than 1%. The results based on the GC content also supplemented the identification of the species of the strain WYJY-01.
The genome of the strain WYJY-01 was compared with the most similar sequence in the Nr database by BLAST. Thus, the species distribution of the matched sequence was calculated (Additional file 4: Fig. S3A). The statistical chart showed that 85.71% of the sequences in the genome of the strain WYJY-01 were most similar to C. pyruviciproducens, and 4.34% of the sequences were similar to C.glucuronolyticum. Based on the genome components predicted by the database and the gene COG function classification predicted by the eggNOG database (Additional file 4: Fig. S3B), the location of annotated loci between the genome components of the strain WYJY-01 was presented by Circos v0.66 (Fig. 5).
The ANI value based on the strain WYJY-01 and ATCC BAA-1742Â T was 97.44%. And the ANI value based on the strain WYJY-01 and UMB0763 was 98.56%. The dDDH value based on the strain WYJY-01 and ATCC BAA-1742Â T was 76.70%. And the dDDH value based on the strain WYJY-01 and UMB0763 was 87.00%. Venn diagram of gene families of the strain WYJY-01, ATCC BAA-1742Â T, and UMB0763 genome showed (Additional file 5: Fig. S4) that there were 19 unique gene families in the genome of the strain WYJY-01, and 128 gene families different from strain ATCC BAA-1742Â T genome. Furthermore, the number of unique genes of the strain WYJY-01 was 314, and the number of differential genes different from strain ATCC BAA-1742Â T was 441.
Virulence factors and antimicrobial resistance genes
By searching the Comprehensive Antibiotic Resistance Database (CARD) using BLSAT, the predicted resistance genes in the genome of the strain WYJY-01 were ermX, aph(3')-Ia,aph(3″)-Ib,aph(6)-Id, cmx. ErmX is an rRNA methyltransferase that protected the ribosome from inactivation due to antibiotic binding. Therefore, ErmX made bacteria resistant to macrolide antibiotics, lincosamide antibiotics, and streptogramin antibiotics by changing the target of antibiotics. APH(3')-Ia, APH(3″)-Ib, APH(6)-Id were aminoglycoside phosphotransferases, which can enable bacteria to acquire resistance against aminoglycoside antibiotics by inactivating antibiotics. Cmx was a chloramphenicol exporter, which can make bacteria resistant to phenicol antibiotics and chloramphenicol by exuding antibiotics.
By searching the Virulence Factor Database (VFDB) using BLAST, it was predicted that the genes in the genome of the strain WYJY-01 might be involved in encoding the virulence factors, including capsule, intercellular adhesion proteins, hemolysin, type III secretion system (TTSS), type IV secretion system, type VI secretion system, type VII secretion system, rhamnolipid, mucoid exopolysaccharide, SpaA-type pilus, polar flagella, biofilm, adhesin, metal ions transporter (mainly iron ions), P60 extracellular protein, etc.
Horizontal gene transfer is the exchange of genetic information between different species of bacteria through the mobile components. Mobile components can carry virulence factor genes and resistance genes, including genomic islands, prophages, insertion sequences, integrons, transposons, plasmids, etc. [16, 17].
The prophage region of the genome of the strain WYJY-01 predicted by PhiSpy v2.3 was 69,911 bp in length (starting at 814,919 bp and ending at 884,829 bp). Combined with the VFDB prediction results, the virulence factor genes in the prophage region of the genome include farB, clpP, hspX, and flpF (Fig. 6A). Among them, FarB was a fatty acid efflux system protein, and it mediated the resistance to antimicrobial long-chain fatty acids. HspX stabilized cell structure during long-term survival and permitted the bacilli to survive within the low-oxygen environment of the granuloma.
There were 13 genomic islands predicted by IslandPath-DIMOB v0.1 in the genome of the strain WYJY-01. Genomic island 4 (Fig. 6B) started at 613,513 bp and ended at 647,861 bp. The virulence factor genes in this genomic island include irp2, fbpC, irtB, irtA, pchC, rhlB, cap8J, and cylA. Meanwhile, these genes were the different genes that were different from the genome of the strain ATCC BAA-1742 T. Genomic island 6 (Fig. 6C) started at 885,274 bp and ended at 928,269 bp. The genomic island contained virulence factor genes, including irtA, iroC, phoR, phoP, and aha1389. It was worth noting that the five resistance genes in the genome of the strain WYJY-01 searched by the CARD were all located on Genomic island 6. Genomic island 8 (Fig. 6D) started at 1,494,249 bp and ended at 1,505,409 bp, and Genomic island 10 (Fig. 6E) started at 1,916,529 bp and ended at 1,928,888 bp. The virulence factor genes contained in these two genomic islands included bopD, fbpC, wzt, cpsE, and ugd, and these genes were also differential genes from the genome of the strain ATCC BAA-1742 T.
The results of gene function annotated by the KEGG database (Fig. 7) showed the functions of the differential genes between the strain WYJY-01 and the strain ATCC BAA-1742 T were mainly clustered in environmental information processing (Fig. 7B). The differential genes clustered into ABC transporter and predicted by VFDB included fbpC on Genomic island 4, fbpC on Genomic island 8, and wzt on Genomic island 10, shuU, and fepC. One of the differential genes clustered into the two-component system was the phoR of Genomic island 6. The results of gene function annotated by the GO database showed (Additional file 6: Fig. S5) that the proportion of differential genes in the two categories of the extracellular region and electron carrier activity in each functional category was significantly higher than that of the genome of the strain WYJY-01.
Secondary metabolite gene clusters
The secondary metabolite gene clusters of the strain WYJY-01 predicted by antiSMASH v5.0.0 consists of two regions. The non-ribosomal polypeptides (NRPS) gene cluster (Fig. 8A) in the genome started at 585,737 bp and ended at 641,084 bp, and the secondary metabolites were roughly predicted as (cys) + (cys—cys), peptide bonds linked three Cysteines. The polyketides (T1PKS) gene cluster (Fig. 8B) in the genome started at 1,167,876 bp and ended at 1,212,552 bp, and the secondary metabolite was roughly predicted to be pyruvic acid (CH3COCOOH). It was worth noting that the core structure region of the NRPS gene cluster had overlapping genes with Genomic island 4 (Fig. 6B), and most of the genes in the core structure region were the differential genes from the genomic of the strain ATCC BAA-1742 T.
Discussion
C. pyruviciproducens is a lipogenic Corynebacterium capable of producing pyruvate, which was named C. pyruviciproducens strain ATCC BAA-1742 T in 2006 by Tong J [10]. A clinical isolate of C. pyruviciproducens involved in this reaearch was named C. pyruviciproducens strain WYJY-01. The colony morphology of the two C. pyruviciproducens strains was consistent, and the morphology was consistent under the microscope after Gram staining, and both were lipophilic bacteria. In addition, the GC content and ANI and DDH analysis of the genomes of both strains indicated that they were of the same species [18]. However, there were some differences in other general microbiological phenotypes of the two strains, such as bacterial biochemical reactions, antimicrobial susceptibility tests, etc. It is worth noting that, in our research on animal models with different infection routes, we found that under the guidance of Koch’s postulates [19,20,21], the strain WYJY-01 may be pathogenic. The strain WYJY-01 can infect BALB/c mice through intraperitoneal and subcutaneous routes, and cause significant organic changes in BALB/c mice. As C. pyruviciproducens has been discovered for a short time and there are few relevant studies, these data are of great significance for public health monitoring and provide a reference for the phenotypic identification of this bacterium in future clinical work.
At present, studies on the genomic characteristics of C. pyruviciproducens are scarce. We extracted the genomic data of C. pyruviciproducens strain ATCC BAA-1742Â T and UMB0763 from the NCBI Genome database. These two C. pyruviciproducens strains are clinical isolates. However, the microbiological phenotype description of strain UMB0763 is not clear [22]. Therefore, we focus on the results of the genomic analysis of C. pyruviciproducens strain WYJY-01 and ATCC BAA-1742Â T.
Horizontal gene transfer can lead to DNA exchange between different species of bacteria, thereby enhancing their adaptability to the environment, and is one of the important pathways for bacterial evolution. Genomic islands are very active carriers in horizontal gene transfer, and play a crucial role in the plasticity and evolutionary evolution of bacteria [23]. Genomic islands can carry virulence factor genes, resistance genes, metabolic regulation genes, etc., which can lead to the production of pathogenic multi-drug resistant bacteria [24, 25].
Firstly, through CARD resistance gene prediction and comparative genomics analysis, it was concluded that the resistance phenotype and genotype of the strain WYJY-01 against aminoglycoside antibiotic, macrolide antibiotic, and lincosamide antibiotic are consistent. However, the genotype associated with the drug resistance phenotype of the strain WYJY-01 against quinolones antibiotic was not predicted. Moreover, the genotype of the strain WYJY-01 resistant to chloramphenicol was predicted. Combined with the genomic islands prediction results of the strain WYJY-01 genome, it can be seen that five resistance genes of the strain WYJY-01 were located in Genomic island 6 (Fig. 6C). Therefore, Genomic island 6 is the resistance island of the strain WYJY-01, and it may play a decisive role in the process of bacterial resistance.
Secondly, pathogenic bacteria have gone through many processes from adsorption into host cells to successful infection, such as adsorption, colonization, escape from host immune response, local disease, and other stages. According to the prediction results of virulence factors of VFDB, the genome of the strain WYJY-01 contains genes related to virulence factors at various stages of bacterial infection of the host. Combined with the genomic islands prediction results of the strain WYJY-01 genome, it can be seen that the virulence factor genes in Genomic island 4 (Fig. 6B), such as irp2 [26], fbpC [27], irtA, irtB [28], pchC [29, 30], were mainly involved in the biological process of bacterial iron ion transport. This process plays an essential role in regulating the local pathogenic effect stage of bacteria, because the ability of pathogenic bacteria to compete with the host for iron ions plays an important role in bacterial survival and virulence [31]. Virulence factor genes in Genomic island 10 (Fig. 6E), such as cpsE [32] and ugd [33], were mainly involved in the formation process of the capsule, assisting in evading the host immune system by protecting bacteria from opsonophagocytosis and serum killing. Finally, it was predicted through the KEGG database that the genome of the strain WYJY-01 and the differential genes of the strain WYJY-01 and ATCC BAA-1742 T were mainly clustered into the ABC transporters module of environmental information processing. Further research found that the differential genes clustered to ABC transporters were mainly involved in the transport of iron ions [28], such as fbpC, wzt, shuU, fepC, etc. This further proved that the strain WYJY-01 showed a possible molecular mechanism at the gene level of the pathogenic phenotype.
Furthermore, we also predicted the gene cluster of the secondary metabolite of the strain WYJY-01. Through the analysis of the T1PKS database, C. pyruviciproducens had the core gene cluster of pyruvate (Fig. 8B), and the gene sequence of the T1PKS gene cluster also appeared in the genome of the strain ATCC BAA-1742 T. The core biological gene of the T1PKS gene cluster was also annotated as wcbR by VFDB, encoding capsular polysaccharide biosynthesis fatty acid synthase, which was the critical virulence determinant [34]. Through the NRPS database, we predicted that the strain WYJY-01 might metabolize into a compound consisting of three cysteines linked by a peptide bond. At the same time, the NRPS gene cluster overlaps with Genomic island 4 (Fig. 6B). This showed that Genomic island 4 was not only the core pathogenicity island of the strain WYJY-01, but also the metabolism island of the strain WYJY-01.
Conclusion
In general, this study supplemented the microbiological phenotypic characteristics of C. pyruviciproducens, and especially explored its pathogenic phenotype. Through whole-genome sequence analysis, this study revealed the possible molecular mechanism of the pathogenicity of this clinical isolate at the gene level, and also provided a basis for further research on the molecular basis of the relevant microbiological phenotypic characteristics of this species. This study expanded the comprehension for identifying, analyzing, and treating C. pyruviciproducens in clinical.
Methods
Identification of clinical isolates and general microbiological phenotypes
The clinical isolates of C. pyruviciproducens involved in this study were from a female patient with sebaceous gland cysts on the skin near the sternum in the Affiliated Hospital of Weifang Medical University. The Medical Ethics Committee of Weifang Medical University approved the study (No. 2022YX65) and waived the informed consent requirement because this study used only bacterial isolate and did not have any negative impact on patients. After the abscess was cut and drained, the clinical isolates were obtained, and the pus samples were isolated and cultured in our hospital. First, the isolate was identified and analyzed by VITEK 2 Compact full-automatic microbial analyzer [35]. Then the isolated bacteria were observed under the light microscope through Gram staining. Due to the inconsistency between the theoretical bacterial morphology predicted by the analyzer and the observed bacterial morphology under the light microscope after Gram staining, the isolate was identified by 16S rRNA and rpoB gene sequencing, and was named C. pyruviciproducens strain WYJY-01 [36, 37]. The 16S rRNA gene sequence of C. pyruviciproducens strain WYJY-01 was uploaded to the GenBank database. Preliminary sequence analysis of the 16S rRNA gene was conducted using the NCBI database, and strains with high homology in the NCBI database were selected for phylogenetic tree analysis. The phylogenetic tree was constructed by using Mega10.0 and the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches.
Next, the drug resistance phenotype of the isolate was determined by the Etest, a reference method for testing antimicrobial susceptibility to diverse drugs in accordance with CLSI guidelines [15]. In addition, the isolate was cultured in NB medium supplemented with 1% Tween 80, which further verified its lipid requirement and measured its growth curve [38, 39].
Study on pathogenic phenotypes based on Koch’s postulates
All experiments involving mice were approved by the Institutional Animal Care and Use Committee of Weifang Medical University (No. 2022SDL162). SPF grade experimental mice were BALB/c, male, weighing about 25–30 g, purchased from Qingdao Da Ren Fucheng Livestock Co., Ltd. C. pyruviciproducens strain WYJY-01 was resuscitated in Columbia blood agar medium. The colony was scraped and resuspended in normal saline, and adjusted to about 3 × 109 CFU/mL. Twenty-four BALB/c mice were randomly divided into four groups. In the first group, 6 mice were intraperitoneally injected with 1 mL of bacterial fluid. In the second group, 6 mice were injected with 100 μL of bacterial fluid subcutaneously. The third group was the control group. 6 mice were intraperitoneally injected with 1 mL sterile normal saline. The fourth group was the control group. 6 mice were subcutaneously injected with 100 μL sterile normal saline. After 3 days of injection, 6 mice in each experimental group were dissected, and taken out the relevant diseased organs to make ordinary paraffin sections. After HE staining, the changes in tissue structure were observed. The subcutaneous abscess drainage fluid of the infected model was isolated and cultured again. The isolated strains were identified.
Whole genome sequencing and assembly
C. pyruviciproducens strain WYJY-01 stored at -80 °C was resuscitated on Columbia blood agar medium and grew at 37 °C for 48 h. The fresh isolates were resuspended in molecular biology-grade water without nuclease. According to instructions, the collected C. pyruviciproducens strain WYJY-01 was extracted with a QIAGEN Genomic-tip DNA kit. The harvested DNA was electrophoresed through agarose gel and quantified by Qubit assays (Invitrogen). The concentration of the extracted genome was 70.2 ng/μL. OD260/280 of the gemome was 1.87, and OD260/230 of the gemome was 2.04. Whole genome sequencing was done by Oxford Nanopore Technologies DNA sequencing platform (Biomarker Technologies, Beijing, China). The reads after the quality control were assembled by Canu v1.5 [40]. Finally, we polished the assemblies using Pilon v1.22 [41]. The assembled C. pyruviciproducens strain WYJY-01 genome was uploaded to the NCBI database.
Genomic analysis
Prodigal v2.6.3 [42] was used to predict coding sequences (CDSs) in the genome of the strain. RepeatMasker v4.0.5 [43] was used to search for repetitive sequences. Rfam v12.0 [44], and tRNAscan-SE v1.3.1 [45] were used to predict tRNA, and other ncRNAs. And genome information obtained by assembly and prediction, such as information on tRNAs, rRNAs, repeat sequences, GC contents, and gene functions, was used to draw a circular genome map with the software Circos v0.66 [46].
The genomic islands and the prophage were predicted by IslandPath-DIMOB v0.1 [47] and PhiSpy v2.3 [48]. The gene sequence was compared with eggNOG [49], KEGG [50], Nr, and other functional databases by BLAST v2.2.29 to provide the functional annotation. Based on the comparison results of the Nr database, the gene sequence was annotated by Blast2GO v2.5 [51] for the functions of the GO database. The protein sequences of the strain WYJY-01 were blasted by the CARD [52] and VFDB [53], and the corresponding annotation results were obtained. The secondary metabolite gene cluster was predicted by antiSMASH v5.0.0 [54].
The Average Nucleotide Identity (ANI) and DNA–DNA hybridization (dDDH) of the genomic sequences of C. pyruviciproducens strain WYJY-01, ATCC BAA-1742 T, and UMB0763 were analyzed by fastANIv1.33 [55] and GGDC v3.0 [56]. The protein sequences of C. pyruviciproducens strain WYJY-01, ATCC BAA-1742 T, and UMB0763 were classified by the OrthoMCL [57]. Then, gene families were analyzed, including the specific gene family of the strain WYJY-01, the gene family shared by the three strains, and the gene family with a single copy of each strain. The differential genes of the strain WYJY-01 between the strain ATCC BAA-1742 T were annotated by KEGG metabolic pathway and GO annotation.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files. The GenBank accession numbers for the whole-genome and 16S rRNA gene sequences of C. pyruviciproducens strain WYJY-01 are SRX14997086 and OP799734.1, respectively.
References
Dover LG, Thompson AR, Sutcliffe IC, Sangal V. Phylogenomic reappraisal of fatty acid biosynthesis, mycolic acid biosynthesis and clinical relevance among members of the genus Corynebacterium. Front Microbiol. 2021;12:802532.
Burkovski A. Corynebacterium pseudodiphtheriticum : Putative probiotic, opportunistic infector, emerging pathogen. Virulence. 2015;6:673–4.
Coyle MB, Lipsky BA. Coryneform bacteria in infectious diseases: clinical and laboratory aspects. Clin Microbiol Rev. 1990;3:20.
Pappenheimer AM, Gill DM. Diphtheria: recent studies have clarified the molecular mechanisms involved in its pathogenesis. Science. 1973;182:353–8.
Dias FF, Bhat JV. Nutritional properties of Corynebacterium laevaniformans. Antonie Van Leeuwenhoek. 1964;30:176–84.
Sabbe LJM, Van De Merwe D, Schouls L, Bergmans A, Vaneechoutte M, Vandamme P. Clinical Spectrum of Infections Due to the Newly Described Actinomyces Species A. turicensis, A. radingae, and A. europaeus. J Clin Microbiol. 1999;37:8–13.
Fernandez LV, Fortuny AS, Rodriguez EF. Mastitis por Corynebacterium pyruviciproducens y Corynebacterium amycolatum en mujeres inmunocompetentes no lactantes. Rev Argent Microbiol. 2021;53:39–42.
Mathelin C, Riegel P, Chenard M-P, Brettes J-P. Association of Corynebacteria with granulomatous mastitis. Eur J Obstet Gynecol Reprod Biol. 2005;119:260–1.
Asuncion P, Liu C, Castro R, Yon V, Rosas M, Hooshmand S, et al. The effects of fresh mango consumption on gut health and microbiome – Randomized controlled trial. Food Sci Nutr. 2023;11:2069–78.
Tong J, Liu C, Summanen PH, Xu H, Finegold SM. Corynebacterium pyruviciproducens sp. nov., a pyruvic acid producer. Int J Syst Evol Microbiol. 2010;60:1135–40.
Zhang Y, Li X-X, Ma Y, Xu J, Zhao L-N, Qian X-F, et al. CP and CP-PGN protect mice against MRSA infection by inducing M1 macrophages. Sci Rep-uk. 2017;7:16877.
Tong J, Han Q, Wang S, Su Z, Zheng D, Shen P, et al. Corynebacterium pyruviciproducens, as an immune modulator, can promote the activity of macrophages and up-regulate antibody response to particulate antigen. Exp Biol Med (Maywood). 2012;237:1322–30.
Goldenberger D, Hinić V, Turan S, Schultheiss E, Pacheco AL, Frei R, et al. Extended characterization of Corynebacterium pyruviciproducens based on clinical strains from Canada and Switzerland. J Clin Microbiol. 2014;52:3180–3.
Stackebrandt E, Koch C, Gvozdiak O, Schumann P. Taxonomic Dissection of the Genus Micrococcus: Kocuria gen nov, Nesterenkonia gen. nov., Kytococcus gen. nov., Dermacoccus gen. nov., and Micrococcus Cohn 1872 gen. emend. Int J Syst Bacteriol. 1995;45:682–92.
CLSI. Methods for antimicrobial dilution and disk susceptibility testing of infrequently isolated or fastidious bacteria; M45-A2, 2nd ed. Wayne: Clinical and Laboratory Standards Institute; 2010.
Liu Y-Y, Wang Y, Walsh TR, Yi L-X, Zhang R, Spencer J, et al. Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect Dis. 2016;16:161–8.
Kumarasamy KK, Toleman MA, Walsh TR, Bagaria J, Butt F, Balakrishnan R, et al. Emergence of a new antibiotic resistance mechanism in India, Pakistan, and the UK: a molecular, biological, and epidemiological study. Lancet Infect Dis. 2010;10:597–602.
Kujur RRA, Das SK. Pseudomonas phenolilytica sp. nov., a novel phenol-degrading bacterium. Arch Microbiol. 2022;204:320.
Walker L, LeVine H, Jucker M. Koch’s postulates and infectious proteins. Acta Neuropathol. 2006;112:1–4.
Ramakrishnan L. Mycobacterium tuberculosis pathogenicity viewed through the lens of molecular Koch’s postulates. Curr Opin Microbiol. 2020;54:103–10.
Hosainzadegan H, Khalilov R, Gholizadeh P. The necessity to revise Koch’s postulates and its application to infectious and non-infectious diseases: a mini-review. Eur J Clin Microbiol Infect Dis. 2020;39:215–8.
Thomas-White KJ, Kliethermes S, Rickey L, Lukacz ES, Richter HE, Moalli P, et al. Evaluation of the urinary microbiota of women with uncomplicated stress urinary incontinence. Am J Obstet Gynecol. 2017;216:55.e1-55.e16.
Juhas M, Van Der Meer JR, Gaillard M, Harding RM, Hood DW, Crook DW. Genomic islands: tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol Rev. 2009;33:376–93.
Bardaji L, EcheverrÃa M, RodrÃguez-Palenzuela P, MartÃnez-GarcÃa PM, Murillo J. Four genes essential for recombination define GInts, a new type of mobile genomic island widespread in bacteria. Sci Rep-uk. 2017;7:46254.
Coetzee JN, Datta N, Hedges RW. R Factors from Proteus rettgeri. J Gen Microbiol. 1972;72:543–52.
Santos MCFD, Anderson CP, Neschen S, Zumbrennen-Bullough KB, Romney SJ, Kahle-Stephan M, et al. Irp2 regulates insulin production through iron-mediated Cdkal1-catalyzed tRNA modification. Nat Commun. 2020;11:296.
Li X, He J, Fu W, Cao P, Zhang S, Jiang T. Effect of Mycobacterium tuberculosis Rv3717 on cell division and cell adhesion. Microb Pathog. 2018;117:184–90.
Rodriguez GM, Smith I. Identification of an ABC transporter required for iron acquisition and virulence in Mycobacterium tuberculosis. J Bacteriol. 2006;188:424–30.
Ronnebaum TA, Lamb AL. Nonribosomal peptides for iron acquisition: pyochelin biosynthesis as a case study. Curr Opin Struc Biol. 2018;53:1–11.
Reimmann C, Patel HM, Walsh CT, Haas D. PchC Thioesterase optimizes nonribosomal biosynthesis of the peptide siderophore pyochelin in Pseudomonas aeruginosa. J Bacteriol. 2004;186:6367–73.
Teng T, Xi B, Chen K, Pan L, Xie J, Xu P. Comparative transcriptomic and proteomic analyses reveal upregulated expression of virulence and iron transport factors of Aeromonas hydrophila under iron limitation. Bmc Microbiol. 2018;18:52.
Cartee RT, Forsee WT, Bender MH, Ambrose KD, Yother J. CpsE from type 2 Streptococcus pneumoniae catalyzes the reversible addition of glucose-1-phosphate to a polyprenyl phosphate acceptor, initiating type 2 capsule repeat unit formation. J Bacteriol. 2005;187:7425–33.
Mouslim C, Groisman EA. Control of the Salmonella ugd gene by three two-component regulatory systems: Complex regulation by two-component systems. Mol Microbiol. 2003;47:335–44.
Price EP, Sarovich DS, Mayo M, Tuanyok A, Drees KP, Kaestli M, et al. Within-Host Evolution of Burkholderia pseudomallei over a Twelve-Year Chronic Carriage Infection. mBio. 2013;4:e00388-13.
Robinson A, McCarter YS, Tetreault J. Comparison of crystal enteric/nonfermenter system, API 20E system, and Vitek autoMicrobic system for identification of gram-negative bacilli. J Clin Microbiol. 1995;33:364–70.
Khamis A, Raoult D, La Scola B. rpoB Gene sequencing for Identification of Corynebacterium Species. J Clin Microbiol. 2004;42:3925–31.
Janda JM, Abbott SL. 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: pluses, perils, and pitfalls. J Clin Microbiol. 2007;45:2761–4.
Kropinski AM. Practical Advice on the One-Step Growth Curve. In: Clokie MRJ, Kropinski AM, Lavigne R, editors. Bacteriophages. Springer, New York: New York, NY; 2018. p. 41–7.
Govantes F. Serial dilution-based growth curves and growth curve synchronization for high-resolution time series of bacterial biofilm growth. In: Medina C, López-Baena FJ, editors. Host-pathogen interactions. Springer, New York: New York, NY; 2018. p. 159–69.
Koren S, Walenz BP, Berlin K, Miller JR, Bergman NH, Phillippy AM. Canu: scalable and accurate long-read assembly via adaptive k -mer weighting and repeat separation. Genome Res. 2017;27:722–36.
Walker BJ, Abeel T, Shea T, Priest M, Abouelliel A, Sakthikumar S, et al. Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement. PLoS One. 2014;9:e112963.
Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.
Tarailo-Graovac M, Chen N. Using repeat masker to identify repetitive elements in genomic sequences. Curr Protoc Bioinformatics. 2009;4:4.10.1-14.
Nawrocki EP, Burge SW, Bateman A, Daub J, Eberhardt RY, Eddy SR, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 2015;43:D130-7.
Chan PP, Lowe TM. tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences. In: Kollmar M, editor. Gene Prediction. Springer, New York: New York, NY; 2019. p. 1–14.
Krzywinski M, Schein J, Birol İ, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.
Bertelli C, Brinkman FSL. Improved genomic island predictions with IslandPath-DIMOB. Bioinformatics. 2018;34:2161–7.
Akhter S, Aziz RK, Edwards RA. PhiSpy: a novel algorithm for finding prophages in bacterial genomes that combines similarity- and composition-based strategies. Nucleic Acids Res. 2012;40:e126–e126.
Powell S, Forslund K, Szklarczyk D, Trachana K, Roth A, Huerta-Cepas J, et al. eggNOG v4.0: nested orthology inference across 3686 organisms. Nucl Acids Res. 2014;42:D231-9.
Kanehisa M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32:277D-280.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45:D566–73.
Chen L. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2004;33 Database issue:D325-8.
Blin K, Shaw S, Steinke K, Villebro R, Ziemert N, Lee SY, et al. antiSMASH 5.0: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019;47:W81-7.
Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun. 2018;9:5114.
Meier-Kolthoff JP, Carbasse JS, Peinado-Olarte RL, Göker M. TYGS and LPSN: a database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic Acids Res. 2022;50:D801–7.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for Eukaryotic genomes. Genome Res. 2003;13:2178–89.
Acknowledgements
Not applicable.
Funding
This work was supported in part by awards to the National Natural Science Foundation of China (No. 81802054); the support program for Youth Innovation Technology in Colleges and Universities of Shandong Province of China (Grant No. 2021KJ106); the Tai-Shan Scholar Program from Shandong Province (Grant No. tsqn202103116); the Shandong Provincial Natural Science Foundation (Grant No. ZR2022MH061); the Public Sponsored Domestic Visiting Program of Weifang Medical University (No.20237-11).
Author information
Authors and Affiliations
Contributions
QL, YYT, and JQW conceived and designed the experiments; JQW, JJF, WJ, TXY, XYH, and QQW performed the experiments; QL, JQW, FJP, ZFY, and WG analyzed the data; JQW and JJF prepared figures and/or tables; JQW drafted manuscripts; QL, JQW, and YYT revised the manuscript. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of Weifang Medical University (No. 2022YX65). Written informed consents were waived by the Medical Ethics Committee of Weifang Medical University because the study involved no more than minimal risk of harm to subjects. All experiments involving mice were performed in accordance with the Chinese National Laboratory Animal Guideline for Ethical Review of Animal Welfare and approved by the Institutional Animal Care and Use Committee of Weifang Medical University (No. 2022SDL162). All applicable international, national, and institutional guidelines for the care and use of animals were strictly followed.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
 Table S1. Biochemical traits of C. pyruviciproducens isolates.
Additional file
 2: Fig. S1. Subcutaneous abscesses of BALB/c mice after subcutaneous infection with C. pyruviciproducens strain WYJY-01 were re-isolated, cultured, and identified by Gram staining.
Additional file 3: Fig. S2.
Quality control of whole genome extraction and sequencing data of C. pyruviciproducens strain WYJY-01.
Additional file 4: Fig. S3.
Genome annotation of C. pyruviciproducens strain WYJY-01.
Additional file 5: Fig. S4.
Venn diagram based on gene families clustered from the Pfam database of the genomes of three strains of C. pyruviciproducens strain WYJY-01, ATCC BAA-1742T, and UMB0763.
Additional file 6: Fig. S5.
The GO database annotations of the genome of C. pyruviciproducens strain WYJY-01 and its differential genes with that of ATCC BAA-1742T.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, J., Feng, J., Jia, W. et al. Genomic and phenotypic analysis of a novel clinical isolate of Corynebacterium pyruviciproducens. BMC Microbiol 23, 385 (2023). https://doi.org/10.1186/s12866-023-03075-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12866-023-03075-6