Iron Starvation Conditions Upregulate Ehrlichia ruminantium Type IV Secretion System, tr1 Transcription Factor and map1 Genes Family through the Master Regulatory Protein ErxR

Introduction

Ehrlichia ruminantium, the causal agent of heartwater, a fatal disease of ruminants in sub-Saharan African and other tropical regions, belongs to the Anaplasmataceae family and is transmitted by ticks of the genus Amblyomma (Dumler et al., 2001; Allsopp, 2010). In the mammalian host, E. ruminantium mainly infects brain capillary endothelial cells and replicates inside membrane-bound vacuoles (Zweygarth and Josemans, 2001). E. ruminantium has a biphasic developmental cycle in which elementary bodies (EB), the infectious form of the bacterium, first adhere and enter host cells. After internalization, EB differentiate into reticulate bodies (RB), the vegetative and non-infectious form, which divide by binary fission. Within 4–5 days, RB reorganize into EB, which are released from the vacuole by the lysis of the host cell to initiate a new infectious cycle (Moumène and Meyer, 2015).

Intracellular pathogenic bacteria belonging to the Anaplasmataceae family such as Ehrlichia, use a dedicated system, the type IV secretion system (T4SS), to inject some bacterial proteins, named effectors, to evade the host's immune responses and to hijack host cell processes in order to survive and proliferate in a safe replicative niche (Moumène and Meyer, 2015). T4SS is well conserved in the Anaplasmataceae family and several T4SS effectors (T4Es) have been shown to be crucial for the pathogenicity of Anaplasma phagocytophilum, Anaplasma marginale, and Ehrlichia chaffeensis (as reviewed by Rikihisa, 2010). Little is known about T4SS in E. ruminantium and no T4Es have yet been characterized (Collins et al., 2005; Frutos et al., 2007).

In several bacteria, it has been shown that the expression of T4SS is tightly regulated by transcription factors (Li and Carlow, 2012; Martín-Martín et al., 2012). Cheng et al. (2008) showed that the five virB/D4 genetic loci of E. chaffeensis T4SS are co-regulated by the transcription factor EcxR to allow specific expression, depending on the developmental stage of the bacteria (Cheng et al., 2008). Previous studies have also demonstrated developmental regulation of the expression of T4SS components during the intracellular life cycle of A. phagocytophilum in human peripheral blood neutrophils (Niu et al., 2006). A ecxR ortholog, apxR, has been found in several A. phagocytophilum strains (Wang et al., 2007a,b). ApxR regulates the expression of the transcription factor tr1 (Wang et al., 2007b) and of the downstream p44E locus (Wang et al., 2007a). Interestingly, the virB/D4 loci of A. phagocytophilum, E. chaffensis, and E. canis are flanked by genes encoding outer membrane proteins (OMP) belonging to the p44/msp2 family, which are paralogs of Ehrlichia ruminantium map1 genes (Dunning Hotopp et al., 2006; Rikihisa, 2011). The exact environmental cues, which stimulate the expression of apxR and ecxR, are not known.

Considering all these observations, we analyzed the genome of E. ruminantium and found an ortholog of ecxR and apxR in E. ruminantium, hereafter termed erxR, for the Ehrlichia ruminantium expression regulator. Moreover, as we observed in A. phagocytophilum, tr1 is present in E. ruminantium upstream of the cluster of map1 OMPs (Wang et al., 2007b). In E. ruminantium, the 16 paralogs of the map1 multigene family are expressed in bovine endothelial cells and some are preferentially transcribed in the tick or in the mammalian host (van Heerden et al., 2004). However, whether or not ErxR is a regulatory protein, which drives the expression of tr1 gene, the map1 genes family, and the virB/D4 loci in E. ruminantium and the triggering stimuli are currently unknown.

Microorganisms have evolved sensory mechanisms to regulate their cellular activities in response to environmental changes. This is particularly true for bacterial pathogens whose expression of virulence factors is tightly regulated in response to host and non-host environments (Hyytiäinen et al., 2003). Thus, the regulation of T4SS in response to host cues enables the efficient use of bacterial resources and facilitates colonization, leading to full infection (Abromaitis et al., 2013). One such environmental signal is iron, which is an essential cofactor in various enzymatic reactions like respiration, DNA replication, oxygen transport, response to oxidative stress, but can be toxic at excessive intracellular concentrations (Andrews et al., 2003). Therefore, iron scavenging from the limited sources of free iron available in the host is a crucial determinant of bacterial pathogenicity (Ratledge and Dover, 2000). Pathogens have evolved several ways to scavenge iron from the host, including the expression of iron acquisition genes in response to low iron concentrations (Brickman et al., 2011; Portier et al., 2014). In Gram negative bacteria, TonB-dependent outer membrane receptors (TBDR) are required to transfer iron chelates and heme into the periplasm under poor iron conditions, with subsequent transport to the cytoplasm (Shultis et al., 2006). Whether Map1 proteins play a similar role in E. ruminantium's iron uptake or sensing is currently not known.

In this study, we consequently investigated the effect of the newly identified ErxR protein and iron conditions on the regulation of the virB/D4, tr1 and map genes expression in E. ruminantium. For the first time in Ehrlichia, we identified (i) a master regulatory protein responsible for the coordinated regulation of the expression of virB/D4, tr1, and map genes, and (ii) one triggering environmental signal i.e., iron depletion, which activates this regulation cascade. Our work also demonstrates that map genes share at least one regulatory pathway with genes encoding the T4SS and may therefore be important pathogenicity determinants for iron acquisition or host cell infection.

Materials and Methods

Culture Conditions

E. ruminantium Gardel strain (passages 30 to 52) was routinely propagated in bovine aortic endothelial (BAE) cells as previously described (Marcelino et al., 2005). To evaluate the growth characteristics of E. ruminantium under iron depletion conditions, the strain was grown in TC25 cm² flasks in BHK-21 cell medium supplemented with 2 mM glutamine, 10% heat inactivated fotal bovine serum (FBS), penicillin (100 IU/ml), streptomycin (100 mg/ml). For iron response experiments, the medium was supplemented with iron (100 μM FeSO₄) or an iron chelator (100 μM 2,2′-bipyridyl; BPD) as described in Breuer et al. (1995) and Romeo et al. (2001). The cells were kept in a humidified atmosphere supplemented with 5% CO₂ at 37°C. FeSO₄ or 2,2′-bipyridyl was added when 80% cell lysis was observed, at 120 h post-inoculation (hpi). The E. ruminantium infected cell monolayer (1 TC25 cm² flask) was harvested by trypsinisation 24 h after chemical inoculation. From the 6 ml of infected cell supernatant, a 600 μl sample was collected by centrifugation at 20,000 × g for 10 min. The pellet was stored at −80°C until DNA extraction. The remaining 5,400 μl of infected cells were centrifuged at 20,000 × g for 10 min. The pellet was resuspended in TRIzol reagent (Invitrogen) and stored at −80°C until RNA extraction (Pruneau et al., 2012).

Quantitative Detection of E. ruminantium

Genomic DNA was extracted from the 600 μl samples described above using the QIAamp DNA Mini Kit (Qiagen, France). The number of bacteria per sample was quantified by q-PCR, targeting the single copy of map1 gene encoding a major antigenic protein. The primer sequences are shown in Supplementary Table 1 (Pruneau et al., 2012). A standard curve was established using gDNA of Gardel serially diluted from 7 × 10⁶ to 7 × 10¹ copies μL⁻1, to determine the number of bacteria per microliter (Pruneau et al., 2012). Four microliters were added to Taqman master mix (Applied Biosystems, France), following the manufacturer's instructions. PCR conditions were as follows: 2 min at 50°C, 10 min at 95°C, and 40 cycles with 15 s at 95°C and 1 min at 60°C.

Relative Gene Expression: RNA Preparation and qRT-PCR

Total RNA was extracted using TRIzol reagent. RNA pellets were dissolved in 100 μl of DEPC water and treated with turboDNAse (Ambion, France). The purity and concentration of the isolated RNA were assessed using a NanoDrop 2000c (Thermo Scientific, France). RNA samples were diluted in RNase-free water to obtain a final concentration of 0.5 μg/μL. RNA samples were reverse-transcribed with the SuperScript VILO cDNA Synthesis Kit (Invitrogen, France), according to the manufacturer's instructions. Quantitative PCR was performed in a 7500 Real-Time PCR System (Applied Biosystems, France) using a Power SYBR Green PCR Master Mix (Applied Biosystems, France) and the primers listed in Supplementary Table 2. Reactions were performed in 25 μl volume with 5 ng template cDNA and 5 μM of each primer. The amplification conditions were as follows: 2 min at 50°C, 10 min at 95°C and 40 PCR cycles (30 s at 95°C and 1 min at 60°C). An additional dissociation step of 15 s at 95°C, 20 s at 60°C and 15 s at 95°C was added to assess non-specific amplification. A negative control without cDNA template was included for each primer combination. Amplifications were performed in technical replicates consisting of independent cDNA syntheses derived from the same RNA sample and in three independent biological replicates. The relative expression of erxR was calculated by dividing the number of transcripts by the total number of bacteria at each time point. Fold change was then calculated by comparing the relative expression at each time point and the relative expression at 96 hpi, the stationary phase. Ratios were calculated from the number of transcripts and normalized to recA as described in Gonzalez-Rizzo et al. (2006).

Identification of Orthologs of EcxR in the E. ruminantium Genome

To identify orthologs of E. chaffeensis ecxR in the genome of E. ruminantium, we used the same strategy as that described in Li and Carlow (2012). The protein sequence of ECH_0795 (YP_507593) was used as a query to search the genome of E. ruminantium. Multiple sequence alignment was performed using ClustalW (http://www.ebi.ac.uk/Tools/msa/clustalw2/) (Larkin et al., 2007). Sequence identity values between the two sequences were generated using BlastP. A Helix-Turn-Helix motif was determined using Pfam (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan) to find domains and motifs present in the protein sequence.

Cloning and Expression of erxR

Full-length erxR was PCR amplified using the primers listed in Supplementary Table 2, and ligated into the NdeI and XhoI sites of the pET29a(+) vector (Novagen). The resulting plasmid, defined herein as pErxR, was cloned into E. coli DH5α (Invitrogen) for amplification, purified using the QIAGEN Plasmid Maxi Kit (Qiagen, France), and cloned into E. coli BL21 (DE3) (Invitrogen) for protein expression. Protein expression was induced with 4 mM isopro-pyl-β-D-thiogalactopyranoside (IPTG) in 250 ml terrific broth. The protein then was purified using the Ni-NTA Fast Start Kit (Qiagen, France). ErxR expression was determined by Western blot analysis using anti-His tag antibody (Qiagen, France).

Construction of pUA66-Derived Promoter Plasmids

The pUA66 plasmid was used to analyse promoter activity based on the expression of a green fluorescent protein (GFP). The promoter regions were PCR amplified from the genomic DNA of Ehrlichia ruminantium Gardel strain, using the primers listed in Supplementary Table 2. Forward and reverse primers contained HindIII and BamHI restriction sites for cloning into the pET29a(+) plasmid. After cloning in pET29a(+), the promoters were digested with XhoI and BamHI for directed cloning into the pUA66 plasmid (Castaño-Cerezo et al., 2011). BL21 (DE3) cells were co-transformed with pErxR and each of the GFP reporter constructs, individually. The pET29a(+) vector alone was used as a negative control. Cotransformants were grown in LB medium supplemented with 50 μg/ml kanamycin at 37°C for 2 h, followed by induction with 1 mM IPTG for 4 h. Induced bacteria were visualized as described below.

Microscopy

A 6-μl drop of E. coli BL21 (DE3) co-transformed with the pUA66 promoter containing one of the different promoters and the pET29a-erxR plasmid (Supplementary Table 3), resuspended in LB growth medium, was spotted onto a Superfrost Plus slide (Fisher Scientific Ltd, UK) and visualized using a Nikon Eclipse 80i epifluorescence microscope (Nikon, France). Fluorescent images were acquired with a Nikon DXM1200F digital camera (Nikon, France), using Nikon ACT-1 software (Nikon, France). Fluorescence intensity was calculated by measuring the area, integrated intensity and mean gray value of the fluorescent bacteria and the background with ImageJ (National Institute of Health, USA). Corrected total cell fluorescence (CTCF) was calculated using the following formula: integrated density—(area of the cell × mean background readings). The average and statistical differences between the bacteria containing the plasmids with the different promoters and controls were calculated using the CTCF values of all the bacteria in four different fields of view. Images were processed to size, and brightness and contrast were adjusted after the measurements, using Adobe Photoshop cs5 (Adobe Systems Inc., California, USA).

Statistical Analyses

Statistical analyses used Student's t-test and a P < 0.05 was considered significant.

Results

Identification of One ecxR Ortholog in the E. ruminantium Genome

We used the ecxR sequence (YP_507593.1) from E. chaffeensis to search NCBI databases using the BLAST tool. With this approach, we identified ERGA_CDS_03000 (YP_196226.1) in the genome of E. ruminantium as the closest ortholog to ecxR. The results of the BLAST search revealed a putative conserved domain belonging to the HXT_XRE superfamily of DNA binding proteins (cl17200). We identified a helix-turn-helix structure by comparing it with the structure found in Wolbachia (Li and Carlow, 2012) and further confirmed by Pfam (https://www.ebi.ac.uk/Tools/hmmer/search/hmmscan). This protein comprising 124 amino acids has a predicted molecular mass of 14.25 kDa. Alignment of the deduced amino acid sequences of the various orthologs is shown in Figure 1A. By homology with ApxR and EcxR, we named this protein E. ruminantium expression regulator, ErxR. Comparison of the sequence identities of these proteins revealed a high degree of conservation (82% identity) between EcxR and ErxR. ApxR and ErxR showed 40% identity at the amino acid level. Structural analyses indicated that all orthologs shared a conserved helix-turn-helix domain that may function as a sequence specific DNA binding domain, such as in transcription regulators (Aravind et al., 2005).

FIGURE 1

Figure 1. Potential ortholog of the Ehrlichia type IV secretion system regulator EcxR in E. ruminantium and organization of virBD loci and map1 cluster. (A) Amino acid alignment of EcxR and its orthologs from E. ruminantium (ERGA_CDS_03000) and A. phagocytophilum (ApxR). Identical and conserved residues are indicated by stars and dots, respectively. The positions of amino acids predicted to form a helix structure and beta-strand are shown. (B) The genomic map of E. ruminantium is represented by a circle. The origin of replication (ori) is indicated by a black box. Gray arrows show the virBD gene loci and cluster map1. T4SS operon 1 consists in virB8a, virB9a, virB10, virB11, and virD4 and operon 2 consists in sodB, virB3, virB4a, virB6a, virB6b, virB6c, and virB6d. Four copies of virB2 were identified upstream of virB4b. map1 genes belong to a multigenic family organized in one cluster of 16 paralogs located downstream of the transcriptional factor tr1. The length of the arrows is proportional to the length of the gene. (C) Cluster structure of virBD gene loci and map1 genes. virBD and map1 genes whose expression under iron starvation condition was tested and are represented by gray arrows. The names and length of the genes are indicated above and below the arrows. The upstream regions amplified by PCR for regulation assays (hatched boxes) are indicated.

Architecture of T4SS and map1 Genes Clusters of E. ruminantium

We compared the genetic arrangement of E. ruminantium to that of E. chaffeensis. The five virBD loci are represented in Figures 1B,C. In E. ruminantium, the genome sequence revealed the presence of two major operons. virD4, virB11, virB10, virB9a, and virB8a were located in operon 1. Operon 2 was seen to be located in the negative strand and contained four copies of virB6, along with one copy of virB4a and virB3, all located downstream of sodB. Four duplicated versions of virB2 were also located upstream of virB4b, while copies of virB8 and virB9 (namely virBb8 and virB9b) were scattered along the genome. We also present the arrangement of the map1 gene family previously reported by Postigo et al. (2007) (Figure 1C).

Analysis of erxR Expression during the Life Cycle of E. ruminantium

To determine the relative expression of erxR throughout the developmental cycle of E. ruminantium, erxR mRNA expression was analyzed by qRT-PCR. Compared to 96 hpi, the expression of erxR decreased at 24, 48, and 72 hpi, and peaked at 120 hpi, having increased 4 fold (Figure 2A), which corresponds to the time of lysis (early time point in the following round of infection). Data from previous independent microarray experiments (Pruneau et al., 2012) also demonstrated that the maximal expression of erxR in the post-exponential growth phase is significant and reproducible (data not shown). These results suggest that the up-regulation expression of erxR correlates with early stages in the development cycle in vitro before the bacteria enter the host cell.

FIGURE 2

Figure 2. (A) erxR expression peaks at 120 hpi in E. ruminantium. Quantitative RT-PCR was used to determine the temporal expression of erxR. Relative expression at different developmental stages was normalized by dividing the transcripts number by the number of bacteria. Fold differences were evaluated by comparing each time point to 96 hpi (the stationary phase). Data were obtained from triplicate samples and are expressed as means + standard deviation. (B) ErxR activates the transcription of virBD, tr1 and map genes. Fluorescence intensity was used to measure the transcriptional activities of gfp reporter constructs. The values are means + standard deviations for three specimens, and measurements were taken from four different view fields. An asterisk indicates that values differ significantly (P < 0.001) from the controls. Western blot analyses were performed of samples from the fluorescence assays using an anti-His antibody to verify the expression of rErxR. Arrowheads indicate the position of rErxR.

rErxR Activates gfp Reporter Fusions

The promoter regions (hatched boxes in Figure 1) for virBD genes, map1, tr1, and erxR were cloned into a gfp reporter plasmid and transformed into E. coli BL21 (DE3) carrying pErxR or an empty pET29a(+) vector to investigate if ErxR regulatory protein drove their expression. The virB3-gfp reporter constructs presented a significant increase in fluorescence intensity after IPTG induction (270,000 units of fluorescence intensity) compared to samples lacking IPTG (10,000 units of fluorescence) or compared to the control (50,000 units of fluorescence) (Figure 2B). Similarly, the reporter construct showed an induction of virB2a by the recombinant ErxR protein. Activation was also observed for tr1 (200,000 units), erxR promoter (250,000), and map1 promoter (~45,000) (Figure 2B). Western blot experiments confirmed that the expression of ErxR (14 kDa band) was only detected following induction with IPTG (Figure 3).

FIGURE 3

Figure 3. Production of rErxR in E. coli. E. ruminantium erxR was cloned into the pET29a(+) vector, expressed, and purified using nickel chelate chromatography. The purified protein was subjected to SDS- PAGE analysis, followed by Coomassie blue staining (A) and Western blot analysis using an anti-His tag antibody (B). We used pre-stained protein size standards. Each lane contained 8 μg of recombinant protein.

The Expression of T4SS, erxR, tr1, and map Genes Is Induced by Iron Depletion in E. ruminantium

As iron uptake mechanisms are closely associated with bacterial pathogenesis and may be connected to the expression of certain virulence determinants in E. ruminantium, we investigated the expression of T4SS, tr1 and map genes clusters in response to iron starvation. We incubated the bacteria in media containing iron or an iron chelator (BPD) for 24 h after lysis. The five T4SS transcription units of E. ruminantium were significantly up-regulated under iron depletion (25 to 190 fold increase as shown for virB3, virB4a, virB4b, virB8a, virB8b, and virB9b) (Figure 4A). virB2a was the only gene that did not present a significant change in expression (Figure 4A). Moreover, the erxR gene was up-regulated 25 fold during iron starvation (Figure 4B). We tested the expression of the map1 gene cluster under iron starvation in the two genes at the border of the cluster (tr1 and map1+1) as well as in the archetypal map1 and map1-6, which are located in a central position in the cluster. Culture under iron limitation strongly increased the expression of all genes (Figure 4C), however map1+1 showed the highest change in expression with 400-fold up-regulation (Figure 4C), suggesting that this gene may play a key role during iron starvation.

FIGURE 4

Figure 4. virBD, tr1 and map genes as well as erxR are up-regulated under iron-depletion. The expression of the virB (A), erxR (B), tr1 and map1 (C) genes were measured during the lysis phase of infection under iron repletion or iron depletion using quantitative real-time PCR. The data represent the mean + SD of 2 or 3 biological replicates, each of which comprised 3 technical replicates. Ratios were calculated from the transcript numbers and normalized to recA. The resol gene was used as negative control.

Discussion

The preferential expression of virulence factors with diverse functions in response to host and environmental cues has already been characterized in several bacteria (Oogai et al., 2011; Weber et al., 2014). Thus, several Ehrlichia and other members of the Anaplasmataceae family have been shown to differentially express certain genes in response to the host cells (Singu et al., 2006; Nelson et al., 2008).

These genes include genes encoding components of the T4SS, which have been shown to play an essential role in pathogenicity (Rikihisa, 2010; Moumène and Meyer, 2015). Many intracellular bacterial pathogens use T4SS to deliver effector molecules, which subvert the eukaryotic host cell defenses and other cellular processes to their own advantage (Trokter et al., 2014). The genetic arrangement of the 18 genes encoding the T4SS apparatus in E. ruminantium T4SS genes resembles that of E. chaffeensis (Figure 1).

Likewise, OMPs are often regulated by environmental signals and play an important role in bacterial pathogenesis by enhancing the ability of the bacteria to adapt to different environments (Blanvillain et al., 2007). Interestingly, T4SS affects outer membrane properties, which might be important for the adaptation of Brucella to growth both in vitro and in vivo (Wang et al., 2010). Herein, we have described the effects of one environmental cue, iron concentration, on the expression of the transcriptional regulator erxR, the components of the T4SS apparatus and members of the Map1 family during E. ruminantium infection cycle of BAE cells.

Regulatory proteins are known to play an important role in the survival and persistence of intracellular pathogens in their host (Cheng et al., 2008). EcxR (ECH_0795) is the only transcriptional regulator to be associated with the expression of T4SS components in E. chaffeensis (Cheng et al., 2008). In E. ruminantium, we identified erxR, an orthologous gene of ecxR (Figure 1A), by sequence homology. ErxR (Ehrlichia ruminantium expression regulator) binds and regulates its own promoter and the promoters of some virBD genes. These results show that ErxR is the regulatory protein of the T4SS of E. ruminantium. Using qRT-PCR at different time points, we also show that erxR is strongly expressed at an elementary body stage (120 hpi), like that demonstrated in E. chaffeensis (Cheng et al., 2008).

ErxR orthologs in other Anaplasmataceae also appear to be associated with the expression of important antigenic OMPs. Thus, ApxR regulates the transcription of p44 transcription by binding to the tr1 promoter during A. phagocytophilum infection of mammalian host cells (Wang et al., 2007a,b). p44E encodes the immunodominant pleomorphic 44-kDa major surface protein, which shows homology with the Map1 family in E. ruminantium as well as the P30 family in E. canis and the P28 family in E. chaffeensis (Dunning Hotopp et al., 2006). The major antigenic protein Map1 is part of a multigene family containing 16 paralogs tandemly organized in a head to tail arrangement that are located downstream of a hypothetical transcriptional regulator gene (tr1) (Postigo et al., 2007; Figure 1C), a similar arrangement to that reported for p44, p30, and p28 in the other Anaplasmataceae (Dunning Hotopp et al., 2006). tr1 is one of the three promoters identified in the p44 expression locus and shown to be the strongest promoter driving the expression of a polycistronic mRNA containing OMP1, p44ESup, and p44 (Barbet et al., 2005). tr1 harbors a winged helix-turn-helix and a DNA binding motif and its part of the xenobiotic response element family of transcriptional regulators. However, the function of tr1 remains unclear (Nelson et al., 2008), and whether or not tr1 drives the expression of a polycistronic tandem mRNA containing several map homologs is still not known. According to our results, ErxR binds to tr1 and to map1 promoters (Figure 3). It is thus possible that it regulates the expression of the map1 members in a similar way to tr1 in A. phagocytophilum for p44 expression.

Taken together, our results show that ErxR binds to tr1, map1, and virBD promoters, suggesting for the first time coordinated regulation of T4SS and OMP in Anaplasmataceae. Interestingly, in the endosymbiotic bacterium Wolbachia (wBm) wBmxR1 and wBmxR2, which are orthologs of ErxR, have been shown to co-regulate genes of the T4SS and riboflavin biosynthesis pathway (Li and Carlow, 2012). Riboflavin is an important co-factor for the survival of the endosymbiont's host, the filarial parasite Brugia malayi (Li and Carlow, 2012).

Next, we investigated whether the regulation of genes encoding the T4SS and members of the Map1 family by ErxR could be triggered by environmental and nutritional cues in the host cell. Our transcriptional analysis showed that erxR, virBD, tr1, and map1 genes were upregulated in response to iron starvation (Figure 4). Interestingly, the fact that virB2a is not expressed under iron depletion could be due to a functional redundancy depending on the environment (e.g. mammalian host, vector cell, etc.) as previously shown for PopF1 and PopF2 proteins of the T3SS of Ralstonia solanacearum (Meyer et al., 2006). Similarly, A. phagocytophilum appears to express specific virB2 paralogs in a host cell dependent manner as well as differential expression of virB2 gene in tick cells or in human cells (Nelson et al., 2008). Many bacterial pathogens sense iron depletion as a signal indicating that they are within a vertebrate host (Skaar, 2010).

Although the sensing mechanisms of changes in iron concentration are not known in E. ruminantium, three pathways are possible. One possibility is that the regulator ErxR is activated by an unidentified sensor kinase which responds to this environmental signal and up-regulates expression of the T4SS genes and certain map1 genes in response to iron starvation (Figure 5). Sensor kinases are part of two component systems (TCS), which regulate the differential expression of genes in bacteria in response to environmental cues. Three TCS composed of three response regulators and three sensor kinases have been identified in A. phagocytophilum and E. chaffeensis (Cheng et al., 2006). All three sets of TCS have orthologs in E. ruminantium and one of these TCS could thus be involved in the regulation of pathogenic factors in response to changes in iron abundance.

FIGURE 5

Figure 5. Putative regulation model for ErxR, virBD, tr1 and map genes in E. ruminantium under iron starvation conditions.

Another possibility is the presence of an unknown Fur repressor in the Anaplasmataceae family. This sensing involves transcriptional control mediated by the transcriptional repressor, Fur (Escolar et al., 1998). But many organisms, including E. coli, Campylobacter jejuni, and Vibrio cholera, have been shown to use Fur to negatively regulate gene expression with increasing concentrations of iron (O'Sullivan et al., 1994; Palyada et al., 2004; Mey et al., 2005). For example, Fur activates sodB, an iron superoxide dismutase, under iron repletion. SODs are metalloproteins, which play an important role in protection against oxidative stress by catalyzing dismutation of the superoxide radical (${\text{O}}_2^{-}$). Interestingly, in E. ruminantium, sodB is located upstream of operon 2 of the T4SS and is co-transcribed along these genes in Ehrlichia and Anaplasma, suggesting an effect of iron in the expression of virBD genes. Moreover, it has been demonstrated that the regulator ApxR binds to the promoter regions upstream of sodB (Wang et al., 2007b). We searched for putative Fur boxes in the genome of E. ruminantium and found one located upstream of erxR (Supplementary Figure 1). The 19 bp sequence consisted in two repeated hexamers (nATWAT) flanking a 7 nt sequence, which is commonly found upstream of iron regulated enzymes such as succinate dehydrogenase iron-sulfur subunits, major ferric iron binding protein precursors, and adenosine tRNA methylthiotransferase (Escolar et al., 1998; Grifantini et al., 2003). Thus, it is possible that E. ruminantium is capable of sensing low iron concentrations in the environment and of regulating the expression of ErxR through this putative Fur box.

Finally, it is possible that some Map1 proteins play a role like that of TBDRs in the perception of environmental cues and in iron uptake (Blanvillain et al., 2007). As mentioned above, Map1 proteins are orthologs of members of the p28 family in E. chaffeensis, which have been shown to function as porins and possibly act in nutrient uptake during intracellular infection (Kumagai et al., 2008). It has been shown that porins, such as OmpA and OmpC, bind to transferrin and act in iron uptake in the enteropathogenic strains of E. coli, Salmonella typhimurium, and several Shigella species (Sandrini et al., 2013). Similarly, Mycobacterium smegmatis is able to acquire ferric ions through members of the Msp family of porins (Jones and Niederweis, 2010). The up-regulation of map1, map1+1, and map1-6 showed under iron starvation suggests that these three porins play a role in iron acquisition and raises the possibility that Map1 proteins may also act as sensors, although more evidence is needed before concluding. These results suggest that the Map1 proteins may fulfil several functions during E. ruminantium infection. Characterizing these functions could advance our understanding of the adaptation of E. ruminantium to its host, as done by Blanvillain et al. (2007).

In conclusion, we have demonstrated that exposure of E. ruminantium to iron limitation induces ErxR-dependent expression of the T4SS apparatus and map1 genes. These findings reveal an important degree of coordination between T4SS and map1 genes at the transcriptional level and raise the possibility of the involvement of Map proteins in environmental sensing and in the infection process. The data presented herein enables us to propose a model for the regulation of E. ruminantium T4SS and map1 genes in which ErxR acts as a global regulator integrating iron as an important triggering environmental signal (Figure 5). Understanding Ehrlichia gene regulation in response to environmental signals provides valuable cues for the development of alternative treatments.

Author Contributions

DM, Conceived and designed the experiments. AM, SG-R, and DM, Performed the experiments. AM, SG-R, TL, NV, and DM, Analyzed the data. AM and DM, 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.

Acknowledgments

We are grateful to. S. Castaño-Cerezo at the Universidad de Murcia for providing pUA66 plasmid. We thank E. Albina for critical reading of the manuscript and A.O. Chavez for technical assistance in the regulation experiments. AM acknowledges financial support from FEDER grant FED 1/1.4-30305, 2007–2013, “Risque en santé animale et végétale” (PhD grant to AM), and from EPIGENESIS RegPot European project (n°315988).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2017.00535/full#supplementary-material

Supplementary Figure 1 (A). Nucleotide sequence of the predicted erxR promoter. The initiation codon of erxR is shown in capital letters at the end of the sequence. The sequence in red represents the nATWAT motif. The−35 and−10 boxes of the erxR promoter are capitalized. (B) The potential erxR Fur box is compared with the consensus E. coli Fur box. Conserved residues are shown in red.

Supplementary Table 1. Oligonucleotide primers used for qPCR and qRT-PCR.

Supplementary Table 2. Oligonucleotide primers used to construct plasmids.

Supplementary Table 3. Strains and plasmids used in this work.

References

Abromaitis, S., Nelson, C. S., Previte, D., Yoon, K. S., Clark, J. M., DeRisi, J. L., et al. (2013). Bartonella quintana deploys host and vector temperature-specific transcriptomes. PLoS ONE 8:e58773. doi: 10.1371/journal.pone.0058773

PubMed Abstract | CrossRef Full Text | Google Scholar

Allsopp, B. A. (2010). Natural history of Ehrlichia ruminantium. Vet. Parasitol. 167, 123–135. doi: 10.1016/j.vetpar.2009.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrews, S. C., Robinson, A. K., and Rodríguez-Qui-ones, F. (2003). Bacterial iron homeostasis. FEMS Microbiol. Rev. 27, 215–237. doi: 10.1016/S0168-6445(03)00055-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Aravind, L., Anantharaman, V., Balaji, S., Babu, M. M., and Iyer, L. M. (2005). The many faces of the helix-turn-helix domain: transcription regulation and beyond. FEMS Microbiol. Rev. 29, 231–262. doi: 10.1016/j.fmrre.2004.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Barbet, A. F., Agnes, J. T., Moreland, A. L., Lundgren, A. M., Alleman, A. R., Noh, S. M., et al. (2005). Identification of functional promoters in the msp2 expression loci of Anaplasma marginale and Anaplasma phagocytophilum. Gene 353, 89–97. doi: 10.1016/j.gene.2005.03.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanvillain, S., Meyer, D., Boulanger, A., Lautier, M., Guynet, C., Denancé, N., et al. (2007). Plant carbohydrate scavenging through tonB-dependent receptors: a feature shared by phytopathogenic and aquatic bacteria. PLoS ONE 2:e224. doi: 10.1371/journal.pone.0000224

PubMed Abstract | CrossRef Full Text | Google Scholar

Breuer, W., Epsztejn, S., and Cabantchik, Z. I. (1995). Iron acquired from transferrin by K562 cells is delivered into a cytoplasmic pool of chelatable iron(II). J. Biol. Chem. 270, 24209–24215. doi: 10.1074/jbc.270.41.24209

PubMed Abstract | CrossRef Full Text | Google Scholar

Brickman, T. J., Cummings, C. A., Liew, S.-Y., Relman, D. A., and Armstrong, S. K. (2011). Transcriptional profiling of the iron starvation response in Bordetella pertussis provides new insights into siderophore utilization and virulence gene expression. J. Bacteriol. 193, 4798–4812. doi: 10.1128/JB.05136-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Castaño-Cerezo, S., Bernal, V., Blanco-Catalá, J., Iborra, J. L., and Cánovas, M. (2011). cAMP-CRP co-ordinates the expression of the protein acetylation pathway with central metabolism in Escherichia coli. Mol. Microbiol. 82, 1110–1128. doi: 10.1111/j.1365-2958.2011.07873.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, Z., Kumagai, Y., Lin, M., Zhang, C., and Rikihisa, Y. (2006). Intra-leukocyte expression of two-component systems in Ehrlichia chaffeensis and Anaplasma phagocytophilum and effects of the histidine kinase inhibitor closantel. Cell. Microbiol. 8, 1241–1252. doi: 10.1111/j.1462-5822.2006.00704.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, Z., Wang, X., and Rikihisa, Y. (2008). Regulation of type IV secretion apparatus genes during Ehrlichia chaffeensis intracellular development by a previously unidentified protein. J. Bacteriol. 190, 2096–2105. doi: 10.1128/JB.01813-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Collins, N. E., Liebenberg, J., de Villiers, E. P., Brayton, K. A., Louw, E., Pretorius, A., et al. (2005). The genome of the heartwater agent Ehrlichia ruminantium contains multiple tandem repeats of actively variable copy number. Proc. Natl. Acad. Sci. U.S.A. 102, 838–843. doi: 10.1073/pnas.0406633102

PubMed Abstract | CrossRef Full Text | Google Scholar

Dumler, J. S., Barbet, A. F., Bekker, C. P., Dasch, G. A., Palmer, G. H., Ray, S. C., et al. (2001). Reorganization of genera in the families Rickettsiaceae and Anaplasmataceae in the order Rickettsiales: unification of some species of Ehrlichia with Anaplasma, Cowdria with Ehrlichia and Ehrlichia with Neorickettsia, descriptions of six new species combinations and designation of Ehrlichia equi and “HGE agent” as subjective synonyms of Ehrlichia phagocytophila. Int. J. Syst. Evol. Microbiol. 51, 2145–2165. doi: 10.1099/00207713-51-6-2145

CrossRef Full Text | Google Scholar

Dunning Hotopp, J. C., Lin, M., Madupu, R., Crabtree, J., Angiuoli, S. V., Eisen, J. A., et al. (2006). Comparative genomics of emerging human ehrlichiosis agents. PLoS Genet. 2:e21. doi: 10.1371/journal.pgen.0020021

PubMed Abstract | CrossRef Full Text | Google Scholar

Escolar, L., Pérez-Martín, J., and de Lorenzo, V. (1998). Binding of the fur (ferric uptake regulator) repressor of Escherichia coli to arrays of the GATAAT sequence. J. Mol. Biol. 283, 537–547.

PubMed Abstract | Google Scholar

Frutos, R., Viari, A., Vachiery, N., Boyer, F., and Martinez, D. (2007). Ehrlichia ruminantium: genomic and evolutionary features. Trends Parasitol. 23, 414–419. doi: 10.1016/j.pt.2007.07.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Rizzo, S., Crespi, M., and Frugier, F. (2006). The Medicago truncatula CRE1 cytokinin receptor regulates lateral root development and early symbiotic interaction with Sinorhizobium meliloti. Plant Cell 18, 2680–2693. doi: 10.1105/tpc.106.043778

PubMed Abstract | CrossRef Full Text | Google Scholar

Grifantini, R., Sebastian, S., Frigimelica, E., Draghi, M., Bartolini, E., Muzzi, A., et al. (2003). Identification of iron-activated and -repressed Fur-dependent genes by transcriptome analysis of Neisseria meningitidis group B. Proc. Natl. Acad. Sci. U.S.A. 100, 9542–9547. doi: 10.1073/pnas.1033001100

CrossRef Full Text | Google Scholar

Hyytiäinen, H., Sjöblom, S., Palomäki, T., Tuikkala, A., and Tapio Palva, E. (2003). The PmrA-PmrB two-component system responding to acidic pH and iron controls virulence in the plant pathogen Erwinia carotovora ssp. carotovora. Mol. Microbiol. 50, 795–807. doi: 10.1046/j.1365-2958.2003.03729.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, C. M., and Niederweis, M. (2010). Role of porins in iron uptake by Mycobacterium smegmatis. J. Bacteriol. 192, 6411–6417. doi: 10.1128/JB.00986-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumagai, Y., Huang, H., and Rikihisa, Y. (2008). Expression and porin activity of P28 and OMP-1F during intracellular Ehrlichia chaffeensis development. J. Bacteriol. 190, 3597–3605. doi: 10.1128/JB.02017-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Larkin, M. A., Blackshields, G., Brown, N. P., Chenna, R., McGettigan, P. A., McWilliam, H., et al. (2007). Clustal, W., and Clustal X version 2.0. Bioinformatics 23, 2947–2948. doi: 10.1093/bioinformatics/btm404

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., and Carlow, C. K. S. (2012). Characterization of transcription factors that regulate the type IV secretion system and riboflavin biosynthesis in Wolbachia of Brugia malayi. PLoS ONE 7:e51597. doi: 10.1371/journal.pone.0051597

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcelino, I., Veríssimo, C., Sousa, M. F. Q., Carrondo, M. J. T., and Alves, P. M. (2005). Characterization of Ehrlichia ruminantium replication and release kinetics in endothelial cell cultures. Vet. Microbiol. 110, 87–96. doi: 10.1016/j.vetmic.2005.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Martín-Martín, A. I., Sancho, P., de Miguel, M. J., Fernández-Lago, L., and Vizcaíno, N. (2012). Quorum-sensing and BvrR/BvrS regulation, the type IV secretion system, cyclic glucans, and BacA in the virulence of Brucella ovis: similarities to and differences from smooth brucellae. Infect. Immun. 80, 1783–1793. doi: 10.1128/IAI.06257-11

PubMed Abstract | CrossRef Full Text | Google Scholar

Mey, A. R., Wyckoff, E. E., Kanukurthy, V., Fisher, C. R., and Payne, S. M. (2005). Iron and fur regulation in Vibrio cholerae and the role of fur in virulence. Infect. Immun. 73, 8167–8178. doi: 10.1128/IAI.73.12.8167-8178.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Meyer, D. F., Cunnac, S., Guéneron, M., Declercq, C., Van Gijsegem, F., Lauber, E., et al. (2006). PopF1 and PopF2, two proteins secreted by the type III protein secretion system of Ralstonia solanacearum, are translocators belonging to the HrpF/NopX family. J. Bacteriol. 88, 4903–4917 doi: 10.1128/JB.00180-06

CrossRef Full Text | Google Scholar

Moumène, A., and Meyer, D. F. (2015). Ehrlichia's molecular tricks to manipulate their host cells. Microbes Infect. 18, 172–179. doi: 10.1016/j.micinf.2015.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, C. M., Herron, M. J., Felsheim, R. F., Schloeder, B. R., Grindle, S. M., Chavez, A. O., et al. (2008). Whole genome transcription profiling of Anaplasma phagocytophilum in human and tick host cells by tiling array analysis. BMC Genomics 9:364. doi: 10.1186/1471-2164-9-364

PubMed Abstract | CrossRef Full Text | Google Scholar

Niu, H., Rikihisa, Y., Yamaguchi, M., and Ohashi, N. (2006). Differential expression of VirB9 and VirB6 during the life cycle of Anaplasma phagocytophilum in human leucocytes is associated with differential binding and avoidance of lysosome pathway. Cell. Microbiol. 8, 523–534. doi: 10.1111/j.1462-5822.2005.00643.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Oogai, Y., Matsuo, M., Hashimoto, M., Kato, F., Sugai, M., and Komatsuzawa, H. (2011). Expression of virulence factors by Staphylococcus aureus grown in serum. Appl. Environ. Microbiol. 77, 8097–8105. doi: 10.1128/AEM.05316-11

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Sullivan, D. J., Dowling, D. N., deLorenzo, V., and O'Gara, F. (1994). Escherichia coli ferric uptake regulator (Fur) can mediate regulation of a pseudomonad iron-regulated promoter. FEMS Microbiol. Lett. 117, 327–332. doi: 10.1111/j.1574-6968.1994.tb06787.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Palyada, K., Threadgill, D., and Stintzi, A. (2004). Iron acquisition and regulation in Campylobacter jejuni. J. Bacteriol. 186, 4714–4729. doi: 10.1128/JB.186.14.4714-4729.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Portier, E., Zheng, H., Sahr, T., Burnside, D. M., Mallama, C., Buchrieser, C., et al. (2014). IroT/mavN, a new iron-regulated gene involved in Legionella pneumophila virulence against amoebae and macrophages. Environ. Microbiol. 17, 1338–13350. doi: 10.1111/1462-2920.12604

PubMed Abstract | CrossRef Full Text | Google Scholar

Postigo, M., Taoufik, A., Bell-Sakyi, L., de Vries, E., Morrison, W. I., and Jongejan, F. (2007). Differential transcription of the major antigenic protein 1 multigene family of Ehrlichia ruminantium in Amblyomma variegatum ticks. Vet. Microbiol. 122, 298–305. doi: 10.1016/j.vetmic.2007.01.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruneau, L., Emboulé, L., Gely, P., Marcelino, I., Mari, B., Pinarello, V., et al. (2012). Global gene expression profiling of Ehrlichia ruminantium at different stages of development. FEMS Immunol. Med. Microbiol. 64, 66–73. doi: 10.1111/j.1574-695X.2011.00901.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ratledge, C., and Dover, L. G. (2000). Iron metabolism in pathogenic bacteria. Annu. Rev. Microbiol. 54, 881–941. doi: 10.1146/annurev.micro.54.1.881

PubMed Abstract | CrossRef Full Text | Google Scholar

Rikihisa, Y. (2010). Anaplasma phagocytophilum and Ehrlichia chaffeensis: subversive manipulators of host cells. Nat. Rev. Microbiol. 8, 328–339. doi: 10.1038/nrmicro2318

PubMed Abstract | CrossRef Full Text | Google Scholar

Rikihisa, Y. (2011). Mechanisms of obligatory intracellular infection with Anaplasma phagocytophilum. Clin. Microbiol. Rev. 24, 469–489. doi: 10.1128/CMR.00064-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Romeo, A. M., Christen, L., Niles, E. G., and Kosman, D. J. (2001). Intracellular chelation of iron by bipyridyl inhibits DNA virus replication: ribonucleotide reductase maturation as a probe of intracellular iron pools. J. Biol. Chem. 276, 24301–24308. doi: 10.1074/jbc.M010806200

PubMed Abstract | CrossRef Full Text | Google Scholar

Sandrini, S., Masania, R., Zia, F., Haigh, R., and Freestone, P. (2013). Role of porin proteins in acquisition of transferrin iron by enteropathogens. Microbiology 159, 2639–2650. doi: 10.1099/mic.0.071928-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Shultis, D. D., Purdy, M. D., Banchs, C. N., and Wiener, M. C. (2006). Outer membrane active transport: structure of the BtuB:TonB complex. Science 312, 1396–1399. doi: 10.1126/science.1127694

PubMed Abstract | CrossRef Full Text | Google Scholar

Singu, V., Peddireddi, L., Sirigireddy, K. R., Cheng, C., Munderloh, U., and Ganta, R. R. (2006). Unique macrophage and tick cell-specific protein expression from the p28/p30-outer membrane protein multigene locus in Ehrlichia chaffeensis and Ehrlichia canis. Cell. Microbiol. 8, 1475–1487. doi: 10.1111/j.1462-5822.2006.00727.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Skaar, E. P. (2010). The battle for iron between bacterial pathogens and their vertebrate hosts. PLoS Pathog. 6:e1000949. doi: 10.1371/journal.ppat.1000949

PubMed Abstract | CrossRef Full Text | Google Scholar

Trokter, M., Felisberto-Rodrigues, C., Christie, P. J., and Waksman, G. (2014). Recent advances in the structural and molecular biology of type IV secretion systems. Curr. Opin. Struct. Biol. 27, 16–23. doi: 10.1016/j.sbi.2014.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

van Heerden, H., Collins, N. E., Brayton, K. A., Rademeyer, C., and Allsopp, B. A. (2004). Characterization of a major outer membrane protein multigene family in Ehrlichia ruminantium. Gene 330, 159–168. doi: 10.1016/j.gene.2004.01.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Cheng, Z., Zhang, C., Kikuchi, T., and Rikihisa, Y. (2007a). Anaplasma phagocytophilum p44 mRNA expression is differentially regulated in mammalian and tick host cells: involvement of the DNA binding protein ApxR. J. Bacteriol. 189, 8651–8659. doi: 10.1128/JB.00881-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Kikuchi, T., and Rikihisa, Y. (2007b). Proteomic identification of a novel Anaplasma phagocytophilum DNA binding protein that regulates a putative transcription factor. J. Bacteriol. 189, 4880–4886. doi: 10.1128/JB.00318-07

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Chen, Z., Qiao, F., Zhong, Z., Xu, J., Wang, Z., et al. (2010). The type IV secretion system affects the expression of Omp25/Omp31 and the outer membrane properties of Brucella melitensis. FEMS Microbiol. Lett. 303, 92–100. doi: 10.1111/j.1574-6968.2009.01866.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Weber, G. G., Kortmann, J., Narberhaus, F., and Klose, K. E. (2014). RNA thermometer controls temperature-dependent virulence factor expression in Vibrio cholerae. Proc. Natl. Acad. Sci. U.S.A. 111, 14241–14246. doi: 10.1073/pnas.1411570111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zweygarth, E., and Josemans, A. I. (2001). Continuous in vitro propagation of Cowdria ruminantium (Welgevonden stock) in a canine macrophage-monocyte cell line. Onderstepoort J. Vet. Res. 68, 155–157. Available online at: http://hdl.handle.net/2263/18399

PubMed Abstract | Google Scholar