New species of Ehrlichia isolated from Rhipicephalus (Boophilus) microplus shows an ortholog of the E. canis major immunogenic glycoprotein gp36 with a new sequence of tandem repeats

Ehrlichia species are the etiological agents of emerging and life-threatening tick-borne human zoonoses that inflict serious and fatal infections in companion animals and livestock. The aim of this paper was to phylogeneticaly characterise a new species of Ehrlichia isolated from Rhipicephalus (Boophilus) microplus from Minas Gerais, Brazil. The agent was isolated from the hemolymph of Rhipicephalus (B.) microplus engorged females that had been collected from naturally infested cattle in a farm in the state of Minas Gerais, Brazil. This agent was then established and cultured in IDE8 tick cells. The molecular and phylogenetic analysis was based on 16S rRNA, groEL, dsb, gltA and gp36 genes. We used the maximum likelihood method to construct the phylogenetic trees. The phylogenetic trees based on 16S rRNA, groEL, dsb and gltA showed that the Ehrlichia spp isolated in this study falls in a clade separated from any previously reported Ehrlichia spp. The molecular analysis of the ortholog of gp36, the major immunoreactive glycoproteins in E. canis and ortholog of the E. chaffeensis gp47, showed a unique tandem repeat of 9 amino acids (VPAASGDAQ) when compared with those reported for E. canis, E. chaffeensis and the related mucin-like protein in E. ruminantium. Based on the molecular and phylogenetic analysis of the 16S rRNA, groEL, dsb and gltA genes we concluded that this tick-derived microorganism isolated in Brazil is a new species, named E. mineirensis (UFMG-EV), with predicted novel antigenic properties in the gp36 ortholog glycoprotein. Further studies on this new Ehrlichia spp should address questions about its transmissibility by ticks and its pathogenicity for mammalian hosts.


Background
The emergence of multiple Ehrlichia species as etiological agents of newly discovered human zoonoses and the previous recognition of these agents as causing serious disease in companion animals and livestock have intensified the interest in these pathogens. Ehrlichiae are tick-transmitted obligate intracellular gram-negative bacteria that are maintained in nature by persistent infection of mammalian hosts [1]. They are microorganisms residing within the cytoplasmic vacuoles of monocytes, granulocytes, or platelets of humans and animals. Ehrlichia species elicit illnesses with fever, headache, leukopenia, and thrombocytopenia [2].
The obligately intracellular alpha-proteobacterial genus Ehrlichia (Rickettsiales: Anaplasmataceae) is spread all over the world and are comprised of five recognized species that are tick-transmitted, with three of the five causing human ehrlichiosis (E. canis, E. chaffeensis, and E. ewingii) [3]. The agent that causes the veterinary disease heartwater (E. ruminantium) can potentially infect humans [2,4] and Ehrlichia muris has not been associated with human infection. In addition, numerous candidate entities have been reported ("E. walkerii", "E. shimanensis", "Ixodes ovatus ehrlichia", "Panola Mountain ehrlichia", etc.), all isolated from hard ticks and mainly characterized by PCR sequencing [3]. To date, only three species of the genus Ehrlichia have been reported in Brazil: E. canis, E. ewingii and E. chaffeensis [5].
Polyphasic taxonomy has been advocated to ensure well-balanced determination of taxonomic relationships [7]. Different genes have been proposed to classify ehrlichial agents. The most widely used are 16S rRNA [8,9], groESL operon [10], groEL gene [11], gltA [7], dsb [12], gp36 and gp19 [13]. The gp36 belong to the group of major immunogenic antigen in E. canis (gp36) and E. chaffeensis (gp47) and both are orthologs to the mucin-like protein in E. ruminantium. These glycoproteins have tandem repeats that contain major B-cell epitopes with carbohydrate determinants, which contribute substantially to the immunoreactivity of these proteins. Only five types of tandem repeats have been characterized [14]. Of these glycoproteins, gp36 is the most divergent gene among E. canis isolates [15]. Nevertheless, the tandem repeat is highly conserved among different isolates, changing only in the number of repeats [13] and in few amino acids among E. canis isolates [15].
Recently, we have isolated an organism from hemolymph of R. (B.) microplus engorged females which had been collected from naturally infested cattle in Brazil (unpublished data). This organism has been propagated continuously in vitro, both in a tick cell line (IDE8) and in a monocyte-macrophage cell line from a dog (DH82), and has been initially characterised as a new genotype of Ehrlichia spp (UFMG-EV strain) [16]. In the present study we report further molecular and phylogenetic analyses focusing on five genes (16S rRNA, groESL, gltA, dsb and gp36) of this new organism, from now on referred as Ehrlichia mineirensis (UFMG-EV).

Organism isolation and in vitro cultivation
Eleven R. (B.) microplus engorged females, larger than 4.5 mm in length, were collected from naturally infested calves (4 to 6 months old) from a farm in Minas Gerais, Brazil. The ticks were washed, blotted dry, and disinfected with Germekil (Johnson, Brazil.), for 30 minutes at room temperature. After several washes in sterile distilled water, the ticks were individually placed into polystyrene plates and were incubated at 27°C and relative humidity over 83%. After a 10-day incubation period hemolymph were collected to provide material for infecting IDE8 cells [17]. Each tick was held with sterile forceps, the cuticula was again sterilized, as previously described, and the leg cut with a sterile scalpel blade. The hemolymph was collected using a capillary tube to gather the draining fluid. Hemolymph from three ticks were pooled in a tube containing 200 μl of culture medium, which constitute the inoculum to infect one culture flask containing an on growing IDE8 cell monolayer.
After infection, the culture flask was monitored daily by examination of cytocentrifuge smears made from 50 μl aliquots taken from the culture suspension. Smears were fixed twice with methanol (for 10 min), stained with an 8% Giemsa solution for 30 min and examined under oil immersion at 1,000x magnification. The first infected cells were detected 28 days after culture initiation.
Maintenance of cultures was carried out with medium changes weekly. Briefly, IDE8 cells were maintained at 32°C in L-15B medium [18], supplemented with 5% heat-inactivated foetal bovine serum, 10% tryptose phosphate broth, 0.1% bovine lipoprotein concentrate (MP Biomedicals, Santa Ana, CA, USA), 100 IU/ml penicillin and 100 μg/ml streptomycin. Infected IDE8 cultures were propagated in a modified L-15B medium as outlined above, further supplemented with 0.1% NaHCO 3 and 10 mM HEPES. The pH of the medium was adjusted to 7.5 with 1 N NaOH. Infected cultures were propagated at 34°C in 25 cm 2 plastic culture flasks in 5 ml of the medium under normal atmospheric conditions.

Genomic DNA isolation
The DNeasy Blood & Tissue Kit (Qiagen Inc. Valencia, Calif.) was used for extraction of DNA from infected IDE8 cells. DNA extraction was performed according to the manufacturer's instructions. The extracted material was eluted from the columns in 100 μl of sterile double distilled H 2 O (ddH2O), and the DNA concentration and purity were determined by measuring the optical density at both 260 and 280 nm with a DNA-RNA calculator (NanoDrop W ND-1000, Peqlab, Erlangen, Germany). Ten-fold dilutions were done with the genomic DNA and separated in aliquot of 10 μl each and kept frozen until their use in a PCR reaction.

PCR
The primers used in this study are shown in (Table 1). The oligonucleotide primers used for the amplification of dsb gene and gltA gene were designed for this study using primer design software (PrimerSelect; DNAStar, USA) and information from the E. canis genome [Gen-Bank: CP000107] [19]. Two independent PCR reactions were performed for each gene. For each PCR amplification, 2 μL of extracted DNA was used as the template in a 25 μL reaction mixture containing 20pmol of each primer and 2X PCR Master Mix (Promega, USA). The reactions were conducted in an Eppendorf thermocycler (Eppendorf Mastercycler personal AG, 22331 Hamburg, Germany) according to the parameters: 2 min at 94°C followed by 40 cycles of 30 sec at 94°C, 1 min at 45°C, and 1.5 min at 72°C with a final extension step of 5 min. The PCR products were stained using an Ethidium bromide free system, 6X Orange DNA Loading Dye (Thermo Scientific, Germany) and visualized in 0.8% agarose minigels.

Cloning and sequencing
The resulting PCR products were electrophoresed on a 0.8% agarose gel. The size of the amplified fragments was checked by comparison to a DNA molecular weight marker (100-bp DNA Ladder; Promega, USA). In each case, the single amplified product of the expected size was column purified using the QIAquick PCR Purification Kit (Qiagen, USA) and then ligated into the TOPO TA Cloning Kit (Invitrogen, USA) for subsequent transformation in Escherichia coli TOP 10 Chemically Competent cells. For each gene, five individuals clones containing the cloned fragment in the TOPO vector were purified using the QIAprep Spin Miniprep Kit (Qiagen, USA) and prepared for sequencing using an ABI 3130 sequencer (Applied Biosystems, USA) and the Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA) with the M13F and M13R vector primer. Both the sense and antisense strands of each PCRamplified product were sequenced, and the sequences were then manually edited to resolve any ambiguities. A consensus sequence was obtained for each amplified PCR product by comparing both the sense and antisense sequences from the five clones.

DNA sequence analysis
To find the homology of our sequences we used the database Nucleotide collection (nr/nt) using Megablast (optimize for highly similar sequences) from the BLAST server [20]. Nucleotide sequences were aligned using BLAST [20] and protein sequences were aligned using the multiple-alignment program CLUSTALW [21]. The homology between sequences was analyzed using MegAlign, DNAStar, USA. Nucleotide sequences were translated to amino acid (aa) sequence by the ExPASy translation tool of the Swiss Institute of Bioinformatics [22].
The phylogenetic analysis was performed as follows: sequences were aligned with MUSCLE (v3.7) configured Table 1 Primers used in this study for the amplification of the 16S rRNA, groESL, gltA, dsb and gp36 genes from E. mineirensis (UFMG-EV) genomic DNA

Target
Primers* Sequence Expected size (Kb) for highest accuracy [23]. After alignment, ambiguous regions (i.e., containing gaps and/or poorly aligned) were removed with Gblocks (v0.91b) [24]. The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.0 aLRT) [25,26]. Reliability for internal branch was assessed using the bootstrapping method (100 bootstrap replicates). Graphical representation and edition of the phylogenetic tree were performed with TreeDyn (v198.3) [27]. The nomenclature used in the trees is according to Dumler et al., [19]. The same analysis of similarity and phylogenetic relationships was performed for the genes 16S rRNA, groEL, gltA and dsb with the exception that the dsb tree is unrooted and the rest are rooted.

Analysis of the glycoprotein gp36 gene and putative aa sequence
The gp36 ortholog was tested for the presence of signal peptide sequences with the computational algorithm Sig-nalP trained on gram-negative bacteria [28]. The gp36 protein sequence was evaluated for potential mucin-type O-linked glycosylation on serines and threonines with the computational algorithm NetOGlyc v3.1 [29] and for N-linked glycosylation was used the NetNGlyc 1.0 Server [30]. The Tandem Repeats Finder database [31] was used to analyze the tandem repeats. The prediction of continuous B cell epitopes was done using the B cells Epitopes Prediction Tool [32] and the 3D structure of the glycoprotein and the predicted epitopes was obtained using the algorithm contained in the ElliPro epitope modeling tool and sequences available in the ElliPro server [33]. As previously reported [14], for the convenience of sequence comparison the gp36 gene orthologs were divided into three regions: 5' end prerepeat region, a tandem repeat region, and 3' end postrepeat region.

Sequences used in this study
The sequences obtained from Ehrlichia mineirensis

Sequence analysis of 16S rRNA
In order to obtain relevant information from 16S rRNA at the species level, the primers 8 F and 1448R were used to isolate a fragment of~1.4Kb. Approximately a 1.4Kb amplicon corresponding to the expected size of targeted 16S rRNA gene fragment was obtained (data not shown). A consensus sequence of 1.384 Kb was obtained from 2 independent PCRs and five clones were sequenced. In total, our sequence had 10 changes of nucleotides when compared with E. canis [GenBank: GU810149] with two insertions and three deletions (data not shown). The percent of identities with all the members of the Ehrlichia genus are shown in the Table 2 upper triangle. Figure 1A shows the tree build using the maximum likelihood method; it shows that E. mineirensis (UFMG-EV) falls in a clade separated from all the previous reported sequences. The tree build with the neighbour joining method using the Kimura 2 parameters substitution model show identical results (data not shown). The gene 16S rRNA has a highly variable region located at the 5' end of the gene [8]. This fragment is useful in identifying Ehrlichia spp [9]. Figure 2 shows three changes in nucleotides in E. mineirensis (UFMG-EV) in comparison with E. canis and seven changes in nucleotides when compare with Ehrlichia. sp. Tibet which was isolated from R. microplus [8].

Sequence analysis of dsb
The amplicon obtained from the PCR set up with the primers dsbF2 and dsbR2 gave a band with the expected size of 0.7 Kb. A fragment of 0.683 Kb of the gene dsb was obtained and sequenced. Dsb gene sequences for available Ehrlichia spp. were aligned using clustalW. The alignment shows that dsb gene is conserved (76.4% -94.7%) within the genus (  (Table 3) The phylogenetic tree shows that E. mineirensis (UFMG-EV) dsb is separated from its homologs in other species of the Ehrlichia genus ( Figure 3).

Sequence analysis of groESL operon
The amplification with primers HS1-HS6 produced a PCR product in the expected size 1.4Kb. The nucleotide sequences of the PCR products amplified from E. mineirensis (UFMG-EV) contained a reading frame corresponding to the 26 aa carboxyl-terminus of groES, 416 aa of the amino-terminal end of groEL, and the spacer between them. The length of the nucleotide sequence of the spacer region in the sequence reported here were 95 bases. Sequence homology analyses were done for each of the nucleotide sequences and the deduced aa sequences from the partial GroES and GroEL reading frames. Nucleotide and aa sequence homologies with other members of the Ehrlichia genus are presented in Table 4. A phylogenetic tree based on multiple sequence alignment of the 1.249 Kb corresponding to groEL is presented in Figure 1B.  after sequencing and consensus analysis. The putative citrate synthase protein predicted using the E. mineirensis (UFMG-EV) gltA gene was 416 aa. Table 5 shows the nucleotide and the aa similarities with other members of the Ehrlichia genus. The gltA gene has been proposed as an alternative tool for the phylogenetic analysis of the genus Ehrlichia [7]. Using the maximum likelihood method we built a phylogenetic tree showing that E. mineirensis (UFMG-EV) falls in a clade apart from any previously reported gltA genes in the family Anaplasmataceae ( Figure 4).
Sequence analysis of the gp36 gene and the putative encoded protein sequence The gp36 based PCR products derived from the isolate reported here had a molecular size of 1000 base pair (bp) (data not shown). Subsequent cloning of the PCR amplicons followed by sequencing showed that our gene was 0.948 Kb encoding a predicted protein with 315 aa and a molecular mass of 31.51 KDa (28.89 KDa without the predicted 23-aa signal peptide). We found that the gp36 protein isolated in our study is a putative glycoprotein. The aa sequence of gp36 in our study has five potential sites of O-glycosylation and two of N-glycosylation. The O-carbohydrates were predicted to be linked to three serines (S) of the tandem repeat region at position 155, 164 and 173 and two threonines (T) present in the post-repeat region at position 286 and 289. We explored as well the possibility to find N-glycosylation on putative glycosylated asparagines (N). Two sequons of N-glycosylation (N-Xaa-T/S) at the pre-repeat region were found: NRS (at position 81) and NFS (at position 106).
Differences found in the Region I (The 5′ end pre-repeat region) Alignment of the gp36 ortholog obtained in this study revealed that our sequence was 422 nucleotides in length encoding for 141 aa ( Table 6). The nucleotide and predicted aa sequences exhibited relatively low identities, ranging from 54.9% to 91.2%, and from 38.0% to 82.0%, respectively, in comparison with related genes previously published for the gp36 orthologs in E. canis, E. chaffeensis and E. ruminantium [14] (Table 6).

Region II (the tandem repeat region)
Region II in E. mineirensis (UFMG-EV) contains 16 tandem repeats of 27 bp, each encoding nine aa. The single   Table 7). The tandem repeat of E. mineirensis (UFMG-EV) is a serine enriched area of the total protein sequence but does not contain threonine. Its glycoprotein gene shows a high C + G percent in the whole gene (42.0%) and in the tandem repeat region (52.1%).

Region III (the 3′ end post-repeat region)
The comparison of region III among the orthologs show that it is a quite variable region, presenting differences in length, nucleotide and aa sequence. It has been widely revised by [14] and [15]. Our sequence was 94-bp length, which differ from any previously reported (data not shown
We found that all epitopes were exposed on the surface of the predicted 3D structure of each protein. The superposition analysis of the epitopes 3D structure showed that they were structurally dissimilar with a root mean square deviation (rmsd) of 5-6 Å between the epitope of E. mineirensis (UFMG-EV) and others three Figure 5A-E. A linear correlation between the rmsd and % (dis) similarities among structure and sequences, respectively, is a valid interpretation for the evolution of homolog proteins [34]. Correlation for the epitopes of E. mineirensis (UFMG-EV) when compared with the other three orthologs gives an R 2 = 0.77.

Discussion
Polyphasic taxonomy has been advocated to ensure wellbalanced determinations of taxonomic relationships [7]. Different genes have been proposed to classify ehrlichial agents, however, the most widely used are 16S rRNA [8,9], groESL operon [10], groEL gene [11], gltA [7], dsb [12], gp36, and gp19 [13]. Sequence comparison of the 16S rRNA gene is recognized as one of the most powerful and precise methods for determining the phylogenetic relationships of bacteria [8,11,35]. Our results were consistent with previous phylogenetic analysis of Ehrlichia spp by using the 16S rRNA gene sequences [9,36]. In this study, our analysis of a relevant fragment of 16S rRNA sequences revealed that the novel agent found in Brazilian R. (B.) microplus ticks was closely related to E. canis [GenBank: GU810149], but was also closely related to E. chaffeensis [GenBank: AF147752] showing 98.3% and 96.9% of homology, respectively. It is worth noting that the hypervariable region 16S rRNA is well conserved in members of the same species (data not shown) and are different among members of Ehrlichia genus [8,9]. However, our hypervariable region of 16S rRNA was different when compared with other members of Ehrlichia genus.
Since the 16S rRNA gene is known to exhibit a high level of structural conservation with a low evolutionary rate, levels of sequence divergence greater than 0.5% in comparisons with nearly complete 16S rRNA gene sequences of members of the genus Ehrlichia have been considered sufficient to classify organisms as different species [8,35]. The levels of divergence of the 16S rRNA sequence between this novel Brazilian ehrlichial agent and the closest member of the Anaplasmataceae, E. canis was 1.7% in pairwise comparisons of 1384 base sequences (data not shown), and this level of difference should be sufficient to classify the novel ehrlichial agent as a new species of the genus Ehrlichia. Furthermore, the 16S rRNA phylogenetic tree constructed with a maximum likelihood method show that E. mineirensis (UFMG-EV) falls in a different clade separated from any previously reported Ehrlichia spp.
The genes groEL [11] and gltA [7] have been proposed as an alternative to 16S rRNA for the phylogenetic analysis of the Anaplasmatacaea family as they are less  conserved than 16S rRNA among the family members [7] and dsb gene has been previously used to classified members of the Ehrlichia genus [12]. It is important to note that the spacer of the groESL operon was 95 bp in E. mineirensis (UFMG-EV), which differs from the reported for E. canis, E. chaffeensis, E. ruminantium with 93, 100 and 96 bp, respectively [10]. The gp36 orthologs are a divergent gene in E. canis, E. chaffeensis and E. ruminantium due to their high evolutionary pressure [14,15]. This gene has been used to differentiate new isolates of E. canis where 16S rRNA was not well suited to discriminate between E. canis isolates [13].
In our study the level of similarity among ehrlichial gltA and dsb were lower than that of 16S rRNA and groEL gene sequences in the genus Ehrlichia. E. canis was the closest Ehrlichia species to E. mineirensis (UFMG-EV) in all the studied genes. Similar phylogenetic relationships are observed between other members of the Ehrlichia genusi.e., E. chaffeensis/E. muris, N. risticii/N. sennetsu and A. marginale/A. platys.
The architecture of gltA, groEL and dsb based phylogenetic trees were similar to that of the tree derived from the 16S rRNA gene sequences. However, the trees constructed from gltA and dsb show more divergence Table 6 Length and percent of nucleotide and aa homology of the 5' end pre-repeat region between the orthologs of gp36 in Ehrlichia mineirensis (UFMG-EV) and related genes  than that from the 16S rRNA and groEL gene. The difference of E. canis and E. mineirensis (UFMG-EV) was well established in all the four trees based on nucleotide sequences. E. mineirensis (UFMG-EV) was well defined, with higher bootstrap values in the gltA (100) and dsb (100) based trees than for those of the 16S rRNA (97) and groEL (93) based tree. Based on aa homology and genomic synteny analyses, it has been determined that the mucin-like protein of Ehrlichia ruminantium, gp36 of E. canis and gp47 of E. chaffeensis are orthologs [14]. Identity of 87.2% has been found in the pre-repeat region among geographically distant E. canis isolates [13]. The single tandem repeat was highly conserved among isolates (TEDSVSAPA) with variations in the number of repeats [13][14][15] and few conservative changes in amino acid sequences [15]. The tandem repeat genetic unit varies in length (from 27 bp -99 bp) among the different orthologs, number of repeats (from 3.4 -56) and the homology of the nucleotide and the aa sequence encoded in the repeat (Table 7). Our sequence contains a tandem repeat that shares an extremely low homology with the gp36 orthologs reported until now ranging from 22% (E. ruminantium and E. canis) to 33% (E. chaffeensis). Doyle et al. [14] describes gp36 and gp47 as glycoprotein sharing O-glycosylation predicted sites in the serines and threonines of the tandem repeat. It is noteworthy that the tandem repeat of our sequence does not contain threonine; nevertheless, we predicted three sites of O-glycosylation in the serines of the tandem repeat and two in threonines of the postrepeat region. Two N-glycosylation sites were found in our aa sequence. The analysis for N-glycosylation was done for E. ruminantium, E. canis and for E. chaffeensis ortholog sequences (data not shown) and potential sites of N-glycosylation were found as well for these sequences. Glycosylation plays a crucial role in the immunogenicity of these glycoproteins [14,15]. Deglycosylation of the gp36 tandem repeat drastically reduces its immunogenicity [14]. Both gp36 and gp47 are described as the major immunoreactive protein of E. canis and E. chaffeensis and the tandem repeats contain the major antibody epitope [14,15]. It was found that the tandem repeat of gp36 from E. mineirensis (UFMG-EV) contain the major B cell epitope previously reported for the glycoprotein orthologs. The prediction of the 3D structure of the B cell epitopes present in the tandem repeat shows a high structural divergence among the closest gp36 orthologs in E. mineirensis (UFMG-EV), E. canis and E. chaffeensis. These structural differences may explain the results obtained by Doyle et al. [14] in which neither gp36 nor gp47 reacted with heterologous antisera. The C + G content of the gp36 gene of E. mineirensis (UFMG-EV) is higher than the rest of the orthologs previously reported (data not shown). The C + G content in specific genes have been used in systematics as support for the classification of organisms [7], and it is known that recombination significantly increases the silent C + G content of a genome in a selectively neutral manner [37].
Although it is well known that Babesia bovis, B. bigemina and Anaplasma marginale are the most common etiological agents transmitted by R. (B.) microplus ticks [38], the detection of any species of Ehrlichia in R. (B.) microplus ticks has been infrequently reported. The first two reports were in China in the Guangxi Autonomous Region in 1999 [39] and Tibet in 2002 [8]; the second in Thailand in 2003 [36] and the latest one in Xiamen, China in 2011 [40]. Except the isolate from Guangxi, E. canis [39], the rest share, based on 16S rRNA, a 99.9% of homology [36,40] and differ from the ehrlichial species previously reported and classified as Ehrlichia spp strain Tibet [8]. In the present study, determined by pairwise alignment, the E. mineirensis (UFMG-EV) isolated from R. (B.) microplus shares 97% of similarity with the 16S rRNA sequences of the referred species (data not shown). This is the second report of a new Ehrlichia spp isolated from R. (B.) microplus, but the first to be reported in the American continent. The identification of E. mineirensis (UFMG-EV) in R. (B.) microplus ticks suggests a potential of infection and transmission of this agent to cattle in the area where infected ticks are present.

Conclusions
Based on the molecular and phylogenetic analysis of the genes 16S rRNA, groEL, dsb and gltA we concluded that the new microorganism isolated from the hemolymph of R. (B.) microplus is a new species of Ehrichia with new predicted antigenic properties in the gp36 glycoprotein ortholog. Complementary analysis of C + G content in the gp36 orthologs, distant of groESL spacer and hypervariable region of 16S rRNA supports the fact that E. mineirensis (UFMG-EV) is a separate phylogenetic entity.
Further studies should address the question whether R. (B.) microplus is a competent vector for this and other Ehrlichia species and whether this new organism is an emerging pathogen for cattle or an endosymbiont of R. (B.) microplus.