A new multiple-locus variable-number tandem repeat analysis reveals different clusters for Anaplasma phagocytophilum circulating in domestic and wild ruminants
Parasites & Vectors volume 7, Article number: 439 (2014)
Anaplasma phagocytophilum is a tick-borne intragranulocytic alpha-proteobacterium. It is the causative agent of tick-borne fever in ruminants, and of human granulocytic anaplasmosis in humans, two diseases which are becoming increasingly recognized in Europe and the USA. However, while several molecular typing tools have been developed over the last years, few of them are appropriate for in-depth exploration of the epidemiological cycle of this bacterium. Therefore we have developed a Multiple-Locus Variable number tandem repeat (VNTR) Analysis typing technique for A. phagocytophilum.
Five VNTRs were selected based on the HZ human-derived strain genome, and were tested on the Webster human-derived strain and on 123 DNA samples: 67 from cattle, 7 from sheep, 15 from roe deer, 4 from red deer, 1 from a reindeer, 2 from horses, 1 from a dog, and 26 from ticks.
From these samples, we obtained 84 different profiles, with a diversity index of 0.96 (0.99 for vertebrate samples, i.e. without tick samples). Our technique confirmed that A. phagocytophilum from roe deer or domestic ruminants belong to two different clusters, while A. phagocytophilum from red deer and domestic ruminants locate within the same cluster, questioning the respective roles of roe vs red deer as reservoir hosts for domestic ruminant strains in Europe. As expected, greater diversity was obtained between rather than within cattle herds.
Our technique has great potential to provide detailed information on A. phagocytophilum isolates, improving both epidemiological and phylogenic investigations, thereby helping in the development of relevant prevention and control measures.
It is now recognized that a pathogen’s ability to infect multiple, rather than just a single host species, is more the rule than the exception[1, 2]. However for each pathogenic species, it remains difficult to determine whether different animal species are randomly infected, or whether certain strains circulate in more or less distinct epidemiological cycles[3, 4]. Such knowledge is crucial in order to understand pathogenic epidemiological cycles and to formulate and implement relevant prevention and control measures. Nowadays, improved pathogen detection using molecular tools has opened up new opportunities to address these issues. Molecular tools are invaluable for studying the circulation of isolates between hosts, as well as the respective roles of the different animal species involved.
Many pathogens transmitted by Ixodes ricinus typically exploit multiple different host species, as vector ticks can feed on diverse numbers of vertebrates. This multiplicity of hosts is one of the main reasons why it is so difficult to identify the precise role of hosts involved in the epidemiological cycle(s) of certain tick-borne pathogens.
Anaplasma phagocytophilum is a tick-borne intragranulocytic alpha-proteobacterium, mainly transmitted by I. ricinus in Europe, by I. scapularis and I. pacificus in North America, and by I. persulcatus in Asia. It infects a large range of hosts, including wild and domestic ruminants, dogs, horses, and rodents. Animal infection has been detected in Europe, North and South America, Asia, and Africa[8–10].
A. phagocytophilum is the causative agent of granulocytic anaplasmosis in humans, horses, dogs and occasionally cats, and tick-borne fever (TBF) in domestic ruminants. TBF is characterized by anorexia, agalactia, and in some cases, fetal abortion. The epidemiology of A. phagocytophilum infection differs greatly between the US and Europe. In the US, human granulocytic anaplasmosis (HGA) is an increasing public health problem (the CDC reported 1800 human cases in 2010, with a 0.7% fatality rate), whereas no US cases of TBF have been described to date. Conversely, HGA appears to be more rare in Europe (even though the number of reported cases has increased during recent years, probably linked in part to improved surveillance[12, 13]), whereas many TBF cases have been described in both cattle and sheep, causing significant economic losses[14–16]. These discrepancies between continents within one species, suggest A. phagocytophilum variability between and within host species and between continents. Experimental studies support the hypothesis that different epidemiological contexts are associated with considerable strain variation[17, 18]: an American strain infectious for horses was not infectious for ruminants, whereas a European variant pathogenic for cattle did not cause any clinical disease in horses. In the US, the Ap-ha variant, which is pathogenic for humans, can also infect both ruminants and mice under experimental conditions, whereas the Ap-Variant 1, which is not infectious for humans, can infect goats and deer, but not mice[21–24].
As A. phagocytophilum is not transmitted transovarially in ticks, it is thought to be maintained in vertebrate reservoir hosts. However, the reservoir hosts for European cattle and human strains have not yet been identified. Roe deer have been suspected as reservoir hosts for sheep strains in Norway. However to date, their role is still unclear. Isolates from either roe deer, or sheep and cattle, belonged to different clusters based on ankA gene phylogeny, whereas examining the groEL locus showed that isolates from domestic goats (Capra hircus), belonged to the same cluster as those from roe deer.
To understand the circulation of A. phagocytophilum, and to identify its reservoir hosts, the genetic diversity of this pathogen needs to be explored, requiring discriminant genetic markers. Currently, there are only a few molecular typing techniques available for A. phagocytophilum, which are currently unable to discriminate between all strains. These techniques include pulsed-field gel electrophoresis (PFGE), and single locus typing[26, 27, 29]. Two potentially more discriminant multi-locus approaches have been developed. The first approach is based on sequencing four loci including three genes and a 16S-RNA locus segment, and the second, called MLST (Multi-Locus Sequence Typing), relies on the sequencing of seven loci[30–33]. MLST is well adapted for phylogenetic analysis, however, this technique does not have enough discriminatory power for traceability studies, which are better performed by Multiple Locus VNTR (Variable-Number Tandem Repeat) Analysis (MLVA)[34–37]. MLVA determines the number of tandemly repeat sequences at different polymorphic VNTR loci within a bacterial genome. In a growing number of prokaryotes, which includes those displaying low genetic heterogeneity, VNTR typing has proved to be discriminatory, simple and transferable, with excellent marker stability. This type of epidemiological tool has now been used for many pathogenic bacteria, in particular for those bacteria with clinical impact on animals and/or with zoonotic impact[34, 39–43]. In a previous study, Bown et al. developed an MLVA technique for A. phagocytophilum, based on intergenic microsatellite VNTRs. Paradoxically, this technique was too discriminatory for use in epidemiology, when analyzing epidemiological links between isolates from different host species. Moreover, as the basic VNTR units were too short to discriminate between the different alleles by gel electrophoresis, any results had to be systematically determined or confirmed by sequencing or capillary electrophoresis.
In the present study we have developed an MLVA technique that would be more appropriate for A. phagocytophilum typing. Our objective was to explore the genetic diversity of A. phagocytophilum obtained from domestic and wild ruminants, in order to better understand the circulation of this pathogen between these two populations. This information will aid in subsequent efforts to control cattle infection. To address our objective, five new intragenic VNTR minisatellites were selected. Using these new MLVA markers, 97 samples from different domestic and wild animal species, 26 tick samples, the HZ genome, and the Webster human strain were analyzed. Resultant sample diversity was then analyzed and compared to current A. phagocytophilum diversity data.
The domestic animals used in this study met the definition of "farm animals" or "pets", which are not currently covered by French regulations (Décret 2013–118 dated the 1st February 2013, from the French Ministry of Agriculture). The owners of the animals provided permission for studies using samples obtained from their animals.The wild ruminant (roe deer specimens from French departments 12 and 21, Figure 1) and red deer specimens (French departments 1, 8 and 19, Figure 1) were killed by gunshot by licensed hunters, during the legal hunting season, which fell under official annual hunting quotas delineated by the county prefect. Specific accreditation to collect such samples was granted by the National Agency for Hunting and Wildlife (ONCFS, accreditation number 2009–2014). Specimens were obtained from animals which had been killed for hunting reasons, and were not specifically killed for this study. Thus, ethics committee approval was not required for this part of the study.
The wild roe deer from French department 8 were live-captured during the winter of 2013 using large-scale drives and long-nets. Blood samples were obtained from the jugular vein, and all roe deer were released after sample collection. The UR 035 team was granted administrative authorization by the Direction Départementale des Territoires from department 8 to catch and inventory wild animals. Blood sampling was performed by team veterinarians who were authorized to perform animal experimentation.
Samples and DNA extraction
DNA samples were obtained from animal blood or spleens, and from ticks. A. phagocytophilum was detected in these samples by PCR as previously described. Hereafter, the term "samples" applies to field samples that gave positive PCR results for A. phagocytophilum. The HZ strain genome and 124 A. phagocytophilum-positive samples were analyzed: 67 cattle (Bos taurus), 7 sheep (Ovis aries), 15 roe deer (Capreolus capreolus), 4 red deer (Cervus elaphus), 1 reindeer (Rangifer tarandus), 2 horses (Equus caballus), 1 dog (Canis lupus familiaris), 26 ticks (19 Rhipicephalus bursa, 5 Rhipicephalus sanguineus, 1 Rhipicephalus turanicus, and 1 Ixodes scapularis), and the North American human Webster strain. Sample characteristics are summarized in Table 1 (for more details see Additional file1). DNA was subsequently extracted from all samples.
For DNA extraction from blood or ticks, the NucleoSpin® Blood QuickPure kit (Macherey-Nagel, Bethlehem, USA), or the NucleoSpin® Tissue kit (Macherey-Nagel) were used respectively, according to manufacturer’s instructions. DNA extracts were then stored at -20°C prior to testing.
Computer analysis of repetitive DNA sequences for use as VNTR candidates
The genomic DNA sequence of the human A. phagocytophilum HZ strain (Reference Sequence NC_007797.1) was screened for repetitive DNA sequences using the tandem repeats database developed by Le Flèche et al..
The following criteria were applied to potential VNTR candidates: i/ Total length between 9 and 5000 bp; ii/ unit length between 9 and 500 bp; iii/ basic unit (BU) copy number between 1 and 2000; iv/ 70 to 100% BU identity.
BLASTN analysis of the repeat sequences excluded VNTRs present in other available genomes, especially those present in other Anaplasmataceae.
Primer design for VNTR candidates
The tandem repeats database described by Le Flèche et al. also provided 500 bp of flanking sequences both upstream and downstream of each VNTR candidate locus. BLASTN analysis of these flanking sequences was used to design forward and reverse primers. Primers were designed ensuring that no annealing would occur with; i) other A. phagocytophilum genome regions, ii) other available pathogen genomes and iii) available host genomes (i.e. human, both wild and domestic ruminants, canine, equine, rodent and tick), iv) the paired primer, or itself. Each VNTR candidate was designated as an APV (A naplasma p hagocytophilum VNTR), followed by a letter. The location of each APV in the A. phagocytophilum genome is described as the ‘locus’. An ‘allele’ corresponds to a given number of repeat units for a particular APV or locus. APV characteristics are summarized in Table 2.
Amplification, analysis and selection of VNTR candidates
Protocol optimization was carried out on five different samples, and the ability of each VNTR to explore diversity was then tested on all samples. VNTR amplification was conducted in a volume of 25 μl, containing 5 μl purified DNA, 4 μl 5X high fidelity amplification buffer, 200 μM each dNTP, 0.5 μM each primer (Eurofins MWG Operon, Ebersberg, Germany) and 0.4 units of Phusion DNA polymerase (Fisher Scientific, Waltham, USA). An initial denaturation step at 98°C for 30 s was followed by the subsequent thermocycling protocol: DNA was denatured for 10 s at 98°C, and primers were annealed for 30 s at the optimal temperature (56°C for APVA, APVB, APVC and APVD, and 58°C for APV-E), and extended at 72°C for 1 min. After 35 cycles, there was a final extension step at 72°C for 10 min.
PCR products were separated by electrophoresis on 1% SeaKem LE agarose gels (Ozyme, Saint Quentin en Yvelines, France) in TBE buffer (Lonza, Basel, Switzerland), and stained with ethidium bromide for imaging.
PCR products (20 μl) were sequenced by Sanger sequencing (Eurofins MWG Operon) to characterize the sequence of each new allele. Results were analyzed using Bioedit software version 7.2.5 (Ibis Biosciences, Carlsbad, USA).
GenBank accession numbers
For each VNTR, the alleles are listed in order.
[GenBank: KM216273 - KM216278]
APVA_4.5rep; APVA_3rep; APVA_2rep; APVA_1.5rep; APVA_1rep; APVA_0.5rep
[GenBank: KM216279 - KM216287]
APVB_9rep; APVB_8rep; APVB_7rep; APVB_6rep; APVB_5rep; APVB_4rep; APVB_3.5rep; APVB_2.5rep; APVB_2rep
[GenBank: M216288 - KM216296]
APVC_7.5rep; APVC_5rep; APVC_4rep; APVC_3.5; APVC_3rep; APVC_2rep; APVC_1.5rep; APVC_0.7rep; APVC_0.5rep
[GenBank: KM216297 - KM216300]
APVD_3rep; APVD_4rep; APVD_5rep; APVD_6rep
[GenBank: KM216301 - KM216311]
APVE_45rep; APVE_32rep; APVE_26rep; APVE_23rep; APVE_21rep; APVE_18rep; APVE_15rep; APVE_11rep; APVE_7rep; APVE_4rep; APVE_0.5rep
The Hunter and Gaston discrimination index (DI) was used to evaluate the discriminatory power of the selected APVs,, as recommended by the European Society of Clinical Microbiology and Infectious Diseases Study Group on Epidemiological Markers. This index measures the probability that two samples or strains, randomly chosen, will have different types. It is defined by:
N: number of samples or strains; S: total number of profiles or alleles and nj: number of samples or strains with the profile or allele j. Polymorphism rate is considered high when the index is higher than 95%.
Diversity of MLVA profiles obtained for cattle from the same herd vs cattle from different herds
The Fisher’s exact test was used to compare the proportion of A. phagocytophilum from cattle sharing the same profile among two groups: animals originating from distinct herds (one bovine sample per herd), and animals having common origins (at least two bovine per herd). Statistical analysis was conducted at 95% confidence level. A p value less than 0.05 was considered statistically significant.
Similarity of MLVA profiles obtained within the same species
The similarities of MLVA profiles obtained from the same animal species were analyzed using a bootstrap method as previously described. Only animal species for which two or more MLVA profiles had been observed were taken into account. The null hypothesis was that two pairs of MLVA profiles obtained from the same species, should be as similar as two pairs of MLVA profiles obtained from different species. The distance between two MLVA profiles was the number of VNTRs with differential numbers of repeats (distances thus varied between 0 and 5). The analysis was focused on the average of within-species mean distances between MLVA profiles. The observed value of this statistic of interest was first computed. A resampling procedure was then used to simulate bootstrap samples from the data, under the null hypothesis. The resampling procedure was based on an adjacency matrix that links species (matrix rows) and MLVA profiles (matrix columns). Each bootstrap sample was obtained by generating random permutations of the matrix column headers (the MLVA profiles). This procedure guarantees that the marginal sums of the adjacency matrix are constant (i.e. number of profiles per species, number of species per profile), whereas the profiles are randomly redistributed. Ten thousand bootstrap samples were thus generated and, for each of these, the statistic of interest was calculated. The corresponding distribution was finally examined to determine the empirical p-value of the null hypothesis test: this p-value was the proportion of the samples (simulated under the null hypothesis) for which the statistic of interest was below the actual value (computed from the data).
MLVA clustering was performed using the BioNumerics software package version 6.6 (Applied-Maths, Sint-Martens-Latem, Belgium), with the categorical distance coefficient and UPGMA (Un-weighted Pair Group Method with Arithmetic mean) clustering method. Data were analyzed as a character dataset. Based on this categorical distance coefficient analysis, the Minimum Spanning Tree (MST) graphing algorithm was used to represent the relationships between strains. The priority rule for constructing MSTs was set so that the type which had the highest number of single-locus variants would be linked first. A cut-off value of maximum differences of 1 VNTR was applied to define clonal complexes under the MST method.
Using the tandem repeats database developed by Le Flèche et al., we pre-selected 721 VNTR candidates according to the criteria described above. APV candidates corresponding to the selection criteria were further tested by evaluating their polymorphism. Following amplification optimization (amplification cycle and temperature gradient modifications), only those VNTRs that were amplified and appeared to be polymorphic when tested on five of our samples, were then selected and used in this study. These VNTRs corresponded to those listed as APV-A, APV-B, APV-C, APV-D and APV-E. Except for APV-E (15 bp), all VNTR BUs are longer than 100 bp. The VNTR names, BU lengths, primer sequences, gene and genome localization, and allele size ranges are shown in Table 2.
Specificity of the MLVA technique
In silico analysis of all selected primers in BLASTN confirmed that they were unable to hybridize with either the genomes of other known animal pathogens, or those of microorganisms found in ticks. This feature was experimentally confirmed, as no PCR amplification was observed for DNA extracts from these microorganisms.
The length of each allele was confirmed by sequencing. The different allele sizes are detailed in Table 2.
The discriminatory power of each APV locus was estimated from the genetic DI values based on the numbers of alleles and their frequencies in different samples. The high individual DI values for the 124 samples, as well as the HZ genome included in our study, reflect their excellent potential as markers of genetic diversity: 0.80 for APV-A, 0.74 for APV-B, 0.74 for APV-C, 0.50 for APV-D, and 0.87 for APV-E (and without tick samples: 0.76 for APV-A, 0.69 for APV-B, 0.76 for APV-C, 0.60 for APV-D, and 0.88 for APV-E). Moreover, the global DI value, when combining the five APVs, was 0.96 (and 0.99 without tick samples). In total, we observed 84 different profiles for 125 samples.
Sixty-seven cattle samples tested positive for A. phagocytophilum. Of these, twenty-one were obtained from seven farms located in various French departments, indicated in Figure 1: three samples from a Côtes-d’Armor farm (department 4 in Figure 1), two samples from a Finistère farm (department 6), three samples from a Haute-Vienne farm (department 23), four samples from an Ille-et-Vilaine farm (department 10), three samples from a Mayenne farm (department 15), and six samples from two Morbihan farms (three samples per farm, department 16). Different samples from the same farms harbored similar or identical profiles. For this set of 21 samples, 38% (n = 8) shared their profile with at least one other sample. For the 46 remaining cattle samples from different herds, our technique identified 44 distinct profiles. Within this group, 7% (n = 3) shared their profile with at least one other sample, the proportion of which was significantly higher for the first group (≥2 cattle per farm) than for the second group (1 cow per farm) (p = 0.003). This difference was also significant for cattle samples from Bretagne (French departments 4, 6, 10, 14, 15 and 16), where the majority of farms with multiple samples were located (exact Fisher’s test, p = 0.01347).
Two different but very similar profiles were observed for the two horse samples from the same area, as only the APV-E allele varied between those samples. Seven different profiles were obtained for seven sheep samples from one Pyrénées-Atlantiques herd. Three different but close profiles were observed for 10 roe deer samples from French department 8. These three MLVA profiles differed only by their APV-C alleles. Four different profiles were obtained for five roe deer samples from Franche-Comté (French departments 12 and 21), which differed according to their APV-A, APV-C and APV-E alleles. And finally, four different profiles were obtained from four red deer samples (French departments (1, 8 and 19). According to Figure 2, these profiles are very close to those from both cattle and sheep.
All 25 tick samples from the three different French departments of Camargue (15 from the Bouches-du-Rhône, five from Gard and five from the Hérault department; Figure 1, Table 1 and Additional file2), a 250-km2 area in southern France, all harbored the same profile. Conversely, the tick sample aph1 (North American strain) provided a profile that was different from that observed in Camargue ticks.
A total of 125 samples (84 VNTR profiles) were utilized for this analysis. Bootstrap analysis showed that MLVA profiles obtained from the same host species were significantly closer than the null hypothesis (according to which, the average distance between MLVA profiles is the same between species and within species) (p < 0.0001). MST analysis clearly differentiated A. phagocytophilum MLVA profiles from cattle, sheep, and red deer, compared to A. phagocytophilum from roe deer, positioning them in two different clusters (Figure 2). Moreover, A. phagocytophilum from domestic ruminants shared alleles at the APV-B locus with red deer, but not with A. phagocytophilum from roe deer.
MLVA has become a popular typing approach for studying the epidemiology of pathogenic agents, due to several advantages: high discriminatory power, robustness, repeatability, inter-laboratory portability and speed[53, 54]. In this study, we developed an MLVA technique for A. phagocytophilum, an intracellular tick-borne bacterium which can be isolated from both humans and various animal species. To the best of our knowledge, our work is the first study to be performed in France which typed A. phagocytophium from domestic and wild ruminants. Our MLVA technique has been shown to be effective for the typing of human, cattle, sheep, roe deer, red deer, reindeer, dog, horse, and tick isolates. This technique is based on five VNTRs with a global Hunter and Gaston discrimination index of 0.96 (84 profiles/125 samples), or 0.99 without tick samples, emphasizing its high discriminatory power.
In a previous study, Bown et al. developed an MLVA technique based on intergenic minisatellite VNTRs, but paradoxically, their method was too discriminatory. The 20 strains tested in their study shared only a few alleles at each of the four loci (VNTR1 and 4: 17 different alleles/20 isolates; VNTR5 and 8: 14 different alleles/20 isolates, where each given allele was either shared by one isolate – a majority of cases – or up to 5 isolates). For this reason, their technique does not seem well adapted to studying either A. phagocytophilum transmission between species, or the circulation of isolates within cattle herds, or between wild animals – such as rodents – and domestic ruminants. This was an additional reason why we based our MLVA technique on intragenic minisatellite VNTRs, which are probably under higher selective pressure, and could be less variable, than those selected by Bown et al. as already shown for other bacteria. Our results confirmed this hypothesis, as we obtained, among the 125 samples tested, six different alleles for APV-A, 10 different alleles for APV-B and APV-C, five different alleles for APV-D, and 14 different alleles for APV-E. Moreover, our technique appeared to efficiently discriminate related strains isolated from the same area and/or within the same animal species, and even regardless of the species. The two horse samples, both collected in Lyon (France), shared four alleles. APV-E, which is the most variable of our markers, was the only VNTR that differentiated between these two samples, differing by 30 nucleotides. In addition, whereas cattle samples isolated from different herds gave 44 profiles from 45 samples, samples obtained from cattle belonging to the same herds gave identical or very similar MLVA profiles. These results were confirmed even for cattle originating from the same region (Bretagne), which therefore excluded any links between geographical location and profile distribution. Conversely, we obtained highly diverse MLVA profiles within the sheep herd that was tested. If these divergent results were confirmed in more herds under similar conditions, they could perhaps be explained by differences in farming methods between the two animal species. For example, whereas cattle are confined to small pastures, sheep have extensive grazing areas, and are also subject to seasonal transhumance across a variety of territories, thus increasing their likelihood of exposure to different ticks, and consequently, diverse A. phagocytophilum strains. In addition, cattle are not persistently infected, in contrast to sheep. Thus genetic shifts and/or inter-strain recombinations within the same sheep flock cannot be excluded, whereas the isolates obtained for cattle belonging to the same farm are most likely possibly derived from the same parental strain. This last result suggests that in vivo variations could occur under the selective pressure of the host immune system. However, no variation was detected after eight passages of the Webster strain in HL-60 cells, corresponding to approximately one year of culture (results not shown), suggesting that the observed lack of VNTR variations could be due to the absence of in vitro selective factors. Unfortunately, as artificially infecting animals with A. phagocytophilum requires specific animal facilities, we could not confirm whether in vivo variations are likely to occur in the field.
Taken together, our results strongly suggest that this particular MLVA technique has a good concordance with epidemiological contexts. The in-depth resolution obtained will be of great value for studying A. phagocytophilum circulation within cattle, or between wildlife and domestic ruminants. Epidemiological studies with suitable samples could validate this assertion, as well as the applicability of our MLVA technique for such studies.
We also noticed that the APV-B VNTR was systematically different in roe deer from other ruminant species samples. As our sampling is limited, further analysis of additional samples is needed in order to confirm or not if this VNTR could be considered as a marker of roe deer strains.
Our results suggest that red deer are reservoir hosts for domestic ruminant strains, in contrast to roe deer. Nevertheless, this observation has to be taken with caution as we could only test four red deer samples that were all collected from the same region, which does not permit the examination of geographical effects. However, red deer are becoming increasingly suspected as having an important role in the epidemiological cycle of A. phagocytophilum, as reservoir hosts. Different arguments support this hypothesis. Firstly, sheep experimentally infected with A. phagocytophilum red deer isolates developed TBF symptoms. Secondly, high prevalence of A. phagocytophilum infection has been detected in red deer in many European countries.[27, 56–61]. In France, 80% of red deer (28/35) recently tested by Zehnter et al. from one French region (Correze), were reported to be seropositive for A. phagocytophilum. According to official hunting bag data (ONCFS: the French national agency for wildlife), the average red deer population has increased three-fold in France over a twenty year period, and is continuing to increase[63, 64]. Unfortunately, no data are available which could either confirm or deny whether a parallel increase in French BGA incidence has occurred during the same period.
Thus, to date, the circulation of isolates between wild and domestic ruminants still remains unclear. Within this context, our MLVA technique could be a very useful tool with which to address this problem. For a reliable approach, more wild and domestic animal samples originating from the same regions will be required.
In addition, our results suggest the existence of an alternative epidemiological cycle in the French Camargue, involving Rhipicephalus ticks as potential vectors.
In Europe, A. phagocytophilum is mainly transmitted by I. ricinus, but A. phagocytophilum DNA has also been detected in other tick genera, such as Dermacentor and Rhipicephalus[65, 66]. Nevertheless, the ability of those ticks to transmit A. phagocytophilum has not yet been proven. In contrast to Rhipicephalus ticks, I. ricinus is considered to be rare, or even absent in the Camargue. It means that in this region, A. phagocytophilum is, at least predominantly, transmitted by another vector, which could be Rhipicephalus ticks.
Chastagner et al. have found only one A. phagocytophilum MLST genotype in Rhipicephalus ticks in the Camargue. Using our MLVA technique on the same samples, we have shown that all presented the same profile. Moreover, the fact that one North American tick sample provided a profile that was completely different from the profile obtained from Camargue tick samples, confirms that our MLVA technique is effective in revealing diversity at the level of tick isolates. As the ticks were collected at different sites, over an area of 250 km2, our results suggest the circulation of one ecotype of A. phagocytophilum in the French Camargue, which is adapted to Rhipicephalus ticks, representing a very particular ecosystem. In order to test this hypothesis, the competence of Rhipicephalus ticks for A. phagocytophilum transmission should be explored. Using our MLVA technique, we would be able to determine which variant(s) are detected, and if they could be multiplied and transmitted by Rhipicephalus ticks. Typing a large number of Camargue samples collected from different animal species would be an excellent complementary method to further characterize the epidemiological cycle of A. phagocytophilum in this very distinct French region, and in particular, the role of Rhipicephalus ticks.
In conclusion, our results potentially reveal the existence of at least two different epidemiological transmission cycles of A. phagocytophilum, as has been previously described in the United Kingdom. The first cycle may involve red deer as reservoir hosts, and possibly domestic ruminants, as either accidental or longer-term hosts, whereas the second might involve roe deer. In addition, our study reveals the existence of a potentially alternative epidemiological cycle in the Camargue region, which could involve Rhipicephalus ticks as vectors. These hypotheses should be further explored, and the MLVA technique described here would be suitable for such studies. Our MLVA technique could also be used in order to compare the molecular profiles of A. phagocytophilum from different areas and/or countries, as has already been performed for B. henselae in our laboratory[52, 68, 69]. In the case of A. phagocytophilum, this would be an important contribution to the in-depth characterization of European epidemiological cycles, and aid in the comparison with the rather different cycles identified in the USA.
Anaplasma phagocytophilum VNTR
Bovine granulocytic anaplasmosis
Basic local alignment search tool
Centers for disease control and prevention
Human granulocytic anaplasmosis
Multi-locus sequence typing
Multiple-locus variable-number tandem repeat analysis
Minimum spanning tree
French national agency for wildlife
Pulsed-field gel electrophoresis
Unweighted pair group method with arithmetic mean
Variable number tandem repeat.
Pedersen AB, Altizer S, Poss M, Cunningham AA, Nunn CL: Patterns of host specificity and transmission among parasites of wild primates. Int J Parasitol. 2005, 35: 647-657. 10.1016/j.ijpara.2005.01.005.
Cleaveland S, Laurenson MK, Taylor LH: Diseases of humans and their domestic mammals: pathogen characteristics, host range and the risk of emergence. Philos Trans R Soc Lond B Biol Sci. 2001, 356: 991-999. 10.1098/rstb.2001.0889.
Woolhouse ME, Taylor LH, Haydon DT: Population biology of multihost pathogens. Science. 2001, 292: 1109-1112. 10.1126/science.1059026.
Bordes F, Morand S: The impact of multiple infections on wild animal hosts: a review. Infect Ecol Epidemiol. 2011, 1: 1-10.
Pedersen AB, Fenton A: Emphasizing the ecology in parasite community ecology. Trends Ecol Evol. 2007, 22: 133-139. 10.1016/j.tree.2006.11.005.
Salman MD, Tarrés-Call J, Estrada-Peña A: Ticks and Tick-Borne Diseases: Geographical Distribution and Control Strategies in the Euro-Asia Region. 2013, Colorado State University Press: Fort Collins
Dumler JS, Barbet AF, Bekker CP, Dasch GA, Palmer GH, Ray SC, Rikihisa Y, Rurangirwa FR: Reorganization of genera in the families Rickettsiaceae and Anaplasmataceae in the order Rickettsiales: unification of some species of Ehrlichia with Anaplasma, Cowdria with Ehrlichia and Ehrlichia with Neorickettsia, descriptions of six new species combinations and designation of Ehrlichia equi and "HGE agent" as subjective synonyms of Ehrlichia phagocytophila. Int J Syst Evol Microbiol. 2001, 51: 2145-2165. 10.1099/00207713-51-6-2145.
Stuen S, Granquist EG, Silaghi C: Anaplasma phagocytophilum–a widespread multi-host pathogen with highly adaptive strategies. Front Cell Infect Microbiol. 2013, 3: 1-33.
M’ghirbi Y, Yaïch H, Ghorbel A, Bouattour A: Anaplasma phagocytophilum in horses and ticks in Tunisia. Parasit Vectors. 2012, 5: 180-10.1186/1756-3305-5-180.
Djiba ML, Mediannikov O, Mbengue M, Thiongane Y, Molez J-F, Seck MT, Fenollar F, Raoult D, Ndiaye M: Survey of Anaplasmataceae bacteria in sheep from Senegal. Trop Anim Health Prod. 2013, 45: 1557-1561. 10.1007/s11250-013-0399-y.
Centers for Disease Control and Prevention, Statistics and Epidemiology of Anaplasmosis. [http://www.cdc.gov/anaplasmosis/stats/],
Edouard S, Koebel C, Goehringer F, Socolovschi C, Jaulhac B, Raoult D, Brouqui P: Emergence of human granulocytic anaplasmosis in France. Ticks Tick-Borne Dis. 2012, 3: 403-405. 10.1016/j.ttbdis.2012.10.002.
Cochez C, Ducoffre G, Vandenvelde C, Luyasu V, Heyman P: Human anaplasmosis in Belgium: a 10-year seroepidemiological study. Ticks Tick-Borne Dis. 2011, 2: 156-159. 10.1016/j.ttbdis.2011.06.004.
Woldehiwet Z: Anaplasma phagocytophilum in ruminants in Europe. Ann N Y Acad Sci. 2006, 1078: 446-460. 10.1196/annals.1374.084.
Pusterla N, Braun U: Clinical findings in cows after experimental infection with Ehrlichia phagocytophila. Zentralblatt Für Veterinärmedizin Reihe A. 1997, 44: 385-390.
Brodie TA, Holmes PH, Urquhart GM: Some aspects of tick-borne diseases of British sheep. Vet Rec. 1986, 118: 415-418. 10.1136/vr.118.15.415.
Foley J, Nieto NC, Madigan J, Sykes J: Possible differential host tropism in Anaplasma phagocytophilum strains in the Western United States. Ann N Y Acad Sci. 2008, 1149: 94-97. 10.1196/annals.1428.066.
Morissette E, Massung RF, Foley JE, Alleman AR, Foley P, Barbet AF: Diversity of Anaplasma phagocytophilum Strains, USA. Emerg Infect Dis. 2009, 15: 928-931. 10.3201/eid1506.081610.
Stannard AA, Gribble DH, Smith RS: Equine ehrlichiosis: a disease with similarities to tick-borne fever and bovine petechial fever. Vet Rec. 1969, 84: 149-150. 10.1136/vr.84.6.149.
Pusterla N, Lutz H, Braun U: Experimental infection of four horses with Ehrlichia phagocytophila. Vet Rec. 1998, 143: 303-305. 10.1136/vr.143.11.303.
Massung RF, Priestley RA, Miller NJ, Mather TN, Levin ML: Inability of a variant strain of Anaplasma phagocytophilum to infect mice. J Infect Dis. 2003, 188: 1757-1763. 10.1086/379725.
Massung RF, Mather TN, Levin ML: Reservoir competency of goats for the Ap-variant 1 strain of Anaplasma phagocytophilum. Infect Immun. 2006, 74: 1373-1375. 10.1128/IAI.74.2.1373-1375.2006.
Massung RF, Courtney JW, Hiratzka SL, Pitzer VE, Smith G, Dryden RL: Anaplasma phagocytophilum in white-tailed deer. Emerg Infect Dis. 2005, 11: 1604-1606. 10.3201/eid1110.041329.
Tate CM, Mead DG, Luttrell MP, Howerth EW, Dugan VG, Munderloh UG, Davidson WR: Experimental infection of white-tailed deer with Anaplasma phagocytophilum, etiologic agent of human granulocytic anaplasmosis. J Clin Microbiol. 2005, 43: 3595-3601. 10.1128/JCM.43.8.3595-3601.2005.
Stuen S, Scharf W, Schauer S, Freyburger F, Bergström K, von Loewenich FD: Experimental infection in lambs with a red deer (Cervus elaphus) isolate of Anaplasma phagocytophilum. J Wildl Dis. 2010, 46: 803-809. 10.7589/0090-3558-46.3.803.
Scharf W, Schauer S, Freyburger F, Petrovec M, Schaarschmidt-Kiener D, Liebisch G, Runge M, Ganter M, Kehl A, Dumler JS, Garcia-Perez AL, Jensen J, Fingerle V, Meli ML, Ensser A, Stuen S, von Loewenich FD: Distinct host species correlate with Anaplasma phagocytophilum ankA gene clusters▿. J Clin Microbiol. 2011, 49: 790-796. 10.1128/JCM.02051-10.
Silaghi C, Hamel D, Thiel C, Pfister K, Passos LMF, Rehbein S: Genetic variants of Anaplasma phagocytophilum in wild caprine and cervid ungulates from the Alps in Tyrol, Austria. Vector Borne Zoonotic Dis Larchmt N. 2011, 11: 355-362. 10.1089/vbz.2010.0051.
Dumler JS, Asanovich KM, Bakken JS: Analysis of genetic identity of North American Anaplasma phagocytophilum strains by pulsed-field gel electrophoresis. J Clin Microbiol. 2003, 41: 3392-3394. 10.1128/JCM.41.7.3392-3394.2003.
Bown KJ, Lambin X, Ogden NH, Begon M, Telford G, Woldehiwet Z, Birtles RJ: Delineating Anaplasma phagocytophilum ecotypes in coexisting, discrete enzootic cycles. Emerg Infect Dis. 2009, 15: 1948-1954. 10.3201/eid1512.090178.
Shukla SK, Aswani V, Stockwell PJ, Reed KD: Contribution of polymorphisms in ankA, gltA, and groESL in defining genetic variants of Anaplasma phagocytophilum. J Clin Microbiol. 2007, 45: 2312-2315. 10.1128/JCM.00704-07.
Zeman P, Jahn P: An entropy-optimized multilocus approach for characterizing the strains of Anaplasma phagocytophilum infecting horses in the Czech Republic. J Med Microbiol. 2009, 58 (Pt 4): 423-429.
Huhn C, Winter C, Wolfsperger T, Wüppenhorst N, Strašek Smrdel K, Skuballa J, Pfäffle M, Petney T, Silaghi C, Dyachenko V, Pantchev N, Straubinger RK, Schaarschmidt-Kiener D, Ganter M, Aardema ML, von Loewenich FD: Analysis of the population structure of Anaplasma phagocytophilum using multilocus sequence typing. PLoS ONE. 2014, 9: e93725-10.1371/journal.pone.0093725.
Chastagner A, Bailly X, Leblond A, Pradier S, Vourc’h G: Single genotype of Anaplasma phagocytophilum identified from ticks, Camargue, France. Emerg Infect Dis. 2013, 19: 825-826.
Bouchouicha R, Boulouis H-J, Berrich M, Monteil M, Chomel B, Haddad N: Comparison of the performances of MLVA vs. the main other typing techniques for Bartonella henselae. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2009, 15 (Suppl 2): 104-105.
Elberse KEM, Nunes S, Sa-Leao R, van der Heide HGJ, Schouls LM: Multiple-locus variable number tandem repeat analysis for Streptococcus pneumoniae: comparison with PFGE and MLST. PLoS ONE. 2011, 6: 1-8.
Maatallah M, Bakhrouf A, Habeeb MA, Turlej-Rogacka A, Iversen A, Pourcel C, Sioud O, Giske CG: Four genotyping schemes for phylogenetic analysis of Pseudomonas aeruginosa: comparison of their congruence with multi-locus sequence typing. PLoS ONE. 2013, 8: 1-10.
Riegler M, Iturbe-Ormaetxe I, Woolfit M, Miller WJ, O’Neill SL: Tandem repeat markers as novel diagnostic tools for high resolution fingerprinting of Wolbachia. BMC Microbiol. 2012, 12 (Suppl 1): S12-10.1186/1471-2180-12-S1-S12.
Lindstedt B-A: Multiple-locus variable number tandem repeats analysis for genetic fingerprinting of pathogenic bacteria. Electrophoresis. 2005, 26: 2567-2582. 10.1002/elps.200500096.
Beyer W, Turnbull PCB: Co-infection of an animal with more than one genotype can occur in anthrax. Lett Appl Microbiol. 2013, 57: 380-384.
Hill V, Zozio T, Sadikalay S, Viegas S, Streit E, Kallenius G, Rastogi N: MLVA based classification of Mycobacterium tuberculosis complex lineages for a robust phylogeographic snapshot of its worldwide molecular diversity. PLoS ONE. 2012, 7: 1-16.
Li ZJ, Cui BY, Chen H, Chen JD, Zhao HY, Piao DR, Jiang H, Zhang L, Tang X, Ke CW, Yao Z, Tian GZ: Molecular typing of Brucella suis collected from 1960s to 2010s in China by MLVA and PFGE. Biomed Environ Sci BES. 2013, 26: 504-508.
Pilet H, Vachiéry N, Berrich M, Bouchouicha R, Durand B, Pruneau L, Pinarello V, Saldana A, Carasco-Lacombe C, Lefrançois T, Meyer DF, Martinez D, Boulouis H-J, Haddad N: A new typing technique for the Rickettsiales Ehrlichia ruminantium: Multiple-locus variable number tandem repeat analysis. J Microbiol Methods. 2012, 88: 205-211. 10.1016/j.mimet.2011.11.011.
Saleh-Lakha S, Allen VG, Li J, Pagotto F, Odumeru J, Taboada E, Lombos M, Tabing KC, Blais B, Ogunremi D, Downing G, Lee S, Gao A, Nadon C, Chen S: Subtyping of a large collection of historical Listeria monocytogenes strains from Ontario, Canada, by an improved multilocus variable-number tandem-repeat analysis (MLVA). Appl Environ Microbiol. 2013, 79: 6472-6480. 10.1128/AEM.00759-13.
Bown KJ, Lambin X, Ogden NH, Petrovec M, Shaw SE, Woldehiwet Z, Birtles RJ: High-resolution genetic fingerprinting of European strains of Anaplasma phagocytophilum by use of multilocus variable-number tandem-repeat analysis. J Clin Microbiol. 2007, 45: 1771-1776. 10.1128/JCM.00365-07.
Drazenovich N, Foley J, Brown RN: Use of real-time quantitative PCR targeting the msp2 protein gene to identify cryptic Anaplasma phagocytophilum infections in wildlife and domestic animals. Vector Borne Zoonotic Dis Larchmt N. 2006, 6: 83-90. 10.1089/vbz.2006.6.83.
Chen SM, Dumler JS, Bakken JS, Walker DH: Identification of a granulocytotropic Ehrlichia species as the etiologic agent of human disease. J Clin Microbiol. 1994, 32: 589-595.
Dunning Hotopp JC, Lin M, Madupu R, Crabtree J, Angiuoli SV, Eisen J, Seshadri R, Ren Q, Wu M, Utterback TR, Smith S, Lewis M, Khouri H, Zhang C, Niu H, Lin Q, Ohashi N, Zhi N, Nelson W, Brinkac LM, Dodson RJ, Rosovitz MJ, Sundaram J, Daugherty SC, Davidsen T, Durkin AS, Gwinn M, Haft DH, Selengut JD, Sullivan SA, Zafar N, Zhou L: Comparative genomics of emerging human ehrlichiosis agents. PLoS Genet. 2006, 2: 208-223. 10.1371/journal.pgen.0020208.
Le Fleche P, Hauck Y, Onteniente L, Prieur A, Denoeud F, Ramisse V, Sylvestre P, Benson G, Ramisse F, Vergnaud G: A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis. BMC Microbiol. 2001, 1: 2-10.1186/1471-2180-1-2.
Hunter PR, Gaston MA: Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988, 26: 2465-2466.
Struelens MJ: Consensus guidelines for appropriate use and evaluation of microbial epidemiologic typing systems. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 1996, 2: 2-11.
Felsenstein J: Mathematics vs. evolution: mathematical evolutionary theory. Science. 1989, 246: 941-942. 10.1126/science.246.4932.941.
Azzag N, Haddad N, Durand B, Petit E, Ammouche A, Chomel B, Boulouis H-J: Population structure of Bartonella henselae in Algerian urban stray cats. PLoS ONE. 2012, 7: 1-13.
Pourcel C, Minandri F, Hauck Y, D’Arezzo S, Imperi F, Vergnaud G, Visca P: Identification of Variable-Number Tandem-Repeat (VNTR) sequences in Acinetobacter baumannii and interlaboratory validation of an optimized multiple-locus VNTR analysis typing scheme. J Clin Microbiol. 2011, 49: 539-548. 10.1128/JCM.02003-10.
Vergnaud G, Pourcel C: Multiple locus variable number of tandem repeats analysis. Methods Mol Biol Clifton NJ. 2009, 551: 141-158. 10.1007/978-1-60327-999-4_12.
Kimura M, Sakamuri RM, Groathouse NA, Rivoire BL, Gingrich D, Krueger-Koplin S, Cho S-N, Brennan PJ, Vissa V: Rapid variable-number tandem-repeat genotyping for Mycobacterium leprae clinical specimens. J Clin Microbiol. 2009, 47: 1757-1766. 10.1128/JCM.02019-08.
Alberdi MP, Walker AR, Urquhart KA: Field evidence that roe deer (Capreolus capreolus) are a natural host for Ehrlichia phagocytophila. Epidemiol Infect. 2000, 124: 315-323. 10.1017/S0950268899003684.
Groen J, Koraka P, Nur YA, Avsic-Zupanc T, Goessens WHF, Ott A, Osterhaus ADME: Serologic evidence of ehrlichiosis among humans and wild animals in The Netherlands. Eur J Clin Microbiol Infect Dis Off Publ Eur Soc Clin Microbiol. 2002, 21: 46-49. 10.1007/s10096-001-0659-z.
Liz JS: Ehrlichiosis in Ixodes ricinus and wild mammals. Int J Med Microbiol IJMM. 2002, 291 (Suppl 33): 104-105.
De la Fuente J, Ruiz-Fons F, Naranjo V, Torina A, Rodríguez O, Gortázar C: Evidence of Anaplasma infections in European roe deer (Capreolus capreolus) from southern Spain. Res Vet Sci. 2008, 84: 382-386. 10.1016/j.rvsc.2007.05.018.
Zele D, Avberšek J, Gruntar I, Ocepek M, Vengušt G: Evidence of Anaplasma phagocytophilum in game animals from Slovenia. Acta Vet Hung. 2012, 60: 441-448. 10.1556/AVet.2012.038.
Stuen S, Pettersen KS, Granquist EG, Bergström K, Bown KJ, Birtles RJ: Anaplasma phagocytophilum variants in sympatric red deer (Cervus elaphus) and sheep in southern Norway. Ticks Tick-Borne Dis. 2013, 4: 197-201. 10.1016/j.ttbdis.2012.11.014.
Zehnter A: Etude épidémiologique de maladies bactériennes vectorisées sur la faune sauvage et les bovins de Corrèze. DVM thesis. 2013, Ecole Nationale Vétérinaire d’Alfort: Maisons-Alfort
ONCFS RO sauvages: Tableaux de chasse ongulés sauvages saison 2011-2012 et évolution départementale par décennie. Lettre d’information du réseau Ongulés sauvages. 2013, 17: 2-3.
ONCFS RO sauvages: Tableaux de chasse ongulés sauvages saison 2012–2013. Lettre d’information du réseau Ongulés sauvages. 2013, 301: 1-8.
Bonnet S, de la Fuente J, Nicollet P, Liu X, Madani N, Blanchard B, Maingourd C, Alongi A, Torina A, Fernándezde Mera IG, Vicente J, George J-C, Vayssier-Taussat M, Joncour G: Prevalence of tick-borne pathogens in adult Dermacentor spp. ticks from nine collection sites in France. Vector Borne Zoonotic Dis Larchmt N. 2013, 13: 226-236. 10.1089/vbz.2011.0933.
de la Fuente J, Naranjo V, Ruiz-Fons F, Vicente J, Estrada-Peña A, Almazán C, Kocan KM, Martín MP, Gortázar C: Prevalence of tick-borne pathogens in ixodid ticks (Acari: Ixodidae) collected from European wild boar (Sus scrofa) and Iberian red deer (Cervus elaphus hispanicus) in central Spain. Eur J Wildl Res. 2004, 50: 187-196. 10.1007/s10344-004-0060-1.
Pérez-Eid C: Les tiques, Identification, Biologie, Importance Médicale Et Vétérinaire. 2007, Cachan: TEC&DOC, Lavoissier Press
Bouchouicha R, Durand B, Monteil M, Chomel BB, Berrich M, Arvand M, Birtles RJ, Breitschwerdt EB, Koehler JE, Maggi R, Maruyama S, Kasten R, Petit E, Boulouis H-J, Haddad N: Molecular epidemiology of feline and human Bartonella henselae isolates. Emerg Infect Dis. 2009, 15: 813-816. 10.3201/eid1505.080995.
Podsiadly E, Haddad N, Berrich M, Bouchouicha R, Durand B, Monteil M, Boulouis H-J, Tylewska-Wierzbanowska S: Characterization of Polish feline B. henselae isolates by multiple-locus tandem repeat analysis and pulse-field gel electrophoresis. Ann Agric Environ Med AAEM. 2012, 19: 39-43.
The authors thank Prof. Stephen Dumler (The Johns Hopkins University School of Medicine, Baltimore, USA.); Prof. Philippe Brouqui and Prof. Didier Raoult (Faculty of medicine, URMIT Joint Unit, Marseille) for providing the Webster Strain; Dr Michael L. Levin (Center for Disease Control and Prevention, Atlanta, USA) for providing DNA from American tick sample aph-1; Charlette Chandosné and the French Hunters’ Federation of Franche-Comté; researchers, technicians and volunteers from the INRA-CEFS of Toulouse that provided help to trap the roe deer; Dr. Nathalie Vassallo and the Laboratory of Côtes d’Armor LDA22 for providing bovine blood samples, and Dr. Martine Monteil and Guillaume Girault for their technical support and helpful advice. This work was performed within the context of the Laboratory of Excellence (Labex) of Integrative Biology of Emerging Infectious Diseases (IBEID).
The authors declare that they have no competing interests.
TD contributed to experimental design, to MLVA data analysis, was responsible for laboratory work, and drafted the manuscript. AC, ACL and EP participated in laboratory work as well as critical revision of the manuscript. BD was responsible for statistical analysis, contributed to MLVA data analysis and to drafting the manuscript. SM was responsible for MST analysis and helped to draft the manuscript. XB, FC, LC, AL, GV and HV were involved in sample collection and in critical revision of the manuscript. HJB, NH and RM were responsible for the conception and contributed to experimental design, to MLVA data analysis and to drafting the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Origin (host species, herd, and geographical area), and year of sampling for the A. phagocytophilum field samples, strain and genome used in this study.(XLS 34 KB)
About this article
Cite this article
Dugat, T., Chastagner, A., Lagrée, AC. et al. A new multiple-locus variable-number tandem repeat analysis reveals different clusters for Anaplasma phagocytophilum circulating in domestic and wild ruminants. Parasites Vectors 7, 439 (2014). https://doi.org/10.1186/1756-3305-7-439