First report of tomato chlorotic dwarf viroid infecting litchi tomato (Solanum sisymbriifolium).

Litchi tomato (LT) (Solanum sisymbriifolium) is a solanaceous weed that is considered a biological control tool to manage potato cyst nematode (PCN) in Europe and is being explored for use in Idaho. Two Several LT lines were clonally maintained as stocks in the university greenhouse since 2013 and were also established in tissue culture at the same time. In 2018, tomato (Solanum lycopersicum cv. Alisa Craig) scions were grafted onto two LT rootstocks originating either from healthy-looking greenhouse stocks or from tissue culture-maintained plants. Unexpectedly, tomatoes grafted onto the greenhouse-maintained rootstocks of LT displayed severe symptoms of stunting, foliar deformation, and chlorosis, while grafts onto the same LT lines from tissue culture produced healthy-looking tomato plants. Tests for the presence of several viruses known to infect solanaceous plants were conducted on symptomatic tomato scion tissues using ImmunoStrips (Agdia, Elkhard, IN) and RT-PCR (Elwan et al. 2017) but yielded negative results. High throughput sequencing (HTS) was then used to identify possible pathogens that could have been responsible for the symptoms observed in tomato scions. Samples from two symptomatic tomato scions, two asymptomatic scions grafted onto the tissue culture-derived plants, and two greenhouse-maintained rootstocks were subjected to HTS. Total RNA from the four tomato and two LT samples was depleted of ribosomal RNA and subjected to HTS on an Illumina MiSeq platform producing 300-bp paired-end reads and raw reads were adapter and quality cleaned. For the tomato samples, the clean reads were mapped against the S. lycopersicum L. reference genome, and unmapped paired reads were assembled producing between 4,368 and 8,645 contigs. For the LT samples, all clean reads were directly assembled, producing 13,982 and 18,595 contigs. In the symptomatic tomato scions and the two LT rootstock samples, a 487-nt contig was found, comprising an ~1.35 tomato chlorotic dwarf viroid (TCDVd) genome and exhibiting 99.7% identity with it (GenBank accession AF162131; Singh et al. 1999). No other virus-related or viroid contigs were identified. RT-PCR analysis using a pospiviroid primer set Pospi1-FW/RE (Verhoeven et al. 2004), and a TCDVd-specific primer set TCDVd-Fw/TCDVd-Rev (Olmedo-Velarde et al. 2019) produced 198-nt and 218-nt bands, respectively, thus confirming the presence of TCDVd in tomato and LT samples. These PCR products were Sanger sequenced and confirmed to be TCDVd-specific; the complete sequence of the Idaho isolate of TCDVd was deposited in GenBank under the accession number OQ679776. Presence of TCDVd in LT plant tissue was confirmed by the APHIS PPQ Laboratory in Laurel, MD. Asymptomatic tomatoes and LT plants from tissue culture were found negative for TCDVd. Previously, TCDVd was reported to affect greenhouse tomatoes in Arizona and Hawaii (Ling et al. et al. 2009; Olmedo-Velarde et al. 2019), however, this is the first report of TCDVd infecting litchi tomato (S. sisymbriifolium). Five additional greenhouse-maintained LT lines were found TCDVd-positive using RT-PCR and Sanger sequencing. Given the very mild or asymptomatic infection of TCDVd in this host, molecular diagnostic methods should be used to screen LT lines for the presence of this viroid to avoid inadvertent spread of TCDVd. Another viroid, potato spindle tuber viroid, was reported to be transmitted through LT seed (Fowkes et al. 2021), and transmission of TCDVd through LT seed may also be responsible for this TCDVd outbreak in the university greenhouse, although no direct evidence was collected. To the best of our knowledge, this is the first report of TCDVd infection in S. sisymbriifolium and also the first report of the TCDVd occurrence in Idaho.

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.

Target specific-trisubstituted pyrrole inhibits Babesia bovis erythrocytic growth.

Babesiosis, a significant veterinary disease and an emerging zoonotic human infection, is caused by certain species of the protozoan parasite, Babesia. Here we report that a trisubstituted pyrrole is a potent inhibitor of Babesia bovis, a bovine parasite. Furthermore, B. bovis expresses the known target of the compound, the cGMP dependent protein kinase. Target conservation and the in vitro efficacy support further investigation of this compound and validation of Babesia cGMP dependent protein kinase as its in vivo target.

Babesia bovis: a bipartite signal directs the glutamyl-tRNA synthetase to the apicoplast.

Babesia bovis contains a prokaryotic derived organelle known as the apicoplast. Many participants of the metabolic pathways within the apicoplast are encoded in the nuclear genome and post-translationally imported with the help of a bipartite signal. Recently, an all encompassing algorithm was derived to predict apicoplast targeted proteins for many non-Plasmodium apicomplexans in which it reported the presence of 260 apicoplast targeted proteins in Babesia. One of these proteins is glutamyl tRNA synthetase (GltX). This study investigates if the putative bipartite signal of GltX alone is sufficient to direct proteins into the apicoplast. Using a transient transfection system consisting of a green fluorescent protein as the reporter, we tested the signal and transit portions of the bipartite signal in apicoplastic transport. We first identified the transcript of gltX to be expressed during the asexual blood stages and subsequently confirmed that the complete bipartite signal is responsible for directing the reporter protein into a compartment distinct from the nucleus and the mitochondrion. As GltX bipartite signal successfully guided the reporter protein into the apicoplast, our finding implies that it also directs native GltX into the same organelle.

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.

Structure and evolution of a proviral locus of Glyptapanteles indiensis bracovirus.

BACKGROUND: Bracoviruses (BVs), a group of double-stranded DNA viruses with segmented genomes, are mutualistic endosymbionts of parasitoid wasps. Virus particles are replication deficient and are produced only by female wasps from proviral sequences integrated into the wasp genome. Virus particles are injected along with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and therefore perpetuation of proviral DNA. Here we describe a 223 kbp region of Glyptapanteles indiensis genomic DNA which contains a part of the G. indiensis bracovirus (GiBV) proviral genome. RESULTS: Eighteen of ~24 GiBV viral segment sequences are encoded by 7 non-overlapping sets of BAC clones, revealing that some proviral segment sequences are separated by long stretches of intervening DNA. Two overlapping BACs, which contain a locus of 8 tandemly arrayed proviral segments flanked on either side by ~35 kbp of non-packaged DNA, were sequenced and annotated. Structural and compositional analyses of this cluster revealed it exhibits a G+C and nucleotide composition distinct from the flanking DNA. By analyzing sequence polymorphisms in the 8 GiBV viral segment sequences, we found evidence for widespread selection acting on both protein-coding and non-coding DNA. Comparative analysis of viral and proviral segment sequences revealed a sequence motif involved in the excision of proviral genome segments which is highly conserved in two other bracoviruses. CONCLUSION: Contrary to current concepts of bracovirus proviral genome organization our results demonstrate that some but not all GiBV proviral segment sequences exist in a tandem array. Unexpectedly, non-coding DNA in the 8 proviral genome segments which typically occupies ~70% of BV viral genomes is under selection pressure suggesting it serves some function(s). We hypothesize that selection acting on GiBV proviral sequences maintains the genetic island-like nature of the cluster of proviral genome segments described herein. In contrast to large differences in the predicted gene composition of BV genomes, sequences that appear to mediate processes of viral segment formation, such as proviral segment excision and circularization, appear to be highly conserved, supporting the hypothesis of a single origin for BVs.

Bumped kinase inhibitor prohibits egression in Babesia bovis.

Babesiosis is a global zoonotic disease acquired by the bite of a Babesia-infected Ixodes tick or through blood transfusion with clinical relevance affecting humans and animals. In this study, we evaluated a series of small molecule compounds that have previously been shown to target specific apicomplexan enzymes in Plasmodium, Toxoplasma and Cryptosporidium. The compounds, bumped kinase inhibitors (BKIs), have strong therapeutic potential targeting apicomplexa-specific calcium dependent protein kinases (CDPKs). We investigated if BKIs also show inhibitory activities against piroplasms such as Babesia. Using a subset of BKIs that have promising inhibitory activities to Plasmodium and Toxoplasma, we determined that their actions ranged from 100% and no inhibition against Babesia bovis blood stages. One specific BKI, RM-1-152, showed complete inhibition against B. bovis within 48h and was the only BKI that showed noticeable phenotypic changes to the parasites. Focusing our study on this BKI, we further demonstrated that RM-1-152 has Babesia-static activity and involves the prohibition of merozoite egress while replication and re-invasion of host cells are unaffected. The distinct, abnormal phenotype induced by RM-1-152 suggests that this BKI can be used to investigate less studied cellular processes such as egression in piroplasm.

Characterization of acyl carrier protein and LytB in Babesia bovis apicoplast.

The apicoplast is a highly specialized organelle that mediates required functions in the growth and replication of apicomplexan parasites. Despite structural conservation of the apicoplast among different parasite genera and species, there are also critical differences in the metabolic requirements of different parasites and at different stages of the life cycle. To specifically compare apicoplast pathways between parasites that have both common and unique stages, we characterized the apicoplast in Babesia bovis, which has only intraerythrocytic asexual stages in the mammalian host, and compared it to that of Plasmodium falciparum, which has both asexual intraerythrocytic and hepatic stages. Specifically focusing on the type II fatty acid (FASII) and isoprenoid (MEP) biosynthesis pathways, we searched for pathway components and retention of active sites within the genome, localized key components [acyl carrier protein (ACP) and 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (LytB)] to the apicoplast, and demonstrated that the N-terminal bipartite signals of both proteins are required and sufficient for trafficking to the apicoplast lumen. Using specific pharmacologic inhibition, we demonstrated that MEP biosynthesis may be disrupted and its presence is required for intraerythrocytic growth of B. bovis asexual stages, consistent with the genomic pathway analysis and with its requirement in the asexual erythrocytic stages of P. falciparum. In contrast, FASII biosynthesis may or may not be present and specific drug targets did not have any inhibitory effect to B. bovis intraerythrocytic growth, which is consistent with the lack of requirement for P. falciparum intraerythrocytic growth. However, genomic analysis revealed the loss of FASII pathway components in B. bovis whereas the pathway is intact for P. falciparum but regulated to be expressed when needed (hepatic stages) and silent when not (intraerythrocytic stages). The results indicate specialized molding of apicoplast biosynthetic pathways to meet the requirements of individual apicomplexan parasites and their unique intracellular niches.

Genome-wide analysis of peptidase content and expression in a virulent and attenuated Babesia bovis strain pair.

Identifying virulence determinants in Apicomplexan parasites remains a major gap in knowledge for members within this phylum. We hypothesized that peptidases would segregate with virulence between a virulent parent Babesia bovis strain and an attenuated daughter strain derived by rapid in vivo passage. Using the complete genome sequence of the virulent T2Bo strain, 66 peptidases were identified and active sites confirmed. The presence, sequence identity and expression levels were tested for each of the 66 peptidases in the virulent parent and attenuated daughter T2Bo strains using whole genome, targeted sequencing approaches and microarrays analyses. Quantitative PCR revealed that there was no significant difference in peptidase expression between the virulent and attenuated strains. We conclude that while peptidases may well play a required role in B. bovis pathogenesis, neither loss of peptidase gene content nor reduced gene expression underlies the loss of virulence associated with in vivo passage and attenuation.

Genotypic diversity of merozoite surface antigen 1 of Babesia bovis within an endemic population.

Multiple genetically distinct strains of a pathogen circulate and compete for dominance within populations of animal reservoir hosts. Understanding the basis for genotypic strain structure is critical for predicting how pathogens respond to selective pressures and how shifts in pathogen population structure can lead to disease outbreaks. Evidence from related Apicomplexans such as Plasmodium, Toxoplasma, Cryptosporidium and Theileria suggests that various patterns of population dynamics exist, including but not limited to clonal, oligoclonal, panmictic and epidemic genotypic strain structures. In Babesia bovis, genetic diversity of variable merozoite surface antigen (VMSA) genes has been associated with disease outbreaks, including in previously vaccinated animals. However, the extent of VMSA diversity within a defined population in an endemic area has not been examined. We analyzed genotypic diversity and temporal change of MSA-1, a member of the VMSA family, in individual infected animals within a reservoir host population. Twenty-eight distinct MSA-1 genotypes were identified within the herd. All genotypically distinct MSA-1 sequences clustered into three groups based on sequence similarity. Two thirds of the animals tested changed their dominant MSA-1 genotypes during a 6-month period. Five animals within the population contained multiple genotypes. Interestingly, the predominant genotypes within those five animals also changed over the 6-month sampling period, suggesting ongoing transmission or emergence of variant MSA-1 genotypes within the herd. This study demonstrated an unexpected level of diversity for a single copy gene in a haploid genome, and illustrates the dynamic genotype structure of B. bovis within an individual animal in an endemic region. Co-infection with multiple diverse MSA-1 genotypes provides a basis for more extensive genotypic shifts that characterizes outbreak strains.

Comparative genomics of mutualistic viruses of Glyptapanteles parasitic wasps.

BACKGROUND: Polydnaviruses, double-stranded DNA viruses with segmented genomes, have evolved as obligate endosymbionts of parasitoid wasps. Virus particles are replication deficient and produced by female wasps from proviral sequences integrated into the wasp genome. These particles are co-injected with eggs into caterpillar hosts, where viral gene expression facilitates parasitoid survival and, thereby, survival of proviral DNA. Here we characterize and compare the encapsidated viral genome sequences of bracoviruses in the family Polydnaviridae associated with Glyptapanteles gypsy moth parasitoids, along with near complete proviral sequences from which both viral genomes are derived. RESULTS: The encapsidated Glyptapanteles indiensis and Glyptapanteles flavicoxis bracoviral genomes, each composed of 29 different size segments, total approximately 517 and 594 kbp, respectively. They are generated from a minimum of seven distinct loci in the wasp genome. Annotation of these sequences revealed numerous novel features for polydnaviruses, including insect-like sugar transporter genes and transposable elements. Evolutionary analyses suggest that positive selection is widespread among bracoviral genes. CONCLUSIONS: The structure and organization of G. indiensis and G. flavicoxis bracovirus proviral segments as multiple loci containing one to many viral segments, flanked and separated by wasp gene-encoding DNA, is confirmed. Rapid evolution of bracovirus genes supports the hypothesis of bracovirus genes in an 'arms race' between bracovirus and caterpillar. Phylogenetic analyses of the bracoviral genes encoding sugar transporters provides the first robust evidence of a wasp origin for some polydnavirus genes. We hypothesize transposable elements, such as those described here, could facilitate transfer of genes between proviral segments and host DNA.