Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2003, p. 1609-1616, Vol. 41, No. 4
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.4.1609-1616.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Genotypes Evolved under Positive Selection Pressure but Are Not Markers for Geographic Isolates
Department of Veterinary Pathobiology,1 Veterinary Teaching Hospital, College of Veterinary Medicine,3 Department of Zoology and Collection of Vertebrates, Oklahoma State University, Stillwater, Oklahoma 747082
Received 15 July 2002/ Returned for modification 16 December 2002/ Accepted 23 January 2003
|
|
|---|
is not a marker for the characterization of A. marginale geographic isolates and suggest that the genetic heterogeneity observed among isolates of A. marginale within Oklahoma could be explained by cattle movement and the maintenance of different genotypes by independent transmission events. |
|
|---|
The organisms in the order Rickettsiales have been reclassified into two families, the Anaplasmataceae and Rickettsiaceae, on the basis of genetic analyses of the 16S rRNA, groesl, and surface protein genes (12). Recent research has focused on MSPs that are involved in interactions with vertebrate and invertebrate host cells (7, 8, 20, 21) and have been used to elucidate the phylogeographic patterns of A. marginale (5, 6, 9, 16, 19). These MSPs are involved in host-pathogen interactions and may evolve more rapidly than other nuclear genes because of selective pressures exerted by host immune systems. Of the six A. marginale MSPs that have been identified and characterized, only three (MSP1a, MSP4, and MSP5) are encoded by single genes. Because these MSPs do not appear to undergo antigenic variation in cattle or ticks (4), they were posited to be more stable genes for phylogenetic studies. MSP1a, encoded by msp1
, has been reported to be an adhesin for bovine erythrocytes and tick cells and to effect adhesion, infection, and transmission of A. marginale by ticks of the genus Dermacentor (7, 8, 20, 21). Although the specific function of MSP4 is not known, our previous analysis of the msp4 gene from A. marginale isolates demonstrated that its sequence varies sufficiently among isolates to support its use in phylogeographic studies (6, 9).
Preliminary analysis of the phylogenetic relationships of A. marginale isolates on the basis of the MSP1a and MSP4 sequences demonstrated phylogeographic partitioning of A. marginale isolates in the United States (6). This information, together with the fact that for many isolates infection with geographic isolates is not cross protective, supports the inclusion of several isolates of A. marginale in the development of vaccine formulations for the United States (17). However, recent data support the existence of genetic heterogeneity in the structures of the msp1
sequences of isolates recovered from infected animals within a state (9) and even in a single herd of cattle in an area of endemicity (23), thus questioning the use of msp1
sequences to identify geographic isolates of A. marginale. Although single msp1
genotypes have been identified in individual cattle that were naturally or experimentally infected and sampled at different stages of infection (4, 23), a finding that we demonstrated could be explained by infection exclusion of A. marginale isolates in infected cattle (10), these findings suggest a rapid evolution of msp1
sequences, resulting in the maintenance of A. marginale genotypes in the population through different transmission events (9, 10, 23).
To test these hypotheses, we characterized the msp1
and msp4 genotypes of A. marginale isolates obtained during 2001 from infected cattle in east-central Oklahoma, where A. marginale is endemic, and conducted phylogenetic studies with two other isolates previously obtained from Oklahoma and isolates from Latin America and the United States.
|
|
|---|
sequence of an A. marginale isolate from Israel (GenBank accession no. AF352560) were used for outgroup comparisons. |
View this table: [in a new window] |
TABLE 1. Amplified sequences of A. marginale isolates from Oklahoma and their characteristicsa
|
![]() View larger version (32K): [in a new window] |
FIG. 1. Distribution by county of Oklahoma isolates of A. marginale used in the study.
|
and msp4 PCR and sequencing.
A. marginale DNA was extracted from 0.5 to 1 ml of erythrocytes prepared from blood collected in EDTA-treated tubes as reported previously (6). The msp1
gene was amplified from 1 µl (1 to 10 ng) of DNA by PCR with 10 pmol of each primer (primer MSP1
P [5'-GCATTACAACGCAACGCTTGAG-3'] and primer MSP1
3 [5'-GCTTTACGCCGCCGCCTGCGCC-3']) in a 50-µl volume (1.5 mM MgSO4, 0.2 mM deoxynucleoside triphosphates, 1x avian myeloblastosis virus-Tfl reaction buffer, 5 U of Tfl DNA polymerase) with the Access RT-PCR system (Promega, Madison, Wis.). The reactions were performed in an automated DNA thermal cycler (Mastercycler personal; Eppendorf, Westbury, N.Y.) for 35 cycles. After an initial denaturation step at 94°C for 30 s, each cycle consisted of denaturation at 94°C for 30 s and annealing-extension at 68°C for 2.5 min. The program ended by storage of the reaction mixtures at 4°C. The msp4 gene was amplified as described above but with oligonucleotides MSP45 (5'-GGGAGCTCCTATGAATTACAGAGAATTGTTTAC-3') and MSP43 (5'-CCGGATCCTTAGCTGAACAGGAATCTTGC-3') and a PCR profile of denaturation at 94°C for 30 s, annealing at 60°C for 30 s, and extension at 68°C for 1 min. The PCR products were electrophoresed on 1% agarose gels to check the sizes of the amplified fragments.
The amplified fragments were purified with resin (Wizard; Promega) and were cloned into the vector pGEM-T (Promega) or used directly for sequencing of both strands by double-stranded dye-termination cycle sequencing (Core Sequencing Facility, Department of Biochemistry and Molecular Biology, Noble Research Center, Oklahoma State University). The msp4 gene was completely sequenced. For the msp1
gene, only the fragment containing the upstream and variable regions was sequenced for all isolates (6, 9). When clones were recovered, at least two clones from each PCR were sequenced.
Sequence alignment and phylogenetic analysis.
The A. marginale msp1
gene fragment and the derived protein sequence containing the upstream and variable regions and the entire msp4-coding region were used for sequence alignment and phylogenetic analysis. Multiple-sequence alignment was performed with the program AlignX (Vector NTI Suite, version 5.5; InforMax, North Bethesda, Md.) with an engine based on the Clustal W algorithm (28). Nucleotides were coded as unordered, discrete characters with five possible character states (A, C, G, T, or N), and gaps were coded as missing data. Maximum-parsimony (MP) analyses were conducted by using equal weights for all characters and substitutions, heuristic searches with 25 random additions of input taxa, and tree bisection-reconnection (TBR) branch swapping by using the PAUP*4.0b4a program (27). The stability or the accuracy of the inferred topology(ies) was assessed by bootstrap analysis (14) of 200 iterations of heuristic searches with 25 random additions of input taxa and TBR branch swapping. Character-state changes were polarized by designating A. centrale and A. ovis as outgroups for the msp4 data set and by designating an A. marginale isolate from Israel as the outgroup for the msp1
and MSP1
data sets. To examine the effect of the analysis method on the resulting phylogenetic tree, we constructed a neighbor-joining phylogenetic tree (26) based on the distances obtained by maximum-likelihood (ML) analysis, with the model of DNA substitution chosen with MODELTEST (24).
To examine the phylogeographic patterns among the Oklahoma isolates of A. marginale, unrooted phylogenetic trees were constructed by ML analysis with the appropriate model of DNA sequence evolution, as determined with MODELTEST (24). The stabilities of the clades on the resulting msp4 and msp1
or MSP1a phylogenetic trees were estimated via 200 bootstrap iterations.
Selection pressure on msp1
and msp4.
The relative frequencies of nonsynonymous substitutions (dNs) and synonymous substitutions (dSs) among all pairwise comparisons of msp1
and msp4 sequences were estimated with MEGA (Molecular Evolutionary Genetics Analysis, version 1.01) software (18) by the method of Nei and Gojobori (22) and by application of the correction for multiple substitutions of Jukes and Cantor (15).
Nucleotide sequence accession numbers.
The GenBank accession numbers for the msp1
sequences of the Missouri and Oklahoma isolates of A. marginale are AY127052 to AY127064; and those for the msp4 sequences of the Oregon, Missouri, and Oklahoma A. marginale isolates are AY127065 to AY127078.
|
|
|---|
genes (512 to 861 bp; mean, 766 bp) and msp4 genes (851 bp) of 11 isolates of A. marginale collected from cattle in east-central Oklahoma during 2001. For comparative purposes, we included the DNA sequences of these two genes from two isolates of A. marginale previously collected in Wetumka, Okla., in 1997 (Oklahoma isolate) (5-10) and Pawhuska, Okla., in the 1960s (Table 1). All A. marginale isolates from Oklahoma possessed different msp1
and MSP1a sequences, with the number of tandem repeats varying from 1 (isolate Oklahoma City) to 5 (isolates Glencoe 1, Cushing 2, and Stillwater 1) (Fig. 2 and Table 1). Isolates of A. marginale from Glencoe and New Castle and one of the two isolates from Stillwater possessed identical msp4 sequences. Moreover, both isolates from Cushing possessed the same msp4 sequence as the isolate from Pawhuska that was collected in the 1960s. The remaining five isolates possessed unique msp4 sequences. All Oklahoma isolates of A. marginale had identical MSP4 sequences.
![]() View larger version (59K): [in a new window] |
FIG. 2. Comparison of the deduced amino acid sequences of the MSP1a variable region of Oklahoma isolates of A. marginale. The amino-terminal decapeptide and the beginning of the conserved region after the tandem repeats are underlined. The one-letter amino acid code is used. Asterisks indicate identical amino acids. Gaps indicate deletions or insertions.
|
and msp4 sequences of the Oklahoma A. marginale isolates, we included sequences from 7 Latin American isolates and 13 (for msp1
) or 11 (for msp4) U.S. isolates of A. marginale in subsequent analyses. The mean number of differences in the msp1
sequences for all pairwise comparisons of A. marginale isolates from Oklahoma was 14.5 (Table 2), whereas the mean number of differences in the msp1
sequences for all pairwise comparisons of the New World isolates examined was 19.7. In contrast, only seven unique msp4 sequences were detected among the isolates from Oklahoma, with the mean number of DNA substitutions for all pairwise comparisons being 3.8 (Table 2), which is only slightly less than the mean number of differences for all pairwise comparisons of New World isolates examined (n = 4.6). |
View this table: [in a new window] |
TABLE 2. Absolute number of DNA sequence differences of msp1 and msp4 for all pairwise comparisons of Oklahoma A. marginale isolates
|
(where GTR is general time reversible, I is the proportion of invariable sites, and
is the shape parameter of the gamma distribution) model of DNA sequence evolution as most appropriate. The model parameters used in the ML analysis that resulted in a single tree were as follows: base frequencies = 0.25, 0.25, and 0.25; NST (number of substitution types) = 6; Revmat (rate matrix) = 1.0000, 16.9376, 1.0000, 1.0000, and 1.0000; rates = gamma; alpha parameter of the gamma distribution (
) = 0.9197; and proportion of invariant sites (Pinvar) = 0.939 (Fig. 4). Strong bootstrap analysis support (73%) was detected for east-central (isolates Oklahoma, Stigler, Okmulgee, New Castle, and Glencoe 1 to 3) and north-central (isolates Oklahoma City, Cushing 1 and 2, and Pawhuska) clades within Oklahoma (Fig. 4). Both clades contained isolates from Stillwater (Fig. 4).
![]() View larger version (26K): [in a new window] |
FIG. 3. Topology of the strict consensus of 15 equally parsimonious trees of 51 steps (consistency index = 0.7843; retention index = 0.9209) based on unweighted MP analysis of msp4 sequence data. The numbers above the branches indicate the percent bootstrap analysis support for those clades supported in more than 50% of 200 bootstrap iterations with 25 random additions of input taxa and TBR branch swapping.
|
![]() View larger version (14K): [in a new window] |
FIG. 4. Unrooted phylogram obtained by ML analysis depicting the phylogenetic relationships among Oklahoma A. marginale isolates based on msp4 sequence data. The numbers above the branches indicate the percent bootstrap analysis support for those clades supported in more than 50% of 200 bootstrap iterations.
|
sequences demonstrated greater variability, which resulted in less phylogenetic and phylogeographic resolution than that afforded by use of msp4 sequences. The msp1
sequences of New World A. marginale isolates and one isolate from Israel provided 72 phylogenetically informative characters. MP analysis of these data resulted in 603 equally parsimonious trees of 163 steps (consistency index = 0.4908; retention index = 0.7552) and little support for phylogeographic clustering with the exception of a clade containing the isolates from Morelos and Veracruz, Mexico (Fig. 5). MODELTEST chose a modification of the GTR + I +
model of DNA substitution as the one that best fit the evolution of msp1
sequences of Oklahoma A. marginale isolates. The model parameters used in the ML analysis that resulted in a single tree were as follows: base frequencies = 0.2311, 0.2011, and 0.2986; NST = 6; Revmat = 1.0000, 7.8706, 1.0000, 1.0000, and 3.0680; rates = gamma;
= 0.7835; and Pinvar = 0.7001 (Fig. 6). In agreement with the results of MP analysis, ML analysis found little support for any phylogeographic structuring of the msp1
sequences of isolates from Oklahoma (Fig. 6). Similar results were obtained by MP and ML analyses of the MSP1a protein sequence data set (data not shown).
![]() View larger version (24K): [in a new window] |
FIG. 5. Topology of the strict consensus of 630 equally parsimonious trees of 163 steps (consistency index = 0.4908; retention index = 0.7552) based on unweighted MP analysis of msp1 sequence data. The numbers above the branches indicate the percent bootstrap analysis support for those clades supported in more than 50% of 200 bootstrap iterations with 25 random additions of input taxa and TBR branch swapping.
|
![]() View larger version (12K): [in a new window] |
FIG. 6. Unrooted phylogram obtained by ML analysis depicting the phylogenetic relationships among Oklahoma A. marginale isolates based on msp1 sequence data. The numbers above the branches indicate the percent bootstrap analysis support for those clades supported in more than 50% of 200 bootstrap iterations.
|
, dS was equal to 0.0286 ± 0.0157, whereas dN was equal to 0.0468 ± 0.0119, resulting in a dN/dS of 1.636. Thus, the selective pressures acting on msp1
are moderately positive and differ from those acting on msp4. |
|
|---|
genotypes circulate in some geographic regions (e.g., an area of endemicity in Oregon [23]), questioned the use of MSP1
sequences for identification of geographic isolates of A. marginale and suggested that msp1
sequences may be rapidly evolving (9). In this study we tested these hypotheses by analyzing the msp1
and msp4 genotypes of A. marginale isolates obtained during 2001from infected cattle in east-central Oklahoma, where A. marginale is endemic, and comparing them with those of New World isolates. As in previous studies (6, 9), msp4 sequences provided phylogeographic patterns for A. marginale isolates on a broad geographic scale. Additionally, DNA sequence characters further subdivided the U.S. isolates into west-central and southeastern clades. Both clades contained isolates from Oklahoma, suggesting extensive movement of cattle. Cattle movement during the 19th and early 20th centuries was responsible for the dissemination of tick-borne diseases in the United States (29, 30). On a smaller geographic scale with only the Oklahoma isolates, the msp4 sequences provided some phylogeographic information. However, the east-central and north-central clades included isolates from Stillwater, again suggesting the movement of cattle across this area. In contrast, analysis of the MSP1a DNA and protein sequences demonstrated extensive genotypic variation among Oklahoma isolates of A. marginale and did not provide phylogeographic information.
These results suggest that even if msp1
sequences are rapidly evolving, msp1
genotypes probably reflect the history of cattle movement rather than the geographic distribution of A. marginale isolates. Recent results suggest that different A. marginale genotypes are maintained within a herd in an area of endemicity by independent transmission events and that infection with more than one genotype per host is prevented by a mechanism of infection exclusion (10, 23). Therefore, if cattle movement introduces a new A. marginale genotype, it could be established by mechanical and/or biological transmission to susceptible cattle. In regions with few introductions of A. marginale isolates such as Australia, genotypic variation was found to be minimal (3, 19). In regions like Oklahoma, where the movement of cattle is extensive, a highly heterogeneous population of A. marginale isolates would be expected.
MSP1a is an immunodominant surface ligand involved in the interaction of A. marginale with host cell receptors (5, 7, 8). Furthermore, the tandem repeats that constitute the variable region of the protein mediate the interactions of the rickettsia with host cells and codify the capacity of the isolate to be transmitted by Dermacentor ticks (5, 8, 11). The support for positive selection pressure on msp1
is consistent with the biological function of MSP1a and suggests that this protein could be useful as a vaccine candidate. The sequence of MSP4 is highly conserved among A. marginale isolates and between A. marginale isolates and other Anaplasma spp. (9), but its biological function is unknown. Therefore, the lack of support for positive selection pressure on msp4 is difficult to interpret at present. Recently, Allsopp et al. (2) obtained similar results for the major antigenic protein (map1) of Ehrlichia ruminantium, which shows sequence identity with msp4, and suggested that map1 may not be a good candidate antigen for vaccination against E. ruminantium.
In summary, as we have demonstrated in previous studies (6, 9), msp4 sequences provided sufficient variation to detect the phylogeographic patterns of A. marginale isolates on different geographic scales. The results support the use of msp4 as a marker for elucidation of the phylogeographic patterns and the phylogenetic relationships among A. marginale isolates. In contrast to the results obtained with msp4, the DNA and deduced amino acid sequences of MSP1a failed to resolve either the phylogenetic or phylogeographic patterns of A. marginale isolates. The results indicate that MSP1a is not a marker for the characterization of geographic isolates of A. marginale and suggest that cattle movement and the maintenance of different genotypes by independent transmission events could explain the genetic heterogeneity observed among isolates of A. marginale in Oklahoma and, possibly, in other regions of endemicity.
Yu-Wei Chiang and Stephen Whu (Fort Dodge Animal Health, American Home Products, Fort Dodge, Iowa) and Katherine Curtis (University of Missouri) are acknowledged for providing the A. marginale Pawhuska and Missouri isolates, respectively. The contributions of Robert Mihalovich (Large Animal Resident, Oklahoma State University), field veterinarians, and cattle farmers in the collection of blood samples from infected cattle in Oklahoma are acknowledged. Dollie Clawson (Department of Veterinary Pathobiology, Oklahoma State University) is acknowledged for technical assistance. Sue Ann Hudiburg and Janet J. Rogers (Core Sequencing Facility, Department of Biochemistry and Molecular Biology, Noble Research Center, Oklahoma State University) are acknowledged for oligonucleotide synthesis and DNA sequencing, respectively.
|
|
|---|
polymerase chain reaction assay for specific detection and differentiation of Anaplasma marginale isolates. Vet. Microbiol. 86:325-335.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»