On the origin of Mycobacterium ulcerans, the causative agent of Buruli ulcer
© Doig et al.; licensee BioMed Central Ltd. 2012
Received: 25 March 2012
Accepted: 30 May 2012
Published: 19 June 2012
Mycobacterium ulcerans is an unusual bacterial pathogen with elusive origins. While closely related to the aquatic dwelling M. marinum, M. ulcerans has evolved the ability to produce the immunosuppressive polyketide toxin mycolactone and cause the neglected tropical disease Buruli ulcer. Other mycolactone-producing mycobacteria (MPM) have been identified in fish and frogs and given distinct species designations (M. pseudoshottsii, M. shinshuense, M. liflandii and M. marinum), however the evolution of M. ulcerans and its relationship to other MPM has not been defined. Here we report the comparative analysis of whole genome sequences from 30 MPM and five M. marinum.
A high-resolution phylogeny based on genome-wide single nucleotide polymorphisms (SNPs) showed that M. ulcerans and all other MPM represent a single clonal group that evolved from a common M. marinum progenitor. The emergence of the MPM was driven by the acquisition of the pMUM plasmid encoding genes for the biosynthesis of mycolactones. This change was accompanied by the loss of at least 185 genes, with a significant overrepresentation of genes associated with cell wall functions. Cell wall associated genes also showed evidence of substantial adaptive selection, suggesting cell wall remodeling has been critical for the survival of MPM. Fine-grain analysis of the MPM complex revealed at least three distinct lineages, one of which comprised a highly clonal group, responsible for Buruli ulcer in Africa and Australia. This indicates relatively recent transfer of M. ulcerans between these continents, which represent the vast majority of the global Buruli ulcer burden. Our data provide SNPs and gene sequences that can differentiate M. ulcerans lineages, suitable for use in the diagnosis and surveillance of Buruli ulcer.
M. ulcerans and all mycolactone-producing mycobacteria are specialized variants of a common Mycobacterium marinum progenitor that have adapted to live in restricted environments. Examination of genes lost or retained and now under selective pressure suggests these environments might be aerobic, and extracellular, where slow growth, production of an immune suppressor, cell wall remodeling, loss or modification of cell wall antigens, and biofilm-forming ability provide a survival advantage. These insights will guide our efforts to find the elusive reservoir(s) of M. ulcerans and to understand transmission of Buruli ulcer.
The genomes of Mycobacterium ulcerans and Mycobacterium marinum are closely related, sharing >97% overall nucleotide identity , but cause very different kinds of infections in humans. M. marinum causes minor skin infections, characterised by intracellular bacteria and the granulomatous lesions that are features of infection with many mycobacterial pathogens, notably Mycobacterium tuberculosis. In contrast, M. ulcerans causes Buruli ulcer (BU), a slowly progressing, ulcerative disease characterized by necrosis of subcutaneous tissue. BU has a characteristic histopathology with large numbers of extracellular bacteria during the acute phase of the infection, with a marked lack of inflammatory response in advanced infection. This unusual pathology is principally mediated by an immunosuppressive polyketide called mycolactone, which is not produced by M. marinum or M. tuberculosis. In BU patients, mycolactone is present in cutaneous lesions but also diffuses and can be detected in serum [3, 4].
Two main features differentiate the genomes of M. ulcerans and M. marinum. The first is the pMUM megaplasmid, found in M. ulcerans but absent from M. marinum[5, 6]. This plasmid harbours three large genes (mlsA1: 51 kb, mlsA2: 7.2 kb, mlsB: 43 kb) that encode the polyketide synthases (PKSs) required for mycolactone synthesis . The second is the insertion sequence (IS) IS2404 that is absent from M. marinum but present in high copy number (>200) in M. ulcerans genomes. IS2404 expansion in the M. ulcerans genome has led to the inactivation of many genes through disruption of coding and promoter sequences and has mediated the deletion of approximately 1 Mbp of DNA from M. ulcerans compared with M. marinum. Together with loss of DNA, there is also evidence of extensive loss of gene function in M. ulcerans - the genome of M. ulcerans isolate Agy99 harbours 771 pseudogenes (inactivated genes), while the M. marinum genome harbours only 65. Acquisition of foreign DNA, IS expansion, pseudogene accumulation and genome reduction are features in common with bacterial populations that have passed through an evolutionary bottleneck [7–12], suggesting there has been constriction of population size during adaptation to a new, niche environment. Analysis of the M. ulcerans Agy99 genome showed deletion or inactivation of genes expressing potent T-cell antigens, and genes required for pigment biosynthesis, anaerobiosis, and intracellular growth . This profile suggested a bacterium that has adapted to a dark, extracellular environment where slow growth, loss of immunogens and production of an immunosuppressive molecule provide a selective advantage [1, 13]. In contrast, its progenitor, M. marinum, has the characteristics of both a specialist bacterium that can persist within an intracellular environment as well as a generalist that can survive in extracellular conditions. A niche environment for M. ulcerans has not yet been demonstrated although the recent discovery that Australian native possums inhabiting BU endemic areas appear to harbour the bacteria in their gastrointestinal tracts raises some interesting possibilities .
The species definition of M. ulcerans has recently been challenged by the discovery of variously named mycobacteria that also make mycolactones but are not always associated with BU. These mycobacteria, isolated from humans, fish and frogs in diverse geographic locations (including Japan, the Mediterranean, Israel, Belgium and the United States), have been variously named Mycobacterium shinshuense, Mycobacterium marinum Mycobacterium pseudoshottsii, and “Mycobacterium liflandii”, the latter an unofficial species name with no standing in nomenclature [15–19]. Studies of these isolates using an 8-gene multi locus sequence typing (MLST) scheme and patterns of genome deletion suggest that they evolved from a M. marinum common ancestor and subsequently diverged into two principal lineages [20, 21]. These mycolactone-producing Mycobacteria (MPM) lineages have been termed “classical” and “ancestral”. The “classical” MPM lineage includes M. ulcerans isolates associated with BU from Australia, South East Asia and Africa while the “ancestral” lineage includes the fish and frog isolates as well as M. ulcerans BU isolates from Japan, China, Mexico, Surinam and French Guiana [22, 23]. We recently proposed that based on their genetic coherence all MPM should be renamed M. ulcerans, a proposition we will revisit in this paper .
Comparative genomic analysis of closely related bacteria has facilitated dramatic improvements in our understanding of bacterial pathogen evolution . Low-coverage 454 and Illumina sequencing of three M. ulcerans isolates from Ghana and a single isolate from Japan demonstrated the capacity of high throughput sequencing to differentiate M. ulcerans co-circulating in a geographic region . Here we report the sequencing and analysis of whole genomes from a diversity of isolates, which are representative of what we define as the M. ulcerans M. marinum complex (abbreviated hereafter as MuMC). That is, mycobacteria previously identified as M. marinum M. ulcerans and other MPM that share >97% nucleotide identity based on MLST [23, 24]. This collection includes mycobacteria isolated from humans, possums, a bilby (small Australian mammal), fish, frogs, an insect and an armadillo from diverse geographic locations. We used whole genome comparisons to reveal the phylogenomic relationship among these mycobacteria, enabling us to investigate the hypothesis that all MPM (including M. ulcerans) are derived from a common ancestor. Our genome comparisons strongly supported this hypothesis and showed that the niche-adapted genomic signature of M. ulcerans caused by reductive evolution was established in the progenitor of all MPM before their intercontinental dispersal. We also identify DNA segments that could be used to develop molecular diagnostic tools to distinguish MPM causing BU from other members of the MuMC.
Isolate selection and sequence analysis strategy
Isolates used in this study
M. ulcerans isolates
Houedja, Ouinhi, Zou, Benin
Lalo, Couffo, Benin
Adjohoun, Ouémé, Benin
Cayenne, French Guiana
Ga District, Ghana
Ga District, Ghana
Amansie West District, Ghana
Ga District, Ghana
Ga District, Ghana
Ga District, Ghana
*Mu_DL045 (M. marinum)
Wokon, Ouinhi, Zou, Benin
Djigbé, Zé, Atlantique, Benin
Hwegoudo, Zé, Atlantique, Benin
Sagon, Ouinhi, Zou, Benin
*Mu_CC240299 (M. marinum)
*Mu_06-3844 (M. marinum)
*Mu_8765 (M. shinshuense)
*Mu_JKD8071 (M. liflandii)
*Mu_L15 (M. pseudoshottsii)
Point Lonsdale, Australia
Point Lonsdale, Australia
Point Lonsdale, Australia
St Leonards, Australia
East Gippsland, Australia
Point Lonsdale, Australia
Point Lonsdale, Australia
Point Lonsdale, Australia
M. marinum isolates
Bilby (Macrotis lagotis)
Mycolactone-producing mycobacteria form a monophyletic group
The most sensitive method to identify a core set of conserved nucleotide positions among a set of genomes is via read mapping to a reference genome. The major advantages of this approach, compared to de novo assembly and ortholog clustering, are its independence from gene annotation and orthology and high (single nucleotide) resolution. Hence we mapped sequencing reads from each isolate to the M. marinum M reference and defined the MuMC core genome as the set of M. marinum M nucleotide positions that were covered by at least three reads from every isolate. Repeat regions and regions spanning PE/PPE genes were excluded from this analysis due to ambiguous read mapping and poor coverage, respectively. We found a core of 4,362,138 bp (65.7%) and 3,318 CDS of M. marinum M were present in all 35 mycobacteria of the complex.
Geographical restriction of M. ulcerans over short time scales
We used the same approach to construct a phylogenetic tree for the mycolactone plasmid pMUM001. No pMUM001 sequences were detected among the four M. marinum isolates or M. ulcerans isolate Mu_8765 (Lineage 2), indicating that the plasmid was absent from these isolates. Mu_8765 has previously been reported to have lost its plasmid during laboratory passage . All other MPM isolates carried the plasmid, confirming the central role of pMUM001 in the evolution of this complex. A total of 315 pMUM001 SNPs were identified among the MPM isolates, permitting construction of a distance-based phylogenetic tree for the plasmid (Figure 4, inset). The plasmid phylogenetic tree closely matches both the topology and relative branch lengths of the core chromosome tree, consistent with co-evolution of the plasmid and chromosome in each lineage . By comparing read depths for plasmid and chromosome (see Methods), we estimated an average pMUM copy number of 1.8 copies per cell (range 1.3 for Mu_CC240299 to 3.1 for Mu_JKD8071).
Inferred genomic characteristics of M. ulcerans MRCA
First, we identified ancestral SNPs that differentiate the most recent common ancestor (MRCA) of M. ulcerans from the rest of M. marinum, which may have played a role in evolution of the distinctive phenotypic characteristics of M. ulcerans. Using M. marinum “M” as a reference, we identified 4,170 such SNPs (including small indels), spanning the entire chromosome and comprising 607 intergenic, 2,254 synonymous, 1,301 non-synonymous and eight SNPs occurring within pseudogenes (Figure 1, Additional file 3: Table S2). This list of mutations, which includes the recently reported hspR_2 regulatory mutation , will be important for research exploring functional characteristics of M. ulcerans. The identification of these SNPs is also the foundation for evaluation of evolutionary forces acting on specific genes (see later section).
We also explored differences in gene content between M. ulcerans isolated from different geographic locations or from different hosts (Figure 6B, 6C, Additional file 5: Table S4). The core genome of M. ulcerans isolates known to infect humans comprised a strict subset of the core genome of fish and frog isolates. Nine CDS were found only in the fish/frog isolates and may be assumed to be superfluous for M. ulcerans infections in humans, although care needs to be exercised as there were only five fish and frog isolates on which to base this analysis. The products of the nine CDS include five hypothetical proteins, three putative hydrogenases and a putative oxononanoate synthase (BioF2_5). A total of 76 CDS were conserved in M. ulcerans lineages 1 and 2 but absent from lineage 3 (which contains only African and Australian BU isolates) (Figure 6C). These data suggest that the African and Australian BU isolates, which represent the majority of the global burden of BU infection in humans, have undergone further reductive evolution. This may be associated with additional adaptation to a more specialised niche, or genetic drift associated with passing through an evolutionary bottleneck. The 76 genes lost (Additional file 5: Table S4) include some involved in metabolic and respiratory processes, which may be no longer required in the more restricted environment occupied by lineage 3 isolates . We also investigated novel M. ulcerans genes not present in the reference genome (Agy99, an African BU isolate from Lineage 3) via annotation and comparison of de novo assemblies of the novel isolates. The African and Australian isolates, also from Lineage 3, had no novel CDS apart from a putative integrase in six of the ten Australian isolates. The remaining M. ulcerans isolates from other lineages had many more novel CDS as follows; Mu_JKD8071 (128), Mu_06-3844 (93), Mu_CC240299 (93), Mu_L15 (13), Mu_1G897 (8), Mu_8765 (3). The numbers of novel CDS largely reflects the isolates’ similarity to the reference genomes as shown in Additional file 2: Figure S2.
Figure 6 also shows the largest segments of DNA that could be used to discriminate between isolates by species, host and region (genomic locations are listed in Additional file 5: Table S4). For example, to discriminate between a M. ulcerans and a M. marinum isolate, there is a 2,688 bp segment that is present in all M. ulcerans isolates and absent from all other M. marinum isolates. Together with previously described deletions, these DNA segments could provide the basis of a set of DNA diagnostic tests to identify the species, host or region of an unknown isolate.
Noteworthy genes or loci absent in the M. ulcerans MRCA
selenocysteine synthase, required for the synthesis of proteins containing selenocysteine.
Formate dehydrogenase alpha and beta subunits, selenocysteine-containing. Likely role in anaerobic growth.
Anaerobic dehydrogenase (putative nitrate reductase), possible role in anaerobic growth.
Lgt acylates prolipoproteins, thus key role in lipoprotein synthesis. Lgt mutant of S. aureus is less immunogenic.
M. marinum has six plc genes. Four of these are absent in all M. ulcerans. Phospholipase C enzymes can cause direct or indirect enzymatic hydrolysis of host cell membrane phospholipids and appear important for mycobacterial intracellular survival.
Carbon starvation protein, CstA. In E. coli, cstA encodes a peptide transporter and is induced by carbon starvation. Maybe part of a redundant stress response system in M. ulcerans.
Multicopper oxidases protect against oxidative stress. This enzyme has functions in tolerance to copper, and iron and manganese oxidation in a range of bacteria. It catalyses the oxidation of cuprous copper, ferrous iron and diphenolic compounds. In Salmonella, a cueO deletion mutant is less virulent in mouse model. Divergent transcriptional arrangement with cstA.
Catalyzes the trans-addition of three molecules of IPP onto DMAPP to form geranylgeranyl pyrophosphate which is a precursor of the ether-linked lipids. Impact here of diverting all isoprenoid biosynthesis to the non-mevalonate pathway is not known. Possibly more favourable energetically if spending cellular resources on mycolactone synthesis.
Impact of IS2404 and IS2606 on M. ulcerans genome architecture
M. ulcerans genes under selective pressure
Identification of CDS that are under selective pressure can give additional insights into the specific nature of the environment(s) an organism encounters. The ratio dN/dS is the number of non-synonymous substitutions per non-synonymous site (dN) to the number of synonymous substitutions per synonymous site (dS); hence dN/dS > 1 implies diversifying selection, while dN/dS < 1 implies purifying selection through loss of mutations that cause changes at the protein level. To compare selection in the M. marinum isolates against the M. ulcerans isolates, two likelihood function (LF) models were generated that compared a constant dN/dS across all isolates against a model where dN/dS differed across the tree. The dN/dS value across a core set of 2,379 CDS for the full tree was 0.181 (constant model) while the values for the differing sub-branches (varying model) were 0.093 for M. marinum isolates, 0.484 for African Lineage 3 M. ulcerans isolates, 0.515 for Australian Lineage 3 M. ulcerans isolates and 0.459 for M. ulcerans lineages 1–2. This is consistent with relaxation of purifying selection in M. ulcerans compared to M. marinum (log likelihood ratio statistic = 8988.6, degrees of freedom = 63, p < 1e-6), however the younger age of M. ulcerans could also explain this result . This dN/dS variation was consistent across all annotated functional groups: (M. ulcerans isolates dN/dS, M. marinum isolates dN/dS) conserved hypotheticals (0.538, 0.191), cell wall and cell processes (0.524, 0.165), lipid metabolism (0.484, 0.160), regulatory proteins (0.470, 0.160), intermediate metabolism (0.469, 0.148), information pathways (0.402, 0.128). Next we compared dN/dS at the level of individual CDS. Our model parameter start values were taken from the LF and the tree previously used to compare regions using the full set of CDS. This improved the ability of the new LF to be more rapidly and reliably optimised. Of the 2,379 CDS examined, 172 CDS had a dN/dS > 1 for the M. ulcerans isolates (Additional file 7: Table S6). This set was dominated by two functional groups, ‘cell wall and cell processes’ (44 CDS) and ‘conserved hypotheticals’ (56 CDS). Notable CDS in this set of genes under diversifying selection included esxE and esxF that encode WXG100 domains and are orthologs of the secreted T-cell antigens EsxA and EsxB. The stress response sigma factor sigM was also identified as under diversifying selection, perhaps indicating remodelling of the SigM regulon in M. ulcerans (Figure 1). In M. tuberculosis, sigM appears to positively regulate esxE and esxF expression while negatively regulating certain non-ribosomal peptide synthetases and polyketide synthase genes involved in biosynthesis of cell wall lipids [52, 53]. Lipoproteins, another class of cell wall-associated molecule, also featured among the genes displaying evidence of diversifying selection (lppN, lppD, lprA) and perhaps links to the predicted loss of Lgt activity and the ability to acylate prolipoproteins (Table 2).
We also examined CDS that contained only non-synonymous SNPs, prohibiting calculation of a dN/dS ratio, or were not core to M. ulcerans, and compiled a list of additional CDS with an overabundance of non-synonymous SNPs, sorted by the difference between non-synonymous and synonymous SNPs (see Methods). The 50 most variable proteins mirrored the dominant functional groupings of the high dN/dS analysis (Additional file 7: Table S6), including esxE/F as well as enoyl-CoA hydratase echA4_2 and ssp. Ssp (signal peptide protease IV) is a serine protease required for digestion of cleaved signal peptides, indicating that M. ulcerans may be under pressure to modify their secretome. Also of interest is mymT, encoding a copper-binding protein, a cytoplasmic metallothionein and in M. tuberculosis, sequesters excess copper in the bacterial cell .
The close relationship between M. tuberculosis and the MuMC led us to investigate whether the orthologous human T cell epitopes for M. tuberculosis were conserved in M. ulcerans. The genes encoding antigens in human pathogens tend to be under diversifying selection in order to evade host immunity and it has been suggested that this is not the case with M. tuberculosis. We found 20 MuMC orthologs (>80% AA identity) of the 86 human T-cell antigens in M. tuberculosis H37Rv (Additional file 8: Table S7). These orthologs all had very low dN/dS ratios (M. marinum isolates median 0.05, range 0.15, M. ulcerans isolates median 0.17, range 0.95) suggesting that these CDS are under purifying selection in the MuMC as they are in M tuberculosis.
Whole genome sequencing and comparison of this large collection of isolates has shown that all mycolactone-producing mycobacteria (collectively referred to in this study as M. ulcerans) evolved from a common M. marinum progenitor by a combination of horizontal gene transfer (pMUM and phage), IS-mediated deletion and point mutation (Figure 4). These evolutionary processes resulted in extensive gene loss (185 CDS), changes in gene expression (e.g. hsp18 and possibly sigM) and likely changes in gene function through positive selection, establishing M. ulcerans as a highly specialized, niche-adapted mycobacterium. Our data indicate these changes were followed by global dispersal of M. ulcerans and further diversification and adaptive evolution.
Scrutiny of the types of genes gained, lost or modified in conjunction with other experimental evidence provides some clues regarding the nature of the niche environment in which the M. ulcerans MRCA was able to flourish and in which today’s isolates still survive. Foremost among the DNA gained was in the pMUM plasmid. The acquisition of this plasmid occurred early in the evolution of M. ulcerans as demonstrated by the congruent tree topologies inferred from chromosome and plasmid sequence alignments (Figure 4). The synthesis of mycolactone, a potent immunosuppressive small molecule encoded in the pMUM plasmid, presumably gave a population of M. marinum cells the ability to persist in a place their generalist relatives could not. In previous research we have described how recombination and gene conversion has shaped the unusually repetitive gene structure of the mycolactone PKS, such that the 110 kb, three-gene PKS locus on pMUM comprises only 10 kb of unique sequence. This unusual gene structure and its resultant instability are strongly suggestive of intense selection acting on M. ulcerans populations to maintain mycolactone production [5, 56, 57].
The loss of the mevalonate pathway for synthesis of isoprenoid lipids was originally observed in the genome of M. ulcerans Agy99 . In the current study we show that this trait was established in the M. ulcerans MRCA and was probably a key adaptive response of the bacterium following the acquisition of pMUM (Table 2). It is possible that loss of this metabolic capacity freed essential resources for critical mycolactone synthesis or alternatively that these pathways and metabolites were redundant in the niche environment occupied by M. ulcerans. Other potentially significant gene losses in M. ulcerans include selenocysteine synthase (selA), required for the synthesis of proteins containing selenocysteine, and the linked genes encoding the alpha and beta subunits of a putative selenocysteine-containing formate dehydrogenase, with a possible role in anaerobic growth; mutations that suggest the M. ulcerans MRCA lost the ability to grow anaerobically.
It is also striking that genes associated with the intracellular lifestyle of mycobacterial pathogens such as M. marinum and M. tuberculosis have been lost in M. ulcerans. Genes predicted to be inactive include four phospholipase enzymes (PlcB_2,3,5&6) and cueO (Table 2). As well as playing roles in intracellular replication in other bacteria, these CDS are also components of the bacterial cell wall. There appears to have been significant selective pressure on M. ulcerans to reduce or change its cell wall and cell surface antigenic profile. In this respect another noteworthy pseudogene conserved in M. ulcerans is lgt (MUL_1594). Lgt is a prolipoprotein diacylglyceryl transferase that acylates prolipoproteins at a conserved N-terminal cysteine . An lgt mutation in Staphylococcus aureus causes growth rate attenuation, an accumulation of prolipoproteins in the culture supernatant, and reduced activation of innate immune responses [46, 58]. The loss of lgt in M. ulcerans might therefore lead to aberrantly or non-acylated lipoproteins with reduced immunogenicity, like the S. aureus mutant.
The modification of the cell wall appears to have continued in the BU-associated M. ulcerans lineage 3 genomes, which contain an additional 589 pseudogenes or deleted regions of which 30% are predicted to have encoded antigens or cell wall associated proteins, including EsxA_2, EsxA_3, and Hspx_1. The deep branching lineage and clonal nature of the African and Australian lineage 3 isolates, which are most commonly involved in human infections, have the signature of passing through a second evolutionary bottleneck: gene deletions, further loss of gene function, chromosomal rearrangements and the expansion of another IS (IS2606 from the pMUM plasmid). Each of the M. ulcerans lineages probably represents different ecotypes, reflecting adaptation to related but distinct niche environments. It may be that each lineage is best described as an M. ulcerans ecovar.
We tried to identify a temporal signal in our sequence data to estimate divergence dates for particular lineages of the MuMC. However, while our phylogenetic inferences were highly robust (100% bootstrap values for major branches of our tree, Figure 4) no linear correlation between branch length and year of isolation was observed. We suspect there is variation in the effective number of generations per year across the complex, perhaps related to different niches, reservoirs or modes of living. Such variation has recently been observed in M. tuberculosis. The lack of a temporal signal in these data raises doubts around previous estimates of divergence time that have assumed a constant molecular clock rate among the lineages of the complex [26, 34].
However, sequence comparisons did reveal a compelling correlation between genes undergoing positive selection - as revealed by dN/dS analysis - and those CDS inactivated or deleted (Additional file 6: Table S5, Additional file 7: Table S6), with mutations in all these groups skewed towards CDS involved in cell wall and lipid biosynthesis. These patterns point to significant selective pressures acting on M. ulcerans populations to devote resources (substrate and energy) towards the synthesis of mycolactones and modification of cell wall structures. Intriguingly, many of the cell wall metabolites lost via mutation are known to be highly antigenic in other mycobacteria. One interpretation of these observations is that the bacteria are responding to pressures from a host immune system, a point argued in a previous study . Furthermore, when one considers that mycolactone is a potent immune suppressor with reported specificity for a mammalian microRNA that controls T-cell chemotaxis , this in turn leads to the idea that the niche occupied by M. ulcerans is a higher organism with a complex immune system. The discovery that Australian possums inhabiting BU endemic areas are susceptible to BU disease and harbour large number of M. ulcerans in their gastrointestinal tracts is consistent with this idea . Nevertheless, these arguments are not consistent with the significant lack of variation seen among M. ulcerans proteins with putative T-cell epitopes (Additional file 8: Table S7), where an immune escape hypothesis would predict hypervariability not hyperconservation in these regions.
This study has also reinforced the close relationship between isolate origin and genotype for M. ulcerans strains that cause Buruli ulcer, notably those from Africa and Australia, where multiple isolates from one region were sequenced. The complete resolution of strain differences afforded by whole genome sequencing has shown how the genotype of M. ulcerans strains from two African countries correlate with place of origin. Isolates from the east of Benin are distinct from isolates in the West of the country or from a different country (Ghana, Figure 5). These findings suggest - as also demonstrated in a previous study in Ghana - that M. ulcerans transmission and microevolution generally occurs at a local level and therefore the source of the bacterium is somewhat fixed within a local region . This observation should guide our thinking regarding the source of the bacteria in BU endemic areas, indicating that animal reservoirs of M. ulcerans are unlikely to be highly mobile. Yet, one should also consider that the relative paucity of genomic differences between isolates from Ghana and Benin also reflects the relatively recent spread of the bacteria across this entire region. Efforts to establish the rate of mutation of these isolates or genome analysis of a more temporally and spatially diverse collection of isolates from this region might help estimate the amount of time M. ulcerans has been extant in West Africa.
This study has given us our most comprehensive insight yet into the M. ulcerans-M. marinum complex. The isolates examined in the study covered a wide geographical diversity, and while future studies should include M. ulcerans isolates from other endemic countries together with more M. marinum isolates, whole genome comparisons of 35 isolates across the genetic diversity of the complex have yielded nucleotide-level granular detail of each isolate and their relatedness to each other. We have shown that all MPM are a single lineage whose divergence from M. marinum was characterised by acquisition of the pMUM plasmid, which conferred the ability to synthesise mycolactones (largely responsible for the immunomodulatory properties of M. ulcerans) and also introduced IS2404 and IS2606 into the chromosome; this was then followed by IS2404 and IS2606 expansion and extensive loss and modification of gene function, consistent with an evolutionary bottleneck and adaptation to a new niche. Based on these shared features we propose that all members of this lineage should be considered M. ulcerans. Furthermore, we suggest that M. ulcerans sublineages (such as the three we have identified in this study) should be considered M. ulcerans ecovars.
Examination of the classes of genes lost and modified in the M. ulcerans MRCA suggests a bacterial population occupying a niche environment that is aerobic, osmotically stable, dark (at least for lineage 3 isolates as shown by loss of UV-protecting pigment genes in a previous study) and possibly extracellular given the number of genes known to be involved in intracellular survival lost from M. ulcerans[1, 13, 21]. Support for an extracellular niche is also found in two separate studies showing M. ulcerans elaborates a mycolactone-rich extracellular matrix and specifically expresses a surface protein that promotes adherence during initial stages of biofilm formation [42, 62]. Furthermore, it is difficult to ignore the striking depletion of immunogens and modification of cell-wall proteins in M. ulcerans that may be subject to interactions with the immune defences of a host organism.
This study has clarified our understanding of the origins of an emerging human pathogen, provided important insights into the search for the reservoir(s) of this pathogen and generated a significant resource for future research. It is clear that M. ulcerans is not a generalist, saprophytic mycobacterium but a highly specialized occupant of a protected niche. This concept should be considered in ongoing research efforts to pinpoint the reservoir(s) of M. ulcerans and to understand the transmission of BU.
Several bioinformatics approaches were used to analyse the sequence data. The workflow we developed is summarised in Additional file 9: Figure S1 and involved analysis of sequence differences and similarities between all strains. A mapping-to-reference approach was applied as well as de novo assembly of the sequenced reads. A summary of the isolate data is presented in Additional file 1: Table S1, including the number of reads and the de novo assembly statistics for each strain.
Bacterial isolates and sequencing
The isolates of M. ulcerans and M. marinum used in this study are described in Table 1. Genomic DNA was prepared from at least 50 mg (wet weight) cell pellets, harvested from mycobacteria that were cultivated on egg-yolk agar media as previously described .
Genome sequencing of the isolates (except for Mu_DL045) was performed using an Illumina GAIIx DNA sequencer with 36, 76, or 100 cycle paired-end chemistry, generating read lengths of 35–101 nucleotides and 16-240x mean read depth per isolate. Sequencing statistics are provided in Additional file 1: Table S1. Mu_DL045 was sequenced using the Ion Torrent platform to generate single-ended reads. Read data for the study isolates (Mu_06-3845, Mu_06-3846, Mu_07-1082, Mu_113, Mu_119, Mu_74, Mu_81, Mu_93, Mu_05142109, Mu_JKD8170, Mu_08009899, Mu_JKD8049) have been deposited in the NCBI Sequence Read Archive (SRA) under accession SRP004497 .
Prior to further analysis, reads were filtered to remove those containing ambiguous base calls, any reads <24 nucleotides in length and reads containing only homopolymers. All reads were trimmed from the 3’ end removing bases with the Illumina Read Segment Quality Control Indicator “B”.
An in-house Python utility called Nesoni, which uses SHRiMP2  for read mapping, was used to map novel sequence reads to reference genomes . Nesoni identified SNPs and indels up to ~10 bp and predicted the consequences of SNPs and indels on CDS, e.g. frameshifts, premature stop codons and non-synonymous codon changes. Five reference genomes and fragments were used for read mapping, (sequence name, Genbank accession number) (M. ulcerans Agy99, CP000325), (M. marinum strain “M”, CP000854), (pMUM001 plasmid, BX649209), (pMUM002 plasmid, EU271968) and (pMUM003 plasmid, EU271967). To facilitate comparison of M. marinum and M. ulcerans reference genomes using Nesoni, short read data for M. marinum “M” was simulated via in silico shredding of the reference genome to generate a set of 100 bp short reads (read depth x15 insert length 280 bp) using a utility function available within Nesoni. Average read depth and the number of SNPs detected by mapping to M. ulcerans Agy99 and M. marinum M are given in Additional file 1: Table S1.
The number of reads mapping to unique regions of the plasmids was used to infer plasmid copy number. The ratio of the mode of read depth for the reference plasmid (pMUM002 for Mu_06-3844, Mu_CC240299, Mu_JKD8071, Mu_L15, Mu_1G897 and pMUM001 for all other M. ulcerans strains) was compared to the mode of read depth for the chromosome.
Variant sites (SNPs and indels) identified using Nesoni were concatenated to form a multiple alignment, from which a phylogenetic tree was constructed using SplitsTree4 (neighbour-joining tree based on uncorrected P distances) . All major bifurcations had >95% bootstrap support based on 1,000 runs. The tree shown in Figure 4 is based on SNPs identified using M. marinum M as the reference, however the same tree was recovered when using M. ulcerans Agy99 as the reference. SNPs were classified by their position on the annotated references of M. ulcerans Agy99 and M. marinum “M” and identified as ‘intergenic’, ‘coding’ or ‘pseudogene’; ‘coding’ SNPs were further divided into synonymous, non-synonymous or nonsense SNPs. All SNPs called within annotated ISs or PE/PPE genes were excluded from SNP analysis due to the unreliability of read mapping in repetitive regions [8, 67].
De novo assembly of short-read sequences
Genomes were assembled using Velvet  with the help of VelvetOptimiser  to select appropriate parameter values for k-mer size, expected coverage and coverage cutoff. The Velvet assemblies yielded an average N50 of 18,399 bp; assembly statistics are given in Additional file 1: Table S1. Repetition of hundreds of copies of IS2404 and IS2606 within the M. ulcerans isolates was problematic and restricted the contig lengths, as the size of sequenced DNA fragments (≤240 bp) was much smaller than the length of the ISs (1365–1444 bp). To identify IS disrupted contigs, paired reads were mapped to consensus IS2404 and IS2606 sequences. Where a read was mapped to the IS but its partner read was not, the unmapped partner reads were pooled and de novo assembled to identify the DNA sequences that flank the IS sequences. This analysis demonstrated that 51% (11% s.d.) of the contig breakpoints terminated at either IS2404 or IS2606. Mycobacterial genomes contain hundreds of PE/PPE genes, which pose problems for sequencing, read mapping and assembly due to their repetitive, low G + C content 3’ sequences . Hence PE/PPE genes were excluded from most downstream analysis including SNP calling and gene content analysis. Contigs were annotated by Prokka, an in-house bacterial genome annotation pipeline using Glimmer, EasyGene, Orpheus and GeneMark.
Core and pan genome analysis
The core genome was defined as all positions in the reference genomes that were covered by ≥3 reads in all isolates; this information was derived from depth plots generated by Nesoni. The set of genomic material conserved in one set of isolates, A, but absent from a second disjoint set of isolates, B, was obtained by taking the core genome of A minus (all positions conserved among A) and removing any bases present in B.
To define the pan genome, CDS annotated in the reference or de novo assembly contigs were subjected to ortholog clustering. CDS were translated to amino acid sequences and clustered using CD-Hit  with a 70% amino acid sequence identity threshold. A core set of M. ulcerans CDS was defined by identifying all clusters containing a M. ulcerans reference CDS and one or more CDS fragments, with a combined length within 15% of the reference sequence, from each novel M. ulcerans genome.
To assess the completeness of the set of study isolates, we calculated the median number of new genes added to the pan-genome for every additional genome, using 1,000 permutations of the study isolates (M. ulcerans and M. marinum). A trend was modelled by fitting a power curve to estimate the exponent, indicating whether the pan genome is ‘open’ or ‘closed’ [37, 38].
The relative rates of change at synonymous and non-synonymous sites (dN/dS) within protein coding regions was calculated using the BioPerl module Bio::Align::DNAStatistics and the PyCogent toolkit . The dN/dS analysis was performed for both individual CDS (using BioPerl and PyCogent) and for aggregate concatenated CDS (using PyCogent). To test for selection across multiple isolates using multiple CDS, analysis was restricted to an alignment of DNA sequences from all CDS that were present, without a premature stop codon, in all isolates (2,379 CDS). This alignment, together with an unweighted neighbour-joining tree of the isolates, was provided as input to PyCogent modules that calculated the dN/dS values across subsets of the isolates. The CNF codon substitution model was used for all analyses . A representative set of study isolates was used: M. marinum Mm_1726, Mm_99/84, Mm_99/87, Mm_99/89; M. ulcerans lineage 1 Mu_CC40299, Mu_JKD8071, Mu_DL045, Mu_1G897; M. ulcerans lineage 2 Mu_8765; M. ulcerans lineage 3 Mu_05142109, Mu_Agy99. This alignment, together with the corresponding subset of the tree and a codon substitution model, was used to create a likelihood function (LF) which models selection using the PyCogent framework . By specifying how the dN/dS parameter values of the model vary, optimisation of the LF can estimate dN/dS across the whole tree, a sub branch of isolates or individual isolates. The dN/dS metric cannot identify CDS that contained only non-synonymous SNPs but have no synonymous SNPs. As an alternative method of identifying CDS with an overabundance of non-synonymous SNPs, we calculated the average difference between the number of non-synonymous and synonymous SNPs normalised for the length of CDS. Functional groups used to classify CDS were defined by the reference genome annotations and based on COG v1.0 / CDD groupings.
In order to identify the pseudogenes in the M. ulcerans MRCA we examined the CDS SNP changes that were common in all M. ulcerans isolates but not present in any M. marinum isolates. All PE/PPE, ISs and M. marinum pseudogenes were excluded from this analysis. The remaining M. marinum CDSs were examined to determine whether the SNPs called by read mapping the M. ulcerans isolates onto M. marinum resulted in a change to the CDS that would likely render it inactive within the M. ulcerans isolates; namely (i) indels that result in frameshifts, (ii) premature stop codons truncating encoded protein, (iii) a change to the start codon, (iv) an IS insertion or (v) a deletion of part or whole of the CDS.
Average nucleotide identity
Protein-coding DNA sequence
Insertion or deletion in genome
Multi-locus sequence typing
Mycolactone producing mycobacteria
Most Recent Common Ancestor
Mycobacterium tuberculosis complex
Mycobacterium ulcerans-Mycobacterium marinum complex
The length of the contig such that 50% of the assembly is contained in contigs of length N50 or greater.
This research was supported in part by the National Health and Medical Research Council of Australia (628930 & 65195), the Stop Buruli Initiative of the UBS-Optimus Foundation and the Victorian Government Department of Health. We are grateful to Scott Coutts and Jessica Porter for sequencing support, and Matthew Wakefield and Helen Lindsay for their expert assistance regarding PyCogent.
- Stinear TP, Seemann T, Pidot S, Frigui W, Reysset G, Garnier T, Meurice G, Simon D, Bouchier C, Ma L: Reductive evolution and niche adaptation inferred from the genome ofMycobacterium ulcerans, the causative agent of Buruli ulcer. Genome Res. 2007, 17 (2): 192-200. 10.1101/gr.5942807.PubMed CentralView ArticlePubMedGoogle Scholar
- Stinear TP, Seemann T, Harrison PF, Jenkin GA, Davies JK, Johnson PDR, Abdellah Z, Arrowsmith C, Chillingworth T, Churcher C: Insights from the complete genome sequence ofMycobacterium marinumon the evolution ofMycobacterium tuberculosis. Genome Res. 2008, 18 (5): 729-741. 10.1101/gr.075069.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Sarfo FS, Le Chevalier F, Aka N, Phillips RO, Amoako Y, Boneca IG, Lenormand P, Dosso M, Wansbrough-Jones M, Veyron-Churlet R: Mycolactone diffuses into the peripheral blood of Buruli ulcer patients–implications for diagnosis and disease monitoring. PLoS Negl Trop Dis. 2011, 5 (7): e1237-10.1371/journal.pntd.0001237.PubMed CentralView ArticlePubMedGoogle Scholar
- Hong H, Coutanceau E, Leclerc M, Caleechurn L, Leadlay PF, Demangel C: Mycolactone Diffuses fromMycobacterium ulcerans-Infected Tissues and Targets Mononuclear Cells in Peripheral Blood and Lymphoid Organs. PLoS Negl Trop Dis. 2008, 2 (10): e325-10.1371/journal.pntd.0000325.PubMed CentralView ArticlePubMedGoogle Scholar
- Stinear TP, Mve-Obiang A, Small PL, Frigui W, Pryor MJ, Brosch R, Jenkin GA, Johnson PD, Davies JK, Lee RE: Giant plasmid-encoded polyketide synthases produce the macrolide toxin ofMycobacterium ulcerans. Proc Natl Acad Sci U S A. 2004, 101 (5): 1345-1349. 10.1073/pnas.0305877101.PubMed CentralView ArticlePubMedGoogle Scholar
- Stinear TP, Hong H, Frigui W, Pryor MJ, Brosch R, Garnier T, Leadlay PF, Cole ST: Common evolutionary origin for the unstable virulence plasmid pMUM found in geographically diverse strains ofMycobacterium ulcerans. J Bacteriol. 2005, 187 (5): 1668-1676. 10.1128/JB.187.5.1668-1676.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MT, Churcher CM, Bentley SD, Mungall KL: Comparative analysis of the genome sequences ofBordetella pertussis,Bordetella parapertussisandBordetella bronchiseptica. Nat Genet. 2003, 35 (1): 32-40. 10.1038/ng1227.View ArticlePubMedGoogle Scholar
- Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, Goodhead I, Rance R, Baker S, Maskell DJ, Wain J: High-throughput sequencing provides insights into genome variation and evolution inSalmonellaTyphi. Nat Genet. 2008, 40 (8): 987-993. 10.1038/ng.195.PubMed CentralView ArticlePubMedGoogle Scholar
- Hiller NL, Janto B, Hogg JS, Boissy R, Yu S, Powell E, Keefe R, Ehrlich NE, Shen K, Hayes J: Comparative genomic analyses of seventeenStreptococcus pneumoniaestrains: insights into the pneumococcal supragenome. J Bacteriol. 2007, 189 (22): 8186-8195. 10.1128/JB.00690-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersson JO, Andersson SG: Genome degradation is an ongoing process in Rickettsia. Mol Biol Evol. 1999, 16 (9): 1178-1191. 10.1093/oxfordjournals.molbev.a026208.View ArticlePubMedGoogle Scholar
- Cole ST, Eiglmeier K, Parkhill J, James KD, Thomson NR, Wheeler PR, Honore N, Garnier T, Churcher C, Harris D: Massive gene decay in the leprosy bacillus. Nature. 2001, 409 (6823): 1007-1011. 10.1038/35059006.View ArticlePubMedGoogle Scholar
- Thomson NR, Howard S, Wren BW, Holden MT, Crossman L, Challis GL, Churcher C, Mungall K, Brooks K, Chillingworth T: The complete genome sequence and comparative genome analysis of the high pathogenicityYersinia enterocoliticastrain 8081. PLoS Genet. 2006, 2 (12): e206-10.1371/journal.pgen.0020206.PubMed CentralView ArticlePubMedGoogle Scholar
- Demangel C, Stinear TP, Cole ST: Buruli ulcer: reductive evolution enhances pathogenicity ofMycobacterium ulcerans. Nat Rev Microbiol. 2009, 7 (1): 50-60. 10.1038/nrmicro2077.View ArticlePubMedGoogle Scholar
- Fyfe JAM, Lavender CJ, Handasyde KA, Legione AR, O'Brien CR, Stinear TP, Pidot SJ, Seemann T, Benbow ME, Wallace JR: A Major Role for Mammals in the Ecology ofMycobacterium ulcerans. PLoS Negl Trop Dis. 2010, 4 (8): e791-10.1371/journal.pntd.0000791.PubMed CentralView ArticlePubMedGoogle Scholar
- Kazumi Y, Ohtomo K, Takahashi M, Mitarai S, Sugawara I, Izumi J, Andoh A, Hasegawa H: Mycobacterium shinshuenseisolated from cutaneous ulcer lesion of right lower extremity in a 37-year-old woman. Kekkaku. 2004, 79: 437-441.PubMedGoogle Scholar
- Mve-Obiang A, Lee RE, Portaels F, Small PL: Heterogeneity of mycolactones produced by clinical isolates ofMycobacterium ulcerans: implications for virulence. Infect Immun. 2003, 71 (2): 774-783. 10.1128/IAI.71.2.774-783.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Rhodes MW, Kator H, McNabb A, Deshayes C, Reyrat JM, Brown-Elliott BA, Wallace R, Trott KA, Parker JM, Lifland B: Mycobacterium pseudoshottsiisp. nov., a slowly growing chromogenic species isolated from Chesapeake Bay striped bass (Morone saxatilis). Int J Syst Evol Microbiol. 2005, 55 (Pt 3): 1139-1147.View ArticlePubMedGoogle Scholar
- Stragier P, Ablordey A, Bayonne L, Lugor Y, Sindani I, Suykerbuyk P, Wabinga H, Meyers W, Portaels F: Heterogeneity amongMycobacterium ulceransisolates from Africa. Emerg Infect Dis. 2006, 12: 844-847. 10.3201/eid1205.051191.PubMed CentralView ArticlePubMedGoogle Scholar
- Ucko M, Colorni A, Kvitt H, Diamant A, Zlotkin A, Knibb WR: Strain variation in Mycobacterium marinum fish isolates. Appl Environ Microbiol. 2002, 68 (11): 5281-5287. 10.1128/AEM.68.11.5281-5287.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaser M, Rondini S, Naegeli M, Stinear T, Portaels F, Certa U, Pluschke G: Evolution of two distinct phylogenetic lineages of the emerging human pathogenMycobacterium ulcerans. BMC Evol Biol. 2007, 7: 177-10.1186/1471-2148-7-177.PubMed CentralView ArticlePubMedGoogle Scholar
- Yip MJ, Porter JL, Fyfe JA, Lavender CJ, Portaels F, Rhodes M, Kator H, Colorni A, Jenkin GA, Stinear T: Evolution ofMycobacterium ulceransand other mycolactone-producing mycobacteria from a commonMycobacterium marinumprogenitor. J Bacteriol. 2007, 189 (5): 2021-2029. 10.1128/JB.01442-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Röltgen K, Stinear TP, Portaels F, Certa U, Pluschke G: The genome, evolution and diversity ofMycobacterium ulcerans. Evol. 2012, 12 (3): 522-529.Google Scholar
- Kaser M, Hauser J, Small P, Pluschke G: Large sequence polymorphisms unveil the phylogenetic relationship of environmental and pathogenic mycobacteria related toMycobacterium ulcerans. Appl Environ Microbiol. 2009, 75 (17): 5667-5675. 10.1128/AEM.00446-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Pidot S, Asiedu K, Kaser M, Fyfe J, Stinear T: Mycobacterium ulceransand Other Mycolactone-Producing Mycobacteria Should Be Considered a Single Species. PLoS Negl Trop Dis. 2010, 4 (7): e663-10.1371/journal.pntd.0000663.PubMed CentralView ArticlePubMedGoogle Scholar
- Parkhill J, Wren BW: Bacterial epidemiology and biology - lessons from genome sequencing. Genome Biol. 2011, 12 (10): 230-10.1186/gb-2011-12-10-230.PubMed CentralView ArticlePubMedGoogle Scholar
- Qi W, Kaser M, Roltgen K, Yeboah-Manu D, Pluschke G: Genomic diversity and evolution ofMycobacterium ulceransrevealed by next-generation sequencing. PLoS Pathog. 2009, 5 (9): e1000580-10.1371/journal.ppat.1000580.PubMed CentralView ArticlePubMedGoogle Scholar
- Pearson T, Busch JD, Ravel J, Read TD, Rhoton SD, U'Ren JM, Simonson TS, Kachur SM, Leadem RR, Cardon ML: Phylogenetic discovery bias inBacillus anthracisusing single-nucleotide polymorphisms from whole-genome sequencing. Proc Natl Acad Sci U S A. 2004, 101 (37): 13536-13541. 10.1073/pnas.0403844101.PubMed CentralView ArticlePubMedGoogle Scholar
- Portaels F, Meyers WM, Ablordey A, Castro AG, Chemlal K, de Rijk P, Elsen P, Fissette K, Fraga AG, Lee R: First Cultivation and Characterization ofMycobacterium ulceransfrom the Environment. PLoS Negl Trop Dis. 2008, 2 (3): e178-10.1371/journal.pntd.0000178.PubMed CentralView ArticlePubMedGoogle Scholar
- De Gentile PL, Mahaza C, Rolland F, Carbonnelle B, Verret JL, Chabasse D: Cutaneous ulcer from Mycobacterium ulcerans. Apropos of 1 case in French Guiana. Bull Soc Pathol Exot. 1992, 85 (3): 212-214.PubMedGoogle Scholar
- Roltgen K, Qi W, Ruf M-T, Mensah-Quainoo E, Pidot SJ, Seemann T, Stinear TP, Kaser M, Yeboah-Manu D, Pluschke G: Single Nucleotide Polymorphism Typing ofMycobacterium ulceransReveals Focal Transmission of Buruli Ulcer in a Highly Endemic Region of Ghana. PLoS Negl Trop Dis. 2010, 4 (7): e751-10.1371/journal.pntd.0000751.PubMed CentralView ArticlePubMedGoogle Scholar
- Stragier P, Hermans K, Stinear T, Portaels F: First report of a mycolactone-producingMycobacteriuminfection in fish agriculture in Belgium. FEMS Microbiol Lett. 2008, 286 (1): 93-95. 10.1111/j.1574-6968.2008.01264.x.View ArticlePubMedGoogle Scholar
- Tsukamura M, Kaneda K, Imaeda T, Mikoshiba H: A taxonomic study on a mycobacterium which caused a skin ulcer in a Japanese girl and resembledMycobacterium ulcerans. Kekkaku. 1989, 64 (11): 691-697.PubMedGoogle Scholar
- Trott KA, Stacy BA, Lifland BD, Diggs HE, Harland RM, Khokha MK, Grammer TC, Parker JM: Characterization of aMycobacterium ulcerans-like infection in a colony of African tropical clawed frogs (Xenopus tropicalis). Comp Med. 2004, 54 (3): 309-317.PubMedGoogle Scholar
- Stinear TP, Jenkin GA, Johnson PD, Davies JK: Comparative genetic analysis ofMycobacterium ulceransandMycobacterium marinumreveals evidence of recent divergence. J Bacteriol. 2000, 182 (22): 6322-6330. 10.1128/JB.182.22.6322-6330.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Stragier P, Ablordey A, Meyers W, Portaels F: GenotypingMycobacterium ulceransandMycobacterium marinumby using mycobacterial interspersed repetitive units. J Bacteriol. 2005, 187: 1639-1647. 10.1128/JB.187.5.1639-1647.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645. 10.1101/gr.092759.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Tettelin H, Riley D, Cattuto C, Medini D: Comparative genomics: the bacterial pan-genome. Curr Opin Microbiol. 2008, 11 (5): 472-477. 10.1016/j.mib.2008.09.006.View ArticlePubMedGoogle Scholar
- Mira A, Martin-Cuadrado AB, D'Auria G, Rodriguez-Valera F: The bacterial pan-genome:a new paradigm in microbiology. Int Microbiol. 2010, 13 (2): 45-57.PubMedGoogle Scholar
- Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM: DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007, 57 (Pt 1): 81-91.View ArticlePubMedGoogle Scholar
- Yang F, Yang J, Zhang X, Chen L, Jiang Y, Yan Y, Tang X, Wang J, Xiong Z, Dong J: Genome dynamics and diversity ofShigellaspecies, the etiologic agents of bacillary dysentery. Nucleic Acids Res. 2005, 33 (19): 6445-6458. 10.1093/nar/gki954.PubMed CentralView ArticlePubMedGoogle Scholar
- Huson DH: SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics. 1998, 14 (1): 68-73. 10.1093/bioinformatics/14.1.68.View ArticlePubMedGoogle Scholar
- Pidot SJ, Porter JL, Tobias NJ, Anderson J, Catmull D, Seemann T, Kidd S, Davies JK, Reynolds E, Dashper S: Regulation of the 18 kDa heat shock protein inMycobacterium ulcerans: an alpha-crystallin orthologue that promotes biofilm formation. Mol Microbiol. 2010, 78 (5): 1216-1231. 10.1111/j.1365-2958.2010.07401.x.View ArticlePubMedGoogle Scholar
- Wiedenbeck J, Cohan FM: Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches. FEMS Microbiol Rev. 2011, 35 (5): 957-976. 10.1111/j.1574-6976.2011.00292.x.View ArticlePubMedGoogle Scholar
- Liu Y, Harrison PM, Kunin V, Gerstein M: Comprehensive analysis of pseudogenes in prokaryotes: widespread gene decay and failure of putative horizontally transferred genes. Genome Biol. 2004, 5 (9): R64-10.1186/gb-2004-5-9-r64.PubMed CentralView ArticlePubMedGoogle Scholar
- Rezwan M, Grau T, Tschumi A, Sander P: Lipoprotein synthesis in mycobacteria. Microbiol. 2007, 153 (Pt 3): 652-658.View ArticleGoogle Scholar
- Stoll H, Dengjel J, Nerz C, Gotz F: Staphylococcus aureusdeficient in lipidation of prelipoproteins is attenuated in growth and immune activation. Infect Immun. 2005, 73 (4): 2411-2423. 10.1128/IAI.73.4.2411-2423.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Sander P, Rezwan M, Walker B, Rampini SK, Kroppenstedt RM, Ehlers S, Keller C, Keeble JR, Hagemeier M, Colston MJ: Lipoprotein processing is required for virulence ofMycobacterium tuberculosis. Mol Microbiol. 2004, 52 (6): 1543-1552. 10.1111/j.1365-2958.2004.04041.x.View ArticlePubMedGoogle Scholar
- Bakala N'goma JC, Schue M, Carriere F, Geerlof A, Canaan S: Evidence for the cytotoxic effects ofMycobacterium tuberculosisphospholipase C towards macrophages. Biochim Biophys Act. 2010, 1801 (12): 1305-1313. 10.1016/j.bbalip.2010.08.007.View ArticleGoogle Scholar
- Dubey AK, Baker CS, Suzuki K, Jones AD, Pandit P, Romeo T, Babitzke P: CsrA regulates translation of theEscherichia colicarbon starvation gene,cstA, by blocking ribosome access to thecstAtranscript. J Bacteriol. 2003, 185 (15): 4450-4460. 10.1128/JB.185.15.4450-4460.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Achard ME, Tree JJ, Holden JA, Simpfendorfer KR, Wijburg OL, Strugnell RA, Schembri MA, Sweet MJ, Jennings MP, McEwan AG: The multi-copper-ion oxidase CueO ofSalmonella entericaserovar Typhimurium is required for systemic virulence. Infect Immun. 2010, 78 (5): 2312-2319. 10.1128/IAI.01208-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Rocha EP, Smith JM, Hurst LD, Holden MT, Cooper JE, Smith NH, Feil EJ: Comparisons of dN/dS are time dependent for closely related bacterial genomes. J Theor Biol. 2006, 239 (2): 226-235. 10.1016/j.jtbi.2005.08.037.View ArticlePubMedGoogle Scholar
- Agarwal N, Woolwine SC, Tyagi S, Bishai WR: Characterization of theMycobacterium tuberculosissigma factor SigM by assessment of virulence and identification of SigM-dependent genes. Infect Immun. 2007, 75 (1): 452-461. 10.1128/IAI.01395-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Raman S, Puyang X, Cheng TY, Young DC, Moody DB, Husson RN: Mycobacterium tuberculosisSigM positively regulates Esx secreted protein and nonribosomal peptide synthetase genes and down regulates virulence-associated surface lipid synthesis. J Bacteriol. 2006, 188 (24): 8460-8468. 10.1128/JB.01212-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Gold B, Deng H, Bryk R, Vargas D, Eliezer D, Roberts J, Jiang X, Nathan C: Identification of a copper-binding metallothionein in pathogenic mycobacteria. Nat Chem Biol. 2008, 4 (10): 609-616. 10.1038/nchembio.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Comas I, Chakravartti J, Small PM, Galagan J, Niemann S, Kremer K, Ernst JD, Gagneux S: Human T cell epitopes ofMycobacterium tuberculosisare evolutionarily hyperconserved. Nat Genet. 2010, 42 (6): 498-503. 10.1038/ng.590.PubMed CentralView ArticlePubMedGoogle Scholar
- Hong H, Demangel C, Pidot SJ, Leadlay PF, Stinear T: Mycolactones: immunosuppressive and cytotoxic polyketides produced by aquatic mycobacteria. Nat Prod Rep. 2008, 25 (3): 447-454. 10.1039/b803101k.PubMed CentralView ArticlePubMedGoogle Scholar
- Pidot SJ, Hong H, Seemann T, Porter JL, Yip MJ, Men A, Johnson M, Wilson P, Davies JK, Leadlay PF: Deciphering the genetic basis for polyketide variation among mycobacteria producing mycolactones. BMC Genomics. 2008, 9: 462-10.1186/1471-2164-9-462.PubMed CentralView ArticlePubMedGoogle Scholar
- Hashimoto M, Tawaratsumida K, Kariya H, Kiyohara A, Suda Y, Krikae F, Kirikae T, Gotz F: Not lipoteichoic acid but lipoproteins appear to be the dominant immunobiologically active compounds inStaphylococcus aureus. J Immunol. 2006, 177 (5): 3162-3169.View ArticlePubMedGoogle Scholar
- Ford CB, Lin PL, Chase MR, Shah RR, Iartchouk O, Galagan J, Mohaideen N, Ioerger TR, Sacchettini JC, Lipsitch M: Use of whole genome sequencing to estimate the mutation rate ofMycobacterium tuberculosisduring latent infection. Nat Genet. 2011, 43 (5): 482-486. 10.1038/ng.811.PubMed CentralView ArticlePubMedGoogle Scholar
- Huber CA, Ruf MT, Pluschke G, Kaser M: Independent loss of immunogenic proteins inMycobacterium ulceranssuggests immune evasion. Clin Vaccine Immunol. 2008, 15 (4): 598-606. 10.1128/CVI.00472-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Guenin-Mace L, Carrette F, Asperti-Boursin F, Le Bon A, Caleechurn L, Di Bartolo V, Fontanet A, Bismuth G, Demangel C: Mycolactone impairs T cell homing by suppressing microRNA control of L-selectin expression. Proc Nat Acad Sci. 2011, 108 (31): 12833-12838. 10.1073/pnas.1016496108.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsollier L, Brodin P, Jackson M, Kordulokova J, Tafalmeyer P, Carbonelle E, Aubry J, Milon G, Legras P, Saint Andre JP: mpact ofMycobacterium ulceransbiofilm on transmissibility to ecological niches and Buruli ulcer pathogenesis. PLoS Pathogens. 2007, 3 (5): e62-10.1371/journal.ppat.0030062.PubMed CentralView ArticlePubMedGoogle Scholar
- NCBI Sequence Read Archive. http://www.ncbi.nlm.nih.gov/sra.
- David M, Dzamba M, Lister D, Ilie L, Brudno M: SHRiMP2: Sensitive yet Practical Short Read Mapping. Bioinformatics. 2011, 27 (7): 1011-1012. 10.1093/bioinformatics/btr046.View ArticlePubMedGoogle Scholar
- Harrison P, Seemann T: Nesoni. http://www.vicbioinformatics.com/nesoni.shtml.
- Huson DH, Bryant D: Application of Phylogenetic Networks in Evolutionary Studies. Mol Biol Evol. 2006, 23 (2): 254-267.View ArticlePubMedGoogle Scholar
- Harris SR, Feil EJ, Holden MT, Quail MA, Nickerson EK, Chantratita N, Gardete S, Tavares A, Day N, Lindsay JA: Evolution of MRSA during hospital transmission and intercontinental spread. Science. 2010, 327 (5964): 469-474. 10.1126/science.1182395.PubMed CentralView ArticlePubMedGoogle Scholar
- Zerbino D, Birney E: Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18 (5): 821-10.1101/gr.074492.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Gladman S: VelvetOptimiser. http://bioinformatics.net.au/software.velvetoptimiser.shtml.
- Sampson SL: Mycobacterial PE/PPE Proteins at the Host-Pathogen Interface. Clin Dev Immunol. 2011, 2011: 11-View ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- Knight R, Maxwell P, Birmingham A, Carnes J, Caporaso JG, Easton BC, Eaton M, Hamady M, Lindsay H, Liu Z: PyCogent: a toolkit for making sense from sequence. Genome Biol. 2007, 8 (8): R171-10.1186/gb-2007-8-8-r171.PubMed CentralView ArticlePubMedGoogle Scholar
- Yap VB, Lindsay H, Easteal S, Huttley G: Estimates of the effect of natural selection on protein-coding content. Mol Biol Evol. 2010, 27 (3): 726-734. 10.1093/molbev/msp232.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.