Antigen Diversity in the Parasitic Bacterium Anaplasma phagocytophilum Arises from Selectively-Represented,Spatially Clustered Functional Pseudogenes

Anaplasma phagocytophilum is a tick-transmitted bacterial pathogen of humans and other animals, and is an obligate intracellular parasite. Throughout the course of infection, hosts acquire temporary resistance to granulocytic anaplasmosis as they develop immunity specific for the major antigen, major surface protein 2 (Msp2). However, the bacterium then utilizes a novel recombination mechanism shuffling functional pseudogenes sequentially into an expression cassette with conserved 5′ and 3′ ends, bypassing host immunity. Approximately 100 pseudogenes are present in the only fully sequenced human-origin HZ genome, representing the possibility for almost unlimited antigenic diversity. In the present study, we identified a select group of 20% of the A. phagocytophilum HZ msp2 pseudogenes that have matched preferentially to human, canine, and equine expression cassettes. Pseudogenes cluster predominantly in one spatial run limited to a single genomic island in less than 50% of the genome but phylogenetically related pseudogenes are neither necessarily located in close proximity on the genome nor share similar percent identity with expression cassettes. Pseudogenes near the expression cassette (and the origin) are more likely to be expressed than those farther away. Taken together, these findings suggest that there may be natural selection pressure to retain pseudogenes in one cluster near the putative origin of replication, even though global recombination shuffles pseudogenes around the genome, separating pseudogenes that share genetic origins as well as those with similar identities.     

Conservation of a gene conversion mechanism in two distantly related paralogues of Anaplasma marginale

Anaplasmataceae, the causative agents of anaplasmosis and ehrlichiosis, persist in the bloodstream of their mammalian hosts, allowing acquisition and transmission by tick vectors. Anaplasma marginale establishes persistent infection characterized by sequential cycles of rickettsaemia in which new antigenic variants emerge. The two most immunodominant outer membrane proteins, MSP2 and MSP3, are paralogues, each encoded by a distinct family of related genes. This study demonstrates that, although the two gene families have diverged substantially, each has maintained a similar mechanism to generate structurally and antigenically polymorphic surface antigens. Like MSP2, MSP3 is expressed from a single locus in which variation of the expressed msp3 gene is generated by recombination using msp3 pseudogenes. Each of the msp3 pseudogenes encodes a unique central variable region (CVR) flanked by conserved 5' and 3' regions. Changes in the CVR of the expressed msp3, concomitant with invariance of the pseudogenes, indicate that expression site variation is generated using gene conversion. A. marginale thus maintains two large, separate systems within its small genome to generate antigenic variation of its surface proteins, while analogous structural elements indicate a common mechanism.     

Structural basis for segmental gene conversion in generation of Anaplasma marginale outer membrane protein variants

Bacterial pathogens in the genus Anaplasma generate surface coat variants by gene conversion of chromosomal pseudogenes into single-expression sites. These pseudogenes encode unique surface-exposed hypervariable regions flanked by conserved domains, which are identical to the expression site flanking domains. In addition, Anaplasma marginale generates variants by recombination of oligonucleotide segments derived from the pseudogenes into the existing expression site copy, resulting in a combinatorial increase in variant diversity. Using the A. marginale genome sequence to track the origin of sequences recombined into the msp2 expression site, we demonstrated that the complexity of the expressed msp2 increases during infection, reflecting a shift from recombination of the complete hypervariable region of a given pseudogene to complex mosaics with segments derived from hypervariable regions of different pseudogenes. Examination of the complete set of 1183 variants with segmental changes revealed that 99% could be explained by one of the recombination sites occurring in the conserved flanking domains and the other within the hypervariable region. Consequently, we propose an 'anchoring' model for segmental gene conversion whereby the conserved flanking sequences tightly align and anchor the expression site sequence to the pseudogene. Associated with the recombination sites were deletions, insertions and substitutions; however, these are a relatively minor contribution to variant generation as these occurred in less than 2% of the variants. Importantly, the anchoring model, which can account for more variants than a strict segmental sequence identity mechanism, is consistent with the number of msp2 variants predicted and empirically identified during persistent infection.     

Antigenic variation of Anaplasma marginale msp2 occurs by combinatorial gene conversion

The rickettsial pathogen Anaplasma marginale establishes lifelong persistent infection in the mammalian reservoir host, during which time immune escape variants continually arise in part because of variation in the expressed copy of the immunodominant outer membrane protein MSP2. A key question is how the small 1.2 Mb A. marginale genome generates sufficient variants to allow long-term persistence in an immunocompetent reservoir host. The recombination of whole pseudogenes into the single msp2 expression site has been previously identified as one method of generating variants, but is inadequate to generate the number of variants required for persistent infection. In the present study, we demonstrate that recombination of a whole pseudogene is followed by a second level of variation in which small segments of pseudogenes recombine into the expression site by gene conversion. Evidence for four short sequential changes in the hypervariable region of msp2 coupled with the identification of nine pseudogenes from a single strain of A. marginale provides for a combinatorial number of possible expressed MSP2 variants sufficient for lifelong persistence.     

The hypervariable region of Anaplasma marginale major surface protein 2 (MSP2) contains multiple immunodominant CD4+ T lymphocyte epitopes that elicit variant-specific proliferative and IFN-gamma responses in MSP2 vaccinates

Major surface protein 2 (MSP2) is an immunodominant outer membrane protein of Anaplasma marginale and Anaplasma phagocytophilum pathogens that cause bovine anaplasmosis and human granulocytic ehrlichiosis, respectively. MSP2 has a central hypervariable region (HVR) flanked by highly conserved amino and carboxyl termini. During A. marginale infection, dynamic and extensive amino acid sequence variation in MSP2 occurs through recombination of msp2 pseudogenes into the msp2 expression site, followed by sequential segmental gene conversions to generate additional variants. We hypothesized that MSP2 variation leads to significant changes in Th cell recognition of epitopes in the HVR. T cell epitopes were mapped using T cells from native MSP2-immunized cattle and overlapping peptides spanning the most abundant of five different MSP2 HVRs in the immunogen. Several epitopes elicited potent effector/memory Th cell proliferative and IFN-gamma responses, including those in three discreet blocks of sequence that undergo segmental gene conversion. Th cell clones specific for an epitope in the block 1 region of the predominant MSP2 variant type failed to respond to naturally occurring variants. However, some of these variants were recognized by oligoclonal T cell lines from MSP2 vaccinates, indicating that the variant sequences contain immunogenic CD4(+) T cell epitopes. In competition/antagonism assays, the nonstimulatory variants were not inhibitory for CD4(+) T cells specific for the agonist peptide. Dynamic amino acid sequence variation in MSP2 results in escape from recognition by some effector/memory MSP2-specific Th cells. Antigenic variation in MSP2 Th cell and B cell epitopes may contribute to immune evasion that allows long-term persistence of A. marginale in the mammalian reservoir.     

Anaplasma phagocytophilum in small mammals and ticks in northeast Florida

Human anaplasmosis is an emerging tick-borne disease in the United States, but few studies of the causative agent, Anaplasma phagocytophilum, have been conducted in southeastern states. The aim of this study was to determine if A. phagocytophilum is present in small mammals and ticks in northeast Florida. Polymerase chain reaction assays designed to amplify portions of the major surface protein 2 gene (p44), 16S rDNA, and groESL operons were used to test rodent blood and tick DNA samples for the presence of A. phagocytophilum. Positive samples were confirmed by DNA sequence analysis. Anaplasma phagocytophilum DNA was detected in less than 5% of cotton mice and 45% of cotton rats from two sites in northeast Florida. Anaplasma phagocytophilum DNA was also confirmed in 1.3% of host-seeking adult Ixodes scapularis tested and 2.7% of host-seeking adult Amblyomma americanum. This report describes the first DNA sequence data confirming strains of A. phagocytophilum in rodents and ticks in Florida. The DNA sequences of the msp2, 16S rDNA, and groESL gene fragments obtained in this study were highly similar to reference strains of human pathogenic strains of A. phagocytophilum. These findings suggest that A. phagocytophilum is present and established among some small mammal species in northeast Florida. Although the infection prevalence was low in the total number of ticks tested, the presence of A. phagocytophilum in two human biting tick species, one of which is a known competent vector, suggests that humans in this region may be at risk of granulocytic anaplasmosis caused by this pathogen.     

PCR detection of Anaplasma phagocytophilum in goat flocks in an area endemic for tick-borne fever in Switzerland

Central Switzerland is a highly endemic region for tick-borne fever (TBF) in cattle, however, little is known about A. phagocytophilum in goats. In the present study, 72 animals from six goat flocks (373 EDTA blood-samples) in Central Switzerland were analysed for A. phagocytophilum DNA. A real-time PCR targeting the msp2 gene of A. phagocytophilum was performed and in positive samples the partial 165 rRNA, groEL and msp4 gene were amplified for sequence analysis. Four DNA extracts were positive. Different sequence types on basis of the amplified genes were found. For comparison, sequences of A. phagocytophilum from 12 cattle (originating from Switzerland and Southern Germany) were analysed. The 165 rRNA gene sequences from cattle were all identical amongst each other, but the groEL and msp4 gene differed depending on the origin of the cattle samples and differed from the variants from goats. This study clearly provides molecular evidence for the presence of different types of A. phagocytophilum in goat flocks in Switzerland, a fact which deserves more thorough attention in clinical studies.     

Coexistence of emerging bacterial pathogens in Ixodes ricinus ticks in Serbia

The list of tick-borne pathogens is long, varied and includes viruses, bacteria, protozoa and nematodes. As all of these agents can exist in ticks, their co-infections have been previously reported. We studied co-infections of emerging bacterial pathogens (Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Francisella tularensis) in Ixodes ricinus ticks in Serbia. Using PCR technique, we detected species-specific sequences, rrf-rrl rDNA intergenic spacer for B. burgdorferi s.l., p44/msp2 paralogs for A. phagocytophilum, and the 17 kDa lipoprotein gene, TUL4, for F. tularensis, respectively, in total DNA extracted from the ticks. Common infections with more than one pathogen were detected in 42 (28.8%) of 146 infected I. ricinus ticks. Co-infections with two pathogens were present in 39 (26.7%) of infected ticks. Simultaneous presence of A. phagocytophilum and different genospecies of B. burgdorferi s.l. complex was recorded in 16 ticks, co-infection with different B. burgdorferi s. l. genospecies was found in 15 ticks and eight ticks harbored mixed infections with F. tularensis and B. burgdorferi s.l. genospecies. Less common were triple pathogen species infections, detected in three ticks, one infected with A. phagocytophilum / B. burgdorferi s.s. / B. lusitaniae and two infected with F. tularensis / B. burgdorferi s.s. / B. lusitaniae. No mixed infections of A. phagocytophilum and F. tularensis were detected.     

Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum, Francisella tularensis and their co-infections in host-seeking Ixodes ricinus ticks collected in Serbia

To evaluate the prevalence rate of tick-borne bacterial pathogens, unfed adult Ixodes ricinus ticks were collected from vegetation in 2001, 2003, and 2004 at 18 localities throughout Serbia. A total of 287 ticks were examined by PCR technique for the presence of Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum, and Francisella tularensis. The highest prevalence rate was that for B. burgdorferi sensu lato (42.5%), followed by A. phagocytophilum (13.9%) and F. tularensis (3.8%). The presence of five B. burgdorferi sensu lato genospecies, namely, B. burgdorferi sensu stricto, B. afzelii, B. garinii, B. lusitaniae, and B. valaisiana was identified by restriction fragment length polymorphism (RFLP) analysis. The most frequent B. burgdorferi sensu lato genospecies was B. lusitaniae, followed by B. burgdorferi sensu stricto. Co-infection by B. burgdorferi sensu stricto and B. lusitaniae was frequently observed. Co-infection by B. burgdorferi sensu lato and A. phagocytophilum and co-infection by B. burgdorferi sensu lato and F. tularensis appeared in 24 ticks. Sequencing of p44/msp2 paralogs of Serbian A. phagocytophilum showed that they were unique and distinct from those of A. phagocytophilum in US and UK. This is the first report of B. garinii, B. lusitaniae, B. valaisiana, as well as A. phagocytophilum and F. tularensis infected ticks in Serbia. These findings indicate a public health threat in Serbia of tick-borne diseases caused by B. burgdorferi sensu lato, A. phagocytophilum and F. tularensis.     

The influence of population structure on gene expression and flowering time variation in the ubiquitous weed Capsella bursa‐pastoris (Brassicaceae)

Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa‐pastoris. We used RNA‐Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping‐by‐sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.     

Cloning and modeling of CD8 beta in the amphibian ambystoma Mexicanum. Evolutionary conserved structures for interactions with major histocompatibility complex (MHC) class I molecules

Bovine anaplasmosis is a rickettsial disease of world-wide economic importance caused by Anaplasma marginale. Several major surface proteins with conserved gene sequences have been examined as potential candidates for vaccines and/or diagnostic assays. Major surface protein 1 (MSP1) is composed of polypeptides MSP1a and MSP1b. MSP1a is expressed from the single copy gene msp1 alpha and MSP1b is expressed by members of the msp1 beta multigene family. In order to determine if the msp1 genes are conserved, primers specific for msp1 alpha, msp1 beta(1), and msp1 beta(2) genes were synthesized and used to amplify msp1 sequences of A. marginale from tick cell cultures, from cattle during acute and chronic infections and from salivary glands of Dermacentor variabilis. Protein sequences of MSP1a, MSP1b(1) and MSP1b(2) were conserved during the life cycle of the parasite. No amino acid changes were observed in MSP1a. However, small variations were observed in the MSP1b(1) and MSP1b(2) protein sequences, which could be attributed to recombination, selection for sub-populations of A. marginale in the vertebrate host and/or PCR errors. Several isolate-specific sequences were also observed. Based on the information obtained in this study, the MSP1 protein appears to be fairly well conserved and a potential vaccine candidate.     

PERSPECTIVE: THE THEORIES OF FISHER AND WRIGHT IN THE CONTEXT OF METAPOPULATIONS: WHEN NATURE DOES MANY SMALL EXPERIMENTS

We critically review the two major theories of adaptive evolution developed early in this century, Wright's shifting balance theory and Fisher's large population size theory, in light of novel findings from field observations, laboratory experiments, and theoretical research conducted over the past 15 years. Ecological studies of metapopulations have established that the processes of local extinction and colonization of demes are relatively common in natural populations of many species and theoretical population genetic models have shown that these ecological processes have genetic consequences within and among local demes. Within demes, random genetic drift converts nonadditive genetic variance into additive genetic variance, increasing, rather than limiting, the potential for adaptation to local environments. For this reason, the genetic differences that arise by drift among demes, can be augmented by local selection. The resulting adaptive differences in gene combinations potentially contribute to the genetic origin of new species. These and other recent findings were not discussed by either Wright or Fisher. For example, although Wright emphasized epistatic genetic variance, he did not discuss the conversion process. Similarly, Fisher did not discuss how the average additive effect of a gene varies among demes across a metapopulation whenever there is epistasis. We discuss the implications of such recent findings for the Wright-Fisher controversy and identify some critical open questions that require additional empirical and theoretical study.     

Allelic dimorphism-associated restriction of recombination in Plasmodium falciparum msp1

Allelic dimorphism is a characteristic feature of the Plasmodium falciparum msp1 gene encoding the merozoite surface protein 1, a strong malaria vaccine candidate. Meiotic recombination is a major mechanism for the generation of msp1 allelic diversity. Potential recombination sites have previously been mapped to specific regions within msp1 (a 5' 1-kb region and a 3' 0.4-kb region) with no evidence for recombination events in a central 3.5-kb region. However, evidence for the lack of recombination events is circumstantial and inconclusive because the number of msp1 sequences analysed is limited, and the frequency of recombination events has not been addressed previously in a high transmission area, where the frequency of meiotic recombination is expected to be high. In the present study, we have mapped potential allelic recombination sites in 34 full-length msp1 sequences, including 24 new sequences, from various geographic origins. We also investigated recombination events in blocks 6 to 16 by population genetic analysis of P. falciparum populations in Tanzania, where malaria transmission is intense. The results clearly provide no evidence of recombination events occurring between the two major msp1 allelic types, K1-type and Mad20-type, in the central region, but do show recombination events occurring throughout the entire gene within sequences of the Mad20-type. Thus, the present study indicates that allelic dimorphism of msp1 greatly affects inter-allelic recombination events, highlighting a unique feature of allelic diversity of P. falciparum msp1.     

Evaluating the role of natural selection in the evolution of gene regulation

Surveys of gene expression reveal extensive variability both within and between a wide range of species. Compelling cases have been made for adaptive changes in gene regulation, but the proportion of expression divergence attributable to natural selection remains unclear. Distinguishing adaptive changes driven by positive selection from neutral divergence resulting from mutation and genetic drift is critical for understanding the evolution of gene expression. Here, we review the various methods that have been used to test for signs of selection in genomic expression data. We also discuss properties of regulatory systems relevant to neutral models of gene expression. Despite some potential caveats, published studies provide considerable evidence for adaptive changes in gene expression. Future challenges for studies of regulatory evolution will be to quantify the frequency of adaptive changes, identify the genetic basis of expression divergence and associate changes in gene expression with specific organismal phenotypes.     

Differential expression of VirB9 and VirB6 during the life cycle of Anaplasma phagocytophilum in human leucocytes is associated with differential binding and avoidance of lysosome pathway

Anaplasma phagocytophilum, an obligate intracellular bacterium, is the aetiologic agent of human granulocytic anaplasmosis (HGA). A. phagocytophilum virB/D operons encoding type IV secretion system are expressed in cell culture and in the blood of HGA patients. In the present study, their expression across the A. phagocytophilum intracellular developmental cycle was investigated. We found that mRNA levels of both virB9 and virB6 were upregulated during infection of human neutrophils in vitro. The antibody against the recombinant VirB9 protein was prepared and immunogold and immunofluorescence labelling were used to determine the VirB9 protein expression by individual organisms. Majority of A. phagocytophilum spontaneously released from the infected host cells poorly expressed VirB9. At 1 h post infection, VirB9 was not detectable on most bacteria associated with neutrophils. However, VirB9 was strongly expressed by A. phagocytophilum during proliferation in neutrophils. In contrast, with HL-60 cells, approximately 80% of A. phagocytophilum organisms associated at 1 h post infection expressed VirB9 protein and were colocalized with lysosome-associated membrane protein-1 (LAMP-1), whereas, VirB9-undetectable bacteria were not colocalized with LAMP-1. These results indicate developmental regulation of expression of components of type IV secretion system during A. phagocytophilum intracellular life cycle and suggest that bacterial developmental stages influence the nature of binding to the hosts and early avoidance of late endosome-lysosome pathway.     

Swimming against the current: genetic structure, host mobility and the drift paradox in trematode parasites

Life-cycle characteristics and habitat processes can potentially interact to determine gene flow and genetic structuring of parasitic species. In this comparative study, we analysed the genetic structure of two freshwater trematode species with different life histories using cytochrome c oxidase I gene (COI) sequences and examined the effect of a unidirectional river current on their genetic diversity at 10 sites along the river. We found moderate genetic structure consistent with an isolation-by-distance pattern among subpopulations of Coitocaecum parvum but not in Stegodexamene anguillae. These contrasting parasite population structures were consistent with the relative dispersal abilities of their most mobile hosts (i.e. their definitive hosts). Genetic diversity decreased, as a likely consequence of unidirectional river flow, with increasing distance upstream in C. parvum, which utilizes a definitive host with only restricted mobility. The absence of such a pattern in S. anguillae suggests that unidirectional river flow affects parasite species differently depending on the dispersal abilities of their most mobile host. In conclusion, genetic structure, genetic diversity loss and drift are stronger in parasites whose most mobile hosts have low dispersal abilities and small home ranges. An additional prediction can be made for parasites under unidirectional drift: those parasites that stay longer in their benthic intermediate host or have more than one benthic intermediate hosts would have relatively high local recruitment and hence increased retention of upstream genetic diversity.