|
|
||||||||
1 Department of Pathology and Laboratory Medicine, University of British Columbia, 950 West 28th Avenue, Vancouver, BC V5Z 4H4, Canada
2 Department of Pathology and Laboratory Medicine, University of British Columbia, Department of Cancer Genetics, BC Cancer Research Centre, 601 West 10th Avenue, Vancouver, BC V5Z 1L3, Canada
3 Departments of Paediatrics and Pathology and Laboratory Medicine, University of British Columbia, Division of Infectious and Immunological Diseases, BC Children's Hospital, 950 West 28th Avenue, Vancouver, BC V5Z 4H4, Canada
Correspondence
Richard W. Stokes
rstokes{at}interchange.ubc.ca
| ABSTRACT |
|---|
|
|
|---|
, bone-marrow-derived macrophages; FRD, fumarate reductase; MNPO, mercaptopyridine-N-oxide; qPCR, quantitative real-time PCR; SDH, succinate dehydrogenase; TB, tuberculosis| INTRODUCTION |
|---|
|
|
|---|
The bacterial model applied in this study is that of M. tuberculosis H37Ra and H37Rv, sibling strains derived from the parental H37 (Oatway & Steenken, 1936
; Steenken et al., 1934
), a laboratory strain originally isolated from the sputum of a pulmonary TB patient in 1906 (Steenken & Gardner, 1946
). Strains H37Ra and H37Rv are highly similar, as studies have shown little genomic difference between them (Bhargava et al., 1990
; Collins & De Lisle, 1984
; Imaeda, 1985
), save for the RvD2 region of difference (Brosch et al., 1999
; Lari et al., 2001
), which has not been conclusively linked to pathogenesis (Lari et al., 2001
). However, in spite of this lack of genomic differences between the strains, H37Rv has been shown to display a marked advantage in survival and growth in animal hosts over that of H37Ra (Collins & Smith, 1969
; Jung et al., 2002
). Thus, dissecting the molecular aspects of the phenotypic differences between these highly related strains may lead to a better understanding of mycobacterial virulence.
Studies comparing the expression differences between broth-grown H37Ra and H37Rv and/or intracellular and broth-grown H37Ra/Rv bacteria have elucidated genes differentially expressed between the strains, some of which have been shown to have roles in mycobacterial virulence (Graham & Clark-Curtiss, 1999
; Kinger & Tyagi, 1993
; Rindi et al., 1999
, 2001
; Rivera-Marrero et al., 1998
; Wei et al., 2000
; Gao et al., 2004
). For example, Gao et al. (2004)
examined, via microarray, genes involved in the cording and non-cording phenotypes of the respective H37Rv and H37Ra. Graham & Clark-Curtiss (1999)
developed the novel SCOTS technique to minimize the noise introduced by rRNA populations and, in their comparisons of the transcriptomes of broth-grown and intracellular M. tuberculosis H37Rv, they were able to isolate several genes upregulated in intracellular M. tuberculosis. One of these genes was pks2, a putative polyketide synthase which has been found to play a role in the synthesis of sulpholipid (Sirakova et al., 2001
), an important component of the mycobacterial cell wall. Furthermore, Kinger & Tyagi (1993)
examined, via differential display, genes differentially expressed between aerobic cultures of M. tuberculosis H37Ra and H37Rv. From this study, the devR/devS (also known as dosR/dosS) two-component system was identified, which has subsequently been found to have a role in the survival of M. tuberculosis under anoxic conditions (Dasgupta et al., 2000
; Guinn et al., 2004
; Park et al., 2003
; Saini et al., 2004
; Wei et al., 2000
). These studies highlight the value of the H37Rv/H37Ra model in the search for genes that may have a role in M. tuberculosis pathogenesis.
Our study also utilized the H37Rv/H37Ra model; however, we reasoned that gene expression changes that would affect bacterial survival and growth within host cells are more likely to be discovered by focusing on differences between H37Ra and H37Rv expression when both strains are in interaction with macrophages rather than by comparing broth-grown with intracellular bacteria. To assess these transcriptome changes between strains, we used a novel technique – bacterial artificial chromosome fingerprint array (BACFA) – that has been successfully utilized for comparative genomic studies (Brosch et al., 1998
, 1999
; Gordon et al., 1999a
). We also provide evidence that, through use of BACFA, we have identified an enzyme complex, fumarate reductase, that aids the survival of intracellular M. tuberculosis.
| METHODS |
|---|
|
|
|---|
Isolation and culturing of bone-marrow-derived macrophages (BM-M
).
BM-M
were obtained from the femora, tibiae and humeri of CD-1 mice as previously described (Rooyakkers & Stokes, 2005
; Stokes et al., 1993
). Briefly, the bones were dissected out from euthanized, 7–10-week-old female CD-1 mice, and the marrow was flushed from the bones using BM-M
medium [RPMI containing 10 % heat-inactivated fetal calf serum, 10 % L-929 cell medium (Rooyakkers & Stokes, 2005
), 2 mM L-glutamine and 1 mM sodium pyruvate]. The cell suspension was left to adhere in a 175 cm2 flask for 3 h to deplete non-stem cells. Non-adherent cells (stem cells) were collected and cultured for 7 days in BM-M
medium at 37 °C with 5 % CO2 to allow differentiation into macrophages in 150 cm2 tissue culture flasks for RNA studies, or on 13 mm diameter coverslips for c.f.u. determinations (Rooyakkers & Stokes, 2005
).
Infection and determination of macrophage-associated growth of M. tuberculosis.
To assess bacterial expression profiles 96 h and 168 h post-infection, monolayers were incubated at an m.o.i. of 10 bacteria to 1 macrophage, which resulted in an average infection rate (average c.f.u. at 4 h per macrophage plated in the well) of 0.1 bacteria per macrophage. After a 4 h incubation with binding medium (138 mM NaCl, 8.1 mM Na2HPO4, 1.5 mM KH2PO4, 2.7 mM KCl, 0.6 mM CaCl2, 1.0 mM MgCl2 and 5.5 mM D-glucose) containing the requisite number of bacteria, the monolayers were washed three times with pre-warmed medium to rinse off unbound bacteria, before being submerged in cRPMI (RPMI 1640 with 10 % fetal calf serum, 2 mM L-glutamate and 1 mM sodium pyruvate) until they were to be processed for RNA at days 4 and 7 post-infection. To assess bacterial expression at 4 h post-infection, monolayers were infected at an m.o.i. of 60 bacteria to 1 macrophage, resulting in an average infection rate of 0.6 bacteria per macrophage. This higher m.o.i. was used to allow for the harvest of a sufficient amount of bacterial RNA and was comparable to the bacterial load in macrophages 4 days post-infection at an m.o.i. of 10 : 1 (data not shown). At each time point, coverslips and supernatants were briefly sonicated (10 s) to release and disperse bacteria, and plated on 7H10 agar plates supplemented with 10 % OADC to determine c.f.u. per ml (Rooyakkers & Stokes, 2005
). For all macrophage experiments, three replicate coverslips were assessed at each time point for each of three independent experiments.
Extraction and purification of bacterial RNA.
RNA was extracted from broth cultures as previously described (Mangan et al., 1997
; Stewart et al., 2002
). For intracellular bacteria used to examine M. tuberculosis expression profiles at 96 and 168 h post-infection, 50 ml GTC lysis buffer (Monahan et al., 2001
) was added to each flask at the designated time and rocked gently until Schlieren lines disappeared, indicating homogeneous mixing of the solutions. For intracellular bacteria used to examine mycobacterial expression at 4 h post-infection, flasks were washed three times with pre-warmed medium, and 50 ml GTC lysis buffer was poured directly onto the monolayer. The entire contents of flasks were poured into 50 ml conical tubes, and spun at 3700 r.p.m. (3500 g) for 15 min. Pellets were resuspended and pooled with 1 ml GTC lysis buffer. RNA was extracted from the pellets and then purified as previously described (Mangan et al., 1997
; Stewart et al., 2002
; Wernisch et al., 2003
).
Reverse transcription and cDNA synthesis.
For BAC arrays, 1–5 µg total RNA was reverse transcribed with Superscript II (Invitrogen) and a universal arbitrary primer (5'-GCCGGAGCTCTGCAGAATTC-3'), henceforth called Uniprime, was used to generate single-stranded cDNA (Mangan & Butcher, 1998
). Eight units of Klenow fragment of DNA polymerase I (Invitrogen) was then used to synthesize second-strand cDNA according to previously published protocols (Mangan & Butcher, 1998
).
DIG-labelling of cDNA for use in BACFA.
A 2 µl aliquot of double-stranded cDNA was added to 8 µl water, boiled for 10 min and added to the following: 2 µl 10x PCR buffer (Invitrogen), 2.5 µl 2 mM DIG-dUTP : dTTP (3 : 1 ratio), 2.5 µl 2 mM dNTPs (dATP, dCTP, dGTP), 0.8 µl 50 µM Uniprime stock solution, 0.6 µl 50 mM MgCl2 and 1 U Taq Polymerase (Invitrogen), to a final volume of 20 µl. DIG labelling of cDNA was then done via PCR (94 °C for 2 min followed by 30 cycles of 94 °C, 55 °C and 72 °C for 1, 2 and 3 min, respectively).
Generation of bacterial artificial chromosome fingerprint arrays (BACFAs).
The M. tuberculosis H37Rv BAC library was kindly provided by Drs Stewart Cole and Roland Brosch at the Pasteur Insitute (Paris, France), who used it to isolate and identify genomic differences between H37Ra and H37Rv (Brosch et al., 1998
, 1999
; Gordon et al., 1999b
). This library contains 78 BACs, each containing an average of 68 kb of genomic sequence from H37Rv. For preliminary BACFA assessment of differential expression between H37Ra and H37Rv, all BACs were digested with 10 U each of PvuII and StuI (New England Biolabs) overnight at 37 °C. After preliminary screening, some BACs were digested with StuI and SalI. The BACs of interest were first digested with SalI overnight at 37 °C; the products were cleaned with a Qiagen PCR purification column and eluted with 20 µl water. Eluted volumes were then digested overnight at 37 °C with 10 U StuI.
Digestion products were run on a 1.5 % TAE-agarose gel for 2 h at 80 V. Southern transfer of fragments to nylon membrane was done by the manufacturer's protocol (Roche). Fragments were immobilized onto nylon membranes by UV cross-linking at 260 nm for 5 min followed by baking at 80 °C for 2 h.
Generation of BACFA hybridization profiles.
Nylon membranes were hybridized according to the manufacturer's instructions with 2.5 µl DIG-labelled PCR reaction product (above) per ml hybridization buffer. Membranes were hybridized for 16–18 h at 50 °C. Immunological detection of blots was carried out by the manufacturer's protocols (Roche). Profiles were generated in duplicate for each of three independent pools of cDNA from H37Ra and H37Rv. Differences were only singled out if they appeared in all six profiles.
Quantitative real-time PCR (qPCR) analysis of expression differences observed with BACFA.
As hybridization protocols could potentially allow for non-specific binding resulting in a false positive, all differences identified via BACFA were assessed via qPCR. Using gene-specific primers, this technique allows for unambiguous expression analysis of the candidates. Second-strand cDNA generated with the universal primer was used in qPCRs with SYBR Green as the indicator dye. For qPCRs, primers were designed with PrimerQuest, a web-based program freely available from the IDT website (http://www.idtdna.com/Scitools/Applications/Primerquest/), as well as Primer Software from Molecular Biology Tools. Primers were designed with an annealing temperature of 57 °C, and used at a final concentration of 300 nM in a reaction volume of 20 µl. qPCR cycling conditions were as follows: 95 °C for 10 min followed by 35 cycles of 94 °C for 30 s, 57 °C for 20 s, and 72 °C for 30 s. qPCR data are presented as fold difference of expression in H37Ra over that in H37Rv using the 2–
Ct method (Schmittgen et al., 2000
), with rrnAP1(Menendez Mdel et al., 2005
) and dnaK used as normalizing genes for the 96 h and 168 h transcripts, and rrnAP1 and 16S used as normalizing genes for 4 h transcripts. 16S was used in addition to rrnAP1 at 4 h post-infection as it was found to be more reliable at the earlier time point compared to the later time point. Two normalizing genes were used for each time point as previous reports have shown that normalizing genes do not necessarily stay constant over time and that more than one is needed for reliable expression analysis (Vandesompele et al., 2002
).
Treatment of extracellular and intracellular bacteria with mercaptopyridine-N-oxide (MPNO).
We assessed bacterial viability in the presence of MPNO, a putative inhibitor of fumarate reductase. Frozen aliquots of M. tuberculosis H37Ra and H37Rv were thawed, and diluted to 6x106 ml–1 in 7H9. Aliquots (200 µl) of the bacterial suspension were then dispensed into 96-well plates and incubated at 37 °C for 72 h to produce actively growing bacteria. At 72 h, stock MPNO (10 mM in water) was added to wells to give a final concentration ranging from 0 to 4.8 µM and incubated for a further 72 h. Then 50 µl Alamar Blue reagent was added to each well; the plate was left at 37 °C for 24 h and read at both 550 nm and 570 nm on a Bio-Rad microplate reader. Analysis of Alamar Blue reduction as a measure of bacterial viability was done according to the manufacturer's protocol (TREK Diagnostic Systems, Fisher Biosciences).
To monitor the effects of MPNO on intracellular mycobacteria, monolayers of BM-M
were infected as described above. After infection, the monolayers were washed thoroughly to remove unbound bacteria and three replicate coverslips were removed to assess the numbers of bacteria bound/internalized by macrophages. Half of the remaining coverslips were overlaid with cRPMI supplemented with 2.4 µM MPNO whilst the other half were overlaid with unsupplemented cRPMI. Coverslips and supernatants were processed at 96 and 168 h post-infection as described above to monitor intracellular bacterial growth.
Statistics.
Statistical significance of comparisons between H37Rv and H37Ra, or between two time points, was determined with the two-tailed, unpaired Student's t test. P values <0.05 were considered significant.
| RESULTS |
|---|
|
|
|---|
(Fig. 1c
the doubling time of H37Rv was 31.2±1.3 h, whereas the doubling time of H37Ra was 51.1±1.4 h (P=0.014). While not significant, a marked difference in growth between strains was observed on day 4 (96 h) post-infection (for H37Ra and for H37Rv, P=0.1058, Fig. 1c
|
|
|
|
|
Assessment of candidate gene expression profiles in broth cultures and at 4 h and 96 h post-infection
The 168 h post-infection time point was chosen primarily because that was the time at which sufficient RNA could be harvested for BACFA analysis. However, there remained the possibility that the significant difference we saw between H37Ra and H37Rv at 168 h post-infection was independent of events at earlier time points. Taking candidate genes identified as having expression differences at 168 h post-infection, the expression profiles of these genes were also assessed at 4 and 96 h post-infection (Fig. 4a
). Expression profiles of the genes were also assessed for PB+T broth cultures to determine whether the inter-strain differences are inherent or are more pronounced in interactions with the host macrophage (Fig. 4b, c
).
|
Expression of all frd operon genes was elevated in bacteria interacting with macrophages versus bacteria grown in enriched broth. This was manifested as a substantial increase in expression at 4 h that declined over the 168 h period of the experiment for both H37Rv and H37Ra (Fig. 4b, c
).
The upregulation of pks2 in macrophage-associated M. tuberculosis over that of broth-grown bacteria was maintained in H37Rv from initial interaction throughout the infection (Fig. 4c
). This trend was not detected in macrophage-associated H37Ra at 4 h post-infection, but the expression of pks2 was observed to increase over the duration of the infection (Fig. 4b
).
aceE displayed higher levels of expression in macrophage-associated bacteria versus broth-grown bacteria only at 168 h post-infection for both H37Ra and H37Rv (Fig. 4b, c
); however, levels of aceE transcript were higher in intracellular H37Rv versus intracellular H37Ra at all time points assessed during the infection (Fig. 4a
).
Expression of Rv1571 was upregulated in intracellular H37Rv over H37Ra at both 4 and 168 h post-infection, but in broth, expression was 16-fold higher in H37Ra (Fig. 4a
). Interestingly, the expression of Rv1571 was greatly enhanced in both H37Rv and H37Ra on exposure to the intracellular environment (Fig. 4b, c
).
Lastly, statistical evaluation of gene expression revealed a significant difference between H37Ra and H37Rv for frdC at 96 h post-infection (P=0.0068). However, statistical significance was not noted for the other genes.
Growth of M. tuberculosis under oxygen-limiting conditions
As FRD is an important enzyme in fumarate respiration, an energy-production pathway relied upon by bacteria exposed to anoxic environments, we assessed the growth of H37Ra and H37Rv under oxygen-limiting conditions similar to those described by Wayne (1976)
. Doubling times of H37Ra and H37Rv grown under these conditions (40.7±1.3 h for H37Ra; 45.9±7.9 h for H37Rv) differed significantly from either strain grown in aerated roller bottle cultures (23.49±2.4 h for H37Ra P=0.0033; 20.3±1.1 h for H37Rv, P=0.0325). However, the growth rates of the oxygen-limited cultures of H37Ra and H37Rv were not significantly different (P=0.297). All doubling times were calculated at day 5 (120 h).
qPCR analysis of frd transcripts in oxygen-limited cultures of M. tuberculosis
Concurrent with growth rate assessments of oxygen-limited cultures, RNA from these and their aerated culture counterparts were also harvested and levels of the frd transcripts were assessed. Transcripts from static cultures used to inoculate the low-oxygen and aerophilic cultures –henceforth referred to as seed cultures – were also assessed to determine if changes after inoculation were pre-existing and strain specific. As seen in the previous experiment (Fig. 4a
), seed cultures of M. tuberculosis H37Ra displayed higher levels of the frd transcripts compared to seed cultures of H37Rv (Fig. 5a
). Comparable to the results with intracellular bacteria, frdA and frdB were expressed at a higher level in H37Rv shortly after the oxygen constraints were initiated but, thereafter, H37Ra expression of these genes increased to levels higher than in H37Rv (Fig. 5a
). Unlike the results obtained for intracellular bacteria, transcription of frdC and frdD was higher in H37Ra over the period of this study (Fig. 5a
). Comparing the frd transcripts in oxygen-limited cultures versus their respective aerated culture counterparts, frd genes were expressed rapidly in oxygen-limited cultures of the virulent H37Rv, with a subsequent fall-off in expression levels over time (Fig. 5c
). Expression of frdC in H37Ra followed a similar pattern but, in contrast, expression of frdA, B and D was initially low, gradually increasing over the 96 h of the experiment. Using Student's t-test, only the expression of frdA and frdB at 4 h was significantly different between strains (P=0.036, 0.013, respectively) whereas all other transcripts were not significantly different between strains.
|
|
| DISCUSSION |
|---|
|
|
|---|
Firstly, the primer selected for the generation of DIG-labelled probes was of a defined sequence, which has the potential to limit the differences detected. However, in our experience random primers that would enable the screening of a more complex pool of transcripts did not allow for the generation of probes that could be used in expression analysis. Thus a compromise would be to use a defined primer of similar length, but different sequence, allowing for a different snapshot of the transcriptome. A further possibility would be to use genome-directed primers (Talaat et al., 2000
) in separate RT and PCR applications, with pools of these probes enabling a more complex analysis of the differences between strains. For the present study, the arbitrary primers were used to elucidate the ability to generate meaningful BACFAs for M. tuberculosis; subsequent studies will make use of genome-directed primers.
A second consideration that affects the use of BACFAs in expression analyses is that one signal may be the result of several genes. It is this complication that will require resolution by generating alternative BACFAs with different sets of enzymes, which may be time-consuming. Conversely, this limitation is also an advantage of the BACFA method as it will allow for the identification of genes present on the same fragment that may be co-regulated. As in the case of the frd operon isolated in this study, each gene on its own may not have given a signal intense enough to warrant further analysis. However, as three genes of the upregulated operon were present in one fragment, their combined signals led to the isolation of this operon. Ultimately, BACFA analyses can be a viable alternative to, if not a substitute for, other expression analysis techniques, particularly if cost and specialized facilities are of concern.
Using the BACFA technique, both previously reported and novel differences were found. Differences in expression of pks2 (a gene whose product is thought to be involved in sulpholipid synthesis as well as other lipid synthesis roles: Sirakova et al., 2001
) have been previously reported in transcriptome comparisons of broth-grown versus macrophage-associated H37Rv (Graham & Clark-Curtiss, 1999
). Similarly, we found that pks2 expression was upregulated in intracellular H37Rv and H37Ra versus their broth-grown counterparts. Additionally, we report, apparently for the first time, that pks2 expression is higher in intracellular H37Rv versus H37Ra.
Previously unreported differences in gene expression identified using BACFA and subsequently confirmed via qPCR analyses were those of the genes Rv1571, aceE and frdD. aceE, seen to be expressed more readily in the virulent M. tuberculosis H37Rv throughout the course of infection, encodes pyruvate decarboxylase, otherwise known as the E1 subunit for the pyruvate dehydrogenase complex (PDC). The PDC catalyses the conversion of pyruvate to acetyl-CoA, which feeds into the tricarboxylic acid cycle, one of the main pathways of cellular respiration and biosynthesis in both eukaryotes and prokaryotes (Madigan et al., 1997a
). However, this E1 subunit can act independently of the PDC in glycolysis, converting pyruvate into acetaldehyde and a molecule of CO2 (Madigan et al., 1997a
). Previous studies have shown that H37Rv displays a higher respiration rate than H37Ra (Heplar et al., 1954
) and that H37Ra and H37Rv rely on both glycolytic and oxidative means of glucose metabolism, although H37Rv was seen to rely more heavily on glycolytic pathways (Ramakrishnan et al., 1962
). aceE upregulation could help to support both aerobic and glycolytic pathways, acting either in concert with the rest of the PDC in aerobic respiration or independently as pyruvate decarboxylase under low-oxygen growth.
The dramatic increase in Rv1571 expression on transfer of both H37Rv and H37Ra from broth to the intracellular environment is of interest and warrants further investigation. However, there was no consistent pattern of expression differences for this gene between intracellular H37Rv and H37Ra.
An additional difference found with BACFA analysis and confirmed with subsequent qPCR analysis was that of the components of the fumarate reductase (FRD) complex (frdA, B, C and D). This enzyme complex catalyses the conversion of fumarate into succinate, and is composed of four subunits: FRD-A and FRD-B, which constitute the catalytic domain, and FRD-C and FRD-D, which constitute the anchoring domain. Fumarate respiration is an alternative means to acquire energy by utilizing fumarate as the terminal electron acceptor when oxygen or
is absent. BACFA analysis of 168 h post-infection transcripts revealed an upregulation of frdB, frdC and frdD in H37Ra versus H37Rv. It was previously found in E. coli that FRD can partially compensate for missing succinate dehydrogenase (SDH) activity (Guest, 1981
), a critical enzyme of the tricarboxylic acid cycle; however, PCR analysis of the sdh operon in both H37Ra and H37Rv showed no genomic or expression differences between strains (Fig. 3b
).
In mycobacterial studies, frdA has been found to be upregulated in M. tuberculosis interacting with the macrophage versus M. tuberculosis grown in broth cultures (Schnappinger et al., 2003
) as well as in studies that examined the behaviour of M. tuberculosis grown under carbon starvation (Betts et al., 2002
). Microarray studies examining the transcriptome in stationary-phase M. tuberculosis H37Rv cultures found an increase of frdB and frdC transcripts versus bacteria in exponential growth (Raman et al., 2001
). Additionally, investigations into the respiratory behaviour of Mycobacterium phlei found that FRD activity increased fourfold when bacteria were grown under low-oxygen conditions (Gillespie et al., 1988
). Our studies have taken these initial observations further by including a direct comparison between H37Rv and H37Ra in the intracellular environment as well as in broth culture. The recently published genomic sequence of M. tuberculosis H37Ra allowed the direct comparison of frd gene sequences of the two strains (Wang et al., 2007
). However, no differences were noted in either the gene or promoter sequences of the frd operon in H37Ra and H37Rv (Cole et al., 1998
; Wang et al., 2007
). When we examined the expression of the frd operon over the duration of the infection period at 4, 96 and 168 h post-infection, the trend observed was one where expression was initiated and climaxed earlier in the virulent H37Rv strain versus that in the attenuated H37Ra. These results support a hypothesis where H37Rv is able to respond more quickly to the challenge of the intracellular environment – possibly in response to the oxygen-limited conditions encountered inside the macrophage (Schnappinger et al., 2003
), thereby resulting in enhanced replication within the macrophage.
To address the role of FRD under anoxic conditions, we sought to characterize FRD activity and frd transcripts in cultures of M. tuberculosis H37Rv and H37Ra grown as both unaerated static (oxygen-limited) and aerated roller-bottle cultures. The attenuated strain grown under low-oxygen conditions displayed a lag in gene expression with regard to the catalytic domain of FRD (frdA, frdB) but not the anchoring domain (frdC, frdD). The data do not exactly reproduce the trends seen with the intracellular bacteria, probably due to the multiple challenges encountered by the bacterium inside host cells. As such, although H37Ra is observed to transcribe the genes encoding the anchoring domain at even a higher level than H37Rv, the lag in transcription of the catalytic subdomain could limit the function of the complex as a whole.
None of the frd genes were found in the TraSH (transposon site hybridization) screens used to identify genes essential for survival of the bacterium in macrophages and in mice (Rengarajan et al., 2005
; Sassetti & Rubin, 2003
). However, as the authors state when comparing their results to expression analysis, their findings emphasize the complementary nature of the two experimental approaches and that TraSH may fail to pick up mutants that are complemented in trans by either bacterial or host factors (Rengarajan et al., 2005
). In addition, TraSH would not identify differences such as we found, where differences in the kinetics of frd gene expression correlated with differences between H37Ra and H37Rv in their intracellular growth and virulence. Thus, the absence of the frd genes in the TraSH screens does not minimize the importance of FRD in M. tuberculosis infections. Studies of directed knockout mutants of all frd genes in addition to knockout mutants of individual frd genes have been initiated and are ongoing.
The transcriptomic data from the infection study suggested a role for FRD in the intracellular survival of M. tuberculosis. We therefore investigated the effect of a putative FRD inhibitor, MPNO, on the growth of intracellular M. tuberculosis. The FRD enzyme complex is not found in mammalian cells (Ge, 2002
; Phillips et al., 1987
; Van Hellemond & Tielens, 1994
), so FRD inhibitors are expected to have low to negligible adverse effects upon the viability of mammalian cells (Ge, 2002
). With regard to other pathogenic organisms, FRD has been considered a target in the treatment of Helicobacter pylori, as it was found to be an essential gene in the establishment of H. pylori colonization of the mouse stomach (Birkholz et al., 1994
; Ge et al., 2000
; Wang et al., 2000
). Furthermore, FRD has been reported to be a successful target in the treatment of protozoan and helminth infections using a variety of compounds (Bryant & Bennet, 1983
; Chen et al., 2001
; Omura et al., 2001
; Prichard, 1973
; Turrens et al., 1996
, 1999
), although it should be emphasized that the specificity of these compounds for FRD has not been established and the mode of action of at least some of the compounds is unclear.
MPNO is a common antimicrobial and antifungal ingredient in household products. It has been characterized to be an inhibitor of protozoa and is the precursor to L921-021, a drug developed and implemented for the treatment of protozoan infections. However, it should be stressed that although MPNO has been reported to inhibit FRD activity (Turrens et al., 1999
), its mode of action is not understood and the specificity of MPNO for FRD has not been demonstrated.
For our purposes of examining MPNO effects on the intracellular growth of M. tuberculosis H37Ra and H37Rv, we used the concentration of the inhibitor previously used to inhibit the growth of the intracellular parasite Trypanasoma cruzi without adverse effects on mammalian cells (Turrens et al., 1999
). The significant reduction in intracellular bacteria in MPNO-treated macrophages throughout the study was possibly due to inhibitory effects of the inhibitor on bacterial FRD. However, as the literature provides no real understanding of the mode of action of MPNO, it is not yet possible to determine whether the MPNO was acting on other bacterial systems to produce its mycobactericidal effect. Indeed, the specificity of MPNO and other reported FRD inhibitors has not been demonstrated (Bryant & Bennet, 1983
; Chen et al., 2001
; Omura et al., 2001
; Prichard, 1973
; Turrens et al., 1996
, 1999
). Further work will be necessary before we can be sure that other putative FRD inhibitors have activity against extracellular and intracellular M. tuberculosis and that any observed activity is indeed due solely to inhibition of FRD. Nevertheless, our observations indicate that MPNO and other putative FRD inhibitors warrant further investigation as anti-mycobacterials.
In this study we have outlined an alternative means to study transcriptomes of intracellular bacteria. Using BACFA, we have identified differences between virulent and attenuated strains of M. tuberculosis that may ultimately help to explain why these highly related bacteria have such different phenotypes when interacting with the host. We have also characterized the expression profile of the genes encoding the FRD enzyme complex, which in other organisms provides energy and substrates under suboptimal conditions, and which appears to play a supportive role in the survival of M. tuberculosis in macrophages.
| ACKNOWLEDGEMENTS |
|---|
Edited by: S. V. Gordon
| REFERENCES |
|---|
|
|
|---|
Betts, J. C., Lukey, P. T., Robb, L. C., McAdam, R. A. & Duncan, K. (2002). Evaluation of a nutrient starvation model of Mycobacterium tuberculosis persistence by gene and protein expression profiling. Mol Microbiol 43, 717–731.[CrossRef][Medline]
Bhargava, S., Tyagi, A. K. & Tyagi, J. S. (1990). tRNA genes in mycobacteria: organization and molecular cloning. J Bacteriol 172, 2930–2934.
Birkholz, S., Knipp, U., Lemma, E., Kroger, A. & Opferkuch, W. (1994). Fumarate reductase of Helicobacter pylori – an immunogenic protein. J Med Microbiol 41, 56–62.
Brosch, R., Gordon, S. V., Billault, A., Garnier, T., Eiglmeier, K., Soravito, C., Barrell, B. G. & Cole, S. T. (1998). Use of a Mycobacterium tuberculosis H37Rv bacterial artificial chromosome library for genome mapping, sequencing, and comparative genomics. Infect Immun 66, 2221–2229.
Brosch, R., Philipp, W. J., Stavropoulos, E., Colston, M. J., Cole, S. T. & Gordon, S. V. (1999). Genomic analysis reveals variation between Mycobacterium tuberculosis H37Rv and the attenuated M. tuberculosis H37Ra strain. Infect Immun 67, 5768–5774.
Bryant, C. & Bennet, E. M. (1983). Observations on the fumarate reductase system in Haemonchus contortus and their relevance to anthelmintic resistance and to strain variations of energy metabolism. Mol Biochem Parasitol 7, 281–292.[CrossRef][Medline]
Chen, M., Zhai, L., Christensen, S. B., Theander, T. G. & Kharazmi, A. (2001). Inhibition of fumarate reductase in Leishmania major and L. donovani by chalcones. Antimicrob Agents Chemother 45, 2023–2029.
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]
Collins, D. M. & De Lisle, G. W. (1984). DNA restriction endonuclease analysis of Mycobacterium tuberculosis and Mycobacterium bovis BCG. J Gen Microbiol 130, 1019–1021.
Collins, F. M. & Smith, M. M. (1969). A comparative study of the virulence of Mycobacterium tuberculosis measured in mice and guinea pigs. Am Rev Respir Dis 100, 631–639.[Medline]
Dasgupta, N., Kapur, V., Singh, K. K., Das, T. K., Sachdeva, S., Jyothisri, K. & Tyagi, J. S. (2000). Characterization of a two-component system, devR-devS, of Mycobacterium tuberculosis. Tuber Lung Dis 80, 141–159.[CrossRef][Medline]
Dolin, P. J., Raviglione, M. C. & Kochi, A. (1994). Global tuberculosis incidence and mortality during 1990–2000. Bull World Health Organ 72, 213–220.[Medline]
Gao, Q., Kripke, K., Arinc, Z., Voskuil, M. & Small, P. (2004). Comparative expression studies of a complex phenotype: cord formation in Mycobacterium tuberculosis. Tuberculosis (Edinb) 84, 188–196.[CrossRef][Medline]
Ge, Z. (2002). Potential of fumarate reductase as a novel therapeutic target in Helicobacter pylori infection. Expert Opin Ther Targets 6, 135–146.[CrossRef][Medline]
Ge, Z., Feng, Y., Dangler, C. A., Xu, S., Taylor, N. S. & Fox, J. G. (2000). Fumarate reductase is essential for Helicobacter pylori colonization of the mouse stomach. Microb Pathog 29, 279–287.[CrossRef][Medline]
Gillespie, J., Barton, L. L. & Rypka, E. W. (1988). Influence of oxygen tension on the respiratory activity of Mycobacterium phlei. J Gen Microbiol 134, 247–252.
Goldberg, I., Lonberg-Holm, K., Bagley, E. A. & Stieglitz, B. (1983). Improved conversion of fumarate to succinate by Escherichia coli strains amplified for fumarate reductase. Appl Environ Microbiol 45, 1838–1847.
Gordon, S. V., Brosch, R., Billault, A., Garnier, T., Eiglmeier, K. & Cole, S. T. (1999a). Identification of variable regions in the genomes of tubercle bacilli using bacterial artificial chromosome arrays. Mol Microbiol 32, 643–655.[CrossRef][Medline]
Gordon, S. V., Heym, B., Parkhill, J., Barrell, B. & Cole, S. T. (1999b). New insertion sequences and a novel repeated sequence in the genome of Mycobacterium tuberculosis H37Rv. Microbiology 145, 881–892.
Graham, J. E. & Clark-Curtiss, J. E. (1999). Identification of Mycobacterium tuberculosis RNAs synthesized in response to phagocytosis by human macrophages by selective capture of transcribed sequences (SCOTS). Proc Natl Acad Sci U S A 96, 11554–11559.
Guest, J. R. (1981). Partial replacement of succinate dehydrogenase function by phage- and plasmid-specified fumarate reductase in Escherichia coli. J Gen Microbiol 122, 171–179.
Guinn, K. M., Hickey, M. J., Mathur, S. K., Zakel, K. L., Grotzke, J. E., Lewinsohn, D. M., Smith, S. & Sherman, D. R. (2004). Individual RD1-region genes are required for export of ESAT-6/CFP-10 and for virulence of Mycobacterium tuberculosis. Mol Microbiol 51, 359–370.[CrossRef][Medline]
Heplar, J. Q., Clifton, C. E., Raffel, S. & Futrelle, C. M. (1954). Virulence of the tubercle bacillus. I. Effect of oxygen tension upon respiration of virulent and avirulent bacilli. J Infect Dis 94, 90–98.[Medline]
Imaeda, T. (1985). Deoxyribonucleic acid relatedness among strains of Mycobacterium tuberculosis, Mycobacterium bovis BCG, Mycobacterium microti and Mycobacterium africanum. Int J Syst Bacteriol 35, 147–150.
Jannasch, H. W. (1969). Estimations of bacterial growth rates in natural waters. J Bacteriol 99, 156–160.
Jung, Y. J., LaCourse, R., Ryan, L. & North, R. J. (2002). Virulent but not avirulent Mycobacterium tuberculosis can evade the growth inhibitory action of a T helper 1-dependent, nitric oxide synthase 2-independent defense in mice. J Exp Med 196, 991–998.
Kinger, A. K. & Tyagi, J. S. (1993). Identification and cloning of genes differentially expressed in the virulent strain of Mycobacterium tuberculosis. Gene 131, 113–117.[CrossRef][Medline]
Lari, N., Rindi, L. & Garzelli, C. (2001). Identification of one insertion site of IS6110 in Mycobacterium tuberculosis H37Ra and analysis of the RvD2 deletion in M. tuberculosis clinical isolates. J Med Microbiol 50, 805–811.
Li, M. S., Monahan, I. M., Waddell, S. J., Mangan, J. A., Martin, S. L., Everett, M. J. & Butcher, P. D. (2001). cDNA-RNA subtractive hybridization reveals increased expression of mycocerosic acid synthase in intracellular Mycobacterium bovis BCG. Microbiology 147, 2293–2305.
Madigan, M. T., Martinko, J. M. & Parker, J. (1997a). Nutrition and metabolism. In Brock Biology of Microorganisms, 8th edn, pp. 111–150. Upper Saddle River, NJ: Prentice-Hall.
Madigan, M. T., Martinko, J. M. & Parker, J. (1997b). Microbial growth. In Brock Biology of Microorganisms, 8th edn, pp. 151–179. Upper Saddle River, NJ: Prentice-Hall.
Mangan, J. A. & Butcher, P. D. (1998). Analysis of mycobacterial differential gene expression by RAP-PCR. Methods Mol Biol 101, 307–322.[Medline]
Mangan, J. A., Sole, K. M., Mitchison, D. A. & Butcher, P. D. (1997). An effective method of RNA extraction from bacteria refractory to disruption, including mycobacteria. Nucleic Acids Res 25, 675–676.
Menendez Mdel, C., Rebollo, M. J., Nunez Mdel, C., Cox, R. A. & Garcia, M. J. (2005). Analysis of the precursor rRNA fractions of rapidly growing mycobacteria: quantification by methods that include the use of a promoter (rrnA P1) as a novel standard. J Bacteriol 187, 534–543.
Monahan, I. M., Betts, J., Banerjee, D. K. & Butcher, P. D. (2001). Differential expression of mycobacterial proteins following phagocytosis by macrophages. Microbiology 147, 459–471.
Monod, J. (1949). The growth of bacterial cultures. Annu Rev Microbiol 3, 371–394.[Medline]
Mostowy, S., Cleto, C., Sherman, D. R. & Behr, M. A. (2004). The Mycobacterium tuberculosis complex transcriptome of attenuation. Tuberculosis (Edinb) 84, 197–204.[CrossRef][Medline]
Narain, J. P., Raviglione, M. C. & Kochi, A. (1992). HIV-associated tuberculosis in developing countries: epidemiology and strategies for prevention. Tuber Lung Dis 73, 311–321.[CrossRef][Medline]
Oatway, W. H. & Steenken, W (1936). The pathogenesis and fate of tubercle produced by dissociate variants of tubercle bacilli. J Infect Dis 59, 306–325.
Omura, S., Miyadera, H., Ui, H., Shiomi, K., Yamaguchi, Y., Masama, R., Nagamitsu, T., Takano, D., Sunazuka, T. & other authors (2001). An anthelmintic compound, nafuredin, shows selective inhibition of complex I in helminth mitochondria. Proc Natl Acad Sci U S A 98, 60–62.
Park, H. D., Guinn, K. M., Harrell, M. I., Liao, R., Voskuil, M. I., Tompa, M., Schoolnik, G. K. & Sherman, D. R. (2003). Rv3133c/dosR is a transcription factor that mediates the hypoxic response of Mycobacterium tuberculosis. Mol Microbiol 48, 833–843.[CrossRef][Medline]
Phillips, M. K., Hederstedt, L., Hasnain, S., Rutberg, L. & Guest, J. R. (1987). Nucleotide sequence encoding the flavoprotein and iron-sulfur protein subunits of the Bacillus subtilis PY79 succinate dehydrogenase complex. J Bacteriol 169, 864–873.
Prichard, R. K. (1973). The fumarate reductase reaction of Haemonchus contortus and the mode of action of some anthelmintics. Int J Parasitol 3, 409–417.[CrossRef][Medline]
Rahn, O. (1930). The non-logarithmic order of death of some bacteria. J Gen Physiol 13, 395–407.
Ramakrishnan, T., Indira, M. & Maller, R. K. (1962). Evaluation of the routes of glucose utilization in virulent and avirulent strains of Mycobacterium tuberculosis. Biochim Biophys Acta 59, 529–532.[Medline]
Raman, S., Song, T., Puyang, X., Bardarov, S., Jacobs, W. R., Jr & Husson, R. N. (2001). The alternative sigma factor SigH regulates major components of oxidative and heat stress responses in Mycobacterium tuberculosis. J Bacteriol 183, 6119–6125.
Rengarajan, J., Bloom, B. R. & Rubin, E. J. (2005). Genome-wide requirements for Mycobacterium tuberculosis adaptation and survival in macrophages. Proc Natl Acad Sci U S A 102, 8327–8332.
Rindi, L., Lari, N. & Garzelli, C. (1999). Search for genes potentially involved in Mycobacterium tuberculosis virulence by mRNA differential display. Biochem Biophys Res Commun 258, 94–101.[CrossRef][Medline]
Rindi, L., Lari, N. & Garzelli, C. (2001). Genes of Mycobacterium tuberculosis H37Rv downregulated in the attenuated strain H37Ra are restricted to M. tuberculosis complex species. New Microbiol 24, 289–294.[Medline]
Rivera-Marrero, C. A., Burroughs, M. A., Masse, R. A., Vannberg, F. O., Leimbach, D. L., Roman, J. & Murtagh, J. J., Jr (1998). Identification of genes differentially expressed in Mycobacterium tuberculosis by differential display PCR. Microb Pathog 25, 307–316.[CrossRef][Medline]
Rooyakkers, A. W. & Stokes, R. W. (2005). Absence of complement receptor 3 results in reduced binding and ingestion of Mycobacterium tuberculosis but has no significant effect on the induction of reactive oxygen and nitrogen intermediates or on the survival of the bacteria in resident and interferon-gamma activated macrophages. Microb Pathog 39, 57–67.[CrossRef][Medline]
Saini, D. K., Malhotra, V., Dey, D., Pant, N., Das, T. K. & Tyagi, J. S. (2004). DevR-DevS is a bona fide two-component system of Mycobacterium tuberculosis that is hypoxia-responsive in the absence of the DNA-binding domain of DevR. Microbiology 150, 865–875.
Sassetti, C. M. & Rubin, E. J. (2003). Genetic requirements for mycobacterial survival during infection. Proc Natl Acad Sci U S A 100, 12989–12994.
Schmittgen, T. D., Zakrajsek, B. A., Mills, A. G., Gorn, V., Singer, M. J. & Reed, M. W. (2000). Quantitative reverse transcription-polymerase chain reaction to study mRNA decay: comparison of endpoint and real-time methods. Anal Biochem 285, 194–204.[CrossRef][Medline]
Schnappinger, D., Ehrt, S., Voskuil, M. I., Liu, Y., Mangan, J. A., Monahan, I. M., Dolganov, G., Efron, B., Butcher, P. D. & other authors (2003). Transcriptional adaptation of Mycobacterium tuberculosis within macrophages: insights into the phagosomal environment. J Exp Med 198, 693–704.
Selwyn, P. A., Hartel, D., Lewis, V. A., Schoenbaum, E. E., Vermund, S. H., Klein, R. S., Walker, A. T. & Friedland, G. H. (1989). A prospective study of the risk of tuberculosis among intravenous drug users with human immunodeficiency virus infection. N Engl J Med 320, 545–550.[Abstract]
Sirakova, T. D., Thirumala, A. K., Dubey, V. S., Sprecher, H. & Kolattukudy, P. E. (2001). The Mycobacterium tuberculosis pks2 gene encodes the synthase for the hepta- and octamethyl-branched fatty acids required for sulfolipid synthesis. J Biol Chem 276, 16833–16839.
Steenken, W. & Gardner, L. U. (1946). History of H37 strain of tubercle bacillus. Am Rev Respir Dis 54, 62–66.
Steenken, W., Jr, Oatway, W. H. & Petroff, S. A. (1934). Biological studies of the tubercle bacillus. III. Dissociation and pathogenicity of the R and S variants of the human tubercle bacillus (H37). J Exp Med 60, 515–540.[Abstract]
Stewart, G. R., Wernisch, L., Stabler, R., Mangan, J. A., Hinds, J., Laing, K. G., Young, D. B. & Butcher, P. D. (2002). Dissection of the heat-shock response in Mycobacterium tuberculosis using mutants and microarrays. Microbiology 148, 3129–3138.
Stokes, R. W., Haidl, I. D., Jefferies, W. A. & Speert, D. P. (1993). Mycobacteria–macrophage interactions. Macrophage phenotype determines the nonopsonic binding of Mycobacterium tuberculosis to murine macrophages. J Immunol 151, 7067–7076.[Abstract]
Talaat, A. M., Hunter, P. & Johnston, S. A. (2000). Genome-directed primers for selective labeling of bacterial transcripts for DNA microarray analysis. Nat Biotechnol 18, 679–682.[CrossRef][Medline]
Turrens, J. F., Watts, B. P., Jr, Zhong, L. & Docampo, R. (1996). Inhibition of Trypanosoma cruzi and T. brucei NADH fumarate reductase by benznidazole and anthelmintic imidazole derivatives. Mol Biochem Parasitol 82, 125–129.[CrossRef][Medline]
Turrens, J. F., Newton, C. L., Zhong, L., Hernandez, F. R., Whitfield, J. & Docampo, R. (1999). Mercaptopyridine-N-oxide, an NADH-fumarate reductase inhibitor, blocks Trypanosoma cruzi growth in culture and in infected myoblasts. FEMS Microbiol Lett 175, 217–221.[CrossRef][Medline]
Van Hellemond, J. J. & Tielens, A. G. (1994). Expression and functional properties of fumarate reductase. Biochem J 304, 321–331.[Medline]
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A. & Speleman, F. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3, RESEARCH0034[Medline]
Waddell, S. J., Chung, G. A., Gibson, K. J., Everett, M. J., Minnikin, D. E., Besra, G. S. & Butcher, P. D. (2005). Inactivation of polyketide synthase and related genes results in the loss of complex lipids in Mycobacterium tuberculosis H37Rv. Lett Appl Microbiol 40, 201–206.[CrossRef][Medline]
Wang, S., Liu, F. & Zhao, B. (2000). Study on mRNA differential display in Mycobacterium tuberculosis H37Rv and H37Ra. Zhonghua Jie He He Hu Xi Za Zhi 23, 669–671.[Medline]
Wang, S. Y., Zheng, H. J., Wang, B. F., Zhang, X. L., Pu, S. Y., Zhu, G. F. & Zhao, G. P. (2007). Complete genomic sequence of Mycobacterium tuberculosis strain H37Ra, a non-pathogenic variant closely related to the well-characterized pathogenic strain H37Rv. GenBank Accession Number: CP000611, http://www.ncbi.nlm.nih.gov/entrez/viewer.fcgi?db=nuccore&id=148503909.
Wayne, L. G. (1976). Dynamics of submerged growth of Mycobacterium tuberculosis under aerobic and microaerophilic conditions. Am Rev Respir Dis 114, 807–811.[Medline]
Wei, J., Dahl, J. L., Moulder, J. W., Roberts, E. A., O'Gaora, P., Young, D. B. & Friedman, R. L. (2000). Identification of a Mycobacterium tuberculosis gene that enhances mycobacterial survival in macrophages. J Bacteriol 182, 377–384.
Wernisch, L., Kendall, S. L., Soneji, S., Wietzorrek, A., Parish, T., Hinds, J., Butcher, P. D. & Stoker, N. G. (2003). Analysis of whole-genome microarray replicates using mixed models. Bioinformatics 19, 53–61.
Received 4 January 2008;
revised 16 April 2008;
accepted 16 April 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 | |