Comparative proteomic profiling of Ehrlichia ruminantium pathogenic strain and its high-passaged attenuated Strain reveals virulence and attenuation-associated proteins

Isabel Marcelino, Miguel Ventosa, Elisabete Pires, Markus Müller, Frédérique Lisacek, Thierry Lefrancois, Nathalie Vachery, Ana Varela Coelho

<figure><img src="image_1.png" /></figure>

► To cite this version:

Isabel Marcelino, Miguel Ventosa, Elisabete Pires, Markus Müller, Frédérique Lisacek, et al.. Comparative proteomic profiling of Ehrlichia ruminantium pathogenic strain and its high-passaged attenuated Strain reveals virulence and attenuation-associated proteins. PLoS ONE, 2015, 10 (12), 23 p. 10.1371/journal.pone.0145328. [hal-01595370](https://hal.science/hal-01595370)

HAL Id: hal-01595370

https://hal.science/hal-01595370

Submitted on 26 Sep 2017

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Comparative Proteomic Profiling of *Ehrlichia ruminantium* Pathogenic Strain and Its High-Passaged Attenuated Strain Reveals Virulence and Attenuation-Associated Proteins

Isabel Marcelino¹,²,³,⁴*, Miguel Ventosa¹,², Elisabete Pires², Markus Müller⁵,
Frédérique Lisacek⁵, Thierry Lefrançois⁶, Nathalie Vachieri³,⁴, Ana Varela Coelho²

¹ Instituto de Biologia Experimental e Tecnológica (iBET), Oeiras, Portugal, ² Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB-UNL), Oeiras, Portugal, ³ CIRAD, UMR CMAEE, F-97170, Petit-Bourg, Guadeloupe, France, ⁴ INRA, UMR1309 CMAEE, F-34398, Montpellier, France, ⁵ Proteome Informatics, Swiss Institute of Bioinformatics (SIB), Geneva, Switzerland, ⁶ CIRAD, UMR CMAEE, F-34398, Montpellier, France

* Isabel.Marcelino@cirad.fr

<figure><img src="image_3.png" /></figure>

CrossMark
click for updates

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Marcelino I, Ventosa M, Pires E, Müller M, Lisacek F, Lefrançois T, et al. (2015) Comparative Proteomic Profiling of *Ehrlichia ruminantium* Pathogenic Strain and Its High-Passaged Attenuated Strain Reveals Virulence and Attenuation-Associated Proteins. PLoS ONE 10(12): e0145328. doi:10.1371/journal.pone.0145328

**Editor:** Xue-jie Yu, University of Texas Medical Branch, UNITED STATES

**Received:** August 27, 2015

**Accepted:** December 1, 2015

**Published:** December 21, 2015

**Copyright:** © 2015 Marcelino et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the Fundação para a Ciência e Tecnologia (FCT) through the PTDC/CVT/114118/2009 research project, grant # Pest-OE/EQB/LA0004/2011 and the National Research Program for "Rede Nacional de Espectrometria de Massa - RNEM" (REDE/1504/RNEM/2005). The authors MV and IM also acknowledge funding from the EU projects COST action FA-1002 (for networking opportunities) and

## Abstract

The obligate intracellular bacterium *Ehrlichia ruminantium* (ER) causes heartwater, a fatal tick-borne disease in livestock. In the field, ER strains present different levels of virulence, limiting vaccine efficacy, for which the molecular basis remains unknown. Moreover, there are no genetic tools currently available for ER manipulation, thus limiting the knowledge of the genes/proteins that are essential for ER pathogenesis and biology. As such, to identify proteins and/or mechanisms involved in ER virulence, we performed the first exhaustive comparative proteomic analysis between a virulent strain (ERGvir) and its high-passaged attenuated strain (ERGatt). Despite their different behaviors *in vivo* and *in vitro*, our results from 1DE-nanoLC-MS/MS showed that ERGvir and ERGatt share 80% of their proteins; this core proteome includes chaperones, proteins involved in metabolism, protein-DNA-RNA biosynthesis and processing, and bacterial effectors. Conventional 2DE revealed that 85% of the identified proteins are proteoforms, suggesting that post-translational modifications (namely glycosylation) are important in ER biology. Strain-specific proteins were also identified: while ERGatt has an increased number and overexpression of proteins involved in cell division, metabolism, transport and protein processing, ERGvir shows an overexpression of proteins and proteoforms (DIGE experiments) involved in pathogenesis such as Lpd, AnkA, VirB9 and B10, providing molecular evidence for its increased virulence *in vivo* and *in vitro*. Overall, our work reveals that ERGvir and ERGatt proteomes are streamlined to fulfill their biological function (maximum virulence for ERGvir and replicative capacity for ERGatt), and we provide both pioneering data and novel insights into the pathogenesis of this obligate intracellular bacterium.
FEDER 2007-2013 project (FED 1/1.4-30305). IM acknowledges financial support from Post-doc grant SFRH/BPD/45978/2008 from FCT (Portugal) and from EPIGENESIS RegPot European project (n° 315988). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

**Abbreviations:** ER, *Ehrlichia ruminantium*; ERGvir, *Ehrlichia ruminantium* Gardel virulent strain; ERGatt, *Ehrlichia ruminantium* Gardel attenuated strain; 1DE, standard SDS-PAGE; nanoLC, nano-scale liquid chromatography; 2DE, Bidimensional electrophoresis; PTM, post-translational modification; SPG, sucrose-phosphate-glutamate buffer; DIGE, Differential In-Gel Electrophoresis.

## Introduction

*Ehrlichia ruminantium* (ER), an obligate intracellular bacterium of the order *Rickettsiales*, causes heartwater, a fatal and economically important disease of domestic and wild ruminants. This tick-borne disease occurs throughout sub-Saharan Africa and some islands in the Indian Ocean and the Caribbean, from where it threatens to invade the Americas, posing a serious threat to the livestock industry [1].

Several candidate vaccines such as inactivated, attenuated, DNA and subunit vaccines are available [2–6], but the development of a fully effective vaccine has been hindered by the difficulty of identifying protective antigens. This is due to the limited knowledge of the protective immune responses against heartwater and the lack of accurate *in vitro* antigen screening tests. Moreover, ER has a high antigenic diversity inducing low cross-protection between strains [1,7–12]. Several virulent ER strains have been reported [7–12] and at the moment, only two strains have been recognized as being attenuated by serial passage in bovine endothelial cells or canine macrophages: Senegal (from West Africa [13]) and Welgevonden (from South Africa [4]), respectively. A third one, Gardel strain (originated from Guadeloupe, French West Indies) has been suggested to be attenuated after roughly 200 passages in host endothelial cells [14] but no *in vivo* nor *in vitro* results have been so far published. Despite the recent advances in the study of ER biology [15–19], the lack of genetic tools for ER manipulation hampers knowledge of the genes that are actually expressed in live bacteria, namely those associated with ER virulence and attenuation.

Disparity in the pathogenicity or other biological aspects observed between strains can be related to important differences at the genomic level [20–22]. In ER, only 64% of the genome is predicted to be coding sequence (CDS), encoding 950 proteins and 39 stable RNA species [23]; 30% of these CDSs have unknown function. Genome analysis also indicates the presence of several classes of predictive virulence factors, such as genes/proteins involved in secretion and/or the trafficking of molecules between the pathogen and host cells, or evasion and/or modulation of the host immune system [16,17,24,25]. Still, comparative genomic analyses between the genomes of Gardel and Welgevonden strains indicate few genomic differences [25]. The gain of pathogenicity in bacteria can result from horizontal gene transfer, either directly or through mobile genetic elements [26]. In the case of an intracellular parasite such as ER, only one species of intracellular parasite inhabits a host cell, which restricts the parasite’s access to new genes, so few genes acquired by horizontal gene transfer are to be expected [24].

In this context, it can be assumed that differences between ER virulent and attenuated strains might not be predominantly at the genomic level but may be associated with variations in expression levels of genes/proteins encoding virulence factors, as previously observed for *Mycoplasma hyopneumoniae* [20].

Recent transcriptomic studies of ERGvir during its development cycle revealed that the extracellular and infectious form of the bacterium (elementary bodies, EBs) presented an over-expression of DksA [17], a transcription factor known to be related to virulence in other bacteria [27,28]. Preliminary proteomic studies on EBs also highlighted the expression of proteins related to virulence namely from the Type IV secretion system (T4SS, VirB9 and VirB11), and T4SS effectors AnkA, DksA and ElbB [29,30]. Still, no transcriptomic or proteomic studies comparing virulent and attenuated ER strains have been published. In fact, global “omics” studies on *Rickettsiales* species comparing virulent and avirulent strains are so far only available for *Rickettsia prowazekii* [22,31,32] and for *Rickettsia peacockii* [33].

In this study, we first confirmed the attenuation of Gardel after 200 passages *in vivo* and characterized ERGatt growth kinetics *in vitro*. Afterwards, we compared the proteomes of the ERGvir and ERGatt strains using gel-based separation approaches associated with
MALDI-TOF/TOF mass spectrometry. Proteome analysis using 1DE-nLC-MS/MS revealed that both strains share 80% of their proteins (292 non-redundant protein), which represent approximately 31% (292/950) of the genome coding capacity of this bacterium. Surprisingly, 85% of the proteins identified by conventional 2DE electrophoresis are proteoforms. We also performed a quantitative proteomic analysis using DIGE to assess differentially expressed proteins between the two strains. This study is the first comprehensive and comparative proteomic analysis of virulent and attenuated ER strains and it provides an important proteomic basis for ER pathogenesis and suggests an important role of proteoforms in ER biology.

## Materials and Methods

### Experimental Design and Statistical Analyses

For the 1DE-nanoLC experiments and conventional 2DE gel analyses, three independent biological replicates per strain were used. Results from 1-DE-nanoLC-MS/MS experiments were analysed using PEAKS search engine tool (PEAKS Studio 5.3; Bioinformatics Solutions Inc., Waterloo, ON, Canada) [34], combining three search engines MASCOT, X!Tandem and Peaks DB. Results from conventional 2DE gels were analyzed and compiled using MDM software [35]. For DIGE experiments, four independent biological replicates per strain were used; image analyses were performed using Progenesis SameSpots v3.0 software and spot-normalized volume was used to select statistically significant (fold-change, ANOVA, false discovery rate and power value) differentiated spots between ERG strains analyzed in the experiment. All these parameters are described below in more detail.

### *E. ruminantium* cultivation and growth kinetics

Infected blood from ERG-infected goats was loaded on finite cultures of bovine aortic endothelial cells (BAE) to isolate the ERGardel virulent strain (ERGvir). ERGvir was then routinely propagated in BAE as described elsewhere [36,37]; ERGvir samples from passages up to passage 44 were used throughout this study, as they have been proven to be highly infectious in previous studies [38]. In order to obtain an attenuated ERG strain, ERGvir was cultivated over 230 passages in BAE cells, as suggested by Martinez (1987) [14]. When 80% cell lysis was observed (at 120hpi for ERGvir and 96hpi for ERGatt, S1 Fig), supernatant and cellular debris containing infectious ER elementary bodies were harvested and then used to (i) infect a freshly confluent monolayer or (ii) to be purified using a multistep centrifugation methodology [39]. Purified ERs were stored in SPG [39] at -80°C with a "Complete EDTA-free" anti-protease cocktail (Roche, Germany) prior to proteomic analysis. Prior to *in vivo* assays, purified ER were stored in SPG in liquid nitrogen [38].

Growth of ERGvir and ERGatt was monitored by phase-contrast light microscopy and quantified as previously described [17].

### Animal infection studies

To evaluate the virulent and attenuated status of ERGvir and ERGatt respectively, *in vivo* experiments were performed using naive goats from Les Saintes Island (Guadeloupe, FWI), a heartwater-free region. The ER inoculum was thawed and diluted in cold culture medium for immediate use on animals. Two goats (0541 and 0614) were infected intravenously with 72 x 10⁴ and 72 x 10⁵ viable elementary bodies from ERG passage 230 (ERGatt); the doses were calibrated as previously described [38] and correspond, respectively, to 2.4 and 24 times the lethal dose for the ERGvir strain [38–40]. As a positive control, one goat (0212) was infected with 72 x 10⁴ viable ERG passage 27. This part of the study was conducted in 2011–2012 according to
internationally approved World Organization for Animal Health (OIE) standards, and CIRAD institution was authorized by the director of the veterinary services of Guadeloupe on behalf of the Prefect of Guadeloupe on August 2006 (authorization number: A-971-18-01). During these experiments, the intensity of the disease was monitored daily using a clinical reaction index, assigned to the different clinical symptoms [38,40]; daily clinical scores (S2 Fig) values correspond to the sum of clinical reaction indices *per* animal and *per* day. Goat 0212 (the positive control infected with ERGvir) was treated with antibiotics (tetracycline) at day 15 to prevent animal death (S2 Fig); no animals were sacrificed in this study.

## Protein extracts preparation

Purified ER were washed in TSB (Tris-sucrose buffer: 33mM Trizma, 250mM sucrose, pH 7.4) and centrifuged at 20,000 x g, for 30 minutes at 4°C. Pellets were snap frozen in liquid nitrogen and then dissolved in DIGE solubilization buffer (CHAPS (4%, w/v), urea (7M), thiourea (2M), Trizma (30mM)) using a vortex. Four cycles of freeze and thaw in liquid nitrogen were performed and samples were sonicated (using a microtip) on ice for 30 seconds (10% amplitude, 1 second ON, 0.5 second OFF) (Branson, USA). Samples were then centrifuged at 20,000 x g for 30 min at 4°C to pellet contaminants. Protein in the supernatant was then quantified using a 2D Quant kit (GE Healthcare, Sweden).

## 1DE-nanoLC experiments

**1D-SDS-PAGE.** 20 µg of ERGvir and ERGatt protein samples were loaded on pre-cast NuPage 4–12% Novex Tris-Glycine Gels, as described elsewhere [29]. Briefly, samples were solubilized in Laemmli buffer, under reducing conditions, and subjected to electrophoresis for 35 minutes in NuPage MES Buffer (all from Invitrogen, UK). Each gel included pre-stained molecular mass markers (SeeBlue 2, Invitrogen, UK). Proteins were stained with Coomassie/ Colloidal Brilliant Blue R-250 and 14 gel bands were cut across the three ERGvir and three ERGatt biological replicates lanes (S3 Fig).

**In-gel digestion.** Each band was excised into smaller pieces, placed into a 1.5mL microcentrifuge tube and washed with 100µL of MilliQ water twice for 15 min. Gel pieces were washed three times each with 50% acetonitrile and water and the samples were reduced at 56°C for 45 min with 10mM DTT in 100mM NH₄HCO₃ and alkylated for 30 min at RT with 55mM iodoacetamide in 100mM NH4HCO3. Gel pieces were then dehydrated twice with 100% acetonitrile for 15 min. Following drying in a SpeedVac, the gel pieces were mixed with 6.7µg/mL of trypsin (Promega, Madison, WI, USA) in 50mM NH4HCO3 and incubated on ice for 45 min. Tryptic digestion was carried out at 37°C overnight and stopped with 0.5% formic acid (Sigma). Tryptic peptides were extracted from the gel using 50% acetonitrile and water. Following 20 min of vortexing, the supernatant was collected and saved. The extraction was repeated once and all of the supernatants were combined. After evaporation of acetonitrile and water in a SpeedVac, the samples were dissolved in water with 0.1%AF (12µl) and loaded onto the Ther-moplate for chromatographic separation [41].

**NanoLC experiments.** Chromatographic peptide separation was performed on a Thermo EASY-nLC 1000 with a pre-column Acclaim PepMap 100 C18 (75 µm x 2 cm) used as peptrap and an Acclaim PepMap RSLC C18 (50µm x 15cm) as the chromatographic separation column, as previously described [41]. Briefly, a chromatographic gradient was established using mixed volumes of 0.1% formic acid in water (buffer A) and 0.1% formic acid in ACN (buffer B, all LC-MS grade, from MERCK); the peptides were eluted in 5–40% buffer A for 40min according to their hydrophilic/hydrophobic properties. Peptide fractions were spotted onto MALDI plates with alpha-Cyano-4-hydroxycinnamic acid (5mg/mL, Sigma) at a constant rate
of 2ml/min using a micro-spotter (Sunchrom, Germany). Maldi plates were then analyzed on a 4800 MALDI TOF/TOF instrument (Applied Biosystems, MA, USA) with 4000 series explorer v3.5 software.

## Conventional 2DE

**Isoelectric focussing (IEF).** Protein separation of three biological replicates *per* strain was performed using 600µg of protein. This amount of protein was diluted in DIGE buffer (as described above and supplemented with IPG buffer pH 3–10 NL (0.8%, v/v) and 60mM DTT) to a final volume of 450µl. IEF was then performed using the IPGphor system (using cup loading and manifold) and 24cm immobiline drystrips with a non-linear pH gradient from 3 to 10 (all from GE Healthcare, Sweden). The protocol consisted of a sequence of 7 steps as follows: rehydration of the strips (50µA/strip at 20°C) was carried out for 12 h at 30V, followed by a step-and-hold running condition at 100Vh (3h), step from 300V to 600Vhr, step from 500V to 500Vhr, gradient at 3500V (2h), step from 3500V to 7000Vhr, gradient at 10000Vh (3h), and a final step at 10000Vh (5h), with a total of 77,7kVh.

**SDS-PAGE.** After IEF, the IPG strips were equilibrated: first the samples were reduced in equilibration buffer (50mM Tris pH 8.8, 6M urea, 30% (v/v) glycerol and 2% (w/v) SDS) supplemented with DTT (10mg/ml), followed by alkylation in equilibration buffer supplemented with iodocetamide (25 mg/ml). Each equilibration step lasted 15 min under slow agitation at room temperature [42]. Electrophoresis was then performed as previously described [29]. The IPG strips were then embedded in a precast gel (GE Healthcare, Sweden) and sealed into place using 0.5% (w/v) agarose sealing solution. The SDS-PAGE was performed using an Ettan six DALT system (GE Healthcare, Sweden) with a discontinuous buffer system of SDS electrophoresis buffer (25 mM Tris-HCl, pH 8.3, 192 mM glycine, 0.1% (w/v) SDS) in the bottom chamber and SDS electrophoresis buffer (50 mM Tris-HC pH 8.3, 384 mM glycine, 0.2% (w/v) SDS) in the top chamber overnight at 150 V, 12 mA/gel, 2 W/gel. Gels were stained using Colloidal Coomassie Blue, according to methodology described by Neuhoff et al. [43]. Briefly, gels were stained for 48h and subsequently washed three times in double distilled water. Gels were stored at 4°C in a 20% (w/v) ammonium sulphate solution until image acquiring and spot excision. Digital images of the gels were acquired using a laser-based scanner FLA-5100 (Fuji Inc., Japan). Three biological replicates were used per strain.

## 2D-DIGE

Each protein sample (50 µg) was labeled with 400 pmol of Cy3 or Cy5, and Cy2 was used as an internal calibrator. After incubating on ice for 30 min in the dark, the labelling reaction was stopped with 10 mM lysine. For each gel, Cy3- and Cy5-labeled proteins were mixed with 130µL rehydration buffer (7M urea, 2M thiourea, 2% (w/v) CHAPS, 130mM DTT, 2% IPG buffer pH 3–10). The labeled protein mixture was applied to Immobiline DryStrip strips (24cm, pH 3-10NL; GE Healthcare). Isoelectric focusing (IEF) and SDS-PAGE was performed as described above. Gel images were acquired on a laser-based scanner FLA-5100 (Fuji Inc., Japan) using 532nm and 635nm excitation lasers (DGR1double filter) for Cy3 and Cy5 respectively, and 473nm excitation laser (LPB filter) for Cy2 under Image Reader FLA 500 version 1.0 (FujiFilm). The gels were scanned using low-fluorescence glass plates at a resolution of 100µm. Four biological replicates per strain were used.

## Image analysis

Coomassie and DIGE images were analyzed using the Progenesis SameSpots v3.0 software (Non-linear Dynamics, Newcastle, UK). First, images were aligned. Next, prominent spots
were used to manually assign vectors to digitized images within each gel and then the auto-
matic vector tool was used to add additional vectors, which were manually revised and edited
for correction if necessary. These vectors were used to warp and align gel images with a refer-
ence image of one internal standard across and within each gel. After automatic spot detection,
spots were manually revised with editing tools for correct detection.

For DIGE analysis, gel groups were established according to the experimental design and spot-normalized volume was used to select statistically significant (fold-change, ANOVA, false discovery rate and power value) differentiated spots between ERG strains analyzed in the experiment. A value of 1.5-fold increase or decrease was used as a cut-off and statistically significant differences in spot intensities were identified using a Student's t-test with *p*<0.05 (with visual inspection of the results).

A protein spot list was generated and the spots with significant changes were excised from the Coomassie gels. Excision of protein spots was done as follows: with a set of paper reference circles attached to each side of the glass plate, the ordinance information for each protein of interest was translated and transferred to an automatic spot picker (GE Ettan Spot Handling Work Station, Sweden) through the pick list. Spots then were excised by the picker and transferred to a 96-well collecting plate containing 120µL of a 10% ethanol solution. Plates with excised spots were stored at -20°C until further use.

In-gel digestion for 2DE-spots

The proteins of interest were digested into their component peptides with trypsin and eluted
from the gel plugs, as described above. Briefly, spots were excised, destained, reduced with
DTT, alkylated with iodoacetamide, and dried in a speedvac. Gel pieces were rehydrated with
digestion buffer (50 mM NH4HCO3) containing trypsin (6.7 ng/µl) (Promega, Madison, WI,
USA) and incubated overnight at 37°C. The buffered peptides were acidified with formic acid,
desalted and concentrated using homemade reversed phase microcolumns (POROS R2,
Applied Biosystems, Foster City, CA, USA). The peptides were eluted onto a MALDI plate
using a matrix solution that contained 5 mg/ml α-cyano-4- hydroxycinnamic acid dissolved in
50% (v/v) ACN/0.1% (v/v) formic acid.

Protein identification by MALDI-TOF/TOF

Protein identification was done on a 4800 MALDI-TOF/TOF instrument (Applied Biosystems, Framingham, MA, USA) with 4000 series explorer v3.5 software in both MS and MS/MS mode. Each MS spectrum was obtained in a result independent acquisition mode with a total of 800 laser shots per spectrum, with internal calibration using Angiotensin II (1046.2 Da), Angiotensin I (1296.5 Da), Neurotensin (1672.9 Da), ACTH [1–17] (2093.5 Da), ACTH [18–39] (2465.7 Da) (PepMix1, from Laserbio labs). Fifteen s/n best precursors from each MS spectrum were selected for MS/MS analysis, starting with the most intense peak and ending with the least intense peak. MS/MS analyses were performed using CID (Collision Induced Dissociation) assisted with air, using a collision energy of 1 kV and a gas pressure of 1 × 10⁻⁶¹ Torr. A total of 1400 laser shots were collected for each MS/MS spectrum with the peak detection minimum signal to noise (S/N) setting at 20.

The MS/MS spectra obtained for the samples processed by 1DE-nanoLC were searched
against the database search using PEAKS search engine tool (PEAKS Studio 5.3; Bioinformatics
Solutions Inc., Waterloo, ON, Canada) [34], combining 3 search engines MASCOT, X!Tandem
and Peaks DB (S1, S2 and S3 Tables). The search parameters for the MS/MS spectra are pre-
sented in S1 File.
For protein identification in 2D gels, Mascot Generic Files combining MS and MS/MS spectra were automatically created and used to interrogate a non-redundant protein database using a local version of Mascot v2.2 from Matrix Science through the Global Protein Server (GPS) v3.6 (Applied Biosystems). The search parameters for the MS/MS spectra were described in S1 File. MDM [41] was used to organize Mascot files into Tables (S4 and S5 Tables).

## Protein functional annotation

The identified proteins were classified according to Collins et al [24], Frutos et al ([25], Moumene et al [19], Meyer et al [16]), the Clusters of Orthologous genes (COG) functional categories [44,45], the online pathway tools on the Kyoto Encyclopedia of Genes and Genomes web site (KEGG, (http://www.genome.jp/kegg-bin/show_pathway?org_name = erg&mapno=01100&mapscale=0.35&show_description = hide&show_module_list) [46] and Uniprot database (http://www.uniprot.org/proteomes/UP000000533). For proteins with no assigned function, homology searches were performed using the BlastP program against all non-redundant protein sequences deposited in the NCBI database http://blast.ncbi.nlm.nih.gov/Blast.cgi and the Uniprot database website. Protein hits were considered to be significant when their e-value was below $10^{-5}$.

## *In silico* PTM analysis

An adapted version of the QuickMod software [47] was used to search for unidentified modifications in the MS/MS data of MAP1 (ERGA_CDS_09160) and the Porin_05140 (ERGA_CDS_05140) protein spots. First, all MAP1 or Porin_05140 spectra were clustered in order to increase the spectrum quality by forming consensus spectra. Then we searched all MS/MS spectra against all consensus spectra in Open Modification Search (OMS) mode [48]. This approach relies on the alignment of modified query spectra to database spectra, where some fragment peaks are shifted in accordance with the mass of the modification(s) [49].

## Results

### Validation of the attenuated status of ERGatt strain

ERGvir has previously been shown to be highly virulent *in vivo*, leading to animal death, unless it is treated with antibiotics (tetracycline) [1,38–40]. Its growth kinetics *in vitro* have also been previously described [37], presenting a complete life cycle within 120 hours post-infection (hpi). Herein, we showed for the first time that ERG virulence was attenuated after 230 passages *in vitro*. Indeed, naive goats injected with high doses of ERGp230 (goats 0541 and 0614) survived the infection with low clinical signs (only mild hyperthermia occurred during 6 days). The goat infected with a lethal dose of ERGvir (goat 0212) presented strong clinical symptoms (high fever and neurological symptoms associated in general with fatality) and was treated with antibiotics to avoid animal death (S2 Fig). Another independent experiment with 4 goats infected with ERGatt resulted in similar results, with low clinical symptoms for naive goats infected with ERGatt (data not shown). *In vitro*, we observed that ERGatt presented a shorter life cycle of 96 hpi instead of 120hpi for ERGvir (S1A Fig). Interestingly, we could use a high multiplicity of infection for ERGatt sub-culturing without any cytotoxic effect on host cells (data not shown). On the contrary, for ERGvir, the inocula should not be above 400 ER/host cell (1/20 of the starting volume) as it can lead to cell death at the time of infection [50]. Moreover, a higher number of morula *per* host cell is also observed in ERGatt-infected cells (S1B Fig).
## General features of the *E.ruminantium* Gardel proteome

Whole proteome analysis of ERGvir and ERGatt strains using 1DE-nanoLC-MALDI-TOF/TOF and Peaks software analysis resulted in the identification of 341 and 364 non-redundant proteins, respectively. Of these, 292 proteins (31% (292/950) of the CDS in ERG genome [25] (S1 Table, Fig 1) were found to be common between ERGvir and ERGatt, which corresponds to approximately 80% of the identified ERG proteome. This low variability in the identified proteins confirms the few observed differences in the SDS-PAGE migration pattern (S3 Fig). ERGvir and ERGatt expressed 49 and 72 strain-specific proteins respectively (Fig 1, S2 and S3 Tables).

To complement the data from 1DE-nanoLC, we also used conventional bidimensional electrophoresis (2DE). S4 Fig is a representative example of a 2DE gel obtained for the three biological replicates used within this study. It shows that proteins exhibit isoelectric points (pI) ranging from 4 to 10 with most of them having an acidic pI and molecular mass (MM) range from 10 to 157 kDa (S4 Fig). Of all spots analyzed per gel, only about 50% were identified as ER proteins (S4 Table), the others being proteins of host cell origin (data not shown). The 2DE data were generally in agreement with the data obtained by 1DE-nanoLC for both strains, although there were some proteins (such as AnkC, PurE, IhfA, Def, PdxJ, ThiC, RibH, Dut, Rsfs, PpnK, PdhA, GatC and Rho) that were exclusively detected by 2DE and DIGE (S4 and S5 Tables). The sequence coverage of the identified proteins ranged from 1.1% (spot 506, protein PutA) to 78.5% (spot 2468, protein TsAA) (S4 and S5 Tables). Interestingly, we also observed that 71 spots (approx 15% (71/483 spots) of the whole conventional 2DE proteome) corresponded to non-redundant ERG proteins while all the other proteins were found in multiple spots (from 2 to 36); this suggests that ER single genes resulted in several protein molecules

<figure><img src="image_5.png" /><figcaption>**Fig 1.** Venn diagram representing the number of non-redundant proteins identified in ERGvir and ERGatt and in both (core proteome) by 1DE-nanoLC-MALDI-TOF/TOF analysis. Three biological replicates per strain were independently used and the results analyzed with Peaks software (using using simultaneously Peaks DB, MASCOT and X!Tandem algorithms).</figcaption></figure>

doi:10.1371/journal.pone.0145328.g001
("proteoforms") due to genetic variation or post-translational modifications [51] (Table 1). MAP1 (ERGA_CDS_09160, with 36 spots) and Porin_05140 (ERGA_CDS_05140, with 27 spots) were found to be the proteins with the highest number of proteoforms in the ERG proteome map. In an attempt to identify the PTMs associated to each protein spot, we used an adapted version of the QuickMod software [47]. The analyses revealed many modification mass shifts of +14Da, +16Da, +28Da, +48Da, +91Da and +283Da between different MAP1 spots, as well as +14Da, +15Da, +16Da, +91Da and +283Da between different Porin_5140 spots. These mass shifts could correspond to methylation (+14Da, +28Da), oxidation (+16Da) and GlcNAc (+283Da) or sequence variations. The best OMS scores revealed that some proteoforms of MAP1 (spots 2372, 3464 and 2581, S4 Fig) and Porin (spots 1765, 1815, 1885 and 1814, S4 Fig) correspond to N-glycosylated proteins.

## Analysis of the functional distribution of the proteome

The core proteome of ERG consists of 292 proteins common to both strains, as detected by 1DE-nanoLC whose functional distribution is represented in Fig 2 (blue color code). As predicted from the ER genome [24,25] and previous studies of ER proteomics and transcriptomics [18], ER has a high number of uncharacterized proteins (30% of CDS). We identified bacterial effectors homologous to AnkA, B and C, Ats-1, ApxR and a Bax-1 related protein by manual reannotation of uncharacterized proteins. Proteins exclusively expressed in ERGvir or ERGatt ERG reflect the proteome diversity and the proteins that may be involved in the transition from virulence to attenuation. We observed that ERGvir has a higher number of proteins related to the MAP-1 family proteins but also a higher number related to central intermediary metabolism, regulatory function, the ATP-synthase complex, chaperones, pyruvate dehydrogenase, the TCA cycle and accurate chromosome replication (Fig 2, S2 Table). ERGatt expresses a higher number of proteins involved in the biosynthesis of co-factors, cell division,

<table><caption>**Table 1. Maximum number of spots per protein (proteoforms) detected by conventional 2DE-MALDI-TOF/TOF (N = 6).**</caption><thead><tr><td>Protein name</td><td>Maximum nb of proteins spots detected in 2DE gels (N = 6)</td></tr></thead><tbody><tr><td>Map1</td><td>36</td></tr><tr><td>Porin_05140</td><td>27</td></tr><tr><td>Tuf1/Tuf2</td><td>15</td></tr><tr><td>GroES</td><td>12</td></tr><tr><td>TsaA, DnaK, ClpB, GroEL, Q5FGU8</td><td>9</td></tr><tr><td>ArgD, Pal, GlyA</td><td>8</td></tr><tr><td>AtpA, Q5FH79</td><td>7</td></tr><tr><td>RpoC, Bcp, FtsZ</td><td>6</td></tr><tr><td>Q5FGD3, VirB9</td><td>5</td></tr><tr><td>DapE, ClpX, HtpG, Ssb, IscS, Q5FFV7, PepA, SdhB, Q5FGB7, TolC</td><td>4</td></tr><tr><td>RibB, GrpE, Tme, SucC, Tig, RpoD, ArgG, AtpD, SodB, FolP, FolK, PpdK, SdhA, SecB, PurD, Q5FGI8, Ndk, Map1-1, Map1+1, Q5FGQ4, Q5FGQ2, Def, RibB, RpoA, SucB</td><td>3</td></tr><tr><td>VirB11, GlnA, PdhB, SucD, FusA, RplK, RplA, RplL, FabF, DapA, Q5FHE9, HupB, GatA, Mdh, HscA, ArgB, NusA, FumC, FolD, RpsF, Trx, GltX1, Rho, Map1-14, ElbB, Q5FF93, Efp, GyrB, Icd, Lpd, Pnp, PurH, PutA</td><td>2</td></tr></tbody></table>

doi:10.1371/journal.pone.0145328.t001
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Functional distribution of the identified proteins in ERGvir (green) and ERGatt (red) by 1DE-nanoLC-MALDI-TOF/TOF and after Peaks software analysis (using simultaneously Peaks DB, MASCOT and X!Tandem algorithms). The proteins constituting the core proteome are depicted in blue. The number of identified proteins associated with each COG functional category is shown in the X axis (total number) and in the graph bars (number per strain, n = 3).</figcaption></figure>

doi:10.1371/journal.pone.0145328.g002

metabolism (lipid, fatty acid, aminoacids), membrane-associated proteins (including transporters), DNA-RNA-protein synthesis and processing (Fig 2, S3 Table).

Alterations in protein abundance between virulent and attenuated ERG strains

To quantitatively assess proteins differentially expressed between ERGvir and ERGatt strains we used a DIGE strategy. Image analysis with Samespots software, revealed that the expression levels of 117 proteins were found to be significantly differentially expressed between ERGvir and ERGatt (p<0.05), 67 being overexpressed in ERGvir and 48 in ERGatt (S5 Fig and S5 Table). In ERGvir, some of the proteins with higher fold-change are known virulence factors in other bacteria, or related to cell redox homeostasis (such as AnkA and one proteoform of Lpd) (Fig 3A). Interestingly, amongst the proteins overexpressed by ERGatt, there are 4 proteoforms of the Porin_05140 (Fig 3B).
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig 3. Proteins overexpressed (fold change ≥ 3) in ERGvir (A) and ERGatt (B) according to DIGE analyses and MALDI-TOF/TOF.**</figcaption></figure>

doi:10.1371/journal.pone.0145328.g003

**Discussion**

ER was first identified more than 90 years ago but little data on its biology and pathogenesis are currently available. The lack of genetic manipulation techniques for ER has hampered an understanding of the role of numerous uncharacterized genes and the comprehension of virulence/attenuation mechanisms.
*In vitro* and *in vivo* infection assays revealed that the ER virulent and attenuated Gardel strains behave differently: the virulent strain is highly pathogenic to goats and can be toxic to host cells *in vitro* immediately after infection [50] while the attenuated strain undergoes a shorter life cycle and does not induce death in goats. The aim of our work was to identify virulence and attenuation-associated proteins and find proteins and/or biological processes that may have impacted the biology of ER during its long term-passaging *in vitro*, contributing to the attenuation process (Fig 4).

## General features of *Ehrlichia ruminantium* Gardel proteomes

Despite the *in vivo* and *in vitro* differences, ERGvir and ERGatt share 80% of the identified proteins, constituting a core proteome (Figs 1 and 2, S1 Table). This is in agreement with other

<figure><img src="image_9.png" /><figcaption>**Fig 4.** Schematic overview of metabolic pathways and membrane proteins found in ERG, based on genomic (KEGG database) and proteomic information obtained in this work. Nodes correspond to substrates and edges to enzymatic reactions and the 11 major metabolic pathways are color-coded (for example, light orange for amino acid metabolism). The proteins that were identified on both ERGvir and ERGatt are indicated in dark blue, the proteins and pathways found exclusively in ERGatt are in red and those detected only in ERGvir are highlighted in green. The dashed lines correspond to metabolic pathways found in both strains.</figcaption></figure>

doi:10.1371/journal.pone.0145328.g004
proteome studies comparing bacterial strains, that have regularly found the core proteome to consist of 70 to 90% of the identified proteins [29,52].

The ERG core proteome includes proteins from several functional categories, namely proteins involved in DNA-RNA-protein processing; aminoacid/carbohydrates biosynthesis, metabolism and transport; energy production; intracellular trafficking and secretion and proteins involved in homeostasis control (Fig 2). This global set of proteins is essential for the biology of ERG. These types of proteins were also shown to be crucial for other *Rickettsiales* species such as *Anaplasma phagocytophilum*, *Anaplasma marginale* and *Ehrlichia chaffensis* [18,53]. Proteins with uncharacterized function were one of the major groups of protein detected due to the lack of homology with other organisms. As these proteins are within the core proteome they must also be essential for ER biology.

Whole comparative proteomic profiling between the two Gardel strains also allowed us to detect strain-specific proteins, which might be related to virulent and attenuated status. For instance, ERGvir expresses the well conserved ParA/ParB system [54], which is known to guarantee an accurate partitioning of chromosomes between bacterial daughter cells prior to cell division. The ParA protein was not detected in ERGatt. We thus suggest that long-term passaging of ERGvir *in vitro* could have resulted in the loss of ParA protein affecting chromosome segregation and eventually bacteria growth rate, as previously observed for *Pseudomonas aeruginosa* [55]. Additionally, ERGatt was found to have a higher number of proteins involved in cellular replication (ERGA_CDS_06380, ERGA_CDS_06680, ERGA_CDS_0882, ERGA_CDS_08930) which supports its increased growth rate *in vitro* compared to ERGvir. On the other hand, quantitative differential proteomics indicate that differences in virulence might not only be associated with strain-specific proteins (S2 and S3 Tables) but also related to differential regulation of common biochemical processes (Fig 3, S5 Table). These topics will be discussed in more detail below.

An interesting feature of the ER proteomes is the high number of protein proteoforms expressed in both strains. In a previous work performed by our group with ERGvir [29], we discovered that 25% of the proteins identified were proteoforms. Herein, after optimizing protein extract preparation and 2DE conditions we found that up to 85% of proteins were proteoforms. Moreover, our results in DIGE experiments clearly revealed that different proteoforms are differently expressed between virulent and attenuated ER strains. Both results clearly indicate the biological importance of PTMs in ER, although we do not know yet their real impact on the bacterium biology and pathogenesis. In pathogenic bacteria, post-translationally modified proteins can promote bacterial survival, replication, and evolution from the host immune system [56]. For instance, PTMs such as lipoylation, glycosylation, phosphorylation and SUMOylation can have a high impact on host immune modulation during *Ehrlichia muris* [57], *Ehrlichia chaffeensis* [58] and *Anaplasma phagocytophilum* [59] infections and, multimethylation can lead to different levels of virulence in *R. prowazekii* strains [22,60]. *In silico* analysis of MS/MS data from two major proteins (MAP1 and Porin_05140) revealed that some ER protein proteoforms have N-glycosylated moieties. This corroborates with the preliminary results obtained by Postigo and co-workers [30] regarding the glycosylation of MAP1 protein proteoforms. Preliminary assays of ER PTM mapping performed by our group revealed that more than 50% of ER proteome is composed of glycoproteins [61].

### *Ehrlichia ruminantium* basic metabolic activities

Functional studies of obligate intracellular bacteria from the *Rickettsiales* order reveal the presence of several genes/proteins involved in (i) energy production and conversion and (ii) the
transport and metabolism of nucleotides, aminoacids, inorganic ions, carbohydrate and co-
enzymes [18]. Our results summarized in Fig 4 are in general agreement with these findings.

More specifically, ERGvir and ERGatt strains express enzymes involved in the Embden-Meyerhof-Parnas pathway, aminoacid metabolism, including one Proline/Betaine transporter for aminoacids (Fig 4). All predicted enzymes of the TCA cycle [24] were also found (Fig 4). We also detected the proteins involved in a partial gluconeogenesis pathway and a complete non-oxidative pentose-phosphate pathway (Fig 4).

In the virulent strain, we detected a higher number of proteins involved in energy conversion and production. This energy could be use to fuel specific processes, namely those involved in virulence. Interestingly, a higher number of proteins related to metabolism of lipids and amino acids, protein processing and biosynthesis of co-factors was detected in ERGatt (Fig 4); this could relate to a higher growth rate of ERGatt *in vitro*. This could also suggest that ERGatt is metabolically more efficient than ERGvir and that it might not need to compete with bovine endothelial host cells for essential vitamins and nucleotides, and may even supply them to the host cells. Similar processes have been proposed to occur between the obligatory intracellular bacterium *Wigglesworthia glossinidia* and its insect host, tsetse fly [62] and for *Anaplasma phagocytophilum* and *Ehrlichia chaffeensis* [52].

Expression of ERG membrane proteins

Membrane proteins are an important group of proteins as they perform a variety of functions vital to the survival of organisms (host invasion, transport, immune response, adhesion, etc). As predicted from the ER genome annotation [24,25] and previous proteomic analyses [19,29], ER has many membrane proteins. Here, we found 21 membrane proteins in the ERG core proteome, including the members from the Major Antigenic Protein 1 (MAP1) family proteins (Figs 2 and 4). Indeed, 8 out of the sixteen MAP1 paralogs (MAP1, MAP1+1, MAP1-13, MAP1-6, MAP1-1, MAP1-11, MAP1-12, MAP1-14) were detected in both ERG strains (Fig 4). Some of these proteins (MAP1, MAP1+1, MAP1-6, MAP1-13 and MAP1-14) were previously detected in host endothelial cell cultures [29,30]. Here we report for the first time that MAP1-1 is found in ER cultivated in BAE cells, and not only in tick cell lines as previously mentioned [30,63]. It has also been previously observed that the ER *map1* gene cluster can be differentially expressed according to the microenvironment (host-tick-extracellular milieu) but also according to the ER strain [63]. Our results identified a total of 12 MAP1 related proteins in ERGvir and nine MAP1 family proteins in ERGatt (Fig 4). Additionally, we found that some of these proteins have several proteoforms: 2 for MAP1-14, 3 for MAP1+1 and MAP1-1 and 36 for MAP1 (known to be the most abundant protein in ERG [29]). Although the role of MAP1-family proteins in ER biology has not yet been established, they are known to be highly immunogenic; we suggest that they could then be used by the virulent strain as bait to confound the host immune system. From the results presented above, it is clear that they must be relevant for ER biology and further investigation into their role is necessary.

We also detected the porin ERG_CDS_05140, which appears to be the second most abundant protein in ERG after MAP1. This protein was found in both strains but four proteoforms were mainly expressed in ERGatt. Although we do not know the role of this protein in ERG biology, we suggest that this porin could contribute to the exchange of small metabolites (sugar, aminoacids, etc) between ER and its environment/host cell and thereby contribute to the increased growth rate in ERGatt. Additional studies are currently being performed by our group to expand the knowledge of the function of this outer membrane protein.

Another interesting membrane protein common to both strains is ERGA_CDS_05800,
which is homologous to the *A.phagocytophilum* invasin OmpA [64]. OmpA (outer membrane
protein A) is conserved among most Gram-negative bacteria and is an important virulence factors for several Gram-negative pathogens [64–66]. It could be thus be used by both ERG strains during host cell invasion.

In ERGvir, we detected two additional proteins that could be used by this strain to invade the host cells: the enolase (ERGA_CDS_04960, Fig 3) and the protein ERGA_CDS_00060 which is homologous to *A.phagocytophilum* hemolysin (S2 Table). Enolase is a prototypic moonlighting protein in both prokaryotes and eukaryotes [67,68] and has been recently recognized to have a significant impact in a variety of pathophysiological processes. Hemolysins are membranolytic enzymes that form pores of varying diameters in the membrane of cells. This protein could thus aid ER invasion, and it has been previously observed in other bacterial infections [69–71]. Its presence uniquely in ERGvir could also support the cytotoxic effect in host cell cultivated *in vitro* culture conditions and its invasion capability *in vivo*.

## Proteins involved in *E.ruminantium*-host interactions

In order to be able to infect and successfully colonize the host, ER must rapidly modulate host cell gene transcription and function after adhesion/invasion. In this study, several transporters and molecules involved in ER-host cell cross talk were detected in both strains, including proteins from the Sec pathway (involved in both the secretion of unfolded proteins across the cytoplasmic membrane and the insertion of membrane proteins into the cytoplasmic membrane [72]) and also those related to the Type Four Secretion System (T4SS). The role of the T4SS in the pathogenicity or parasitic lifestyle of bacterial pathogens is well established, including in other *Rickettsiales* such as *A. marginale*, *E. canis* and *E. chaffeensis* [18,73–75]. Here, we detected 8 out the 10 “building blocks” of the ER T4SS [76] (VirB10, B11, B2, B4, B6, B8, B9 and VirD4, S1 Table). Manual reannotation of proteins allowed us to identify ER bacterial effectors homologous to *A.phagocytophilum* proteins, such as T4SS-secreted proteins AnkA (previously predicted in ER by S4TE software [16] and detected in ERGvir [29]), Ank–B-C, Ats-1 and AprX proteins. We also identified a BAX inhibitor (BI)-1 like protein. Ats-1 and BI-1 like protein are known to interfere with apoptosis [77], and ApxR and Ank-related proteins regulate gene expression [78–80].

Differential expression analyses using DIGE revealed that DksA (known as a virulent factor in *Salmonella typhimurium* [28] and enterohemorrhagic *Escherichia coli* [27]) and AnkA are overexpressed in ERGvir (Fig 3A). As above mentioned, AnkA alters host cell gene expression; an increased level of this protein in ERGvir would facilitate host manipulation, interfering with cellular responses in inflammatory diseases [81]. VirB9 and B10 proteins were also found to be overexpressed in ERGvir (fold change > 3, Fig 3A). Increased amounts of these T4SS core complex building blocks could be used as a support for intracellular development and therefore contribute for increased virulence of the bacteria.

In the ERGatt proteome, we detected a Patatin like-protein (ERGA_CDS_01780, homologous to *E.chaffeensis* EchaDRAFT_0464). This intracellular cytotoxin is known to perturb membrane trafficking and modulate intracellular bacterial growth [82]; it could thus promote intracellular replication inside host cells, again contributing to the increased growth rate of ERGatt *in vitro*.

## Proteins involved in *E.ruminantium* survival in the extracellular milieu

During its period outside the host cells as an infectious elementary body, ER must interact with its surroundings and protect itself from adverse conditions. Two-Component Signal transduction systems (TCS) are known to be important for these interactions as they allow organisms (especially prokaryotes) to sense and respond to changes in many different environmental
conditions [83,84]. They typically consist of a membrane-bound histidine kinase that senses a specific environmental stimulus and a corresponding response regulator that mediates the cellular response, mostly through differential expression of target genes. Several *Rickettsiales* (including ER) are known to have the genes coding for three histidine kinases (homologs of NtrY, PleC, and CckA) and their corresponding response regulators (homologs of NtrX, PleD, and CtrA) [83]. In both virulent and attenuated ERG strains, we detected the response regulators CtrA and PleD, while the sensor protein NtrY was only identified in ERGatt (S1-S3 Tables). In comparison, the three potential pairs of TCSs were detected in *E.chaffeensis* and *A. phagocytophilum* using double immunofluorescence labeling and Western blot analysis (using polyclonal antibodies) but only during their intracellular development, suggesting that TCSs might be active mostly during the bacterial intracellular development [83].

Several chaperones (such as DnaK, DnaJ, HslU, HslV, GroEL, GroES, and HtpG proteins) and other key proteins involved in cell homeostasis/oxidative stress response (such as PepA, ClpP, ClpB, DnaK, SurE, ElbB, TsA) were also identified in both strains. Some of these proteins were previously detected in ERG [29] and in other *Rickettsiales* [53,85,86]. Interestingly, one of the overexpressed proteins in ERGvir is the protein Lpd (Fig 3A). Although it is generally considered to be involved in cell homeostasis, in *A. phagocytophilum*, this protein was found to act as an immunopathological molecule, affecting cytokine and chemokine production [87]. It is an important virulent factor in several bacteria such as *Mycobacterium tuberculosis* [88], *Mycoplasma gallisepticum* [89] and *Pseudomonas aeruginosa* [90]. As host endothelial cells are able to produce cytokines and chemokines, high Lpd expression could explain the strong clinical signs observed *in vivo* during the infection with ERGvir [39] and their eventual toxicity *in vitro* [50]. Apart from the proteins found in the core proteome, no additional known proteins related to pathogen-extracellular milieu interaction were detected specifically in ERGatt.

## From virulence to attenuation: what happens?

From all the points discussed above, we propose the following hypothesis for the conversion of ERGvir to ERGatt: *In vivo*, ERGvir needs to have the “tools” to protect itself from the external environment/immune system and rapidly infect new target cells. To avoid the dangers of the surroundings, ERGvir expresses not only a high number of chaperones and proteins to avoid oxidative stress, but also a higher number of MAP1-related proteins that induce a high but unprotective antibody response. To efficiently infect target cells, ERGvir expresses more proteins related to central intermediary metabolism to fuel the higher number of proteins related to virulence (AnkA, Hemolysin, Lpd, Enolase, VirB-proteins, etc). The high number of proteoforms would also provide an efficient strategy to perform the tasks required for infection and survival. Indeed, different proteoforms with slight difference in their PTM could “sidetrack” the immune system while others could be essential for ER pathogenesis. After successive long-term passaging *in vitro* with no major selective pressure from the immune system, ERGvir adapts itself to the less constraining *in vitro* culture conditions, and eventually suffers from mutations (eventually due to the loss of ParA protein). This adaptation could also be coupled to the overexpression of proteins related to cell division, biosynthesis of co-factors, electron transport, membrane associated proteins, ribosomal proteins, proteins associated to protein production and processing and transporters in ERGatt, all resulting in rapid bacterial replication. Interestingly, ERGatt has 4 major proteoforms of a porin that could also contribute to a more efficient metabolism. On the other hand, the disappearance or lower expression of virulence-associated proteins is also observed in ERGatt, which could explain why animals survive without major clinical signs when ERGatt is administrated.
**Conclusion**

By establishing the most complete proteome profiling of two ER strains with different levels of virulence, we helped to answer some questions raised by ER genome sequences and to highlight the role and the structure of key proteins involved in ER survival, infectivity and attenuation. This is of major importance as at the moment no genetic tools are currently available to transform ER. Our data open a window of opportunity for further *in vivo* or *in vitro* studies to investigate the role of specific proteins. By using bidimensional electrophoresis, we also highlight the importance of protein proteoforms (and PTMs) in ER biology, namely in virulence. We are currently performing PTM mapping using several ER strains with different levels of virulence. As these results and outcomes constitute a post-genome reference study on ER, we believe they could be of general relevance for the biology and the mechanistic basis of pathogenesis and attenuation phenomena in other *Rickettsiales* with high impact in human and animal health.

**Supporting Information**

S1 Fig. Representative growth kinetics of ERGvir and ERGatt obtained by (A) real time PCR targeting map-1 gene (dashed arrows represent the time of total medium exchange) and (B) reverse phase microscopy (N stands for nuclei, M for morula and EBs for elementary bodies).
(PDF)

S2 Fig. Daily clinical scores of the animals infected with ERGatt (0541 and 0614) and ERGvir (0212) strains.
(PDF)

S3 Fig. 1DE-SDS-PAGE protein migration profiles of the four biological replicates used for ERGvir and ERGatt strains within this work. The molecular marker profile is presented in the first lane.
(PDF)

S4 Fig. Conventional 2D electrophoretic map of ERGatt proteins expressed at 96hpi in BAE cells. The crude extract of EBs were separated using a non-linear pH 3–10 IPG strip in the first dimension, followed by a pre-cast 12% SDS–PAGE in the second dimension. A representative gel of 3 biological replicates assays is shown.
(PDF)

S5 Fig. Spots differentially expressed between ERGvir and ERGatt strains according to DIGE experiments (using 4 biological replicates per strain) and after image analysis using Samespots software analysis (A) and PCA results (B).
(PDF)

S1 File. Protein identification by MALDI-TOF/TOF.
(PDF)

S1 Table. Proteins common to both strains (core proteome) detected by 1DE-nanoLC–MALDI-TOF/TOF. Three biological replicates per strain were independently used and the results analyzed with Peaks software (using using simultaneously Peaks DB, MASCOT and X! Tandem algorithms). For each protein, columns denote the accession number, species, corresponding accession number in ERGardel strain, protein name, ordered locus name, gene name, function, Peaks software protein score (%), sequence coverage (%), number of peptides matched per identified protein, and the number of unique peptides.
(XLSX)
**S2 Table.** Strain-specific proteins identified in ERGvir strain by 1DE-nanoLC-MALDI-TOF/TOF. Three biological replicates per strain were independently used and the results analyzed with Peaks software (using using simultaneously Peaks DB, MASCOT and X!Tandem algorithms). For each protein, columns denote the accession number, species, corresponding accession number in ERGardel strain, protein name, ordered locus name, gene name, function, Peaks software protein score (%), sequence coverage (%), number of peptides matched per identified protein, and the number of unique peptides.
(XLSX)

**S3 Table.** List of strain-specific proteins detected in ERGatt strain by 1DE-nanoLC-MALDI-TOF/TOF. Three biological replicates per strain were independently used and the results analyzed with Peaks software (using using simultaneously Peaks DB, MASCOT and X!Tandem algorithms). For each protein, columns denote the accession number, species, corresponding accession number in ERGardel strain, protein name, ordered locus name, gene name, function, Peaks software protein score (%), sequence coverage (%), number of peptides matched per identified protein, and the number of unique peptides.
(XLSX)

**S4 Table.** ERG proteins identified by conventional 2DE-MALDI-TOF/TOF (data compiled using MDM software). Spot numbers correspond to those indicated in S3 Fig. For each protein spot, columns denote the accession number in ERGardel strain, protein name, ordered locus name, gene name, function, theoretical isoelectric point (pI) and molecular masses (MM), protein scores, sequence coverage (%) and the number of peptides matched.
(XLSX)

**S5 Table.** Proteins differentially 15344040expressed between ERGvir (A) and ERGatt (B) strains according to DIGE experiments. Spot numbers correspond to those indicated in S4 Fig. For each protein spot, columns denote the accession number in ERGardel, protein name, ordered locus name, gene name, function, protein scores, sequence coverage (%), the number of peptides matched, the average normalized volume for ERGvir and ERGatt replicates, the p-value and the fold change.
(XLSX)

**Acknowledgments**

The authors also thank Joana Martins and Jonathan Gordon for paper reviewing.

**Author Contributions**

Conceived and designed the experiments: IM MV AVC NV. Performed the experiments: IM MV EP MM NV TL. Analyzed the data: IM MV EP MM NV AVC. Contributed reagents/materials/analysis tools: IM AVC NV MM FL. Wrote the paper: IM FL MM NV TL AVC.

**References**

1. Vachiery N, Marcelino I, Martinez D, Lefrançois T. Opportunities in diagnostic and vaccine approaches to mitigate potential heartwater spreading and impact on the american mainland. In Dev Biol (Basel), Vol. 135, Switzerland, 2013. 191–200.

2. Martinez D, Maillard JC, Coisne S, Sheikboudou C, Bensaid A (1994) Protection of goats against heartwater acquired by immunisation with inactivated elementary bodies of Cowdria ruminantium. *Vet Immunol Immunopathol* 41, 153–163. PMID: 8066991

3. Allsopp BA (2009) Trends in the control of heartwater. *Onderstepoort J Vet Res* 76, 81–88. PMID: 19967932
4. Zweygarth E, Josemans AI, Van Strijp MF, Lopez-Rebollar L, Van Kleef M, Allsopp BA (2005) An attenuated Ehrlichia ruminantium (Welgevonden stock) vaccine protects small ruminants against virulent heartwater challenge. *Vaccine* 23, 1695-1702. PMID: 15705474

5. Sebatjane SI, Pretorius A, Liebenberg J, Steyn H, Van Kleef M In vitro and in vivo evaluation of five low molecular weight proteins of Ehrlichia ruminantium as potential vaccine components. *Vet Immunol Immunopathol* 137, 217-225. doi: 10.1016/j.vetimm.2010.05.011 PMID: 20566221

6. Pretorius A, Liebenberg J, Louw E, Collins NE, Allsopp BA Studies of a polymorphic Ehrlichia ruminantium gene for use as a component of a recombinant vaccine against heartwater. *Vaccine* 28, 3531-3539. doi: 10.1016/j.vaccine.2010.03.017 PMID: 20338214

7. Du Plessis JL, Van Gas L, Olivier JA, Bezuidenhout JD (1989) The heterogenicity of Cowdria ruminantium stocks: cross-immunity and serology in sheep and pathogenicity to mice. *Onderstepoort J Vet Res* 56, 195-201. PMID: 2812704

8. Jongejan F, Thielemans MJ, Briere C, Uilenberg G (1991) Antigenic diversity of Cowdria ruminantium isolates determined by cross-immunity. *Res Vet Sci* 51, 24-28. PMID: 1896626

9. Martinez D, Vachieri N, Stachurski F, Kandassamy Y, Ralinaina M, Aprelon R et al. (2004) Nested PCR for detection and genotyping of Ehrlichia ruminantium: use in genetic diversity analysis. *Ann N Y Acad Sci* 1026, 106-113. PMID: 15604477

10. Perez JM, Martinez D, Sheikboudou C, Jongejan F, Bensaid A (1998) Characterization of variable immunodominant antigens of Cowdria ruminantium by ELISA and immunoblots. *Parasite Immunol* 20, 613-622. PMID: 9990646

11. Uilenberg G, Dobbelaere DA, de Gee AL, Koch HT (1993) Progress in research on tick-borne diseases: theileriosis and heartwater. *Vet Q* 15, 48-54. PMID: 8372422

12. Jongejan F, Uilenberg G, Franssen FF, Gueye A, Nieuwenhuijs J (1988) Antigenic differences between stocks of Cowdria ruminantium. *Res Vet Sci* 44, 186-189. PMID: 3387670

13. Jongejan F, Vogel SW, Gueye A, Uilenberg G (1993) Vaccination against heartwater using in vitro attenuated Cowdria ruminantium organisms. *Rev Elev Med Vet Pays Trop* 46, 223-227. PMID: 8134636

14. Martinez D. Ph.D. thesis, Utrecht University 1997.

15. Marcelino I, de Almeida AM, Ventosa M, Pruneau L, Meyer DF, Martinez D et al. Tick-borne diseases in cattle: applications of proteomics to develop new generation vaccines. In *J Proteomics*, Vol. 75, Netherlands, 2012. 4232-4250.

16. Meyer DF, Noroy C, Moumene A, Raffaele S, Albina E, Vachieri N. Searching algorithm for type IV secretion system effectors 1.0: a tool for predicting type IV effectors and exploring their genomic context. In *Nucleic Acids Res*, Vol. 41, England, 2013. 9218-9229.

17. Pruneau L, Emboule 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 PMID: 22098128

18. Pruneau L, Moumene A, Meyer DF, Marcelino I, Lefrancois T, Vachieri N (2014) Understanding Anaplasmataceae pathogenesis using "Omics" approaches. Front Cell Infect Microbiol 4, 86. doi: 10.3389/fcimb.2014.00086 PMID: 25072029

19. Moumene A, Marcelino I, Ventosa M, Gros O, Lefrancois T, Vachieri N et al. Proteomic Profiling of the Outer Membrane Fraction of the Obligate Intracellular Bacterial Pathogen Ehrlichia ruminantium. In PLoS One, Vol. 10, United States, 2015. e0116758.

20. Pinto P, Klein C, Zaha A, Ferreira H (2009) Comparative proteomic analysis of pathogenic and non-pathogenic strains from the swine pathogen Mycoplasma hyopneumoniae. *Proteome Science* 7, 45. doi: 10.1186/1477-5956-7-45 PMID: 20025764

21. Dettman JR, Rodrigue N, Kassen R. Genome-wide patterns of recombination in the opportunistic human pathogen Pseudomonas aeruginosa. In *Genome Biol Evol*, Vol. 7, England, 2014. 18-34.

22. Bechah Y, El Karkouri K, Mediannikov O, Leroy Q, Pelletier N, Robert C Genomic, proteomic, and transcriptomic analysis of virulent and avirulent Rickettsia prowazekii reveals its adaptive mutation capabilities. In *Genome Res*, Vol. 20, United States, 2010. 655-663.

23. Frutos R, Viari A, Ferraz C, Bensaid A, Morgat A, Boyer F et al. (2006) Comparative Genomics of Three Strains of Ehrlichia ruminantium: A Review. *Ann N Y Acad Sci* 1081, 417-433. PMID: 17135545

24. Collins NE, Liebenberg J, de Villiers EP, Brayton KA, 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. PMID: 15637156

25. Frutos R, Viari A, Ferraz C, Morgat A, Eychenie S, Kandassamy Y et al. (2006) Comparative Genomic Analysis of Three Strains of Ehrlichia ruminantium Reveals an Active Process of Genome Size Plasticity. *J Bacteriol* 188, 2533-2542. PMID: 16547041
26. Fournier PE, El Karkouri K, Leroy Q, Robert C, Giumelli B, Renesto P et al. Analysis of the Rickettsia africana genome reveals that virulence acquisition in Rickettsia species may be explained by genome reduction. In *BMC Genomics*, Vol. 10, England, 2009. 166.

27. Nakanishi N, Abe H, Ogura Y, Hayashi T, Tashiro K, Kuhara S et al. ppGpp with DksA controls gene expression in the locus of enterocyte effacement (LEE) pathogenicity island of enterohaemorrhagic Escherichia coli through activation of two virulence regulatory genes. In *Mol Microbiol*, Vol. 61, England, 2006. 194–205.

28. Webb C, Moreno M, Wilmes-Riesenberg M, Curtiss R 3rd, Foster JW. Effects of DksA and ClpP protease on sigma S production and virulence in Salmonella typhimurium. In *Mol Microbiol*, Vol. 34, England, 1999. 112–123.

29. Marcelino I, de Almeida AM, Brito C, Meyer DF, Barreto M, Sheikboudou C et al. (2012) Proteomic analyses of Ehrlichia ruminantium highlight differential expression of MAP1-family proteins. *Vet Microbiol* 156, 305–314. doi: 10.1016/j.vetmic.2011.11.022 PMID: 22204792

30. Postigo M, Taoufik A, Bell-Sakyi L, Bekker CP, de Vries E, Morrison WI et al. (2008) Host cell-specific protein expression in vitro in Ehrlichia ruminantium. *Vet Microbiol* 128, 136–147. PMID: 18006251

31. Ge H, Chuang YY, Zhao S, Temenak JJ, Ching WM (2003) Genomic studies of Rickettsia prowazekii virulent and avirulent strains. *Ann N Y Acad Sci* 990, 671–677. PMID: 12860705

32. Chao CC, Chelius D, Zhang T, Mutumanje E, Ching WM. Insight into the virulence of Rickettsia prowazekii by proteomic analysis and comparison with an avirulent strain. In *Biochim Biophys Acta*, Vol. 1774, Netherlands, 2007. 373–381.

33. Felsheim RF, Kurtti TJ, Munderloh UG (2009) Genome sequence of the endosymbiont Rickettsia peacockii and comparison with virulent Rickettsia rickettsii: identification of virulence factors. *PLoS One* 4, e8361. doi: 10.1371/journal.pone.0008361 PMID: 20027221

34. Zhang J, Xin L, Shan B, Chen W, Xie M, Yuen D et al. PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identification. In *Mol Cell Proteomics*, Vol. 11, United States, 2012. M111 010587.

35. Dyrlund TF, Poulsen ET, Scavenius C, Sanggaard KW, Enghild JJ (2012) MS Data Miner: a web-based software tool to analyze, compare, and share mass spectrometry protein identifications. *Proteomics* 12, 2792–2796. doi: 10.1002/pmic.201200109 PMID: 22833312

36. Emboule L, Daigle F, Meyer DF, Mari B, Pinarello V, Sheikboudou C et al. (2009) Innovative approach for transcriptomic analysis of obligate intracellular pathogen: selective capture of transcribed sequences of Ehrlichia ruminantium. *BMC Mol Biol* 10, 111. doi: 10.1186/1471-2199-10-111 PMID: 20034374

37. Marcelino I, Verissimo C, Sousa MF, Carrondo MJ, Alves PM (2005) Characterization of Ehrlichia ruminantium replication and release kinetics in endothelial cell cultures. *Vet Microbiol* 110, 87–96. PMID: 16139967

38. Vachieri N, Lefrancois T, Esteves I, Molia S, Sheikboudou C, Kandassamy Y et al. (2006) Optimisation of the inactivated vaccine dose against heartwater and in vitro quantification of Ehrlichia ruminantium challenge material. *Vaccine* 24, 4747–4756. PMID: 16621174

39. Marcelino I, Vachieri N, Amaral AI, Roldao A, Lefrancois T, Carrondo MJ et al. (2007) Effect of the purification process and the storage conditions on the efficacy of an inactivated vaccine against heartwater. *Vaccine* 25, 4903–4913. PMID: 17531356

40. Marcelino I, Lefrancois T, Martinez D, Giraud-Girard K, Aprelon R, Mandonnet N et al. (2015) A user-friendly and scalable process to prepare a ready-to-use inactivated vaccine: The example of heartwater in ruminants under tropical conditions. *Vaccine* 33, 678–685. doi: 10.1016/j.vaccine.2014.11.059 PMID: 25514207

41. Moumene A, Marcelino I, Ventosa M, Gros O, Lefrancois T, Vachieri N et al. (2014) OM proteome.

42. Gorg A, Obermaier C, Boguth G, Harder A, Scheibe B, Wildgruber R et al. (2000) The current state of two-dimensional electrophoresis with immobilized pH gradients. *Electrophoresis* 21, 1037–1053. PMID: 10786879

43. Neuhoff V, Arold N, Taube D, Ehrhardt W (1988) Improved staining of proteins in polyacrylamide gels including isoelectric focusing gels with clear background at nanogram sensitivity using Coomassie Brilliant Blue G-250 and R-250. *Electrophoresis* 9, 255–262. PMID: 2466658

44. Tatusov RL, Galperin MY, Natale DA, Koonin EV (2000) The COG database: a tool for genome-scale analysis of protein functions and evolution. *Nucleic Acids Res* 28, 33–36. PMID: 10592175

45. Tatusov RL, Natale DA, Garkavtsev IV, Tatusova TA, Shankavaram UT, Rao BS et al. (2001) The COG database: new developments in phylogenetic classification of proteins from complete genomes. *Nucleic Acids Res* 29, 22–28. PMID: 11125040
46. Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M (1999) KEGG: Kyoto Encyclopedia of Genes and Genomes. *Nucleic Acids Res* 27, 29–34. PMID: 9847135

47. Ahrne E, Nikitin F, Lisacek F, Muller M (2011) QuickMod: A tool for open modification spectrum library searches. *J Proteome Res* 10, 2913–2921. doi: 10.1021/pr200152g PMID: 21500769

48. Ahrne E, Muller M, Lisacek F (2010) Unrestricted identification of modified proteins using MS/MS. *Proteomics* 10, 671–686. doi: 10.1002/pmic.200900502 PMID: 20029840

49. Pevzner PA, Dancik V, Tang CL (2000) Mutation-tolerant protein identification by mass spectrometry. *J Comput Biol* 7, 777–787. PMID: 11382361

50. Marcelino I, Sousa MF, Verissimo C, Cunha AE, Carrondo MJ, Alves PM (2006) Process development for the mass production of Ehrlichia ruminantium. *Vaccine* 24, 1716–1725. PMID: 16257480

51. Smith LM, Kelleher NL. Proteoform: a single term describing protein complexity. In *Nat Methods*, Vol. 10, United States, 2013. 186–187.

52. Lin M, Kikuchi T, Brewer HM, Norbeck AD, Rikihisa Y (2011) Global Proteomic Analysis of Two Tick-borne Emerging Zoonotic Agents: Anaplasma phagocytophilum and Ehrlichia chaffeensis. *Front. Microbio*. 2, doi: 10.3389/fmicb.2011.00024

53. Lin M, Kikuchi T, Brewer HM, Norbeck AD, Rikihisa Y (2011) Global proteomic analysis of two tick-borne emerging zoonotic agents: anaplasma phagocytophilum and ehrlichia chaffeensis. *Front Microbiol* 2, 24. doi: 10.3389/fmicb.2011.00024 PMID: 21687416

54. Mierzejewska J, Jagura-Burdzy G. Prokaryotic ParA-ParB-parS system links bacterial chromosome segregation with the cell cycle. In *Plasmid*, Vol. 67, United States, 2012. 1–14.

55. Lasocki K, Bartosik AA, Mierzejewska J, Thomas CM, Jagura-Burdzy G. Deletion of the parA (soj) homologue in Pseudomonas aeruginosa causes ParB instability and affects growth rate, chromosome segregation, and motility. In *J Bacteriol*, Vol. 189, United States, 2007. 5762–5772.

56. Cain JA, Solis N, Cordwell SJ (2014) Beyond gene expression: the impact of protein post-translational modifications in bacteria. *J Proteomics* 97, 265–286. doi: 10.1016/j.jprot.2013.08.012 PMID: 23994099

57. Thomas S, Thirumalapura N, Crossley EC, Ismail N, Walker DH. Antigenic protein modifications in Ehrlichia. In *Parasite Immunol*, Vol. 31, England, 2009. 296–303.

58. Dunphy PS, Luo T, McBride JW. Ehrlichia chaffeensis exploits host SUMOylation pathways to mediate effector-host interactions and promote intracellular survival. In *Infect Immun*, Vol. 82, United States, 2014. 4154–4168.

59. Beyer AR, Truchan HK, May LJ, Walker NJ, Borjesson DL, Carlyon JA (2014) The Anaplasma phagocytophilum effector AmpA hijacks host cell SUMOylation. *Cell Microbiol* 17, 504–519. doi: 10.1111/cmi. 12380 PMID: 25308709

60. Abeykoon A, Wang G, Chao CC, Chock PB, Gucek M, Ching WM et al. Multimethylation of Rickettsia OmpB catalyzed by lysine methyltransferases. In *J Biol Chem*, Vol. 289, United States, 2014. 7691–7701.

61. Marcelino I, Colomé-Calls N, Laires R, Lefrançois T, Bassols A, Vachiéry N et al. PTM mapping of Ehrlichia ruminantium proteins using mass spectrometry: successful development of methods for the identification of glyco- and phosphoproteins. In *Proceedings of the 5th Management Committee Meeting and 3rd Meeting of Working Groups 1, 2 & 3 of COST Action FA1002 Košice*, Slovakia 25–26 April 2013., pp112-115 Wageningen Academic Publishers, 2014.

62. Zientz E, Dandekar T, Gross R. Metabolic interdependence of obligate intracellular bacteria and their insect hosts. In *Microbiol Mol Biol Rev*, Vol. 68, United States, 2004. 745–770.

63. Bekker CP, Bell-Sakyi L, Paxton EA, Martinez D, Bensaid A, Jongejan F (2002) Transcriptional analysis of the major antigenic protein 1 multigene family of Cowdria ruminantium. *Gene* 285, 193–201. PMID: 12039046

64. Ojogun N, Kahlon A, Ragland SA, Troese MJ, Mastronunzio JE, Walker NJ et al. Anaplasma phagocytophilum outer membrane protein A interacts with sialylated glycoproteins to promote infection of mammalian host cells. In *Infect Immun*, Vol. 80, United States, 2012. 3748–3760.

65. Smith SG, Mahon V, Lambert MA, Fagan RP. A molecular Swiss army knife: OmpA structure, function and expression. In *FEMS Microbiol Lett*, Vol. 273, England, 2007. 1–11.

66. Martinez E, Cantet F, Fava L, Norville I, Bonazzi M. Identification of OmpA, a Coxiella burnetii protein involved in host cell invasion, by multi-phenotypic high-content screening. In *PLoS Pathog*, Vol. 10, United States. e1004013.

67. Diaz-Ramos A, Roig-Borrellas A, Garcia-Melero A, Lopez-Alemany R alpha-Enolase, a multifunctional protein: its role on pathophysiological situations. J Biomed Biotechnol 2012, 156795. doi: 10.1155/2012/156795 PMID: 23118496
68. Pancholi V (2001) Multifunctional alpha-enolase: its role in diseases. *Cell Mol Life Sci* 58, 902–920.
PMID: 11497239

69. Gillespie JJ, Kaur SJ, Rahman MS, Rennoll-Bankert K, Sears KT, Beier-Sexton M et al. (2014) Secretome of obligate intracellular Rickettsia. *FEMS Microbiol Rev*

70. Wang R, Zhong Y, Gu X, Yuan J, Saeed AF, Wang S (2015) The pathogenesis, detection, and prevention of Vibrio parahaemolyticus. *Front Microbiol* 6, 144. doi: 10.3389/fmicb.2015.00144 PMID: 25798132

71. Goebel W, Chakraborty T, Kreft J (1988) Bacterial hemolysins as virulence factors. *Antonie Van Leeuwenhoek* 54, 453–463. PMID: 3144241

72. Natale P, Bruser T, Driessen AJ. Sec- and Tat-mediated protein secretion across the bacterial cytoplasmic membrane—distinct translocases and mechanisms. In *Biochim Biophys Acta*, Vol. 1778, Netherlands, 2008. 1735–1756.

73. Rikihisa Y, Lin M, Niu H, Cheng Z (2009) Type IV secretion system of Anaplasma phagocytophilum and Ehrlichia chaffeensis. *Ann N Y Acad Sci* 1166, 106–111. doi: 10.1111/j.1749-6632.2009.04527.x
PMID: 19538269

74. Lopez JE, Palmer GH, Brayton KA, Dark MJ, Leach SE, Brown WC (2007) Immunogenicity of Anaplasma marginale type IV secretion system proteins in a protective outer membrane vaccine. *Infect Immun* 75, 2333–2342. PMID: 17339347

75. Felek S, Huang H, Rikihisa Y (2003) Sequence and expression analysis of virB9 of the type IV secretion system of Ehrlichia canis strains in ticks, dogs, and cultured cells. *Infect Immun* 71, 6063–6067. PMID: 14500531

76. Gillespie JJ, Brayton KA, Williams KP, Diaz MA, Brown WC, Azad AF et al. Phylogenomics reveals a diverse Rickettsiales type IV secretion system. In *Infect Immun*, Vol. 78, United States, 2010. 1809–1823.

77. Niu H, Kozjak-Pavlovic V, Rudel T, Rikihisa Y (2010) Anaplasma phagocytophilum Ats-1 is imported into host cell mitochondria and interferes with apoptosis induction. *PLoS Pathog* 6, e1000774. doi: 10.1371/journal.ppat.1000774 PMID: 20174550

78. Rikihisa Y, Lin M (2010) Anaplasma phagocytophilum and Ehrlichia chaffeensis type IV secretion and Ank proteins. *Curr Opin Microbiol* 13, 59–66. doi: 10.1016/j.mib.2009.12.008 PMID: 20053580

79. Garcia-Garcia JC, Rennoll-Bankert KE, Pelly S, Milstone AM, Dumler JS. Silencing of host cell CYBB gene expression by the nuclear effector AnkA of the intracellular pathogen Anaplasma phagocytophilum. In *Infect Immun*, Vol. 77, United States, 2009. 2385–2391.

80. Wang X, Kikuchi T, Rikihisa Y. Proteomic identification of a novel Anaplasma phagocytophilum DNA binding protein that regulates a putative transcription factor. In *J Bacteriol*, Vol. 189, United States, 2007. 4880–4886.

81. Sinclair SH, Rennoll-Bankert KE, Dumler JS (2014) Effector bottleneck: microbial reprogramming of parasitized host cell transcription by epigenetic remodeling of chromatin structure. *Front Genet* 5, 274. doi: 10.3389/fgene.2014.00274 PMID: 25177343

82. Mansueto P, Vitale G, Cascio A, Seidita A, Pepe I, Carroccio A et al. (2012) New insight into immunity and immunopathology of Rickettsial diseases. *Clin Dev Immunol* 2012, 967852. doi: 10.1155/2012/967852 PMID: 21912565

83. Cheng Z, Kumagai Y, Lin M, Zhang C, Rikihisa Y. Intra-leukocyte expression of two-component systems in Ehrlichia chaffeensis and Anaplasma phagocytophilum and effects of the histidine kinase inhibitor closantel. In *Cell Microbiol*, Vol. 8, England, 2006. 1241–1252.

84. Beceiro A, Tomas M, Bou G. Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? In *Clin Microbiol Rev*, Vol. 26, United States, 2013. 185–230.

85. Renesto P, Azza S, Dolla A, Fourquet P, Vestris G, Gorvel JP et al. (2005) Rickettsia conorii and R. prowazekii proteome analysis by 2DE-MS: a step toward functional analysis of rickettsial genomes. *Ann N Y Acad Sci* 1063, 90–93. PMID: 16481497

86. Hajem N, Weintraub A, Nimtz M, Romling U, Pahlson C (2009) A study of the antigenicity of Rickettsia helvetica proteins using two-dimensional gel electrophoresis. *Apmis* 117, 253–262. doi: 10.1111/j. 1600-0463.2009.02435.x PMID: 19338513

87. Chen G, Severo MS, Sakhon OS, Choy A, Herron MJ, Felsheim RF et al. Anaplasma phagocytophilum dihydrolipoamide dehydrogenase 1 affects host-derived immunopathology during microbial colonization. In *Infect Immun*, Vol. 80, United States, 2012. 3194–3205.

88. Venugopal A, Bryk R, Shi S, Rhee K, Rath P, Schnappinger D et al. Virulence of Mycobacterium tuberculosis depends on lipoamide dehydrogenase, a member of three multienzyme complexes. In *Cell Host Microbe*, Vol. 9, United States. 21–31.
89. Hudson P, Gorton TS, Papazisi L, Cecchini K, Frasca S Jr, Geary SJ. Identification of a virulence-associated determinant, dihydrolipoamide dehydrogenase (Ipd), in Mycoplasma gallisepticum through in vivo screening of transposon mutants. In *Infect Immun*, Vol. 74, United States, 2006. 931–939.

90. Hallstrom T, Morgelin M, Barthel D, Raguse M, Kunert A, Hoffmann R et al. Dihydrolipoamide dehydrogenase of Pseudomonas aeruginosa is a surface-exposed immune evasion protein that binds three members of the factor H family and plasminogen. In *J Immunol*, Vol. 189, United States. 4939–4950.