Comparative transcriptome analysis of geographically distinct virulent and attenuated Babesia bovis strains reveals similar gene expression changes through attenuation

Background: Loss of virulence is a phenotypic adaptation commonly seen in prokaryotic and eukaryotic pathogens. This mechanism is not well studied, especially in organisms with multiple host and life cycle stages such as Babesia, a tick-transmitted hemoparasite of humans and animals. B. bovis, which infects cattle, has naturally occurring virulent strains that can be reliably attenuated in vivo. Previous studies suggest the virulence loss mechanism may involve post-genomic modification. We investigated the transcriptome profiles of two geographically distinct B. bovis virulent and attenuated strain pairs to better understand virulence loss and to gain insight into pathogen adaptation strategies. Results: Expression microarray and RNA-sequencing approaches were employed to compare transcriptome profiles of two B. bovis strain pairs, with each pair consisting of a virulent parental and its attenuated derivative strain. Differentially regulated transcripts were identified within each strain pair. These included genes encoding for VESA1, SmORFs, undefined membrane and hypothetical proteins. The majority of individual specific gene transcripts differentially regulated within a strain were not shared between the two strains. There was a disproportionately greater number of ves genes upregulated in the virulent parental strains. When compared with their attenuated derivatives, divergently oriented ves genes were included among the upregulated ves genes in the virulent strains, while none of the upregulated ves genes in the attenuated derivatives were oriented head to head. One gene family whose specific members were consistently and significantly upregulated in expression in both attenuated strains was spherical body protein (SBP) 2 encoding gene where SBP2 truncated copies 7, 9 and 11 transcripts were all upregulated. Conclusions: We conclude that ves heterodimer pair upregulation and overall higher frequency of ves gene expressions in the virulent strains is consistent with the involvement of this gene family in virulence. This is logical given the role of VESA1 proteins in cytoadherence of infected cells to endothelial cells. However, upregulation of some ves genes in the attenuated derivatives suggests that the consequence of upregulation is gene-specific. Furthermore, upregulation of the spherical body protein 2 gene family may play a role in the attenuated phenotype. Exactly how these two gene families may contribute to the loss or gain of virulence is discussed.EEA RafaelaFil: Pedroni, Monica J. Washington State University. College of Veterinary Medicine. Department of Veterinary Microbiology & Pathology. Program of Genomics; Estados UnidosFil: Sondgeroth, Kerry S. Washington State University. College of Veterinary Medicine. Department of Veterinary Microbiology & Pathology. Program of Genomics; Estados UnidosFil: Gallego-Lopez, Gina M. Washington State University. College of Veterinary Medicine. Department of Veterinary Microbiology & Pathology. Program of Genomics; Estados UnidosFil: Echaide, Ignacio Eduardo. Instituto Nacional de Tecnología Agropecuaria (INTA). Estación Experimental Agropecuaria Rafaela. Laboratorio de Inmunología y Parasitología; ArgentinaFil: Lau, Audrey OT. Washington State University. College of Veterinary Medicine. Department of Veterinary Microbiology & Pathology. Program of Genomics; Estados Unidos. Washington State University. College of Veterinary Medicine. Paul G. Allen School for Global Animal Health; Estados Unido

Theileria equi isolates vary in susceptibility to imidocarb dipropionate but demonstrate uniform in vitro susceptibility to a bumped kinase inhibitor

Abstract Background The apicomplexan hemoparasite Theileria equi is a causative agent of equine piroplasmosis, eradicated from the United States in 1988. However, recent outbreaks have sparked renewed interest in treatment options for infected horses. Imidocarb dipropionate is the current drug of choice, however variation in clinical response to therapy has been observed. Methods We quantified the in vitro susceptibility of two T. equi isolates and a lab generated variant to both imidocarb dipropionate and a bumped kinase inhibitor compound 1294. We also evaluated the capacity of in vitro imidocarb dipropionate exposure to decrease susceptibility to that drug. The efficacy of imidocarb dipropionate for clearing infection in four T. equi infected ponies was also assessed. Results We observed an almost four-fold difference in imidocarb dipropionate susceptibility between two distinct isolates of T. equi. Four ponies infected with the less susceptible USDA Florida strain failed to clear the parasite despite two rounds of treatment. Importantly, a further 15-fold decrease in susceptibility was produced in this strain by continuous in vitro imidocarb dipropionate exposure. Despite a demonstrated difference in imidocarb dipropionate susceptibility, there was no difference in the susceptibility of two T. equi isolates to bumped kinase inhibitor 1294. Conclusions The observed variation in imidocarb dipropionate susceptibility, further reduction in susceptibility caused by drug exposure in vitro, and failure to clear T. equi infection in vivo, raises concern for the emergence of drug resistance in clinical cases undergoing treatment. Bumped kinase inhibitors may be effective as alternative drugs for the treatment of resistant T. equi parasites

Attenuation of virulence in an apicomplexan hemoparasite results in reduced genome diversity at the population level

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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

Comparative genomic analysis and phylogenetic position of Theileria equi

BACKGROUND: Transmission of arthropod-borne apicomplexan parasites that cause disease and result in death or persistent infection represents a major challenge to global human and animal health. First described in 1901 as Piroplasma equi, this re-emergent apicomplexan parasite was renamed Babesia equi and subsequently Theileria equi, reflecting an uncertain taxonomy. Understanding mechanisms by which apicomplexan parasites evade immune or chemotherapeutic elimination is required for development of effective vaccines or chemotherapeutics. The continued risk of transmission of T. equi from clinically silent, persistently infected equids impedes the goal of returning the U. S. to non-endemic status. Therefore comparative genomic analysis of T. equi was undertaken to: 1) identify genes contributing to immune evasion and persistence in equid hosts, 2) identify genes involved in PBMC infection biology and 3) define the phylogenetic position of T. equi relative to sequenced apicomplexan parasites. RESULTS: The known immunodominant proteins, EMA1, 2 and 3 were discovered to belong to a ten member gene family with a mean amino acid identity, in pairwise comparisons, of 39%. Importantly, the amino acid diversity of EMAs is distributed throughout the length of the proteins. Eight of the EMA genes were simultaneously transcribed. As the agents that cause bovine theileriosis infect and transform host cell PBMCs, we confirmed that T. equi infects equine PBMCs, however, there is no evidence of host cell transformation. Indeed, a number of genes identified as potential manipulators of the host cell phenotype are absent from the T. equi genome. Comparative genomic analysis of T. equi revealed the phylogenetic positioning relative to seven apicomplexan parasites using deduced amino acid sequences from 150 genes placed it as a sister taxon to Theileria spp. CONCLUSIONS: The EMA family does not fit the paradigm for classical antigenic variation, and we propose a novel model describing the role of the EMA family in persistence. T. equi has lost the putative genes for host cell transformation, or the genes were acquired by T. parva and T. annulata after divergence from T. equi. Our analysis identified 50 genes that will be useful for definitive phylogenetic classification of T. equi and closely related organisms