Jin Zhao
Xinli Shen†
Lulu Jin
Songjun Ji
Xinling Pan
94% of researchers rate our articles as excellent or good
Learn more about the work of our research integrity team to safeguard the quality of each article we publish.
Find out more
- Department of Biomedical Sciences Laboratory, Affiliated Dongyang Hospital of Wenzhou Medical University, Dongyang, Zhejiang, China
Background: Mycobacterium seoulense (M. seoulense) is an emerging pathogen increasingly associated with infections; however, its resistance phenotypes and genomic characteristics remain largely unknown.
Methods: Seven M. seoulense isolates were collected from clinical samples. Drug susceptibility testing was conducted using Sensititre™ SLOMYCO2 susceptibility plates. Whole genome sequencing and supporting bioinformatics analyses were performed to analyze the genomic features.
Results: All M. seoulense isolates (n=7) exhibited growth on 7H10 agar medium containing thiophenecarboxylic acid hydrazide or p-Nitrobenzoic acid, with marked diversity in growth rates in liquid culture. All strains exhibited high minimum inhibitor concentrations (MICs) for minocycline (>8 μg/mL), doxycycline (>8 μg/mL), and amikacin (16-32 μg/mL). The MICs for linezolid, rifabutin, moxifloxacin, ciprofloxacin, streptomycin, clarithromycin, and rifampicin varied among the isolates. High levels of genomic diversity were noted among these strains concerning genome-called single nucleotide polymorphisms and average nucleotide identity. In total, 4,282 genes were shared across all genomes, while 315 unique genes were restricted to one strain. Comparative genomic analysis identified two unique virulence genes encoding a catalase enzyme and a protein involved in capsule biosynthesis and transport. Additionally, all M. seoulense strains demonstrated the ability to survive within macrophages.
Conclusion: The clinical M. seoulense isolates analyzed in this study exhibited varying levels of antibiotic susceptibility, suggesting the potential need for susceptibility testing to guide clinical treatment. Genomic features of these isolates indicated that they are likely pathogenic non-tuberculous mycobacterium, highlighting a need for closer epidemiological monitoring.
Introduction
Most non-tuberculous mycobacteria (NTM) are widely distributed in the environment, including dust and water sources. Although typically nonpathogenic, they can cause infections in immunocompetent individuals, resulting in pulmonary, soft tissue, and lymphatic infections (Nguyen et al., 2024; Sun et al., 2024). Among immunosuppressed patients, NTM can also cause disseminated infections (Kitazawa et al., 2024; Nithirungruang et al., 2024), posing an immense therapeutic challenge to those with preexisting diseases. With the rapid development of molecular detecting techniques, over 200 species of NTM have been identified to date (https://lpsn.dsmz.de/genus/mycobacterium) (Armstrong et al., 2023), many of which have been implicated in clinical infections. Several treatment guidelines have been proposed for NTM infections (Haworth et al., 2017; Kurz et al., 2020), but there remains limited evidence regarding the most appropriate therapeutic regimens when combatting infections caused by less common NTM species (Lange et al., 2022). Therefore, phenotypic and genomic investigations are warranted for those novel NTM species with rising clinical prevalence (Zhang et al., 2013).
Mycobacterium seoulense (M. seoulense), which was first reported in Korea in 2007 (Mun et al., 2007), has also been found to cause pulmonary infection in Japan and China (Zhao et al., 2022; Nakamoto et al., 2024). Four cases of M. seoulense infection were identified in our previous study through radiological examination. However, anti-mycobacterial treatment was not initiated due to limited clinical experience, a lack of public data on M seoulense infections, and patients’ reluctance to receive antibiotic therapy. Only the case from Japan received antibiotic therapy, and had a optimal outcome. While the growth features and biochemical features of this NTM have been reported (Mun et al., 2007), the resistance phenotype for M. seoulense remains unknown. Like most NTM, the addition of thiophene carboxylic acid hydrazide (TCH) and p-nitrobenzoic acid (PNB) to the culture medium has no impact on the growth of M. seoulense, potentially providing a means of differentiating it from M. nebraskense, a biochemically similar species. Although two genomes for M. seoulense are available (https://www.ncbi.nlm.nih.gov/datasets/genome/?taxon=386911), no genomic analysis has been conducted, leaving some uncertainty regarding the relevant virulence factors and core genes for this species. Therefore, in this study, seven strains isolated from clinical samples in our hospital were sequenced, and drug susceptibility test was performed to provide clinicians with insights to help guide the management of patients infected by M. seoulense.
Methods
Bacterial culture
Strains were recovered from -80°C by subculturing them in Middlebrook 7H9 broth medium (BD Bioscience, USA) containing 10% ADC broth enrichment medium (Shanghai Yiyan Bio-technology Co. Ltd, Shanghai, China) and 0.05% Tween-80. Isolates were incubated at 37°C until reaching an optical density at 600 nm (OD600) of 0.6-1.0 as measured with an Implen (Germany) instrument. A shaking culture was then performed at 37°C and 250 rpm, and growth curves were obtained by testing the OD600 every three days before day 21 and every six days after day 21. In addition, the medium was imaged after an 18-day incubation.
To assess the morphology of these isolates on the Middlebrook 7H10 agar medium, cultures with an OD600 between 0.6 and 1.0 were diluted to an OD600 of 0.02 using a complete 7H9 medium, and 5 µL of the culture was spotted onto a 7H10 complete agar plate supplemented with 10% glycerol. The 7H10 agar plates were supplemented with 5 μg/mL of TCH or 500 μg/mL of PNB. The plates were incubated in an incubator at 37°C and imaged until the colonies had formed after 21 days.
Drug susceptibility test
Bacterial cultures with an OD600 between 0.6 and 1.0 were obtained as described above, followed by dilution with phosphate-buffered saline (PBS) to a final OD600 of 0.2 (about 0.5 McFarland) (Penuelas-Urquides et al., 2013). The bacterial suspension was then diluted 100-fold using Mueller-Hinton medium, and 100 μL was transferred into the Sensititre™ SLOMYCO2 drug susceptibility test plates (Thermo Fisher, USA). These plates were sealed with parafilm and cultured in an incubator at 37°C. Readings were taken after 7 days, with the minimum inhibitory concentration defined as the lowest drug concentration at which no visible bacterial sediment was observed at the bottom of the well. Two replicate wells were established for this drug susceptibility assay. M. intracellulare ATCC 13950 was used as the reference strain.
Genomic DNA extraction and next-generation sequencing
Bacteria in the early exponential phase were harvested via centrifugation, and cell pellets were used to extract genomic DNA based on the provided instructions (Tiangen, Beijing, China). High-quality DNA was used to prepare a library with a NEXTFLEX Rapid DNA-Seq Kit (Revvity, MA, USA). Sequencing was performed on a NovaSeq™ X Plus (Illumina, CA, USA), and raw data was analyzed to remove low-quality reads using the fastp software (Chen et al., 2018). The high-quality reads were assembled with Unicycler software (Wick et al., 2017), and the sequences were uploaded to the NCBI database (PRJNA1196008). Genes were annotated using Prokka (Seemann, 2014), and pan-genome analyses were conducted using Roary (Page et al., 2015).
Genomic analysis of the assembled sequences
M. seoulense strain JCM16018 was the reference genome (Accession No: GCA_010731595.1) for calling core single nucleotide polymorphisms (coreSNPs) of clinical isolates using the Snippy software as described previously (Seemann, 2015). Following the removal of SNPs in recombinant regions using Gubbins (Croucher et al., 2015), the coreSNP distribution matrix was used to construct a Phylogenic tree with IQ-TREE (Nguyen et al., 2015). Genomic alignment was performed via BLAST and visualized in the form of rings with BRIG (Alikhan et al., 2011). The identity percentages between any two strains were evaluated based on average nucleotide identity (ANI) (Yoon et al., 2017).
Third-generation sequencing
Genomic DNA was fragmented into sequences ~10 kb in length, after which an SMRT Bell library was constructed with an SMRT Bell Prep Kit 3.0 (Pacific Biosciences, CA). The raw data was sequenced using the PacBio Sequel IIe sequencing platform. Genomic assembly was performed with the Unicycler software, and the results were uploaded to the NCBI database (PRJNA1195676). The whole genome was aligned to the JCM16018 strain, and the arranged genomes were analyzed using Mauve (Darling et al., 2004).
Analyses of intracellular survival in macrophages
THP-1 cells were purchased from the China Center for Type Culture Collection (Shanghai, China) and cultured in RPMI 1640 medium (Gibco, USA) containing 10% fetal bovine serum (Vivacell, China). Cells were seeded into 24-well plates (1 x 105/well) and differentiated into macrophages via treatment for 24 h with 60 ng/mL phorbol 12-myristate 13-acetate. Bacterial cultures with an OD600 between 0.6 and 1.0 were harvested by centrifugation (2,000 rpm, 5 min). After washing thrice with PBS, the bacterial pellet was resuspended as a single-cell suspension in RPMI 1640 medium via repeated pipetting. The concentration was adjusted into an OD600 of 0.6 (equivalent to 1×108/mL) and co-cultured with the macrophage at a multiplicity of infection of 10. After three hours of incubation, the extracellular bacteria were washed away twice with PBS and killed by adding 30 µg/mL gentamycin for one hour. The infected cells were incubated in a complete RPMI 1640 medium supplemented with 5 µg/mL gentamycin in a 37°C 5% CO2 incubator. After incubation for 24 h, the cells were lysed with sterile deionized water and serially diluted with PBS. The diluted lysed cell solutions were spread on the 7H10 agar medium plates with 10% glycerol. After incubation for 20 days, the colony-forming units (CFUs) on these plates were counted. The survival rate was calculated as the ratio of CFU obtained after infection for 24 h to the initial CFU obtained after infection. This experiment was repeated three times.
Statistical analysis
The survival rates of M. seoulense within macrophages were compared to that of M. intracellulare ATCC13950 strain using a t-test in GraphPad Prism (version 9.3.1). Statistical significance was defined as *, P<0.05; **, P<0.01.
Results
Biological features of M. seoulense in culture
A total of seven M. seoulense strains were isolated from six patients (Supplementary Table 1), including four previously described cases. The two additional cases in this study had M. seoulense in their bronchoalveolar lavage fluid specimens. In the 7H9 broth medium, these clinical isolates exhibited different growth rates (Figure 1A). Specifically, the 16NTM017 and 20NTM121 isolates grew more rapidly than the 23NTM055 and 21NTM074 isolates, with these four strains entering the stationary phase after about 21 days, similar to M. intracellulare ATCC 13950. The 20NTM058, 20NTM062, and 22NTM124 strains grew gradually and entered the stationary phase until 40 days. The color of the liquid culture medium for all M. seoulense isolates was brown-yellow, while that for M. intracellular was opalescent (Figure 1B).
Figure 1
Figure 1. The growth features of Mycobacterium seoulense in culture medium. (A) The growth curve in the 7H9 Middlebrook medium with 10% ADC and 0.05% Tween 80; (B) The color of 7H9 medium liquid medium; (C) The morphology of Mycobacterium seoulense on the 7H10 agar medium with 10% glycerol. 1, 16NTM017; 2, 20NTM058; 3, 20NTM062; 4, 20NTM121; 5, 21NTM074; 6, 21NTM160 (M. colombiense); 7, 22NTM124; 8, 23NTM055; 9, ATCC13950.
On solid medium, all isolates were able to grow on the 7H10 complete medium with or without TCH or PNB (Figure 1C). However, the clones grown on PNB were white, whereas those grown on other mediums were yellow. Moreover, strains 16NTM017 and 20NTM121 exhibited a darker color compared to other strains.
Drug susceptibility test of M. seoulense isolates
All M. seoulense strains exhibited high MICs for Amikacin, ranging from 16-32 μg/mL (Table 1). For doxycycline and minocycline, all the strains exhibited MIC levels higher than the upper limit of detection. These clinical isolates presented varying levels of susceptibility to moxifloxacin, ciprofloxacin, streptomycin, and clarithromycin, as indicated by at least four MICs for these seven strains.
Table 1
Table 1. Distribution of MICs for M. seoulense strains.
Phylogenetic analyses of M. seoulense isolates
Based on the coreSNPs detected in individual strains, an established phylogenetic tree revealed the clustering of the JCM16017 and DSM 44998 strains (Figure 2A), with an average nucleotide identity of 99.99% (Supplementary Figure 1). The two strains derived from the same case (20NTM062 and 20NTM058) were also clustered with ANI of 99.98%. Two strains that were isolated from different cases were close to each other in the phylogenetic tree. The differing SNP distributions among the analyzed strains provided support for the phylogenetic tree results (Figure 2B). For the clustered strains, the strains derived from the same case exhibited eight coreSNPs that differed from one another, while those from different cases exhibited 181-1643 differing coreSNPs. In addition, 22NTM124 and 23NTM055 were distant from the other strains (Figure 2A), with approximately 50,000 coreSNPs differing from the other strains (Figure 2B). The diversity of coreSNPs was similar to the ANI distributions between any two strains (Supplementary Figure 1).
Figure 2
Figure 2. The phylogenetic analysis of Mycobacterium seoulense strains. (A) phylogenetic tree based on the coreSNPs; (B) distribution matrix of coreSNPs among Mycobacterium seoulense strains.
Pan-genome analyses of M. seoulense isolates
In order to establish the core genes of M. seoulense, nine isolates, including our clinical isolates and reference strains, were subjected to alignment and pan-genome analyses. Relative to the 16NTM017 strain, diversity hot spots were observed across the genome, including from 1965.5-2391 kb and 3600-3700 kb (Figure 3A). These regions exhibited a higher percentage of GC content than other regions. For the genes in these strains, 4,282 conserved genes were present in all genomes, while 315 unique genes were present in a single strain (Figure 3B, Supplementary Table 2). There were 552 genes associated with lipid transport and metabolism, accounting for 15.4% of the 3,595 genes with available COG annotation. This category was associated with the most diversity in terms of the number of related genes among these clinical isolates, while the numbers of genes associated with “cell motility”, “extracellular structure”, and “intracellular trafficking, secretion, and vesicular transport” were largely conserved (Figure 3C).
Figure 3
Figure 3. The pan-genome analysis of Mycobacterium seoulense strains. (A) Genomic alignment by circular ring; (B) A summary of genes in the pan-genome analysis; (C) The number of genes with COG annotation of genes in the pan-genome analysis.
Genomic analyses of M. seoulense
In order to comprehensively analyze the genomic features of 16NTM017, assembled genomic data was obtained via third-generation sequencing. Relative to strain JCM16018, 16NTM017 exhibited an ANI of 99.03% (Figure 4A), with several long insertion sequences and one rearranged region in the genome (Figure 4B). The longest inserted sequence corresponded to position 1965.5 kb to 2391 kb, encoding 434 CDS.
Figure 4
Figure 4. The genomic alignment of 16NTM017 and JCM16018 based on the assembled whole genome sequence. (A) The global view of genomic alignment; (B) the rearrangement of genomic regions between the aligned genomes. The red regions indicated for inserted regions in 16NTM017 comparing to JCM16018.
Virulence factors in M. seoulense
When compared with three reference strains (Mycobacterium tuberculosis, Mycobacterium bovis BCG strain, and Mycobacterium avium complex (MAC)), the reference M. seoulense strain and our clinical isolates exhibited a range of common virulence factors (Figure 5). Specifically, the M. seoulense strains were found to harbor two unique genes encoding a catalase enzyme and a protein involved in capsule biosynthesis and transport. In addition, the M. seoulense strains were found to lack rmlA, mbtL, fad23, and eis, which are associated with the GPL locus, a trehalose-recycling ABC transporter, mycobactin, and the enhanced intracellular survival protein. The distribution of several genes also varied among clinical isolates (Supplementary Table 3).
Figure 5
Figure 5. The virulence factor-related genes for Mycobacterium seoulense. The genes with yellow background were specifically absent in Mycobacterium seoulense strains comparing to other mycobacterial species; The genes with green background were uniquely present in the Mycobacterium seoulense strains comparing to other mycobacterial species.
Analyses of intracellular survival in macrophages
The M. intraellulare reference strain ATCC13950 exhibited a survival rate of ~60% in macrophages for 24 h (Figure 6). Relative to M. intracellulare, the clinical M. seoulense isolates demonstrated enhanced intracellular survival in macrophages, albeit with variations among these strains.
Figure 6
Figure 6. The intracellular survival of Mycobacterium seoulense strains within the macrophage. The survival rate was defined as the ratio of CFU after 24 hours infection to CFU at the time of infection. *, P<0.05; **, P<0.01.
Discussion
In this study, whole genome sequencing was performed to characterize this novel mycobacterial species and to identify virulence factors related to pathogenicity. These efforts provided an opportunity to analyze the genomic features of clinical M. seoulense isolates. Compared to the commonly isolated NTM species, such as the Mycobacterium avium complex, M. seoulense strains exhibited a yellow pigmentation in culture, indicating carotenoid production (Robledo et al., 2011; Tran et al., 2020). In mycobacteria, the carotenoid production is used mainly for taxonomic identification. However, in other bacterial species, carotenoids possess antioxidant, anti-inflammatory, anti-cancer, and antimicrobial properties (Tran et al., 2020). Although most NTM can be isolated from living environments and medical devices, M. seoulense has not been reported in the soil or water (Abbas et al., 2024; Zhang et al., 2024). One explanation is that it was overlooked in epidemiological studies, as M. seoulense infections tend to be less common in clinical practice relative to other NTM species (Zhao et al., 2022; Nakamoto et al., 2024). Owing to the absence of data regarding the recommended threshold of SNPs in defining clonal transmission in this species, evidence of transmission between these cases could not be obtained based on the available epidemiological data. Based on the SNPs that differed between the strains isolated from the same patient (8 SNPs), the genomic diversity among clonal M. seoulense strains may be greater than for other species belonging to the MAC (Namkoong et al., 2021; Wetzstein et al., 2024). Different drug susceptibility phenotypes for several antibiotic agents were indicated for high heterogeneity of these two strains. Interestingly, two M. seoulense strains with available genomes (DSMD44998 and JCM16018) shared highly similar genomes despite their submission by different individuals. Moreover, their identity percentages were higher than those from the same case in this study, indicating the possibility of a common source. This finding warrants further study, as these microbes were isolated from different countries, potentially indicating interregional transmission. Unfortunately, the possible epidemiological connections for these two strains could not be established based on the information submitted to the database.
The MICs for M. seoulense were first reported in this study, but the resistance phenotypes could not be established owing to limited data on these MICs and recommended breakpoints (Brown-Elliott and Woods, 2019). Based on its growth feature, M. seoulense is classified as a slow NTM (Mun et al., 2007). When using the breakpoints for slow-growth NTM species other than for MAC or M. kansasii used in a previous study (Liu et al., 2021), M. seoulense strains could be not efficiently inhibited by doxycycline, minocycline, or ciprofloxacin. Based on susceptibility results, a multi-drug regimen containing clofazimine or clarithromycin in combination with rifampicin and moxifloxacin may be effective for M. seoulense infections, offering a potential alternative drug regimen in clinical practice (Nakamoto et al., 2024). In cases of suboptimal treatment response, drug susceptibility test for rifampicin and clarithromycin may be necessary, as the MICs vary among isolates. However, further research is required to establish resistance breakpoints, evaluate the efficacy of different drug regimens in clinical settings, and assess treatment outcomes for M. seoulense infections (Huh et al., 2019; Kwon et al., 2019).
There was substantial diversity among the analyzed M. seoulense genomes in this study based on identity and annotated genes. Despite two strains isolated from the same region, a higher level of diversity was noted among strains. Pan-genomic analyses allowed us to focus on conserved genes in molecular diagnosis and pathogenicity investigations while assessing accessory genes helps to study intra-species diversity (Zekic et al., 2018; Zakham et al., 2021). Further research revealed that the diversity among these clinical strains was related to the number of genes, especially for those annotated in the category “lipid transport and metabolism”. Variations in lipid metabolism-related genes among clinical isolates are also observed in different lineages of Mycobacterium tuberculosis, contributing to differences in virulence and transmissibility (O’Neill et al., 2015; Moopanar and Mvubu, 2020). In addition, pathogenic mycobacteria exhibit a contraction in genes related to lipid and secondary metabolite biosynthesis and metabolism (Zhu et al., 2023). Given that lipids are a critical source of energy supporting bacterial survival under nutrient-limiting conditions, the diversity of lipid-related genes may play a significant role in mycobacterial virulence and could have implications for drug development (Bailo et al., 2015). The inconsistent number of transport-associated genes may have contributed to the diverse pathogenicity in the cell infection model and the varying levels of drug resistance (Samuels et al., 2022). In contrast, the genes associated with nucleotide transport and metabolism were more conserved than those related to carbohydrate, amino acid, and lipid transport and metabolism. However, further investigation will be necessary owing to the limited data on the genomic analyses of metabolic activity for NTM (Karmakar and Sur, 2024).
The distribution of virulence factor-related genes in M. seoulense may provide insight into the mechanisms that may be responsible for infection, helping to guide the selection of molecular targets for drug design efforts. Mycobactin is important for supporting mycobacterial survival under iron-stress conditions in host macrophages through siderophore formation (Luo et al., 2005; Shyam et al., 2021; Shankar and Akhter, 2024). There are 14 genes in the mbt gene clusters that are involved in the biosynthesis of mycobactin, with three of these genes (mbtL, mbtN, and mbtJ) being absent in the M. seoulense strains relative to the MTB reference strain (H37Rv). The synthesis of the core mycobactin scaffold is mediated by enzymes encoded by mbt A-J, while the transfer of the lipophilic aliphatic chain to the ϵ-amino group of the lysine fragment is facilitated by enzymes encoded by mbtK-N, resulting in the potential loss of function of mycobactin in M. seoulense (Shyam et al., 2021). In addition, faD33 (mbtM) was absent in M. seoulense, resulting in the absence of inhibitory effects by the antituberculosis agents (C16-AMS nucleoside) targeting mbtM (Guillet et al., 2016; Shyam et al., 2021). Moreover, pks2 and fad23 were absent in M. seoulense, which would impair the biosynthesis of Sulfolipid-1 (Yan et al., 2023). This protein can restrict the intracellular survival of MTB in human macrophages in a species-specific manner (Gilmore et al., 2012), and it can also serve as a nociceptive molecule that activates nociceptive neurons and induces coughing in infected animal models (Ruhl et al., 2020). The rmlA gene, which was not detected in the genome of M. seoulense, has been identified as a virulence factor enhancing the intracellular growth of MTB (Sassetti et al., 2003) and may be a target for antituberculosis agents (Xiao et al., 2021). Eis, a secreted protein not present in M. seoulense but present in other mycobacterial species, suppresses host innate immune defenses by modulating autophagy, inflammation, and cell death in a redox-dependent manner (Shin et al., 2010). M. seoulense strains lacking certain genes mentioned above were able to survive within macrophages, other virulence factors may have played a compensatory role in enhancing intracellular survival, such as those associated with the ESX-1 secretion system (Guo et al., 2021; Bo et al., 2023). Notably, katA is uniquely present in M. seoulense, while KatG is conserved across all the mycobacterial species (Supplementary Table 3). Both katA and katG encode catalase enzymes that are essential for detoxifying antimicrobial hydrogen peroxide (Stacy et al., 2014; Santander et al., 2018), with katG also contributing to isoniazid resistance (Ofori-Anyinam et al., 2024).
This study has several limitations. While certain virulence factor-encoding genes were uniquely detected or absent in the genomes of these M. seoulense strains, their precise roles in intracellular survival and pathogenicity remain unclear. The current in vitro macrophage infection data serve as an initial reference, but future studies utilizing gene-editing techniques will be necessary to elucidate their functional significance. Additionally, information regarding possible interactions between cases could not be obtained owing to the prolonged period over which these clinical isolates were obtained, although their living addresses were not geographically linked. Furthermore, clinical follow-up after species identification was not conducted, resulting in missing data on patient prognosis following M. seoulense infection.
Conclusion
Based on the MICs measured, M. seoulense strains can be efficiently inhibited by rifampicin, rifabutin, linezolid, and clarithromycin but not by doxycycline, minocycline, or ciprofloxacin. Pan-genomic analyses revealed the core genes across all the strains, and the distribution of virulence factor-related genes differed from those of other common NTM species. Our findings may not only help clinicians select appropriate drug therapies for M. seoulense-related infections but may also provide possible targets for studies of the mechanisms underlying the pathogenicity of these mycobacteria.
Data availability statement
All the genomic sequences could be accessible from NCBI data with No. PRJNA1196008 and No. PRJNA1195676.
Ethics statement
All the enrolled cases gave written consent to provide the strains and clinical information in an anonymous format at their admission to the hospital. This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee and Institutional Review Board of Dongyang People’s Hospital (2023-YX-272).
Author contributions
JZ: Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing. XS: Data curation, Validation, Writing – original draft, Writing – review & editing. LJ: Funding acquisition, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. SJ: Data curation, Resources, Validation, Writing – original draft, Writing – review & editing. XP: Conceptualization, Funding acquisition, Investigation, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Medical Health Science and Technology Project of Zhejiang Provincial Health Commission (2023KY386) and the Science and Technology Bureau of Jinhua (2022-3-013 and 2022-3-011).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2025.1553591/full#supplementary-material
Supplementary Figure 1 | The distribution of average nucleotide identity among Mycobacterium seoulense strains.
Abbreviations
MTB, Mycobacterium tuberculosis; NTM, non-tuberculous mycobacterium; MIC, minimum inhabitation concentration; SNP, single nucleotide polymorphism; CDS, coding sequence; CFU, colony forming units; TCH, thiophene carboxylic acid hydrazide; PNB, p-Nitrobenzoic acid.
References
Abbas, M., Khan, M. T., Iqbal, Z., Ali, A., Eddine, B. T., Yousaf, N., et al. (2024). Sources, transmission and hospital-associated outbreaks of nontuberculous mycobacteria: a review. Future Microbiol. 19, 715–740. doi: 10.2217/fmb-2023-0279
Alikhan, N. F., Petty, N. K., Ben Zakour, N. L., Beatson, S. A. (2011). BLAST Ring Image Generator (BRIG): simple prokaryote genome comparisons. BMC Genomics 12, 402. doi: 10.1186/1471-2164-12-402
Armstrong, D. T., Eisemann, E., Parrish, N. (2023). A brief update on mycobacterial taxonomy, 2020 to 2022. J. Clin. Microbiol. 61, e0033122. doi: 10.1128/jcm.00331-22
Bailo, R., Bhatt, A., Ainsa, J. A. (2015). Lipid transport in Mycobacterium tuberculosis and its implications in virulence and drug development. Biochem. Pharmacol. 96, 159–167. doi: 10.1016/j.bcp.2015.05.001
Bo, H., Moure, U. A. E., Yang, Y., Pan, J., Li, L., Wang, M., et al. (2023). Mycobacterium tuberculosis-macrophage interaction: Molecular updates. Front. Cell Infect. Microbiol. 13, 1062963. doi: 10.3389/fcimb.2023.1062963
Brown-Elliott, B. A., Woods, G. L. (2019). Antimycobacterial susceptibility testing of nontuberculous mycobacteria. J. Clin. Microbiol. 57, e00834-19. doi: 10.1128/JCM.00834-19
Chen, S., Zhou, Y., Chen, Y., Gu, J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–ii90. doi: 10.1093/bioinformatics/bty560
Croucher, N. J., Page, A. J., Connor, T. R., Delaney, A. J., Keane, J. A., Bentley, S. D., et al. (2015). Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 43, e15. doi: 10.1093/nar/gku1196
Darling, A. C., Mau, B., Blattner, F. R., Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704
Gilmore, S. A., Schelle, M. W., Holsclaw, C. M., Leigh, C. D., Jain, M., Cox, J. S., et al. (2012). Sulfolipid-1 biosynthesis restricts Mycobacterium tuberculosis growth in human macrophages. ACS Chem. Biol. 7, 863–870. doi: 10.1021/cb200311s
Guillet, V., Galandrin, S., Maveyraud, L., Ladeveze, S., Mariaule, V., Bon, C., et al. (2016). Insight into structure-function relationships and inhibition of the fatty acyl-AMP ligase (FadD32) orthologs from mycobacteria. J. Biol. Chem. 291, 7973–7989. doi: 10.1074/jbc.M115.712612
Guo, Q., Bi, J., Wang, H., Zhang, X. (2021). Mycobacterium tuberculosis ESX-1-secreted substrate protein EspC promotes mycobacterial survival through endoplasmic reticulum stress-mediated apoptosis. Emerg. Microbes Infect. 10, 19–36. doi: 10.1080/22221751.2020.1861913
Haworth, C. S., Banks, J., Capstick, T., Fisher, A. J., Gorsuch, T., Laurenson, I. F., et al. (2017). British Thoracic Society guidelines for the management of non-tuberculous mycobacterial pulmonary disease (NTM-PD). Thorax 72, ii1–ii64. doi: 10.1136/thoraxjnl-2017-210929
Huh, H. J., Kim, S. Y., Jhun, B. W., Shin, S. J., Koh, W. J. (2019). Recent advances in molecular diagnostics and understanding mechanisms of drug resistance in nontuberculous mycobacterial diseases. Infect. Genet. Evol. 72, 169–182. doi: 10.1016/j.meegid.2018.10.003
Karmakar, M., Sur, S. (2024). Unlocking the Mycobacteroides abscessus pan-genome using computational tools: insights into evolutionary dynamics and lifestyle. Antonie Van Leeuwenhoek 118, 30. doi: 10.1007/s10482-024-02042-z
Kitazawa, T., Yamamoto, A., Nakayama, S., Kawase, K., Wakabayashi, Y. (2024). Disseminated nontuberculous mycobacterium infection during treatment of multiple myeloma: A case report and review of the literature. Intern. Med. doi: 10.2169/internalmedicine.4234-24
Kurz, S. G., Zha, B. S., Herman, D. D., Holt, M. R., Daley, C. L., Ruminjo, J. K., et al. (2020). Summary for clinicians: 2020 clinical practice guideline summary for the treatment of nontuberculous mycobacterial pulmonary disease. Ann. Am. Thorac. Soc. 17, 1033–1039. doi: 10.1513/AnnalsATS.202003-222CME
Kwon, Y. S., Daley, C. L., Koh, W. J. (2019). Managing antibiotic resistance in nontuberculous mycobacterial pulmonary disease: challenges and new approaches. Expert Rev. Respir. Med. 13, 851–861. doi: 10.1080/17476348.2019.1638765
Lange, C., Bottger, E. C., Cambau, E., Griffith, D. E., Guglielmetti, L., van Ingen, J., et al. (2022). Consensus management recommendations for less common non-tuberculous mycobacterial pulmonary diseases. Lancet Infect. Dis. 22, e178–ee90. doi: 10.1016/S1473-3099(21)00586-7
Liu, C. F., Song, Y. M., He, W. C., Liu, D. X., He, P., Bao, J. J., et al. (2021). Nontuberculous mycobacteria in China: incidence and antimicrobial resistance spectrum from a nationwide survey. Infect. Dis. Poverty 10, 59. doi: 10.1186/s40249-021-00844-1
Luo, M., Fadeev, E. A., Groves, J. T. (2005). Mycobactin-mediated iron acquisition within macrophages. Nat. Chem. Biol. 1, 149–153. doi: 10.1038/nchembio717
Moopanar, K., Mvubu, N. E. (2020). Lineage-specific differences in lipid metabolism and its impact on clinical strains of Mycobacterium tuberculosis. Microb. Pathog. 146, 104250. doi: 10.1016/j.micpath.2020.104250
Mun, H. S., Kim, H. J., Oh, E. J., Kim, H., Bai, G. H., Yu, H. K., et al. (2007). Mycobacterium seoulense sp. nov., a slowly growing scotochromogenic species. Int. J. Syst. Evol. Microbiol. 57, 594–599. doi: 10.1099/ijs.0.64744-0
Nakamoto, K., Fujiwara, K., Morimoto, K., Ohta, K. (2024). Nontuberculous mycobacterial pulmonary disease caused by Mycobacterium seoulense. Respirol Case Rep. 12, e70049. doi: 10.1002/rcr2.v12.10
Namkoong, H., Omae, Y., Asakura, T., Ishii, M., Suzuki, S., Morimoto, K., et al. (2021). Genome-wide association study in patients with pulmonary Mycobacterium avium complex disease. Eur. Respir. J. 58, 1902269. doi: 10.1183/13993003.02269-2019
Nguyen, M. H., Haas, M. K., Kasperbauer, S. H., Calado Nogueira de Moura, V., Eddy, J. J., Mitchell, J. D., et al. (2024). Nontuberculous mycobacterial pulmonary disease: patients, principles, and prospects. Clin. Infect. Dis. 79, e27–e47. doi: 10.1093/cid/ciae421
Nguyen, L. T., Schmidt, H. A., von Haeseler, A., Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Nithirungruang, P., Tanpowpong, P., Getsuwan, S., Boonsathorn, S. (2024). Disseminated macrolide-resistant mycobacterium intracellulare infection in a child with autoimmune lymphoproliferative disorder: A case report and literature review. Am. J. Trop. Med. Hyg 111, 1051–1055. doi: 10.4269/ajtmh.24-0076
O’Neill, M. B., Mortimer, T. D., Pepperell, C. S. (2015). Diversity of Mycobacterium tuberculosis across Evolutionary Scales. PloS Pathog. 11, e1005257. doi: 10.1371/journal.ppat.1005257
Ofori-Anyinam, B., Hamblin, M., Coldren, M. L., Li, B., Mereddy, G., Shaikh, M., et al. (2024). Catalase activity deficiency sensitizes multidrug-resistant Mycobacterium tuberculosis to the ATP synthase inhibitor bedaquiline. Nat. Commun. 15, 9792. doi: 10.1038/s41467-024-53933-8
Page, A. J., Cummins, C. A., Hunt, M., Wong, V. K., Reuter, S., Holden, M. T., et al. (2015). Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics 31, 3691–3693. doi: 10.1093/bioinformatics/btv421
Penuelas-Urquides, K., Villarreal-Trevino, L., Silva-Ramirez, B., Rivadeneyra-Espinoza, L., Said-Fernandez, S., de Leon, M. B. (2013). Measuring of Mycobacterium tuberculosis growth. A correlation of the optical measurements with colony forming units. Braz. J. Microbiol. 44, 287–289. doi: 10.1590/S1517-83822013000100042
Robledo, J. A., Murillo, A. M., Rouzaud, F. (2011). Physiological role and potential clinical interest of mycobacterial pigments. IUBMB Life 63, 71–78. doi: 10.1002/iub.v63.2
Ruhl, C. R., Pasko, B. L., Khan, H. S., Kindt, L. M., Stamm, C. E., Franco, L. H., et al. (2020). Mycobacterium tuberculosis sulfolipid-1 activates nociceptive neurons and induces cough. Cell 181, 293–305 e11. doi: 10.1016/j.cell.2020.02.026
Samuels, A. N., Wang, E. R., Harrison, G. A., Valenta, J. C., Stallings, C. L. (2022). Understanding the contribution of metabolism to Mycobacterium tuberculosis drug tolerance. Front. Cell Infect. Microbiol. 12, 958555. doi: 10.3389/fcimb.2022.958555
Santander, R. D., Figas-Segura, A., Biosca, E. G. (2018). Erwinia amylovora catalases KatA and KatG are virulence factors and delay the starvation-induced viable but non-culturable (VBNC) response. Mol. Plant Pathol. 19, 922–934. doi: 10.1111/mpp.2018.19.issue-4
Sassetti, C. M., Boyd, D. H., Rubin, E. J. (2003). Genes required for mycobacterial growth defined by high density mutagenesis. Mol. Microbiol. 48, 77–84. doi: 10.1046/j.1365-2958.2003.03425.x
Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153
Seemann, T. (2015). “Snippy: rapid haploid variant calling and core SNP phylogeny,” in GitHub. Available at: github.com/tseemann/snippy.
Shankar, G., Akhter, Y. (2024). Stealing survival: Iron acquisition strategies of Mycobacteriumtuberculosis. Biochimie. 227, 37–60. doi: 10.1016/j.biochi.2024.06.006
Shin, D. M., Jeon, B. Y., Lee, H. M., Jin, H. S., Yuk, J. M., Song, C. H., et al. (2010). Mycobacterium tuberculosis eis regulates autophagy, inflammation, and cell death through redox-dependent signaling. PloS Pathog. 6, e1001230. doi: 10.1371/journal.ppat.1001230
Shyam, M., Shilkar, D., Verma, H., Dev, A., Sinha, B. N., Brucoli, F., et al. (2021). The mycobactin biosynthesis pathway: A prospective therapeutic target in the battle against tuberculosis. J. Med. Chem. 64, 71–100. doi: 10.1021/acs.jmedchem.0c01176
Stacy, A., Everett, J., Jorth, P., Trivedi, U., Rumbaugh, K. P., Whiteley, M. (2014). Bacterial fight-and-flight responses enhance virulence in a polymicrobial infection. Proc. Natl. Acad. Sci. U.S.A. 111, 7819–7824. doi: 10.1073/pnas.1400586111
Sun, Y., Zhang, C., Lu, B., Chen, J., Pan, X. (2024). Case report: Mycobacterium chimaera-induced lymph node infection in a patient with chronic myeloproliferative neoplasm misdiagnosed as tuberculous lymphadenitis. Front. Public Health 12, 1387722. doi: 10.3389/fpubh.2024.1387722
Tran, T., Dawrs, S. N., Norton, G. J., Virdi, R., Honda, J. R. (2020). Brought to you courtesy of the red, white, and blue-pigments of nontuberculous mycobacteria. AIMS Microbiol. 6, 434–450. doi: 10.3934/microbiol.2020026
Wetzstein, N., Diricks, M., Anton, T. B., Andres, S., Kuhns, M., Kohl, T. A., et al. (2024). Clinical and genomic features of Mycobacterium avium complex: a multi-national European study. Genome Med. 16, 86. doi: 10.1186/s13073-024-01359-8
Wick, R. R., Judd, L. M., Gorrie, C. L., Holt, K. E. (2017). Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PloS Comput. Biol. 13, e1005595. doi: 10.1371/journal.pcbi.1005595
Xiao, G., Alphey, M. S., Tran, F., Pirrie, L., Milbeo, P., Zhou, Y., et al. (2021). Next generation Glucose-1-phosphate thymidylyltransferase (RmlA) inhibitors: An extended SAR study to direct future design. Bioorg Med. Chem. 50, 116477. doi: 10.1016/j.bmc.2021.116477
Yan, M., Cao, L., Zhao, L., Zhou, W., Liu, X., Zhang, W., et al. (2023). The key roles of mycobacterium tuberculosis fadD23 C-terminal domain in catalytic mechanisms. Front. Microbiol. 14, 1090534. doi: 10.3389/fmicb.2023.1090534
Yoon, S. H., Ha, S. M., Lim, J., Kwon, S., Chun, J. (2017). A large-scale evaluation of algorithms to calculate average nucleotide identity. Antonie Van Leeuwenhoek 110, 1281–1286. doi: 10.1007/s10482-017-0844-4
Zakham, F., Sironen, T., Vapalahti, O., Kant, R. (2021). Pan and Core Genome Analysis of 183 Mycobacterium tuberculosis Strains Revealed a High Inter-Species Diversity among the Human Adapted Strains. Antibiotics (Basel) 10, 500. doi: 10.3390/antibiotics10050500
Zekic, T., Holley, G., Stoye, J. (2018). Pan-genome storage and analysis techniques. Methods Mol. Biol. 1704, 29–53. doi: 10.1007/978-1-4939-7463-4_2
Zhang, L., Lin, T. Y., Liu, W. T., Ling, F. (2024). Toward characterizing environmental sources of non-tuberculous mycobacteria (NTM) at the species level: A tutorial review of NTM phylogeny and phylogenetic classification. ACS Environ. Au. 4, 127–141. doi: 10.1021/acsenvironau.3c00074
Zhang, Z. Y., Sun, Z. Q., Wang, Z. L., Hu, H. R., Wen, Z. L., Song, Y. Z., et al. (2013). Identification and pathogenicity analysis of a novel non-tuberculous mycobacterium clinical isolate with nine-antibiotic resistance. Clin. Microbiol. Infect. 19, 91–96. doi: 10.1111/j.1469-0691.2012.03818.x
Zhao, J., Lu, B., Zhou, Y., Wang, S., Pan, X. (2022). Molecular identification and clinical significance of Mycobacterium seoulense strains from patients with nontuberculous mycobacterium infections. Diagn. Microbiol. Infect. Dis. 104, 115756. doi: 10.1016/j.diagmicrobio.2022.115756
Keywords: M. seoulense, whole genome sequencing, virulence factor, drug susceptibility, intracellular survival
Citation: Zhao J, Shen X, Jin L, Ji S and Pan X (2025) Resistance phenotypes and genomic features of Mycobacterium seoulense isolates. Front. Cell. Infect. Microbiol. 15:1553591. doi: 10.3389/fcimb.2025.1553591
Received: 31 December 2024; Accepted: 13 March 2025;
Published: 07 April 2025.
Edited by:
Mattia Pirolo, University of Copenhagen, DenmarkReviewed by:
Elias Eger, Helmholtz Institute for One Health, GermanyStacy Julisa Carrasco Aliaga, University of Milan, Italy
Copyright © 2025 Zhao, Shen, Jin, Ji and Pan. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinling Pan, cGFuZmVuZ3l1d3V6dUAxNjMuY29t
†These authors have contributed equally to this work