Systems analysis of host-parasite interactions.

Parasitic diseases caused by protozoan pathogens lead to hundreds of thousands of deaths per year in addition to substantial suffering and socioeconomic decline for millions of people worldwide. The lack of effective vaccines coupled with the widespread emergence of drug-resistant parasites necessitates that the research community take an active role in understanding host-parasite infection biology in order to develop improved therapeutics. Recent advances in next-generation sequencing and the rapid development of publicly accessible genomic databases for many human pathogens have facilitated the application of systems biology to the study of host-parasite interactions. Over the past decade, these technologies have led to the discovery of many important biological processes governing parasitic disease. The integration and interpretation of high-throughput -omic data will undoubtedly generate extraordinary insight into host-parasite interaction networks essential to navigate the intricacies of these complex systems. As systems analysis continues to build the foundation for our understanding of host-parasite biology, this will provide the framework necessary to drive drug discovery research forward and accelerate the development of new antiparasitic therapies

CACTI: an in silico chemical analysis tool through the integration of chemogenomic data and clustering analysis.

It is well-accepted that knowledge of a small molecules target can accelerate optimization. Although chemogenomic databases are helpful resources for predicting or finding compound interaction partners, they tend to be limited and poorly annotated. Furthermore, unlike genes, compound identifiers are often not standardized, and many synonyms may exist, especially in the biological literature, making batch analysis of compounds difficult. Here, we constructed an open-source annotation and target hypothesis prediction tool that explores some of the largest chemical and biological databases, mining these for both common name, synonyms, and structurally similar molecules. We used this Chemical Analysis and Clustering for Target Identification (CACTI) tool to analyze the Pathogen Box collection, an open-source set of 400 drug-like compounds active against a variety of microbial pathogens. Our analysis resulted in 4,315 new synonyms, 35,963 pieces of new information and target prediction hints for 58 members.Scientific contributionsWith the employment of this tool, a comprehensive report with known evidence, close analogs and drug-target prediction can be obtained for large-scale chemical libraries that will facilitate their evaluation and future target validation and optimization efforts

The effects of pH on the torsional flexibility of DNA bound to a nucleosome core particle

The effects of pH on the torsional flexibility of DNA bound to a nucleosome core particle were investigated by studying the time-resolved fluorescence anisotropy decays of ethidium bromide intercalated into the DNA of the core particle. As the torsional flexibility of DNA is affected by the presence of an intercalating dye, the decays were studied at different ethidium bromide to core particle binding ratios. The anisotropy decays were collected using the method of time-resolved single-photon counting and were fit to a model developed by J. M. Schurr (Schurr, 1984) using a non-linear least squares fitting algorithm developed by the author for this purpose. It was shown that below a binding ratio of 0.1 there was no demonstrable change in the anisotropy as a function of binding ratio. Our results show, that the apparent torsional flexibility of DNA of to a nucleosome core particle is dependent on the number of base pairs of the DNA between points of attachment to the histone core. If this number is as high as 30 base pairs, then the torsional flexibility of DNA on a nucleosome core particle is as high or higher than DNA free in solution. Also, for reasonable values of N, the friction felt by the DNA on a core particle is much higher than that felt by free DNA. This indicates that the DNA on a core particle is highly constrained in its motions. The hydrogen ion concentration was shown to have a substantial effect on the fluorescent anisotropy decays, particularly in the early regions of the decay. These analyses indicated that the observed change could be attributed to either a loosening of the contacts between the DNA and the histone core, or a relaxing of the torsional flexibility of the DNA

Genome-wide nucleosome mapping of Plasmodium falciparum reveals histone-rich coding and histone-poor intergenic regions and chromatin remodeling of core and subtelomeric genes

<p>Abstract</p> <p>Background</p> <p>Epigenetic modifications of histones and regulation of chromatin structure have been implicated in regulation of virulence gene families in <it>P. falciparum</it>. To better understand chromatin-mediated gene regulation, we used a high-density oligonucleotide microarray to map the position and enrichment of nucleosomes across the entire genome of <it>P. falciparum </it>at three time points of the intra-erythrocytic developmental cycle (IDC) in vitro. We used an unmodified histone H4 antibody for chromatin immunoprecipitation of nucleosome-bound DNA.</p> <p>Results</p> <p>We observed generally low nucleosomal occupancy of intergenic regions and higher occupancy of protein coding regions. In contract to the overall small fluctuation of nucleosomal occupancy in most coding regions throughout the IDC, subtelomeric genes encoding surface proteins such as <it>var </it>and <it>rif</it>, as well as some core chromosomal genes such as transcription factors, showed large changes in chromatin structure. Telomeres harbored a region with the highest nucleosomal occupancy of the genome and also exhibited large changes with higher nucleosomal occupancy at schizont stages. While many of these subtelomeric genes were previously shown to be modified by H3K9 trimethylation, we also identified some housekeeping genes in core chromosome regions that showed extensive changes in chromatin structure but do not contain this modification. tRNA and basal transcription factor genes showed low nucleosomal occupancy at all times, suggesting of an open chromatin structure that might be permissive for constitutively high levels of expression. Generally, nucleosomal occupancy was not correlated with the steady-state mRNA levels. Several <it>var </it>genes were exceptions: the <it>var </it>gene with the highest expression level showed the lowest nucleosomal occupancy, and selection of parasites for <it>var2CSA </it>expression resulted in lower nucleosomal occupancy at the <it>var2CSA </it>locus. We identified nucleosome-free regions in intergenic regions that may serve as transcription start sites or transcription factor binding sites. Using the nucleosomal occupancy data as the baseline, we further mapped the genome-wide enrichment of H3K9 acetylation and detected general enrichment of this mark in intergenic regions.</p> <p>Conclusions</p> <p>These data on nucleosome enrichment changes add to our understanding of the influence of chromatin structure on the regulation of gene expression. Histones are generally enriched in coding regions, and relatively poor in intergenic regions. Histone enrichment patterns allow for identification of new putative gene-coding regions. Most genes do not show correlation between chromatin structure and steady-state mRNA levels, indicating the dominant roles of other regulatory mechanisms. We present a genome-wide nucleosomal occupancy map, which can be used as a reference for future experiments of histone modification mapping.</p

Leveraging two-way probe-level block design for identifying differential gene expression with high-density oligonucleotide arrays

BACKGROUND: To identify differentially expressed genes across experimental conditions in oligonucleotide microarray experiments, existing statistical methods commonly use a summary of probe-level expression data for each probe set and compare replicates of these values across conditions using a form of the t-test or rank sum test. Here we propose the use of a statistical method that takes advantage of the built-in redundancy architecture of high-density oligonucleotide arrays. RESULTS: We employ parametric and nonparametric variants of two-way analysis of variance (ANOVA) on probe-level data to account for probe-level variation, and use the false-discovery rate (FDR) to account for simultaneous testing on thousands of genes (multiple testing problem). Using publicly available data sets, we systematically compared the performance of parametric two-way ANOVA and the nonparametric Mack-Skillings test to the t-test and Wilcoxon rank-sum test for detecting differentially expressed genes at varying levels of fold change, concentration, and sample size. Using receiver operating characteristic (ROC) curve comparisons, we observed that two-way methods with FDR control on sample sizes with 2–3 replicates exhibits the same high sensitivity and specificity as a t-test with FDR control on sample sizes with 6–9 replicates in detecting at least two-fold change. CONCLUSIONS: Our results suggest that the two-way ANOVA methods using probe-level data are substantially more powerful tests for detecting differential gene expression than corresponding methods for probe-set level data

Infection of laboratory-colonized Anopheles darlingi mosquitoes by Plasmodium vivax.

Anopheles darlingi Root is the most important malaria vector in the Amazonia region of South America. However, continuous propagation of An. darlingi in the laboratory has been elusive, limiting entomological, genetic/genomic, and vector-pathogen interaction studies of this mosquito species. Here, we report the establishment of an An. darlingi colony derived from wild-caught mosquitoes obtained in the northeastern Peruvian Amazon region of Iquitos in the Loreto Department. We show that the numbers of eggs, larvae, pupae, and adults continue to rise at least to the F6 generation. Comparison of feeding Plasmodium vivax ex vivo of F4 and F5 to F1 generation mosquitoes showed the comparable presence of oocysts and sporozoites, with numbers that corresponded to blood-stage asexual parasitemia and gametocytemia, confirming P. vivax vectorial capacity in the colonized mosquitoes. These results provide new avenues for research on An. darlingi biology and study of An. darlingi-Plasmodium interactions

Esterase mutation is a mechanism of resistance to antimalarial compounds

Pepstatin is a potent peptidyl inhibitor of various malarial aspartic proteases, and also has parasiticidal activity. Activity of pepstatin against cultured Plasmodium falciparum is highly variable depending on the commercial source. Here we identify a minor contaminant (pepstatin butyl ester) as the active anti-parasitic principle. We synthesize a series of derivatives and characterize an analogue (pepstatin hexyl ester) with low nanomolar activity. By selecting resistant parasite mutants, we find that a parasite esterase, PfPARE (P. falciparum Prodrug Activation and Resistance Esterase) is required for activation of esterified pepstatin. Parasites with esterase mutations are resistant to pepstatin esters and to an open source antimalarial compound, MMV011438. Recombinant PfPARE hydrolyses pepstatin esters and de-esterifies MMV011438. We conclude that (1) pepstatin is a potent but poorly bioavailable antimalarial; (2) PfPARE is a functional esterase that is capable of activating prodrugs; (3) Mutations in PfPARE constitute a mechanism of antimalarial resistance

In silico discovery of transcription regulatory elements in Plasmodium falciparum

<p>Abstract</p> <p>Background</p> <p>With the sequence of the <it>Plasmodium falciparum </it>genome and several global mRNA and protein life cycle expression profiling projects now completed, elucidating the underlying networks of transcriptional control important for the progression of the parasite life cycle is highly pertinent to the development of new anti-malarials. To date, relatively little is known regarding the specific mechanisms the parasite employs to regulate gene expression at the mRNA level, with studies of the <it>P. falciparum </it>genome sequence having revealed few <it>cis</it>-regulatory elements and associated transcription factors. Although it is possible the parasite may evoke mechanisms of transcriptional control drastically different from those used by other eukaryotic organisms, the extreme AT-rich nature of <it>P. falciparum </it>intergenic regions (~90% AT) presents significant challenges to <it>in silico cis</it>-regulatory element discovery.</p> <p>Results</p> <p>We have developed an algorithm called Gene Enrichment Motif Searching (GEMS) that uses a hypergeometric-based scoring function and a position-weight matrix optimization routine to identify with high-confidence regulatory elements in the nucleotide-biased and repeat sequence-rich <it>P. falciparum </it>genome. When applied to promoter regions of genes contained within 21 co-expression gene clusters generated from <it>P. falciparum </it>life cycle microarray data using the semi-supervised clustering algorithm Ontology-based Pattern Identification, GEMS identified 34 putative <it>cis</it>-regulatory elements associated with a variety of parasite processes including sexual development, cell invasion, antigenic variation and protein biosynthesis. Among these candidates were novel motifs, as well as many of the elements for which biological experimental evidence already exists in the <it>Plasmodium </it>literature. To provide evidence for the biological relevance of a cell invasion-related element predicted by GEMS, reporter gene and electrophoretic mobility shift assays were conducted.</p> <p>Conclusion</p> <p>This GEMS analysis demonstrates that <it>in silico </it>regulatory element discovery can be successfully applied to challenging repeat-sequence-rich, base-biased genomes such as that of <it>P. falciparum</it>. The fact that regulatory elements were predicted from a diverse range of functional gene clusters supports the hypothesis that <it>cis</it>-regulatory elements play a role in the transcriptional control of many <it>P. falciparum </it>biological processes. The putative regulatory elements described represent promising candidates for future biological investigation into the underlying transcriptional control mechanisms of gene regulation in malaria parasites.</p

Identification of pathogen genomic variants through an integrated pipeline

Background: Whole-genome sequencing represents a powerful experimental tool for pathogen research. We present methods for the analysis of small eukaryotic genomes, including a streamlined system (called Platypus) for finding single nucleotide and copy number variants as well as recombination events. Results: We have validated our pipeline using four sets of Plasmodium falciparum drug resistant data containing 26 clones from 3D7 and Dd2 background strains, identifying an average of 11 single nucleotide variants per clone. We also identify 8 copy number variants with contributions to resistance, and report for the first time that all analyzed amplification events are in tandem. Conclusions: The Platypus pipeline provides malaria researchers with a powerful tool to analyze short read sequencing data. It provides an accurate way to detect SNVs using known software packages, and a novel methodology for detection of CNVs, though it does not currently support detection of small indels. We have validated that the pipeline detects known SNVs in a variety of samples while filtering out spurious data. We bundle the methods into a freely available package