FurA contributes to the oxidative stress response regulation of Mycobacterium avium ssp. paratuberculosis.

The ferric uptake regulator A (FurA) is known to be involved in iron homeostasis and stress response in many bacteria. In mycobacteria the precise role of FurA is still unclear. In the presented study, we addressed the functional role of FurA in the ruminant pathogen Mycobacterium avium ssp. paratuberculosis (MAP) by construction of a furA deletion strain (MAPΔfurA). RNA deep sequencing revealed that the FurA regulon consists of repressed and activated genes associated to stress response or intracellular survival. Not a single gene related to metal homeostasis was affected by furA deletion. A decisive role of FurA during intracellular survival in macrophages was shown by significantly enhanced survival of MAPΔfurA compared to the wildtype, indicating that a principal task of mycobacterial FurA is oxidative stress response regulation in macrophages. This resistance was not associated with altered survival of mice after long term infection with MAP. Our results demonstrate for the first time, that mycobacterial FurA is not involved in the regulation of iron homeostasis. However, they provide strong evidence that FurA contributes to intracellular survival as an oxidative stress sensing regulator.

INTRODUCTION

For efficient infection and colonization, pathogenic bacteria need to adapt their metabolism to the host environment and to combat innate antimicrobial host defense mechanisms such as iron limitation by iron sequestering molecules such as ferritin, transferrin and lactoferrin, and the production of toxic reactive oxygen (ROS) and nitrogen species (RNS).

Reactive oxygen and nitrogen species comprise highly reactive molecules such as the superoxide anion (O2-), hydroxyl radical (OH-), and peroxynitrite (ONOO-) as well as more stable oxidants like hydrogen peroxide (H2O2) and nitric oxide (NO). These products are generated during auto oxidative metabolic processes and play important roles in prokaryotic and eukaryotic homeostasis and cell signaling (Yoshikawa and Naito, 2000; Devasagayam et al., 2004; Gough and Cotter, 2011). However, when enriched in infected phagocytes by the activity of phagosomal myeloperoxidases and NO-oxidases, they represent effective antimicrobial measures. Accordingly, intracellular persisting mycobacteria have developed host defense escape mechanisms and the regulation of these is often linked to the expression of virulence factors.

Like other pathogenic mycobacteria, Mycobacterium avium ssp. paratuberculosis (MAP), the causative agent of Johne’s disease, a chronic, non-curable granulomatous inflammation of the ruminant intestine (Harris and Barletta, 2001; Lombard, 2011), is able to persist and multiply intracellularly in the phagosome of macrophages (Kuehnel et al., 2001). Specifically, MAP adapts to the intestinal environment by adjusting its metabolism to the host environment and responding to antimicrobial defense mechanisms indicated by enhanced expression of KatG and SodA (Weigoldt et al., 2011, 2013). Hence, MAP is clearly able to withstand the host innate defense response, but its enabling mechanisms are still unclear.

In general, oxidative stress response in bacteria is regulated by the major transcription factors OxyR and PerR (Zheng et al., 1999; Herbig and Helmann, 2001). OxyR belongs to the family of LysR-regulators and is not only found in Gram negative bacteria but also in some Gram positive bacteria (Oh et al., 2007). OxyR is a transcriptional activator under oxidizing conditions acting through direct interaction with the RNA polymerase α subunit (Kullik et al., 1995; Tao et al., 1995). PerR is a member of the ferric uptake regulator (FUR) family of metallo-regulators and functions as a peroxide-responsive repressor (Herbig and Helmann, 2001). PerR is found mostly in Gram positive bacteria where it seems to substitute for OxyR (Mongkolsuk and Helmann, 2002). Nevertheless, in some bacterial species PerR is present together with OxyR (van Vliet et al., 1999; Hahn et al., 2002; Tseng et al., 2003). To prevent damage by ROS, OxyR, and/or PerR regulate, amongst others, the expression of a catalase-peroxidase (kat) and an alkyl hydroperoxide reductase (ahp) in many bacteria (Christman et al., 1989; Storz and Altuvia, 1994; Bsat et al., 1996, 1998).

Iron plays a particular role in bacterial metabolism. On the one hand, iron is an essential structural and catalytic cofactor for many metabolic enzymes. On the other hand, iron excess in the cell is harmful due to the spontaneous generation of ROS by the Fenton reaction. Therefore, the bacterial intracellular iron homeostasis must be tightly regulated and often is closely linked to the response to oxidative stresses (Cornelis et al., 2011). In Mycobacterium tuberculosis (Mtb), three genes encoding putative iron-dependent regulatory proteins have been described (Cole et al., 1998): IdeR (Schmitt et al., 1995) and SirR (Hill et al., 1998), both belonging to the diptheria-toxin repressor (DtxR) family, as well as ferric uptake regulator A (FurA), a member of the FUR family (Lee and Helmann, 2007; Fillat, 2014).

IdeR is the global iron-dependent regulator in mycobacteria. As an iron-loaded homodimer, it controls the expression of about 50 genes in Mtb which are mostly involved in iron homeostasis, and in addition a number of virulence associated factors (Rodriguez et al., 2002). The ideR gene is essential in pathogenic mycobacteria like Mtb (Rodriguez et al., 2002) or MAP (Janagama et al., 2009), but not in saprophytic species such as M. smegmatis (Dussurget et al., 1996). Using a conditional ideR mutant, Pandey and Rodriguez (2014) recently demonstrated that IdeR is a factor important for the survival of Mtb in macrophages as well as for the infection of mice. SirR was originally identified in Staphylococcus epidermidis with a sequence identity to Mtb IdeR of 33% (Hill et al., 1998). The biological role of SirR homologs in pathogenic mycobacteria is still unclear.

The FUR family of transcriptional metallo-regulator proteins is widespread in bacteria and archaea (Fillat, 2014). The best characterized member of the FUR family, the Escherichia coli Fur protein, employs Fe2+ as co-repressor and acts as a global regulator, influencing the expression of more than 100 genes involved in iron homeostasis, intermediary metabolism, and oxidative stress response (McHugh et al., 2003). The classical mechanism of Fur regulation includes Fur-binding of Fe2+, formation of a metal-containing dimer which subsequently binds to a specific Fur-box found in the promoter sequences of the corresponding target genes (Bagg and Neilands, 1987). Recent studies indicated that the regulatory mechanisms governed by Fur exceed this classical repressor model. Thus, Fur homologs in different bacterial species can also act as metal-dependent positive regulators or even repress gene expression in the absence of the iron cofactor [apo-Fur repression; see review by Carpenter et al. (2009)]. In addition, a metal-independent activation of genes by apo-Fur was recently described (Butcher et al., 2012).

In all mycobacterial species sequenced so far furA is located upstream of katG. Furthermore, for some mycobacteria species, co-expression of furA and katG indicated an involvement of FurA in mycobacterial oxidative stress response. Indeed, it was proposed that, in mycobacteria, FurA could functionally replace the peroxide repressor PerR (Lucarelli et al., 2008). Nevertheless, the precise regulatory role of FurA and the regulation of stress response in mycobacteria are still not completely clear. One explanation for this might be the absence of a functional OxyR regulator in some mycobacterial species like those of the M. tuberculosis complex and M. smegmatis. In these species, it was found that FurA is involved in the oxidative stress response by regulating katG (Pym et al., 2001; Zahrt et al., 2001) but there was no contribution of FurA to ahpC expression. In contrast, many other mycobacterial species, including the pathogens M. avium and M. marinum, harbor a functional oxyR gene (Deretic et al., 1995; Sherman et al., 1995). In these species it was shown that ahpC but not katG is under control of OxyR (Dhandayuthapani et al., 1997; Pagan-Ramos et al., 1998). However, the regulation of katG in these species has not been addressed. Overall, these data suggest a more complex regulatory network of oxidative stress response regulation in mycobacteria.

In the present work we aimed to analyze the relevance of mycobacterial FurA in iron homeostasis and stress response. We characterized the FurA regulon in MAP by construction and analysis of an isogenic MAPΔfurA mutant. By this, we were able to clarify the influence of FurAMAP in iron metabolism, stress response, and intracellular survival.

MATERIALS AND METHODS

CHEMICALS, BACTERIAL STRAINS, AND GROWTH CONDITIONS

All chemicals were purchased from Sigma–Aldrich (Munich, Germany) if not stated otherwise. Bacterial strains, oligonucleotides, plasmids, and phages are listed in Tables S1 and S2. E. coli strains were cultivated at 37°C in Luria-Bertani (LB) broth or on LB-agar supplemented with 100 μg/ml ampicillin or 100 μg/ml hygromycin, as required.

Mycobacterium smegmatis mc2 155 and transformants were grown in LB media or on Middlebrook phage (MBP) agar (Supplementary Methods) supplemented with 100 μg/ml hygromycin if necessary. Liquid cultures were incubated at 37°C in a shaking incubator at 100 rpm. MAP DSM44135 wildtype (wt), MAPΔfurA, and the complemented strain MAPΔfurAC were grown in DifcoTM Middlebrook 7H9 broth or on solid Middlebrook 7H10 agar (Beckton Dickinson, Franklin Lakes, NJ, USA) supplemented with mycobactin J (2 mg/l, Allied Monitor, Fayette, MO, USA), 2.5% glycerol and 10% OADC (0.06% oleic acid, 5% albumin, 2% dextrose, 0.085% NaCl, 0.003% catalase; MB-complete broth) and 50 μg/ml hygromycin and/or 50 μg/ml kanamycin, as required. Agar plates were sealed in plastic bags after inoculation and incubated at 37°C for 8–10 weeks. Liquid cultures were grown in Duran® laboratory glass bottles on a magnetic stirrer at 100 rpm and 37°C until they reached an OD600 of 1.0 and used for further experiments.

For growth kinetics, a starter culture was grown to an OD600 of 1.0 in MB-complete broth. Bacteria were harvested by centrifugation at 3,800 ×g for 15 min at 4°C, resuspended in MB-complete broth, singularized with glass beads (Ø 3 mm) by vortexing and inoculated into fresh medium to obtain an initial OD600 of 0.2. The OD600 of cultures was measured twice a week and bacteria were grown until they entered the stationary phase.

The response of MAPwt and MAPΔfurA to oxidative stress was examined upon exposure of bacteria to hydrogen peroxide (H2O2). Cells were grown in MB-complete broth and harvested at an OD600 of 1.0 as described above. The bacterial pellet was washed twice with PBS and singularized cells were transferred to 50 ml Middlebrook 7H9 broth supplemented with 2 mg/l mycobactin, 2.5% glycerol, and 10% ADC lacking catalase (5% albumin, 2% dextrose, 0.085% NaCl; MB-cat broth) to obtain an OD600 of 1.0. H2O2 was added to a final concentration of 10 mM and cultures were exposed for 2 h at 37°C with gentle shaking (100 rpm) in Erlenmeyer flasks. Following, bacteria were harvested and immediately used for RNA-extraction.

CONSTRUCTION OF A MAPΔfurA STRAIN AND COMPLEMENTATION OF THE MUTANT STRAIN

Construction of a furA (map1669c) deletion mutant of MAPwt was performed by specialized transduction according to Park et al. (2008) with minor modifications. A detailed protocol is given in the Supplementary Methods file.

For complementation of the mutant strain the furA gene, including 88 bp of the upstream UTR, was amplified by PCR using primers oMAP-furA-K1/oMAP-furA-K2. The PCR product was digested with XbaI and HindIII and ligated into the corresponding sites of the mycobacterial integrative vector pMV306 (Stover et al., 1991), carrying a kanamycin resistance. The resulting plasmid pMAP-furA1101 was transformed into electro-competent MAPΔfurA cells and plated on MB-complete agar with 50 μg/ml hygromycin and kanamycin. Successful deletion and integration of furA was confirmed by PCR.

CELL CULTURE OF MACROPHAGES AND SURVIVAL ASSAY OF INTRACELLULAR MYCOBACTERIA

The intracellular survival of MAPwt, MAPΔfurA, and MAPΔfurAC was determined in the mouse macrophage cell line J774A.1 (Ralph et al., 1975). Macrophage cells were seeded in tissue culture dishes and were maintained in Dulbecco’s modified Eagle medium (DMEM) supplemented with 10% fetal calf serum (FCS), 1% glutamine, 100 units/ml penicillin, 100 μg/ml streptomycin at 37°C and 8% CO2. Twenty-four hours prior to infection, medium was changed to antibiotic-free DMEM. For infection experiments, MAP cells were grown to an OD600 of 1.0, harvested, singularized and stored in MB-complete broth containing 10% glycerol at -80°C until further usage. For infection, MAP cells were thawed on ice, diluted in antibiotic-free DMEM to an OD600 of 0.15 and incubated with the macrophages (2.0 × 106 cells per cell culture dish) as described earlier (Kuehnel et al., 2001), resulting in a MOI of 10. The mean values of the initial inocula (five experiments) of all strains (MAPwt, MAPΔfurA, MAPΔfurAC) are shown in the Supplementary data files. At indicated time points, infected macrophages were washed twice with PBS and scraped off the plates in 1 ml PBS containing 0.1% SDS. J774A.1 cells were disrupted by passage through a 24-gauge needle. The homogenates as well as the infection inocula were 10-fold serial diluted in PBS and 25 μl of the 10-4 and 10-5 dilutions were plated on supplemented MB-complete agar plates in duplicate. Colony forming units (CFU) were determined after incubation for up to 10 weeks at 37°C and normalized to the initial inoculum.

DETERMINATION OF REACTIVE OXYGEN SPECIES IN MACROPHAGE INFECTION EXPERIMENTS

Induction of ROS production by macrophages upon infection with MAPwt, MAPΔfurA, and MAPΔfurAC was determined by use of CellRox® Deep Red (Life Technologies), according to the manufacturer’s protocol. In brief, macrophages were cultured and infected as described above and incubated for 2 h at 37°C/8% CO2 or treated with the chemical menadione as a positive control at a final concentration of 100 μM (Malorni et al., 1993). CellRox® Deep Red was added to a final concentration of 2.5 μM 30 min before harvesting and detection. Signals were measured at 660 nm by flow cytometry (Guava® easyCyte 8HT, Millipore, Billerica, MA, USA) and shown as mean fluorescence intensity (MFI) of four independent experiments.

IMMUNOFLUORESCENCE AND CONFOCAL MICROSCOPY

Macrophages were seeded on coverslips and infected as described above. After 30 min and 2 h of infection, cells were fixed with 3% formaldehyde in PBS for 10 min, washed twice with PBS and treated with blocking buffer (PBS/1% BSA/10% FCS) for 20 min at room temperature. For staining of extracellular mycobacteria, coverslips were incubated with a 1:100 dilution of a polyclonal rabbit anti-MAP-HBHA serum (Supplementary Methods) in blocking buffer at room temperature in a humid chamber for 45 min. After washing with PBS, coverslips were incubated with goat anti-rabbit IgG coupled to Alexa Fluor® 488 (Life Technologies, Darmstadt, Germany) for 30 min and subsequently washed with PBS. To label intracellular mycobacteria, cells were permeabilized with Triton X-100 (0.1%) for 5 min at room temperature, washed in PBS, followed by incubation with a 1:100 dilution of polyclonal rabbit anti-MAP-HBHA serum in blocking buffer for 45 min at room temperature. After washing, coverslips were treated with goat anti-rabbit IgG coupled to Alexa Fluor® 568 (Life Technologies, Darmstadt, Germany) for 30 min. Coverslips were washed three times in PBS and, after one brief washing step with ddH2O, samples were mounted using ProLong® Gold with DAPI (Life Technologies, Darmstadt, Germany). Mounted samples were examined using a TCS SP5 confocal laser scanning microscope equipped with a 63×/1.4–0.6 NA HCX PL APO objective (Leica, Wetzlar, Germany). Image stacks were acquired using 1 Airy unit pinhole diameter in sequential imaging mode to avoid bleed through. Image stacks were deconvolved using Huygens® Essential (Scientific Volume Imaging, Hilversum, The Netherlands), maximum intensity projections were calculated for display purposes and adjusted for brightness and contrast using ImageJ/Fiji (Schindelin et al., 2012). After this labeling procedure, extracellular mycobacteria appear yellow, whereas intracellular mycobacteria are stained red.

ANIMAL EXPERIMENT

Mouse infection experiments were approved by the Lower Saxony Federal State Office for Consumer Protection and Food Safety, Germany (reference number 08/1504). In each group, 9 female C57BL/6 mice aged 8 weeks (Charles River, Erkrath, Germany) were challenged intraperitoneally with an infection dose of 1 × 108 bacteria of MAPwt and MAPΔfurA in 200 μl Dulbecco’s Phosphate-Buffered Saline (DPBS, Life Technologies GmbH, Darmstadt, Germany). DPBS was used as negative control. Mice were sacrificed after 4 weeks, and liver and spleen were weighed. Bacteria were quantified by plating of homogenized tissue on MB-complete agar plates. Histology was performed in the Mouse Pathology Platform at HZI Braunschweig as recently described (Meißner et al., 2014; Suwandi et al., 2014).

EXTRACTION OF NUCLEIC ACIDS AND QUANTITATIVE REAL-TIME PCR (qRT-PCR)

Genomic DNA was prepared as previously described (Stratmann et al., 2004). Plasmids were prepared using NucleoBond® AX kit (Macherey Nagel GmbH, Düren, Germany) according to the manufacturer’s protocol. Total RNA from bacteria was isolated according to Rustad et al. (2009) with minor modifications. In brief, bacterial pellets were resuspended in TRIzol® reagent, mechanically disrupted in a FastPrep® instrument (ThermoSavant, Carlsbad, CA, USA) and RNA was separated by chloroform, chloroform-isoamylalcohol (49:1 v/v) extraction, and precipitation with 2-propanol. Total RNA was treated twice with 50 U DNase I (Roche, Mannheim, Germany) according to the manufacturer’s protocol. Integrity and quality of RNA was determined by gel electrophoresis and spectrophotometric analysis using an Epoch instrument (Biotek, Bad Friedrichshall, Germany) at 260 nm.

For cDNA synthesis, 4 μg of DNA-depleted RNA were diluted with RNase free water up to a volume of 20 μl and incubated for 10 min at 70°C with 0.4 μg random primers (Promega, Madison, WI, USA). After 5 min cooling on ice, samples were split in two and mixed with 5x RT buffer, 10 mM dNTP’s, and either 200 U MMLV-superscript transcriptase (Promega) or buffer without transcriptase as a negative control. The mix was incubated for 1 h at 42°C, following incubation for 5 min at 85°C. Samples were diluted with 90 μl ddH2O and stored at -20°C for further analysis. 2.5 μl of each cDNA sample was mixed with 0.4 μM primer and 10 μl SYBR-Green Mix (Qiagen, Düsseldorf, Germany) in a total volume of 20 μl. Quantitative real-time PCR (qRT-PCR) experiments were performed using a Mx3005P qPCR system (Agilent Technologies, Santa Clara, CA, USA) with a thermal cycling profile as follows: segment 1, 20 min at 95°C, 1 cycle; segment 2, 45 s at 95°C, 1 min at 58°C, 1 min at 72°C, 45 cycles; segment 3, 1 min 95°C, 30 s 55°C, 30 s 95°C, 1 cycle. Each sample was analyzed in duplicate. Ct values were shown either as absolute data or normalized to the housekeeping gene gap and expressed as fold change to the untreated control (ΔΔct). All oligonucleotide primer pairs (Table S1) were tested for efficacy with a serial dilution of genomic DNA.

RNA DEEP SEQUENCING AND ANALYSIS

RNA deep sequencing was used to determine the regulon of FurA in MAP. Sequencing libraries of MAPwt and MAPΔfurA RNA, of three independent samples each, were prepared and sequenced using 36 bp single-ends sequencing on a Genome Analyzer IIx (Illumina, San Diego, CA, USA) or 50 bp single-ends sequencing on a HiSeq2500 instrument (Illumina), respectively. In brief, libraries of 300 bp were prepared according the manufacturer’s instructions “Preparing Samples for Sequencing of mRNA” (Illumina) and “Protocol for ScriptSeqTM v2 RNA-Seq Library Preparation” (Epicentre Biotechnologies, Madison, WI, USA). Quality of the libraries was validated using an Agilent Bioanalyzer (Agilent Technologies) following the manufacturer’s instruction. Cluster generation was performed using the Illumina cluster station; sequencing on the Genome Analyzer IIx or HiSeq2500 followed a standard protocol. The fluorescent images were processed to sequences and transformed to FastQ format using the Genome Analyzer Pipeline Analysis software 1.8.2 (Illumina). The sequence output was controlled for general quality features using the fastq-mcf tool of ea-utils (Aronesty, 2013) and was mapped against the genome sequence of the reference strain MAP-K10 (AE016958) using BWA version 0.7.5 (Li and Durbin, 2009) and SAMtools (Li et al., 2009) for storing nucleotide sequence alignments. Strain MAP-K10 was used as a reference strain as the genome of MAP DSM44135 has not been sequenced yet. Raw data sets are available in the European Nucleotide Archive Repository [1]. For data analyses, reads of the 50 bp sequencing were clipped to 36 bp reads. Following, all sequences were computed with Rockhopper tool (McClure et al., 2013). Genes with a q-value ≤0.01 were considered as significantly differentially expressed and genes with raw counts in all replicates were included for further analysis (Table S3). The fold change was calculated as the expression value of MAPwt/MAPΔfurA.

BIOINFORMATICS AND STATISTICS

Differentially expressed genes identified by Rockhopper analysis were further processed with Blast2Go tool (Conesa et al., 2005) and NCBI blastx [2] as well as tuberculosis database [3] to investigate functions of putative proteins. Additionally, the putative proteins were clustered into orthologous groups (COG), based on the genome of MAP-K10 [4]. Data are expressed as mean ± SEM, and statistical analyses were performed using GraphPad Prism 5.03 (GraphPad, San Diego, CA, USA). Depending on the experiment, either the non-parametric t-test (Mann–Whitney) or one-way ANOVA test (Kruskal–Wallis) were used. A p-value of <0.05 between samples and controls was considered as statistically significant.

Prediction of putative FurA binding sites in MAP (NCBI accession No. NC_002944) was performed by FIMO analysis (Bailey et al., 2009). Initially, we generated a consensus sequence for FurA binding by comparison of sequences located -70/+30 bp up-/downstream of predicted FurA translation start sites in Mtb (NC_000962, Rv1909c), M. bovis BCG pasteur (NC_008769, BCG_1948c), M. avium ssp. avium 104 (NC_008595, MAV_2752), MAP (NC_002944, MAP_1669c), and M. smegmatis (NC_008596, MSMEG_6383) with MEME Suite (Bailey et al., 2009). The resulting binding motif T[CT]TTGACT[CG][AG]TTCCA[GA]A[AT]AA[GT][GT][GC][ACG][GC][TG]CAT[AT] was subsequently submitted to FIMO analysis. Nucleotides in brackets are variable, single nucleotides are conserved.

RESULTS

CHARACTERIZATION OF A MAPΔfurA DELETION MUTANT

Ferric uptake regulator A (FurA) has been described as a global regulator in E. coli and was proposed to be involved in oxidative stress response and regulation of iron homeostasis (Zheng et al., 1999). In mycobacteria, iron homeostasis is maintained by IdeR, and FurA is assumed to take part in general stress response due to its genomic co-localization and its (partly) co-transcription with katG. However, the regulatory role of FurA has not yet been fully understood.

To gain insight into the role of FurA in MAP, we generated a furA (map1669c) deletion mutant MAPΔfurA by specialized transduction. In addition, a complemented strain was constructed by introducing the integrative plasmid pFurA-MAP1101, harboring the furA gene under control of its own promoter, into the ΔfurA mutant strain, resulting in strain MAPΔfurAC. Deletion of the furA gene in MAPΔfurA and complementation of MAPΔfurAC were confirmed by PCR (Figure ). Expression of furA and the adjacent gene katG were analyzed by qRT-PCR. furA transcripts were detected in MAP wildtype (MAPwt) and MAPΔfurAC but not in MAPΔfurA. The expression levels of katG were similar in all strains (Figure ), indicating that the genetic manipulation had no polar effect on the expression of downstream located genes. Next we compared growth of MAPwt, MAPΔfurA, and MAPΔfurAC in MB-complete broth. As shown in Figure , all strains exhibited similar growth kinetics. Growth of MAPΔfurA and MAPΔfurAC was slightly reduced in comparison to MAPwt, and both strains reached stationary growth phase at slightly lower OD600 than MAPwt. These effects might be due to the antibiotic supplements in the growth media.

Characterization of a MAPΔ Confirmation of deletion and complementation of furA. Chromosomal DNA was prepared from MAPwt (lanes 1), MAPΔfurA (lanes 2), and MAPΔfurAC (lanes 3) grown in MB-complete broth to an OD600 of 1.0 and was subsequently tested in PCR for the presence of furA (left panel) with primers oRTfurA1 fw/rev and hygromycin (right panel) with primers oTKfurAKatGfw/oRTlppsrev. A negative control with water only is shown in lanes 4. M1 = 100 bp DNA ladder (Invitrogen), M2 = 1 kb plus DNA ladder. Samples were run on a 1.5% agarose gel. (B) Total RNA was extracted from all strains grown to an OD600 of 1.0 in MB-complete broth. Expression of furA (black bars) and the adjacent gene katG (white bars) was analyzed by qRT-PCR. Results are expressed as mean Ct values of three independent experiments. The level of the housekeeping gene gap for all samples was 19.27 ± 0.4060. (C) Growth of MAPwt (filled triangles), MAPΔfurA (open rectangle), and the complemented strain MAPΔfurAC (filled rectangle) was monitored for 60 days by measurement of the optical density (OD600) in MB-complete broth. Each growth experiment was conducted in triplicate for each strain. Shown are representative results of one out of two independent experiments.

DETERMINATION OF THE FurA REGULON OF MAP

To analyze the regulatory role of FurA in MAP, we prepared RNA from MAPwt and MAPΔfurA grown to an OD600 of 1.0 in MB-complete broth and performed RNA deep sequencing. Gene expression profiles of three independent replicates of each strain were analyzed with the Rockhopper analysis tool (McClure et al., 2013), and genes with a q-value of ≤0.01 were considered as significantly differentially expressed. Differences in gene expression were calculated as fold change of expression values calculated by Rockhopper analysis (MAPΔfurA vs. MAPwt).

In total, we identified 48 genes, which were differentially expressed in the mutant strain. Compared to the wt strain, 13 genes were higher (Table ) and 35 lower (Table ) expressed in MAPΔfurA. Clustering of putative proteins encoded by the differentially expressed genes into orthologous groups (COG) revealed that almost 40% of the differentially expressed genes were assigned to metabolism, 20.75% to cellular processes and signaling, 9.43% belong to the group of poorly characterized proteins, and 30.19% were not assigned to any COG (Tables and ). No genes related to iron homeostasis, such as siderophores (mycobactin) and siderophore uptake systems (Neilands, 1992; Escolar et al., 1998; McHugh et al., 2003) or iron storage proteins bfr, fntA (Abdul-Tehrani et al., 1999; Wildermann et al., 2004), were affected by the furA deletion in MAP.

Among the 13 higher expressed genes, 11 were organized in 4 putative operons. The majority of these genes encodes for proteins predicted to be involved in stress response or redox processes. Within this group, the gene cluster map1589c–1587c, encoding for two alkyl hydroperoxide reductases involved in oxidative stress and drug resistance as well as intracellular survival in Mtb (Sherman et al., 1999; Master et al., 2002), and a NAD(P)H nitroreductase (map1743c), which has been shown to be associated with detoxification of nitroaromatic components in Mtb (Purkayastha et al., 2002), exhibited the most significant expression changes (Table ).

The group of lower expressed genes is comprised of 35 genes. Aside from furA, this group consists of two putative transporters, five genes predicted to be involved in metabolism, three stress response-associated genes and 13 genes were of unknown function. The remaining 11 genes can be associated to virulence which is represented by two putative invasion proteins (map1203/1204), two putative antigens (map0047c/map2168c), one Pro-Pro-Glu protein (PPE; map1003c), one PE-PGRS glycin-rich protein (map3868), one nlpc p60 protein (map1272c), and one mammalian cell entry (mce) family protein (map2114c) as well as three other enzymes (map1706, map1967c, map3421c) for cell entry and persistence. Altogether, map0047c, map0847 (conserved protein), and map4206c (putative efflux ABC transporter) showed the most significant differences in expression compared to the wt (Table ). In addition, the 5′ UTRs of map1589c (ahpC), map0081 (putative oxidoreductase), map0847, and map3421c (putative diguanylate cyclase phosphodiesterase) featured a predicted FurA binding site (Tables and ). Overall, these analyses indicated a considerable involvement of MAP FurA in the adaptation of MAP to the host cell environment and general stress response.

To confirm the results of our RNA deep sequencing analyses and to further specify the involvement of FurA in regulation, we analyzed the transcription of ahpC and ahpD (map1589c/1588c, higher expressed in MAPΔfurA), and of map0847 and map0047c (lower expressed in MAPΔfurA) by qRT-PCR. MAPwt, MAPΔfurA, and the complemented strain MAPΔfurAC were grown in MB-complete broth to an OD600 of 1.0, RNA was extracted and analyzed by qRT-PCR. As a control, we included the mbtB gene (map2177c) which has been shown to be regulated iron-dependently by IdeR in Mtb (Gold et al., 2001) and MAP (Janagama et al., 2009). As expected, the expression of mbtB was similar in all strains. In agreement with our RNA sequencing analyses, ahpC, ahpD, map0847, and map0047c were up- and down-regulated in the mutant, respectively. In the complemented mutant MAPΔfurAC expression levels were mostly restored to wt levels (Figure ). The differences to the wt strain might be due to the slightly lower expression level of furA in MAPΔfurAC compared to MAPwt (Figure ). Overall, these results confirm the regulatory influence of FurA on expression of ahpC, ahpD, map0847, and map0047c.

Functional complementation of FurA-regulated genes. MAPwt, MAPΔfurA and the complemented strain MAPΔfurAc were grown in MB-complete broth to an OD600 of 1.0. Total RNA was extracted, followed by cDNA synthesis and qRT-PCR analysis. Expression of genes mbtB, ahpC (map1589c), ahpD (map1588c), map0847, and map0047c in MAPΔfurA (white bars) and MAPΔfurAc (black bars) was compared to transcription levels in MAPwt. Three independent experiments were analyzed, qRT-PCR was performed in duplicate for each sample, normalized to the housekeeping gene gap and expressed as fold change compared to the MAPwt (mean ± SEM). Statistical analysis was performed using Mann–Whitney test with *p < 0.05, **p < 0.005.

ahp EXPRESSION IN MAP IS HYDROGEN PEROXIDE DEPENDENT

AhpC and AhpD are metal dependent alkyl hydroperoxide reductases and constitute an antioxidant system important in reduction and neutralization of a wide range of organic hydroperoxides and peroxynitrite (Wood et al., 2003). Genes encoding for AhpC/D have been shown to be induced by peroxide stress mediated by OxyR in other mycobacteria (Dhandayuthapani et al., 1997). However, in MAP these genes seem to be regulated by FurA which resembles PerR mediated activation in Bacillus and Campylobacter (van Vliet et al., 1999; Helmann et al., 2003). Thus, we exposed MAPwt and MAPΔfurA cultures without catalase supplement to oxidative stress by addition of 10 mM H2O2 for 2 h and analyzed ahpD expression by qRT-PCR. As a negative control we used the FurB dependent but FurA independent regulated mptA (map3736c) gene (Eckelt et al., 2014).

As shown in Figure , expression of ahpD in MAPΔfurA was increased compared to the wt and supports RNA-sequencing data. Upon peroxide treatment, both strains responded with an induction of ahpD expression. However, the relative fold change in MAPΔfurA was considerably higher than in MAPwt compared to the untreated control.

Response of the FurA regulated alkyl hydroperoxide reductase system AhpCD to oxidative stress. MAPwt and MAPΔfurA cultures were grown in MB-complete broth, harvested at OD600 of 1.0, washed, transferred to MB-cat broth and incubated for 2 h with or without 10 mM H2O2. RNA was isolated and used for cDNA synthesis. Gene expression analysis by qRT-PCR was performed for the genes ahpD (map1588c) and mptA (map3736c) in untreated (MAPwt, white bars; ΔfurA, dark gray bars) and H2O2 treated samples (MAPwt, light gray bars; ΔfurA, black bars). Results were normalized to the housekeeping gene gap and are expressed as fold change compared to the untreated MAPwt control (mean ± SEM, n = 3).

In contrast mptA expression was similar in all strains at both conditions (standard/oxidative stress).

These data on the one hand side indicate that repression of ahpC/D by FurA can be relieved by peroxide. On the other hand they demonstrate that an additional peroxide sensing regulatory activity, most probably by OxyR or a OxyR-like regulator, co-regulates ahpC/D expression in MAP.

MAPΔfurA EXHIBITS HIGHER RESISTANCE TO INTRACELLULAR KILLING IN MACROPHAGES BUT FurA IS DISPENSABLE DURING CHRONIC INFECTION

It has been shown that ahpC/D are necessary for defense against ROS and RNS in the mycobacterial phagosome (Sherman et al., 1995; Wilson et al., 1998; Master et al., 2002). By this pathogenic mycobacteria are able to survive within macrophages, enabling metabolic adaptation to the phagosomal compartment which is necessary for persistence. The above data showed that FurA is involved in the regulation of MAP response to oxidative stress, as indicated by the regulation of enzymes involved in neutralization of oxygen and nitrogen radicals [ahpC/D, NAD(P)H nitroreductase]. Therefore, we aimed to address the role of FurA for intracellular survival of MAPwt, MAPΔfurA, and MAPΔfurAC in macrophages. First we confirmed that the MAP strains are exposed to ROS after being taken up by the macrophages. For this, J774A.1 macrophages were infected with MAPwt, MAPΔfurA, and MAPΔfurAC. After 1.5 h of infection, cells were treated with the reagent CellROX Deep Red for 30 min. This dye fluoresces upon oxidation and allows the detection of ROS. As positive control we included macrophages treated (+) with menadione (Figure ) which is known to induce the generation of ROS in eukaryotic cells (Malorni et al., 1993). As shown in Figures , the MFI measured by FACS analysis in macrophages infected with MAPwt, MAPΔfurA and the complemented mutant strain was higher than in untreated (-) macrophages (Figure ), indicating the generation of ROS induced by all MAP strains. Next, cells were infected with MAPwt, MAPΔfurA, and MAPΔfurAC as described above and bacterial survival was determined after 2 h, 2 and 7 days by CFU-counting after serial dilution plating. After 2 h and 2 days infection we found significantly higher numbers of viable MAPΔfurA cells compared to MAPwt (Figure ). Complementation of MAPΔfurA almost completely restored the phenotype to wt level. Confocal analyses of macrophages infected with the different MAP strains for 30 min or 2 h after inside-outside staining of bacteria showed no differences in mycobacterial uptake by and adhesion to the macrophages (Figure ). These findings clearly exclude a better phagocytosis of MAPΔfurA which could account for the initial higher survival of MAPΔfurA in the macrophages. Interestingly, the percentage of viable MAPΔfurA cells decreased over time from ∼60% (of the initial inoculum) at 2 h to ∼40% at day 2 and to ∼18% at day 7. In contrast, no change in survival was visible in the wt and complement strain, indicating that the mutant was also being killed (Figure ). Calculating the survival rates at day 2 and 7 relative to the number of bacteria taken up by the macrophages after the 2 h infection period revealed that MAPΔfurA even though being more resistant to initial killing, is much more efficiently killed during the time course of macrophage infection (data not shown). Overall, these data indicate that MAPΔfurA can better resist the initial killing by the macrophage most probably caused by oxidative stress but is susceptible to macrophage defense mechanisms present at later time points of infection or less adapted to the phagosomal life style.

Induction of oxidative burst in macrophages and survival of MAP strains in macrophages. J774A.1 macrophages (Mø) were grown in DMEM and infected with an MOI of 10 with MAPwt, MAPΔfurA and the complemented strain MAPΔfurAC grown in MB-complete broth to an OD600 of 1.0. Mø were infected as described above for 1.5 h or treated with 100 μM menadione (+). CellROX Deep Red was then added to a final concentration of 2.5 mM for 30 min. (A) Histograms showing the production of ROS measured by CellROX Deep Red staining of unstimulated cells (green histogram) or infected Mø with MAPwt (dark blue histogram), MAPΔfurA (light blue histogram), and the complemented strain MAPΔfurAC (turquoise histogram). Mø treated with 100 μM menadione were used as a positive control (red histogram). Unstained cells are shown in the gray histogram. One representative experiment out of four is shown. (B) Fluorescence of untreated (-, green bar), menadione treated (+, red bar), and with MAPwt (dark blue bar), MAPΔfurA (light blue bar) or MAPΔfurAC (turquoise bar) infected macrophages was determined by FACS analysis. Shown are the results of four independent experiments (mean fluorescence intensity, MFI ± SEM). (C) Mø were infected as described above and incubated for 2 h, 2 and 7 days. Cells were disrupted and the lysate was diluted prior plating. 25 μl of the dilutions were plated in duplicate on MB-complete agar and incubated for up to 10 weeks. Five independent experiments were conducted, CFU of lysates and of the respective inocula used for macrophage infection were determined, and mycobacterial survival was expressed as percent to the inoculum. Statistical analysis was performed using one-way ANOVA (Kruskal–Wallis) with **p < 0.005. (D) Confocal maximum intensity projections of a double immunofluorescence staining of J774A.1 macrophages infected for 30 min or 2 h with MAPwt, MAPΔfurA, and MAPΔfurAC. Extracellular mycobacteria were labeled with anti-MAP-HBHA and green fluorescent Alexa Fluor® 488, intracellular bacteria with anti-MAP-HBHA and orange fluorescent Alexa Fluor® 568. Nuclei were stained blue with DAPI. Extracellular bacteria appear yellow and intracellular appear red. Bar represents 15 μm.

Hence, we analyzed the role of FurA for survival in the host. For this, C57BL/6 mice were infected intraperitoneally with 1 × 108 bacteria of exponentially grown cultures (OD600 of 1.0) and were sacrificed 4 weeks post infection. As shown in Figure , no significant differences in weight of liver and spleen, as well as number of bacteria or granuloma in the liver upon infection with MAPwt and MAPΔfurA were detected. However, as a trend we observed lower CFU numbers in the livers of MAPΔfurA infected cells. Overall, the lack of FurA in MAP did not alter survival rates compared to the wt and suggests that FurA might be dispensable for MAP-infection in the mouse model.

Biological fitness of MAPwt and MAPΔ Eight weeks old female C57BL/6 mice were infected with MAPwt (circles) and MAPΔfurA (triangle), grown to an OD600 of 1.0, with an infection dose of 1 × 108 cells and sacrificed 4 weeks later. The control group was treated with DPBS buffer (squares). Weight of spleen (A), liver (B), CFU detection of bacteria from liver (C), and the amount of granuloma in the liver (D) were chosen as read-out parameters. The results represent the mean ± SEM of one animal experiment with 9 mice in each group.

DISCUSSION

Bacterial iron homeostasis is tightly regulated and closely linked to the cellular response against cytotoxic effects due to the generation of ROS by free iron ions via the Fenton reaction (Cornelis et al., 2011). Accordingly, FurA was first identified in E. coli as global transcriptional regulator influencing expression of genes associated with iron acquisition and storage, but was later also shown to be involved in the regulation of genes activated by cellular stress (Dubrac and Touati, 2000). In mycobacteria, furA is located immediately upstream of katG encoding a catalase-peroxidase. Both genes are co-induced by peroxide treatment (Master et al., 2001; Milano et al., 2001; Jang et al., 2009). From this and the higher resistance of a furA mutant of the fast growing M. smegmatis to peroxide stress (Zahrt et al., 2001), it was proposed that in mycobacterial species, FurA is involved in the oxidative stress response regulation (Lucarelli et al., 2008). However, the precise role of FurA in pathogenic mycobacteria has not yet been elucidated.

In the present study, we characterized the role of furA in MAP. By RNAseq analyses, we were able to dissect a specific involvement of FurA in the regulation of the MAP response to reactive nitrogen and oxygen radicals, in intracellular survival, and in the expression of predicted virulence-associated genes. Most interestingly, unlike in many other bacterial species, where Fur regulates the expression of iron uptake and storage systems (Bagg and Neilands, 1987; Bsat et al., 1998; Abdul-Tehrani et al., 1999), FurA did not participate in the regulation of iron homeostasis in MAP. The involvement of FurA in modulating stress response and mechanisms for evading host cell defense as well as for invasion and infection is emphasized as 30 of the 48 regulated genes have been shown by works of other groups to be differentially expressed in MAP after infection of macrophages or in vivo infection of cattle (9 M/O, 21 in vivo; Wu et al., 2007; Cossu et al., 2012; Ghosh et al., 2013) and 26 of these genes have been found to be associated to stress response, as they were significantly regulated when MAP was exposed to acidic pH, nitrosamine, hydrogen peroxide, or heat (Wu et al., 2007; Kawaji et al., 2010; Cossu et al., 2012).

Our transcriptomic data also revealed that the FurA regulon includes genes with repressed and with activated transcription. Expression of a defined group of 13 genes encoding proteins of general stress response and enzymes necessary for the generation of reductive activity to detoxify oxygen radicals was higher in MAPΔfurA and is therefore repressed by FurA. For instance map1589c (ahpC) and map1588c (ahpD) were found to be induced in MAP upon oxidative stress (Kawaji et al., 2010) and during macrophage infection (Ghosh et al., 2013). Additionally, we found a peroxide dependent regulation of these genes by FurA in MAP. The alkyl hydroperoxide reductases AhpC and AhpD of the peroxireductase family (Prx) form an antioxidant system that is necessary for the reduction and neutralization of a wide range of organic hydroperoxides and peroxynitrite (Wood et al., 2003) to protect structures particularly sensitive to peroxide-mediated damage, such as lipids and nucleic acids (Halliwell and Gutteridge, 1984). A regulation of ahp genes by the Fur-like regulator PerR has been described for other bacteria, e.g. Bacillus subtilis and Campylobacter jejuni (van Vliet et al., 1999; Helmann et al., 2003). The regulatory mechanism bases on the property of mononuclear iron proteins such as FurA to lose repressive activity during both O2- and H2O2 stress by oxidation of Fe2+ to Fe3+. In addition to the missing activation of iron acquisition systems it supports our hypothesis of the role of FurA as ROS sensing regulator in MAP. CysH has been shown to contribute to protection against oxidative and nitrosative stress in Mtb (Senaratne et al., 2006) and many other FurA-repressed genes encode for enzymes involved in reactions liberating reducing activity in form of H+ ions.

Beyond this, FurA actively contributes to the expression of a subset of 35 genes many of which are suggested to be associated to virulence or intracellular survival. A major group of FurA-activated genes comprises genes encoding invasins (MAP1204/1204), antigens (MAP2168c, Webb et al., 1998), proteins of the Pro-Pro-Glu (PPE) family (MAP1003c, Akhter et al., 2012), PE-PGRS glycine-rich proteins, (MAP3868, Iantomasi et al., 2012), nlpc p60 (MAP1272c, Parthasarathy et al., 2012), and a mce family protein (MAP2114c, Gioffre et al., 2005) as well as several enzymes for cell entry, and it has been demonstrated for other bacteria, that homologs to these genes are important virulence factors (Kulasakara et al., 2006; Frederiksen et al., 2013). Interestingly, also some proteins involved in electron transfer were affected which suggest that FurA contributes to the maintenance of the canonical endogenous redox balance. Finally we found, that katG expression was not affected by furA deletion. Despite the fact, that a co-transcription of both genes and regulation of katG by FurA was suggested for other mycobacteria species (Master et al., 2001; Pym et al., 2001; Zahrt et al., 2001), we could not confirm this for MAP. Hence, katG might be regulated independently of FurA as proposed for Mycobacterium sp. strain JC1 DSM 3803 (Lee et al., 2009) or other bacteria (Italiani et al., 2011).

The particular role of FurA for virulence in mycobacteria has not been shown yet. Like other pathogenic mycobacteria, MAP is able to infect and survive in macrophages. Therefore, MAP must be able to react to the antimicrobial macrophage response such as reactive oxygen radicals. Indeed, our macrophage infection experiments with MAPwt, MAPΔfurA, and MAPΔfurAC showed the induction of ROS formation in the host cell. When monitoring the intracellular fate of MAPwt, MAPΔfurA, and MAPΔfurAC, we observed a significantly higher percentage of viable MAPΔfurA cells compared to the wt and MAPΔfurAC during the initial phase of infection. The higher percentage resulted from higher survival most probably given by an enhanced resistance of MAPΔfurA to oxidative stress mediated by the expression of FurA-repressed genes. Enhanced resistance to oxidative stress was also described for C. jejuni and M. smegmatis: in these species, deletion of perR and furA, respectively, resulted in a hyper resistancy to peroxide stress (van Vliet et al., 1999; Pym et al., 2001). In contrast, the survival rate of MAPΔfurA was considerably lower compared to the wt and MAPΔfurAC, indicated by killing of the mutant over time. This suggests that genes activated by FurA are needed for the adaptation of MAP to the macrophage phagosomal microenvironment.

In the intestinal tissue of the infected host, metabolic adaptation, and oxidative stress defense dominate the metabolism of MAP (Weigoldt et al., 2011, 2013). In contrast to the macrophage infections, the furA deletion seems not to affect MAP survival after 4 weeks intraperitoneal MAP infection of mice. Histological staining did not show significant differences in the number of granuloma and granuloma size as well as composition (data not shown). This suggests that the phenotype of MAPΔfurA is not advantageous when MAP has to counteract a complex immune response or has to metabolically adapt to the host environment. Nevertheless, we cannot exclude that FurA might play a role in a later stage of infection or effects of the FurA deletion might only be notable very early after infection.

CONCLUSION

In conclusion, in the presented study, we determined for the first time the regulon of FurA in a pathogenic mycobacterial species. Our data provide novel insights into the function of FurA in mycobacteria. They suggest two different physiological roles executed by FurA. Firstly, FurA acts as a repressor for a selective group of genes involved in the response to oxidative stresses. During oxidative stress it might operate as an iron-based sensor of ROS and, therefore, might be functionally closely related to PerR. Secondly, FurA contributes to activation of a second group of genes, many of which are suggested to contribute to MAP long term survival in macrophages. To conclusively investigate the role of FurA in virulence, however, infection experiments in animals developing clinical disease are required.

AUTHOR CONTRIBUTIONS

Ralph Goethe, Gerald-F. Gerlach, and Jochen Meens designed the experiments; Thorsten Meißner constructed the mutant, Thorsten Meißner and Kristin Laarmann accomplished mutant complementation, Michael Jarek performed RNA sequencing; Thorsten Meißner, Elke Eckelt, and Kristin Laarmann characterized the mutant and the complemented strain; Elke Eckelt performed stress experiments, Andreas Nerlich and Kristin Laarmann performed the macrophage experiments, Ralph Goethe and Siegfried Weiss designed the mouse infection experiments, Thorsten Meißner performed the mouse infections. Thorsten Meißner, Elke Eckelt, Jochen Meens, Ralph Goethe analyzed data; and Elke Eckelt, Ralph Goethe, Gerald-F. Gerlach, and Jochen Meens wrote the paper.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Table 1

Genes higher expressed in MAPΔfurA compared to MAPwt.

		Orthologous genes (% similarity)
RCNa	Annotation	Mtb	MAA	q-Valueb	Fold changec	Putative functiond	COGe
–	MAP0771	–	MAV_0960 (93.8)	<0.0001	8.23	Membrane protein	–
–	MAP0772	–	MAV_0961 (100)	<0.0001	7.00	Membrane protein	–
–	MAP0773	Rv2193 (45.88)	MAV_0962 (93.22)	<0.0001	11.63	Cytochrome c oxidase subunit iii	C
–	MAP1587c	Rv2471 (45.07)	MAV_2842 (91.81)	<0.0001	21.34	Alpha-amylase	G
ahpD	MAP1588c	Rv2429 (84.57)	MAV_2840 (99.43)	<0.0001	21.70	Alkyl hydroperoxide reductase	S
ahpC	MAP1589c*	Rv2428 (93.33)	MAV_2839 (100)	<0.0001	29.87	Alkyl hydroperoxide reductase c protein ahpc	O
–	MAP1741c	Rv2026c (68.58)	MAV_2507 (94.21)	0.00145	4.11	Universal stress protein	T
–	MAP1742c	Rv2026c (84.49)	MAV_2506 (74.4)	<0.0001	8.33	Universal stress protein	T
–	MAP1743c	Rv2032 (64.64)	MAV_2505 (90.09)	<0.0001	24.33	NAD(P)H nitroreductase	–
cysH_1	MAP2036	Rv2392 (83.85)	MAV_2153 (98.78)	0.00058	6.43	Phosphoadenosine phosphosulfate reductase	EH
–	MAP2037	Rv2393 (73.24)	MAV_2152 (91.19)	<0.0001	6.74	Cobalamin (vitamin B12) biosynthesis CbiX protein	S
fdxC_1	MAP2039	Rv1177 (89.62)	MAV_2150 (99.07)	<0.0001	7.72	Ferredoxin	C
–	MAP3677	Rv2242 (51.28)	MAV_4933 (81.93)	0.00069	4.75	PucR family transcriptional regulator	TQ

Table 2

Genes lower expressed in MAPΔfurA compared to MAPwt.

		Orthologous genes (% similarity)
RCNa	Annotation	Mtb	MAA	q-Valueb	Fold changec	Putative functiond	COGe
–	MAP0047c	Rv0040c (90.9)	MAV_0054 (73.79)	<0.0001	-12.91	Proline-rich 28 kDa antigen	–
–	MAP0079	Rv2002 (50.81)	MAV_0087 (95.12)	<0.0001	-8.00	Oxidoreductase	IQR
–	MAP0081*	Rv0063 (76.76)	MAV_0089 (85.61)	<0.0001	-7.67	Oxidoreductase	C
–	MAP0130	–	MAV_0124 (91.48)	<0.0001	-8.37	Hypothetical protein MAP0130	T
–	MAP0337	Rv1648 (52.83)	MAV_0364 (89.63)	<0.0001	-6.00	Membrane protein	–
–	MAP0847*	Rv1754c (83.96)	MAV_1035 (100)	<0.0001	-19.77	Conserved protein	–
PE	MAP1003c	Rv1040c (75.55)	MAV_1179 (54.03)	<0.0001	-12.40	PPE family protein	N
–	MAP1203	Rv1477 (76.48)	MAV_3301 (91.13)	<0.0001	-11.63	Invasion protein	M
–	MAP1204	Rv1478 (72.76)	MAV_3300 (86.88)	<0.0001	-9.58	Invasin 1	M
–	MAP1272c	Rv1566c (83.58)	MAV_3208 (76.04)	0.00059	-6.19	nlpc p60 family protein	M
–	MAP1388	–	MAV_3090 (74.19)	<0.0001	-4.67	Hypothetical protein MAP1388	–
–	MAP1418c	Rv3821 (56.25)	MAV_3059 (84.67)	0.00950	-4.00	Hypothetical protein MAP1418c	–
–	MAP1570	Rv1861 (50.49)	MAV_2858 (65.16)	<0.0001	-5.52	Membrane protein	–
–	MAP1639c	Rv0854 (75.86)	MAV_2785 (99.32)	0.00074	-4.21	Cyclase/dehydratase	I
furA	MAP1669c*	Rv1909c (88.02)	MAV_2752 (100)	<0.0001	-52.80	FurA	P
–	MAP1706	Rv1987 (73.1)	MAV_2710 (100)	<0.0001	-10.38	Chitinase	G
–	MAP1967c	Rv2223c (78.81)	MAV_2243 (88.43)	0.00014	-5.86	Exported protease	–
–	MAP2114c	Rv1968 (61.83)	MAV_2063 (90.27)	<0.0001	-6.67	Virulence factor mce family protein	Q
–	MAP2117c	Rv3501c (66.0)	MAV_0946 (66.53)	0.00103	-4.69	ABC transporter permease	Q
–	MAP2168c	Rv2376c (55.82)	MAV_2017 (90.35)	<0.0001	-11.59	Low molecular weight antigen mtb12	–
fdxC_2	MAP2607c	Rv1177 (96.22)	MAV_1316 (100)	<0.0001	-11.90	Ferredoxin	C
–	MAP2642	Rv1660 (64.75)	MAV_1280 (96.27)	<0.0001	-11.50	Polyketide synthase	Q
–	MAP3123c	–	MAV_3943 (74.1)	<0.0001	-5.20	Hypothetical protein	–
–	MAP3310	Rv3209 (80.68)	MAV_4157 (100)	0.00471	-4.73	Membrane protein	–
–	MAP3421c*	–	MAV_4276 (99)	0.00114	-4.00	Diguanylate cyclase phosphodiesterase	–
–	MAP3796	Rv0309 (78.99)	MAV_4853 (96.95)	<0.0001	-9.24	Conserved exported protein	S
–	MAP3800	Rv0312 (68.43)	MAV_2177 (52.59)	<0.0001	-10.89	Hypothetical protein	O
–	MAP3812c	Rv0116c (80.39)	MAV_4834 (100)	0.00030	-4.60	Hypothetical protein MAP3812c	S
–	MAP3868	Rv0754 (64.7)	MAV_4778 (91.69)	0.00011	-4.40	PE-PGRS family protein	G
–	MAP4204	Rv1888c (68.75)	MAV_4428 (88.82)	0.00382	-3.78	Transmembrane protein	–
–	MAP4206c	–	MAV_4427 (75.39)	<0.0001	-24.00	Efflux ABC transporter, permease protein	M
–	MAP4208	Rv3003c (40.8)	MAV_4424 (88.85)	<0.0001	-6.71	Acetolactate synthase	EH
–	MAP4209	–	MAV_4423 (100)	<0.0001	-6.73	Hypothetical protein MAP4209	–
–	MAP4210	–	MAV_4422 (99.37)	<0.0001	-10.00	3-oxoacyl-acp synthase	I
–	MAP4211	Rv3772 (46.15)	MAV_4421 (89.77)	<0.0001	-8.75	Histidinol-phosphate aminotransferase	E