Next Article in Journal
A Genomics-Based Discovery of Secondary Metabolite Biosynthetic Gene Clusters in the Potential Novel Strain Streptomyces sp. 21So2-11 Isolated from Antarctic Soil
Previous Article in Journal
Periodontal Inflammation and Dysbiosis Relate to Microbial Changes in the Gut
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain

1
Department of Medical and Biological Sciences, The Catholic University of Korea, Bucheon 14662, Republic of Korea
2
Department of Biotechnology, The Catholic University of Korea, Bucheon 14662, Republic of Korea
3
Bio-Resources Bank Division, Nakdonggang National Institute of Biological Resources, Sangju 37242, Republic of Korea
4
Performance Innovation Division, Nakdonggang National Institute of Biological Resources, Sangju 37242, Republic of Korea
5
Department of Pathology, Yale University School of Medicine, New Haven, CT 06510, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Microorganisms 2024, 12(6), 1227; https://doi.org/10.3390/microorganisms12061227
Submission received: 3 June 2024 / Revised: 14 June 2024 / Accepted: 15 June 2024 / Published: 18 June 2024
(This article belongs to the Section Microbial Biotechnology)

Abstract

:
Dermacoccus barathri is the first reported pathogen within the Dermacoccus genus to cause a catheter-related bloodstream infection, which occurred in 2015. In this study, the complete genome assembly of Dermacoccus barathri was constructed, and the complete genome of Dermacoccus barathri FBCC-B549 consists of a single chromosome (3,137,745 bp) without plasmids. The constructed genome of D. barathri was compared with those of two closely related species within the Dermacoccus genus. D. barathri exhibited a pattern similar to Dermacoccus abyssi in terms of gene clusters and synteny analysis. Contrary to previous studies, biosynthetic gene cluster (BGC) analysis for predicting secondary metabolites revealed the presence of the LAP biosynthesis pathway in the complete genome of D. barathri, predicting the potential synthesis of the secondary metabolite plantazolicin. Furthermore, an analysis to investigate the potential pathogenicity of D. barathri did not reveal any antibiotic resistance genes; however, nine virulence factors were identified in the Virulence Factor Database (VFDB). According to these matching results in the VFDB, despite identifying a few factors involved in biofilm formation, further research is required to determine the actual impact of D. barathri on pathogenicity. The complete genome of D. barathri is expected to serve as a valuable resource for future studies on D. barathri, which currently lack sufficient genomic sequence information.

1. Introduction

The Dermacoccus genus comprises coccoid-shaped bacteria characterized as Gram-positive, lacking capsules, non-spore-forming, aerobic, and commonly found in soil, water, and cured meat products [1]. In addition, some strains are part of the human skin microbiome, and their decrease has been associated with Atopic dermatitis, according to Wongpiyabovorn et al. in 2019 [2]. Despite inhabiting a wide range of terrestrial, aquatic, and biological habitats, only four fully identified species are known within the Dermacoccus genus: Dermacoccus barathri, Dermacoccus abyssi, Dermacoccus nishinomiyaensis, and Dermacoccus profundi [3,4]. Dermacoccus species have traditionally been considered nonpathogenic commensals that live on human skin. However, D. barathri was identified as the first potential pathogen in Dermacoccus species because it has been implicated in catheter-related bloodstream infections (CRBSIs) [5]. Beginning with D. barathri, a few species of the Dermacoccus genus have caused human infections. Specifically, there have been five reported cases of infections caused by D. nishinomiyaensis [6,7,8,9,10] and one reported case of infection caused by Dermacoccus sp. [11].
CRBSIs often result from bacterial colonization near the catheter insertion site, where biofilm formation protects bacteria from antibiotics and immune responses [12]. Gram-positive bacteria, including D. barathri, are more frequently associated with catheter infections than Gram-negative bacteria [12]. Notably, D. barathri was reported in 2015 as the first Dermacoccus species to cause a CRBSI in a 7-year-old boy with an implantable central venous access device [5]. Furthermore, D. barathri isolated from this patient showed the ability to form biofilms after approximately 10 days of culture [5]. Despite being recognized as a novel pathogen within the Dermacoccus genus, the complete genome of and genetic information on D. barathri are lacking. Therefore, constructing the complete genome of D. barathri can address this limitation and aid in further research.
In this study, whole-genome sequencing of Dermacoccus barathri FBCC-B549 was conducted using Nanopore Flongle, resulting in the assembly of the complete genome of D. barathri. Through comparative genomic analysis with two closely related species, D. abyssi, and D. nishinomiyaensis, the evolutionary relationships and unique genetic features of D. barathri were elucidated. In addition, genes associated with antibiotic resistance and virulence factors, particularly those linked to biofilm formation, were scrutinized to evaluate the potential pathogenicity of D. barathri. Furthermore, secondary metabolites that had not been previously discovered before constructing the complete genome were identified in this study. The comprehensive analysis of the complete genome of D. barathri sheds light on its genetic characteristics and potential pathogenic traits. This study promises to enrich our molecular-level understanding and serve as a valuable resource for further research in related fields.

2. Materials and Methods

2.1. Whole-Genome Sequencing

The DNA of D. barathri FBCC-B549 was extracted using the Omega Bio Mag-Bind Universal Pathogen Kit (Omega Bio-tek, Norcress, GA, USA) following the manufacturer’s protocol. The Nanopore Flongle sequencing library was constructed using a Nanopore SQK-LSK114 kit (Oxford Nanopore Technologies, Oxford, UK). Base calling for the generated signal from the nanopore was conducted using Guppy v6.5.7 [13] with a super-accuracy model and CUDA acceleration. Trimming of adapter sequences in the generated reads was performed using Porechop_ABI v0.5.0 [14], and the trimmed reads were assembled using FLYE v2.9.1 [15] with the nano-hq parameter. For a more accurate genome sequence, Medaka v1.8.0 (https://github.com/nanoporetech/medaka, accessed on 14 August 2023) was employed for polishing with a super-accuracy calling model. The quality of the assembled genome was evaluated using BUSCO v5.4.7 [16] with the micrococcales_odb10 database.

2.2. Gene prediction and Functional Annotation

The gene annotation and characteristics of D. barathri FBCC-B549 were analyzed using Prokka, COGclassifier, and the Rapid Annotation based on Subsystem Technology (RAST) server. Prokka v1.14.5 [17] and Proksee [18] were utilized for gene prediction, annotation, and the construction of a circular map of the assembled genome. COGclassifier v1.0.5 [19] was employed to annotate genes and classify them according to COG. RAST [20] was also utilized for gene annotation and subsystem category-wise analyses.

2.3. Comparative Genome Analysis

Comparative genomic analyses were conducted using the complete genomes of closely related species within the Dermacoccus genus, namely, D. abyssi and D. nishinomiyaensis, which are the species that have available complete genomes. Whole-genome comparisons were performed using Mauve v2.4.0 [21]. To align the starting point of the circular genome, each genome of D. barathri, D. abyssi, and D. nishinomiyaensis was manually modified. The locally collinear block (LCB) weight was set to 1625. Gene clusters of the three species were compared and annotated using OrthoVenn3 [22] with default parameters (E-value 0.01; inflation value 1.50). Additional comparative genome analyses were performed using Roary v3.13.0 [23]. Biosynthetic gene clusters (BGCs) were predicted using antiSMASH v7.1.0 [24] with default parameters.

2.4. Prediction of Antimicrobial Resistance Gene and Potential Virulence Factors

The potentially harmful antibiotic resistance and virulence of D. barathri were predicted. Resfinder v4.4.2 [25] and AMRFinderPlus v3.11.26 [26] were used to identify antibiotic resistance factors using the default parameters. Potential virulence prediction and toxin factor analysis were performed against the VFDB [27] using the DIAMOND v2.1.8 tool [28]. Toxicity factors were filtered under relaxed conditions with an identity above 70% and query coverage above 60% [29].

3. Results and Discussion

3.1. Complete Genome Construction and Gene Annotation

Using Nanopore Flongle sequencing data, the complete genome of D. barathri was successfully assembled using Flye and Medaka. The genome of D. barathri is composed of a single chromosome with no plasmids, totaling 3,137,745 base pairs and possessing a GC content of 68.2%, which is typical for Actinomycetes (Table 1). High-GC regions are prone to PCR bias, resulting in reduced amplification efficiency under standard PCR conditions. This bias likely contributed to the challenges in obtaining the complete genome of D. barathri using Illumina sequencing, which relies on PCR. However, nanopore sequencing, based on single-molecule sequencing, is believed to be less affected by PCR bias, facilitating the assembly of the high-GC-content genome of D. barathri. To evaluate the quality of the assembled genome, BUSCO was used with the micrococcales_odb10 database, revealing a high level of completeness at 98.5%. This included 528 complete and single-copy BUSCOs, one complete and duplicated BUSCO, and eight missing BUSCOs (Table 2). Previous studies have suggested that nanopore sequencing can solely produce qualified genomes without additional short-read data for specific research purposes, and the polishing step is critical [30]. While tools like Homopolish [31] have been effective in enhancing the base accuracy of nanopore sequencing by utilizing multiple public genomes of conspecific microorganisms, this specific polishing tool could not be utilized for D. barathri due to the limited availability of sequences in the public database. Nevertheless, the accuracy enhancement achieved using nanopore protein pores, associated reagents, base calling (super-accuracy model), and polishing algorithms, particularly with the default Medaka polishing tool provided by Nanopore, was deemed sufficient through BUSCO analysis. As technology advances, it is increasingly expected that the construction of qualified genomes suitable for various research purposes will become more feasible solely through the use of nanopore technology.
For the gene annotation of the complete genome, Prokka annotation was performed, identifying 2924 CDSs, 56 tRNAs, and 9 rRNAs (Table 1). These values are similar to those of the closely related species D. abyssi and D. nishinomiyaensis. A circular genome map was constructed based on the gene information of the complete D. barathri genome using Proksee (Figure 1).
Additionally, functional annotation of the genes was conducted using the COG classifier and RAST annotation. In the COG classifier results, COG annotation was completed for 2256 of the 2924 CDS (77.15%). Categories related to basic nutrient metabolism, such as “Translation, Ribosomal structure, and biogenesis” (212), “Amino Acid Transport and Metabolism” (209), “Carbohydrate Transport and Metabolism” (171), “Lipid Transport and Metabolism” (152), and “Coenzyme Transport and Metabolism” (152), accounted for a significant portion (Figure 2). Considering the requirement for the secretion of polysaccharides, proteins, amino acids, and various RNAs and DNAs outside the cell for biofilm matrix formation, these basic nutrient-metabolism-related categories are expected to play important roles in biofilm formation. In addition to COG annotation, subsystem-based annotation using the RAST server was performed, with 25% of the annotated genes assigned to subsystems and 75% unassigned. Similar to the COG annotation results, a considerable number of genes involved in basic metabolic processes, such as “Amino Acids and Derivatives” (228), “Protein Metabolism” (170), and “Carbohydrates” (138), were identified to be part of subsystems (Figure 2). Furthermore, through RAST, 29 genes were identified in the “Virulence, Disease, and Defense” category, with 20 genes classified under the “Resistance to Antibiotics and Toxic Compounds” subcategory and 9 genes under the “Invasion and Intracellular Resistance” subcategory. This suggests that the species possesses genes that may contribute to its pathogenicity.

3.2. Comparative Genome Analysis

Using Mauve, a comparison of the LCBs within the genomes of the closely related species D. abyssi and D. nishinomiyaensis was conducted (Figure 3). The regions depicted as blocks in Figure 3 represent LCBs, denoting conserved sequence regions devoid of internal genomic rearrangements. The overall GC content and genome length of the three species were similar, and most gene sequences and orders were preserved in the major LCB arrays, indicating good collinearity. While the major LCB arrays of the three species were similar, the crossing patterns of the lines linking the conserved regions in each genome appeared to be more complex in D. nishinomiyaensis, suggesting a closer similarity to D. abyssi in terms of overall genome structure. According to the phylogenetic tree of the Dermaoccus genus reported in the literature, D. barathri is phylogenetically closer to D. abyssi than D. nishinomiyaensis [1].
This phylogenetic relationship is consistent with the results reported by Roary. D. barathri shared 2176 genes with D. abyssi and 38 with D. nishinomiyaensis. Analyses using OrthoVenn3 showed a similar pattern (Figure 4). As illustrated in the Venn diagram in Figure 4, the total number of proteins in D. barathri, D. abyssi, and D. nishinomiyaensis was 8444, with 2679 clusters. Furthermore, 2056 clusters were shared among the three species. These clusters constituted a significant proportion of individual clusters and primarily consisted of genes related to fundamental bacterial biological functions. Similar to previous results, D. barathri shared more clusters with D. abyssi, sharing 2469 clusters. D. barathri shared 2170 clusters with D. nishinomiyaensis, which was fewer than the number of clusters shared by D. abyssi. Although the patterns between the results of Roary and OrthoVenn3 appear similar, variations in the values suggest differences in the default parameters for each tool.
By predicting secondary metabolites using antiSMASH, biosynthetic gene clusters (BGCs) in D. barathri and related species were investigated. In the complete genome of D. barathri FBCC-B549, six BGCs for secondary metabolites were identified: the NI-siderophore, NAPAA, terpene, RiPP-like, LAP, and ectoine. Among these, the NAPAA, terpene, and RiPP-like clusters were found to be common to D. nishinomiyaensis, D. abyssi, and D. barathri. Additionally, D. abyssi and D. barathri shared two additional clusters: the NI-siderophore and ectoine. In contrast to previous studies [32], the analysis of the newly constructed complete genome of D. barathri FBCC-B549 using antiSMASH revealed the presence of the LAP (Linear azol(in)e-containing peptide) cluster, which was not previously reported in D. barathri literature. LAPs belong to a specific subgroup of ribosomally synthesized and post-translationally modified peptides (RiPPs) characterized by the presence of azole or azoline heterocycles within the linear peptides [33]. RiPPs and LAPs undergo pathways in which precursor peptides mature and become active through the actions of modifying enzymes and transporters [33]. Among these modifying enzymes, YcaO serves as a biosynthetic enzyme that forms azole and azoline heterocycles [33]. The LAP cluster in D. barathri also contains the YcaO protein gene and shares the same precursor peptide sequence with that of Dermacoccus sp. DE3. Previous studies suggested that Dermacoccus sp. DE3 synthesizes plantazolicin (PZN) as a secondary metabolite [33], indicating the potential for D. barathri to also synthesize PZN. PZN was initially discovered in the plant growth-promoting bacterium Bacillus amyloliquefaciens FZB42 [34]. PZN acts as a ribosome-targeting inhibitor and inhibits translation by binding to bacterial 70S ribosomes. It has been reported to function as an antibiotic with a very narrow spectrum, targeting Bacillus anthracis, the causative agent of anthrax [33,35]. Therefore, D. barathri appears to have the potential to synthesize antimicrobial substances effective against anthrax-causing bacteria. This characteristic is expected to be further explored in subsequent research.

3.3. Potential Pathogenicity of D. barathri

To explore the potential pathogenicity of D. barathri (FBCC-B549), an analysis of antibiotic resistance genes and virulence factors within its complete genome was conducted. Unfortunately, our initial attempt to identify antimicrobial resistance genes using Resfinder and AMRfinder with default parameters did not yield matches to known antibiotic resistance genes. Subsequently, for virulence factor analysis, protein sequences obtained from the Prokka annotation of the complete genome of D. barathri FBCC-B549 were employed and compared against VFDB using DIAMOND. Our criteria for matches were set at identity ≥70% and query coverage ≥60%, in line with established parameters from previous research [29]. Through this analysis, nine virulence factors were identified and categorized into nutritional/metabolic factors, regulation, adherence, stress survival, and others (Table 3). As previously documented, in the first case of a catheter-related bloodstream infection, D. barathri was found to possess the ability to form biofilms [5]. Assuming this trait to be common among D. barathri strains, potential factors contributing to biofilm formation among the detected virulence factors were sought. Factors within the adherence category, such as GroEL and EF-Tu, are known for their involvement in the initial adhesion phase of biofilm formation across various pathogens [36]. EF-Tu serves as a prominent constituent of bacterial cell walls, exhibiting properties similar to those of adhesion factors, inducing inflammatory responses, and participating in energy metabolism, functioning similarly to moonlighting proteins [37]. Moreover, EF-Tu plays a role similar to that of amyloid-like proteins discovered in Gallibacterium during biofilm formation [38], contributing to immune evasion, initiation, maturation, and attachment stages of biofilm development [39]. Furthermore, GroEL participates in early aggregation [40], and beyond its essential role in mycolate synthesis during Mycobacterium biofilm maturation [41], it has also been implicated in biofilm formation in Hemophilus influenzae [35]. LeuD, present in the nutritional/metabolic factor category, was upregulated in E. coli biofilms, and experimental evidence showed that LeuD functions as a major regulator of biofilm growth as it becomes S-nitrosylated [42,43]. Finally, KatA, within the stress survival category, safeguards biofilms and organisms from UVA radiation [44]. In addition to these biofilm-related factors, aceA, a gene encoding isocitrate lyase, regulates the TCA cycle and synthesizes carbohydrates from fatty acids via the initial enzyme of the glyoxylate shunt, contributing to sustaining infection; therefore, it is classified as a virulent factor [45]. Additionally, SigA, SigH, and RegX3 within the regulation category are known to regulate virulence and stress resistance operons [26,46,47]. Biofilm formation is associated not only with host immune responses and antibiotic therapy but also with the recurrence and persistence of chronic diseases. This also suggests its involvement in the pathogenicity of opportunistic bacteria [48]. According to the National Institutes of Health, up to 80% of human bacterial infections are associated with biofilm-associated microorganisms [49,50,51]. D. barathri has also been reported to form biofilms and cause catheter-related bloodstream infections [4]. Additionally, despite not being extracted from the reported case, the strain used in this study possessed genetic factors that are potentially involved in biofilm formation. Further research is required to ascertain the biofilm-forming capabilities of these genes. Furthermore, subsequent studies should investigate whether biofilm formation is a common phenotypic trait in D. barathri. The complete genome of D. barathri constructed in this study is expected to serve as a valuable resource for enhancing our understanding of the characteristics and genetic background of D. barathri species lacking an available genome.

4. Conclusions

In this study, the complete genome of D. barathri (strain FBCC-B549) was successfully assembled using Nanopore Flongle technology, thereby filling a significant gap in the available public genomic information for this species. The genome of D. barathri has a notably high GC content, and our assembly reveals a single chromosome spanning 3,137,745 base pairs. From this complete genome, fundamental genetic information was extracted, and genes involved in secondary metabolite synthesis were identified. These findings were then compared with sequences from other species within the same genus. Furthermore, antibiotic resistance genes and virulence factors were analyzed to evaluate the potential pathogenicity of D. barathri at the genetic level. The establishment of the complete genome sequence of D. barathri in this study is expected to significantly advance our molecular-level understanding of this microbial species. Moreover, it will serve as a valuable resource for further research in related fields.

Author Contributions

W.K. designed the study and directed all analysis processes. Y.K. and H.K. conducted the entire analysis and wrote the main manuscript. J.K., S.W.N., J.-H.H., E.J.C. and M.S. reviewed and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the Nakdonggang National Institute of Biological Resources (NNIBR), funded by the Ministry of Environment (MOE) of the Republic of Korea (grant number: NNIBR20241108) and the Korea Environment Industry & Technology Institute (KEITI), through the project to make multi-ministerial national biological research resources a more advanced program funded by the MOE (grant number: 2021003420001).

Data Availability Statement

The original data presented in the study are openly available in the NCBI database Bioproject under accession number PRJNA1118424.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Pathom-Aree, W.; Nogi, Y.; Ward, A.C.; Horikoshi, K.; Bull, A.T.; Goodfellow, M. Dermacoccus barathri sp. nov. and Dermacoccus profundi sp. nov., novel actinomycetes isolated from deep-sea mud of the Mariana Trench. Int. J. Syst. Evol. Microbiol. 2006, 56, 2303–2307. [Google Scholar] [CrossRef] [PubMed]
  2. Wongpiyabovorn, J.; Soonthornchai, W.; Wilantho, A.; Palasuk, M.; Payungporn, S.; Sodsai, P.; Poomipak, W.; Weschawalit, S.; Ruchusatsawat, K.; Baillie, G.S. Effect of tarcolimus on skin microbiome in atopic dermatitis. Allergy 2019, 74, 1400–1406. [Google Scholar] [CrossRef]
  3. Schoch, C.L.; Ciufo, S.; Domrachev, M.; Hotton, C.L.; Kannan, S.; Khovanskaya, R.; Leipe, D.; Mcveigh, R.; O’Neill, K.; Robbertse, B. NCBI Taxonomy: A comprehensive update on curation, resources and tools. Database 2020, 2020, baaa062. [Google Scholar] [CrossRef] [PubMed]
  4. Parte, A.C.; Sardà Carbasse, J.; Meier-Kolthoff, J.P.; Reimer, L.C.; Göker, M. List of Prokaryotic names with Standing in Nomenclature (LPSN) moves to the DSMZ. Int. J. Syst. Evol. Microbiol. 2020, 70, 5607–5612. [Google Scholar] [CrossRef] [PubMed]
  5. Takahashi, N.; Shinjoh, M.; Tomita, H.; Fujino, A.; Sugita, K.; Katohno, Y.; Kuroda, T.; Kikuchi, K. Catheter-related blood stream infection caused by Dermacoccus barathri, representing the first case of Dermacoccus infection in humans. J. Infect. Chemother. 2015, 21, 613–616. [Google Scholar] [CrossRef] [PubMed]
  6. Joron, C.; Roméo, B.; Le Fleche-Mateos, A.; Rames, C.; El Samad, Y.; Hamdad, F. Dermacoccus nishinomiyaensis as a cause of persistent paediatric catheter-related bacteraemia. Clin. Microbiol. Infect. 2019, 25, 1054–1055. [Google Scholar] [CrossRef] [PubMed]
  7. Seifu, W.D.; Gebissa, A.D. Prevalence and antibiotic susceptibility of Uropathogens from cases of urinary tract infections (UTI) in Shashemene referral hospital, Ethiopia. BMC Infect. Dis. 2018, 18, 30. [Google Scholar] [CrossRef] [PubMed]
  8. Ogbac, M.K.E. Community Acquired Pneumonia Secondary to Dermacoccus nishinomiyaensis Infection: A Case Report. Am. J. Respir. Crit. Care Med. 2021, 203, A3978. [Google Scholar] [CrossRef]
  9. Tanaka, A.; Watanabe, Y.; Ito, C.; Murata, M.; Shinjo, H.; Otsuka, Y.; Takeda, A. Successful treatment of peritoneal dialysis-related peritonitis caused by Dermacoccus nishinomiyaensis. CEN Case Rep. 2019, 8, 183–187. [Google Scholar] [CrossRef] [PubMed]
  10. Triguya, J.; Tolina, A.L.; Peñaa, S.; Ballestera, C.; Gallardoa, Á.; Fixmanb, M.J.A. Dermacoccus nishinomiyaensis brain abscesses as the first manifestation of chronic granulomatous disease. Arch. Argent Pediatr. 2023, 121, e202202804. [Google Scholar]
  11. Ivaska, L.; Alyazidi, R.; Hoang, L.; Goldfarb, D.M. Dermacoccus sp. isolated from a brain abscess in a 4-year-old child. J. Infect. Chemother. 2019, 25, 1070–1073. [Google Scholar] [CrossRef] [PubMed]
  12. Gahlot, R.; Nigam, C.; Kumar, V.; Yadav, G.; Anupurba, S. Catheter-related bloodstream infections. Int. J. Crit. Illn. Inj. Sci. 2014, 4, 162–167. [Google Scholar] [CrossRef] [PubMed]
  13. Wick, R.R.; Judd, L.M.; Holt, K.E. Performance of neural network basecalling tools for Oxford Nanopore sequencing. Genome Biol. 2019, 20, 129. [Google Scholar] [CrossRef] [PubMed]
  14. Bonenfant, Q.; Noé, L.; Touzet, H. Porechop_ABI: Discovering unknown adapters in ONT sequencing reads for downstream trimming. BioRxiv 2022, 2007, 499093. [Google Scholar] [CrossRef]
  15. Kolmogorov, M.; Yuan, J.; Lin, Y.; Pevzner, P.A. Assembly of long, error-prone reads using repeat graphs. Nat. Biotechnol. 2019, 37, 540–546. [Google Scholar] [CrossRef] [PubMed]
  16. Manni, M.; Berkeley, M.R.; Seppey, M.; Simão, F.A.; Zdobnov, E.M. BUSCO update: Novel and streamlined workflows along with broader and deeper phylogenetic coverage for scoring of eukaryotic, prokaryotic, and viral genomes. Mol. Biol. Evol. 2021, 38, 4647–4654. [Google Scholar] [CrossRef] [PubMed]
  17. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef] [PubMed]
  18. Grant, J.R.; Enns, E.; Marinier, E.; Mandal, A.; Herman, E.K.; Chen, C.-y.; Graham, M.; Van Domselaar, G.; Stothard, P. Proksee: In-depth characterization and visualization of bacterial genomes. Nucleic Acids Res. 2023, 51, W484–W492. [Google Scholar] [CrossRef] [PubMed]
  19. Shimoyama, Y. COGclassifier: A Tool for Classifying Prokaryote Protein Sequences into COG Functional Category. 2022. Available online: https://github.com/moshi4/COGclassifier (accessed on 9 January 2024).
  20. Overbeek, R.; Olson, R.; Pusch, G.D.; Olsen, G.J.; Davis, J.J.; Disz, T.; Edwards, R.A.; Gerdes, S.; Parrello, B.; Shukla, M. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2014, 42, D206–D214. [Google Scholar] [CrossRef] [PubMed]
  21. Darling, A.C.; Mau, B.; Blattner, F.R.; Perna, N.T. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14, 1394–1403. [Google Scholar] [CrossRef] [PubMed]
  22. Sun, J.; Lu, F.; Luo, Y.; Bie, L.; Xu, L.; Wang, Y. OrthoVenn3: An integrated platform for exploring and visualizing orthologous data across genomes. Nucleic Acids Res. 2023, 51, W397–W403. [Google Scholar] [CrossRef] [PubMed]
  23. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef] [PubMed]
  24. Medema, M.H.; Blin, K.; Cimermancic, P.; De Jager, V.; Zakrzewski, P.; Fischbach, M.A.; Weber, T.; Takano, E.; Breitling, R. antiSMASH: Rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011, 39, W339–W346. [Google Scholar] [CrossRef] [PubMed]
  25. Bortolaia, V.; Kaas, R.S.; Ruppe, E.; Roberts, M.C.; Schwarz, S.; Cattoir, V.; Philippon, A.; Allesoe, R.L.; Rebelo, A.R.; Florensa, A.F. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 2020, 75, 3491–3500. [Google Scholar] [CrossRef] [PubMed]
  26. Feldgarden, M.; Brover, V.; Gonzalez-Escalona, N.; Frye, J.G.; Haendiges, J.; Haft, D.H.; Hoffmann, M.; Pettengill, J.B.; Prasad, A.B.; Tillman, G.E. AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Sci. Rep. 2021, 11, 12728. [Google Scholar] [CrossRef] [PubMed]
  27. Chen, L.; Yang, J.; Yu, J.; Yao, Z.; Sun, L.; Shen, Y.; Jin, Q. VFDB: A reference database for bacterial virulence factors. Nucleic Acids Res. 2005, 33, D325–D328. [Google Scholar] [CrossRef] [PubMed]
  28. Buchfink, B.; Reuter, K.; Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 2021, 18, 366–368. [Google Scholar] [CrossRef] [PubMed]
  29. Umanets, A.; Surono, I.S.; Venema, K. I am better than I look: Genome based safety assessment of the probiotic Lactiplantibacillus plantarum IS-10506. BMC Genom. 2023, 24, 518. [Google Scholar] [CrossRef] [PubMed]
  30. Lee, J.Y.; Kong, M.; Oh, J.; Lim, J.; Chung, S.H.; Kim, J.-M.; Kim, J.-S.; Kim, K.-H.; Yoo, J.-C.; Kwak, W. Comparative evaluation of Nanopore polishing tools for microbial genome assembly and polishing strategies for downstream analysis. Sci. Rep. 2021, 11, 20740. [Google Scholar] [CrossRef] [PubMed]
  31. Huang, Y.-T.; Liu, P.-Y.; Shih, P.-W. Homopolish: A method for the removal of systematic errors in nanopore sequencing by homologous polishing. Genome Biol. 2021, 22, 95. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, X.; Ruan, Y.; Liu, W.; Chen, Q.; Gu, L.; Guo, A. Whole genome sequencing and genome annotation of Dermacoccus abyssi strain HZAU 226 isolated from spoiled eggs. Genomics 2021, 113, 1199–1206. [Google Scholar] [CrossRef] [PubMed]
  33. Burkhart, B.J.; Schwalen, C.J.; Mann, G.; Naismith, J.H.; Mitchell, D.A. YcaO-dependent posttranslational amide activation: Biosynthesis, structure, and function. Chem. Rev. 2017, 117, 5389–5456. [Google Scholar] [CrossRef] [PubMed]
  34. Scholz, R.; Molohon, K.J.; Nachtigall, J.; Vater, J.; Markley, A.L.; Süssmuth, R.D.; Mitchell, D.A.; Borriss, R. Plantazolicin, a novel microcin B17/streptolysin S-like natural product from Bacillus amyloliquefaciens FZB42. J. Bacteriol. 2011, 193, 215–224. [Google Scholar] [CrossRef] [PubMed]
  35. Molohon, K.J.; Saint-Vincent, P.M.; Park, S.; Doroghazi, J.R.; Maxson, T.; Hershfield, J.R.; Flatt, K.M.; Schroeder, N.E.; Ha, T.; Mitchell, D.A. Plantazolicin is an ultranarrow-spectrum antibiotic that targets the Bacillus anthracis membrane. ACS Infect. Dis. 2015, 2, 207–220. [Google Scholar] [CrossRef] [PubMed]
  36. Gallaher, T.K.; Wu, S.; Webster, P.; Aguilera, R. Identification of biofilm proteins in non-typeable Haemophilus influenzae. BMC Microbiol. 2006, 6, 65. [Google Scholar] [CrossRef] [PubMed]
  37. Gannesen, A.V.; Zdorovenko, E.L.; Botchkova, E.A.; Hardouin, J.; Massier, S.; Kopitsyn, D.S.; Gorbachevskii, M.V.; Kadykova, A.A.; Shashkov, A.S.; Zhurina, M.V. Composition of the biofilm matrix of Cutibacterium acnes acneic strain RT5. Front. Microbiol. 2019, 10, 456374. [Google Scholar] [CrossRef] [PubMed]
  38. López-Ochoa, J.; Montes-García, J.F.; Vázquez, C.; Sánchez-Alonso, P.; Pérez-Márquez, V.M.; Blackall, P.J.; Vaca, S.; Negrete-Abascal, E. Gallibacterium elongation factor-Tu possesses amyloid-like protein characteristics, participates in cell adhesion, and is present in biofilms. J. Microbiol. 2017, 55, 745–752. [Google Scholar] [CrossRef] [PubMed]
  39. Erskine, E.; MacPhee, C.E.; Stanley-Wall, N.R. Functional amyloid and other protein fibers in the biofilm matrix. J. Mol. Biol. 2018, 430, 3642–3656. [Google Scholar] [CrossRef] [PubMed]
  40. Arai, T.; Ochiai, K.; Senpuku, H. Actinomyces naeslundii GroEL-dependent initial attachment and biofilm formation in a flow cell system. J. Microbiol. Methods 2015, 109, 160–166. [Google Scholar] [CrossRef] [PubMed]
  41. Ojha, A.; Anand, M.; Bhatt, A.; Kremer, L.; Jacobs, W.R.; Hatfull, G.F. GroEL1: A dedicated chaperone involved in mycolic acid biosynthesis during biofilm formation in mycobacteria. Cell 2005, 123, 861–873. [Google Scholar] [CrossRef] [PubMed]
  42. Barraud, N.; Létoffé, S.; Beloin, C.; Vinh, J.; Chiappetta, G.; Ghigo, J.-M. Lifestyle-specific S-nitrosylation of protein cysteine thiols regulates Escherichia coli biofilm formation and resistance to oxidative stress. Npj Biofilms Microbiomes 2021, 7, 34. [Google Scholar] [CrossRef] [PubMed]
  43. Jones, C.J.; Grotewold, N.; Wozniak, D.J.; Gloag, E.S. Pseudomonas aeruginosa initiates a rapid and specific transcriptional response during surface attachment. J. Bacteriol. 2022, 204, e00086-22. [Google Scholar] [CrossRef]
  44. Pezzoni, M.; Pizarro, R.A.; Costa, C.S. Protective role of extracellular catalase (KatA) against UVA radiation in Pseudomonas aeruginosa biofilms. J. Photochem. Photobiol. B Biol. 2014, 131, 53–64. [Google Scholar] [CrossRef] [PubMed]
  45. Chung, J.C.; Rzhepishevska, O.; Ramstedt, M.; Welch, M. Type III secretion system expression in oxygen-limited Pseudomonas aeruginosa cultures is stimulated by isocitrate lyase activity. Open Biol. 2013, 3, 120131. [Google Scholar] [CrossRef] [PubMed]
  46. Gomez, M.; Doukhan, L.; Nair, G.; Smith, I. sigA is an essential gene in Mycobacterium smegmatis. Mol. Microbiol. 1998, 29, 617–628. [Google Scholar] [CrossRef] [PubMed]
  47. Raman, S.; Song, T.; Puyang, X.; Bardarov, S.; Jacobs, W.R., Jr.; Husson, R.N. The alternative sigma factor SigH regulates major components of oxidative and heat stress responses in Mycobacterium tuberculosis. J. Bacteriol. 2001, 183, 6119–6125. [Google Scholar] [CrossRef] [PubMed]
  48. Hall-Stoodley, L.; Stoodley, P. Evolving concepts in biofilm infections. Cell. Microbiol. 2009, 11, 1034–1043. [Google Scholar] [CrossRef] [PubMed]
  49. Römling, U.; Balsalobre, C. Biofilm infections, their resilience to therapy and innovative treatment strategies. J. Intern. Med. 2012, 272, 541–561. [Google Scholar] [CrossRef] [PubMed]
  50. da Silva, R.A.; Afonina, I.; Kline, K.A. Eradicating biofilm infections: An update on current and prospective approaches. Curr. Opin. Microbiol. 2021, 63, 117–125. [Google Scholar] [CrossRef] [PubMed]
  51. Singh, A.; Amod, A.; Pandey, P.; Bose, P.; Pingali, M.S.; Shivalkar, S.; Varadwaj, P.K.; Sahoo, A.K.; Samanta, S.K. Bacterial biofilm infections, their resistance to antibiotics therapy and current treatment strategies. Biomed. Mater. 2022, 17, 022003. [Google Scholar] [CrossRef] [PubMed]
Figure 1. A circular genome map of D. barathri generated using Prokka and Proksee, illustrating the locations of coding sequences (CDSs) within the genome.
Figure 1. A circular genome map of D. barathri generated using Prokka and Proksee, illustrating the locations of coding sequences (CDSs) within the genome.
Microorganisms 12 01227 g001
Figure 2. (A) COG annotation results for D. barathri, with each letter representing a class of COG categories. (B) RAST annotation results for D. barathri, with the pie chart displaying Subsystem Feature Counts and the number of genes associated with specific functions.
Figure 2. (A) COG annotation results for D. barathri, with each letter representing a class of COG categories. (B) RAST annotation results for D. barathri, with the pie chart displaying Subsystem Feature Counts and the number of genes associated with specific functions.
Microorganisms 12 01227 g002
Figure 3. Locally collinear block (LCB) comparison among complete genomes of D. barathri, D. abyssi (A), and D. nishinomiyaensis (B). Colored blocks indicate aligned genome sequences, highlighting homologous and conserved regions.
Figure 3. Locally collinear block (LCB) comparison among complete genomes of D. barathri, D. abyssi (A), and D. nishinomiyaensis (B). Colored blocks indicate aligned genome sequences, highlighting homologous and conserved regions.
Microorganisms 12 01227 g003
Figure 4. Orthologous gene cluster analysis using OrthoVenn3, depicting shared orthologous clusters among three Dermacoccus species. Numbers adjacent to species indicate the total clusters in each list.
Figure 4. Orthologous gene cluster analysis using OrthoVenn3, depicting shared orthologous clusters among three Dermacoccus species. Numbers adjacent to species indicate the total clusters in each list.
Microorganisms 12 01227 g004
Table 1. Summary of assembly statistics for the D. barathri FBCC-B549 whole genome. “bp” means “base pair”, which is the unit used to measure the length of DNA.
Table 1. Summary of assembly statistics for the D. barathri FBCC-B549 whole genome. “bp” means “base pair”, which is the unit used to measure the length of DNA.
Species NameDermacoccus barathri FBCC-B549
NCBI taxon ID322601
TaxonomyBacteria; Actinomycetota; Actinomycetes; Micrococcales; Dermacoccaceae; Dermacoccus
Genome Size (bp)3,137,745
GC content in the DNA68.2 mol% G+C
Number of Genome Sequences1 Circular
Number of Plasmids0
Number of Coding Sequences2924
Number of rRNAs9
Number of tRNAs56
Table 2. Summary of BUSCO analysis for D. barathri genome.
Table 2. Summary of BUSCO analysis for D. barathri genome.
CategoriesNo. GenesPercentage
Complete52998.5
Complete single-copy BUSCOs52898.3
Complete duplicate BUSCOs10.2
Fragmented BUSCOs00
Missing BUSCOs81.5
Table 3. Virulence factors identified in D. barathri using blastp with the VFDB database.
Table 3. Virulence factors identified in D. barathri using blastp with the VFDB database.
CategoryVirulence Factors
Nutritional/metabolic factorsLeuD, SugC
RegulationSigA, SigH, RegX3
AdherenceEF-Tu, GroEL
Stress survivalKatA
OthersAceA
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, Y.; Kim, H.; Kim, J.; Han, J.-H.; Chung, E.J.; Nam, S.W.; Shin, M.; Kwak, W. The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain. Microorganisms 2024, 12, 1227. https://doi.org/10.3390/microorganisms12061227

AMA Style

Kim Y, Kim H, Kim J, Han J-H, Chung EJ, Nam SW, Shin M, Kwak W. The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain. Microorganisms. 2024; 12(6):1227. https://doi.org/10.3390/microorganisms12061227

Chicago/Turabian Style

Kim, Yeha, Hyaekang Kim, Jina Kim, Ji-Hye Han, Eu Jin Chung, Seung Won Nam, Miyoung Shin, and Woori Kwak. 2024. "The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain" Microorganisms 12, no. 6: 1227. https://doi.org/10.3390/microorganisms12061227

APA Style

Kim, Y., Kim, H., Kim, J., Han, J. -H., Chung, E. J., Nam, S. W., Shin, M., & Kwak, W. (2024). The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain. Microorganisms, 12(6), 1227. https://doi.org/10.3390/microorganisms12061227

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop