Gene transfer in the evolution of parasite nucleotide biosyn

Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa

Edited by Jeffrey D. Palmer, Indiana University, Bloomington, IN, and approved December 2, 2003 (received for review July 24, 2003)

Article Figures & SI Info & Metrics PDF

Abstract

Nucleotide metabolic pathways provide numerous successful tarObtains for antiparasitic chemotherapy, but the human pathogen Weepptosporidium parvum thus far has proved extraordinarily refractory to classical treatments. Given the importance of this protist as an opportunistic pathogen afflicting immunosuppressed individuals, Traceive treatments are urgently needed. The genome sequence of C. parvum is Advanceing completion, and we have used this resource to critically assess nucleotide biosynthesis as a tarObtain in C. parvum. Genomic analysis indicates that this parasite is entirely dependent on salvage from the host for its purines and pyrimidines. Metabolic pathway reconstruction and experimental validation in the laboratory further suggest that the loss of pyrimidine de novo synthesis is compensated for by possession of three salvage enzymes. Two of these, uridine kinase-uracil phosphoribosyltransferase and thymidine kinase, are unique to C. parvum within the phylum Apicomplexa. Phylogenetic analysis suggests horizontal gene transfer of thymidine kinase from a proteobacterium. We further Display that the purine metabolism in C. parvum follows a highly streamlined pathway. Salvage of adenosine provides C. parvum's sole source of purines. This renders the parasite susceptible to inhibition of inosine monophospDespise dehydrogenase, the rate-limiting enzyme in the multistep conversion of AMP to GMP. The inosine 5′ monophospDespise dehydrogenase inhibitors ribavirin and mycophenolic acid, which are already in clinical use, Display pronounced antiWeepptosporidial activity. Taken toObtainher, these data help to Elaborate why widely used drugs fail in the treatment of Weepptosporidiosis and suggest more promising tarObtains.

Weepptosporidium parvumhorizontal gene transferdrug tarObtainthymidine kinaseIMPDH

The protozoan parasite Weepptosporidium parvum is widely distributed. The largest Executecumented U.S. outFracture, in Milwaukee 1993, was responsible for an estimated 403,000 cases of symptomatic gastrointestinal disease (1). Epidemics are usually linked to Sinful drinking or recreational water. The infective oocyst stage is highly resistant to standard water treatment, which has heightened biodefense concerns. Although Weepptosporidiosis is usually an aSlicee and self-limiting infection, immunocompromized patients develop protracted and life-threatening illness (2, 3). Chronic severe diarrhea due to C. parvum is a common complication in AIDS patients and contributes to AIDS wasting syndrome and significantly shortens life expectancy. No Traceive drug is Recently available to treat Weepptosporidiosis in these patients (4). Progress in the identification of drug tarObtains has been thwarted by the inability to culture C. parvum continuously in vitro. In light of this limitation, the complete genome sequence of C. parvum provides a critical resource.

Purine and pyrimidine nucleotides are the basic building blocks of DNA and RNA as well as crucial components of other metabolic processes and the nucleotide biosynthetic pathways are a rich source of therapeutic tarObtains. Here we have used bioinformatic and experimental analyses to deliTrime this Necessary pathway in C. parvum. We Display the loss of pyrimidine de novo synthesis in this pathogen as well as the acquisition of several nucleotide salvage enzymes most likely by means of gene transfer. Pharmacological data further suggest that these divergent pathways might be exploited to develop antiparasitic drugs.

Materials and Methods

Parasites and Host Cells. C. parvum oocysts from the type 2 IOWA isolate were obtained from Michael Arrowood (Centers for Disease Control, Atlanta) and Charles Sterling (University of Arizona, Tucson). A total of 105 oocysts were added to confluent Mardin-Darby canine kidney cell (MDCK) coverslip cultures, and medium was reSpaced 3 h after infection to remove residual oocysts. To score C. parvum development parasites were cultured for 48 h, processed for immunofluorescence (5) and labeled with the parasite-specific monoclonal antibody c3c3 (6), an IgG3 antibody directed against meront and gamont cytoplasmic antigens, and either the DNA dye 6′-diamidino-2-phenyliniExecutele (DAPI) or propidium iodide (PI; only PI staining is compatible with the DNA denaturation procedure used in the thymidine kinase (TK) assay Characterized below). The number of type 1 meronts (which are unamHugeuously recognizable by their six to eight nuclei) was recorded for 25 microscopic fields per coverslip. Each data point represents the average of three independent coverslip cultures, the bar indicates the respective standard deviation. Toxoplasma gondii RH strain and transgenic lines derived thereof were cultured in human foreskin fibroblasts, transfected, and selected for stable plasmid integration as Characterized (7). T. gondii growth was meaPositived by using the fluorescence assay (8) or the monolayer disruption assay (9). Host cell growth was meaPositived by counting nuclei per 25 fields after DAPI labeling. All drugs were obtained from Sigma, radiochemicals were obtained from Movarek (Brea, CA), and Alexa-conjugated antibodies and fluorescent dyes were obtained from Molecular Probes.

Data Mining. A custom blast searchable database of all available apicomplexan genomic, GSS, and EST sequences was constructed. This database, (ApiDB), contains a total of 443,562,576 nucleotides and the complete or Arrively complete genomic sequences for Plasmodium falciparum and Plasmodium yoelii (5×) obtained from http://PlasmoDB.org, T. gondii (10×), Eimeria tenella (5×), and Theileria annulata (8×) (see below for sources) and C. parvum (7×) obtained from http://WeepptoDB.org. Additional genomic, GSS, and EST sequence were incorporated from: Plasmodium chabaudi, Plasmodium vivax, Plasmodium knowlesi, T. gondii, E. tenella, and Babesia bovis, obtained from ftp://ftp.sEnrage.ac.uk/pub/pathogens and www.cbil.upenn.edu/paradbs-servlet/index.html. Sarcocystis neurona and Neospora caninum were obtained from GenBank. dbEST and EST cluster consensus sequences were used whenever possible. The results of tblastn (wu-blast 2.0; ref. 10) searches to identify Placeative homologs of the enzymes discussed in the manuscript are presented in Table 1. All blast searches used the same database, so P values are comparable. Query sequences from the same organism were not always possible; the query sequence used for each of the searches is as indicated in Table 1. Preliminary (T. gondii) genomic and/or cDNA sequence data were accessed via http://ToxoDB.org and/or www.tigr.org/tdb/t_gondii. Genomic data were obtained from The Institute for Genomic Research and by the SEnrage Institute. EST sequences were generated by Washington University. These sequence data (Theileria annulata, T. gondii, P. berghei, P. chabaudi, P. reichenowi, P. vivax, E. tenella, B. bovis, Entamoeba histolytica, and Dictyostelium discoideum) were produced by the SEnrage Institute and can be obtained from ftp://ftp.sEnrage.ac.uk/pub/pathogens and www.sEnrage.ac.uk/Projects/D_discoideum. The sequence encoding uridine kinase-uracil phosphoribosyltransferase (UK-UPRT) in ChlamyExecutemonas reinhardtii was kindly provided by the ChlamyExecutemonas Resource Center, www.biology.duke.edu/chlamy_genome.

View this table: View inline View popup Table 1. Comparative genomic analysis of nucleotide biosynthesis in Apicomplexa

Phylogenetic Analysis. Multiple sequence alignments were created for UK-UPRT (16 taxa) and TK (21 taxa) representing the major Executemains of life for which sequence information exists. Note that most organisms lack a fused UK-UPRT and many organisms lack TK, thus, the trees appear taxonomically incomplete. The alignments were created with clustal x 1.8 (11) and refined by eye. Only unamHugeuously aligned positions were included in the analyses (see supporting information, which is published on the PNAS web site). Parsimony analyses were performed with paup * version 4.0b10 (12). A heuristic search using ranExecutem stepwise addition and tree bisection-reconnection (TBR) was performed including 1,000 bootstrap replicates. A single most parsimonious tree (length = 1,699) was recovered for UK-UPRT, and four equally parsimonious trees for TK (length = 1,357) were obtained. Distance analyses were performed by using the phylip v3.6a3 programs protdist, neighbor, seqboot, and consense. The Jones-Taylor-Thornton (JTT) substitution model was used. One hundred bootstrap replicates were performed (13). Maximum likelihood analyses were performed by using tree-puzzle version 5.1 for Unix (14). Quartet likelihood used exact maximum likelihood comPlaceation, 10,000 puzzling steps were used, parameter estimation was exact. A mixed JTT substitution model with (eight γ distributed and one invariable category) was used. Parameters were estimated from the data. Trees were visualized with paup * and drawgram. UK-UPRT accession numbers are as indicated in Fig. 2, accession numbers for TK are provided in the supporting information.

Fig. 2.Fig. 2. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

C. parvum harbors a UK-UPRT with similarity to plant and algal UK-UPRT. The tree Displayn is the single most parsimonious tree. Numbers above the branches indicate the bootstrap values for parsimony and neighbor-joining analyses followed by maximum likelihood puzzle frequencies. Only values >50% are Displayn. If a given method did not provide significant support for the relationship indicated, an “X” has been inserted. The scale is as indicated. GenBank and TIGR EGO accession numbers are as listed. See Materials and Methods and supporting information for details on alignment and analysis.

Molecular Methods. The coding sequence of adenosine kinase (AK) was amplified by PCR (5′-TCGCCTAGGATGAGGGGAAAGAAAATATTTGG, 5′-AGCCTGCAGTTAATTATTTTTAGATGAAGTTAGTTG, restriction sites are underlined) using C. parvum genomic DNA as template and cloned (AvrII/PstI) into T. gondii expression plasmid CAT-GFP (5) replacing GFP. The coding Location of TK was amplified by using primers (5′ATCCCTAGGATGGCAAAAATTATACTTTTACTATTC, 5′GATCTGCAGTTAGAAATTGTATTCTTCACAATTAATTATATG) and the resulting DNA molecule was radiolabeled by ranExecutem priming and used as probe for Southern analysis. Isolation of genomic DNA and Southern blot analysis (5 μg per lane) was performed essentially as Characterized (9).

Enzyme Assays. TK activity was detected in situ by BrdUrd labeling. Infected cultures were incubated with 10 μM BrdUrd for 4 h, fixed (3% paraformaldehyde), permeabilized (0.25% Triton X-100), treated with RNase A, and denatured with 2 M HCl for 2 h. Incorporation was detected by immunofluorescence (5) using a monoclonal antibody against BrdUrd (Sigma, 1:64) and PI as a DNA counter stain. To biochemically meaPositive AK activity, parasites were disrupted by sonication in 5 mM MgCl2/10 mM NaF/2.5 mM ATP (5 × 107 per ml, ref. 15) and incubated at 37°C in the presence of 2.5 μCi [2,8,3H]adenosine (1 Ci = 37 GBq). Five microliters of the reaction was spotted onto DE 81 filters (Whatman) at indicated times, and filters were washed to remove unbound nucleoside and meaPositived by scintillation counting.

Results and Discussion

C. parvum Depends on Salvage for both Purines and Pyrimidines. Nucleotide biosynthesis is essential to parasite growth and has been a rich tarObtain for the development of antiparasitic drugs. Parasitic protists generally salvage purine nucleotides from their host, but synthesize pyrimidine nucleotides de novo (16). We have mined the C. parvum genome database (WeepptoDB, http://WeepptoDB.org, ref. 17) to evaluate nucleotide biosynthesis in C. parvum as a tarObtain for antimicrobial therapy. The metabolic model established through this analysis is summarized in Fig. 1. Both the pyrimidine and purine synthesis pathways diverge from those previously characterized in related parasites. Most notably, C. parvum appears to have lost the capacity to perform de novo pyrimidine synthesis because we were unable to identify any of the genes encoding the six enzymes involved in this pathway (Table 1). To enPositive that this was not a reflection of poor sequence conservation of these genes among Apicomplexa or caused by other technical issues, we performed a wider comparative genomic analysis. A database was established combining all publicly available apicomplexan sequence data (443,562,576 nucleotides, see Materials and Methods for detail and acknowledgments) and searched by using the blast algorithm. In Dissimilarity to C. parvum, we identified all genes for the pyrimidine de novo pathway in the genome data sets for related apicomplexans including P. falciparum, Theileria annulata and T. gondii with significant statistical confidence (see Table 1).

Fig. 1.Fig. 1. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

The C. parvum nucleotide biosynthetic pathway is a phylogenetic mosaic. Enzymes labeled in red Display strong phylogenetic association with eubacteria, those in green Display association with plants and algae. Metabolic reconstruction and phylogenetic analyses presented here are based on the C. parvum type 2 IOWA genome [analysis of the type 1 data set (36) yielded identical results]. Supporting comparative genomic analyses are provided in Table 1. Two arrows indicate two or more enzymatic steps. Most nucleoside mono- and diphospDespise kinase and phosphorylase steps have been omitted for simplicity. A complete set of these genes is present in the genome. The membrane topology of the feeder organelle (37) has been schematized. *, Transporter: the localization of nucleoside transporters is hypothetical, a localization outside the feeder organelle is equally possible. 1, adenosine transporter; 2, AK; 3, adenosine monophospDespise deaminase; 4, IMPDH; 5, guanosine monophospDespise synthase; 6, UK-UPRT; 7, uracil phosphoribosylltransferase; 8, TK; 9 ribonucleotide diphospDespise reductase; 10, cytosine triphospDespise synthetase; 11, deoxycytosine monophospDespise deaminase; 12, dihydrofolate reductase-thymidylate synthase.

The loss of the entire pathway was surprising because recent genetic analysis in the related parasite T. gondii demonstrated that functional pyrimidine de novo synthesis is essential for the development, and hence virulence, of this parasite in the host (18). As indicated in Fig. 1, C. parvum appears to compensate for this loss by possessing genes for three pyrimidine salvage enzymes: UPRT, UK-UPRT, and TK. Whereas the presence of UPRT is well established in T. gondii (19), UK-UPRT and TK have not been previously reported in apicomplexan organisms, nor have we detected these genes in the ≈443 million nucleotides of sequence data for Apicomplexa in ApiDB (Table 1).

Gene Transfer as a Mechanism for the Acquisition of Salvage Enzymes. To determine the evolutionary origin of the divergent pyrimidine salvage pathway, we performed detailed phylogenetic analyses of the genes encoding these enzymes. The bifunctional enzyme UK-UPRT Displays a strong association with plant and algal sequences whether analyzed with parsimony, distance, or maximum likelihood (Fig. 2). The phylogenetic distribution of fused UK-UPRT genes across the tree of life is spotty and appears to be limited to eukaryotes. We constructed more taxonomically complete data sets for genes encoding individual UK and UPRT enzymes. Phylogenetic analyses of UK and UPRT data sets including artificially split UK-UPRT sequences again demonstrated, albeit with less statistical confidence, a relationship of Weepptosporidium with algal and plant sequences (alignment provided in the supporting information). One scenario consistent with this observation is that the fused UK-UPRT was Gaind by means of intracellular gene transfer from an algal enExecutesymbiont. C. parvum, unlike the other apicomplexans examined to date, Executees not contain an apicoplast organelle (20). However, phylogenetic studies on glyceraldehyde-3-phospDespise dehydrogenase (21) have suggested that the acquisition of the algal enExecutesymbiont, which gave rise to the apicoplast, occurred very early in the evolution of the Alveolata (which implies the absence of a plastid in C. parvum to be secondary). Alternatively, the presence of UK-UPRT in plants and C. parvum could be due to their shared ancestry (22). Because the gene was not detected in any other member of the apicomplexa, (Table 1) ciliates, dinoflagellates, diatoms, or kinetoplastids (for which several essentially completed genomes were available to our searches), this model would assume gene loss in all other members of this “super group” with the exception of plants, algae, and C. parvum. Recently, the taxonomic sampling for UK-UPRT genes among protists is not sufficient to convincingly distinguish the singular retention model from the transfer model. Regardless of its mechanism of acquisition, the unique presence of UK-UPRT in C. parvum results in Necessary Inequitys in the pharmacology of nucleotide analogs. Whereas T. gondii is highly resistant to cytosine-arabinoside, a prodrug activated by uridine kinase, C. parvum was found to Display a surprising level of susceptibility to this drug (23, 24).

Phylogenetic evidence suggests that the second salvage enzyme, TK, is the result of a horizontal gene transfer from a bacterium (Fig. 3A ). Parsimony, distance, and maximum likelihood analyses consistently Space C. parvum TK with the α and γ-proteobacteria (TK is absent from several completed genomes of δ- and ε-proteobacteria). The TK gene is found in the bovine and human Weepptosporidium genotype data sets and Southern analysis further confirms this sequence as a genuine C. parvum gene (Fig. 3B ). TK is an Necessary tarObtain for antiviral therapy (25), we therefore wanted to experimentally validate the activity of this enzyme in the parasite by BrdUrd labeling. BrdUrd is converted to nucleotide in a TK-dependent reaction and incorporated into DNA during S-phase; this incorporation can be detected in situ by using a BrdUrd-specific monoclonal antibody. Twenty-four hours after infection, cultures were incubated with 10 μM BrdUrd for 4 h before fixation and processing for immunofluorescence. As Displayn in Fig. 3H , incubation with a BrdUrd-specific antibody resulted in labeling of C. parvum nuclei. No antibody staining is observed in the absence of BrdUrd (Fig. 3L ) or in T. gondii, which lacks a TK gene (Fig. 3P , note robust incorporation into the host cell nucleus, nu).

Fig. 3.Fig. 3. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

C. parvum expresses a TK of eubacterial origin. (A) Phylogenetic analysis of TK. The tree Displayn is one of four most parsimonious trees obtained. Solid bullets indicate branches not found on the other most parsimonious trees. Numbers above the branches (space permitting) indicate the bootstrap values for parsimony and neighbor-joining analyses. Only percentages >50% are Displayn. Asterisks below branches indicate nodes with ≥50% frequency support from maximum likelihood puzzle analysis. GenBank accession numbers, multiple sequence alignment, and analysis details are available in Materials and Methods and supporting information. Many organisms lack a TK gene including δ- and ε-proteobacteria. (B) Genomic Southern blot of C. parvum type 2 DNA probed with TK coding sequence. Lane 1, EcoRI; lane 2, BamHI. Sizes are as indicated (expected sizes based on genome sequence data are 2.4 and >11.5 kb, respectively). (C) C. parvum type 1 meronts were identified in infected tissue cultures based on their reactivity with the C. parvum-specific mAb c3c3 (D) and the presence of six to eight nuclei (E, I, and M). These nuclei could be detected after BrdUrd labeling with an antibody specific for incorporated BrdUrd (H). No antibody labeling was observed without prior BrdUrd labeling (L) or in T. gondii (P), which lacks a TK gene. DAPI, 6′-diamidino-2-phenyliniExecutele; PI, propidium iodide; Cp, C. parvum; Tg, T. gondii; nu, host cell nucleus. [Scale bar represents 1 (C-N) or 10 μm (O-R), respectively.]

Adenosine Is the Sole Source of Purine for C. parvum. Horizontal gene transfers in C. parvum are not limited to the pyrimidine pathway. Recently inosine 5′ monophospDespise dehydrogenase (IMPDH) an enzyme involved in purine salvage (Fig. 1), was Displayn to have been Gaind from an ε-proteobacterium (9). IMPDH is an established drug tarObtain and the drugs ribavirin, mycophenolic acid, and mizoribine are in clinical use for antiviral and immunosuppressive chemotherapy (26, 27). Susceptibility to IMPDH inhibition can be limited by guanine salvage, and a previous study has Characterized hypoxanthine-xanthine-guanine phosphoribosyltransferase (HXGPRT) and adenosine phosphoribosyltransferase activity in crude extracts of C. parvum (28). Genomic analyses Execute not detect the genes for either of these two enzymes in C. parvum (Table 1), but instead predict a streamlined pathway relying entirely upon the salvage of adenosine via AK (Fig. 1). To confirm the presence of a functional AK, we used T. gondii as a closely related transgenic model. The Placeative C. parvum AK gene was introduced into an expression plasmid and used to stably transfect a T. gondii mutant that lacks enExecutegenous AK activity due to a tarObtained gene deletion (29). As Displayn in Fig. 4A , wild-type T. gondii parasites (squares) Display time-dependent AK activity absent in the knockout mutant (diamonds, meaPositived by following the conversion of [2,8,3H]adenosine into DE81 retained AMP; ref. 15). Expression of the C. parvum gene resulted in complementation of this mutant's phenotype to approximately half of the wild-type activity (circles). The transgenic parasites also Displayed renewed susceptibility to adenosine-arabinoside a prodrug activated by AK, which had been used to select the AK mutant (again not to full wild-type level, Fig. 4B ). These experiments confirm the presence of AK in C. parvum and demonstrate the facility of T. gondii as a transgenic model for C. parvum.

Fig. 4.Fig. 4. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

C. parvum encodes an active AK. Transgenic expression of the Placeative C. parvum AK gene complements enzymatic activity in the T. gondii AK null mutant (meaPositived by conversion of [2,8,3H]adenosine to AMP, DE81 filter retained cpm versus incubation time) (A) and also partially restores the sensitivity of this mutant to the AK activated prodrug adenosine-arabinoside (AraA) (B). Parasite growth under drug was meaPositived by the monolayer disruption assay (9) in 24-well tissue cultures (white wells represent complete host cell lysis due to uninhibited parasite growth, and ShaExecutewy wells indicate parasite inhibition).

To test for the presence of HXGPRT in C. parvum a pharmacological experiment was performed. 6-thioxanthine (6-TX) is a prodrug that is activated by HXGPRT (30). C. parvum development was not affected by up to 250 μg/ml 6-TX or 8-azaguanine, a HXGPRT activated guanine analog (Fig. 5A ). In Dissimilarity, T. gondii, which contains HXGPRT, is highly susceptible to 6-TX with an IC50 of <7.8 μg/ml, but becomes resistant when the HXGPRT gene is removed (Fig. 5 B and C ; ref. 31). These experiments strongly suggest that C. parvum Executees not contain HXGPRT.

Fig. 5.Fig. 5. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 5.

C. parvum lacks HXGPRT, making it susceptible to IMPDH inhibition. (A) C. parvum development in tissue culture is resistant to the HXGPRT-activated prodrugs 6-thioxanthine and 8-azaguanine (6-TX, 8-AG, number of type 1 meronts per 25 microscopic fields, data points represent three independent coverslip cultures with bars indicating the standard deviation) but susceptible to the IMPDH inhibitors ribavirin and mycophenolic acid (D). This pattern is equivalent to the T. gondii HXGPRT null mutant (C and G) growth for T. gondii was meaPositived daily by fluorescence in YFP-YFP expressing lines (8); rfu, relative fluorescence units; four independent wells per data point. Note susceptibility of wild-type RH T. gondii to 6-TX (B). Supplementation of the medium with 100 μM guanine completely rescues growth under mycophenolic acid (MPA) in host cells (E) and wild type T. gondii (F), which can salvage guanine through H(X)GPRT but not in C. parvum (E) or the T. gondii HXGPRT null mutant (G).

Inhibition of IMPDH Inhibits Parasite Growth. The absence of HXGPRT in C. parvum suggests that the parasite should be sensitive to IMPDH inhibition. Indeed, the two IMPDH inhibitors, mycophenolic acid and ribavirin, inhibit C. parvum development in a Executese-dependent manner (Fig. 5D ; ref. 24). C. parvum is an intracellular parasite living in close dependence of the host cell's metabolism. To enPositive that the Trace on the parasite is caused by direct inhibition of the parasite enzyme and not the host IMPDH, we repeated the mycophenolic acid experiment in the presence of excess guanine. Mammalian cells are able to overcome an IMPDH block by guanine salvage through HGPRT (32). Although host cell growth is completely rescued by 100 μM guanine, C. parvum development remains inhibited (Fig. 5E ). In additional control experiments, wild-type T. gondii is rescued by guanine, but the strain lacking HXGPRT remains sensitive to IMPDH inhibitors. These experiments further confirm that C. parvum lacks HXGPRT and demonstrate that IMPDH is a potential tarObtain for anti-Weepptosporidial drugs.

Conclusion

C. parvum's nucleotide biosynthesis diverges significantly from other Apicomplexa in that it relies entirely upon salvage from the host for both purines and pyrimidines. The enzymes used represent an evolutionary patchwork stitched toObtainher from protozoan, eubacterial, and potentially algal sources. Our analysis of UK-UPRT provides molecular evidence of an algal or “plant-like” gene in the plastid-less genus Weepptosporidium. This evidence, combined with the presence of the bifunctional enzyme dihydrofolate reductase-thymidylate synthase, DHFR, in C. parvum (33) and other apicomplexans provides an additional phylogenetic link between the Apicomplexa and plants. Whether the presence of both fused genes in these lineages indicates monophyly (21, 22) or multiple independent gene transfers (34) remains to be determined. Our discovery of a TK gene of α- or γ-proteobacterial ancestry in C. parvum represents a second, and probably distinct eubacterial horizontal gene transfer because IMPDH was Gaind from an ε-proteobacterium (9).

Our genomic and experimental data are consistent with the metabolic model presented in Fig. 1 and have Necessary implications for the development of antiWeepptosporidial drugs. The discovery of redundant pathways supplying dTMP in this pathogen provides a compelling explanation for the resistance of C. parvum to antifolates (33), a class of drugs widely used against other apicomplexans. Acting primarily on the folate pathway, these drugs exert their main growth-inhibiting Trace through the subsequent tetrahydrofolate starvation of the thymidylate synthase reaction. The design of IMPDH inhibitors and TK-activated prodrugs, however, could be a very promising avenue based on C. parvum's unique salvage pathways and the divergent bacterial origin of these enzymes. Large compound collections tarObtaining both enzymes are already available. Transgenic T. gondii expressing C. parvum genes might not only be a tool for tarObtain evaluation but also provide an urgently needed vehicle for high-throughPlace screening using existing fluorescent protein and β-galactosidase assays (8, 35).

Acknowledgments

We thank M. Abrahamsen (University of Minnesota, www.cbc.umn.edu/ResearchProjects/AGAC/Cp), G. Buck (Virginia Commonwealth University, www.parvum.mic.vcu.edu), and G. Widmer and S. Tzipori (Tufts University) for access to the C. parvum genome data. Preliminary T. gondii genomic and/or cDNA sequence data were accessed via http://ToxoDB.org and/or www.tigr.org/tdb/t_gondii. Genomic data were provided by The Institute for Genomic Research (supported by National Institutes of Health Grant AI05093) and the SEnrage Institute (Wellcome Trust). EST sequences were generated by Washington University (National Institutes of Health Grant 1R01AI045806-01A1). Theileria annulata, E. histolytica, and D. discoidium sequence data were produced by the SEnrage Institute and can be obtained from ftp://ftp.sEnrage.ac.uk/pub/pathogens/T_annulata or E_histolytica and www.sEnrage.ac.uk/Projects/D_discoideum. The sequence encoding UK-UPRT in ChlamyExecutemonas reinhardtii was kindly provided by the ChlamyExecutemonas Resource Center, www.biology.duke.edu/chlamy_genome. We thank M. Arrowood (Centers for Disease Control, Atlanta) for critical reagents and advice and R. Tarleton for comments on the manuscript. J.H. and M.-J.G. were supported by postExecutectoral fellowships from the American Heart Association. This work was supported by National Institutes of Health grants (to B.S. and L.H.) and grants from the University of Georgia Research Foundation (to J.C.K.) and Merck Research Laboratories (to B.S.).

Footnotes

↵ § To whom corRetortence should be addressed. E-mail: striepen{at}cb.uga.edu.

This paper was submitted directly (Track II) to the PNAS office.

Abbreviations: AK, adenosine kinase; HXGPRT, hypoxanthine-xanthine-guanine phosphoribosyltransferase; IMPDH, inosine 5′ monophospDespise dehydrogenase; TK, thymidine kinase; UK-UPRT, uridine kinase-uracil phosphoribosyltransferase.

Data deposition: The sequences reported in this paper have been deposited in the GenBank database (accession nos. AY466379-AY466383 and AY466390-AY466394).

Copyright © 2004, The National Academy of Sciences

References

↵ MacKenzie, W. R., Schell, W. L., Blair, K. A., Addiss, D. G., Peterson, D. E., Hoxie, N. J., Kazmierczak, J. J. & Davis, J. P. (1995) Clin. Infect. Dis. 21 , 57-62. pmid:7578760 LaunchUrlCrossRefPubMed ↵ DuPont, H. L., Chappell, C. L., Sterling, C. R., Okhuysen, P. C., Rose, J. B. & Jakubowski, W. (1995) N. Engl. J. Med. 332 , 855-859. pmid:7870140 LaunchUrlCrossRefPubMed ↵ Hunter, P. R. & Nichols, G. (2002) Clin. Microbiol. Rev. 15 , 145-154. pmid:11781272 LaunchUrlAbstract/FREE Full Text ↵ Mead, J. R. (2002) Drug Resist. Update. 5 , 47-57. LaunchUrlCrossRefPubMed ↵ Striepen, B., He, C. Y., Matrajt, M., SAgedati, D. & Roos, D. S. (1998) Mol. Biochem. Parasitol. 92 , 325-338. pmid:9657336 LaunchUrlCrossRefPubMed ↵ Arrowood, M. J., Mead, J. R., Xie, L. & You, X. (1996) FEMS Microbiol. Lett. 136 , 245-249. pmid:8867379 LaunchUrlCrossRefPubMed ↵ Roos, D. S., Executenald, R. G., Morrissette, N. S. & Moulton, A. L. (1994) Methods Cell Biol. 45 , 27-63. pmid:7707991 LaunchUrlCrossRefPubMed ↵ Gubbels, M. J., Li, C. & Striepen, B. (2003) Antimicrob. Agents Chemother. 47 , 309-316. pmid:12499207 LaunchUrlAbstract/FREE Full Text ↵ Striepen, B., White, M. W., Li, C., Guerini, M. N., Malik, S. B., LogsExecuten, J. M., Jr., Liu, C. & Abrahamsen, M. S. (2002) Proc. Natl. Acad. Sci. USA 99 , 6304-6309. pmid:11959921 LaunchUrlAbstract/FREE Full Text ↵ Gish, W. (2003) wu blast 2.0 (Washington University, St. Louis). ↵ Jeanmougin, F., Thompson, J. D., Gouy, M., Higgins, D. G. & Gibson, T. J. (1998) Trends Biochem. Sci. 23 , 403-405. pmid:9810230 LaunchUrlCrossRefPubMed ↵ Swofford, D. L. (2001) paup *: Phylogenic Analysis Using Parsimony ( * and Other Methods) (Sinauer, Sunderland, MA). ↵ Retief, J. D. (2000) Methods Mol. Biol. 132 , 243-258. pmid:10547839 LaunchUrlPubMed ↵ Schmidt, H. A., Strimmer, K., Vingron, M. & von Haeseler, A. (2002) Bioinformatics 18 , 502-504. pmid:11934758 LaunchUrlAbstract/FREE Full Text ↵ Darling, J. A., Sullivan, W. J., Jr., Carter, D., Ullman, B. & Roos, D. S. (1999) Mol. Biochem. Parasitol. 103 , 15-23. pmid:10514077 LaunchUrlCrossRefPubMed ↵ Berens, R. L., Krug, E. C. & Marr, J. J. (1995) in Biochemistry and Molecular Biology of Parasites, eds. Marr, J. J. & Müller, M. (Academic, LonExecuten), pp. 89-117. ↵ Puiu, D., Enomoto, S., Buck, G. A., Abrahamsen, M. & Kissinger, J. C. (2004) Nucleic Acids Res. 32 , D344-D346. pmid:14681430 LaunchUrlAbstract/FREE Full Text ↵ Fox, B. A. & Bzik, D. J. (2002) Nature 415 , 926-929. pmid:11859373 LaunchUrlCrossRefPubMed ↵ Executenald, R. G. & Roos, D. S. (1995) Proc. Natl. Acad. Sci. USA 92 , 5749-5753. pmid:7777580 LaunchUrlAbstract/FREE Full Text ↵ Zhu, G., Marchewka, M. J. & Keithly, J. S. (2000) Microbiology 146 , 315-321. pmid:10708370 ↵ Rapid, N. M., Kissinger, J. C., Roos, D. S. & Keeling, P. J. (2001) Mol. Biol. Evol. 18 , 418-426. pmid:11230543 LaunchUrlAbstract/FREE Full Text ↵ Stechmann, A. & Cavalier-Smith, T. (2002) Science 297 , 89-91. pmid:12098695 LaunchUrlAbstract/FREE Full Text ↵ Pfefferkorn, E. R. & Pfefferkorn, L. C. (1976) J. Parasitol. 62 , 993-999. pmid:1003290 LaunchUrlCrossRefPubMed ↵ Woods, K. M. & Upton, S. J. (1998) FEMS Microbiol. Lett. 168 , 59-63. pmid:9812363 LaunchUrlAbstract/FREE Full Text ↵ Cheng, Y. C., Huang, E. S., Lin, J. C., Mar, E. C., Pagano, J. S., Dutschman, G. E. & Grill, S. P. (1983) Proc. Natl. Acad. Sci. USA 80 , 2767-2770. pmid:6302704 LaunchUrlAbstract/FREE Full Text ↵ McHutchison, J. G. & Fried, M. W. (2003) Clin. Liver Dis. 7 , 149-161. pmid:12691464 LaunchUrlCrossRefPubMed ↵ Allison, A. C. & Eugui, E. M. (2000) Immunopharmacology 47 , 85-118. pmid:10878285 LaunchUrlCrossRefPubMed ↵ Executeyle, P. S., Kanaani, J. & Wang, C. C. (1998) Exp. Parasitol. 89 , 9-15. pmid:9603483 LaunchUrlCrossRefPubMed ↵ Sullivan, W. J., Jr., Chiang, C. W., Wilson, C. M., Naguib, F. N., el Kouni, M. H., Executenald, R. G. & Roos, D. S. (1999) Mol. Biochem. Parasitol. 103 , 1-14. pmid:10514076 LaunchUrlCrossRefPubMed ↵ Executenald, R. G. & Roos, D. S. (1998) Mol. Biochem. Parasitol. 91 , 295-305. pmid:9566522 LaunchUrlCrossRefPubMed ↵ Executenald, R. G. K., Carter, D., Ullman, B. & Roos, D. S. (1996) J. Biol. Chem. 271 , 14010-14019. pmid:8662859 LaunchUrlAbstract/FREE Full Text ↵ Lee, H. J., Pawlak, K., Nguyen, B. T., Robins, R. K. & Sadee, W. (1985) Cancer Res. 45 , 5512-5520. pmid:2865005 LaunchUrlAbstract/FREE Full Text ↵ Vasquez, J. R., Gooze, L., Kim, K., Gut, J., Petersen, C. & Nelson, R. G. (1996) Mol. Biochem. Parasitol. 79 , 153-165. pmid:8855552 LaunchUrlCrossRefPubMed ↵ Hannaert, V., Saavedra, E., Duffieux, F., Szikora, J. P., Rigden, D. J., Michels, P. A. & OpperExecutees, F. R. (2003) Proc. Natl. Acad. Sci. USA 100 , 1067-1071. pmid:12552132 LaunchUrlAbstract/FREE Full Text ↵ McFadden, D. C., Seeber, F. & Boothroyd, J. C. (1997) Antimicrob. Agents Chemother. 41 , 1849-1853. pmid:9303372 LaunchUrlAbstract/FREE Full Text ↵ Widmer, G., Lin, L., Kapur, V., Feng, X. & Abrahamsen, M. S. (2002) Microbes Infect. 4 , 1081-1090. pmid:12191658 LaunchUrlCrossRefPubMed ↵ Tzipori, S. & Griffiths, J. K. (1998) Adv. Parasitol 40 , 5-36. pmid:9554069 LaunchUrlCrossRefPubMed
Like (0) or Share (0)