|
|
||||||||


1 Veterinary Laboratories Agency (Weybridge), New Haw, Addlestone, Surrey KT15 3NB, UK.
2 Bacterial Microarray Group, Department of Cellular and Molecular Medicine, St George's Hospital Medical School, Cranmer Terrace, London SW17 0RE, UK
Correspondence
Stephen V. Gordon
stephen.gordon{at}ucd.ie
| ABSTRACT |
|---|
|
|
|---|
, gamma-interferon; PPD, purified protein derivative (tuberculin); qRT-PCR, quantitative real-time polymerase chain reaction
Present address: Emergent Biosolutions, 540–545 Eskdale Road, Winnersh Triangle, Wokingham, Berkshire RG41 5TU, UK.
Present address: School of Agriculture, Food Science and Veterinary Medicine, College of Life Sciences, University College Dublin, Belfield, Dublin 4, Republic of Ireland.
The ArrayExpress (and BµG@Sbase) accession number for the microarray data in this paper is A-BUGS-59.
Two supplementary tables of primers are available with the online version of this paper.
| INTRODUCTION |
|---|
|
|
|---|
Antigens have already been identified that are feeding into the design of new control strategies. The potent T-cell antigens ESAT-6 and CFP-10 were originally identified from Mycobacterium tuberculosis culture filtrates (Sørensen et al., 1995
; Skjøt et al., 2000
), with the genes encoding these antigens deleted from the genome of the BCG vaccine strain (Mahairas et al., 1996
). They therefore have obvious application to the differential diagnosis of vaccination versus infection. Similarly, antigens have been described that impart significant levels of protection against challenge when used as subunit vaccines (Ibanga et al., 2006
; Orme, 2006
; Vordermeier et al., 2006
). However, we need to expand our repertoire of known antigens to ensure that the best candidates, or cocktails, are optimized and applied.
One approach to the discovery of new antigens is to exploit the information contained in the M. bovis genome. Sequencing the M. bovis genome revealed
4000 protein-encoding genes (Garnier et al., 2003
); however, it is not feasible to screen all of these proteins for diagnostic or vaccine potential. Instead a biologically relevant filter must be used to sift the genome information to a manageable subset. One subset of obvious interest would be genes that are upregulated in vivo, but there are difficulties in identifying global expression changes of M. bovis in cattle. A parallel approach would be to identify genes that respond in vitro to a surrogate of the in vivo milieu, such as nutrient starvation, hypoxia or acid shock, and then to screen the identified genes for their expression profile responses in vivo using targeted methods. Following this approach we chose acid as our in vitro surrogate and in doing so found 60 M. bovis genes that were upregulated when the bacteria were exposed to acid shock (Golby et al., 2007
).
In this study, we have determined whether in vitro acid-induced genes are also upregulated in vivo, and then screened the resulting candidates for their immunogenicity in M. bovis-infected cattle. Using this approach, we have identified two putative methyltransferase-encoding genes which are highly upregulated in vivo and whose products are immunogenic in cattle. We furthermore identified the regulator of these methyltransferase genes as a route to understanding their upregulation in vivo.
| METHODS |
|---|
|
|
|---|
Extraction of RNA from mouse tissues and broth cultures.
Female BALB/c mice were intranasally infected with 200–300 c.f.u. of M. bovis and euthanized when clinical signs of tuberculosis disease became apparent. Lungs were removed from the mice, and for each lung one half was placed in a vial containing 10 ml Trizol and the other half in 10 ml phosphate-buffered saline solution (PBS) plus 0.05 % Triton X-100; then both vials were stored at –80 °C. When required, the vials were thawed and the contents homogenized. The Trizol-lung homogenate was transferred to Ribolyser tubes and RNA extracted according to the protocol outlined by Bacon et al. (2004)
. The PBS-lung homogenate was serially diluted and plated onto solid growth medium. The plates were incubated at 37 °C for 3–4 weeks before the numbers of M. bovis colonies were counted. Viable counts were approximately 108–109 c.f.u. per lung.
RNA was prepared from broth cultures of M. tuberculosis and M. bovis using the guanidinium thiocyanate procedure (Bacon et al., 2004
). Quantitative real-time PCR (qRT-PCR) experiments were performed as described by Golby et al. (2007)
; the sequences of the primer pairs are given in Supplementary Table S1, available with the online version of this paper.
Antigen discovery.
A set of 469 peptides spanning the lengths of 11 open reading frames were purchased from Mimitopes. The peptides were 20 residues in length, each with a 12 residue overlap with its neighbouring peptide. Peptides were formulated into pools of approximately 8–12 peptides. The gamma-interferon (IFN-
) immunoassay methodology was as described previously (Cockle et al., 2002
). For screening purposes, bloods were taken from 21 reactor cattle naturally infected with M. bovis. The IFN-
concentration was determined using the BOVIGAM ELISA kit (Prionics). An antigen was defined as giving a positive response when the A450 with antigen minus A450 without antigens was
0.1. Absorbance readings were converted to concentration of IFN-
(pg ml–1) using the following equation: (A450x30.9)+0.5658.
Construction of M. tuberculosis H37Rv Rv1403c–Rv1405c deletion strain PG100.
A 0.9 kb fragment containing the 5' ends of the Mb1440c and fmt genes was PCR amplified using M. bovis 2122 chromosomal DNA and primers Mb1440c5'F1 and Mb1440c5'R. (Sequences of primers used in the construction of plasmids and the verification of the mutant are given in Supplementary Table S2.) The PCR product was digested with BamHI and cloned into the ColE1 plasmid pSMT100, which contains a hygR cassette and the sacB gene, to give the construct pPG28. A 0.9 kb fragment containing the 3' ends of Rv1403c and priA was PCR amplified using primers Mb1438c3'F and Mb1438c3'R1, digested with XbaI and PstI and cloned into pPG28 to give pPG34 (Fig. 1
). The plasmid was used to transform M. tuberculosis H37Rv and transformants were plated onto 7H11 medium containing 0.5 % glycerol, OADC, 2 % sucrose and hygromycin. Hygromycin- and sucrose-resistant colonies were screened for loss of the Rv1403c–Rv1405c genes by PCR using primers that anneal to the 5' end of Rv1405c and the 3' end of Rv1403c. Putative
Rv1403c–Rv1405c mutants were verified by PCR using primers that anneal to internal sequences of the hygromycin-resistance cassette and the fmt gene, hygF and fmtR1 respectively.
|
Microarray analysis.
Three independent experiments (biological replicates) were carried out. For each strain in each experiment, two microarrays (technical replicates) were performed, with each microarray having two measurements of every gene. Two colour hybridizations were performed using whole-genome M. bovis/M. tuberculosis microarrays; the array design is available in BµG@Sbase, accession no. A-BUGS-31 (http://bugs.sgul.ac.uk/A-BUGS-31) and also ArrayExpress, accession no. A-BUGS-31. Cy5 and Cy3 fluorescently labelled probes synthesized from RNA and genomic DNA, respectively, were hybridized to microarrays. Further details concerning the design of the microarrays and procedures used for probe manufacture and hybridization can be found in Golby et al. (2007)
.
Scanning and image analysis.
Microarrays were scanned using a GenePix 4000A microarray scanner (Axon Instruments) with the photomultiplier tube set in the range 550–750 V, so that spots with the highest signal intensities were just below the level of saturation. Fluorescent spots on each image were quantified using BlueFuse for Microarrays software (BlueGnome).
Microarray data analysis.
Data from every microarray were normalized by calculating the log ratio of the Cy5 to Cy3 signal for every spot, and then dividing each log ratio by the median of the log ratios of all spots (excluding control spots) on the array. As an additional normalization step, a median absolute deviation (MAD) scale transformation was applied to the normalized data from the previous step. For each microarray, duplicate spots were averaged, and then the average expression of every gene across all technical replicate microarrays was calculated. Averages of the three biological replicates were used to compare gene expression between strains. For each gene, a moderate t-test was applied and those genes with a P-value less than 0.05 were selected. From this gene list, those genes whose average expression differed by more than 3-fold between strains were selected. Fully annotated microarray data have been deposited in BµG@Sbase (accession number: E-BUGS-59; http://bugs.sgul.ac.uk/E-BUGS-59) and also ArrayExpress (accession number: E-BUGS-59).
| RESULTS |
|---|
|
|
|---|
|
immunoassay using blood obtained from cattle naturally infected with M. bovis. Peptides derived from the sequences of the 11 CDS were synthesized and pooled into 46 groups, with each group consisting of 8–12 overlapping peptides (Table 2
production measured after 48 h of culture. Fig. 2
(Fig. 3
|
|
|
Rv1403c–Rv1405c : : hyg mutant
Liquid culture growth experiments using PG100 showed that it grew more slowly, and had a much greater tendency to aggregate, than the wild-type (Fig. 4
). The defective growth phenotype of the mutant could only be partially corrected by transformation with a plasmid (pPG72) that contained the genes Rv1405c, Rv1404 and Rv1403c. Transformation of the mutant with a plasmid overexpressing Rv1404 alone (pPG57) conferred a similar level of complementation to that seen with pPG72 (Fig. 4
). Complementation with Rv1404 had a similar effect to addition of all three deleted genes, suggesting that the PG100 growth defect was due to the loss of Rv1404 and not Rv1405c or Rv1403c. However, failure to achieve full trans complementation using the Rv1404 overexpression construct suggests a complex regulatory circuit that may require cis complementation.
|
Rv1404 background is approximately 70–150-fold higher than in the wild-type. Evidence to support the proposed repressor role of Rv1404 on the expression of Rv1405c was provided by the finding that the expression of Rv1405c in PG100 carrying the plasmid pPG55, which contains both Rv1405c and Rv1404, is 45-fold less than the mutant carrying pPG53. Expression of Rv1403c was also found to be 37.7-fold higher in PG100 carrying the Rv1403c-containing pPG54 compared to the wild-type, but only 3.4-fold higher in PG100 carrying both Rv1403c and Rv1404 (pPG56). This suggests that Rv1404 represses the transcription of both Rv1405c and Rv1403c, but the lower fold induction levels of Rv1403c in the mutant would suggest that the promoter of Rv1403c is weaker than that of Rv1405c.
|
Rv1405c–Rv1403c triple mutant could be partially reversed by complementing with the Rv1404 regulator, suggesting that Rv1404 regulated genes in addition to Rv1405c and Rv1403c, genes that may include the methylation target of Rv1403c and Rv1405c. In order to identify these genes, DNA microarrays were used to compare the transcriptional profiles of the wild-type H37Rv and PG100. Table 4
|
|
| DISCUSSION |
|---|
|
|
|---|
Two of the most highly upregulated genes in response to in vitro acid shock, the methyltransferase-encoding genes Rv1403c/Mb1438c and Rv1405c/Mb1440c, were also found to be highly upregulated in vivo. The expression of Mb1440c in vivo was particularly high, showing an approximately 350-fold higher level of expression compared to an in vitro control. This level of induction is comparable to the 70–150-fold induction of Rv1405c seen in an
Rv1404 mutant background compared to the wild-type, suggesting that Rv1405c/Mb1440c expression was completely derepressed in infected mouse lung tissues.
Expression of Mb1440c/Rv1405c and Mb1438c/Rv1403c was shown to be tightly regulated by the product of the Mb1439/Rv1404 gene, encoding a member of the MarR family of transcriptional regulators (Ellison & Miller, 2006
; Grkovic et al., 2002
). Mb1440c/Rv1405c and Mb1438c/Rv1403c showed very low expression under in vitro non-stressed and ex vivo conditions, but high expression under in vivo and stressed in vitro conditions. It is noteworthy that in Mycobacterium ulcerans a transposon is inserted in the promoter/operator region of the Mb1439/Rv1404 orthologue (Stinear et al., 2007
), presumably resulting in constitutive expression of the Rv1405c orthologue since its protein product was detected in the cytoplasmic fraction of M. ulcerans by LC-MS (annotated at http://genolist.pasteur.fr/BuruList/).
The majority of MarR regulators act as repressors, and their activity can be modulated through binding of an inducer molecule to the MarR regulator (Ellison & Miller, 2006
; Grkovic et al., 2002
). Genes regulated by MarR family members are repressed in the absence of the inducer and derepressed in the presence of inducer. Many MarR repressors are also autoregulatory, repressing their own expression in the absence of inducer. This would explain the observed moderate increase (3–4-fold) in expression of Mb1439/Rv1404 observed under acidic and in vivo conditions, as the regulator would be unable to bind to its operator and repress its own expression. In addition to Rv1405c and Rv1403c, deletion of Rv1404 was also shown to affect the expression of another 25 genes, 10 of which were previously shown to be upregulated in response to acid shock in M. tuberculosis H37Rv. This suggests that Rv1404/Mb1439 is an important regulator in the response of M. tuberculosis and M. bovis to acid shock. Curiously, the three genes that show the highest induction in the Rv1404 mutant (Rv0193c–Rv0195) show no upregulation in response to acid (Golby et al., 2007
), suggesting that the expression of these genes could be regulated by other factors in addition to Rv1404.
In a seminal study, Sassetti and Rubin used a combination of saturation mutagenesis with in vivo selection to identify M. tuberculosis genes that were required for survival in a mouse infection model; Rv1405c was among the 197 genes identified as important for virulence (Sassetti & Rubin, 2003
). Hence, methylation of Rv1405c's target is involved in virulence. There are a number of examples where methylation of mycobacterial cellular components has been implicated in the virulence of M. tuberculosis. Methylation of the mycobacterial heparin-binding haemagglutinin adhesin (HBHA), an extrapulmonary dissemination factor (Pethe et al., 2001
), is required for T-cell immunity and protects the molecule against proteolysis (Pethe et al., 2002
; Temmerman et al., 2004
). Similarly, cyclopropanation of mycolic acids by pcaA and cmaA2 has been shown to be involved in virulence of M. tuberculosis (Glickman et al., 2000
; Rao et al., 2006
). Overexpression of Rv1405c in M. smegmatis mc2155 caused no significant changes in mycolic acid profile (data not shown), suggesting that Rv1405c is not involved in cyclopropanation of mycolic acids. Analysis of the Rv1405c locus reveals genes encoding formyl methionine transferase (fmt) and a putative rRNA methyltransferase (fmu), suggesting that Rv1405c/Rv1403c could be involved in methylating a component of the protein translational machinery; however, how this would play a role in virulence is unknown.
The
Rv1403c-Rv1404-Rv1405c mutant also showed upregulation of distal genes. Rv0193c–Rv0195 showed the highest degree of upregulation (13–52-fold), and their contiguity suggests that their products could be functionally related. Rv0195 encodes a response regulator of the two-component family and, like Rv1403c and Rv1405c, has been shown to be upregulated in response to low oxygen conditions (Muttucumaru et al., 2004
; Voskuil et al., 2004
). The gene Rv0194 encodes a putative ABC-type transporter, but its function and substrate are unknown. The gene is unique amongst other ABC-type transporters of M. tuberculosis in having two membrane-spanning domains and two nucleotide-binding domains in a single polypeptide.
In summary, on an applied level we have defined novel M. bovis antigens that are now being moved forward to larger field trials. On a fundamental level, we have described the regulation of genes encoding two methyltransferases in M. tuberculosis, one of which, Rv1405c, has been implicated as a virulence factor. The next step will be to identify the target(s) of these methyltransferases to elucidate their role in virulence.
| ACKNOWLEDGEMENTS |
|---|
Edited by: G. R. Stewart
| REFERENCES |
|---|
|
|
|---|
Cockle, P. J., Gordon, S. V., Lalvani, A., Buddle, B. M., Hewinson, R. G. & Vordermeier, H. M. (2002). Identification of novel Mycobacterium tuberculosis antigens with potential as diagnostic reagents or subunit vaccine candidates by comparative genomics. Infect Immun 70, 6996–7003.
Cole, S. T., Brosch, R., Parkhill, J., Garnier, T., Churcher, C., Harris, D., Gordon, S. V., Eiglmeier, K., Gas, S. & other authors (1998). Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature 393, 537–544.[CrossRef][Medline]
Donnelly-Wu, M. K., Jacobs, W. R., Jr & Hatfull, G. F. (1993). Superinfection immunity of mycobacteriophage L5: applications for genetic transformation of mycobacteria. Mol Microbiol 7, 407–417.[CrossRef][Medline]
Ellison, D. W. & Miller, V. L. (2006). Regulation of virulence by members of the MarR/SlyA family. Curr Opin Microbiol 9, 153–159.[CrossRef][Medline]
Garnier, T., Eiglmeier, K., Camus, J. C., Medina, N., Mansoor, H., Pryor, M., Duthoy, S., Grondin, S., Lacroix, C. & other authors (2003). The complete genome sequence of Mycobacterium bovis. Proc Natl Acad Sci U S A 100, 7877–7882.
Glickman, M. S., Cox, J. S. & Jacobs, W. R., Jr (2000). A novel mycolic acid cyclopropane synthetase is required for cording, persistence, and virulence of Mycobacterium tuberculosis. Mol Cell 5, 717–727.[CrossRef][Medline]
Golby, P., Hatch, K. A., Bacon, J., Cooney, R., Riley, P., Allnutt, J., Hinds, J., Nunez, J., Marsh, P. D. & other authors (2007). Comparative transcriptomics reveals key gene expression differences between the human and bovine pathogens of the Mycobacterium tuberculosis complex. Microbiology 153, 3323–3336.
Grkovic, S., Brown, M. H. & Skurray, R. A. (2002). Regulation of bacterial drug export systems. Microbiol Mol Biol Rev 66, 671–701.
Ibanga, H. B., Brookes, R. H., Hill, P. C., Owiafe, P. K., Fletcher, H. A., Lienhardt, C., Hill, A. V., Adegbola, R. A. & McShane, H. (2006). Early clinical trials with a new tuberculosis vaccine, MVA85A, in tuberculosis-endemic countries: issues in study design. Lancet Infect Dis 6, 522–528.[CrossRef][Medline]
Mahairas, G. G., Sabo, P. J., Hickey, M. J., Singh, D. C. & Stover, C. K. (1996). Molecular analysis of genetic differences between Mycobacterium bovis BCG and virulent M. bovis. J Bacteriol 178, 1274–1282.
Muttucumaru, D. G., Roberts, G., Hinds, J., Stabler, R. A. & Parish, T. (2004). Gene expression profile of Mycobacterium tuberculosis in a non-replicating state. Tuberculosis (Edinb) 84, 239–246.[CrossRef][Medline]
Orme, I. M. (2006). Preclinical testing of new vaccines for tuberculosis: a comprehensive review. Vaccine 24, 2–19.[CrossRef][Medline]
Pethe, K., Alonso, S., Biet, F., Delogu, G., Brennan, M. J., Locht, C. & Menozzi, F. D. (2001). The heparin-binding haemagglutinin of M. tuberculosis is required for extrapulmonary dissemination. Nature 412, 190–194.[CrossRef][Medline]
Pethe, K., Bifani, P., Drobecq, H., Sergheraert, C., Debrie, A. S., Locht, C. & Menozzi, F. D. (2002). Mycobacterial heparin-binding hemagglutinin and laminin-binding protein share antigenic methyllysines that confer resistance to proteolysis. Proc Natl Acad Sci U S A 99, 10759–10764.
Rao, V., Gao, F., Chen, B., Jacobs, W. R., Jr & Glickman, M. S. (2006). Trans-cyclopropanation of mycolic acids on trehalose dimycolate suppresses Mycobacterium tuberculosis-induced inflammation and virulence. J Clin Invest 116, 1660–1667.[CrossRef][Medline]
Sassetti, C. M. & Rubin, E. J. (2003). Genetic requirements for mycobacterial survival during infection. Proc Natl Acad Sci U S A 100, 12989–12994.
Skjøt, R. L., Oettinger, T., Rosenkrands, I., Ravn, P., Brock, I., Jacobsen, S. & Andersen, P. (2000). Comparative evaluation of low-molecular-mass proteins from Mycobacterium tuberculosis identifies members of the ESAT-6 family as immunodominant T-cell antigens. Infect Immun 68, 214–220.
Sørensen, A. L., Nagai, S., Houen, G., Andersen, P. & Andersen, A. B. (1995). Purification and characterization of a low-molecular-mass T-cell antigen secreted by Mycobacterium tuberculosis. Infect Immun 63, 1710–1717.
Stinear, T. P., Seemann, T., Pidot, S., Frigui, W., Reysset, G., Garnier, T., Meurice, G., Simon, D., Bouchier, C. & other authors (2007). Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans, the causative agent of Buruli ulcer. Genome Res 17, 192–200.
Temmerman, S., Pethe, K., Parra, M., Alonso, S., Rouanet, C., Pickett, T., Drowart, A., Debrie, A. S., Delogu, G. & other authors (2004). Methylation-dependent T cell immunity to Mycobacterium tuberculosis heparin-binding hemagglutinin. Nat Med 10, 935–941.[CrossRef][Medline]
Vordermeier, H. M., Chambers, M. A., Buddle, B. M., Pollock, J. M. & Hewinson, R. G. (2006). Progress in the development of vaccines and diagnostic reagents to control. tuberculosis in cattle. Vet J 171, 229–244.[CrossRef][Medline]
Voskuil, M. I., Visconti, K. C. & Schoolnik, G. K. (2004). Mycobacterium tuberculosis gene expression during adaptation to stationary phase and low-oxygen dormancy. Tuberculosis (Edinb) 84, 218–227.[CrossRef][Medline]
Received 6 November 2007;
revised 18 December 2007;
accepted 3 January 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |