Attenuation of virulence in an apicomplexan hemoparasite results in reduced genome diversity at the population level
© Lau et al; licensee BioMed Central Ltd. 2011
Received: 16 February 2011
Accepted: 12 August 2011
Published: 12 August 2011
Virulence acquisition and loss is a dynamic adaptation of pathogens to thrive in changing milieus. We investigated the mechanisms of virulence loss at the whole genome level using Babesia bovis as a model apicomplexan in which genetically related attenuated parasites can be reliably derived from virulent parental strains in the natural host. We expected virulence loss to be accompanied by consistent changes at the gene level, and that such changes would be shared among attenuated parasites of diverse geographic and genetic background.
Surprisingly, while single nucleotide polymorphisms in 14 genes distinguished all attenuated parasites from their virulent parental strains, all non-synonymous changes resulted in no deleterious amino acid modification that could consistently be associated with attenuation (or virulence) in this hemoparasite. Interestingly, however, attenuation significantly reduced the overall population's genome diversity with 81% of base pairs shared among attenuated strains, compared to only 60% of base pairs common among virulent parental parasites. There were significantly fewer genes that were unique to their geographical origins among the attenuated parasites, resulting in a simplified population structure among the attenuated strains.
This simplified structure includes reduced diversity of the variant erythrocyte surface 1 (ves) multigene family repertoire among attenuated parasites when compared to virulent parental strains, possibly suggesting that overall variance in large protein families such as Variant Erythrocyte Surface Antigens has a critical role in expression of the virulence phenotype. In addition, the results suggest that virulence (or attenuation) mechanisms may not be shared among all populations of parasites at the gene level, but instead may reflect expansion or contraction of the population structure in response to shifting milieus.
Pathogens adapt to maintain selective advantage in their environments. As few environments are themselves stable, pathogen adaptation is a dynamic and continuous process. This principle applies to virulence in which acquisition and loss of virulence is dynamic within a pathogen population, varying with host genetics, host immune status at the individual and population levels, and transmission efficiency. These shifts can be achieved through sexual reproduction where novel recombinations often lead to genome diversity . However, for many multi-stage pathogens, this cannot occur during haploid life stages, and thus, they depend on mutations at the genome level or adopt phenotypic variation as a result of non-mutational mechanisms. Mutations associated with virulence (or attenuation) are frequently documented in eubacteria [2–4] and viruses [5, 6]. However, there is a significant gap in our knowledge of virulence-associated gene mutations in protozoa, including pathogens which have a major impact on global public and animal health. With increased emphasis on development and delivery of attenuated vaccines for hemoparasites such as Babesia, Theileria, and Plasmodium spp. [7–9], the ability to predictably and stably attenuate these pathogens would be a significant step toward vaccine implementation.
In this study, we investigated a natural host-apicomplexan pathogen interaction, infection of cattle with Babesia bovis, to investigate virulence loss at the genome level using multiple strains of diverse geographical origin and genetic background. Virulent B. bovis infection results in a clinical syndrome characterized by severe hemolytic anemia, hyperthermia, and a syndrome of neurovirulence clinically and pathophysiologically similar to cerebral malaria . However, through serial in vivo passages of virulent B. bovis in a splenectomized host, the virulence phenotype is gradually lost, and an attenuated B. bovis derivative is obtained. Animals infected with the attenuated derivative do not exhibit neurovirulence, do not require treatment, and are protected upon virulent parental challenge . Reliable derivation of attenuated parasites is the basis for production of B. bovis live vaccine in Babesia endemic countries such as Australia, Israel and Argentina [9, 12, 13]. Furthermore, attenuated vaccine generated from strains isolated in Australia confer cross protection in cattle in Africa, Latin America and Southeast Asia [11, 13], suggesting that immunogens are shared between geographically divergent strains.
In this report, we took advantage of the babesial attenuation procedure to compare coding regions across the complete genomes of distinct B. bovis strain pairs with distinguishable virulence and attenuation phenotypes in the natural host. With one virulent B. bovis genome sequence already published , we subsequently sequenced its attenuated derivative and two additional geographically distinct strain pairs using pyrosequencing technology. The results of a comprehensive genome analysis among all three virulent/attenuated strain pairs (six strains) are presented and discussed in the context of dynamic virulence acquisition and loss.
Three B. bovis virulent and attenuated strain pairs were used in this study. T2Bo_vir is a strain previously isolated in Texas, USA and its genome sequence has been published . T2Bo_att is the attenuated derivative of T2Bo_vir. L17_vir is a virulent field isolate that originated from an infected animal in Salta province, Northwest Argentina. L17_att is the attenuated derivative of L17_vir. T_vir, also referred to as E61, is a virulent isolate that originated from an infected animal in North Queensland, Australia. T_att, also referred to as F100, is the attenuated derivative of T_vir. For comparative analysis, the fully sequenced T2Bo_vir is used as the reference and all the strains used are designated as "test" strains.
Attenuation of virulent B. bovis strains
Attenuation of the T strain (T_att) was achieved through 28 in vivo rapid passages from its original virulent parent, T_vir as previous described . The subsequent attenuated T strain, T_att, was used as a vaccine in Australia between 1991 and 1993 .
The L17_att and T2Bo_att strains were also attenuated by rapid passage as described  and were evaluated after 24 and 29 in vivo passages, respectively. Before the attenuation procedure, T2Bo_vir and L17_vir strains were passed through non-infected Rhipicephalus microplus ticks to obtain stabilates. Liquid nitrogen preserved stabilates of L17_vir or T2Bo_vir were used to inoculate the first splenectomized calves with 1.2 × 108 parasitized erythrocytes after which the infected blood was collected in heparin (5-10 IU/ml) five days post-infection and transferred intravenously to a second splenectomized calf. Subsequent passages were identical with blood collected and transferred with an average of 5 ± 1 day interval.
Evaluation of attenuated derivatives
In vivo evaluation of Babesia bovis attenuation
# of animals
Average age of animal (mth)
# of animals requiring treatment*
Short term in vitro propagation of B. bovis and genomic DNA extraction
T2Bo_vir or _att and L17_vir or _att strains from splenectomized calves were used to establish short term continuous in vitro propagation (< 30 days of culture) as described previously [16, 17]. Parasitized erythrocytes of each strain were frozen in liquid nitrogen, cryopreserved with 10% polyvinylpyrolidone (PVP) and stored in liquid nitrogen until use.
T_vir was made from pooled blood from 10 animals from a property called 'Reedy Beds' Tharingowa via Townsville in far North Queensland. This Australian virulent strain was passaged twice in splenectomized calves to generate enough parasitized erythrocytes for gDNA extraction. Short term in vitro culturing was not performed for this strain.
Genomic DNA of B. bovis infected erythrocytes was obtained as previously described . Southern dot blot hybridization was carried on test strain gDNA to ensure that all samples contained < 10% bovine host gDNA contamination (data not shown) prior to pyrosequencing.
The B. bovis gDNA for T2Bo_vir was completely sequenced by a classical method (Sanger shotgun) and contains two small assembly gaps and one physical gap . B. bovis gDNA for the remaining five strains was submitted to the 454 pyrosequencing facility (http://www.454.com) for shotgun sequencing (see Additional file 1, table s1). Briefly, genomic DNA was randomly fragmented into small 300 to 800 bp fragments. Adapters were ligated to the generated fragments, thus creating a library of DNA fragments. These fragments were immobilized onto DNA capture beads and individually sequenced on a PicoTiterPlate device. The average read lengths of T_att, T_vir, L17_att, L17_vir and T2Bo_att were 101 bp, 251 bp, 215 bp, 250 bp and 386 bp, respectively. The generated sequences were then assembled into a number of unordered and unoriented contigs using the GS De Novo Assembler Software resulting in a consensus sequence. For more information on pyrosequencing and the assembly tool, please visit http://www.454.com/applications.
Comparative genomic analyses
Comparative genomics analysis was performed on all six strains using the T2Bo_vir strain as the reference strain and the remaining strains as test strains. Blastn was used to assign contigs to chromosome scaffolds (chromosomes 1-4, apicoplastic, mitochondrial) based on the topmost hits. Ordering and orientation of the contigs in the draft assemblies was achieved by mapping along the T2Bo_vir reference genome. Scaffolded contigs were visualized for manual inspection using the Mauve software package [20, 21].
Inter-strain genome comparisons involving two or more genomes from either within or across the virulent and attenuated groups were carried out using the Mauve tool. These genome alignments also served to visualize and quantify the overall quality of the individual draft assemblies for the test strains. A combination of Unix bash and Perl scripts were written to compute genome- and gene-wide coverage of the individual draft assemblies relative to the reference, and to calculate various other measures that relate to the distribution of assembly gaps and the impact of repetitive content (e.g., ves genes) on assembly. These scripts were written specifically for these genomes and will be made freely available upon request.
SNP identification and gene analysis
Initially, gene boundaries were marked along the T2Bo_vir contigs using available GenBank annotation. Subsequently, the subset of T2Bo_att contigs that mapped to the genic portions of T2Bo_vir were identified and isolated. The coordinates for this mapping information were determined using MUMmer program . This subset was further sub-divided into gene groups (≤ e value of 10-5 as an arbitrary cutoff for gene determination), where contigs mapping to the same gene were binned together. Long contigs mapping to multiple genes, although from different parts, were made as part of each of the corresponding gene groups. After binning, a multiple sequence alignment (MSA) was generated separately for each group using the CAP3 fragment assembly tool . Note that using an MSA would restrict the similarity detection process to full-length sequence at the nucleotide level. This highly sensitive strategy was used deliberately as a means to flag the slightest of differences between the two forms of T2Bo. Automated scripts were written to trim the MSAs of their leading and trailing gaps, following which intronic portions were excised, retaining only the coding segments. Output of MSAs was then manually examined by visual inspection using an alignment program (MacVector, vers.11.1). Two types of polymorphic events were identified during the search process - base substitutions such as single nucleotide polymorphisms (SNPs) and sequence gaps. The latter event could have arisen due to true insertions/deletions (indels) or simply lack of sequence data. A list of candidate genes resulting from the two-way (T2Bo_vir and _att only) comparison was generated for further evaluations. MSA of these genes were subsequently generated for the remaining two strain pairs. Genes that were commonly found to be different between all strain pairs were then selected for re-sequencing.
Separately, SNPs and indels in the different test strains relative to the reference were identified by comparing each of the test strains against the reference using the MUMmer tool and its default parameters and a customized-written script (freely available for download upon request).
Re-sequencing and analysis
Gene specific primers were generated and synthesized (http://idtdna.com) for re-sequencing of the final candidate genes. PCR amplicons of genomic DNA from T2Bo virulent and attenuated strains were used in the comparison to confirm SNPs or indels that were identified in silico between this strain pair. If polymorphisms were detected from this comparison, genomic DNA from L17 and T strain pairs were then used for re-sequencing to determine if the observed polymorphisms were also found in the other strains; thereby making the process into a 6-way comparison. High fidelity Taq (Invitrogen) was used in all re-sequencing to minimize PCR-generated errors. Amplicons were cloned into pCR4.1TOPO (Invitrogen) and multiple clones per PCR were sequenced (WSU Sequencing Core Facility). Sequence analysis was performed using MacVector vers.11.1. All putative protein-coding sequences were included in the analysis except for ves1 gene family (n = 119), BBOV_I002740 and BBOV_IV029790 which encode for variant erythrocyte surface antigens, a 200 kDa antigen and a hypothetical protein, respectively. Repetitive sequences within these genes prohibited accurate alignments for comparison.
Pyrosequencing achieved high genomic coverage of all strains
Assembly statistics of the five Babesia bovis genomes in comparison to the reference T2Bo_vir
Number of contigs
N50 contig length (bp)
Longest contig (bp)
Total contig length (bp)
Genomic and genic coverages of the five Babesia bovis draft assemblies over the reference T2Bo_vir
A. Genome-wide coverage
% base coverage by contigs
Total gap length (bp)
Number of gaps
B. Genic coverage
% base coverage by contigs
Gap length in genes (bp)
Genic coverage by chromosome of the five Babesia bovis strains over the reference T2Bo_vir
(Average % coverage)
Sequence divergence decreases with attenuation
No common gene differences are shared among all virulent or attenuated strain pairs
In order to identify common gene differences between all virulent and all attenuated strains, MSA was performed initially between T2Bo_vir reference and T2Bo_att to generate a list of genes with either nucleotide differences or simply the absence of any MSA output between this parent-offspring pair. A total of 226 genes were identified for further investigation (data not shown). Members of this gene list were used to generate MSA between L17_vir and L17_att and T_vir and T_att. A second gene list comprising gene differences shared among the three strain pairs resulted (Additional file 1, table s3). Out of these 14 genes, half of these have annotational information while the remaining genes encode for hypothetical proteins. However, re-sequencing of all 14 genes showed that although common polymorphisms were found between all attenuated/virulent strain pairs, these polymorphisms were insignificant in the translated products using SIFT .
The ability to predictably attenuate B. bovis by multiple passages through splenectomized cattle has been used in Babesia endemic regions to produce the only available vaccine that can protect animals from severe clinical disease [11, 25]. The attenuated parasites are infectious but do not cause clinical disease requiring treatment. A similar loss of virulence has been reported in related organisms such as Theileria, Toxoplasma and Plasmodium[26–28]. We took advantage of this attenuation model system to perform whole genome comparison of three different B. bovis virulent parent and attenuated daughter strain pairs. Our initial focus was to identify gene changes shared among all attenuated strains in an effort to determine if there was a common pathway to attenuation.
Pyrosequencing technology  enabled us to achieve high genome coverage for all the B. bovis test strains (n = 5), with the highest coverage for T2Bo_att, followed by L17 and T strains (Table 3). The quality of the assembly varied from strain to strain. Most notably the N50 contig length was considerably less for the virulent strains (L17_vir and T_vir) than their attenuated counterparts (L17_att and T_att) (Table 2), indicating that the assembly for the virulent forms of these two strains was more fragmented despite the higher depth used in sequencing the virulent parasites (see Additional file 1, table s1). This result is most likely a manifestation of the high level of sequence diversity among the virulent population - i.e., reads emerging from a diverse population are less likely to collapse into long contigs during assembly. Despite the difficulties in assembly encountered in virulent strains, the study was not negatively impacted since genomic coverage of these strains was still substantially high (> 84%) when aligned over the T2Bo_vir reference genome.
Comparison of shared sequences between the three strain pairs showed that there are more shared sequences between the North and South American strains (L17 and T2Bo) (Figure 5). This suggests that L17 and T2Bo are genetically more closely related to each other than to the Australian T strain. More SNPs were identified in both T strains when the orthologues were compared to those in L17 and T2Bo. Many of these SNPs, nonetheless, were synonymous changes. However, this observation suggests that the codon usage is similar between L17 and T2Bo. The greater level of T strain sequence diversity compared to the reference strain may contribute to increased difficulty in recognizing orthologues to the reference genes using our stringent gene identification criteria, subsequently resulting in more and larger sequence gaps.
When genomic- and genic-wide coverages by chromosome (Table 3) in all strains were computed, irrespective of the quality of the coverage, chromosomes 1 and 4 consistently had the poorest and best coverage, respectively. A noticeable difference between these two chromosomes is the distribution of variant erythrocyte surface 1 (ves) gene family members. The ves gene family is the largest gene family in the B. bovis (T2Bo strain) genome and contains multiple repeat sequences [10, 14]. Chromosome 1 has the highest number of known ves, with at least 42 of the 119 total ves in the genome, representing ~16% of the total gene content of chromosome 1. By contrast, only ~3% of the total gene content of chromosome 4 is comprised of ves. Since pyrosequencing is not typically suited to repetitive regions , this could explain the lower gene coverage observed in the assembly of chromosome 1 when compared to chromosome 4 (Table 4). To validate this hypothesis, we mapped all the ves locations in the two chromosomes onto the contigs and investigated if they were located within the gap regions. First, we observed that a high percentage (> 84%) of the ves gene bases fall into assembly gaps consistently across all chromosomes, reflecting the high complexity inherent in sequencing the ves gene family. Next, we compared the distribution of the assembly gaps to ves and non-ves regions in chromosome 1 and 4. Analysis showed that 46% of the gaps in chromosome 1 are attributable to ves genes while the remaining 54% include non-ves genes, intergenic and intronic regions. In contrast, only 21% of the gaps in chromosome 4 are attributable to ves (data not shown). These data support our conclusion that the poor coverage of chromosome 1 in the virulent strains as compared to chromosome 4 is due to the clustering of ves.
The total number of strain-specific sequences in each virulent strain is greater than that of the attenuated strains (Figure 5), explaining in part why there are more gaps in the consensus chromosomal alignments among the virulent strains (Figure 2A). Approximately 60% of the sequences are shared among all virulent strains and 81% among the attenuated strains. This is not a simple artifact of comparing sequences obtained by different techniques, as comparison of the L17_vir and T_vir strains, both of which were pyrosequenced, resulted in the same significant gap differences when compared to the attenuated strains. This reduced diversity in strain-specific sequences may be attributed to a change in population structure during the attenuation process, resulting in a simplified or more uniform attenuated population. Evidence of a similar phenomenon has been reported in Theileria and Plasmodium[28, 31]. However, our study is the first demonstration of the change in population structure at the whole genome level. Among the strain-specific gene sequences, ves predominates in both the virulent and attenuated strains (Figure 6A). Interestingly, 97% of the total strain-specific attenuation-associated ves genes are also present in at least one virulent strain (Figure 6B), suggesting that the attenuation process did not select either or against a specific repertoire or for a unique repertoire of variant genes. Due to their polymorphic nature and repetitive sequences as well as a physical gap of approximately 150 kbp containing ves genes in chromosome 1 of the T2Bo_vir genome sequence, members of this gene family were not individually analyzed in our study. Nonetheless, results from the analysis of the strain-specific ves gene sequences between virulent and attenuated strains imply that there is likely considerable overlap of the known ves repertoire between the virulent parent and its attenuated derivative. This may account for the ability of the attenuated parasites to protect against clinical disease caused by virulent strains after only one infection cycle. In contrast, multiple clinical infections are typically required to generate clinical immune protection against P. falciparum infection [32, 33], possibly reflecting the more diverse repertoire of variant genes such as var (encoding PfEMP1) among virulent isolates .
Analysis of non-ves1 strain-specific gene sequences between virulent and attenuated strains also showed a significant decrease in the number of strain-specific sequences as a result of attenuation. Approximately 60% of these strain-specific non-ves sequences are of unknown function while the remainder in both virulent and attenuated groups includes membrane proteins, merozoite surface antigen-2 and s mall o pen r eading f rames (smorfs) (Figure 7). Non-ves strain-specific sequences in virulent strains also include genes not found in the attenuated non-ves unique sequences group. This suggests that some genes uniquely present in the virulent strains are selected against during the attenuation process, while others are enriched.
There were few non-synonymous differences between coding sequences in all three strain pairs which might account for the shared phenotypic characteristics of virulence or attenuation as all non-synonymous changes resulted in tolerated translated products. This may imply that virulence may be multigenic, involving more than one gene in different strains and depending on the strain, may involve different genes, all belonging to the same pathway. Attenuation may also result from selection of a subpopulation whose virulence factors fail to reach a threshold expression level for virulence, such a mechanism has been reported by Colinet et al.  for Leptopilina boulardi, a parasitic wasp of Drosophila melanogaster.
Our comparative data indicate that ves gene members are repeatedly the dominant strain-specific sequences for each strain. Ves is 70% and 86% of the total strain-specific sequences in the virulent and attenuated groups, respectively. As attenuation was achieved, however, the number of the strain-specific ves sequence decreased. An attractive hypothesis based on the observation of population-based ves gene complexity is that virulence is associated with a unique repertoire of VESA proteins expressed on the infected erythrocyte surface. This is consistent with the mechanism used to generate attenuation. Rapid passage of parasites in cattle without a spleen may result in selection of a non-cytoadherent population that would normally be cleared in spleen intact animals, while cytoadherent parasitized erythrocytes sequester in tissue capillaries. Thus, the small percentage loss of certain ves genes may be responsible for the phenotypic difference between the virulent and attenuated strain pairs. However, an equally plausible hypothesis of attenuation mechanism can be due to the reduction of strain-specific non-ves genes. The reduction of approximately half of all strain-specific non-ves genes post-attenuation is remarkable. Lastly, phenotypic changes exhibited in the infected host may be controlled at the transcriptional or protein level. While non-synonymous polymorphisms were predicted to result in tolerated changes, synonymous single nucleotide polymorphisms differentiating virulent and attenuated parasites in all three strain pairs were also found among the 14 genes. It is plausible that synonymous nucleotide changes may play a role in virulence phenotype, although the mechanism of such is currently unknown.
In summary, this study surveyed and compared genomes of three genetically related B. bovis strain pairs of diverse geographical background and with distinct differences in the ability to cause severe clinical disease (virulence). Common virulence factor(s) at the gene level were not consistently found among all virulent B. bovis strains or their attenuated derivatives. However, the virulent parasite population gene pool was significantly more complex than in attenuated B. bovis parasites. Gene complexity was a result of diversity in both a multi-copy variant gene family and multiple single copy genes. Thus, while there were no common genes consistently associated with attenuation or virulence in these related strain pairs, the distinct and reproducible phenotype of the attenuated derivatives in infected animals may be a result of contraction of virulence-associated genes of varying function during the passage of parasites in the natural host. Virulence mechanisms may also be associated with differences at the transcriptional or protein level, including the ves repertoire, or variations in intergenic or intronic regions . For instance, Iyer et al. recently showed that increased expression of a rhoptry protein in a rodent malaria, P. yoelii (py 235), resulted in the invasion of a wider range of erythrocytes and thus, greater virulence . Alternatively, virulence may be achieved by various means such as epigenetic factors, all of which contribute to the overall phenotype exhibited in infected host [37, 38].
The complete genomic sequences for the five test strains are available at http://www.eecs.wsu.edu/~ananth/WSU_B-bovis_RawData.
We thank Ms. Xiaoya Cheng of Washington State University for performing Southern blot analysis on all genomic DNA samples prior to sequencing and Dr. Kelly Brayton for input on the manuscript. This work is partially funded by the Wellcome Trust GRO075800M (GHP, TFM, AOTL, IE, MBF), USDA-ARS Specific Cooperative Agreement #58-5348-7-528 (TFM), Washington State University, College of Veterinary Medicine intramural grant (AOTL), TCP INTA EEA Rafaela-Asociación Cooperadora 426100 (IE) and National Science Foundation #IIS 0916463 (AK).
- Taylor HM, Kyes SA, Newbold CI: Var gene diversity in Plasmodium falciparum is generated by frequent recombination events. Mol Biochem Parasitol. 2000, 110 (2): 391-397. 10.1016/S0166-6851(00)00286-3.PubMedView ArticleGoogle Scholar
- Rosenberg SM: Mutation for survival. Curr Opin Genet Dev. 1997, 7 (6): 829-834. 10.1016/S0959-437X(97)80047-0.PubMedView ArticleGoogle Scholar
- Taddei F, Matic I, Godelle B, Radman M: To be a mutator, or how pathogenic and commensal bacteria can evolve rapidly. Trends Microbiol. 1997, 5 (11): 427-428. 10.1016/S0966-842X(97)01157-8. discussion 428-429PubMedView ArticleGoogle Scholar
- Taddei F, Radman M, Maynard-Smith J, Toupance B, Gouyon PH, Godelle B: Role of mutator alleles in adaptive evolution. Nature. 1997, 387 (6634): 700-702. 10.1038/42696.PubMedView ArticleGoogle Scholar
- Ping J, Dankar SK, Forbes NE, Keleta L, Zhou Y, Tyler S, Brown EG: PB2 and HA Mutations are Major Determinants of Host Range and Virulence in Mouse-Adapted Influenza A Virus. J Virol.
- Engel AR, Rumyantsev AA, Maximova OA, Speicher JM, Heiss B, Murphy BR, Pletnev AG: The neurovirulence and neuroinvasiveness of chimeric tick-borne encephalitis/dengue virus can be attenuated by introducing defined mutations into the envelope and NS5 protein genes and the 3' non-coding region of the genome. Virology. 405 (1): 243-252.
- Darghouth MA: Review on the experience with live attenuated vaccines against tropical theileriosis in Tunisia: considerations for the present and implications for the future. Vaccine. 2008, 26 (Suppl 6): G4-G10.PubMedView ArticleGoogle Scholar
- Druilhe P, Barnwell JW: Pre-erythrocytic stage malaria vaccines: time for a change in path. Curr Opin Microbiol. 2007, 10 (4): 371-378. 10.1016/j.mib.2007.07.009.PubMedView ArticleGoogle Scholar
- Florin-Christensen M, Schnittger L, Dominguez M, Mesplet M, Rodriguez A, Ferreri L, Asenzo G, Wilkowsky S, Farber M, Echaide I, et al: Search for Babesia bovis vaccine candidates. Parassitologia. 2007, 49 (Suppl 1): 9-12.PubMedGoogle Scholar
- Allred DR: Antigenic variation in Babesia bovis: how similar is it to that in Plasmodium falciparum?. Ann Trop Med Parasitol. 1998, 92 (4): 461-472. 10.1080/00034989859438.PubMedView ArticleGoogle Scholar
- Shkap V, de Vos AJ, Zweygarth E, Jongejan F: Attenuated vaccines for tropical theileriosis, babesiosis and heartwater: the continuing necessity. Trends Parasitol. 2007, 23 (9): 420-426. 10.1016/j.pt.2007.07.003.PubMedView ArticleGoogle Scholar
- Shkap V, Leibovitz B, Krigel Y, Hammerschlag J, Marcovics A, Fish L, Molad T, Savitsky I, Mazuz M: Vaccination of older Bos taurus bulls against bovine babesiosis. Vet Parasitol. 2005, 129 (3-4): 235-242. 10.1016/j.vetpar.2005.01.013.PubMedView ArticleGoogle Scholar
- Bock RE, de Vos AJ: Immunity following use of Australian tick fever vaccine: a review of the evidence. Aust Vet J. 2001, 79 (12): 832-839. 10.1111/j.1751-0813.2001.tb10931.x.PubMedView ArticleGoogle Scholar
- Brayton KA, Lau AO, Herndon DR, Hannick L, Kappmeyer LS, Berens SJ, Bidwell SL, Brown WC, Crabtree J, Fadrosh D, et al: Genome sequence of Babesia bovis and comparative analysis of apicomplexan hemoprotozoa. PLoS Pathog. 2007, 3 (10): 1401-1413.PubMedView ArticleGoogle Scholar
- Bock RE, de Vos AJ, Lew A, Kingston TG, Fraser IR: Studies on failure of T strain live Babesia bovis vaccine. Aust Vet J. 1995, 72 (8): 296-300. 10.1111/j.1751-0813.1995.tb03558.x.PubMedView ArticleGoogle Scholar
- Levy MG, Ristic M: Babesia bovis: continuous cultivation in a microaerophilous stationary phase culture. Science. 1980, 207 (4436): 1218-1220. 10.1126/science.7355284.PubMedView ArticleGoogle Scholar
- Hines SA, McElwain TF, Buening GM, Palmer GH: Molecular characterization of Babesia bovis merozoite surface proteins bearing epitopes immunodominant in protected cattle. Mol Biochem Parasitol. 1989, 37 (1): 1-9. 10.1016/0166-6851(89)90096-0.PubMedView ArticleGoogle Scholar
- Lau AO, Cereceres K, Palmer GH, Fretwell DL, Pedroni MJ, Mosqueda J, McElwain TF: Genotypic diversity of merozoite surface antigen 1 of Babesia bovis within an endemic population. Mol Biochem Parasitol. 2010, 172 (2): 107-112. 10.1016/j.molbiopara.2010.03.017.PubMed CentralPubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Rissman AI, Mau B, Biehl BS, Darling AE, Glasner JD, Perna NT: Reordering contigs of draft genomes using the Mauve aligner. Bioinformatics. 2009, 25 (16): 2071-2073. 10.1093/bioinformatics/btp356.PubMed CentralPubMedView ArticleGoogle Scholar
- Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14 (7): 1394-1403. 10.1101/gr.2289704.PubMed CentralPubMedView ArticleGoogle Scholar
- Delcher AL, Phillippy A, Carlton J, Salzberg SL: Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res. 2002, 30 (11): 2478-2483. 10.1093/nar/30.11.2478.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMed CentralPubMedView ArticleGoogle Scholar
- Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11 (5): 863-874. 10.1101/gr.176601.PubMed CentralPubMedView ArticleGoogle Scholar
- Mangold AJ, Vanzini VR, Echaide IE, de Echaide ST, Volpogni MM, Guglielmone AA: Viability after thawing and dilution of simultaneously cryopreserved vaccinal Babesia bovis and Babesia bigemina strains cultured in vitro. Vet Parasitol. 1996, 61 (3-4): 345-348. 10.1016/0304-4017(95)00839-X.PubMedView ArticleGoogle Scholar
- Somerville RP, Littlebury P, Pipano E, Brown CG, Shkap V, Adamson RE, Oliver RA, Glass EJ, Hall FR: Phenotypic and genotypic alterations associated with the attenuation of a Theileria annulata vaccine cell line from Turkey. Vaccine. 1998, 16 (6): 569-575. 10.1016/S0264-410X(97)84509-4.PubMedView ArticleGoogle Scholar
- Yano K, Nakabayashi T: Attenuation of the virulent RH strain of Toxoplasma gondii by passages in mice immunized with Toxoplasma lysate antigens. Biken J. 1986, 29 (2): 31-37.PubMedGoogle Scholar
- Schmidt LH, Fradkin R, Sesler C, Squires W, Zeyen P: Attenuation of the virulence of the M strain of Plasmodium cynomolgi during prolonged multiplication in splenectomized rhesus monkeys. Am J Trop Med Hyg. 1987, 37 (3): 460-490.PubMedGoogle Scholar
- Van Schaik W, Top J, Riley DR, Boekhorst J, Vrijenhoek JE, Schapendonk CM, Hendrickx AP, Nijman IJ, Bonten MJ, Tettelin H: Pyrosequencing-based comparative genome analysis of the nosocomial pathogen Enterococcus faecium and identification of a large transferable pathogenicity island. BMC Genomics. 11: 239-
- Tauch A, Schneider J, Szczepanowski R, Tilker A, Viehoever P, Gartemann KH, Arnold W, Blom J, Brinkrolf K, Brune I, et al: Ultrafast pyrosequencing of Corynebacterium kroppenstedtii DSM44385 revealed insights into the physiology of a lipophilic corynebacterium that lacks mycolic acids. J Biotechnol. 2008, 136 (1-2): 22-30. 10.1016/j.jbiotec.2008.03.004.PubMedView ArticleGoogle Scholar
- Hall R, Ilhan T, Kirvar E, Wilkie G, Preston PM, Darghouth M, Somerville R, Adamson R: Mechanism(s) of attenuation of Theileria annulata vaccine cell lines. Trop Med Int Health. 1999, 4 (9): A78-84. 10.1046/j.1365-3156.1999.00454.x.PubMedView ArticleGoogle Scholar
- Doolan DL, Dobano C, Baird JK: Acquired immunity to malaria. Clin Microbiol Rev. 2009, 22 (1): 13-36. 10.1128/CMR.00025-08. Table of ContentsPubMed CentralPubMedView ArticleGoogle Scholar
- Langhorne J, Ndungu FM, Sponaas AM, Marsh K: Immunity to malaria: more questions than answers. Nat Immunol. 2008, 9 (7): 725-732. 10.1038/ni.f.205.PubMedView ArticleGoogle Scholar
- Rask TS, Hansen DA, Theander TG, Gorm Pedersen A, Lavstsen T: Plasmodium falciparum erythrocyte membrane protein 1 diversity in seven genomes--divide and conquer. PLoS Comput Biol. 2010, 6 (9):
- Colinet D, Schmitz A, Cazes D, Gatti JL, Poirie M: The origin of intraspecific variation of virulence in an eukaryotic immune suppressive parasite. PLoS Pathog. 6 (11): e1001206-
- Iyer JK, Amaladoss A, Genesan S, Preiser PR: Variable expression of the 235 kDa rhoptry protein of Plasmodium yoelii mediate host cell adaptation and immune evasion. Mol Microbiol. 2007, 65 (2): 333-346. 10.1111/j.1365-2958.2007.05786.x.PubMedView ArticleGoogle Scholar
- May FJ, Li L, Davis CT, Galbraith SE, Barrett AD: Multiple pathways to the attenuation of West Nile virus in South-East Texas in 2003. Virology. 2010, 405 (1): 8-14. 10.1016/j.virol.2010.04.022.PubMed CentralPubMedView ArticleGoogle Scholar
- Broadbent K, Park D, Wolf AR, Van Tyne D, Sims JS, Ribacke U, Volkman S, Duraisingh M, Wirth D, Sabeti PC: A global transcriptional analysis of Plasmodium falciparum malaria reveals a novel family of telomere-associated lncRNAs. Genome Biol. 12 (6): R56-
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.