Description of new species of Mycobacterium terrae complex isolated from sewage at the São Paulo zoological park foundation in Brazil

Five mycobacterial isolates from sewage were classified as members of the genus Mycobacterium but presented inconclusive species assignments. Thus, the isolates (MYC017, MYC098, MYC101, MYC123 and MYC340) were analyzed by phenotypical, biochemical, matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) and genomic features to clarify their taxonomic position. Phenotypic analysis and biochemical tests did not distinguish these isolates from other non-pigmented mycobacteria. In contrast, MALDI-TOF MS analysis showed that isolates were not related to any previously described Mycobacterium species. Comparative genomic analysis showed values of ANI and dDDH between 81.59–85.56% and 24.4–28.8%, respectively, when compared to the genomes of species of this genus. In addition, two (MYC101 and MYC123) presented indistinguishable protein spectra from each other and values of ANI = 98.57% and dDDH = 97.3%, therefore being considered as belonging to the same species. Phylogenetic analysis grouped the five isolates within the Mycobacterium terrae complex (MTC) but in a specific subclade and separated from the species already described and supported by 100% bootstrap value, confirming that they are part of this complex but different from earlier described species. According to these data, we propose the description of four new species belonging to the Mycobacterium genus: (i) Mycobacterium defluvii sp. nov. strain MYC017^T (= ATCC TSD-296^T = JCM 35364^T), (ii) Mycobacterium crassicus sp. nov. strain MYC098^T (= ATCC TSD-297^T = JCM 35365T), (iii) Mycobacterium zoologicum sp. nov. strain MYC101^T (= ATCC TSD-298^T = JCM 35366^T) and MYC123 (= ATCC BAA-3216 = JCM 35367); and (iv) Mycobacterium nativiensis sp. nov. strain MYC340^T (= ATCC TSD-299^T = JCM 35368^T).

1 Introduction

Non-tuberculous mycobacteria (NTM) are species belonging to the genus Mycobacterium unrelated to the Mycobacterium tuberculosis complex and Mycobacterium leprae. Some NTM have similar characteristics and are grouped in complexes. The Mycobacterium terrae complex (MTC) was defined in 1981 by the International Working Group in Mycobacterial Taxonomy (IWGMT) and comprised two species, Mycobacterium terrae and Mycobacterium non-chromogenicum, which are non-pigmented and have intermediate growth rate (requiring 5–15 days to form visible colonies in solid media but phenotypic tests cannot clearly differentiate them (Tortoli, 2003; Tortoli et al., 2013). In the 1990s, a unique genetic signature, a two-nucleotide insertion in helix 18 of the 16S rRNA gene, was described as the most reliable marker for the attribution of mycobacterial species to the MTC (Kirschner et al., 1993; Springer et al., 1996; Tortoli et al., 2013).

Today, the MTC comprises 16 species, 13 of which were validly published,¹ are ubiquitous and considered as potentially opportunistic pathogens. They have been isolated from the environment (water and soil), clinical specimens and animals. Most species are recognized as environmental contaminants of sputum but can also cause pulmonary infections, and tenosynovitis and osteomyelitis, arising from direct inoculation injuries (Vasireddy et al., 2016, 2017; Abudaff and Beam, 2018). Mycobacterium hiberniae, “Mycobacterium icosiumassiliensis” and Mycobacterium minnesotensis have not yet been documented as responsible for human infections (Hannigan et al., 2013; Makovcova et al., 2015; Dibaj et al., 2020).

In a previous study, we assessed a large diversity of bacteria belonging to the Mycobacteriaceae family isolated from different sources of water samples in a zoo in the state of São Paulo - Brazil, as well as a variety of potential new species (Romagnoli et al., 2020). Here, we present five of these isolates as candidates for being new species belonging to MTC, after being characterized by biochemical and cultural tests, matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry (MALDI-TOF-MS) and genomic features. Our results show that the strains represent four new species of the Mycobacterium genus and contribute to the knowledge of the taxonomy of this group of bacteria.

2 Materials and methods

2.1 Bacterial strains and cultivation conditions

Strains MYC017, MYC098, MYC101, MYC123 and MYC340 were isolated from sewage at the São Paulo Zoological Park in Brazil (Romagnoli et al., 2020). Bacteria were cultivated in Middlebrook 7H10 medium supplemented with 0.5% glycerol and 10% OADC (oleic acid, albumin, dextrose, and catalase - Becton Dickinson, Franklin Lakes, USA) at 30°C for up to 10 days.

2.2 Morphology and phenotypic characteristics

The morphological characteristics of the colonies from the different isolates were observed during culturing in 7H10-OADC and after Ziehl-Neelsen (ZN) staining. Growth time, pigment production and temperature growth were evaluated at 30°C, 37°C and 42°C. In addition, classical biochemical tests were performed, including growth in a medium with thiophene 2-carboxylic acid hydrazide (TCH), growth in 5% NaCl, Tween 80 hydrolysis, potassium tellurite reduction, iron uptake, nitrate reduction, hydrolysis of urea, growth in MacConkey, and semiquantitative catalase (Kent and Kubica, 1985).

2.3 Antimicrobial susceptibility test

To perform an antimicrobial susceptibility test (AST), we used the method established by the Clinical Laboratory Standards Institute (CLSI) to determine the Minimum inhibitory concentrations (MIC) (Clinical and Laboratory Standards Institute [CLSI], 2018). From isolated colonies, we prepared a bacterial suspension with an OD of 0.5 on the McFarland scale, adjusted to the final inoculum of approximately 5 × 10⁵ CFU/mL in Muller-Hinton II cation-adjusted medium (Becton Dickinson, Franklin Lakes, USA). The suspensions were inoculated in a 96-well microplate with dilutions indexed to base 2. We tested the following drugs: amikacin, ciprofloxacin, clarithromycin, doxycycline, trimethoprim associated with sulfamethoxazole, moxifloxacin, rifampicin, isoniazid, streptomycin, and ethambutol. All tests were performed in triplicate and incubated at 30°C, and growth was observed after 3 and 7 days. Test interpretation was performed according to CLSI, document M62, which recommends for MTC evaluate the same first- and second-line drugs for M. kansasii using the same breakpoints (Woods et al., 2011; Clinical and Laboratory Standards Institute [CLSI], 2018) and published data (Cloud et al., 2006; Masaki et al., 2006; Mun et al., 2008; Sahraoui et al., 2011; Tortoli et al., 2013; Zhang et al., 2013; Djouadi et al., 2016; Vasireddy et al., 2016) (Supplementary Table 1).

2.4 DNA extraction and whole-genome sequencing

For DNA preparation, the bacterial mass was added in a tube containing 0.06 g of glass beads with 0.5 mm diameter (BioSpec, Bartlesville, USA) and 100 μL of Tris-EDTA (TE) 1X. After vigorous homogenization by vortex for 5 min, the tubes were incubated at 80°C for 10 min and centrifuged at 5000 × g for 2 min. Following, genomic DNA extraction was performed using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s recommendations (Silva et al., 2023). DNA quality was verified using Qubit^® dsDNA HS (High Sensitivity) Assay Kits (Thermo Fisher Scientific, Waltham, United States). The whole-genome sequencing (WGS) was carried out on a MiSeq instrument (Illumina San Diego, United States), using 2 × 150 paired-end chemistry and the Nextera XT library preparation kit (Illumina, San Diego, United States).

2.5 Genome features

A script developed by the authors was used to remove adaptors and primers. The read quality was verified using FastQC v0.11.9 before and after the trimming by Trimmomatic v0.35 (Bolger et al., 2014; Andrews, 2015). To verify possible genomic contamination, we used Kraken2 v2.0.8_beta with the database ftp://ftp.ccb.jhu.edu/pub/data/kraken2_dbs/minikraken_8GB_202003.tgz (Davis et al., 2013). The genomes’ assembly was performed with SPAdes v3.15.5 with k-mer sizes equal to 21, 33, 55, 77, 99, and 127 and the options error correction and carefully activated (Bankevich et al., 2012).

After assembly, the genomic sequences of each isolate were compared with sequences from species of Mycobacterium genus terrae-complex deposited in the GenBank database to determine the Average Nucleotide Identity (ANI) and the digital DNA-DNA hybridisation (dDDH) (Meier-Kolthoff et al., 2013). To determine ANI, the online tool ANI calculator² was used (Konstantinidis and Tiedje, 2005; Goris et al., 2007; Yoon et al., 2017; Tortoli et al., 2019). To calculate dDDH, we used the tool Genome-to-Genome Distance Calculator (GGDC) v2.1,³ using the values provided by the formula 2 (identities/HSP length) (Tortoli et al., 2019). To classify genomes as belonging to a particular species, we used the cut-off ANI > 95% and/or dDDH > 70%, while for classification as a subspecies, the cut-off was ANI ≤ 97% together with dDDH ≤ 80% (Meier-Kolthoff et al., 2013; Tortoli et al., 2019).

2.6 Phylogenetic analysis

We used MEGA-Molecular Evolutionary Genetics Analysis v.11.0.13 software to construct the phylogenetic trees based on partial genes. The hsp65, rpoB, and 16S rRNA gene sequences were first aligned with the CLUSTALW tool integrated into the MEGA program. The Neighbor-joining (N-J) statistical method and the Kimura 2-parameter nucleotide substitution model were used for the analyses. Bootstrap values were calculated from 1,000 replications to determine the branch significances. The support for each branch, derived from 1,000 bootstrap samples, is shown at the base of each node. The bar length (Bar 0.02) corresponds to two nucleotide substitutions per 100 nucleotides. Bootstrap values exceeding 50% were considered statistically significant.

A phylogenetic tree based on WGS was constructed that was based on, in addition to the genomes proposed as new species in the present study (MYC017, MYC098, MYC101, MYC123, and MYC340), 27 additional bacterial genomes representative of each clade of Mycobacterium genus (Supplementary Table 2). Hoyosella altamirensis was used as the outgroup, as described earlier (Tortoli et al., 2019).

All genomic sequences were annotated with Prokka v1.14.6 (employing the parameter “–kingdom Bacteria”) to standardize the protein prediction (Seemann, 2014). We used Ortho Finder v2.5.2 (Emms and Kelly, 2019) (applying the “-M msa” option) to extract the shared single copy orthologs and create the multiple alignments based on 935 genes using MAFFT v7.471 (Katoh et al., 2002). The resulting alignment was fed into RAxML-NG v1.0.1 (Kozlov et al., 2019) to generate a maximum-likelihood (ML) phylogeny, applying the “–model LG + G + F” with 500 bootstrap replicates. The tree was employed, and the midpoint was rooted in iTOL v5.0 (Letunic and Bork, 2021). The bioinformatics pipeline is described at https://doi.org/10.5281/zenodo.4660336.

2.7 Chemotaxonomy

Matrix-assisted laser Desorption/Ionization Time-of-Flight mass spectrometry (MALDI-TOF-MS) analysis was conducted as chemotaxonomy biological classification for the five mycobacterial isolates using a published protocol (El Khéchine et al., 2011). In addition, the reference strains M. fortuitum DSM46621 and M. terrae ATCC15755 (referred to as M.F. and M.T., respectively) were also subjected to MALDI-TOF-MS analysis. After extracting proteins, triplicate samples of 1 μL each were placed in a 96-well desorption plate (MSP; Bruker Daltonics). The samples were then covered with 1 μL of α-cyano-4-hydroxycinnamic acid matrix dissolved in a solution of 50% acetonitrile and 2.5% trifluoroacetic acid, allowing the samples to crystallize for analysis. A MALDI-TOF LT Microflex Bruker mass spectrometer was utilized with a 337 nm nitrogen laser in linear mode, controlled by the Flex Control v3.4 computational program (Bruker Daltonics, Billerica, United States).

The resulting spectra were saved using a laser frequency of 60HzA, with a voltage of 20 kV for ion source 1, 18.6 kV for ion source 2, and 9.28 kV for the lens. The mass range of analysis was set from 2,000 to 20,000 KDa. Multiple laser shots were performed automatically in different areas of each sample, resulting in the collection and analysis of 240 spectra. The equipment was calibrated using Escherichia coli DH5α extract. The generated spectra were analyzed with the BioTyper v10.0 (MBT) software (Bruker Daltonics, Billerica, United States). This software compares the obtained spectra from each sample with the profiles available in the database and generates a score value from 0 to 3. The score ranges from 0.000 to 1.699 when the sample is not identified, from 1.700 to 1.999 when there is a probable identification of the genus, from 2.000 to 2.299, indicating a secure identification of the genus and potential species, and finally from 2.300 to 3.000 representing a secure identification of the species.

3 Results

Initial characterization of the isolates was carried out by sequencing and analyzing the hsp65, rpoB and 16S rRNA genes in addition to phenotypic and biochemical characteristics. The phylogenetic trees clearly separated the newly proposed species from the others already described, uniting them in a clade closely related to MTC (Figures 1A–C, respectively). Furthermore, it is shown that MYC017 is very closely related to MYC101 and MYC123, while MYC098 is related to MYC340 by analysis of the rpoB and 16S rRNA and related to Mycobacterium engbaekii by the hsp65 gene. In addition, all isolates presented the insertion of two nucleotides in helix 18 of the 16S rRNA gene, in comparison with other slow-growing mycobacteria, as well as species of MTC (Supplementary Figure 1).

FIGURE 1

Figure 1. Phylogenetic analysis of representative species of the Mycobacterium terrae complex and the MYC strains (17, 98, 101, 12, 340), using neighbour-joining statistical methods and the Kimura 2-parameter nucleotide substitution model methods, with the following targets/genes/conserved regions: (A) 16S rRNA gene; (B) hsp65 gene; and (C) rpoB gene. (A,D) comparative analysis was performed using complete genome sequencing data from the MYC strains and available genomes in Genbank. Each panel is a representative result among four experiments.

Microscopic examination of isolated colonies after ZN staining showed that all were acid-fast bacilli. All isolates grew in seven days at 30 and 37°C, and were positive for the semiquantitative catalase, TCH growth and nitrate reduction tests, but negative for the MacConkey growth tests, NaCl growth, iron uptake, and growth at 42°C. Variable results were observed for urease tests, Tween 80 hydrolysis and tellurite reduction, in which only the MYC340 isolate was positive for all these tests. Furthermore, none of the isolates produced pigment and were therefore classified as non-chromogenic. The phenotypic and biochemical tests for the new bacterial species are described in Table 1. Variable results from phenotypic tests are also presented by the species belonging to the MTC as well as by the isolates in this study (Supplementary Table 3).

TABLE 1

Table 1. Phenotypic and biochemical results performed for five isolates proposed as new species.

The AST results based on MIC are described in Table 2 which shows that four isolates are susceptible to amikacin, ciprofloxacin, clarithromycin, sulfamethoxazole-trimethoprim, moxifloxacin, and resistant to rifampicin, isoniazid, streptomycin and ethambutol. Isolates MYC017, MYC101, and MYC123 are susceptible to doxycycline, while MYC340 is resistant.

TABLE 2

Table 2. Antimicrobial susceptibility results based on Minimum Inhibitory Concentration (MIC) (μg/mL) for isolates in this study.

The genomes were sequenced, and assembly and annotation were performed (Table 3). Quality control based on taxonomic classification is presented in Supplementary Figure 2. The results of the comparison of ANI and dDDH among the five isolates and their comparison with species of MTC are shown in Supplementary Table 4. Isolates MYC101 and MYC123 presented values for ANI = 98.57 and dDDH = 87.30, therefore belonging to the same species. The other isolates presented values for ANI and dDDH below the cut-off point (<95% for ANI and <70% for dDDH) and, as such, represent different species. Comparative analysis of ANI and dDDH between the genomes of the five isolates in this study and the species of the MTC with available genomes showed values below the established cut-off point, thus configuring four new species of this genus.

TABLE 3

Table 3. Sequencing, assembly and annotation results for five isolates proposed as new species of the Mycobacterium terrae complex.

The ML phylogenetic tree showed that isolates of new species were grouped within the MTC, but in a specific subclade separated from the species already described and supported by 100% bootstrap value, confirming that they are part of the Mycobacterium genus within the terrae clade and different from already established species (Figure 1D).

The five isolates and reference strains of Mycobacterium fortuitum and M. terrae were subjected to MALDI-TOF-MS analysis. The obtained spectra were then compared to those in the database using the Biotyper software, resulting in a value score for each. The reference strains, M. fortuitum and M. terrae were correctly identified with scores of 2.2 and 2.0, respectively, while the analyzed MYC isolates exhibited spectra of good quality but had identification scores ranging from 1.4 to 1.6, therefore not being identified (Figure 2). Additionally, it is noteworthy that MYC101 and MYC123, identified as the same species through genomic analyses, displayed indistinguishable profiles, confirming their similar nature.

FIGURE 2

Figure 2. Comparative Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS) of MYC isolates (17, 98, 101, 123, 340) and reference strains Mycobacterium terrae (M.T) and Mycobacterium fortuitum (M.F) identified different mass-to-charge (m/z) peak patterns, indicating potential novel species within the MYC group.

Supplementary Table 5 details the correlation between culture isolates, species nomenclature, and genomic identifications from biorepositories.

4 Discussion

Description of novel Prokatiotic taxa might have important implications for human health and biotechnological implications. In this sense, studies are essential for advancing our knowledge of microbial diversity, which has implications for fields as diverse as medicine, ecology, industry and environmental sciences. In our previous study, a large sampling of mycobacteria from water samples from the São Paulo Zoological Park was evaluated and the results suggested the presence of a putative new species (Romagnoli et al., 2020). From these, five isolates were selected for further characterization regarding their taxonomic position.

Currently, the Mycobacterium genus is formed up of 196 species (number of child taxa with a validly published and correct name) (LPSN, 2023), the majority of which are found in natural environments such as lakes, rivers, swamps, and soil, as well as water treatment and distribution systems (Falkinham, 2015; Claeys and Robinson, 2018). Phenotypic, biochemical and molecular characteristics are useful for differentiating between species. Growth time is particularly useful for dividing species into three large groups: fast-growing (visible colonies within 7 days), slow-growing (growing of visible colonies after 7 days) and intermediate-growing (growing of visible colonies from 5 to 15 days) organisms. In this work, biochemical and cultural tests were not able to clearly distinguish the five isolates from other non-pigmented mycobacteria but could be classified by growth time as intermediate, characteristic compatible with species of the MTC.

The antimicrobial susceptibility profile presented by the isolates in this study was also observed in other species of the M. terrae complex (Hannigan et al., 2013; Tortoli et al., 2013; Djouadi et al., 2016). It is important to note that the susceptibility variability has also been reported for several isolates of M. arupensis (Cloud et al., 2006). Thus, it is not possible to state that the results observed in the isolates in this study represent the profile of each new species.

Analysis of three housekeeping genes showed that all isolates are related to the MTC and is according to the intermediate growth characteristic. In addition, the isolates have an insertion of two nucleotides in the helix 18 of the 16S rRNA gene. It is consense that this simple two-nucleotide insertion remains one of the most important features to determine whether a mycobacterium is related to the MTC (Tortoli et al., 2013; Abudaff and Beam, 2018).

For further characterization of these isolates, their genomes were sequenced for taxonomic analyses using ANI and dDDH. Genome sequences have gained extensive usage and pertinent software tools are easily accessible both as web services and as standalone applications. Furthermore, commonly adopted threshold values for ANI and dDDH to recognize species are 95–96 and 70%, respectively (Richter and Rosselló-Móra, 2009; Meier-Kolthoff et al., 2013; Chun et al., 2018). Our study showed values of ANI and dDDH between 81.59–85.56 and 24.4–28.8%, respectively, when compared to the genomes of species of the canonical Mycobacterium genus, confirming the classification of these isolates as new species.

Phylogenetic analysis used representative members from the five monophyletic groups, as previously suggested (Gupta et al., 2018), to construct the ML phylogenetic tree. The study was based on genomic sequences on molecular signatures identified in core proteins, dividing the genus Mycobacterium into five monophyletic groups, namely: Mycolicibacterium (Fortuitum-Vaccae clade), Mycolicibacter (Terrae clade), Mycolicibacillus (Triviale clade), Mycobacteroides (Chelonae-Abscessus clade) and Mycobacterium (Tuberculosis-Simiae clade). Thus, representative members from each of the five monophyletic groups were selected for analysis which showed that MYC isolates form a distinct cluster within the MTC (”Terrae Clade” or Mycolicibacter) and exhibit differences from previously established species. However, it is important to mention that the canonical genus Mycobacterium has been reconstituted and was therefore adopted here (Meehan et al., 2021).

Demonstrating whether a strain represents a new species or not can be adequately achieved with a genome sequence alone, without the obligatory generation of various types of chemotaxonomic data. Thus, for the proposal of new species, a comprehensive analysis of genome sequences, a basic assessment of growth characteristics, and the public deposition of type strains should be performed (Vandamme and Sutcliffe, 2021). We therefore have chosen this polyphasic analysis of characteristics such as morphology and growth conditions, physiology, genome and chemotaxonomy. In this context, it is worth highlighting that among the diverse chemotaxonomic information available, MALDI-TOF mass spectrometry analysis has demonstrated usefulness not only for routine diagnostic microbiology but also contributes to taxonomy as the detection of new protein profiles and can therefore readily suggest new species (Rodriguez-Temporal et al., 2017; Alcaide et al., 2018; López Medrano et al., 2021; Oviaño and Rodríguez-Sánchez, 2021; Li et al., 2022). The availability of freely accessible and continuously updated spectral databases would undoubtedly be of great benefit to the systematics community and would further promote the adoption of this technology (Vandamme and Sutcliffe, 2021). In this study, the MALDI-TOF results align with other data indicating that the MYC isolates represent previously unidentified species and also support that isolates MYC101 and MYC123 belong to the same species.

In conclusion, we present evidence to classify the strains isolated (MYC017, MYC098, MYC101, MYC123 and MYC340) from waste samples in a zoo in the state of São Paulo – Brazil as four new species of the genus Mycobacterium, considering MYC101 and MYC123 a single species.

5 Nomenclature

5.1 Description of Mycobacterium vasticus sp. nov

Mycobacterium vasticus (vas.ti.cus. L. masc. adj. vastum). Waste and sewage are characteristics of the origin of this isolate. Growth on solid Middlebrook 7H10-OADC medium and Löwenstein–Jensen medium requires 7 days at 30°C and 37°C, giving rise to shiny, smooth, non-chromogenic colonies with defined edges and creamy consistency. Growth in TCH was observed, neither in NaCl nor in MacConkey. M. vasticus is positive for catalase production and nitrate and tellurite reduction, while it is negative for Tween 80 hydrolysis, urease production and iron uptake. AST revealed susceptibility to amikacin, ciprofloxacin, clarithromycin, sulfamethoxazole-trimethoprim, moxifloxacin, and doxycycline, but resistance to ethambutol, rifampicin, isoniazid, and streptomycin. The type strain is MYC017^T (= ATCC TSD-296^T = JCM 35364^T). The GenBank accession numbers for the draft genome, Biosample and 16S rRNA gene sequence are, respectively, CP084028, SAMN20959233 and MK890459.1.

5.2 Description of Mycobacterium crassicus sp. nov

Mycobacterium crassicus (cras.si.cus. L. masc. adjt. crassa). Crude sewage is characteristic of the environment from where the isolate came. Growth on solid Middlebrook 7H10-OADC and Löwenstein–Jensen medium requires 7 days at 30°C and 37°C, giving rise to shiny, smooth, non-chromogenic colonies with defined edges and creamy consistency. Growth in TCH, but not in NaCl or MacConkey was observed. M. crassicus is positive for catalase production, urease production, Tween 80 hydrolysis and nitrate reduction while negative for tellurite reduction and iron uptake. The type strain is MYC098^T (= ATCC TSD-297^T = JCM 35365^T). The GenBank accession numbers for the draft genome, Biosample and 16S rRNA gene are, respectively: CP084029, SAMN20959234, and MK890478.1.

5.3 Description of Mycobacterium zoologicum sp. nov

Mycobacterium zoologicum (zoo.lo.gi.cum. L. masc. subst. zoo.) was isolated from a zoo’s sewage treatment plant. Growth in solid Middlebrook 7H10 medium and on Löwenstein–Jensen medium requires 7 days at 30°C and 37°C, giving rise to non-chromogenic, shiny, smooth colonies with defined edges and creamy consistency. The specie presents growth in TCH but not in NaCl or MacConkey. M. zoologicum is positive for Tween 80 hydrolysis, catalase production, nitrate, and tellurite reduction, while it is negative for urease production and iron uptake. They are susceptible to amikacin, ciprofloxacin, clarithromycin, sulfamethoxazole-trimethoprim, moxifloxacin, and doxycycline but resistant to ethambutol, rifampicin, isoniazid, and streptomycin. The type strain is MYC101^T (= ATCC TSD-298^T = JCM 35366^T) and an additional isolate MYC123 (= ATCC BAA-3216 = JCM 35367). The GenBank accession numbers for the draft genome, Biosample and 16S rRNA gene sequence are, respectively: MYC101^T (CP084030, SAMN20959235, MK890479.1) and MYC123 (CP083985, SAMN20062777, MK890481.1).

5.4 Description of Mycobacterium nativiensis sp. nov

Mycobacterium nativiensis (na.ti.vi.en.sis. L. masc. subst. nativus.). From “native”, local where bacteria was isolated. Growth in solid Middlebrook 7H10 medium and on Löwenstein–Jensen medium was observed in 7 days at 30°C and 37°C, giving rise to non-chromogenic, shiny, smooth colonies with defined edges and creamy consistency. It shows growth in TCH but not in NaCl or MacConkey. M. nativiensis is positive for Tween 80 hydrolysis, catalase and urease production, and nitrate and tellurite reduction while negative for iron uptake. It is susceptible to amikacin, ciprofloxacin, clarithromycin, sulfamethoxazole-trimethoprim and moxifloxacin but resistant to ethambutol, rifampicin, isoniazid, doxycycline, and streptomycin. The type strain is MYC340^T (= ATCC TSD-299^T = JCM 35368^T). The GenBank accession numbers for the draft genome, Biosample and 16S rRNA gene sequence are, respectively: CP083986, SAMN20062778 and MK890521.1.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Author contributions

CLR: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review and editing. ECC: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing. EM: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing. LB: Investigation, Methodology, Writing – review and editing. AS: Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – original draft, Writing – review and editing. NS: Investigation, Methodology, Writing – review and editing. LM: Investigation, Methodology, Writing – review and editing. MJ: Investigation, Methodology, Writing – review and editing. ML: Data curation, Formal analysis, Investigation, Methodology, Writing – review and editing, Resources, Validation. LD: Writing – review and editing, Data curation, Formal analysis, Investigation, Methodology, Software. PS: Funding acquisition, Resources, Supervision, Writing – review and editing. SL: Conceptualization, Funding acquisition, Supervision, Writing – original draft, Writing – review and editing. CV-N: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Visualization, Writing – original draft, Writing – review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was provided by research grants 2019/14137-4 from Fundação de Apoio à Pesquisa do Estado de São Paulo (FAPESP) and 406710/2018-0 from Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). NMS and CLR were supported by a fellowship from the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) – Finance Code 001.

Acknowledgments

We would like to thank Francisco das Chagas de Sousa, Karolina Rosa Fernandes Beraldo, and Mayla Santos Rodrigues for their valuable technical support in carrying out AST, MALDI-TOF, and English editing tasks, respectively. We would also like to thank the Bioinformatic platform RPT04A (Rede de Plataformas Tecnológicas FIOCRUZ) for their support.

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.

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/fmicb.2024.1335985/full#supplementary-material

Footnotes

1. ^ https://lpsn.dsmz.de/
2. ^ https://www.ezbiocloud.net/tools/ani
3. ^ http://ggdc.dsmz.de/ggdc.php

References

Abudaff, N. N., and Beam, E. (2018). Mycobacterium arupense: A review article on an emerging potential pathogen in the Mycobacterium terrae complex. J. Clin. Tuberc. Mycobact. Dis. 10, 1–5. doi: 10.1016/j.jctube.2017.11.001

PubMed Abstract | Crossref Full Text | Google Scholar

Alcaide, F., Amlerová, J., Bou, G., Ceyssens, P. J., Coll, P., Corcoran, D., et al. (2018). How to: identify non-tuberculous Mycobacterium species using MALDI-TOF mass spectrometry. Clin. Microbiol. Infect. 24, 599–603. doi: 10.1016/j.cmi.2017.11.012

PubMed Abstract | Crossref Full Text | Google Scholar

Andrews, S. (2015). FASTQC A quality control tool for high throughput sequence data. Cambridge: Babraham Institute.

Google Scholar

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021

PubMed Abstract | Crossref Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | Crossref Full Text | Google Scholar

Chun, J., Oren, A., Ventosa, A., Christensen, H., Arahal, D. R., da Costa, M. S., et al. (2018). Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int. J. Syst. Evol. Microbiol. 68, 461–466. doi: 10.1099/ijsem.0.002516

PubMed Abstract | Crossref Full Text | Google Scholar

Claeys, T. A., and Robinson, R. T. (2018). The many lives of nontuberculous mycobacteria. J. Bacteriol. 200, e00739–17. doi: 10.1128/JB.00739-17

PubMed Abstract | Crossref Full Text | Google Scholar

Cloud, J. L., Meyer, J. J., Pounder, J. I., Jost, K. C., Sweeney, A., Carroll, K. C., et al. (2006). Mycobacterium arupense sp. nov., a non-chromogenic bacterium isolated from clinical specimens. Int. J. Syst. Evol. Microbiol. 56, 1413–1418. doi: 10.1099/ijs.0.64194-0

PubMed Abstract | Crossref Full Text | Google Scholar

Clinical and Laboratory Standards Institute [CLSI]. (2018). M24Ed3 | Susceptibility Testing of Mycobacteria, Nocardia spp., and Other Aerobic Actinomycetes, 3rd Edn. Available online at: https://clsi.org/standards/products/microbiology/documents/m24/#:~:text=CLSI (accessed October 22, 2023).

Google Scholar

Davis, M. P. A., van Dongen, S., Abreu-Goodger, C., Bartonicek, N., and Enright, A. J. (2013). Kraken: A set of tools for quality control and analysis of high-throughput sequence data. Methods San Diego Calif. 63, 41–49. doi: 10.1016/j.ymeth.2013.06.027

PubMed Abstract | Crossref Full Text | Google Scholar

Dibaj, R., Shojaei, H., and Narimani, T. (2020). Identification and molecular characterization of mycobacteria isolated from animal sources in a developing country. Acta Trop. 204:105297. doi: 10.1016/j.actatropica.2019.105297

PubMed Abstract | Crossref Full Text | Google Scholar

Djouadi, L. N., Levasseur, A., Khalil, J. B., Blanc-Taileur, C., Asmar, S., Ghiloubi, W., et al. (2016). Mycobacterium icosiumassiliensis sp. nov., a new member in the mycobacterium terrae complex isolated from surface water in Algeria. Curr. Microbiol. 73, 255–264. doi: 10.1007/s00284-016-1062-4

PubMed Abstract | Crossref Full Text | Google Scholar

El Khéchine, A., Couderc, C., Flaudrops, C., Raoult, D., and Drancourt, M. (2011). Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry identification of mycobacteria in routine clinical practice. PLoS One 6:e24720.

Google Scholar

Emms, D. M., and Kelly, S. (2019). OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 20:238. doi: 10.1186/s13059-019-1832-y

PubMed Abstract | Crossref Full Text | Google Scholar

Falkinham, J. O. (2015). Environmental sources of nontuberculous mycobacteria. Clin. Chest Med. 36, 35–41. doi: 10.1016/j.ccm.2014.10.003

PubMed Abstract | Crossref Full Text | Google Scholar

Goris, J., Konstantinidis, K. T., Klappenbach, J. A., Coenye, T., Vandamme, P., and Tiedje, J. M. (2007). DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int. J. Syst. Evol. Microbiol. 57, 81–91. doi: 10.1099/ijs.0.64483-0

PubMed Abstract | Crossref Full Text | Google Scholar

Gupta, R. S., Lo, B., and Son, J. (2018). Phylogenomics and comparative genomic studies robustly support division of the genus Mycobacterium into an emended genus Mycobacterium and four novel genera. Front. Microbiol. 9:67. doi: 10.3389/fmicb.2018.00067

PubMed Abstract | Crossref Full Text | Google Scholar

Hannigan, G. D., Krivogorsky, B., Fordice, D., Welch, J. B., and Dahl, J. L. (2013). Mycobacterium minnesotense sp. nov., a photochromogenic bacterium isolated from sphagnum peat bogs. Int. J. Syst. Evol. Microbiol. 63, 124–128. doi: 10.1099/ijs.0.037291-0

PubMed Abstract | Crossref Full Text | Google Scholar

Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 30, 3059–3066. doi: 10.1093/nar/gkf436

PubMed Abstract | Crossref Full Text | Google Scholar

Kent, P. T., and Kubica, G. P. (1985). Public Health Mycobacteriology: A Guide for the Level III Laboratory. Atlanta, GA: Centers for Disease Control.

Google Scholar

Kirschner, P., Springer, B., Vogel, U., Meier, A., Wrede, A., Kiekenbeck, M., et al. (1993). Genotypic identification of mycobacteria by nucleic acid sequence determination: Report of a 2-year experience in a clinical laboratory. J. Clin. Microbiol. 31, 2882–2889.

Google Scholar

Konstantinidis, K. T., and Tiedje, J. M. (2005). Genomic insights that advance the species definition for prokaryotes. Proc. Natl. Acad. Sci. U.S.A. 102, 2567–2572. doi: 10.1073/pnas.0409727102

PubMed Abstract | Crossref Full Text | Google Scholar

Letunic, I., and Bork, P. (2021). Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, W293–W296. doi: 10.1093/nar/gkab301

PubMed Abstract | Crossref Full Text | Google Scholar

Li, B., Zhu, C., Sun, L., Dong, H., Sun, Y., Cao, S., et al. (2022). Performance evaluation and clinical validation of optimized nucleotide MALDI-TOF-MS for mycobacterial identification. Front. Cell. Infect. Microbiol. 12:1079184. doi: 10.3389/fcimb.2022.1079184

PubMed Abstract | Crossref Full Text | Google Scholar

López Medrano, R., Burgos Asurmendi, I., and Rivero Lezcano, O. M. (2021). A rapid proteomic system (MALDI-TOF) for nontuberculous-mycobacteria identification. Enfermedades Infecc. Microbiol. Clin. 40, 385–387. doi: 10.1016/j.eimc.2021.06.007

PubMed Abstract | Crossref Full Text | Google Scholar

LPSN. (2023). Genus Mycobacterium - List of Prokaryotic names with Standing in Nomenclature (LPSN). Available online at: https://lpsn.dsmz.de/genus/mycobacterium (accessed October 23, 2023).

Google Scholar

Makovcova, J., Babak, V., Slany, M., and Slana, I. (2015). Comparison of methods for the isolation of mycobacteria from water treatment plant sludge. Antonie Van Leeuwenhoek 107, 1165–1179. doi: 10.1007/s10482-015-0408-4

PubMed Abstract | Crossref Full Text | Google Scholar

Masaki, T., Ohkusu, K., Hata, H., Fujiwara, N., Iihara, H., Yamada-Noda, M., et al. (2006). Mycobacterium kumamotonense Sp. Nov. recovered from clinical specimen and the first isolation report of Mycobacterium arupense in Japan: Novel slowly growing, nonchromogenic clinical isolates related to Mycobacterium terrae complex. Microbiol. Immunol. 50, 889–897. doi: 10.1111/j.1348-0421.2006.tb03865.x

PubMed Abstract | Crossref Full Text | Google Scholar

Meehan, C. J., Barco, R. A., Loh, Y.-H. E., Cogneau, S., and Rigouts, L. (2021). Reconstituting the genus Mycobacterium. Int. J. Syst. Evol. Microbiol. 71:004922. doi: 10.1099/ijsem.0.004922

PubMed Abstract | Crossref Full Text | Google Scholar

Meier-Kolthoff, J. P., Auch, A. F., Klenk, H.-P., and Göker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 14:60. doi: 10.1186/1471-2105-14-60

PubMed Abstract | Crossref Full Text | Google Scholar

Mun, H.-S., Park, J.-H., Kim, H., Yu, H.-K., Park, Y.-G., Cha, C.-Y., et al. (2008). Mycobacterium senuense sp. nov., a slowly growing, non-chromogenic species closely related to the Mycobacterium terrae complex. Int. J. Syst. Evol. Microbiol. 58, 641–646. doi: 10.1099/ijs.0.65374-0

PubMed Abstract | Crossref Full Text | Google Scholar

Oviaño, M., and Rodríguez-Sánchez, B. (2021). MALDI-TOF mass spectrometry in the 21st century clinical microbiology laboratory. Enfermedades Infecc. Microbiol. Clin. Engl. Ed 39, 192–200. doi: 10.1016/j.eimc.2020.02.027

PubMed Abstract | Crossref Full Text | Google Scholar

Kozlov, A. M., Darriba, D., Flouri, T., Morel, B., and Stamatakis, A. (2019). RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference | Bioinformatics | Oxford Academic. Available online at: https://academic.oup.com/bioinformatics/article/35/21/4453/5487384 (accessed July 13, 2023).

Google Scholar

Richter, M., and Rosselló-Móra, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. U.S.A. 106, 19126–19131. doi: 10.1073/pnas.0906412106

PubMed Abstract | Crossref Full Text | Google Scholar

Rodriguez-Temporal, D., Perez-Risco, D., Struzka, E. A., Mas, M., and Alcaide, F. (2017). Impact of updating the MALDI-TOF MS database on the identification of nontuberculous mycobacteria. J. Mass Spectrom. 52, 597–602. doi: 10.1002/jms.3944

PubMed Abstract | Crossref Full Text | Google Scholar

Romagnoli, C. L., Pellegrino, K. C. M., Silva, N. M., Brianesi, U. A., Leão, S. C., Rabello, M. C., et al. (2020). Diversity of Mycobacteriaceae from aquatic environment at the São Paulo Zoological Park Foundation in Brazil. PLoS One 15:e0227759. doi: 10.1371/journal.pone.0227759

PubMed Abstract | Crossref Full Text | Google Scholar

Sahraoui, N., Ballif, M., Zelleg, S., Yousfi, N., Ritter, C., Friedel, U., et al. (2011). Mycobacterium algericum sp. nov., a novel rapidly growing species related to the Mycobacterium terrae complex and associated with goat lung lesions. Int. J. Syst. Evol. Microbiol. 61, 1870–1874. doi: 10.1099/ijs.0.024851-0

PubMed Abstract | Crossref Full Text | Google Scholar

Seemann, T. (2014). Prokka: Rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

PubMed Abstract | Crossref Full Text | Google Scholar

Silva, N. M., Romagnoli, C. L., Santiago, C. R., do, N., de Lacerda, J. P. A., Leão, S. C., et al. (2023). Multi-approach characterization of novel pyrene-degrading mycolicibacterium austroafricanum isolates lacking nid genes. Microorganisms 11:1413. doi: 10.3390/microorganisms11061413

PubMed Abstract | Crossref Full Text | Google Scholar

Springer, B., Stockman, L., Teschner, K., Roberts, G. D., and Bottger, E. C. (1996). Two-laboratory collaborative study on identification of mycobacteria: Molecular versus phenotypic methods. J. Clin. Microbiol. 34, 296–303. doi: 10.1128/jcm.34.2.296-303.1996

PubMed Abstract | Crossref Full Text | Google Scholar

Tortoli, E. (2003). Impact of genotypic studies on mycobacterial taxonomy: The new mycobacteria of the 1990s. Clin. Microbiol. Rev. 16, 319–354. doi: 10.1128/CMR.16.2.319-354.2003

PubMed Abstract | Crossref Full Text | Google Scholar

Tortoli, E., Gitti, Z., Klenk, H.-P., Lauria, S., Mannino, R., Mantegani, P., et al. (2013). Survey of 150 strains belonging to the Mycobacterium terrae complex and description of Mycobacterium engbaekii sp. nov., Mycobacterium heraklionense sp. nov. and Mycobacterium longobardum sp. nov. Int. J. Syst. Evol. Microbiol. 63, 401–411. doi: 10.1099/ijs.0.038737-0

PubMed Abstract | Crossref Full Text | Google Scholar

Tortoli, E., Meehan, C. J., Grottola, A., Fregni Serpini, G., Fabio, A., Trovato, A., et al. (2019). Genome-based taxonomic revision detects a number of synonymous taxa in the genus Mycobacterium. Infect. Genet. Evol. 75:103983. doi: 10.1016/j.meegid.2019.103983

PubMed Abstract | Crossref Full Text | Google Scholar

Vandamme, P., and Sutcliffe, I. (2021). Out with the old and in with the new: Time to rethink twentieth century chemotaxonomic practices in bacterial taxonomy. Int. J. Syst. Evol. Microbiol. 71, 005127. doi: 10.1099/ijsem.0.005127

PubMed Abstract | Crossref Full Text | Google Scholar

Vasireddy, R., Vasireddy, S., Brown-Elliott, B. A., Wengenack, N. L., Eke, U. A., Benwill, J. L., et al. (2016). Mycobacterium arupense, Mycobacterium heraklionense, and a newly proposed species, “Mycobacterium virginiense” sp. nov., but Not Mycobacterium nonchromogenicum, as species of the mycobacterium terrae complex causing tenosynovitis and osteomyelitis. J. Clin. Microbiol. 54, 1340–1351. doi: 10.1128/JCM.00198-16

PubMed Abstract | Crossref Full Text | Google Scholar

Vasireddy, R., Vasireddy, S., Brown-Elliott, B. A., Wengenack, N. L., Eke, U. A., Benwill, J. L., et al. (2017). Correction for Vasireddy et al., Mycobacterium arupense, Mycobacterium heraklionense, and a Newly Proposed Species, “Mycobacterium virginiense” sp. nov., but Not Mycobacterium nonchromogenicum, as Species of the Mycobacterium terrae Complex Causing Tenosynovitis and Osteomyelitis. J. Clin. Microbiol. 55:985. doi: 10.1128/JCM.02290-16

PubMed Abstract | Crossref Full Text | Google Scholar

Woods, G. L., Brown-Elliott, B. A., Conville, P. S., Desmond, E. P., Hall, G. S., Lin, G., et al. (2011). Susceptibility testing of mycobacteria, nocardiae, and other aerobic actinomycetes, 2nd Edn. Wayne, PA: Clinical and Laboratory Standards Institute.

Google Scholar

Yoon, S.-H., Ha, S.-M., Lim, J., Kwon, S., and 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

PubMed Abstract | Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar