Iron-sparing Response of Mycobacterium avium subsp. paratuberculosis is strain dependent
© Janagama et al; licensee BioMed Central Ltd. 2010
Received: 23 May 2010
Accepted: 22 October 2010
Published: 22 October 2010
Two genotypically and microbiologically distinct strains of Mycobacterium avium subsp. paratuberculosis (MAP) exist - S and C MAP strains that primarily infect sheep and cattle, respectively. Concentration of iron in the cultivation medium has been suggested as one contributing factor for the observed microbiologic differences. We recently demonstrated that S strains have defective iron storage systems, leading us to propose that these strains might experience iron toxicity when excess iron is provided in the medium. To test this hypothesis, we carried out transcriptional and proteomic profiling of these MAP strains under iron-replete or -deplete conditions.
We first complemented M. smegmatis ΔideR with IdeR of C MAP or that derived from S MAP and compared their transcription profiles using M. smegmatis mc 2 155 microarrays. In the presence of iron, sIdeR repressed expression of bfrA and MAP2073c, a ferritin domain containing protein suggesting that transcriptional control of iron storage may be defective in S strain. We next performed transcriptional and proteomic profiling of the two strain types of MAP under iron-deplete and -replete conditions. Under iron-replete conditions, C strain upregulated iron storage (BfrA), virulence associated (Esx-5 and antigen85 complex), and ribosomal proteins. In striking contrast, S strain downregulated these proteins under iron-replete conditions. iTRAQ (isobaric tag for relative and absolute quantitation) based protein quantitation resulted in the identification of four unannotated proteins. Two of these were upregulated by a C MAP strain in response to iron supplementation. The iron-sparing response to iron limitation was unique to the C strain as evidenced by repression of non-essential iron utilization enzymes (aconitase and succinate dehydrogenase) and upregulation of proteins of essential function (iron transport, [Fe-S] cluster biogenesis and cell division).
Taken together, our study revealed that C and S strains of MAP utilize divergent metabolic pathways to accommodate in vitro iron stress. The knowledge of the metabolic pathways these divergent responses play a role in are important to 1) advance our ability to culture the two different strains of MAP efficiently, 2) aid in diagnosis and control of Johne's disease, and 3) advance our understanding of MAP virulence.
Mycobacterium avium subsp. paratuberculosis (MAP), the causative agent of Johne's disease (JD) of ruminants, often requires eight to sixteen weeks to see colonies in culture - a major hurdle in the diagnosis and therefore in implementation of optimal control measures. Unlike other mycobacteria, which mobilize iron via mycobactins, MAP is unable to produce detectable mycobactin in vitro or in vivo [1–3]. Although the reasons for the in vitro mycobactin dependency of MAP are currently unknown, we have recently shown that the mycobactin (mbt) operon promoter is active and that the mycobactin genes are transcribed by MAP inside macrophages  and in tissues of naturally infected animals (accepted for publication in BMC Genomics).
Pathogenic mycobacteria encounter a wide variety of stressors inside the host cells and their ability to overcome iron deprivation and iron toxicity represents a major virulence determinant . Transcript and protein profiling of MTB and other pathogens in response to in vitro iron stress is well documented [6–9]. While MAP transcriptome or proteome profiles in response to heat shock, pH, oxidative stress, hypoxia, and nutrient starvation have been demonstrated [10–12], stress responses to iron supplementation or starvation are lacking.
Iron dependent regulator (IdeR) has been very well studied as a global regulator involved in maintaining iron homeostasis in Mycobacterium tuberculosis (MTB) . Recently we have demonstrated that IdeR of MAP in the presence of iron recognizes a consensus sequence on the promoter called "iron box" and regulates expression of genes involved in iron acquisition (mbt) and storage (bfrA). More interestingly, we demonstrated that polymorphisms in the promoter of iron storage gene (bfrA) in S MAP strains relative to C MAP strains results in a differential gene regulation . IdeR dependent repression of bfrA in the presence of iron suggests variations in iron storage mechanisms and/or iron requirements in cattle and sheep MAP strains.
Comparative genomic hybridizations, short sequence repeat analysis and single nucleotide polymorphisms of MAP isolates obtained from diverse host species have established and indexed genomic differences between C and S strains of MAP [14–19]. Phylogenetic analysis of sequences has identified C and S strains as separate pathogenic clones that share a common ancestor [20–23]. Furthermore, cellular infection studies show distinctive phenotypes between the two MAP strain types [24, 25]. We also recently demonstrated that S strains have defective iron storage systems, leading us to propose that these strains might experience iron toxicity when excess iron is provided in the medium . Taken together, the literature suggests that MAP strains vary in their iron dependent gene regulation. To test this further, we profiled their transcriptomes and proteomes in response to iron and demonstrated that iron induced metabolic pathways are significantly diverse.
Bacterial strains, DNA manipulations and media
Mycobacterium avium subsp. paratuberculosis strains MAP1018 (C MAP) and MAP7565 (S MAP) were grown in Middlebrook 7H9 supplemented with OADC enrichment medium and mycobactin J (2 mg/mL; Allied Monitor, Fayette, MO).
To test the hypothesis that gene regulation may be dependent on iron availability MAP strains were grown in Middlebrook 7H9 medium without mycobactin J or Sauton medium (0.5 g KH2PO4, 0.5 g MgSO4, 4.0 g L-asparagine, 60 ml glycerol, 0.05 g ferric ammonium citrate, 2.0 g citric acid, 0.1 ml 1% (w/v) ZnSO4 and 2.5 ml 20% Tween 80 in 1 liter). Growth of MAP strains in the absence of mycobactin J took over 6 months to provide sufficient material for proteomics and transcriptional profiling. For iron restriction, 2,2'-dipyridyl (Sigma Aldrich, St. Louis, MO) was added at a concentration of 200 μM.
MAP7565 and MAP1018 have been genotyped by SSR as well as comparative genomics using oligoarrays. They represent the typical genomotypes of sheep and cattle strains, respectively  and show distinct phenotypes in both human and bovine macrophages [24, 25].
M. smegmatis (mc2155) and E. coli TOP10F (Invitrogen Corporation, Carlsbad, CA) competent cells were grown in Luria Bertani (LB) medium and antibiotics (kanamycin (20 μg/ml) or hygromycin (100 μg/ml)) were added when necessary. The open reading frames of ideR (MAP2827) derived from C or S MAP strains were cloned into pSM417 and M. smegmatis ΔideR (SM3) was complemented as previously reported . Briefly, MAP2827 from MAP1018 (cideR) or MAP 7565 (sideR) was amplified via PCR using primers that carried restriction sites for BamHI and HindIII. Amplified products were double digested with BamHI and HindIII and ligated into a pre digested (BamHI and HindIII) expression plasmid pSM417. Accuracy of the ligation and orientation of MAP2827 in pSM417 was verified by sequencing. SM3 was transformed with pSM417 carrying MAP2827 from C or S MAP strains.
A seed stock from logarithmically grown (OD600 = 1.0) cultures were diluted to fresh medium to yield an OD600 = 0.1. These were grown in various aliquots under constant shaking (120 rpm) at 37°C. These cultures were monitored for their growth at weekly intervals by measuring their absorbance at 600 nm wave length using SpectraMax M2 (Molecular Devices, Sunnyvale, CA) until they reached an absorbance of 1.0 (Additional file 1, Figure S1). At this point, the cultures were then pelleted, washed in ice cold 1XPBS and re-suspended in fresh culture medium (with or without the addition of 2,2'-dipyridyl (Sigma Aldrich, St. Louis, MO)). Dipyridyl was added at a concentration of 200 μM. Following three hours of incubation at 37°C under constant shaking, cells were pelleted and washed with ice cold 1X PBS and either used in microarrays or iTRAQ. The detailed experimental design is provided as Additional file 1, Figure S2.
Nucleic acid and protein extraction
Log phase MAP or M. smegmatis cultures were pelleted, washed and re-suspended in fresh culture medium with or without 200 μM of 2,2'-dipyridyl. The cultures were incubated at 37°C with shaking for 3 hr. immediately prior to RNA and protein extraction.
For RNA, cells were homogenized in Mini bead-beater for 4 min. by adding 0.3 ml of 0.1 mm sterile RNase-free zirconium beads followed by extraction using Trizol (Invitrogen, Carlsbad, CA). All samples were treated with RNase-free DNase I (Ambion, Inc., Austin, TX) to eliminate genomic DNA contamination. The purity and yield of total RNA samples was confirmed using Agilent 2100E Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA). RNA was stored at -80 until used in microarrays and real time RT-PCR assays.
For protein, cells were re-suspended in minimal quantity (250 μL) of iTRAQ dissolution buffer (0.5 M TEAB pH 8.5) and 0.1% SDS. The solution was transferred to a 2 ml screw cap tube containing 0.1 mm zirconium beads (Biospec) and disrupted in minibead beater (Biospec) for 4 × 1 minute pulses with samples kept on ice between pulses. The lysate was then centrifuged at 12,000 × g for 10 minutes at 4°C. Supernatant was transferred to a fresh tube without disturbing the pellet and used in iTRAQ labeling for detection of proteome (Additional file 1, Figure S3).
Gene expression profiling of S (1018) and C (7565) MAP strains was performed using MAP K-10 microarrays obtained from Dr. Michael Paustian, NADC, IA. Expression profiling of M. smegmatis ΔideR complemented with c or sideR was carried out using M. smegmatis mc 2 155 arrays provided via Pathogen Functional Genomics Resource Center (PFGRC) at J. Craig Venter Institute (JCVI). Array hybridizations and analyses were performed as described previously and according to the protocols established at PFGRC with minor modifications  and according to MIAME 2.0 guidelines.
Briefly, synthesis of fluorescently labeled cDNA (Cyanine-3 or Cyanine-5) from total RNA and hybridizations of labeled cDNA to MAP K-10 or mc 2 155 oligoarray was performed. Microarray hybridizations were performed from cDNA isolated from two independent experiments. On each independent occasion, bacterial cultures growing under iron-replete or iron-limiting medium were used for RNA extractions, cDNA labeling and array hybridizations. Each slide was competitively hybridized with cDNA obtained from iron-replete (labeled with cy3 or cy5) and iron-limiting growth medium (counter labeled with cy5 or cy3) to reveal relative expressional differences. About 4 μg (2 μg each from iron limitation or sufficient) of cDNA was used to hybridize onto the array. However, if the cDNA yield is low for a sample the RNA from the same sample was used to synthesize more cDNA, pooled and labeled onto the arrays. Hybridized slides were scanned using HP Scan array 5000 (PerkinElmer Inc., Waltham, MA). The images were processed and numerical data was extracted using the microarray image analysis software, BlueFuse (BlueGnome Ltd, Cambridge) and TM4 microarray suite available through JCVI. Genes differentially regulated at a fold change of 1.5 or greater were identified at a false discovery rate of 1% by Statistical Analysis of Microarrays (SAM) program . Genes that showed a fold change 1.5 or greater in all the replicate arrays were retained and reported as being up- or downregulated in the presence of iron.
RNA isolated from MAP strains grown under iron-replete or iron-limiting growth medium was used in real time RT-PCR assays. Genes were selected based on their diverse roles and microarray expression pattern. Selected genes included siderophore transport (MAP2413c, MAP2414c), esx-3 secretion system (MAP3783, MAP3784), aconitase (MAP1201c), fatty acid metabolism (MAP0150c) and virulence (MAP0216, MAP3531c, MAP1122 and MAP0475). RNA was treated with DNaseI (Ambion, Austin, TX) and one step Q-RT PCR was performed using QuantiFast SYBR Green mix (Qiagen, Valencia, CA) and gene specific primers (Additional file 1, Table S1) in a Lightcycler 480 (Roche, Indianapolis, IN).
Protein extracted from the two MAP strains grown in iron-replete or iron-limiting medium was used in iTRAQ analysis (Additional file 1, Figure S3). iTRAQ labeling and protein identification was carried out as described previously with minor modifications .
Briefly, cell lysate was quantified using the bicinchoninic acid (BCA) protein assay (Pierce, Rockford, IL) prior to trypsin digestion. Peptides were labeled with iTRAQ reagents (114 and 115 for MAP 1018 grown in iron-replete and iron-limiting medium respectively; 116 and 117 for MAP 7565 grown in iron-replete and iron-limiting medium respectively) at lysine and arginine amino terminal groups. The labeled peptides were pooled, dried and re-suspended in 0.2% formic acid. The re-suspended peptides were passed through Oasis® MCX 3CC (60 mg) extraction cartridges per manufacturer recommendations (Waters Corporation, Milford, MA) for desalting prior to strong cation exchange (SCX) fractionation.
Eluted peptides were dried and dissolved in SCX buffer A (20% v/v ACN and 5 mM KH2PO4 pH 3.2, with phosphoric acid) and fractionated using a polysulfoethyl A column (150 mm length × 1.0 mm ID, 5 μm particles, 300 Å pore size) (PolyLC Inc., Columbia, MD) on a magic 2002 HPLC system (Michrom BioResources, Inc., Auburn, CA). Peptides were eluted by running a 0-20% buffer B gradient for greater than 55 min. and 20%-100% buffer B (20% v/v ACN, 5 mM KH2PO4 pH 3.2, 500 mM KCL) for 20 min. at a column flow rate of 50 μl/min. Several fractions were collected at frequent intervals and seven fractions that showed mAU280 > 2 were analyzed by LC-MS/MS as previously described .
Fractions were reconstituted in reversed-phase load buffer (10 mM phosphate buffer) and analyzed in a 4800 MALDI TOF/TOF instrument (AB Sciex, Foster City, CA). Protein pilot Software™ 3.0.1 (AB Sciex, Foster city, CA) which utilizes the paragon™ scoring algorithm  was used to identify and quantify the relative abundance of the labeled peptides. Relative abundance of proteins (iron-replete v/s iron-limitation) for each MAP strain was determined by comparing the reporter ion ratios (114/115 for C and 116/117 for S MAP). iTRAQ experiments were repeated on two independent experiments for each treatment of each strain. We searched against the MAP K-10, non redundant (nr) mycobacteria proteins and entire nr protein database deposited in the NCBI along with the contaminants to identify MAP specific peptides at a false discovery rate of 1%.
Transcriptional profiling of MAP IdeR
We recently characterized MAP IdeR and computationally predicted that IdeR in the presence of iron regulates expression of 24 genes . In the current study, we identified that 20 of the 24 previously predicted genes were differentially expressed in response to iron by MAP microarrays. Mycobactin synthesis, transport and fatty acid biosynthesis genes were repressed in the presence of iron by both cattle and sheep MAP strains (Additional file 1, Table S2). However iron storage and oxidoreductase genes were upregulated in the presence of iron only in C MAP (Additional file 1, Figure S4).
We first confirmed if these differences are due to regulation via IdeR. IdeR is essential in MAP and attempts to delete this gene failed . We complemented M. smegmatis ΔideR (SM3) with C or S strain ideR and compared regulational differences in the presence or absence of iron. Genes that showed a log2 fold change of 1.0 in SM3 or SM3 complemented with empty plasmid (negative controls) in the presence or absence of iron while having a fold change >± 1.5 in the complemented strains (test) and plasmid carrying M. smegmatis ideR and mc 2 155 (wild type) (positive control) were considered as being regulated by MAP IdeR. Fourteen of the 20 genes were regulated by IdeRs of both MAP strains in M. smegmatis. Furthermore, our results suggested that sIdeR functions by primarily repressing genes in the presence of iron whereas cIdeR functions both by repressing mycobactin synthesis and de-repressing iron storage genes in the presence of iron (Additional file 1, Table S3). These were further validated by realtime RT-PCR in both M. smegmatis transformants carrying MAP ideRs and MAP genetic background. The data is presented only for MAP (Additional file 1, Table S4).
Transcript profiles under iron-limiting conditions
Under iron-limiting conditions both the MAP strains showed increased transcription of genes belonging to mycobactin synthesis and esx-3, an essential secretory system of mycobactin biosynthesis (Additional file 1, Tables S2 - S5) .
Transcript and protein expression in cattle MAP under iron-limiting (LI) conditions
MAP ORF ID
2.03 ± 0.2
2.87 ± 0.7
FadE33_2 (acyl-coA synthase)
1.79 ± 0.5
1.88 ± 0.8
pdhC alpha-keto acid dehydrogenase
1.68 ± 0.3
2.52 ± 0.4
FadE23 (acyl-CoA dehydrogenase)
2.41 ± 0.2
3.51 ± 1.0
FadE5 (acyl-CoA dehydrogenase)
1.87 ± 0.8
3.15 ± 0.2
heat shock protein
2.18 ± 0.6
2.48 ± 0.3
FeS assembly protein SufD
2.23 ± 1.0
2.73 ± 0.2
FeS assembly ATPase SufC
1.78 ± 0.5
2.03 ± 0.1
heat shock protein HtpX
1.48 ± 0.1
1.66 ± 0.5
Poorly characterized pathways
1.67 ± 0.3
1.56 ± 0.3
iron suphur cluster biosynthesis
1.56 ± 0.9
1.66 ± 0.2
1.84 ± 0.3
2.19 ± 0.8
RES domain containing protein
1.50 ± 0.7
2.40 ± 0.2
Transcript and protein expression in sheep MAP under iron-limiting (LI) conditions
MAP ORF ID
1.54 ± 0.1
1.58 ± 0.6
1.74 ± 0.3
2.05 ± 1.0
thioredoxin domain containing protein
1.82 ± 0.1
2.04 ± 0.3
acyl carrier protein
1.90 ± 0.5
1.68 ± 0.5
1.50 ± 0.4
2.29 ± 0.3
DnaK molecular chaperone
1.63 ± 0.6
3.52 ± 0.5
Information storage and processing
FusA, elongation factor G
1.52 ± 0.2
2.58 ± 0.7
transcriptional regulatory protein
1.52 ± 0.3
1.50 ± 0.1
DNA-directed RNA polymerase alpha subunit
1.56 ± 0.1
1.83 ± 0.3
DNA binding protein, HU
1.60 ± 0.6
1.81 ± 0.5
30S ribosomal protein S5
1.75 ± 0.1
1.55 ± 0.3
1.94 ± 0.3
1.59 ± 0.2
transcription antitermination protein, NusG
1.98 ± 0.3
1.82 ± 0.5
elongation factor Tu
2.08 ± 0.4
2.16 ± 0.1
Poorly characterized pathways
conserved alanine and arginine rich protein
1.54 ± 0.2
2.27 ± 0.5
initiation of DNA replication
1.63 ± 0.1
1.91 ± 0.2
transcriptional regulator like protein
1.75 ± 0.6
1.50 ± 0.2
1.83 ± 1.0
1.52 ± 0.5
Transcript profiles under iron-replete conditions
Transcript and protein expression in cattle MAP under iron-replete (HI) conditions
MAP ORF ID
FadE25_2 (acyl-coA dehydrogenase)
1.72 ± 0.1
1.88 ± 0.2
1.73 ± 0.3
1.56 ± 0.1
1.69 ± 0.2
3.68 ± 0.3
Fas (fatty acid synthase)
1.61 ± 0.5
2.28 ± 0.4
AccA3 (acetyl-/propionyl-coenzyme A)
1.45 ± 0.1
2.18 ± 0.2
1.89 ± 0.3
4.57 ± 0.5
iron regulated conserved protein
1.62 ± 0.2
0.78 ± 0.3
1.79 ± 0.5
2.29 ± 0.2
Information storage and processing
translation initiation factor IF-2
1.57 ± 0.2
1.89 ± 0.2
ribosome releasing factor
1.66 ± 0.3
2.11 ± 0.5
50S ribosomal protein L1
1.61 ± 0.1
1.57 ± 0.2
rplJ 50S ribosomal protein L10
1.52 ± 0.1
1.66 ± 0.5
fusA elongation factor G
2.13 ± 0.4
3.05 ± 0.3
rpsJ 30S ribosomal protein S10
1.68 ± 0.3
2.87 ± 0.4
rpsH 30S ribosomal protein S8
1.79 ± 0.5
2.42 ± 0.1
rpoA DNA-directed RNA polymerase
1.56 ± 0.1
1.65 ± 0.4
Poorly characterized pathways
FbpA antigen 85-A
1.87 ± 0.2
2.16 ± 0.3
mycobacterial integration host factor
1.73 ± 0.3
2.00 ± 0.5
In contrast, we did not document any upregulation (at a log2 fold change of 1.5) in the S MAP under iron-replete conditions. The directionality of transcripts as identified by microarrays under iron-replete conditions by S MAP strain was confirmed by real time RT-PCR (Additional file 1, Table S4).
The following criteria were used for protein identification in each treatment - (1) peptides identified by mass spectrometry were searched against the non-redundant (nr) protein database deposited in NCBI); and (2) MAP specific peptides reported with >95% confidence were used to quantify the relative abundance (iron-replete v/s iron-limitation) of each protein. A peptide with no hits on the MAP genome but with identities with other mycobacterial proteins was considered as unannotated MAP protein.
Protein expression under iron-limiting conditions
Consistent with the transcription profile, the C strain of MAP upregulated proteins belonging to SUF operon involved in Fe-S cluster assembly, fatty acid metabolism and a pyruvate dehydrogenase (MAP2307c). Transporter proteins, two component systems, and cell division associated proteins (MAP1906c, MAP0448 and MAP2997c) were also upregulated by the C strain (Table 1 and Additional file 1, Table S8). The sheep strain also upregulated transporter proteins, fatty acid biosynthesis, DNA replication protein (MAP3433), and stress response proteins (MAP3831c, MAP2764) (Table 2, Additional file 1, Table S9 and Figure S3).
Protein expression under iron-replete conditions
Identification of unannotated MAP proteins
We identified two unique peptides (SSHTPDSPGQQPPKPTPAGK and TPAPAKEPAIGFTR) that originated from the unannotated MAP gene located between MAP0270 (fadE36) and MAP0271 (ABC type transporter). We also identified two peptides (DAVELPFLHK and EYALRPPK) that did not map to any of the annotated MAP proteins but to the amino acid sequence of MAV_2400. Further examination of the MAP genome revealed that the peptides map to the reversed aminoacid sequence of MAP1839. These two unique proteins were not differentially regulated in response to iron. However, two more unique peptides that were translated from other unannotated MAP genes were upregulated (>1.5 fold) under iron-replete conditions in C MAP strain (Figure 3B).
Johne's disease is a major animal health problem of ruminant species worldwide and imposes significant economic losses to the industry. Our ability to culture the causative agent--Mycobacterium avium subsp. paratuberculosis (MAP)--and therefore its rapid diagnosis and our understanding of its virulence is limited. MAP is difficult to culture because of its unusually strict iron requirements. For optimal growth in laboratory media, MAP requires a siderophore (mycobactin) supplementation that makes MAP fastidious ., often requiring eight to sixteen weeks to produce colonies in culture - a major hurdle in the diagnosis and therefore implementation of optimal control measures. Understanding iron regulatory networks in the pathogen invitro is therefore of great importance.
A tale of two strain types of MAP - A case to study iron regulation
Several microbiological and genotyping studies and clinical observations suggest that Johne's in certain hosts such as sheep, goats, deer, and bison is caused by a distinct set of strains that show a relatively high degree of host preference [18, 40]. At least two microbiologically distinct types of MAP have been recognized. A less readily cultivable type is the common, but not invariable, cause of paratuberculosis in sheep (S MAP) [39, 41, 42], while another readily cultivable type is the most common cause of the disease in cattle (C MAP). Cell infection studies have also revealed distinctive host response phenotypes between cattle and sheep MAP strains - the former elicit primarily a pro-inflammatory response while latter strains suppress inflammation and upregulate anti-apoptotic pathways [24, 25]. In addition, since MAP genome sequence was published in 2005, very little research has focused on iron physiology and its contribution to metabolic networks of this fastidious organism.
Based on these classical microbiologic, genotypic, and clinical observations, we addressed the hypothesis that the iron dependent gene regulation is different between cattle and sheep MAP strains using a systems approach.
Iron-sparing response to iron-limitation is unique to cattle MAP strain
Iron is a critical component of several metabolic enzymes . Most bacteria respond to iron starvation with a unique regulatory mechanism called the iron-sparing response . Iron-sparing is a physiological phenomenon used by cells to increase the intracellular iron pool by post-transcriptionally repressing the synthesis of non-essential iron using proteins and sparing iron for essential cellular functions . Therefore, the paradigm is to transcriptionally upregulate all iron uptake systems while repressing non-essential enzymes via post-transcriptional regulatory mechanisms to survive iron-limiting conditions. Both MAP strains upregulated genes involved in siderophore biosynthesis (mbt), ability to acquire iron from synthesized siderophores (esx-3), and to transport iron bound siderophores (irtAB) into the bacterium (Additional file 1, Table S2). Furthermore, cattle MAP strain under iron-limiting conditions upregulated transcription of aconitase (Additional file 1, Table S4) while downregulating its protein expression (Figure 2). It is likely that targets for post-transcriptional repression of these non-essential iron using proteins are mediated via small RNAs . Studies to test this hypothesis in the two MAP strain types are underway.
Differential metabolic responses of cattle and sheep MAP strains to iron-limitation
Under iron-limiting conditions most other bacteria including M. tuberculosis (MTB) upregulate SUF operon [26, 45]. SUF synthesizes [Fe-S] clusters and transports them to iron-sulfur containing proteins involved in diverse cellular functions such as redox balance and gene regulation . This is critical because unlike E. coli, MTB and MAP genomes encode for only one such system to synthesize all the [Fe-S] needed by the cell and free iron and sulfide atoms are toxic to cells . Our data revealed that cattle strain, but not S strain upregulated SUF operon at the transcript and protein level under iron-limiting conditions (Table 1).
Cattle MAP strain upregulated pyruvate dehydrogenase operon involved in catabolism of propionate a key component of lipid biosynthesis under limiting iron conditions . In contrast, sheep strain upregulated isoprenoid synthesis genes involved in cell wall biogenesis . The sheep isolate also upregulated oxidoreductase and stress responses in its transcriptome and proteome during iron-limitation (Table 2). CarD and toxin-antitoxin systems primarily function during unfavorable conditions such as starvation or oxidative stress by arresting cell growth [50, 51]. Sheep strain upregulated transcripts of toxin-antitoxin system involved in arresting cell growth, suggesting a trend toward stringency response (Additional file 1, Table S6). Taken together, our data suggests that cattle strain is able to efficiently modulate its metabolism during iron-limitation - probably a survival advantage.
MAP2325, a hypothetical protein deleted in the sheep strain was found to be upregulated under iron-limiting conditions by the C strain (Additional file 1, Table S5). This is interesting because an ortholog of MAP2325 in MTB called enhanced intracellular survival (eis) interacts with host T cells. Stimulation of recombinant Eis from MTB results in increased production of IL-10 and decreased production of TNF-α thus contributing to mycobacterial survival inside macrophages . We have also demonstrated a similar result in bovine or human macrophages stimulated with diverse MAP strains. Cattle strains produced relatively more IL-10 and less TNF-α and persisted for longer periods of time inside macrophages [24, 25].
There is increased protein synthesis and turn over in response to iron in MTB . Similarly, we observed an increased expression of ribosomal proteins in the transcriptome and proteome in C MAP under iron-replete conditions. In striking contrast, iron-limitation induced a similar theme in sheep strain. Heme degradation is a significant physiological phenomenon where in pathogens recycle iron and gain a survival advantage inside the host . Recently the crystal structure of Rv3592 of MTB was solved and demonstrated its ability as heme degrader . We observed an upregulation of MAP0467c protein (ortholog of Rv3592) under iron-replete conditions in C MAP while it was downregulated in the sheep strain (Figure 4). Similar to our previous reports, iron storage protein, BfrA was upregulated only by C MAP under iron-repletion (Figure 3) . Although the reasons for differential iron storage mechanisms in sheep compared to cattle strains of MAP are currently unknown, differential role of ferritins in bacterial pathogens is not uncommon .
This work was supported in part by a USDA-NRI grant (2005-35204-16106) and Johne's disease Integrated Program (USDA-CSREES 2008-55620-18710) awarded to SS. We would like to thank Microbial and Plant Genomics Institute, Biomedical Genomics Center and Computational Genetics Laboratory at the University of Minnesota for providing resources and services to perform these studies. We would also like to thank JCVI for providing M. smegmatis microarrays.
- Lambrecht RS, Collins MT: Mycobacterium paratuberculosis. Factors that influence mycobactin dependence. Diagn Microbiol Infect Dis. 1992, 15 (3): 239-246. 10.1016/0732-8893(92)90119-E.View ArticlePubMedGoogle Scholar
- Lambrecht RS, Collins MT: Inability to detect mycobactin in mycobacteria-infected tissues suggests an alternative iron acquisition mechanism by mycobacteria in vivo. Microb Pathog. 1993, 14 (3): 229-238. 10.1006/mpat.1993.1022.View ArticlePubMedGoogle Scholar
- Snow GA: Mycobactins: iron-chelating growth factors from mycobacteria. Bacteriol Rev. 1970, 34 (2): 99-125.PubMed CentralPubMedGoogle Scholar
- Janagama HK, Senthilkumar TM, Bannantine JP, Rodriguez GM, Smith I, Paustian ML, McGarvey JA, Sreevatsan S: Identification and functional characterization of the iron-dependent regulator (IdeR) of Mycobacterium avium subsp. paratuberculosis. Microbiology. 2009, 155 (Pt 11): 3683-3690. 10.1099/mic.0.031948-0.PubMed CentralView ArticlePubMedGoogle Scholar
- Waddell SJ, Butcher PD: Microarray analysis of whole genome expression of intracellular Mycobacterium tuberculosis. Curr Mol Med. 2007, 7 (3): 287-296. 10.2174/156652407780598548.PubMed CentralView ArticlePubMedGoogle Scholar
- Rao PK, Li Q: Protein turnover in mycobacterial proteomics. Molecules. 2009, 14 (9): 3237-3258. 10.3390/molecules14093237.View ArticlePubMedGoogle Scholar
- Rao PK, Roxas BA, Li Q: Determination of global protein turnover in stressed mycobacterium cells using hybrid-linear ion trap-fourier transform mass spectrometry. Anal Chem. 2008, 80 (2): 396-406. 10.1021/ac701690d.View ArticlePubMedGoogle Scholar
- Rao PK, Li Q: Principal Component Analysis of Proteome Dynamics in Iron-starved Mycobacterium Tuberculosis. J Proteomics Bioinform. 2009, 2 (1): 19-31. 10.4172/jpb.1000058.PubMed CentralView ArticlePubMedGoogle Scholar
- Hindre T, Bruggemann H, Buchrieser C, Hechard Y: Transcriptional profiling of Legionella pneumophila biofilm cells and the influence of iron on biofilm formation. Microbiology. 2008, 154 (Pt 1): 30-41. 10.1099/mic.0.2007/008698-0.View ArticlePubMedGoogle Scholar
- Gumber S, Whittington RJ: Analysis of the growth pattern, survival and proteome of Mycobacteriumavium subsp. paratuberculosis following exposure to heat. Vet Microbiol. 2009, 136 (1-2): 82-90. 10.1016/j.vetmic.2008.10.003.View ArticlePubMedGoogle Scholar
- Gumber S, Taylor DL, Marsh IB, Whittington RJ: Growth pattern and partial proteome of Mycobacterium avium subsp. paratuberculosis during the stress response to hypoxia and nutrient starvation. Vet Microbiol. 2009, 133 (4): 344-357. 10.1016/j.vetmic.2008.07.021.View ArticlePubMedGoogle Scholar
- Wu CW, Schmoller SK, Shin SJ, Talaat AM: Defining the stressome of Mycobacterium avium subsp. paratuberculosis in vitro and in naturally infected cows. J Bacteriol. 2007, 189 (21): 7877-7886. 10.1128/JB.00780-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez GM: Control of iron metabolism in Mycobacterium tuberculosis. Trends Microbiol. 2006, 14 (7): 320-327. 10.1016/j.tim.2006.05.006.View ArticlePubMedGoogle Scholar
- Motiwala AS, Strother M, Amonsin A, Byrum B, Naser SA, Stabel JR, Shulaw WP, Bannantine JP, Kapur V, Sreevatsan S: Molecular epidemiology of Mycobacterium avium subsp. paratuberculosis: evidence for limited strain diversity, strain sharing, and identification of unique targets for diagnosis. J Clin Microbiol. 2003, 41 (5): 2015-2026. 10.1128/JCM.41.5.2015-2026.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Motiwala AS, Strother M, Theus NE, Stich RW, Byrum B, Shulaw WP, Kapur V, Sreevatsan S: Rapid detection and typing of strains of Mycobacterium avium subsp. paratuberculosis from broth cultures. J Clin Microbiol. 2005, 43 (5): 2111-2117. 10.1128/JCM.43.5.2111-2117.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh IB, Bannantine JP, Paustian ML, Tizard ML, Kapur V, Whittington RJ: Genomic comparison of Mycobacterium avium subsp. paratuberculosis sheep and cattle strains by microarray hybridization. J Bacteriol. 2006, 188 (6): 2290-2293. 10.1128/JB.188.6.2290-2293.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh IB, Whittington RJ: Genomic diversity in Mycobacterium avium: single nucleotide polymorphisms between the S and C strains of M. avium subsp. paratuberculosis and with M. a. avium. Mol Cell Probes. 2007, 21 (1): 66-75. 10.1016/j.mcp.2006.08.002.View ArticlePubMedGoogle Scholar
- Paustian ML, Zhu X, Sreevatsan S, Robbe-Austerman S, Kapur V, Bannantine JP: Comparative genomic analysis of Mycobacterium avium subspecies obtained from multiple host species. BMC Genomics. 2008, 9: 135-10.1186/1471-2164-9-135.PubMed CentralView ArticlePubMedGoogle Scholar
- Paustian ML, Kapur V, Bannantine JP: Comparative genomic hybridizations reveal genetic regions within the Mycobacterium avium complex that are divergent from Mycobacterium avium subsp. paratuberculosis isolates. J Bacteriol. 2005, 187 (7): 2406-2415. 10.1128/JB.187.7.2406-2415.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Turenne CY, Collins DM, Alexander DC, Behr MA: Mycobacterium avium subsp. paratuberculosis and M. avium subsp. avium are independently evolved pathogenic clones of a much broader group of M. avium organisms. J Bacteriol. 2008, 190 (7): 2479-2487. 10.1128/JB.01691-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Turenne CY, Semret M, Cousins DV, Collins DM, Behr MA: Sequencing of hsp65 distinguishes among subsets of the Mycobacterium avium complex. J Clin Microbiol. 2006, 44 (2): 433-440. 10.1128/JCM.44.2.433-440.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Alexander DC, Turenne CY, Behr MA: Insertion and deletion events that define the pathogen Mycobacterium avium subsp. paratuberculosis. J Bacteriol. 2009, 191 (3): 1018-1025. 10.1128/JB.01340-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu CW, Glasner J, Collins M, Naser S, Talaat AM: Whole-genome plasticity among Mycobacterium avium subspecies: insights from comparative genomic hybridizations. J Bacteriol. 2006, 188 (2): 711-723. 10.1128/JB.188.2.711-723.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Motiwala AS, Janagama HK, Paustian ML, Zhu X, Bannantine JP, Kapur V, Sreevatsan S: Comparative transcriptional analysis of human macrophages exposed to animal and human isolates of Mycobacterium avium subspecies paratuberculosis with diverse genotypes. Infect Immun. 2006, 74 (11): 6046-6056. 10.1128/IAI.00326-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Janagama HK, Jeong K, Kapur V, Coussens P, Sreevatsan S: Cytokine responses of bovine macrophages to diverse clinical Mycobacterium avium subspecies paratuberculosis strains. BMC Microbiol. 2006, 6: 10-10.1186/1471-2180-6-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez GM, Voskuil MI, Gold B, Schoolnik GK, Smith I: ideR, An essential gene in mycobacterium tuberculosis: role of IdeR in iron-dependent gene expression, iron metabolism, and oxidative stress response. Infect Immun. 2002, 70 (7): 3371-3381. 10.1128/IAI.70.7.3371-3381.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Seth M, Lamont EA, Janagama HK, Widdel A, Vulchanova L, Stabel JR, Waters WR, Palmer MV, Sreevatsan S: Biomarker discovery in subclinical mycobacterial infections of cattle. PLoS One. 2009, 4 (5): e5478-10.1371/journal.pone.0005478.PubMed CentralView ArticlePubMedGoogle Scholar
- Akkina SK, Zhang Y, Nelsestuen GL, Oetting WS, Ibrahlm HN: Temporal stability of the urinary proteome after kidney transplant: more sensitive than protein composition?. J Proteome Res. 2009, 8 (1): 94-103. 10.1021/pr800646j.PubMed CentralView ArticlePubMedGoogle Scholar
- Shilov IV, Seymour SL, Patel AA, Loboda A, Tang WH, Keating SP, Hunter CL, Nuwaysir LM, Schaeffer DA: The Paragon Algorithm, a next generation search engine that uses sequence temperature values and feature probabilities to identify peptides from tandem mass spectra. Mol Cell Proteomics. 2007, 6 (9): 1638-1655. 10.1074/mcp.T600050-MCP200.View ArticlePubMedGoogle Scholar
- Siegrist MS, Unnikrishnan M, McConnell MJ, Borowsky M, Cheng TY, Siddiqi N, Fortune SM, Moody DB, Rubin EJ: Mycobacterial Esx-3 is required for mycobactin-mediated iron acquisition. Proc Natl Acad Sci USA. 2009, 106 (44): 18792-18797. 10.1073/pnas.0900589106.PubMed CentralView ArticlePubMedGoogle Scholar
- Rao PK, Rodriguez GM, Smith I, Li Q: Protein dynamics in iron-starved Mycobacterium tuberculosis revealed by turnover and abundance measurement using hybrid-linear ion trap-Fourier transform mass spectrometry. Anal Chem. 2008, 80 (18): 6860-6869. 10.1021/ac800288t.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh A, Guidry L, Narasimhulu KV, Mai D, Trombley J, Redding KE, Giles GI, Lancaster JR, Steyn AJ: Mycobacterium tuberculosis WhiB3 responds to O2 and nitric oxide via its [4Fe-4S] cluster and is essential for nutrient starvation survival. Proc Natl Acad Sci USA. 2007, 104 (28): 11562-11567. 10.1073/pnas.0700490104.PubMed CentralView ArticlePubMedGoogle Scholar
- Abdallah AM, Verboom T, Hannes F, Safi M, Strong M, Eisenberg D, Musters RJ, Vandenbroucke-Grauls CM, Appelmelk BJ, Luirink J: A specific secretion system mediates PPE41 transport in pathogenic mycobacteria. Mol Microbiol. 2006, 62 (3): 667-679. 10.1111/j.1365-2958.2006.05409.x.View ArticlePubMedGoogle Scholar
- Gaballa A, Antelmann H, Aguilar C, Khakh SK, Song KB, Smaldone GT, Helmann JD: The Bacillus subtilis iron-sparing response is mediated by a Fur-regulated small RNA and three small, basic proteins. Proc Natl Acad Sci USA. 2008Google Scholar
- Jacques JF, Jang S, Prevost K, Desnoyers G, Desmarais M, Imlay J, Masse E: RyhB small RNA modulates the free intracellular iron pool and is essential for normal growth during iron limitation in Escherichia coli. Mol Microbiol. 2006, 62 (4): 1181-1190. 10.1111/j.1365-2958.2006.05439.x.View ArticlePubMedGoogle Scholar
- Masse E, Gottesman S: A small RNA regulates the expression of genes involved in iron metabolism in Escherichia coli. Proc Natl Acad Sci USA. 2002, 99 (7): 4620-4625. 10.1073/pnas.032066599.PubMed CentralView ArticlePubMedGoogle Scholar
- Bannantine JP, Huntley JF, Miltner E, Stabel JR, Bermudez LE: The Mycobacterium avium subsp. paratuberculosis 35 kDa protein plays a role in invasion of bovine epithelial cells. Microbiology. 2003, 149 (Pt 8): 2061-2069. 10.1099/mic.0.26323-0.View ArticlePubMedGoogle Scholar
- Bannantine JP, Radosevich TJ, Stabel JR, Berger S, Griffin JF, Paustian ML: Production and characterization of monoclonal antibodies against a major membrane protein of Mycobacterium avium subsp. paratuberculosis. Clin Vaccine Immunol. 2007, 14 (3): 312-317. 10.1128/CVI.00353-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Merkal RS, Curran BJ: Growth and metabolic characteristics of Mycobacterium paratuberculosis. Appl Microbiol. 1974, 28 (2): 276-279.PubMed CentralPubMedGoogle Scholar
- Motiwala AS, Li L, Kapur V, Sreevatsan S: Current understanding of the genetic diversity of Mycobacterium avium subsp. paratuberculosis. Microbes Infect. 2006, 8 (5): 1406-1418. 10.1016/j.micinf.2005.12.003.View ArticlePubMedGoogle Scholar
- Harris NB, Robbe-Austerman S, Payeur JB: Effect of egg yolk on the detection of Mycobacterium avium subsp. paratuberculosis using the ESP II liquid culture system. J Vet Diagn Invest. 2005, 17 (6): 554-560.View ArticlePubMedGoogle Scholar
- Whittington RJ, Sergeant ES: Progress towards understanding the spread, detection and control of Mycobacterium avium subsp paratuberculosis in animal populations. Aust Vet J. 2001, 79 (4): 267-278. 10.1111/j.1751-0813.2001.tb11980.x.View ArticlePubMedGoogle Scholar
- Wandersman C, Delepelaire P: Bacterial iron sources: from siderophores to hemophores. Annu Rev Microbiol. 2004, 58: 611-647. 10.1146/annurev.micro.58.030603.123811.View ArticlePubMedGoogle Scholar
- Masse E, Salvail H, Desnoyers G, Arguin M: Small RNAs controlling iron metabolism. Curr Opin Microbiol. 2007, 10 (2): 140-145. 10.1016/j.mib.2007.03.013.View ArticlePubMedGoogle Scholar
- Runyen-Janecky L, Daugherty A, Lloyd B, Wellington C, Eskandarian H, Sagransky M: Role and regulation of iron-sulfur cluster biosynthesis genes in Shigella flexneri virulence. Infect Immun. 2008, 76 (3): 1083-1092. 10.1128/IAI.01211-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Fontecave M, Choudens SO, Py B, Barras F: Mechanisms of iron-sulfur cluster assembly: the SUF machinery. J Biol Inorg Chem. 2005, 10 (7): 713-721. 10.1007/s00775-005-0025-1.View ArticlePubMedGoogle Scholar
- Huet G, Daffe M, Saves I: Identification of the Mycobacterium tuberculosis SUF machinery as the exclusive mycobacterial system of [Fe-S] cluster assembly: evidence for its implication in the pathogen's survival. J Bacteriol. 2005, 187 (17): 6137-6146. 10.1128/JB.187.17.6137-6146.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Savvi S, Warner DF, Kana BD, McKinney JD, Mizrahi V, Dawes SS: Functional characterization of a vitamin B12-dependent methylmalonyl pathway in Mycobacterium tuberculosis: implications for propionate metabolism during growth on fatty acids. J Bacteriol. 2008, 190 (11): 3886-3895. 10.1128/JB.01767-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Eoh H, Brown AC, Buetow L, Hunter WN, Parish T, Kaur D, Brennan PJ, Crick DC: Characterization of the Mycobacterium tuberculosis 4-diphosphocytidyl-2-C-methyl-D-erythritol synthase: potential for drug development. J Bacteriol. 2007, 189 (24): 8922-8927. 10.1128/JB.00925-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Miallau L, Faller M, Chiang J, Arbing M, Guo F, Cascio D, Eisenberg D: Structure and proposed activity of a member of the VapBC family of toxin-antitoxin systems. VapBC-5 from Mycobacterium tuberculosis. J Biol Chem. 2009, 284 (1): 276-283. 10.1074/jbc.M805061200.PubMed CentralView ArticlePubMedGoogle Scholar
- Stallings CL, Stephanou NC, Chu L, Hochschild A, Nickels BE, Glickman MS: CarD is an essential regulator of rRNA transcription required for Mycobacterium tuberculosis persistence. Cell. 2009, 138 (1): 146-159. 10.1016/j.cell.2009.04.041.PubMed CentralView ArticlePubMedGoogle Scholar
- Lella RK, Sharma C: Eis (enhanced intracellular survival) protein of Mycobacterium tuberculosis disturbs the cross regulation of T-cells. J Biol Chem. 2007, 282 (26): 18671-18675. 10.1074/jbc.C600280200.View ArticlePubMedGoogle Scholar
- Frankenberg-Dinkel N: Bacterial heme oxygenases. Antioxid Redox Signal. 2004, 6 (5): 825-834.View ArticlePubMedGoogle Scholar
- Chim N, Iniguez A, Nguyen TQ, Goulding CW: Unusual diheme conformation of the heme-degrading protein from Mycobacterium tuberculosis. J Mol Biol. 395 (3): 595-608. 10.1016/j.jmb.2009.11.025.
- Boughammoura A, Matzanke BF, Bottger L, Reverchon S, Lesuisse E, Expert D, Franza T: Differential role of ferritins in iron metabolism and virulence of the plant-pathogenic bacterium Erwinia chrysanthemi 3937. J Bacteriol. 2008, 190 (5): 1518-1530. 10.1128/JB.01640-07.PubMed CentralView ArticlePubMedGoogle Scholar