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

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

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

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

In vivo transcriptional profiling of Plasmodium falciparum

BACKGROUND: Both host and pathogen factors contribute to disease outcome in Plasmodium falciparum infection. The feasibility of studying the P. falciparum in vivo transcriptome to understand parasite transcriptional response while it resides in the human host is presented. METHODS: A custom made oligonucleotide array with probes based on the P. falciparum 3D7 laboratory strain chromosome 2 sequence was used to detect in vivo P. falciparum transcripts. This study analyzed transcripts from total RNA derived from small blood samples of P. falciparum infected patients and compared the in vivo expression profile to the in vitro cultivated 3D7 strain transcriptome. RESULTS: The data demonstrated that in vivo transcription can be studied from a small blood sample, despite the abundance of human RNA. The in vivo transcriptome is similar to the 3D7 ring stage transcriptome, but there are significant differences in genes encoding a sexual stage antigen and surface proteins. CONCLUSIONS: Whole genome transcription analysis of P. falciparum can be carried out successfully and further studies in selected patient cohorts may provide insight into parasite in vivo biology and defense against host immunity

Combining stage specificity and metabolomic profiling to advance antimalarial drug discovery

We report detailed susceptibility profiling of asexual blood stages of the malaria parasite Plasmodium falciparum to clinical and experimental antimalarials, combined with metabolomic fingerprinting. Results revealed a variety of stage-specific and metabolic profiles that differentiated the modes of action of clinical antimalarials including chloroquine, piperaquine, lumefantrine, and mefloquine, and identified late trophozoite-specific peak activity and stage-specific biphasic dose-responses for the mitochondrial inhibitors DSM265 and atovaquone. We also identified experimental antimalarials hitting previously unexplored druggable pathways as reflected by their unique stage specificity and/or metabolic profiles. These included several ring-active compounds, ones affecting hemoglobin catabolism through distinct pathways, and mitochondrial inhibitors with lower propensities for resistance than either DSM265 or atovaquone. This approach, also applicable to other microbes that undergo multiple differentiation steps, provides an effective tool to prioritize compounds for further development within the context of combination therapies