PlasmoDB: a functional genomic database for malaria parasites

Functional diversification between two related Plasmodium falciparum merozoite invasion ligands is determined by changes in the cytoplasmic domain

The pathogenesis of Plasmodium falciparum depends on efficient invasion into host erythrocytes. Parasite ligands encoded by multi-gene families interact with erythrocyte receptors. P. falciparum reticulocyte binding protein homologues (PfRhs) are expressed at the apical surface of invasive merozoites and have divergent ectodomains that are postulated to bind different erythrocyte receptors. Variant expression of these paralogues results in the use of alternative invasion pathways. Two PfRh proteins, PfRh2a and PfRh2b, are identical for 2700 N-terminal amino acids and differ only in a C-terminal 500 amino acid region, which includes a unique ectodomain, transmembrane domain and cytoplasmic domain. Despite their similarity, PfRh2b is required for a well-defined invasion pathway while PfRh2a is not required or sufficient for this pathway. Mapping the genomic region encoding these proteins revealed a recombinogenic locus with PfRh2a and PfRh2b in a head-to-head orientation. We have generated viable PfRh2a/2b chimeric parasites to identify the regions required for alternative invasion pathway utilization. We find that the differential ability to use these pathways is conferred by the cytoplasmic domains of PfRh2a and PfRh2b, not the ectodomain or transmembrane regions. Our results highlight the importance of the cytoplasmic domain for functional diversification of a major adhesive ligand family in malaria parasites.

Glycerophospholipid acquisition in Plasmodium - a puzzling assembly of biosynthetic pathways

Throughout the Plasmodium life cycle, malaria parasites repeatedly undergo rapid cellular growth and prolific divisions, necessitating intense membrane neogenesis and, in particular, the acquisition of high amounts of phospholipids. At the intraerythrocytic stage, glycerophospholipids are the main parasite membrane constituents, which mostly originate from the Plasmodium-encoded enzymatic machinery. Several proteins and entire pathways have been characterized and their features reported, thereby generating a global view of glycerophospholipid synthesis across Plasmodium spp. The malaria parasite displays a panoply of pathways that are seldom found together in a single organism. The major glycerophospholipids are synthesized via ancestral prokaryotic CDP-diacylglycerol-dependent pathways and eukaryotic-type de novo pathways. The parasite exhibits additional reactions that bridge some of these routes and are otherwise restricted to some organisms, such as plants, while base-exchange mechanisms are largely unexplored in Plasmodium. Marked differences between Plasmodium spp. have also been reported in phosphatidylcholine and phosphatidylethanolamine synthesis. Little is currently known about glycerophospholipid acquisition at non-erythrocytic stages, but recent data reveal that intrahepatocytic parasites, oocysts and sporozoites import various host lipids, and that de novo fatty acid synthesis is only crucial at the late liver stage. More studies on the different Plasmodium developmental stages are needed, to further assemble the different pieces of this glycerophospholipid synthesis puzzle, which contains highly promising therapeutic targets.

Differentiating Plasmodium falciparum alleles by transforming Cartesian X,Y data to polar coordinates

BACKGROUND

Diagnosis of infectious diseases now benefits from advancing technology to perform multiplex analysis of a growing number of variables. These advances enable simultaneous surveillance of markers characterizing species and strain complexity, mutations associated with drug susceptibility, and antigen-based polymorphisms in relation to evaluation of vaccine effectiveness. We have recently developed assays detecting single nucleotide polymorphisms (SNPs) in the P. falciparum genome that take advantage of post-PCR ligation detection reaction and fluorescent microsphere labeling strategies. Data from these assays produce a spectrum of outcomes showing that infections result from single to multiple strains. Traditional methods for distinguishing true positive signal from background can cause false positive diagnoses leading to incorrect interpretation of outcomes associated with disease treatment.

RESULTS

Following analysis of Plasmodium falciparum dihydrofolate reductase SNPs associated with resistance to a commonly used antimalarial drug, Fansidar (Sulfadoxine/pyrimethamine), and presumably neutral SNPs for parasite strain differentiation, we first evaluated our data after setting a background signal based on the mean plus three standard deviations for known negative control samples. Our analysis of single allelic controls suggested that background for the absent allele increased as the concentration of the target allele increased. To address this problem, we introduced a simple change of variables from customary (X,Y) (Cartesian) coordinates to planar polar coordinates (X = rcos(theta), Y = rsin(theta)). Classification of multidimensional fluorescence signals based on histograms of angular and radial data distributions proved more effective than classification based on Cartesian thresholds. Comparison with known diallelic dilution controls suggests that histogram-based classification is effective for major:minor allele concentration ratios as high as 10:1.

CONCLUSION

We have observed that the diallelic SNP data resulting from analysis of P. falciparum mutations is more accurately diagnosed when a simple polar transform of the (X,Y) data into (r,theta) is used. The development of high through-put methods for genotyping P. falciparum SNPs and the refinement of analytical approaches for evaluating these molecular diagnostic results significantly advance the evaluation of parasite population diversity and antimalarial drug resistance.

Thousands of chemical starting points for antimalarial lead identification

Malaria is a devastating infection caused by protozoa of the genus Plasmodium. Drug resistance is widespread, no new chemical class of antimalarials has been introduced into clinical practice since 1996 and there is a recent rise of parasite strains with reduced sensitivity to the newest drugs. We screened nearly 2 million compounds in GlaxoSmithKline's chemical library for inhibitors of P. falciparum, of which 13,533 were confirmed to inhibit parasite growth by at least 80% at 2 microM concentration. More than 8,000 also showed potent activity against the multidrug resistant strain Dd2. Most (82%) compounds originate from internal company projects and are new to the malaria community. Analyses using historic assay data suggest several novel mechanisms of antimalarial action, such as inhibition of protein kinases and host-pathogen interaction related targets. Chemical structures and associated data are hereby made public to encourage additional drug lead identification efforts and further research into this disease.

Cerebral malaria: insights from host-parasite protein-protein interactions

Background

Cerebral malaria is a form of human malaria wherein Plasmodium falciparum-infected red blood cells adhere to the blood capillaries in the brain, potentially leading to coma and death. Interactions between parasite and host proteins are important in understanding the pathogenesis of this deadly form of malaria. It is, therefore, necessary to study available protein-protein interactions to identify lesser known interactions that could throw light on key events of cerebral malaria.

Methods

Sequestration, haemostasis dysfunction, systemic inflammation and neuronal damage are key processes of cerebral malaria. Key events were identified from literature as being crucial to these processes. An integrated interactome was created using available experimental and predicted datasets as well as from literature. Interactions from this interactome were filtered based on Gene Ontology and tissue-specific annotations, and further analysed for relevance to the key events.

Results

PfEMP1 presentation, platelet activation and astrocyte dysfunction were identified as the key events influencing the disease. 48896 host-parasite along with other host-parasite, host-host and parasite-parasite protein-protein interactions obtained from a disease-specific corpus were combined to form an integrated interactome. Filtering of the interactome resulted in five host-parasite PPI, six parasite-parasite and two host-host PPI. The analysis of these interactions revealed the potential significance of apolipoproteins and temperature/Hsp expression on efficient PfEMP1 presentation; role of MSP-1 in platelet activation; effect of parasite proteins in TGF-β regulation and the role of albumin in astrocyte dysfunction.

Conclusions

This work links key host-parasite, parasite-parasite and host-host protein-protein interactions to key processes of cerebral malaria and generates hypotheses for disease pathogenesis based on a filtered interaction dataset. These hypotheses provide novel and significant insights to cerebral malaria.

The Kennedy phospholipid biosynthesis pathways are refractory to genetic disruption in Plasmodium berghei and therefore appear essential in blood stages

Phosphatidylcholine (PC) and phosphatidylethanolamine (PE) are the main membrane phospholipids (PLs) of Plasmodium parasites and can be generated by the de novo (Kennedy) CDP-choline and CDP-ethanolamine pathways and by the CDP-diacylglycerol dependent pathway. The Kennedy pathways initiate from exogenous choline and ethanolamine involving choline kinase (CK) and ethanolamine kinase (EK), followed by the choline-phosphate cytidylyltransferase (CCT) and ethanolamine-phosphate cytidylyltransferase (ECT) that catalyse the formation of CDP-choline and CDP-ethanolamine. Finally, in Plasmodium, PC and PE are apparently synthesized by a common choline/ethanolamine-phosphotransferase (CEPT). Here, we have studied the essential nature of the Kennedy pathways in Plasmodium berghei, a rodent malaria parasite. Sequence analysis of the P. berghei CEPT, CCT, ECT and CK enzymes revealed the presence of all catalytic domains and essential residues and motifs necessary for enzymatic activities. Constructs were designed for the generation of gene knockout and GFP-fusions of the cept, cct, ect and ck genes in P. berghei. We found that all four genes were consistently refractory to knockout attempts. At the same time, successful tagging of these proteins with GFP demonstrated that the loci were targetable and indicated that these genes are essential in P. berghei blood stage parasites. GFP-fusions of CCT, ECT and CK were found in the cytosol whereas the GFP-CEPT mainly localised in the endoplasmic reticulum. These results indicate that both CDP-choline and CDP-ethanolamine de novo pathways are essential for asexual P. berghei development and are non-redundant with other possible sources of PC and PE.

Geographic structure of Plasmodium vivax: microsatellite analysis of parasite populations from Sri Lanka, Myanmar, and Ethiopia

Genetic diversity and population structure of Plasmodium vivax parasites can predict the origin and spread of novel variants within a population enabling population specific malaria control measures. We analyzed the genetic diversity and population structure of 425 P. vivax isolates from Sri Lanka, Myanmar, and Ethiopia using 12 trinucleotide and tetranucleotide microsatellite markers. All three parasite populations were highly polymorphic with 3-44 alleles per locus. Approximately 65% were multiple-clone infections. Mean genetic diversity (H(E)) was 0.7517 in Ethiopia, 0.8450 in Myanmar, and 0.8610 in Sri Lanka. Significant linkage disequilibrium was maintained. Population structure showed two clusters (Asian and African) according to geography and ancestry. Strong clustering of outbreak isolates from Sri Lanka and Ethiopia was observed. Predictive power of ancestry using two-thirds of the isolates as a model identified 78.2% of isolates accurately as being African or Asian. Microsatellite analysis is a useful tool for mapping short-term outbreaks of malaria and for predicting ancestry.

Rodent and nonrodent malaria parasites differ in their phospholipid metabolic pathways

Malaria, a disease affecting humans and other animals, is caused by a protist of the genus Plasmodium. At the intraerythrocytic stage, the parasite synthesizes a high amount of phospholipids through a bewildering number of pathways. In the human Plasmodium falciparum species, a plant-like pathway that relies on serine decarboxylase and phosphoethanolamine N-methyltransferase activities diverts host serine to provide additional phosphatidylcholine and phosphatidylethanolamine to the parasite. This feature of parasitic dependence toward its host was investigated in other Plasmodium species. In silico analyses led to the identification of phosphoethanolamine N-methyltransferase gene orthologs in primate and bird parasite genomes. However, the gene was not detected in the rodent P. berghei, P. yoelii, and P. chabaudi species. Biochemical experiments with labeled choline, ethanolamine, and serine showed marked differences in biosynthetic pathways when comparing rodent P. berghei and P. vinckei, and human P. falciparum species. Notably, in both rodent parasites, ethanolamine and serine were not significantly incorporated into phosphatidylcholine, indicating the absence of phosphoethanolamine N-methyltransferase activity. To our knowledge, this is the first study to highlight a crucial difference in phospholipid metabolism between Plasmodium species. The findings should facilitate efforts to develop more rational approaches to identify and evaluate new targets for antimalarial therapy.

A novel DEAD box helicase Has1p from Plasmodium falciparum: N-terminal is essential for activity

Helicases catalyze the opening of nucleic acid duplexes and are implicated in many nucleic acid metabolic cellular processes that require single stranded DNA or reorganization of RNA structure. Previously we have reported that Plasmodium falciparum genome contains a number of DEAD box helicases. In the present study we report the cloning, expression and characterization of one of the novel members of DEAD box family from P. falciparum. Our results indicate that it is a homologue of Has1p from yeast and it contains DNA and RNA unwinding, nucleic acid-dependent ATPase and RNA binding activities. This enzyme can utilize all the nucleosidetriphosphates (NTPs) and deoxy nucleosidetriphosphates (dNTPs) for its unwinding activity. Using a truncated derivative of this protein we further report that the N-terminal region of the protein is essentially required for its activity. These studies suggest that besides the conserved helicase domain the highly variable N-terminal region also contributes in the activity of the protein.

Transcriptome analysis of the Cryptocaryon irritans tomont stage identifies potential genes for the detection and control of cryptocaryonosis

Background

Cryptocaryon irritans is a parasitic ciliate that causes cryptocaryonosis (white spot disease) in marine fish. Diagnosis of cryptocaryonosis often depends on the appearance of white spots on the surface of the fish, which are usually visible only during later stages of the disease. Identifying suitable biomarkers of this parasite would aid the development of diagnostic tools and control strategies for C. irritans. The C. irritans genome is virtually unexplored; therefore, we generated and analyzed expressed sequence tags (ESTs) of the parasite to identify genes that encode for surface proteins, excretory/secretory proteins and repeat-containing proteins.

Results

ESTs were generated from a cDNA library of C. irritans tomonts isolated from infected Asian sea bass, Lates calcarifer. Clustering of the 5356 ESTs produced 2659 unique transcripts (UTs) containing 1989 singletons and 670 consensi. BLAST analysis showed that 74% of the UTs had significant similarity (E-value < 10^-5) to sequences that are currently available in the GenBank database, with more than 15% of the significant hits showing unknown function. Forty percent of the UTs had significant similarity to ciliates from the genera Tetrahymena and Paramecium. Comparative gene family analysis with related taxa showed that many protein families are conserved among the protozoans. Based on gene ontology annotation, functional groups were successfully assigned to 790 UTs. Genes encoding excretory/secretory proteins and membrane and membrane-associated proteins were identified because these proteins often function as antigens and are good antibody targets. A total of 481 UTs were classified as encoding membrane proteins, 54 were classified as encoding for membrane-bound proteins, and 155 were found to contain excretory/secretory protein-coding sequences. Amino acid repeat-containing proteins and GPI-anchored proteins were also identified as potential candidates for the development of diagnostic and control strategies for C. irritans.

Conclusions

We successfully discovered and examined a large portion of the previously unexplored C. irritans transcriptome and identified potential genes for the development and validation of diagnostic and control strategies for cryptocaryonosis.

Exploring metabolomic approaches to analyse phospholipid biosynthetic pathways in Plasmodium

SUMMARYPlasmodium falciparum, the agent responsible for malaria, is an obligate intracellular protozoan parasite. For proliferation, differentiation and survival, it relies on its own protein-encoding genes, as well as its host cells for nutrient sources. Nutrients and subsequent metabolites are required by the parasites to support their high rate of growth and replication, particularly in the intra-erythrocytic stages of the parasite that are responsible for the clinical symptoms of the disease. Advances in mass spectrometry have improved the analysis of endogenous metabolites and enabled a global approach to identify the parasite's metabolites by the so-called metabolomic analyses. This level of analysis complements the genomic, transcriptomic and proteomic data already available and should allow the identification of novel metabolites, original pathways and networks of regulatory interactions within the parasite, and between the parasite and its hosts. The field of metabolomics is just in its infancy in P. falciparum, hence in this review, we concentrate on the available methodologies and their potential applications for deciphering important biochemical processes of the parasite, such as the astonishingly diverse phospholipid biosynthesis pathways. Elucidating the regulation of the biosynthesis of these crucial metabolites could help design of future anti-malarial drugs.

PlasmoPredict: a gene function prediction website for Plasmodium falciparum

The genome sequence of the malaria parasite Plasmodium falciparum was published in 2002 and revealed that approximately 60% of its genes could not be assigned a function. Eight years later the majority of P. falciparum proteins are still of unknown function. We therefore present PlasmoPredict, an easy-to-use online gene function prediction tool that integrates a wide range of functional genomics data for P. falciparum to aid in the annotation of these genes.

Comparison of the cellular and biochemical properties of Plasmodium falciparum choline and ethanolamine kinases

The proliferation of the malaria-causing parasite Plasmodium falciparum within the erythrocyte is concomitant with massive phosphatidylcholine and phosphatidylethanolamine biosynthesis. Based on pharmacological and genetic data, de novo biosynthesis pathways of both phospholipids appear to be essential for parasite survival. The present study characterizes PfCK (P. falciparum choline kinase) and PfEK (P. falciparum ethanolamine kinase), which catalyse the first enzymatic steps of these essential metabolic pathways. Recombinant PfCK and PfEK were expressed as His6-tagged fusion proteins from overexpressing Escherichia coli strains, then purified to homogeneity and characterized. Using murine polyclonal antibodies against recombinant kinases, PfCK and PfEK were shown to be localized within the parasite cytoplasm. Protein expression levels increased during erythrocytic development. PfCK and PfEK appeared to be specific to their respective substrates and followed Michaelis-Menten kinetics. The Km value of PfCK for choline was 135.3+/-15.5 microM. PfCK was also able to phosphorylate ethanolamine with a very low affinity. PfEK was found to be an ethanolamine-specific kinase (Km=475.7+/-80.2 microM for ethanolamine). The quaternary ammonium compound hemicholinium-3 and an ethanolamine analogue, 2-amino-1-butanol, selectively inhibited PfCK or PfEK. In contrast, the bis-thiazolium compound T3, which was designed as a choline analogue and is currently in clinical trials for antimalarial treatment, affected PfCK and PfEK activities similarly. Inhibition exerted by T3 was competitive for both PfCK and PfEK and correlated with the impairment of cellular phosphatidylcholine biosynthesis. Comparative analyses of sequences and structures for both kinase types gave insights into their specific inhibition profiles and into the dual capacity of T3 to inhibit both PfCK and PfEK.

Pathway Tools version 13.0: integrated software for pathway/genome informatics and systems biology

Pathway Tools is a production-quality software environment for creating a type of model-organism database called a Pathway/Genome Database (PGDB). A PGDB such as EcoCyc integrates the evolving understanding of the genes, proteins, metabolic network and regulatory network of an organism. This article provides an overview of Pathway Tools capabilities. The software performs multiple computational inferences including prediction of metabolic pathways, prediction of metabolic pathway hole fillers and prediction of operons. It enables interactive editing of PGDBs by DB curators. It supports web publishing of PGDBs, and provides a large number of query and visualization tools. The software also supports comparative analyses of PGDBs, and provides several systems biology analyses of PGDBs including reachability analysis of metabolic networks, and interactive tracing of metabolites through a metabolic network. More than 800 PGDBs have been created using Pathway Tools by scientists around the world, many of which are curated DBs for important model organisms. Those PGDBs can be exchanged using a peer-to-peer DB sharing system called the PGDB Registry.

EuPathDB: a portal to eukaryotic pathogen databases

EuPathDB (http://EuPathDB.org; formerly ApiDB) is an integrated database covering the eukaryotic pathogens of the genera Cryptosporidium, Giardia, Leishmania, Neospora, Plasmodium, Toxoplasma, Trichomonas and Trypanosoma. While each of these groups is supported by a taxon-specific database built upon the same infrastructure, the EuPathDB portal offers an entry point to all these resources, and the opportunity to leverage orthology for searches across genera. The most recent release of EuPathDB includes updates and changes affecting data content, infrastructure and the user interface, improving data access and enhancing the user experience. EuPathDB currently supports more than 80 searches and the recently-implemented ‘search strategy’ system enables users to construct complex multi-step searches via a graphical interface. Search results are dynamically displayed as the strategy is constructed or modified, and can be downloaded, saved, revised, or shared with other database users.

Structural determinants of DNA binding by a P. falciparum ApiAP2 transcriptional regulator

Putative transcription factors have only recently been identified in the Plasmodium spp., with the major family of regulators comprising the Apicomplexan Apetala2 (AP2) proteins. To better understand the DNA-binding mechanisms of these transcriptional regulators, we characterized the structure and in vitro function of an AP2 DNA-binding domain from a prototypical Apicomplexan AP2 protein, PF14_0633 from Plasmodium falciparum. The X-ray crystal structure of the PF14_0633 AP2 domain bound to DNA reveals a beta-sheet fold that binds the DNA major groove through base-specific and backbone contacts; a prominent alpha-helix supports the beta-sheet structure. Substitution of predicted DNA-binding residues with alanine weakened or eliminated DNA binding in solution. In contrast to plant AP2 domains, the PF14_0633 AP2 domain dimerizes upon binding to DNA through a domain-swapping mechanism in which the alpha-helices of the AP2 domains pack against the beta-sheets of the dimer mates. DNA-induced dimerization of PF14_0633 may be important for tethering two distal DNA loci together in the nucleus and/or for inducing functional rearrangements of its domains to facilitate transcriptional regulation. Consistent with a multisite binding mode, at least two copies of the consensus sequence recognized by PF14_0633 are present upstream of a previously identified group of sporozoite-stage genes. Taken together, these findings illustrate how Plasmodium has adapted the AP2 DNA-binding domain for genome-wide transcriptional regulation.

InParanoid 7: new algorithms and tools for eukaryotic orthology analysis

The InParanoid project gathers proteomes of completely sequenced eukaryotic species plus Escherichia coli and calculates pairwise ortholog relationships among them. The new release 7.0 of the database has grown by an order of magnitude over the previous version and now includes 100 species and their collective 1.3 million proteins organized into 42.7 million pairwise ortholog groups. The InParanoid algorithm itself has been revised and is now both more specific and sensitive. Based on results from our recent benchmarking of low-complexity filters in homology assignment, a two-pass BLAST approach was developed that makes use of high-precision compositional score matrix adjustment, but avoids the alignment truncation that sometimes follows. We have also updated the InParanoid web site (http://InParanoid.sbc.su.se). Several features have been added, the response times have been improved and the site now sports a new, clearer look. As the number of ortholog databases has grown, it has become difficult to compare among these resources due to a lack of standardized source data and incompatible representations of ortholog relationships. To facilitate data exchange and comparisons among ortholog databases, we have developed and are making available two XML schemas: SeqXML for the input sequences and OrthoXML for the output ortholog clusters.

Ensembl Genomes: extending Ensembl across the taxonomic space

Ensembl Genomes (http://www.ensemblgenomes.org) is a new portal offering integrated access to genome-scale data from non-vertebrate species of scientific interest, developed using the Ensembl genome annotation and visualisation platform. Ensembl Genomes consists of five sub-portals (for bacteria, protists, fungi, plants and invertebrate metazoa) designed to complement the availability of vertebrate genomes in Ensembl. Many of the databases supporting the portal have been built in close collaboration with the scientific community, which we consider as essential for maintaining the accuracy and usefulness of the resource. A common set of user interfaces (which include a graphical genome browser, FTP, BLAST search, a query optimised data warehouse, programmatic access, and a Perl API) is provided for all domains. Data types incorporated include annotation of (protein and non-protein coding) genes, cross references to external resources, and high throughput experimental data (e.g. data from large scale studies of gene expression and polymorphism visualised in their genomic context). Additionally, extensive comparative analysis has been performed, both within defined clades and across the wider taxonomy, and sequence alignments and gene trees resulting from this can be accessed through the site.

A cyclic GMP signalling module that regulates gliding motility in a malaria parasite

The ookinete is a motile stage in the malaria life cycle which forms in the mosquito blood meal from the zygote. Ookinetes use an acto-myosin motor to glide towards and penetrate the midgut wall to establish infection in the vector. The regulation of gliding motility is poorly understood. Through genetic interaction studies we here describe a signalling module that identifies guanosine 3′, 5′-cyclic monophosphate (cGMP) as an important second messenger regulating ookinete differentiation and motility. In ookinetes lacking the cyclic nucleotide degrading phosphodiesterase δ (PDEδ), unregulated signalling through cGMP results in rounding up of the normally banana-shaped cells. This phenotype is suppressed in a double mutant additionally lacking guanylyl cyclase β (GCβ), showing that in ookinetes GCβ is an important source for cGMP, and that PDEδ is the relevant cGMP degrading enzyme. Inhibition of the cGMP-dependent protein kinase, PKG, blocks gliding, whereas enhanced signalling through cGMP restores normal gliding speed in a mutant lacking calcium dependent protein kinase 3, suggesting at least a partial overlap between calcium and cGMP dependent pathways. These data demonstrate an important function for signalling through cGMP, and most likely PKG, in dynamically regulating ookinete gliding during the transmission of malaria to the mosquito.

Malaria parasites are single celled organisms, which must alternate between vertebrate and mosquito hosts to survive and spread. In both hosts, certain parasite stages can glide through tissues and invade cells. Many components of the molecular motor that powers gliding and invasion are known and we have a good idea how these may interact to generate force. It is less well understood how the motor is assembled and how its component parts are regulated to switch it on and off. We have begun to address these questions in the ookinete, a parasite stage, which forms in the blood meal of a mosquito and relies on gliding to penetrate the gut wall. Using a malaria parasite of rodents, we have examined the effect of deleting candidate genes involved in controlling levels of the intracellular signalling molecule cyclic guanosine monophosphate (cGMP). We show that the right balance between cGMP production and degradation is important for ookinetes to glide, while also maintaining their typical cell shape. Overall levels of cGMP are not much affected in the mutants, though, and we therefore believe the messenger exerts its effect either locally within the cell or only while the parasite is gliding.

Suggestive evidence for Darwinian Selection against asparagine-linked glycans of Plasmodium falciparum and Toxoplasma gondii

We are interested in asparagine-linked glycans (N-glycans) of Plasmodium falciparum and Toxoplasma gondii, because their N-glycan structures have been controversial and because we hypothesize that there might be selection against N-glycans in nucleus-encoded proteins that must pass through the endoplasmic reticulum (ER) prior to threading into the apicoplast. In support of our hypothesis, we observed the following. First, in protists with apicoplasts, there is extensive secondary loss of Alg enzymes that make lipid-linked precursors to N-glycans. Theileria makes no N-glycans, and Plasmodium makes a severely truncated N-glycan precursor composed of one or two GlcNAc residues. Second, secreted proteins of Toxoplasma, which uses its own 10-sugar precursor (Glc(3)Man(5)GlcNAc(2)) and the host 14-sugar precursor (Glc(3)Man(9)GlcNAc(2)) to make N-glycans, have very few sites for N glycosylation, and there is additional selection against N-glycan sites in its apicoplast-targeted proteins. Third, while the GlcNAc-binding Griffonia simplicifolia lectin II labels ER, rhoptries, and surface of plasmodia, there is no apicoplast labeling. Similarly, the antiretroviral lectin cyanovirin-N, which binds to N-glycans of Toxoplasma, labels ER and rhoptries, but there is no apicoplast labeling. We conclude that possible selection against N-glycans in protists with apicoplasts occurs by eliminating N-glycans (Theileria), reducing their length (Plasmodium), or reducing the number of N-glycan sites (Toxoplasma). In addition, occupation of N-glycan sites is markedly reduced in apicoplast proteins versus some secretory proteins in both Plasmodium and Toxoplasma.

Apicortin, a unique protein, with a putative cytoskeletal role, shared only by apicomplexan parasites and the placozoan Trichoplax adhaerens

A new protein, termed apicortin, has been identified, which contains a DCX (doublecortin) and a partial p25-alpha domain. The DCX domains of the doublecortin superfamily are responsible for their microtubule binding and stabilizing properties. The p25-alpha domain occurs in TPPPs (Tubulin Polymerization Promoting Proteins) exhibiting Microtubule Associated Protein (MAP)-like functions. TPPP orthologs can be classified as long- and short-type ones. The latter ones do not contain the partial p25-alpha domain, which is the most conservative part of the long-type TPPPs and was evolutionary-preserved independently from the whole domain. Apicortin seems to occur only in the placozoan animal Trichoplax adhaerens and in each of the apicomplexan parasites, the genomes of which are available. The function of the novel protein can be predicted from the fact that it contains two different microtubule-binding domains. Apicomplexans, important pathogens of humans and animals, possess unique cytoskeletal elements, as subpellicular microtubules and apical polar rings, which are specific for this phylum, extremely stabilized by unknown MAPs. One of the candidates can be the apicortin, which is supported by a recent experimental proteomics analysis. The role of the new protein in T. adhaerens may be similar but needs clarification.

A novel family of Apicomplexan glideosome-associated proteins with an inner membrane-anchoring role

The phylum Apicomplexa are a group of obligate intracellular parasites responsible for a wide range of important diseases. Central to the lifecycle of these unicellular parasites is their ability to migrate through animal tissue and invade target host cells. Apicomplexan movement is generated by a unique system of gliding motility in which substrate adhesins and invasion-related proteins are pulled across the plasma membrane by an underlying actin-myosin motor. The myosins of this motor are inserted into a dual membrane layer called the inner membrane complex (IMC) that is sandwiched between the plasma membrane and an underlying cytoskeletal basket. Central to our understanding of gliding motility is the characterization of proteins residing within the IMC, but to date only a few proteins are known. We report here a novel family of six-pass transmembrane proteins, termed the GAPM family, which are highly conserved and specific to Apicomplexa. In Plasmodium falciparum and Toxoplasma gondii the GAPMs localize to the IMC where they form highly SDS-resistant oligomeric complexes. The GAPMs co-purify with the cytoskeletal alveolin proteins and also to some degree with the actin-myosin motor itself. Hence, these proteins are strong candidates for an IMC-anchoring role, either directly or indirectly tethering the motor to the cytoskeleton.

Multiple var2csa-type PfEMP1 genes located at different chromosomal loci occur in many Plasmodium falciparum isolates

Background

The var2csa gene encodes a Plasmodium falciparum adhesion receptor which binds chondroitin sulfate A (CSA). This var gene is more conserved than other PfEMP1/var genes and is found in all P. falciparum isolates. In isolates 3D7, FCR3/It4 and HB3, var2csa is transcribed from a sub-telomeric position on the left arm of chromosome 12, but it is not known if this location is conserved in all parasites. Genome sequencing indicates that the var2csa gene is duplicated in HB3, but whether this is true in natural populations is uncertain.

Methodology/Principal Findings

To assess global variation in the VAR2CSA protein, sequence variation in the DBL2X region of var2csa genes in 54 P.falciparum samples was analyzed. Chromosome mapping of var2csa loci was carried out and a quantitative PCR assay was developed to estimate the number of var2csa genes in P.falciparum isolates from the placenta of pregnant women and from the peripheral circulation of other malaria patients. Sequence analysis, gene mapping and copy number quantitation in P.falciparum isolates indicate that there are at least two loci and that both var2csa-like genes can be transcribed. All VAR2CSA DBL2X domains fall into one of two distinct phylogenetic groups possessing one or the other variant of a large (∼26 amino acid) dimorphic motif, but whether either motif variant is linked to a specific locus is not known.

Conclusions/Significance

Two or more related but distinct var2csa-type PfEMP1/var genes exist in many P. falciparum isolates. One gene is on chromosome 12 but additional var2csa-type genes are on different chromosomes in different isolates. Multiplicity of var2csa genes appears more common in infected placentae than in samples from non-pregnant donors indicating a possible advantage of this genotype in pregnancy associated malaria.

Putative DNA G-quadruplex formation within the promoters of Plasmodium falciparum var genes

Background

Guanine-rich nucleic acid sequences are capable of folding into an intramolecular four-stranded structure called a G-quadruplex. When found in gene promoter regions, G-quadruplexes can downregulate gene expression, possibly by blocking the transcriptional machinery. Here we have used a genome-wide bioinformatic approach to identify Putative G-Quadruplex Sequences (PQS) in the Plasmodium falciparum genome, along with biophysical techniques to examine the physiological stability of P. falciparum PQS in vitro.

Results

We identified 63 PQS in the non-telomeric regions of the P. falciparum clone 3D7. Interestingly, 16 of these PQS occurred in the upstream region of a subset of the P. falciparum var genes (group B var genes). The var gene family encodes PfEMP1, the parasite's major variant antigen and adhesin expressed at the surface of infected erythrocytes, that plays a key role in malaria pathogenesis and immune evasion. The ability of the PQS found in the upstream regions of group B var genes (UpsB-Q) to form stable G-quadruplex structures in vitro was confirmed using ¹H NMR, circular dichroism, UV spectroscopy, and thermal denaturation experiments. Moreover, the synthetic compound BOQ1 that shows a higher affinity for DNA forming quadruplex rather than duplex structures was found to bind with high affinity to the UpsB-Q.

Conclusion

This is the first demonstration of non-telomeric PQS in the genome of P. falciparum that form stable G-quadruplexes under physiological conditions in vitro. These results allow the generation of a novel hypothesis that the G-quadruplex sequences in the upstream regions of var genes have the potential to play a role in the transcriptional control of this major virulence-associated multi-gene family.

The malaria merozoite, forty years on

The invasive blood stage of malaria parasites, merozoites, are complex entities specialized for the capture and entry of red blood cells. Their potential for vaccination and other anti-malaria strategies have attracted much research attention over the last 40 years, and there is now a considerable body of data relating to their biology. In this article some of the major advances over this period and remaining challenges are reviewed.

Redefining the expressed prototype SICAvar gene involved in Plasmodium knowlesi antigenic variation

Background

The SICAvar gene family, expressed at the surface of infected erythrocytes, is critical for antigenic variation in Plasmodium knowlesi. When this family was discovered, a prototypic SICAvar gene was characterized and defined by a 10-exon structure. The predicted 205-kDa protein lacked a convincing signal peptide, but included a series of variable cysteine-rich modules, a transmembrane domain encoded by the penultimate exon, and a cytoplasmic domain encoded by the final highly conserved exon. The 205 SICAvar gene and its family with up to 108 possible family members, was identified prior to the sequencing of the P. knowlesi genome. However, in the published P. knowlesi database this gene remains disjointed in five fragments. This study addresses a number of structural and functional questions that are critical for understanding SICAvar gene expression.

Methods

Database mining, bioinformatics, and traditional genomic and post-genomic experimental methods including proteomic technologies are used here to confirm the genomic context and expressed structure of the prototype 205 SICAvar gene.

Results

This study reveals that the 205 SICAvar gene reported previously to have a 10-exon expressed gene structure has, in fact, 12 exons, with an unusually large and repeat-laden intron separating two newly defined upstream exons and the bona fide 5'UTR from the remainder of the gene sequence. The initial exon encodes a PEXEL motif, which may function to localize the SICA protein in the infected erythrocyte membrane. This newly defined start of the 205 SICAvar sequence is positioned on chromosome 5, over 340 kb upstream from the rest of the telomerically positioned SICAvar gene sequence in the published genome assembly. This study, however, verifies the continuity of these sequences, a 9.5 kb transcript, and provides evidence that the 205 SICAvar gene is located centrally on chromosome 5.

Conclusion

The prototype 205 SICAvar gene has been redefined to have a 12-exon structure. These data are important because they 1) address questions raised in the P. knowlesi genome database regarding SICAvar gene fragments, numbers and structures, 2) show that this prototype gene encodes a PEXEL motif, 3) emphasize the need for further refinement of the P. knowlesi genome data, and 4) retrospectively, provide evidence for recombination within centrally located SICAvar sequences.

Discovery: an interactive resource for the rational selection and comparison of putative drug target proteins in malaria

Background

Up to half a billion human clinical cases of malaria are reported each year, resulting in about 2.7 million deaths, most of which occur in sub-Saharan Africa. Due to the over-and misuse of anti-malarials, widespread resistance to all the known drugs is increasing at an alarming rate. Rational methods to select new drug target proteins and lead compounds are urgently needed. The Discovery system provides data mining functionality on extensive annotations of five malaria species together with the human and mosquito hosts, enabling the selection of new targets based on multiple protein and ligand properties.

Methods

A web-based system was developed where researchers are able to mine information on malaria proteins and predicted ligands, as well as perform comparisons to the human and mosquito host characteristics. Protein features used include: domains, motifs, EC numbers, GO terms, orthologs, protein-protein interactions, protein-ligand interactions and host-pathogen interactions among others. Searching by chemical structure is also available.

Results

An in silico system for the selection of putative drug targets and lead compounds is presented, together with an example study on the bifunctional DHFR-TS from Plasmodium falciparum.

Conclusion

The Discovery system allows for the identification of putative drug targets and lead compounds in Plasmodium species based on the filtering of protein and chemical properties.

Plant-like phosphofructokinase from Plasmodium falciparum belongs to a novel class of ATP-dependent enzymes

Malaria parasite-infected erythrocytes exhibit enhanced glucose utilisation and 6-phospho-1-fructokinase (PFK) is a key enzyme in glycolysis. Here we present the characterisation of PFK from the human malaria parasite Plasmodium falciparum. Of the two putative PFK genes on chromosome 9 (PfPFK9) and 11 (PfPFK11), only the PfPFK9 gene appeared to possess all the catalytic features appropriate for PFK activity. The deduced PfPFK proteins contain domains homologous to the plant-like pyrophosphate (PPi)-dependent PFK beta and alpha subunits, which are quite different from the human erythrocyte PFK protein. The PfPFK9 gene beta and alpha regions were cloned and expressed as His(6)- and GST-tagged proteins in Escherichia coli. Complementation of PFK-deficient E. coli and activity analysis of purified recombinant proteins confirmed that PfPFK9beta possessed catalytic activity. Monoclonal antibodies against the recombinant beta protein confirmed that the PfPFK9 protein has beta and alpha domains fused into a 200 kDa protein, as opposed to the independent subunits found in plants. Despite an overall structural similarity to plant PPi-PFKs, the recombinant protein and the parasite extract exhibited only ATP-dependent enzyme activity, and none with PPi. Unlike host PFK, the Plasmodium PFK was insensitive to fructose-2,6-bisphosphate (F-2,6-bP), phosphoenolpyruvate (PEP) and citrate. A comparison of the deduced PFK proteins from several protozoan PFK genome databases implicates a unique class of ATP-dependent PFK present amongst the apicomplexan protozoans.

An unusual ERAD-like complex is targeted to the apicoplast of Plasmodium falciparum

Many apicomplexan parasites, including Plasmodium falciparum, harbor a so-called apicoplast, a complex plastid of red algal origin which was gained by a secondary endosymbiotic event. The exact molecular mechanisms directing the transport of nuclear-encoded proteins to the apicoplast of P. falciparum are not well understood. Recently, in silico analyses revealed a second copy of proteins homologous to components of the endoplasmic reticulum (ER)-associated protein degradation (ERAD) system in organisms with secondary plastids, including the malaria parasite P. falciparum. These proteins are predicted to be endowed with an apicoplast targeting signal and are suggested to play a role in the transport of nuclear-encoded proteins to the apicoplast. Here, we have studied components of this ERAD-derived putative preprotein translocon complex in malaria parasites. Using transfection technology coupled with fluorescence imaging techniques we can demonstrate that the N terminus of several ERAD-derived components targets green fluorescent protein to the apicoplast. Furthermore, we confirm that full-length PfsDer1-1 and PfsUba1 (homologues of yeast ERAD components) localize to the apicoplast, where PfsDer1-1 tightly associates with membranes. Conversely, PfhDer1-1 (a host-specific copy of the Der1-1 protein) localizes to the ER. Our data suggest that ERAD components have been "rewired" to provide a conduit for protein transport to the apicoplast. Our results are discussed in relation to the nature of the apicoplast protein transport machinery.

Proteome analysis of Plasmodium falciparum extracellular secretory antigens at asexual blood stages reveals a cohort of proteins with possible roles in immune modulation and signaling

The highly co-evolved relationship of parasites and their hosts appears to include modulation of host immune signals, although the molecular mechanisms involved in the host-parasite interplay remain poorly understood. Characterization of these key genes and their cognate proteins related to the host-parasite interplay should lead to a better understanding of this intriguing biological phenomenon. The malaria agent Plasmodium falciparum is predicted to export a cohort of several hundred proteins to remodel the host erythrocyte. However, proteins actively exported by the asexual intracellular parasite beyond the host red blood cell membrane (before merozoite egress) have been poorly investigated so far. Here we used two complementary methodologies, two-dimensional gel electrophoresis/MS and LC-MS/MS, to examine the extracellular secreted antigens at asexual blood stages of P. falciparum. We identified 27 novel antigens exported by P. falciparum in the culture medium of which some showed clustering with highly polymorphic genes on chromosomes, suggesting that they may encode putative antigenic determinants of the parasite. Immunolocalization of four novel secreted proteins confirmed their export beyond the infected red blood cell membrane. Of these, preliminary functional characterization of two novel (Sel1 repeat-containing) parasite proteins, PfSEL1 and PfSEL2 revealed that they down-regulate expression of cell surface Notch signaling molecules in host cells. Also a novel protein kinase (PfEK) and a novel protein phosphatase (PfEP) were found to, respectively, phosphorylate/dephosphorylate parasite-specific proteins in the extracellular culture supernatant. Our study thus sheds new light on malaria parasite extracellular secreted antigens of which some may be essential for parasite development and could constitute promising new drug targets.

Buffer optimization of thermal melt assays of Plasmodium proteins for detection of small-molecule ligands

In the past decade, thermal melt/thermal shift assays have become a common tool for identifying ligands and other factors that stabilize specific proteins. Increased stability is indicated by an increase in the protein's melting temperature (Tm). In optimizing the assays for subsequent screening of compound libraries, it is important to minimize the variability of Tm measurements so as to maximize the assay's ability to detect potential ligands. The authors present an investigation of Tm variability in recombinant proteins from Plasmodium parasites. Ligands of Plasmodium proteins are particularly interesting as potential starting points for drugs for malaria, and new drugs are urgently needed. A single standard buffer (100 mM HEPES [pH 7.5], 150 mM NaCl) permitted estimation of Tm for 58 of 61 Plasmodium proteins tested. However, with several proteins, Tm could not be measured with a consistency suitable for high-throughput screening unless alternative protein-specific buffers were employed. The authors conclude that buffer optimization to minimize variability in Tm measurements increases the success of thermal melt screens involving proteins for which a standard buffer is suboptimal.

The High Throughput Sequence Annotation Service (HT-SAS) - the shortcut from sequence to true Medline words

BACKGROUND

Advances in high-throughput technologies available to modern biology have created an increasing flood of experimentally determined facts. Ordering, managing and describing these raw results is the first step which allows facts to become knowledge. Currently there are limited ways to automatically annotate such data, especially utilizing information deposited in published literature.

RESULTS

To aid researchers in describing results from high-throughput experiments we developed HT-SAS, a web service for automatic annotation of proteins using general English words. For each protein a poll of Medline abstracts connected to homologous proteins is gathered using the UniProt-Medline link. Overrepresented words are detected using binomial statistics approximation. We tested our automatic approach with a protein test set from SGD to determine the accuracy and usefulness of our approach. We also applied the automatic annotation service to improve annotations of proteins from Plasmodium bergei expressed exclusively during the blood stage.

CONCLUSION

Using HT-SAS we created new, or enriched already established annotations for over 20% of proteins from Plasmodium bergei expressed in the blood stage, deposited in PlasmoDB. Our tests show this approach to information extraction provides highly specific keywords, often also when the number of abstracts is limited. Our service should be useful for manual curators, as a complement to manually curated information sources and for researchers working with protein datasets, especially from poorly characterized organisms.

The phosphoproteome of bloodstream form Trypanosoma brucei, causative agent of African sleeping sickness

The protozoan parasite Trypanosoma brucei is the causative agent of human African sleeping sickness and related animal diseases, and it has over 170 predicted protein kinases. Protein phosphorylation is a key regulatory mechanism for cellular function that, thus far, has been studied in T.brucei principally through putative kinase mRNA knockdown and observation of the resulting phenotype. However, despite the relatively large kinome of this organism and the demonstrated essentiality of several T. brucei kinases, very few specific phosphorylation sites have been determined in this organism. Using a gel-free, phosphopeptide enrichment-based proteomics approach we performed the first large scale phosphorylation site analyses for T.brucei. Serine, threonine, and tyrosine phosphorylation sites were determined for a cytosolic protein fraction of the bloodstream form of the parasite, resulting in the identification of 491 phosphoproteins based on the identification of 852 unique phosphopeptides and 1204 phosphorylation sites. The phosphoproteins detected in this study are predicted from their genome annotations to participate in a wide variety of biological processes, including signal transduction, processing of DNA and RNA, protein synthesis, and degradation and to a minor extent in metabolic pathways. The analysis of phosphopeptides and phosphorylation sites was facilitated by in-house developed software, and this automated approach was validated by manual annotation of spectra of the kinase subset of proteins. Analysis of the cytosolic bloodstream form T. brucei kinome revealed the presence of 44 phosphorylated protein kinases in our data set that could be classified into the major eukaryotic protein kinase groups by applying a multilevel hidden Markov model library of the kinase catalytic domain. Identification of the kinase phosphorylation sites showed conserved phosphorylation sequence motifs in several kinase activation segments, supporting the view that phosphorylation-based signaling is a general and fundamental regulatory process that extends to this highly divergent lower eukaryote.

Optimizing comparative genomic hybridization probes for genotyping and SNP detection in Plasmodium falciparum

Microarray-based comparative genomic hybridizations (CGH) interrogate genomic DNA to identify structural differences such as amplifications and deletions that are easily detected as large signal aberrations. Subtle signal deviations caused by single nucleotide polymorphisms (SNPs) can also be detected but is challenged by a high AT content (81%) in P. falciparum. We compared genome-wide CGH signal to sequence polymorphisms between parasite strains 3D7, HB3, and Dd2 using NimbleGen microarrays. From 23,191 SNPs (excluding var/rif/stevor genes), our CGH probe set detected SNPs with >99.9% specificity but low (<10%) sensitivity. Probe length, melting temperature, GC content, SNP location in the probe, mutation type, and hairpin structures affected SNP sensitivity. Previously unrecognized variable number tandem repeats (VNTRs) also were detected by this method. These findings will guide the redesign of a probe set to optimize an openly available CGH microarray platform for high-resolution genotyping suitable for population genomics studies.

Identification of microRNA in the protist Trichomonas vaginalis

MicroRNAs (miRNAs) are a class of small noncoding RNAs that have important regulatory roles in multicellular organisms. However, miRNA has never been identified experimentally in protist. Direct cloning of 438 expressed miRNA tags by microRNA serial analysis of gene expression from the parasitic protist Trichomonas vaginalis identified nine candidate miRNAs. Bioinformatics analysis of the corresponding genomic region revealed that these miRNA candidates contain a classical stem-loop-stem structure of pre-microRNAs. Analysis of the 20 nt long mature tva-miR-001 showed that it is an intergenic miRNA located at the scaffold DS113596. Tva-miR-001 was differentially expressed in the trophozoite, pseudocyst and amoeboid stages. Based on the experimental results of the present study, we provided solid evidence that protist possesses a miRNA regulating network comparable with multicellular organisms for the first time.

Quantitative protein expression profiling reveals extensive post-transcriptional regulation and post-translational modifications in schizont-stage malaria parasites

BACKGROUND

Malaria is a one of the most important infectious diseases and is caused by parasitic protozoa of the genus Plasmodium. Previously, quantitative characterization of the P. falciparum transcriptome demonstrated that the strictly controlled progression of these parasites through their intra-erythrocytic developmental cycle is accompanied by a continuous cascade of gene expression. Although such analyses have proven immensely useful, the correlations between abundance of transcripts and their cognate proteins remain poorly characterized.

RESULTS

Here, we present a quantitative time-course analysis of relative protein abundance for schizont-stage parasites (34 to 46 hours after invasion) based on two-dimensional differential gel electrophoresis of protein samples labeled with fluorescent dyes. For this purpose we analyzed parasite samples taken at 4-hour intervals from a tightly synchronized culture and established more than 500 individual protein abundance profiles with high temporal resolution and quantitative reproducibility. Approximately half of all profiles exhibit a significant change in abundance and 12% display an expression peak during the observed 12-hour time interval. Intriguingly, identification of 54 protein spots by mass spectrometry revealed that 58% of the corresponding proteins--including actin-I, enolase, eukaryotic initiation factor (eIF)4A, eIF5A, and several heat shock proteins--are represented by more than one isoform, presumably caused by post-translational modifications, with the various isoforms of a given protein frequently showing different expression patterns. Furthermore, comparisons with transcriptome data generated from the same parasite samples reveal evidence of significant post-transcriptional gene expression regulation.

CONCLUSIONS

Together, our data indicate that both post-transcriptional and post-translational events are widespread and of presumably great biological significance during the intra-erythrocytic development of P. falciparum.

Nucleic Acids Research annual Database Issue and the NAR online Molecular Biology Database Collection in 2009

The current issue of Nucleic Acids Research includes descriptions of 179 databases, of which 95 are new. These databases (along with several molecular biology databases described in other journals) have been included in the Nucleic Acids Research online Molecular Biology Database Collection, bringing the total number of databases in the collection to 1170. In this introductory comment, we briefly describe some of these new databases and review the principles guiding the selection of databases for inclusion in the Nucleic Acids Research annual Database Issue and the Nucleic Acids Research online Molecular Biology Database Collection. The complete database list and summaries are available online at the Nucleic Acids Research web site (http://nar.oxfordjournals.org/).