Previous Article | Next Article ![]()
Journal of Bacteriology, March 2005, p. 1668-1676, Vol. 187, No. 5
0021-9193/05/$08.00+0 doi:10.1128/JB.187.5.1668-1676.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Roland Brosch,1
Thierry Garnier,1
Peter F. Leadlay,5 and
Stewart T. Cole1
Unité de Génétique Moléculaire Bactérienne,1 Plate-Forme 4-Intégration et Analyse Génomiques, Génopole, Institut Pasteur, Paris, France,4 Department of Microbiology, Monash University, Clayton, Australia,2 Department of Chemistry,3 Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom5
Received 8 October 2004/ Accepted 29 November 2004
|
|
|---|
|
|
|---|
In the original description of mycolactone, George et al. (12) extracted acetone-soluble lipids from the Malaysian M. ulcerans strain 1615 and identified mycolactone A/B ([M + Na]+ at m/z 765) with Z- and E- isomers of a 12-member macrolactone linked via an ester bond to a fatty acyl side chain. The structure of mycolactone A/B has been confirmed by complete chemical synthesis (10) (Fig. 1).
![]() View larger version (10K): [in a new window] |
FIG. 1. Structures of mycolactone A (Z- 4',5') and B (E- 4',5') ([M + Na]+ at m/z 765).
|
One might expect such high sequence homology to promote frequent genetic rearrangements and therefore a diversity of metabolites, but M. ulcerans strains from around the world have thus far been shown to produce a very restricted repertoire of mycolactones. A study of 34 M. ulcerans isolates collected worldwide showed that they all make the same lactone core with minor variations in the acyl side chain (22). These variations have been largely attributed to various degrees of oxidation at C-12' of the side chain (15, 22), and it has been proposed that this is due to the activity (or lack of activity) of a specific P450 monooxygenase (encoded by the plasmid gene MUP053) (15, 30).
In this study we used a large-insert M. ulcerans DNA clone library to examine the stability of plasmid pMUM001. We then explored the distribution and structure of this plasmid in other M. ulcerans strains using PCR, DNA sequencing, pulsed-field gel electrophoresis (PFGE), and Southern hybridization.
|
|
|---|
(mrr-hsdRMS-mcrBC)
80dlacZ
M15
lacX74 deoR recA1 araD139
(ara, leu)7697 galU galK rpsL endA1 nupG] and XL2-Blue (recA1 endA1 gyrA96 thi-1 hsdR17 supE44 relA1 lac [F' proAB lacIqZ
M15 Tn10(Tetr) Amy Camr]) were cultivated in Luria-Bertani broth at 37°C. Mycobacterium marinum strain M was cultivated at 32°C in 7H9 Middlebrook medium (Becton Dickinson) supplemented with oleic acid-albumin-dextrose-catalase (Difco). The following 10 M. ulcerans clinical isolates were used: Agy99 (origin, Ghana, 1999; this strain was used for the M. ulcerans genome sequencing project), Kob (origin, Ivory Coast, 2001; kindly provided by Laurent Marsollier, Institut Pasteur), 1615 (origin, Malaysia, 1963; kindly provided by Pamela Small, University of Tennessee), Chant (origin, southeastern Australia, 1993), IP105425 (from the reference collection of the Institut Pasteur and derived from reference strain ATCC 19428; origin, southeastern Australia, 1948), 01G897 (origin, French Guiana, 1991), ITM-5114 (origin, Mexico, 1958), ITM-941331 (origin, Papua New Guinea, 1994), ITM-98912 (origin, People's Republic of China, 1997), and ITM-941328 (origin, Malaysia, 1994). M. ulcerans isolates were grown as described above for M. marinum. M. ulcerans isolates whose designations begin with ITM were kindly provided by Françoise Portaels (Institute for Tropical Medicine, Antwerp, Belgium). Liquid chromatography-sequential mass spectrometry analysis of mycolactones. Lipid fractions from M. ulcerans were extracted and analyzed for mycolactones as previously described (11, 15)
Oligonucleotides and DNA methods. The oligonucleotides used in this study are shown in Table 1. Standard methods were used for subcloning, PCR, and automated DNA sequencing. DNA sequences were assembled and annotated by using Gap4 and Artemis, respectively (2, 26).
|
View this table: [in a new window] |
TABLE 1. Oligonucleotides used in this study
|
-32P]dCTP. Probes were denatured by heating at 100°C and were then added to hybridization buffer (5x SSC, 0.1% SDS, 1% skim milk) to a final concentration of approximately 10 ng/ml. Hybridization was performed at 68°C for 18 h. The hybridization solution was then removed, and three stringency washes were performed (one wash for 5 min in 2x SSC-0.1% SDS at room temperature and then two washes for 10 min in 0.1x SSC-0.1% SDS at 68°C). Each membrane was then washed in 2x SSC and sealed in clear plastic film before detection with a Storm PhosphorImager (Molecular Dynamics). Probe stripping was performed by washing the membrane twice for 20 min at 68°C with 0.1% SDS-0.2 M NaOH. The sizes of DNA restriction fragments were estimated with the Sigmagel software (Jandel Scientific) by using a lambda low-range DNA size ladder (New England Biolabs) to calibrate the gel and blot images. BAC library construction. A whole-genome M. ulcerans bacterial artificial chromosome (BAC) library was constructed as described previously for Mycobacterium tuberculosis (4). Briefly, genomic DNA from M. ulcerans Agy99 was prepared in agarose plugs as described above and subjected to partial HindIII digestion. The DNA was separated under PFGE conditions. Partially digested DNA in the size range from 40 to 120 kb was cloned into the unique HindIII site of the vector pBeloBAC11 and then used to transform E. coli DH10B by electroporation. The resulting clones were stored in Luria-Bertani broth containing 15% glycerol in a 96-well format at 80°C.
BAC plasmid DNA preparation. BAC DNA for automated sequencing was extracted by the method of Brosch et al. (4). For subcloning of BACs, DNA was prepared from 40-ml overnight E. coli cultures, and the plasmid DNA was extracted as previously described (4).
Phylogenetic analysis. The sequences from the four plasmid loci (repA, parA, mls, MUP045) that were present in all 10 M. ulcerans strains were concatenated in frame to produce a 1,266-bp semantide for each strain. These sequences were then aligned with CLUSTAL W (31). In the same way the plasmid sequences obtained from the seven M. ulcerans strains that contained the repA, parA, MUP011, mls load, mlsAT(II), MUP038, and MUP045 loci were concatenated in frame to produce a 2,208-bp semantide composed of these seven loci. Phylogenetic analysis was performed with the MEGA software, version 2.1 (19). P distances were used throughout this study as the overall level of sequence divergence was small. Values for synonymous (dS) and nonsynonymous (dN) mutation frequencies were calculated by the method of Nei and Gojobori (23), and standard errors of the means of these values were estimated by the method of Nei and Jin (24). The dS and dN calculations were performed by using the dSdNqw program (7).
|
|
|---|
![]() View larger version (40K): [in a new window] |
FIG. 2. Dotter analysis of the pMUM001 DNA sequence, highlighting regions where there are repetitive DNA sequences. Direct repeat sequences are indicated by lines that are parallel to the main diagonal, while inverted repeats are indicated by lines that are perpendicular to the main diagonal. The sites of homologous recombination surrounding the start of mlsA1 and mlsB that led to the creation of plasmid deletion derivatives are indicated by the shaded circles.
|
![]() View larger version (50K): [in a new window] |
FIG. 3. Mapping of the deletion variants of pMUM001. (A) Scaled, circular maps of pMUM001 and the two types of deletion derivatives, with a proposed model for recombination-mediated deletion. The positions of all HindIII sites are indicated. On the outer circles, the solid arrows indicate the locations of several key genes. The sites of recombination are circled and indicated by the crossed dotted lines. The inner grey circles show the sequences spanned by BAC clones. For the deletion derivatives, the HindIII sites where the vector pBeloBAC11 was cloned are also shown. (B) Expanded view of the regions of recombination within pMUM001 surrounding the loading modules at the start of mlsA1 and mlsB that gave rise to the deletion variants. All HindIII and PstI sites are indicated. The grey block between the dotted lines indicates the zone of 100% nucleotide identity that was subject to recombination. The 200-bp sequence hybridizing to probe 74 is also shown. (C) Gel electrophoresis of the preparations resulting from PstI RE digestion of 21 M. ulcerans Agy99 BAC clones, showing the presence of two subfamilies that span the mlsB and mlsA genes. (D) Southern hybridization analysis of the gel shown in panel C, confirming the presence of two copies of the mls loading module sequences in pMUM001 and the presence of single copies in the deletion variants. The different sizes of the hybridizing bands are due to the sites of cloning into pBeloBAC11, which contains three PstI sites.
|
![]() View larger version (37K): [in a new window] |
FIG. 4. Mapping of pMUM in seven M. ulcerans strains. (A) PFGE and Southern hybridization analyses with five selected PCR-derived probes from pMUM001 and undigested and XbaI-digested genomic DNA extracted from M. ulcerans and M. marinum. Lane 1, M. ulcerans Agy99; lane 2, M. ulcerans Kob; lane 3, M. ulcerans 1615; lane 4, M. ulcerans Chant; lane 5, M. ulcerans IP105425; lane 6, M. ulcerans ITM-5114; lane 7, M. ulcerans ITM-941331; lane 8, M. marinum strain M. (B) Physical maps of pMUM for the seven M. ulcerans strains, deduced from the Southern hybridization experiments shown in panel A, showing the plasmid size, the positions of all XbaI sites, and the toxin status of each strain as determined by liquid chromatography-sequential mass spectrometry. Question marks indicate that the exact region deleted from the mls locus could not be determined.
|
|
View this table: [in a new window] |
TABLE 2. PCR analysis of 10 different M. ulcerans strains for the presence of eight plasmid-associated genes
|
PFGE and Southern hybridization were used to study in more detail the structure of the plasmids in 7 of the 10 M. ulcerans strains. M. ulcerans DNA was separated by PFGE. This DNA was then hybridized with a pool of probes derived from five of the plasmid markers described in Table 2. The results are shown in Fig. 4A and demonstrate that there was a considerable difference in plasmid size among the isolates, with the sizes ranging from 59 to 174 kb. M. ulcerans strains harboring plasmids less than 110 kb long would not be expected to produce mycolactone as the Mls biosynthetic cluster is encoded by genes encompassing approximately 110 kb of DNA. Screening of lipid extracts from the seven isolates by liquid chromatography-mass spectrometry confirmed this prediction and that of the PCR analysis, as neither mycolactone nor its cometabolites were detected in extracts from M. ulcerans Kob (a recent West African isolate with a 101-kb plasmid), M. ulcerans ITM-5114 (a Mexican isolate with a 59-kb plasmid), and M. ulcerans IP105425 (an isolate from the culture collection of the Institut Pasteur, derived from reference strain ATCC 19428, with a 76-kb plasmid). Digestion with XbaI and hybridization with the five pooled plasmid markers resulted in a profile consisting of two, three, or four bands. For each strain, the sum of its XbaI fragments was equal to the size of the linear plasmid form in the absence of XbaI digestion (Fig. 4). This demonstrated that none of the plasmids had new, additional XbaI fragments. Hybridization experiments with individual probes then permitted linking of plasmid markers to particular XbaI fragments and construction of low-resolution maps (Fig. 4B). The three mycolactone-minus strains had large deletions of 75, 98, and 115 kb. The hybridization data showing the absence of MUP038 (encoding the type II thioesterase), together with the PCR data showing an absence of the acyltransferase domain of module 5 in mlsA1 and the acyltransferase domain of modules 1 and 2 in mlsB, confirmed that these deletions had occurred, at least in part, within their respective mls loci.
Only the strains with four XbaI fragments (M. ulcerans Agy99, 1616, Chant, and ITM-941331) produced mycolactone, and thus by definition they must all contain an intact mls locus. This fact was supported by the presence of conserved 54- and 13-kb fragments, corresponding to the locus harboring the mlsA genes and MUP038. Therefore, the size variations detected among these four strains occurred in the regions flanking the mls genes.
Plasmid variation correlates with the presence of different mycolactone cometabolites. For M. ulcerans strains Chant and ITM-941331, some of the plasmid size variation could be attributed to the absence of a region that includes the gene MUP053 (encoding a P450 hydroxylase). The product of MUP053 is predicted to hydroxylate the mycolactone side chain at C-12' to produce mycolactone A/B with a mass of [M + Na]+ at m/z 765 (30). Strains lacking the hydroxyl group at C-12' produce a mycolactone with a mass of [M + Na]+ at m/z 749. This metabolite has been called mycolactone C (22), and it is a characteristic of Australian strains. The absence of MUP053 in the Australian M. ulcerans strain Chant correlates well with the presence of mycolactone C and the absence of mycolactone A/B (Fig. 5). However, M. ulcerans ITM-941331 also lacks MUP053, yet this strain produces the same mycolactone profile as M. ulcerans Agy99 (15) (data not shown).
![]() View larger version (17K): [in a new window] |
FIG. 5. Liquid chromatography-mass spectrometry analysis of the lipid extract from Australian isolate M. ulcerans Chant, showing the absence of a mycolactone ([M+Na]+ at m/z 765.5) and the presence of the nonhydroxylated mycolactone ([M+Na]+ at m/z 749.5). (A) Ion trace for m/z 765.5; (B) ion trace for m/z 749.5.
|
![]() View larger version (57K): [in a new window] |
FIG. 6. Phylogenetic analysis of 10 M. ulcerans strains with selected plasmid markers. (A) Alignment of 1,266-bp sequences derived from the four concatenated pMUM protein-coding loci present in all 10 M. ulcerans strains. Only variable nucleotides are shown. A period indicates identity with strain M. ulcerans ITM-941331. (B) Alignment of 2,208-bp sequences derived from the seven concatenated pMUM protein-encoding loci present in six M. ulcerans strains. (C) Neighbor-joining tree of the phylogenetic relationships among the 10 M. ulcerans strains inferred from comparisons of the 1,266-bp sequences. PNG, Papua New Guinea. (D) Neighbor-joining tree of the phylogenetic relationships among the six M. ulcerans strains inferred from comparisons of the 2,208-bp sequences. (E) Neighbor-joining tree of the phylogenetic relationship among six M. ulcerans and five M. marinum genotypes as revealed by a previous sequence analysis of seven chromosomal protein-encoding loci in 18 M. ulcerans isolates and 22 M. marinum isolates (28). (F) CLUSTAL W alignment of the predicted amino acid sequences of a 348-bp region of MUP053 for the five M. ulcerans strains positive for this gene.
|
MUP053, encoding a putative P450 monooxygenase with a possible role in modifying mycolactone, had an uneven distribution in the strains. However, MUP053 was present in strains from Africa, Malaysia, the People's Republic of China, and Mexico, and these strains spanned the known genetic diversity of the species. The levels of DNA and amino acid identity for MUP053 in these strains were 98 and 96%, respectively, which were equal to the values for other plasmid sequences (Fig. 6F). This suggests that MUP053 was present in a progenitor M. ulcerans and was subsequently lost from some strains as the species evolved.
|
|
|---|
One of the unusual features of pMUM001 is the unprecedented DNA homology among the functional domains of the mls genes. While the mls genes occupy 105 kb of pMUM001, this region contains less than 10 kb of unique sequence (30). This extraordinary economy of sequence is reflected in Fig. 2 and suggests that the mls genes were created de novo by successive recombination events, such as in-frame duplications and deletions from a core set of PKS sequences. The precise origin of such a core gene set remains obscure as DNA database searches have revealed no orthologous genes, but the significant levels of amino acid identity to PKS sequences from other species of mycobacteria and streptomyces point to a likely origin among the actinomycetes. In addition to suggesting a recent evolutionary origin for mycolactone biosynthesis, the extended DNA sequence homology also implies that such an arrangement is inherently unstable and acts as a substrate for general recombination. In this study we showed that in M. ulcerans Agy99, pMUM001 is unstable and that recombination between two homologous sequences gave rise to two deletion variants. The larger 109-kb variant, represented by BAC clone 22D03, contains an intact origin of replication and is thus likely to be maintained within a cell population. Cells harboring the 22D03 variant should be unable to produce mycolactone but could theoretically still produce the acyl side chain. However, the smaller 65-kb deletion variant, represented by BAC clone 22A01, should be lost from the population upon cell division as it is not capable of autonomous replication, despite having the genes required for synthesis of the mycolactone core. Spontaneous mycolactone-minus and avirulent M. ulcerans mutants were first reported by George et al. (12) and were used to demonstrate the key role of mycolactone in virulence. Mycolactone confers a pale yellow color to colonies, and mycolactone-minus mutants are readily observed as white colony variants when they are grown on Löwenstein-Jensen medium. We attempted to isolate white colony variants of M. ulcerans Agy99 to try to identify the 109-kb deleted form of pMUM001. While white colonies were readily detected on Löwenstein-Jensen medium, their growth after subculturing was highly impaired, and we were unable to generate the biomass required for additional studies, such as PFGE. Nevertheless, investigation of other M. ulcerans strains revealed deleted forms of pMUM similar to those identified in M. ulcerans Agy99 (in particular, in M. ulcerans Kob), and these deleted forms had corresponding toxin-minus phenotypes. Each strain tested had a different plasmid size, and the mapping data showed that deletions had occurred to various extents and in different regions of pMUM. Recombination between homologous sequences is one explanation for this variety, but given the large number of insertion sequences in pMUM (30), another possibility is that insertion sequences also mediated some of these plasmid rearrangements.
It is probably significant that no pMUM-minus M. ulcerans strains were found. While such mutants may exist, the recent finding that pMUM contains an active partition (par) locus (T. P. Stinear et al., Microbiology, in press) means that spontaneous curing is likely to be an infrequent event. par loci are cis-acting elements that ensure that daughter cells faithfully receive a copy of an episome during cell division.
Based on the assumption that the clinical isolates used in this study were originally mycolactone proficient and thus contained intact pMUM, it appears that spontaneous toxin-minus mutants, caused by deletion of M. ulcerans plasmid DNA, are a common occurrence. We have not been able to calculate the frequency with which deletion mutants arise, but for some strains it appears to be very high. M. ulcerans Agy99 and Kob were recent clinical isolates from West Africa with minimal laboratory passaging. The DNA used for the M. ulcerans Agy99 BAC library was prepared from a liquid culture that was at its fourth passage since primary isolation, and M. ulcerans Kob was at its third passage. One outcome of this work is to highlight the care that researchers must take to continually test the plasmid and mycolactone status of the M. ulcerans strains used in their work. Periodic passaging of lab strains through an animal virulence model (e.g., mice) may be required to ensure maintenance of a homogeneous population of mycolactone-producing bacteria.
Plasmid instability contrasts most strikingly with the fact that M. ulcerans isolates recovered from diverse geographic locations around the world produce a relatively homogeneous range of mycolactones (22). This apparent paradox leads compellingly to the notion that there is strong purifying selection for maintenance of a mycolactone-proficient form of pMUM, presumably because mycolactone has a key function in M. ulcerans in the environment. It is probably unlikely that the cytotoxic properties of mycolactone for human cells are part of the primary survival strategy of the bacterium. However, given the highly episodic and geographically compact epidemiology of Buruli ulcer, where waves of M. ulcerans infection can rapidly appear and then disappear from a given region, one possibility is that deleterious recombination and loss of the plasmid function interrupt the chain of transmission at some point. Perhaps mycolactone is a factor required for colonization or persistence in insect salivary glands (19) or establishment of a biofilm on plant surfaces (20). In other clonal bacterial pathogens, such as Yersinia pestis, a modest number of genetic changes have led to dramatically different routes of transmission and modes of pathogenesis compared with their progenitors. Indeed, despite their radically different disease pathologies, there are many parallels between Y. pestis and M. ulcerans; in the case of the agent of plague, acquisition of the plasmid-encoded genes ymt and hms has conferred the ability to resist digestion in the midgut of fleas and the ability to persist on the surface of spines that line the interior of the proventriculus, respectively, thus facilitating an arthropod-linked mode of transmission (14, 17). Determining the role(s) of mycolactone and its mode(s) of action is now the subject of investigations in our laboratory.
While the repetitive nature of the mls locus has not yet led to heterogeneity among mycolactones, one DNA deletion identified in this study can be linked with the production of variant toxin. The plasmid gene MUP053 encodes a putative P450 monooxygenase, an enzyme thought to be required for hydroxylation of mycolactone at position C-12' of its fatty acid side chain to produce mycolactone A/B (m/z 765). As predicted, the Australian strain M. ulcerans Chant lacks MUP053 and produces a lower-mass metabolite at m/z 749 (mycolactone C) that corresponds with the absence of a hydroxyl group. The fact that M. ulcerans ITM-941331 from Papua New Guinea also lacks MUP053 but still produces oxidized mycolactones suggests that in some strains there may be chromosomal P450 genes encoding hydroxylases that are active against the molecule.
This study showed that there is considerable mutational dynamism in pMUM. It may be that there is constant genetic flux within the Mls genes such that new mycolactones are continuously being created within a given M. ulcerans population. However, if new metabolites do not confer a fitness advantage, then cells with such changes do not persist. Nevertheless, it seems that there is potential to discover further structural heterogeneity among mycolactones by screening a wider selection of isolates, as demonstrated by the recent discovery of a new mycolactone produced by an M. ulcerans strain from the People's Republic of China (H. Hong et al., Chem. Biochem., in press).
Present address: Department of Biochemistry and Molecular Biology, Monash University, Clayton 3800, Australia. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»