Phylogenetic trees were constructed for 150 members of the genus Mycobacterium whose genomes are now sequenced (some characteristics of these genomes are listed in Supplementary Table 1) and six members from the order Corynebacteriales (viz. Corynebacterium diphtheriae NCTC 11397, Gordonia bronchialis DSM 43247, Nocardia farcinica NCTC 11134, Rhodococcus erythropolis PR4, Segniliparus rotundus DSM 44985 and Tsukamurella paurometabola DSM 20162), which served as outgroups. The first of these trees was based on 1941 core proteins from the genomes of Mycobacterium species and its construction was carried out by using a software pipeline (Adeolu et al., 2016). Briefly, the CD-HIT program was used (Li and Godzik, 2006; Fu et al., 2012) to identify protein families sharing a minimum of 50% in sequence identity and sequence length and which were found in at least 80% of the input genomes. The Clustal Omega (Sievers et al., 2011) algorithm was used to generate multiple sequence alignment (MSA) of these protein families. The aligned protein families were trimmed with TrimAl (Capella-Gutierrez et al., 2009) to remove poorly aligned regions (Talavera and Castresana, 2007) before concatenation to the other core proteins. The concatenated sequence alignment of 1941 core proteins consisted of 624,360 aligned amino acids. Another comprehensive phylogenetic tree was constructed based on concatenated sequences for 136 proteins, which comprise the phyloeco markers set for the phylum Actinobacteria (Wang and Wu, 2013). Information regarding these proteins is provided in Supplementary Table 2. The profile Hidden Markov Models of these protein families were used for the identification of members of these protein families in the input genomes using HMMer 3.1 (Eddy, 2011). The sequence alignments were trimmed using TrimAl (Capella-Gutierrez et al., 2009) before their concatenation into a single file. The combined sequence from the phyloeco set of proteins consisted of a total of 44,976 aligned amino acids. Maximum likelihood (ML) trees based on both these sequence alignments were constructed using the Whelan and Goldman model of protein sequence evolution (Whelan and Goldman, 2001) in FastTree 2 (Price et al., 2010) and the Le and Gascuel model of protein sequence evolution (Le and Gascuel, 2008) in RAxML 8 (Stamatakis, 2014). Optimization of the robustness of the tree was completed by conducting SH tests (Guindon et al., 2010) in RAxML 8 (Stamatakis, 2014). The identification of the conserved protein families and the construction of phylogenetic trees were completed using an internal software pipeline (Adeolu et al., 2016).

In addition to these two comprehensive trees, another phylogenetic tree was constructed based on concatenated sequences for 8 conserved housekeeping proteins (viz. RpoA, RpoB, RpoC, GyrA, GyrB, Hsp65, EF-Tu and RecA). After removal of non-conserved regions, the concatenated sequence alignment in this case consisted of 6052 aligned amino acids. A maximum likelihood phylogenetic tree based on this sequence was constructed as described above.

The sequence alignments of the 1941 core proteins identified by the above methods were also used to measure genome relatedness. Using the amino acid sequences from these conserved protein families, the amino acid sequence identity between each pair of Mycobacterium genomes was calculated (Thompson et al., 2013).

Information regarding branching of all type species from the genus Mycobacterium in a tree based on 16S rRNA sequences was obtained from the SILVA All Species Tree of Life Project 128 (Quast et al., 2013).

The identification of CSIs was carried out as described in earlier work (Gao and Gupta, 2005; Bhandari et al., 2012; Gupta, 2014; Naushad et al., 2014; Sawana et al., 2014). All annotated proteins from the genomes of M. tuberculosis H37Rv and M. sinense JDM601 were used in these analyses. BLASTp (Altschul et al., 1997) searches were conducted on all protein sequences >100 amino acids in length against the NCBI non-redundant (nr) database. Multiple sequence alignments were generated by obtaining 15–25 homologs from diverse Mycobacterium species and 8–10 homologs from other groups of bacteria. The alignments were visually inspected for sequence gaps of fixed lengths which were flanked on both sides by at least 5 conserved amino acids in the neighboring 30–40 amino acids, and appeared to be shared by either some or all mycobacterial homologs. Query sequences encompassing the potential indel and flanking regions (60–100 amino acids long) were collected and subjected to a more detailed BLASTp search (500 or more hits) to determine the group specificities of the observed indels. Signature files for all CSIs of interest were created using SIG_CREATE and SIG_STYLE programs in the GLEANS software package (available on Gleans.net). Unless otherwise noted, the described CSIs are specific for the indicated groups of species.

The identification of conserved signature proteins was carried out using the protocol described in earlier work (Gao et al., 2006; Adeolu and Gupta, 2014; Naushad et al., 2014). BLASTp (Altschul et al., 1997) searches were conducted on all sequenced proteins from the genomes of M. tuberculosis H37Rv, M. aurum (LSHTM), M. sinense JDM601 (Zhang et al., 2011), M. triviale DSM 44153 (Fedrizzi et al., 2017), and M. abscessus ATCC 19977 (Ripoll et al., 2009) against the NCBI nr database. Proteins of interest were those where either all significant hits were limited to the genus Mycobacterium or the indicated groups/clades of mycobacteria, or where a large increase in E value was observed from the last hit belonging to these groups and the first hit from any other bacteria, and the E-values for the latter hits were >1e⁻³ (Gao et al., 2006; Gao and Gupta, 2012; Naushad et al., 2014). However, in some cases, a few proteins where an isolated significant hit from an unrelated group of bacteria was observed were also retained as CSPs specific for the group of interest.

In the present work, two comprehensive phylogenomic trees were constructed based on the genome sequences of 150 Mycobacterium species. The first of these trees was a core genome tree of 1941 proteins, whose homologs are present in at least 80% of the input mycobacterial genomes as well as the outgroup species. The second genome sequence tree was based on 136 proteins, which are part of the phyloeco set for the phylum Actinobacteria. The trimmed concatenated sequence alignments for the two sets of core proteins, which were employed for phylogenetic analyses, consisted of 624,360 and 44,976 aligned amino acids, respectively. Although phylogenetic trees based on core genes/proteins for mycobacterial species have also been constructed in earlier studies (Prasanna and Mehra, 2013; Fedrizzi et al., 2017; Wee et al., 2017), they were based only on a small number (between 28 and 47) of Mycobacterium species. In contrast, the trees produced in this work include information for ~80% (150 of the 188) of all known mycobacterial species and thus constitute the most comprehensive phylogenetic trees constructed for the genus Mycobacterium. In addition to the two core genome protein trees, a maximum-likelihood tree was also constructed based on concatenated sequences of 8 conserved housekeeping proteins.

The ML trees based on the core proteins from mycobacterial genomes and for the phylum Actinobacteria are shown in Figures 1A,B, respectively. The tree based on the 8 conserved proteins is provided as Supplementary Figure 1. In all of these phylogenetic trees, which were rooted using the sequences from the Corynebacteriales species, nearly all of the observed nodes were supported with high (100%) bootstrap scores or SH-values. Further, the majority of the interrelationships among the Mycobacterium species were highly similar and consistent in all constructed trees. In all of these trees, members of the genus Mycobacterium consistently grouped into four main clades and a clade consisting of the M. triviale—M. koreense, as indicated in Figure 1. Three of these clades are comprised of the slow-growing species, whereas the other two clades are mostly made up of the fast-growing species. Of the two clades of fast-growing species, the first clade referred to as the “Abscessus-Chelonae” clade, forms the earliest branching lineage within the genus Mycobacterium. The second clade of the fast-growing species designated as the “Fortuitum-Vaccae” clade, encompasses most of the other fast-growing species including those related to M. fortuitum, M. vaccae, M. parafortuitum, and M. mucogenicum (Hartmans et al., 2006; Magee and Ward, 2012; Lory, 2014). Of the three clades of slow-growing mycobacteria, the clade designated as “Tuberculosis-Simiae,” encompasses most of the clinically important Mycobacterium species including those related to M. tuberculosis, M. avium, M. gordonae, M, kansasii and M. simiae (Magee and Ward, 2012). The two other clades of the slow-growing species, often referred to as part of the “M. terrae complex,” group together and they form a sister clade to the “Tuberculosis-Simiae” clade. Of the two clades which form the “M. terrae complex,” most of the species closely related to M. terrae are part of a clade that is designated here as the “Terrae” clade (Magee and Ward, 2012; Tortoli, 2012; Ngeow et al., 2015). Adjacent to the “Terrae” clade, the species M. koreense and M. triviale form a distinct clade (designated here as the “Triviale” clade), which is separated from members of the “Terrae” clade by a long branch. It is important to note that in the phylogenetic trees shown in Figure 1, the two clades of fast-growing species do not form a monophyletic grouping, whereas the clades corresponding to the slow-growing mycobacteria group together and form a monophyletic lineage.

We have also compared the relationships observed in the aforementioned phylogenetic trees with the relationships observed in a tree based on 16S rRNA gene sequences, which was extracted from the SILVA Tree of Life Project 128 (Yarza et al., 2008; Quast et al., 2013). This tree is shown in Supplementary Figure 2 with the analogous groups labeled. Overall, in concordance with the core protein-based phylogenetic trees and the tree based on 8 conserved proteins, the slow-growing mycobacterial species corresponding to the “Tuberculosis-Simiae” clade formed a distinct clade in the 16S rRNA tree. The species corresponding to the “Terrae” clade also branched in the immediate proximity of the “Tuberculosis-Simiae” clade, with members of the “Triviale” clade (viz. M. triviale, M. koreense, and M. parakoreense) forming a deeper-branching lineage. However, in contrast to the different trees based on protein sequences, the rapid-growing Mycobacterium species exhibited extensive polyphyly and their interrelationships were poorly resolved. In particular, the members of the “Abscessus-Chelonae” clade formed a monophyletic lineage within the other rapid-growing Mycobacterium species, whereas the relationships among the other rapid growing species were difficult to discern.

Based on genome sequences, the average amino acid identity between different species can be calculated to determine the overall genome relatedness of the species (Konstantinidis and Tiedje, 2007; Richter and Rossello-Mora, 2009; Thompson et al., 2013; Qin et al., 2014; Yarza et al., 2014). Pairwise amino acid identity was calculated based on the conserved protein families between each genome used in the analysis and the results of these analyses are presented in the form of a matrix in Figure 2. An expanded version of this matrix is provided in Supplementary Figure 3. As seen from the AAI matrix (Figure 2), the members of the four main clades observed in the phylogenetic trees (Figure 1 and Supplementary Figure 1) showed higher amino acid identity to members within each clade than to the other Mycobacterium species. Further, members of the “Triviale” clade could be clearly distinguished from the “Terrae” clade, based on their much lower amino acid identity to the members of this latter clade. In addition, members of the “Abscessus-Chelonae” clade exhibited a high degree of amino acid identity (Avg. 92%) to other members of this clade, but significantly lower similarity to members of the “Fortuitum-Vaccae” or the “Tuberculosis-Simiae” clades (Avg. 62%). The results of the genome relatedness analysis support the existence of the four main clades observed in the phylogenetic trees and also the distinctness of the “Triviale” clade from members of the “Terrae” clade.

The results of phylogenomic studies and genomic similarity analysis indicated that the known mycobacterial species can be divided into five main groups including the “Triviale” clade. However, as the branching of species in phylogenetic trees can be affected by a large number of variables (Stackebrandt, 1992; Ludwig and Klenk, 2005; Klenk and Goker, 2010; Gupta, 2016b), it is important to confirm the genetic cohesiveness of the observed clades by independent means not involving phylogenetic analysis. Rare genetic changes, such as insertions and deletions in genes/proteins as well as novel genes/proteins (viz. CSIs and CSPs) which are uniquely shared by an evolutionary related group of organisms constitute synapomorphic characteristics, whose shared presence in a given group of organisms generally results from the occurrence of the genetic changes in a common ancestor of the group (Gupta, 1998, 2014, 2016b; Rokas and Holland, 2000; Dutilh et al., 2008). In our earlier work on Actinobacteria, we described large numbers of CSIs and CSPs which were distinctive characteristics of either the entire phylum or a number of different clades within this phylum at multiple phylogenetic/taxonomic levels (Gao and Gupta, 2005, 2012; Gao et al., 2006; Gupta et al., 2013b). Although the focus of this earlier study was not on mycobacteria, a limited number of CSIs and CSPs which were then specific for the genus Mycobacterium were also identified (Gao et al., 2006; Gao and Gupta, 2012). Since these earlier studies, genome sequences for a large number of other mycobacterial species have become available (Supplementary Table 1). In the present work, we have carried out comprehensive comparative genomic studies on members of the genus Mycobacterium, to identify molecular markers (CSIs and CSPs) that are specific characteristics of either all mycobacterial species or of the identified main clades within this genus. The results of these analyses have identified 172 molecular markers (CSIs and CSPs) that are uniquely found in either all mycobacteria or by the members of different main clades identified by phylogenomic studies. Brief descriptions of the characteristics of the identified molecular markers and their group specificities are provided below.

Our analysis has identified 10 CSIs in proteins involved in diverse functions that are uniquely found in all available mycobacterial homologs. An example of a CSI that is specific for the genus Mycobacterium is shown in Figure 3. In the partial sequence alignment of the protein EgtB (ergothioneine biosynthesis protein), a two amino acid insertion in a conserved region is exclusively found in all members of the genus Mycobacterium, but it is not present in the top 500 homologs of this protein sequence in other bacteria. Ergothionine is a naturally occurring amino acid (thiourea derivative of histidine), whose synthesis is uniquely carried out by only certain groups of actinobacteria as well as some cyanobacteria and fungi (Fahey, 2001). More detailed sequenced information for this CSI as well as sequence information for 9 other CSIs in important proteins, which are also specific for the genus Mycobacterium is provided in Supplementary Figures 4–13 and their main characteristics are summarized in Table 1. Of the described CSIs, the CSI in the protein orotidine 5'-phosphate decarboxylase (Supplementary Figure 7) was identified in our earlier work (Gao and Gupta, 2012). Although the number of sequenced mycobacterial genomes has increased many folds, this CSI is still found only in members of the genus Mycobacterium.

We have previously described a number of CSPs, whose homologs were uniquely found in the then sequenced mycobacterial species (Gao et al., 2006; Gao and Gupta, 2012). In light of the large increase in the number of sequenced mycobacterial genomes, the group specificities of the previously described CSPs were re-examined. Results of these analyses reveal that despite >20-fold increase in the number of sequenced mycobacterial genomes since these CSPs were first identified (Gao et al., 2006), 9 of the CSPs reported in our earlier work are still specific for members of the genus Mycobacterium and no homologs showing significant similarities to these proteins are present in other bacteria (Table 2). In view of the unique shared presence of these 10 CSIs and 9 CSPs by either all or most members of the genus Mycobacterium (except for an isolated exception), the genetic changes leading to these genetic markers most likely initially occurred in a common ancestor of the genus Mycobacterium and then retained by all descendant species.

The “Abscessus-Chelonae” clade, also referred to as M. chelonae or M. abscessus complex (Adékambi and Drancourt, 2004; Medjahed et al., 2010; Tortoli, 2012; Wee et al., 2017), consists of six members and it has recently gained clinical attention in light of its emerging pathogenicity to humans (Medjahed et al., 2010; Tortoli, 2014). In the phylogenetic trees constructed in our work, members of this clade form a monophyletic grouping which comprises the deepest branching lineage among the Mycobacterium species (Figures 1A,B and Supplementary Figure 1). The deep branching of the “Abscessus-Chelonae” clade in comparison to the other Mycobacterium species is also independently supported by 4 CSIs in four different proteins which are commonly shared by the homologs of all other mycobacterial species except those from the “Abscessus-Chelonae” clade. One example of a CSI depicting this pattern is presented in Figure 4, where in the partial sequence alignment of Nif3-like dinuclear metal center hexameric protein, a two amino acid deletion in a conserved region is present in all members of the genus Mycobacterium except members of the “Abscessus-Chelonae” clade. Additional information for this CSI and the sequence information for the three other CSIs exhibiting similar species distributions is provided in Supplementary Figures 14–17 and their main characteristics are summarized in Table 1. Based upon the species distributions of these CSIs, the genetic changes leading to them have likely occurred in a common ancestor of the other Mycobacterium species after the divergence of the “Abscessus-Chelonae” clade.

Our analyses have also identified 27 CSIs in proteins involved in diverse functions that are uniquely shared by members of the “Abscessus-Chelonae” clade providing strong evidence of the genetic cohesiveness and distinctness of this group of mycobacteria. Two examples of the CSIs specific for the “Abscessus-Chelonae” clade are shown in Figure 5. Figure 5A shows a partial sequence alignment of the protein uracil phosphoribosyltransferase, where a six amino acid insertion in a conserved region is present in all members of the “Abscessus-Chelonae” clade but absent in the homologs from all other Mycobacterium species as well as other groups of bacteria. Likewise, Figure 5B shows a four amino acid deletion in the sequence alignment of protein L-histidine N(alpha)-methyltransferase, which is also specific for the “Abscessus-Chelonae” clade. More detailed information for these CSIs and the 25 other identified CSIs, which are also specific for the “Abscessus-Chelonae” clade, is provided in Supplementary Figures 15, 18–43 and their main characteristics are summarized in Table 3. In addition to these CSIs, our work has also identified 24 CSPs listed in Table 2, for which homologs exhibiting significant similarity are only found in members of the “Abscessus-Chelonae” clade. Thus, the distinctness of the “Abscessus-Chelonae” clade from all other mycobacteria is strongly supported by 51 highly-specific molecular signatures identified in this work.

The “Fortuitum-Vaccae” clade as designated here (see Figure 1) encompasses all rapid-growing mycobacterial species, except those from the “Abscessus-Chelonae” clade. In the present work, 4 CSIs and 10 CSPs have been identified that are specific for either all or most members of the “Fortuitum-Vaccae” clade and support the monophyletic clustering of these species as observed in the phylogenomic trees (Figure 1). One of the identified CSIs, which are specific for the “Fortuitum-Vaccae” clade, is found in the LacI family transcriptional regulator. In the partial sequence alignment of this protein shown in Figure 6, a five amino acid insert in a conserved region is exclusively found in different members of the “Fortuitum-Vaccae” clade but it is not found in any other mycobacteria. Three other CSIs showing similar species specificities are present in three other proteins. Detailed sequence information for all of these CSIs is provided in the Supplementary Figures 44–47 and the main characteristics of all CSIs specific for the “Fortuitum-Vaccae” clade are summarized in Table 4.

BLASTp searches on the protein sequences from the genome of Mycobacterium aurum (LSHTM) have also identified 10 CSPs, whose homologs, except for rare exceptions, are only found in the “Fortuitum-Vaccae” clade of Mycobacterium species. Most of these CSPs are hypothetical proteins and their characteristics are summarized in Table 5. For the first four CSPs listed in Table 5, the homologs are present in different members of the “Fortuitum-Vaccae” clade, while for the remaining six CSPs, although they are specific for the “Fortuitum-Vaccae” clade, homologs were not detected in some members of this clade. In all, our identification of 14 molecular markers (4 CSIs and 10 CSPs), which are uniquely shared by members of the “Fortuitum-Vaccae” clade support its monophyletic origin and genetic cohesiveness.

The slow-growing Mycobacterium species generally form a monophyletic clade in most phylogenetic trees based on protein sequences (see Figure 1 and Supplementary Figure 1) as well as those based on the 16S rRNA gene sequences (see Supplementary Figure 2) (Devulder et al., 2005; Kim et al., 2005; Hartmans et al., 2006; Mignard and Flandrois, 2008; Magee and Ward, 2012; Tortoli, 2012; Quast et al., 2013; Lory, 2014; Wang et al., 2015; Wee et al., 2017). The monophyly of the slow-growing Mycobacterium clade is also supported by 3 CSIs and 4 CSPs that have been identified in this study. One example of a CSI that is largely specific for the slow-growing Mycobacterium clade is shown in Figure 7. In the sequence alignment of alkyl-aryl sulfatase protein, a one amino acid insert in a conserved region is present in all of the homologs from slow-growing Mycobacterium species, but it is not found in the homologs of other Mycobacterium species. Detailed sequence information for this CSI and the two other CSIs showing similar specificities is provided in Supplementary Figures 48–50 and their main characteristics are summarized in Table 4. As noted above, the homologs for four of the identified CSPs (Accession numbers: YP_177721.1, YP_178025.1, WP_011725130.1, WP_003874405.1) are also specifically found in slow-growing Mycobacterium species (Table 5). The last two of these CSPs were identified by our earlier work based on limited number of genomes (Gao and Gupta, 2012) and they continue to be specific for this large clade of mycobacteria. Further, of the identified CSPs, which are specific for the slow-growing mycobacterial clade, three of the CSPs correspond to the PE or PPE family of proteins, which are often involved in mycobacterial virulence (Mukhopadhyay and Balaji, 2011).

In our phylogenetic trees, the slow-growing mycobacterial species form three main clades including a clade consisting of M. triviale and M. koreense (“Triviale” clade). The genetic cohesiveness of these clades of slow-growing mycobacteria is also supported by a large number of molecular signatures that are described below.

The “Tuberculosis-Simiae” clade in our work is comprised of all other slow-growing mycobacteria except those from the “Terrae” and “Triviale” clades. This clade encompasses various pathogenic Mycobacterium species including those from the M. tuberculosis complex, M. avium complex, M. gordonae clade, M. kansasii clade, M. simiae clade, as well as several other slow-growing species (Magee and Ward, 2012; Lory, 2014). We have identified a total of 3 CSIs that are specific for the “Tuberculosis-Simiae” clade (Table 4, Supplementary Figures 51–53). One example of a CSI specific for this clade, which is found in a protein of unknown function is shown in Figure 8, where a single amino acid deletion is found in all members of the “Tuberculosis-Simiae” clade, but it is not present in any other mycobacterial homolog. In addition to these CSIs, BLASTp searches on the proteins found in the genome of Mycobacterium tuberculosis H37Rv have identified 3 CSPs, whose homologs are only found in either all or most members of the “Tuberculosis-Simiae” clade. A summary of the CSPs which are specific for the “Tuberculosis-Simiae” clade is provided in Table 5 and of these CSPs, one protein (Genbank Accession Number NP_218369.1) is annotated as a histone-like protein.

The members of the “M. terrae complex” (Tortoli, 2012; Ngeow et al., 2015) has drawn attention recently as some members of this clade are opportunistic pathogens (Mignard and Flandrois, 2008; Kim et al., 2012, 2013; Tortoli, 2012, 2014; Tortoli et al., 2013; Ngeow et al., 2015; Vasireddy et al., 2016). In the core-genome protein trees and the tree based on 8 conserved proteins, members of the “M. terrae complex” form a monophyletic lineage consisting of two distinct subclades: a larger “Terrae” clade encompassing most of the species from the “M. terrae complex” and a deeper branching “Triviale” clade consisting of M. triviale and M. koreense (M. parakoreense also branches with these species in the 16S rRNA tree, Supplementary Figure 2). The phylogenetic distinctness of this larger “Terrae” + “Triviale” clade is also supported by a number of identified molecular signatures. In this work, we have identified 6 CSIs, which are specific for the larger “Terrae complex” consisting of the “Terrae” + “Triviale” clades (Table 6). Sequence information for one of the CSIs specific for the larger “Terrae complex” is presented in Figure 9A. In this case a four amino acid insertion in the protein ATP-dependent helicase is specifically present in all members of the “Terrae complex,” but it is not present in any other bacteria. Detailed sequence information for this CSI as well as other CSIs specific for this clade is presented in Supplementary Figures 54–59 and summarized in Table 6. In addition to these CSIs, which are commonly shared by the “Terrae” + “Triviale” clades, our analyses have also identified 26 other CSIs listed in Table 6, which are specifically shared by only the members of the “Terrae” clade and not present in M. triviale and M. koreense. An example of such a CSI consisting of a four amino acid insertion found in the protein UDP-N-acetylmuramate–L-alanine ligase is shown in Figure 9B. Sequence information for all the “Terrae” clade CSIs is presented in Supplemntary Figures 35, 60–84 and summarized in Table 6. These CSIs serve to indicate the distinctness of the species from the “Terrae” clade from the deeper branching M. triviale and M. koreense species, which are part of the “Triviale” clade.

Our BLASTp searches on the protein sequences from the genome of M. sinense JDM601 (Zhang et al., 2011) and M. triviale DSM 44153 (Fedrizzi et al., 2017) have also identified many CSPs whose homologs are found specifically in either members of the larger “Terrae complex” or uniquely by species which are part of either the “Terrae” clade or the “Triviale” clade. A summary of these CSPs is provided in Table 7. Of the identified CSPs, two CSPs (viz. accession numbers WP_013830140.1 and WP_013827845.1) are uniquely found in most members of the “Terrae” + “Triviale” clades. However, a large number of the other identified CSPs are specific for only either members of the “Terrae” clade (15 CSPs) or members of the “Triviale” clade (22 CSPs) and their homologs are not detected in other mycobacteria. Four of the CSPs specific for the “Triviale” clade included in Table 7 were also previously identified by Ngeow et al. (2015). The identification of a large number of CSPs, which are uniquely found in either all/most members of the “Terrae” clade or those from the “Triviale” clade again serve to clearly differentiate these two groups of mycobacteria and demarcate them in molecular terms.

Mycobacteriaceae (My.co.bac.te.ri.a.ce´ae. N.L. neut. n. Mycobacterium type genus of the family; suff. -aceae ending to denote a family; N.L. fem. pl. n. Mycobacteriaceae the Mycobacterium family).

The family Mycobacteriaceae contains the type genus Mycobacterium as well as the genera Mycolicibacter gen. nov., Mycolicibacterium gen. nov., Mycolicibacillus gen nov., and Mycobacteroides gen. nov. Additionally, the genus Amycolicoccus is also indicated to be a part of this family (Wang et al., 2010; Parte, 2014). However, the sole type species of this genus, Amycolicoccus subflavus, is now reclassified as Hoyosella subflava (Hamada et al., 2016). The general characteristics of the family Mycobacteriaceae are as described by Magee and Ward (2012) for the genus Mycobacterium. The members of this family are aerobic to microaerophilic, slightly curved or straight rods (0.2–0.6 × 1.0–10 μm), which are acid–alcohol-fast at some stage of growth. Difficult to stain by Gram's-method, but are usually considered Gram-stain-positive. Some species may exhibit filamentous or mycelium-like growth. Cells are nonmotile and asporogenous. Colonies may be white- to cream-colored; some strains produce yellow- or orange-pigmented colonies with or without light stimulation. Whole-organism hydrolysates are rich in meso-diaminopimelic acid, arabinose, and galactose. The peptidoglycan is of the A1g type. Muramic acid moieties are N-glycolated. Cells and cell walls are rich in lipids. These include waxes which have characteristic, chloroform-soluble, mycolic acids with long (60–90 carbon atoms) branched chains. The fatty acid esters released on pyrolysis MS of mycolic acid esters have 22–26 carbon atoms. Cells contain diphosphatidylglycerol, phosphatidyl-ethanolamine, phosphatidylinositol, and phospatidylinositol mannosides as predominant polar lipids, straight-chain saturated, unsaturated, and 10-methyloctadecanoic (tuberculostearic) fatty acids as major fatty acid components, and dihydrogenated menaquinones with nine isoprene units as the predominant isoprenolog. The family includes obligate parasites, saprophytes, and opportunistic forms. The G+C content of genome-sequenced species varies from 57 to 71 (mol %) and genome size ranges from 3.1 to 10.5 Mbp. The members of the family Mycobacteriaceae form a distinct clade in the 16S rRNA tree and they are distinguished from all other members of the order Corynebacteriales by their unique shared presence of conserved signature indels described in this work (Table 1) in the following 10 proteins (viz. serine hydrolase, precorrin-4 C(11)-methyltransferase, NAD(P)H-quinone dehydrogenase, orotidine 5′-phosphate decarboxylase, deoxyribonuclease IV, peptidase C69, SGNH/GDSL hydrolase family protein, succinate dehydrogenase, N-dimethylarginine dimethylaminohydrolase, ergothioneine biosynthesis protein EgtB). Additionally, the homologs of the following nine proteins (accession numbers are in parenthesis) are also uniquely found in members of the family Mycobacteriaceae viz. hypothetical protein (WP_011723520.1), hypothetical protein (WP_011723901.1), MAV_11221(WP_011723955.1), membrane protein (WP_011724283.1), PE-PPE domain-containing protein (WP_011724324.1), DUF2561 domain-containing protein (WP_011724709.1), Membrane protein (WP_009976570.1), hypothetical protein (WP_003876314.1) and hypothetical protein (WP_003874755.1) (see Table 2 in this work).

Mycobacterium (My.co.bac.te´ri.um. Gr. n. mykes a fungus; N.L. neut. n. bacterium, a small rod; N.L. neut. n. Mycobacterium, a fungus rodlet).

The type species is Mycobacterium tuberculosis (Zopf 1883) Lehmann and Neumann 1896 (Approved Lists 1980) (Skerman et al., 1980).

Members of this genus whose are slow-growing bacteria requiring at least 7 days of incubation at optimal temperatures to form colonies. Several species are obligate parasites of human and animals and the genus harbors a number of important human (e.g., Mycobacterium tuberculosis, M. leprae, M. ulcerans) and animal (e.g., Mycobacterium bovis) pathogens. Other phenotypic and chemotaxonomic characteristics of this genus are similar to that for the family Mycobacteriaceae.

Some species from this clade contain a 9–12 nucleotide long insert in helix 18 of the 16S rRNA gene sequence (Supplementary Figure 85; Hartmans et al., 2006; Tortoli, 2014). Species are indicated to generally lack the LivFGMH operon and the shaACDEFG cluster of genes, which encodes respectively for proteins allowing the transportation of leucine, isoleucine and valine into the bacteria and a Na⁺/H⁺ antiporter that is important for the homeostasis of Na⁺ and H⁺ (Wee et al., 2017). Presence of the components of Type VII secretion system has been reported in members of this genus (Wee et al., 2017). The members of this genus form a monophyletic clade in phylogenetic trees constructed based on 16S rRNA gene sequences as well as multiple large datasets of protein sequences described in this work including a tree based on 1941 core mycobacterial proteins, a tree based on 136 core proteins for the phylum Actinobacteria, and a tree based on concatenated sequences for eight conserved housekeeping proteins (viz. RpoA, RpoB, RpoC, GyrA, GyrB, Hsp65, EF-Tu, and RecA). Members of the genus Mycobacterium can be clearly distinguished from other genera within the Mycobacteriaceae family based on conserved signature indels described in this study (Table 4) in the following three proteins, a hypothetical protein, aldehyde dehydrogenase family protein and 23S rRNA (guanosine(2251)-2′-O)-methyltransferase, that are uniquely shared by the members of this genus. In addition, the homologs of the following three proteins (accession numbers are in parenthesis): a histone-like protein HNS (NP_218369.1), a hypothetical protein Rv4010 (YP_004837050.1) and a membrane protein (NP_217322.1), are also unique characteristics of the members of this genus.

The G-C content and genome sizes of the member species ranges from 57.8–69.3 (mol %) to 3.2–7.3 Mbp, respectively.

Mycolicibacter (My.co.li.ci.bac´ter. N.L. n. acidum mycolicum, mycolic acid; N.L. masc. n. bacter, rod; N.L. masc. n. Mycolicibacter, a genus of mycolic acid containing rod-shaped bacteria).

The type species is Mycolicibacter terrae.

The members of the genus Mycolicibacter are commonly referred to as the M. terrae complex. This genus contains species that are slow-growing (more than 7 days) and nonchromogenic with some species that show intermediate growth duration (5–15 days) (Tortoli, 2014; Ngeow et al., 2015). In phylogenetic trees, the Mycolicibacter clade forms a sister clade to a clade comprising of the genus Mycobacterium, harboring other slow-growing mycobacteria. Most members of this genus are non-pathogenic, but some species have been isolated from animal hosts (Tasler and Hartley, 1981) and human patients (Smith et al., 2000). Multiple antibiotic resistance has been reported for many of the isolates (Milne et al., 2009; Zhang et al., 2013b).

The members of this genus form a monophyletic clade in phylogenetic trees based on 16S rRNA gene sequences as well as multiple datasets of gene/protein sequences including a tree based on 1941 core mycobacteria proteins and a tree based on 136 core proteins for the phylum Actinobacteria. The members of the genus Mycolicibacter exhibit a closer relationship to members of the genus Mycolicibacillus in phylogenetic trees, which is also supported by a number of CSIs listed (Table 6) in the proteins ATP-dependent helicase, PDZ domain-containing protein, Ferredoxin reductase, DUF2236 domain-containing protein and two hypothetical protein with the accession number WP_083040170 and DUF4185 domain-containing protein, as well as 2 CSPs (viz. accession numbers WP_013830140.1 and WP_013827845.1) that are commonly shared by the members from these two genera. All of the species from this genus contain a 14 nucleotide insertion in the helix 18 of the 16S rRNA gene (Supplementary Figure 85; Tortoli, 2014). Additionally, the members of this genus are distinguished from members of all other genera within the family Mycobacteriaceae due to their possession of 26 conserved signature indels described in this study (Table 6) present in the following proteins, non-ribosomal peptide synthetase, nucleoside hydrolase, three different indels in TetR family transcriptional regulator, carbon starvation protein A, error-prone DNA polymerase, amidohydrolase, carboxymunconolacton decarboxylase family protein, polyketide cyclase, spirocyclase AveC family protein, TobH protein, UDP-N-acetylmuramate–L-alanine ligase, DUF2236 domain-containing protein, cobaltochelatase subunit CobN, alpha/beta hydrolase, potassium transporter Kef, bifunctional tRNA (adenosine(37)-N6)-threonylcarbamoyltransferase complex dimerization subunit Type 1 TsaB/ribosomal protein alanine acetyltransferase RimI, a membrane protein, DUF222 domain-containing protein, MFS transporter, adenylate/guanylate cyclase domain-containing protein, DUF2029 domain-containing protein and the following hypothetical proteins with the accession numbers (WP_083037591, WP_083040170, WP_083036336 and WP_052618664), that are uniquely found in the members of this genus. In addition, the homologs of the 17 conserved signature proteins, whose accession numbers are as follows (viz. WP_013830140.1, WP_013827845.1, WP_013828100.1, WP_013830932.1, WP_013828443.1, WP_013828919.1, WP_013829267.1, WP_041317168.1, WP_013827978.1, WP_041318963.1, WP_013830185.1, WP_013828762.1, WP_013827315.1, WP_041318191.1, WP_013829648.1, WP_013829864.1, and WP_041317804.1) are also distinctive characteristics of either all or most members of this genus (Table 7).

The members of the genus Mycolicibacter are characterized by high G-C content (66.3–70.3 mol %) and they have relatively short genomes (range 3.87–5.11 Mbp).

The description of Mycolicibacter terrae comb. nov. as well as the descriptions of new name combinations for other species which are part of the genus Mycolicibacter are provided in Table 8.

In addition to the new name combinations for species which are part of the genus Mycolicibacter, we also provide below description of two new species that should also be placed in the genus Mycolicibacter.

Description of Mycolicibacter icosiumassiliensis sp. nov. (i.co.si.u.mas.si.li.en´sis; L. masc. n. icosiumassiliensis, from the combination of Icosium, the Latin name of Algiers where the strain was first isolated and Massilia, the Latin name of Marseille, where the strain was described).

The description of this taxon is as given by Djouadi et al. (2016) for “Mycobacterium icosiumassilensis”. The type strain is 8WA6 (= CSUR P1561 = DSM 100711).

Description of Mycolicibacter sinensis sp. nov. (sin.en´sis. N.L. masc. adj. sinensis means “belonging to China,” indicating the source of the type strain).

The description of this taxon is as given by Zhang et al. (2013b) for “Mycobacterium sinense”. The type strain is JDM601.

Mycolicibacillus (My.co.li.ci.ba.cil´lus. N.L. n. acidum mycolicum, mycolic acid; L. masc. n. bacillus, a small staff or rod; N.L. masc. n. Mycolicibacillus, a genus of mycolic acid containing rod-shaped bacteria).

The type species is Mycolicibacillus trivialis.

The genus Mycolicibacillus is comprised of slow-growing nonchromogenic bacteria requiring more than 7 days of incubation at optimal temperatures to form colonies. In phylogenetic trees, members of this genus form a deep-branching distinct clade that is most closely related to members of the genus Mycolicibacter. A close relationship of the species from the genera Mycolicibacillus and Mycolicibacter is also supported by a number of CSIs listed in Table 6 in the proteins ATP-dependent helicase, PDZ domain-containing protein, ferredoxin reductase, DUF2236 domain-containing protein, non-ribosomal peptide synthetase, hypothetical protein with accession number WP_083040170 and DUF4185 domain-containing protein and CSPs listed in Table 7 (viz. accession numbers WP_013830140.1 and WP_013827845.1) that are commonly shared by these two groups of bacteria. Unlike members of the genus Mycolicibacter, which contain a 14 nucleotide insertion in the helix 18 of the 16S rRNA gene, members of the genus Mycolicibacillus lack an insertion in this position (Tortoli, 2014) (Supplementary Figure 85). In addition, the homologs showing significant sequence similarity for the 22 proteins listed in Table 6 with the accession numbers WP_069390591.1, WP_069390644.1, WP_069390667.1, WP_069390717.1, WP_069391089.1, WP_069391367.1, WP_069391463.1, WP_069391521.1, WP_069391698.1, WP_069391782.1, WP_069391793.1, WP_069392105.1, WP_069392126.1, WP_069392251.1, WP_069392420.1, WP_069392510.1, WP_069392884.1, WP_069392982.1, WP_069392983.1, WP_069393100.1, WP_069393493.1, and WP_069393844.1, are uniquely present in members of this genus. This genus presently contains only three species (M. trivialis, M. koreensis and M. parakoreensis) and their genome sizes (3.89–4.08 Mbp) are among the smallest within the family Mycobacteriaceae. The G+C content of the two sequenced species is 69.4 mol %. Although some members of this genus have been isolated from human patients with pulmonary dysfunction, it is unclear whether they exhibit pathogenicity.

The description of Mycolicibacillus trivialis comb. nov. as well as the descriptions of new name combinations for other species which are part of the genus Mycolicibacillus are provided in Table 9.

Mycobacteroides (My.co.bac.te.ro´i.des. N.L. neut. n. Mycobacterium, a bacterial genus; L. neut. suff. -oides, resembling; N.L. neut. n. Mycobacteroides, a genus resembling Mycobacterium).

The type species is Mycobacteroides abscessus. The genus Mycobacteriodes is comprised of bacteria that are commonly referred to as members of the Abscessus-Chelonae clade. This is another genus within the family Mycobacteriaceae of rapidly-growing bacterial species (besides Mycolicibacterium) which take <7 days to form colonies. Phenotypic characteristics of this genus include a positive 3-day arylsulfatase test, better growth at 30°C than at a 35°C, negative nitrate reductase, negative iron uptake and resistance to polymyxin B (Brown-Elliott and Wallace, 2002). The genome size for the species within this clade ranges from 4.5 to 5.6 Mbp and their G+C content ranges from 63.9 to 64.8 mol %. Phylogenetic studies show that members of the genus Mycobacteriodes form a deep branching monophyletic clade within the family Mycobacteriaceae that is distinct from all other genera within this family. Some members from this genus are known to be involved in causing lung, skin and soft tissue infections (Magee and Ward, 2012; Tortoli, 2014) and some exhibit resistance to multiple antimicrobial drugs (Nessar et al., 2012).

The members of the genus Mycobacteriodes can be reliably distinguished from all other Mycobacteriaceae species as well as other bacteria based upon unique shared presence of 27 CSIs in different proteins listed in Table 3 (viz. uracil phosphoribosyltransferase, L-histidine N(alpha)-methyltransferase, DUF58 domain-containing protein, NADH-quinone oxidoreducatase subunit G, ATP-dependent helicase, tRNA (cytidine(34)-2′-O)-methyltransferase, glutamine-fructose-6-phosphate transaminase (isomerizing), error-prone DNA polymerase, 2-amino-4-hydroxy-6-hydroxymethyldihydropteridine diphosphokinase, DEAD/DEAH box helicase, anion transporter, a membrane protein, nicotinate-nucleotide adenylyltransferase, CoA ester lyase, bifunctional ADP-dependent (S)-NAD(P)H-hydrate dehydratase/NAD(P)H-hydrate epimerase, pyridoxal phosphate-dependent aminotransferase, carotenoid oxygenase, SAM-dependent methyltransferase, phosphoribosylamine-glycine ligase, and hypothetical proteins) and the presence of 24 conserved signature proteins listed in Table 2, (viz. MAB_0188c, MAB_0375, MAB_0601, MAB_2852c, MAB_3058, MAB_3079c, MAB_1107c, MAB_1519, MAB_1642, MAB_0008, MAB_0245c, MAB_2487, MAB_3020c, MAB_1440c, MAB_0014, MAB_0015, MAB_0345, MAB_0448c, MAB_0456, MAB_0460, MAB_2549, MAB_1765, MAB_1767, and MAB_1806) that are also specifically found in these bacteria.

The description of Mycobacteroides abscessus comb. nov. as well as the descriptions of new name combinations for other species which are part of the genus Mycobacteroides are provided in Table 10.

Mycolicibacterium (My.co.li.ci.bac.te´ri.um. N.L. n. acidum mycolicum, mycolic acid; N.L. neut. n. bacterium, a small rod; N.L. neut. n. Mycolicibacterium, a genus of mycolic acid containing rod-shaped bacteria).

The type species Mycolicibacterium fortuitum.

The genus is comprised of rapidly-growing bacterial species, which take <7 days to form colonies upon primary isolation (Parte, 2014). Some other phenotypic characteristics generally common to the members of this genus include absence of pigmentation, positive 3-day arylsulfatase activity (Brown-Elliott and Wallace, 2002), positive for nitrate reductase and iron uptake (Magee and Ward, 2012). Most species are saprophytic and considered non-pathogenic to humans, however some cases of infections and diseases by members of this group have been reported (Stahl and Urbance, 1990; Brown-Elliott and Wallace, 2002; Ripoll et al., 2009). The members of this genus form a monophyletic clade in phylogenetic trees based on concatenated sequences of multiple large datasets of conserved proteins including a tree based on 1941 core proteins from mycobacterial genomes, a tree based on 136 core proteins for the phylum Actinobacteria, and another tree based on concatenated sequences for 8 conserved proteins described in the present study.

The members of the genus Mycolicibacterium can be distinguished from other genera within the family Mycobacteriaceae as well as other bacteria based upon conserved signature indels in the following four proteins viz. LacI family transcriptional regulator, Cyclase, CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferase and CDP-diacylglycerol–serine O-phosphatidyltransferase (Table 4) that are uniquely shared by the members of this genus. Additionally, the homologs of the 10 conserved signature proteins, whose accession numbers are as follows (WP_048630777.1, WP_048632025.1, WP_048632497.1, WP_048634851.1, WP_048633467.1, WP_048633322.1 WP_048631132.1, WP_048634509.1, WP_048630657.1, and WP_048632441.1) are also uniquely found in the members of this genus (Table 5). The genome size for the members of this genus ranges from 3.95 to 8.0 Mbp and their G+C content ranges from 65.4 to 70.3 mol %.

The description of Mycolicibacterium fortuitum comb. nov. as well as the descriptions of new name combinations for other species which are part of the genus Mycolicibacterium are provided in Table 11.

In addition to the new name combinations for species which are part of this genus, we also provide below description of two new species that should also be placed in the genus Mycolicibacterium.

Description of Mycolicibacterium acapulense sp. nov. (a.ce.pul.cen´se. N.L. neut. adj. acapulcense from Acapulco, a town on the Pacific coast of México).

The description of this taxon is as given by Bojalil et al. (1962) for “Mycobacterium acapulensis”. The type strain is AC-103 (= ATCC 14473 = JCM 6402).

Description of Mycolicibacterium komanii sp. nov. (ko.ma´ni.i. N.L. gen. n. komanii named after a town in South Africa where one of the isolates originated from, Komani is the Xhosa name for Queenstown (South Africa)).

The description of this taxon is as given by Gcebe (2015) and Gcebe et al. (2016) for “Mycobacterium komanii”. The type strain is GPK 1020.

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.