{'fraction_non_alphanumeric': 0.03095630484165516, 'fraction_numerical': 0.0232951761614182, 'mean_word_length': 4.883647798742138, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Rhodococcus equi’s Extreme Resistance to Hydrogen Peroxide Is Mainly Conferred by One of Its Four Catalase Genes

Rhodococcus equi is one of the most widespread causes of disease in foals aged from 1 to 6 months. R. equi possesses antioxidant defense mechanisms to protect it from reactive oxygen metabolites such as hydrogen peroxide (H2O2) generated during the respiratory burst of phagocytic cells. These defense mechanisms include enzymes such as catalase, which detoxify hydrogen peroxide. Recently, an analysis of the R. equi 103 genome sequence revealed the presence of four potential catalase genes. We first constructed ΔkatA-, ΔkatB-, ΔkatC-and ΔkatD -deficient mutants to study the ability of R. equi to survive exposure to H2O2 in vitro and within mouse peritoneal macrophages. Results showed that ΔkatA and, to a lesser extent ΔkatC, were affected by 80 mM H2O2. Moreover, katA deletion seems to significantly affect the ability of R. equi to survive within murine macrophages. We finally investigated the expression of the four catalases in response to H2O2 assays with a real time PCR technique. Results showed that katA is overexpressed 367.9 times (±122.6) in response to exposure to 50 mM of H2O2 added in the stationary phase, and 3.11 times (±0.59) when treatment was administered in the exponential phase. In untreated bacteria, katB, katC and katD were overexpressed from 4.3 to 17.5 times in the stationary compared to the exponential phase. Taken together, our results show that KatA is the major catalase involved in the extreme H2O2 resistance capability of R. equi.


Introduction
Rhodococcus equi is a Gram-positive, facultative intracellular pathogen affecting foals up to six months of age. The major manifestation of R. equi infection is chronic suppurative bronchopneumonia, but 50% of contaminated foals also show clinical signs of intestinal disease and occasionally septic arthritis or osteomyelitis [1]. It is now well known that, in addition to being a veterinary concern, R. equi is an important opportunistic zoonotic pathogen mainly identified in immunocompromised people, and particularly AIDS patients [2,3]. Interestingly, different cases of infection by R. equi have also been reported in immunocompetent patients [4,5].
R. equi pathogenicity is closely linked to the bacterium's ability to survive and replicate in the lung macrophages of infected foals [6], where bacteria are generally killed by the combined action of a low pH (4.0-5.0), hydrolytic enzymes (including proteases, lipases, DNAses and RNAses) and the production of reactive oxygen metabolites such as hydrogen peroxide, H 2 O 2 [7]. R. equi has been characterized as being able to survive and multiply within a membrane-bound vacuole inside macrophages. This capability is conferred by the presence of a 80-90 kb virulence plasmid which inhibits i) the maturation of phagosomes, leading to the absence of phagosome acidification and prevents ii) the fusion of phagosomes with the lysosomes containing many hydrolytic enzymes [7]. In addition, R. equi must be able to mount an antioxidant defense system to ensure its persistence and multiplication within the highly oxidative environment of host macrophages generated notably by the reactive oxygen species (ROS) producing phagocyte oxidase Nox2 [8]. It has already been shown that R. equi is highly resistant to the action of H 2 O 2 [9], though plasmid-encoded proteins do not play a role in R. equi's resistance to H 2 O 2 [9]. An analysis of the R. equi 103 genome revealed the presence of numerous genes involved in oxidative stress resistance. Four genes in particular display homology with catalases [10], which are enzymes responsible for detoxifying H 2 O 2 [11]. This was the background for our investigation into the causes of the extreme resistance of R. equi to oxidative stress. Our study therefore focused on an analysis of the four potential catalase genes katA, katB, katC and katD. We studied the role of these proteins in H 2 O 2 resistance and in the ability of R. equi to survive within mouse peritoneal macrophages. These results were confirmed with the expression analysis of the catalase genes in response to H 2 O 2 exposure.

Identification of Four Genes Potentially Involved in H 2 O 2 Resistance
The in silico analysis of the genome sequence of the R. equi 103 strain [10] revealed the presence of four potential catalase genes, designated katA, katB, katC and katD. These four genes are at separate locations on the chromosome and do not seem to share genetic links. Sequence analysis of these proteins shows that KatA, KatB and KatC share the conserved catalytic residues and metal ligand amino acids of members of monofunctional heme catalases [12]. KatD shares the conserved catalytic residues and metal ligand amino acids of members of the nonheme catalases. These nonheme catalases utilize manganese ions instead of ferric heme in their active site and are therefore also known as manganese catalases [13]. No catalase of the bifunctional catalase-peroxidase family was identified. Phylogenetic relationships between the four catalases are represented in Figure 1, which shows that the three heme and manganese catalases constitute a distinct group, and that KatA and KatB are closely related. No signal sequence and/ or transmembrane domain were identified for these four catalases by TMHMM [14] or SignalP [15] analysis.
In order to estimate the involvement of each catalase in R. equi's oxidative stress resistance, we constructed corresponding mutants. To define the physiological roles of the four catalases, we studied the growth characteristics of the four deficient mutants. The characteristics were the same as those of the wild type (WT) strain under aerobic growth in BHI medium (data not shown). This similarity demonstrates that the catalases are not essential to cell viability under aerobic conditions.

DkatA and DkatC are the most Sensitive to H 2 O 2 Treatment
To evaluate the effect of exogenous H 2 O 2 stress on the survival of each mutant, cells were challenged with 80 mM of H 2 O 2 under agitation, and CFUs determined 30 min after treatment ( Figure 2). DkatA was the most sensitive mutant to H 2 O 2 since it did not survive when treated in the stationary phase, and its survival decreased by two orders of magnitude compared to the WT strain when treated in the exponential phase. DkatC survived as well as the WT strain in the stationary phase but in contrast, when the cells were in the exponential phase, DkatC was the most susceptible to H 2 O 2 exposure with a decrease in survival of over three magnitudes compared to the WT strain. Both DkatB and DkatD had the same resistance capability as the WT strain, whether in exponential or stationary phases. The drop in survival rates for DkatA and DkatC after exposure to H 2 O 2 demonstrate that, in these experimental conditions, only KatA and KatC primarily contribute to the oxidative defense system of R. equi.

DkatA is the most Susceptible to Elimination by Macrophages
During phagocytosis, phagocytic cells generate superoxide and other ROS involved in antibacterial activity [16]. This study compared the intracellular survival of the WT and DkatA, DkatB, DkatC and DkatD strains inside infected mouse peritoneal macrophages ( Figure 3

katA is Overexpressed in Response to H 2 O 2 Treatment
We evaluated the effect of a sublethal H 2 O 2 concentration (50 mM) defined previously (data not shown) on the expression of katA, katB, katC and katD transcripts during exponential ( Figure 4A) and stationary phases ( Figure 4B). The expression of catalases was measured at several timepoints: 5, 10 and 20 minutes after treatment in the exponential phase and 5, 10, 20, 30 and 60 minutes after treatment in the stationary phase. Results showed that only the katA gene was overexpressed in response to H 2 O 2 treatment. For example, when H 2 O 2 was added in the exponential phase, katA was overexpressed 3.11 (60.59) times in bacteria treated for 5 min compared to untreated bacteria ( Figure 4A). However, when H 2 O 2 was added in the stationary phase, katA was overexpressed 367.90 (6122.63) times in bacteria treated for 10 min compared to untreated bacteria ( Figure 4B). We did not observe any modifications in the rate of transcripts for katB, katC or katD ( Figures 4A and 4B).

katB, katC and katD are Overexpressed in the Stationary Growth Phase Compared to the Exponential Phase
In untreated bacteria, we observed overexpressions of katB, katC and katD in the stationary phase compared to the exponential phase. Irrespective of exposure to oxidative stress, katB was overexpressed 17.55 (63.22) times in the stationary phase compared to the exponential phase ( Figure 5). katC and katD were overexpressed 4.51 (60.59) and 4.31 (60.50) times in the stationary phase compared to the exponential phase, respectively ( Figure 5).

Discussion
Resistance to oxidative stress is one of the key processes that allow pathogens to survive within macrophages, which is why pathogenic bacteria have developed different antioxidative defenses to eliminate ROS or repair their damage [17]. In silico analysis shows that R. equi 103 contains three monofunctional heme catalases (KatA, KatB and KatC) and a monofunctional manganese catalase (KatD). The present study addressed the issue of the involvement of these potential H 2 O 2 -degrading catalase enzymes in R. equi virulence and demonstrates that survival following exposure to exogenous H 2 O 2 and the efficient proliferation of this bacterium in macrophages mainly rely on the activities of KatA, one of its four catalases.
Benoit et al. showed that the extreme resistance of R. equi to exposure to a high H 2 O 2 concentration (100 mM) did not depend on plasmid-encoded proteins [9]. This was confirmed by the analysis of virulence plasmid-chromosome regulatory crosstalk by Letek et al. [10] which shows that the presence of the virulence plasmid did not interfere with the expression of catalase genes. These results suggest that chromosomal genes are key to this function, and genomic analysis of the R. equi chromosome identified an unusually high number of catalases: three monofunctional heme catalases (katA, katB and katC) and one manganese catalase (katD). Surprisingly, no bifunctional catalase-peroxidase was identified although for numerous pathogens, these enzymes have already been shown to be involved in the infection of macrophages [18][19][20].
An analysis of the physiological functions of the four catalases following exposure to exogenously-provided H 2 O 2 showed that only KatA and KatC are part of the peroxide defense system of R. equi. KatA seemed to be more important than KatC under these stress conditions. The fact that katB and katD mutations did not affect survival following H 2 O 2 exposure suggests that either i) these catalases are non functional, ii) they function as a surrogate in the absence of the primary antioxidant, but their role can only be detected in a double or triple catalase mutant [21] or iii) they have specialized roles for protecting cells against oxidative stress perhaps due to different cellular locations. Although no secretion peptide signal was identified for R. equi catalases, it has already been shown that a catalase from B. subtilis could be secreted independently of the Sec-secretion apparatus [22,23].
Inactivation of catalase-encoding genes had little if any effect on pathogenicity-associated traits of numerous bacteria including Salmonella typhimurium [24], Yersinia pestis [25], Staphylococcus aureus [26], and Neisseria gonorrhoeae [27]. This research demonstrates that significantly lower quantities of viable intracellular bacteria were found in the infected murine macrophages of DkatA-deficient mutants during the 72-h infection period than in the WT strain. This suggests that the katA gene is linked to R. equi virulence mechanisms. This observation could be linked to the oxidative microenvironment present within phagosomes, where ROS and reactive nitrogen species are synthesized in huge quantities following macrophage phagocytosis [17,28].
Studies of physiological functions are complex when multiple catalases are present in a single bacterial cell and differences in regulation are a critical factor in understanding the in vivo roles of these enzymes. Our results show that only katA expression is induced in response to exogenous H 2 O 2 . This fast and transitory overexpression is induced mainly in the stationary phase (overexpression factor of 367.96122.6), which is in accordance with the critical role of katA in stationary-phase survival. In  contrast, despite its role in the exponential phase, katC is not induced in the stationary phase. This phenomenon, already observed for a catalase-peroxidase of Legionella pneumophila [11], is not yet fully understood. In untreated bacteria, katB, katC and katD expressions increase in the stationary growth phase. These overexpressions could prepare the cells for subsequent stresses [29]. In this context, we could hypothesize that the lack of effect of the mutation of these genes in the stationary phase is due to a compensatory catalase mechanism between them or that these catalases play a role in protecting against the oxidative stress generated by growth in a minimal medium as already described [30]. Our results therefore suggest that katA is regulated by a specific stress response transcriptional regulator able to sense the presence of H 2 O 2 , whereas the other three catalases are under the control of a RpoS-like central regulator of gene expression in the stationary phase [31]. Moreover, other unidentified conditions could induce katB, katC and katD expression. We could suggest, for example, that the expression of manganese catalase katD is induced by iron starvation as previously described for the manganese catalase of Pseudomonas aeruginosa [32].

Conclusions
This study addressed two fundamental questions: how are the four R. equi catalases involved in survival i) following external exposure to H 2 O 2 and ii) in mouse peritoneal macrophages? Taken together, our results show that KatA is the major determinant of resistance to exogenous H 2 O 2 exposure in the  . In untreated bacteria, katB, katC and katD are overexpressed in the stationary growth phase. R. equi WT was grown in BHI broth without any H 2 O 2 treatment and cells were collected at the exponential (OD of 0.4) and stationary (OD of 1) growth phases for total RNA extraction. cDNAs derived from total RNA, were used for real time PCR. The overexpression factor was calculated using the 2 2DDCt method [36]. The data are the mean of overexpression factor in stationary phase collected bacteria compared to exponential phase collected bacteria 6 standard deviation of triplicate measurements from four reverse transcriptions of two independent experiments. doi:10.1371/journal.pone.0042396.g005 stationary phase and, to a lesser extent, in the exponential phase. KatA is also crucial to R. equi's survival in mouse peritoneal macrophages. Of the other three R. equi catalases, only KatC plays a role in resistance to external H 2 O 2 exposure, but only in the exponential phase.

Ethics Statement
The animal experiments (Survival assays in mouse peritoneal macrophages) were performed under a protocol approved by the Institutional Animal Use and Care Committee at Università Cattolica del S. Cuore, Rome, Italy (Permit number: N21, 12/05/ 2010) and authorized by the Italian Ministry of Health, according to Legislative Decree 116/92, which implemented the European Directive 86/609/EEC on laboratory animal protection in Italy. Animal welfare was routinely checked by veterinarians of the Service for Animal Welfare. Table 1 lists the bacterial strains and plasmids used in this study. The R. equi 103 strain used is devoid of its virulence plasmid. Wild type (WT) or mutant R. equi strains and the Escherichia coli strain used to construct mutants were routinely grown at 37uC in Brain Heart Infusion (BHI) broth and Luria-Bertani media respectively, with vigorous shaking (200 rpm). When required, antibiotics were added to culture media at the following concentrations: apramycin (30 mg/ml); ampicillin (100 mg/ml) [33]. Dilutions were performed in physiological solution (0.9% NaCl).

Construction of Four Catalase Mutants
A double crossover homologous recombination strategy based on suicide vectors derived from pUC19 containing a cassette consisting of the upstream homologous region, the apramycin resistance gene (aaC4) and the downstream homologous region as previously described [33] was used to inactivate the four catalase genes identified in R. equi: REQ4750, REQ44520, REQ26870 and REQ26290, respectively named DkatA, DkatB, DkatC and DkatD. Briefly, for DkatA construction, a DNA fragment consisting of the katA gene and <600 bp up-and downstream of the start and stop codons was amplified with primers L39 and L79 containing respectively EcoRI and XbaI restriction sites (Table S1). Then, the product was inserted into pUC19 using EcoRI and XbaI restriction sites. To delete the katA gene's internal region, an inverse PCR was performed on the resulting pUC19::katA plasmid using divergent primers L162 and L163 with NsiI restriction site sequence in order to obtain a linear pUC19-DkatA plasmid containing a 70% deleted version of katA. Separately, the aacC4 gene was amplified using primers L118 and L119 each containing an NsiI restriction site. Both pUC19-DkatA and aacC4 amplicons were digested by NsiI and ligated to obtain a pUC19-DkatA::aacC4 plasmid. This plasmid was introduced into R. equi by electroporation, and transformants were selected on BHI agar supplemented with 80 mg/ml apramycin. Allelic exchange double recombinants were selected as previously described [33]. The same procedure was applied for the construction of DkatB, DkatC and DkatD. Construction was verified by Southern blot (data not shown). The primers used are listed in Table S1.

Hydrogen Peroxide Challenge Assays
Ten ml of exponentially growing cultures of R. equi WT, DkatA, DkatB, DkatC or DkatD (Optical Density at 600 nm (OD 600 ) of 0.2) were harvested by centrifugation and pellets were resuspended in 10 ml of 0.9% (w/v) NaCl. OD 600 of stationary growing cultures were measured and the volume necessary for obtaining an OD 600 of 0.2 in a volume of 10 ml was harvested by centrifugation and resuspended in 10 ml of 0.9% (w/v) NaCl. Both cell suspensions were treated with 80 mM H 2 O 2 (Riedel de Haën) and bacterium

Survival Assays in Mouse Peritoneal Macrophages
The survival of R. equi WT, DkatA, DkatB, DkatC and DkatD in mouse peritoneal macrophages was tested by using an in vivo/ in vitro infection model as described previously [34,35]. Briefly, after 16 h of growth, bacterial cells were pelleted and resuspended in an adequate volume of phosphate-buffered saline (PBS) for injection. Male BALB/c mice (10 weeks old) were infected with 10 7 to 10 8 cells of each strain by intraperitoneal injection. After an 8 h infection period, peritoneal macrophages were collected by a peritoneal wash, centrifuged then suspended in Dulbecco's modified Eagle's medium containing 10 mM HEPES, 2 mM glutamine, 10% bovine fetal serum, and 1X nonessential amino acids supplemented with vancomycin (10 mg/ml) and gentamicin (150 mg/ml). The cell suspension was dispensed into 24-well tissue culture plates and incubated at 37uC under 5% CO 2 for 2 h. Duplicate wells of infected macrophages were lysed with detergent at 24, 48, and 72 h post infection, and lysates were plated on BHI agar to monitor bacterial survival. All experiments were performed in triplicate, and the results were subjected to statistical analysis using one-way analysis of variance with Bonferroni's correction post test using GraphPad Prism version 5.00 for Windows (GraphPad Software, San Diego, CA).

Catalase Overexpression Assays
R. equi WT in mid-exponential (OD 600 of 0.4) and stationary (OD 600 of 1) growth phases was treated with 50 mM H 2 O 2 . Previous results have shown that 50 mM H 2 O 2 is a sublethal concentration (data not shown). A control for each growth phase was performed without H 2 O 2 treatment. Bacterial cells were incubated (37uC, 200 rpm) for 5, 10, 20, 30 or 60 min after treatment and total bacterial RNA was isolated: cells were harvested by centrifugation, rinsed in PBS 1X and harvested again by centrifugation. Cells were resuspended in 50 mg/ml lysozyme (Sigma) prepared with Tris/HCl-Ethylenediaminetetraacetic acid (EDTA) buffer, and lysed with a Ribolyser (Hybaid). Samples were incubated for 90 min at 37uC. Total RNA was isolated using the RNeasy RNA mini kit (Qiagen) according to the manufacturer's instructions. A 15-min on-column DNA digestion was also performed. After elution, RNA was treated with Turbo DNA-free (Ambion) according to the manufacturer's instructions. The concentrations and quality of total RNA were evaluated using Experion (Biorad) with the RNA Std Sens kit, according to the manufacturer's instructions. One microgram of total RNA was reverse-transcribed using the Omniscript enzyme (Qiagen), recombinant RNAsin (Promega) and random hexamer primers (Invitrogen) according to the manufacturer's instructions.
The resulting cDNAs were used for subsequent PCR amplification with specific primers, designed using Primer3 software. Primer sequences are listed in Table S1 and were used with the QuantiFast SYBR Green PCR mix (Qiagen) according to the manufacturer's instructions, using a MastercyclerH ep realplex (Eppendorf). The quantification of 16S RNA levels was used as an internal control. Two independent bacterial cultures were used to obtain two different RNA extracts. Each RNA was reversetranscribed twice in order to obtain four cDNAs. Amplification, detection and real-time analysis were performed in triplicate with the four cDNAs. The overexpression factor (OF) in treated bacteria vs control, was calculated using the 2 2DDCt method [36]. To evaluate the efficiency of the amplification, a standard curve was constructed using the cycle threshold (C T ) versus 10-fold dilution (data not shown). The reaction's specificity is given by detection of the melting temperatures (Tm) of the amplification products via the melting curve.

Supporting Information
Table S1 Oligonucleotide primers used in this study. (DOCX)