Applied systems biology and malaria

One of the goals of systems-biology research is to discover networks and interactions by integrating diverse data sets. So far, systems-biology research has focused on model organisms, which are well characterized and therefore suited to testing new methods. Systems biology has great potential for use in the search for therapies for disease. Here, the potential of systems-biology approaches in the search for new drugs and vaccines to treat malaria is examined.

The Plasmodium falciparum sexual development transcriptome: a microarray analysis using ontology-based pattern identification

The sexual stages of malarial parasites are essential for the mosquito transmission of the disease and therefore are the focus of transmission-blocking drug and vaccine development. In order to better understand genes important to the sexual development process, the transcriptomes of high-purity stage I-V Plasmodium falciparum gametocytes were comprehensively profiled using a full-genome high-density oligonucleotide microarray. The interpretation of this transcriptional data was aided by applying a novel knowledge-based data-mining algorithm termed ontology-based pattern identification (OPI) using current information regarding known sexual stage genes as a guide. This analysis resulted in the identification of a sexual development cluster containing 246 genes, of which approximately 75% were hypothetical, exhibiting highly-correlated, gametocyte-specific expression patterns. Inspection of the upstream promoter regions of these 246 genes revealed putative cis-regulatory elements for sexual development transcriptional control mechanisms. Furthermore, OPI analysis was extended using current annotations provided by the Gene Ontology Consortium to identify 380 statistically significant clusters containing genes with expression patterns characteristic of various biological processes, cellular components, and molecular functions. Collectively, these results, available as part of a web-accessible OPI database (http://carrier.gnf.org/publications/Gametocyte), shed light on the components of molecular mechanisms underlying parasite sexual development and other areas of malarial parasite biology.

A comparison of match-only algorithms for the analysis of Plasmodium falciparum oligonucleotide arrays

This study is motivated by two data sets which employ a custom Plasmodium falciparum version of the Affymetrix GeneChip, containing only perfect match (PM) oligonucleotides. A PM-only chip cannot be analysed using the standard Affymetrix-supplied software. We compared the performance of three match-only algorithms on these data: the Match Only Integral Distribution (MOID) algorithm, Robust Multichip Analysis (RMA), and the Model Based Expression Index (MBEI). We validated the differential expression of several genes using quantitative reverse transcriptase-PCR. We also performed a comparison using two publicly available 'benchmarking' data sets: the Latin Square spike-in data set generated by Affymetrix, and the Gene Logic dilution series. Since we know what the true fold changes are in these special data sets, they are helpful for assessment of expression algorithms.

Molecular mechanism for switching of P. falciparum invasion pathways into human erythrocytes

The malaria parasite, Plasmodium falciparum, exploits multiple ligand-receptor interactions, called invasion pathways, to invade the host erythrocyte. Strains of P. falciparum vary in their dependency on sialated red cell receptors for invasion. We show that switching from sialic acid-dependent to -independent invasion is reversible and depends on parasite ligand use. Expression of P. falciparum reticulocyte-binding like homolog 4 (PfRh4) correlates with sialic acid-independent invasion, and PfRh4 is essential for switching invasion pathways. Differential activation of PfRh4 represents a previously unknown mechanism to switch invasion pathways and provides P. falciparum with exquisite adaptability in the face of erythrocyte receptor polymorphisms and host immune responses.

Elucidating genetic diversity with oligonucleotide arrays

DNA microarrays, initially designed to measure gene expression levels, also provide an ideal platform for determining genetic diversity. Oligonucleotide microarrays, predominantly high-density oligonucleotide arrays, have emerged as the principal platforms for performing genome-wide diversity analysis. They have wide-ranging potential applications including comparative genomics, polymorphism discovery and genotyping. The identification of inheritable genetic markers also permits the analysis of quantitative traits, population studies and linkage analysis. In this review, we will discuss the application of oligonucleotide arrays, in particular high-density oligonucleotide arrays for elucidating genetic diversity and highlight some of the directions that the field may take.

In vivo transcriptome of Plasmodium falciparum reveals overexpression of transcripts that encode surface proteins

Infections with the human parasite Plasmodium falciparum continue to present a great challenge to global health. Fundamental questions regarding the molecular basis of virulence and immune evasion in P. falciparum have been only partially answered. Because of the parasite's intracellular location and complex life cycle, standard genetic approaches to the study of the pathogenesis of malaria have been limited. The present study presents a novel approach to the identification of the biological processes involved in host-pathogen interactions, one that is based on the analysis of in vivo P. falciparum transcripts. We demonstrate that a sufficient quantity of P. falciparum RNA transcripts can be derived from a small blood sample from infected patients for whole-genome microarray analysis. Overall, excellent correlation was observed between the transcriptomes derived from in vivo samples and in vitro samples with ring-stage P. falciparum 3D7 reference strain. However, gene families that encode surface proteins are overexpressed in vivo. Moreover, this analysis has identified a new family of hypothetical genes that may encode surface variant antigens. Comparative studies of the transcriptomes derived from in vivo samples and in vitro 3D7 samples may identify important strategies used by the pathogen for survival in the human host and highlight, for vaccine development, new candidate antigens that were not previously identified through the use of in vitro cultures.

Global analysis of transcript and protein levels across the Plasmodium falciparum life cycle

To investigate the role of post-transcriptional controls in the regulation of protein expression for the malaria parasite, Plasmodium falciparum, we have compared mRNA transcript and protein abundance levels for seven different stages of the parasite life cycle. A moderately high positive relationship between mRNA and protein abundance was observed for these stages; the most common discrepancy was a delay between mRNA and protein accumulation. Potentially post-transcriptionally regulated genes are identified, and families of functionally related genes were observed to share similar patterns of mRNA and protein accumulation.

Using expression information to discover new drug and vaccine targets in the malaria parasitePlasmodium falciparum

The recent completion of the malaria parasite Plasmodium falciparum genome has opened the door for applying a variety of genomic-based systems biology approaches that complement existing gene-by-gene methods of investigation. Transcriptomic analyses of P. falciparum using DNA microarrays has allowed for the rapid elucidation of gene function, parasite drug response, and in vivo expression profiles, as well as general mechanisms guiding the parasite life cycle that are vital to disease pathogenesis. The results of these studies have identified promising novel gene targets for the development of new drug and vaccine therapies.

Coordinated functions of WSS1, PSY2 and TOF1 in the DNA damage response

The stabilization and processing of stalled replication forks is required to maintain genome integrity in all organisms. In an effort to identify novel proteins that might be involved in stabilizing stalled replication forks, Saccharomyces cerevisiae mutant wss1Delta was isolated from a high-throughput screening of approximately 5000 deletion strains for genes involved in the response to continuous, low-intensity UV irradiation. Disruption of WSS1 resulted in synergistic increases in UV sensitivity with null mutants of genes involved in recombination (RAD52) and cell cycle control (RAD9 and RAD24). WSS1 was also found to interact genetically with SGS1, TOP3, SRS2 and CTF4, which are involved in recombination, repair of replication forks and the establishment of sister chromatid cohesion. A yeast two-hybrid screen identified a potential physical interaction between Wss1 and both Psy2 and Tof1. Genetic interactions were also detected between PSY2 and TOF1, as well as between each gene and RAD52 and SRS2, and between WSS1 and TOF1. Tof1 is known to be involved in stabilizing stalled replication forks and our data suggest that Wss1 and Psy2 similarly function to stabilize or process stalled or collapsed replication forks.

In silico gene function prediction using ontology-based pattern identification

MOTIVATION

With the emergence of genome-wide expression profiling data sets, the guilt by association (GBA) principle has been a cornerstone for deriving gene functional interpretations in silico. Given the limited success of traditional methods for producing clusters of genes with great amounts of functional similarity, new data-mining algorithms are required to fully exploit the potential of high-throughput genomic approaches.

RESULTS

Ontology-based pattern identification (OPI) is a novel data-mining algorithm that systematically identifies expression patterns that best represent existing knowledge of gene function. Instead of relying on a universal threshold of expression similarity to define functionally related groups of genes, OPI finds the optimal analysis settings that yield gene expression patterns and gene lists that best predict gene function using the principle of GBA. We applied OPI to a publicly available gene expression data set on the life cycle of the malarial parasite Plasmodium falciparum and systematically annotated genes for 320 functional categories based on current Gene Ontology annotations. An ontology-based hierarchical tree of the 320 categories provided a systems-wide biological view of this important malarial parasite.

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.

A Plasmodium gene family encoding Maurer's cleft membrane proteins: structural properties and expression profiling

Upon invasion of the erythrocyte cell, the malaria parasite remodels its environment; in particular, it establishes a complex membrane network, which connects the parasitophorous vacuole to the host plasma membrane and is involved in protein transport and trafficking. We have identified a novel subtelomeric gene family in Plasmodium falciparum that encodes 11 transmembrane proteins localized to the Maurer's clefts. Using coimmunoprecipitation and shotgun proteomics, we were able to enrich specifically for these proteins and detect distinct peptides, allowing us to conclude that four to 10 products were present at a given time. Nearly all of the Pfmc-2tm genes are transcribed during the trophozoite stage; this narrow time frame of transcription overlaps with the specific stevor and rif genes that are differentially expressed during the erythrocyte cycle. The description of the structural properties of the proteins led us to manually reannotate published sequences, and to detect potentially homologous gene families in both P. falciparum and Plasmodium yoelii yoelii, where no orthologs were predicted uniquely based on sequence similarity. These basic proteins with two transmembrane domains belong to a larger superfamily, which includes STEVORs and RIFINs.

Discovery of gene function by expression profiling of the malaria parasite life cycle

The completion of the genome sequence for Plasmodium falciparum, the species responsible for most malaria human deaths, has the potential to reveal hundreds of new drug targets and proteins involved in pathogenesis. However, only approximately 35% of the genes code for proteins with an identifiable function. The absence of routine genetic tools for studying Plasmodium parasites suggests that this number is unlikely to change quickly if conventional serial methods are used to characterize encoded proteins. Here, we use a high-density oligonucleotide array to generate expression profiles of human and mosquito stages of the malaria parasite's life cycle. Genes with highly correlated levels and temporal patterns of expression were often involved in similar functions or cellular processes.

Genetic diversity in yeast assessed with whole-genome oligonucleotide arrays

The availability of a complete genome sequence allows the detailed study of intraspecies variability. Here we use high-density oligonucleotide arrays to discover 11,115 single-feature polymorphisms (SFPs) existing in one or more of 14 different yeast strains. We use these SFPs to define regions of genetic identity between common laboratory strains of yeast. We assess the genome-wide distribution of genetic variation on the basis of this yeast population. We find that genome variability is biased toward the ends of chromosomes and is more likely to be found in genes with roles in fermentation or in transport. This subtelomeric bias may arise through recombination between nonhomologous sequences because full-gene deletions are more common in these regions than in more central regions of the chromosome.

The Ume6 regulon coordinates metabolic and meiotic gene expression in yeast

The Ume6 transcription factor in yeast is known to both repress and activate expression of diverse genes during growth and meiotic development. To obtain a more complete profile of the functions regulated by this protein, microarray analysis was used to examine transcription in wild-type and ume6Delta diploids during vegetative growth in glucose and acetate. Two different genetic backgrounds (W303 and SK1) were examined to identify a core set of strain-independent Ume6-regulated genes. Among genes whose expression is controlled by Ume6 in both backgrounds, 82 contain homologies to the Ume6-binding site (URS1) and are expected to be directly regulated by Ume6. The vast majority of those whose functions are known participate in carbon/nitrogen metabolism and/or meiosis. Approximately half of the Ume6 direct targets are induced during meiosis, with most falling into the early meiotic expression class (cluster 4), and a smaller subset in the middle and later classes (clusters 5-7). Based on these data, we propose that Ume6 serves a unique role in diploid cells, coupling metabolic responses to nutritional cues with the initiation and progression of meiosis. Finally, expression patterns in the two genetic backgrounds suggest that SK1 is better adapted to respiration and W303 to fermentation, which may in part account for the more efficient and synchronous sporulation of SK1.

Excess polymorphisms in genes for membrane proteins in Plasmodium falciparum

The detection of single-nucleotide polymorphisms in pathogenic microorganisms has normally been carried out by trial and error. Here we show that DNA hybridization with high-density oligonucleotide arrays provides rapid and convenient detection of single-nucleotide polymorphisms in Plasmodium falciparum, despite its exceptionally high adenine-thymine (AT) content (82%). A disproportionate number of polymorphisms are found in genes encoding proteins associated with the cell membrane. These genes are targets for only 22% of the oligonucleotide probes but account for 69% of the polymorphisms. Genetic variation is also enriched in subtelomeric regions, which account for 22% of the chromosome but 76% of the polymorphisms.

Large-scale mutagenesis and functional genomics in yeast

One of the major goals for the post-genome era is determining of the function of proteins predicted in the genome sequence. In many organisms functional assignments have been the results of comparative sequencing, proteomics or expression profiling. In the yeast, Saccharomyces cerevisiae, however, the functional role of a gene can be tested directly by disrupting the gene and examining the phenotype of the mutant. Because precise targeted deletions can be easily constructed, it is also possible to systematically delete every gene in the genome. Here we describe recent progress in yeast genome-wide mutagenesis programs and the results produced from analyzing the mutants created by them.

Monitoring the chromosome 2 intraerythrocytic transcriptome of Plasmodium falciparum using oligonucleotide arrays

To test the feasibility of using short oligonucleotide probes to monitor transcript levels in Plasmodium falciparum, a microarray was manufactured containing 4,167 (25 base single-stranded) probes derived from the predicted coding region of P. falciparum chromosome 2. RNA samples from three asexual stages (rings, trophozoites, and schizonts) were labeled and hybridized to the arrays. These results were reproducible, and transcripts were detected for 69% of the 210 genes on chromosome 2. In addition, of the 145 expressed genes, 1/3 appeared to be differentially transcribed during the asexual cycle. Some regions of the chromosome appeared to be transcriptionally silent. Results were confirmed by Northern blot analysis and by quantitative reverse transcriptase-polymerase chain reaction. These data validate the use of relatively short 25-mers for monitoring the expression of a genome that is 82% AT rich.

Treasures and traps in genome-wide data sets: case examples from yeast

Since the publication of the Saccharomyces cerevisiae genome sequence, much effort has been dedicated to developing high-throughput techniques to generate comprehensive information about the function and dynamics of all genes in this yeast's genome. These techniques have generated data sets that typically contain large amounts of reliable and valuable biological information. Nevertheless, there are also uncertainties that are associated with such large-scale studies, which we discuss in this review. These uncertainties increase with the complexity of the organism under study. On the basis of the results from yeast, we should learn much from human and mouse genomic data sets. However, as with yeast data sets, they might also contain misleading results.

Previously uncharacterized genes in the UV- and MMS-induced DNA damage response in yeast

A competitive growth assay has been used to identify yeast genes involved in the repair of UV- or MMS-induced DNA damage. A collection of 2,827 yeast strains was analyzed in which each strain has a single ORF replaced with a cassette containing two unique sequence tags, allowing for its detection by hybridization to a high-density oligonucleotide array. The hybridization data identify a high percentage of the deletion strains present in the collection that were previously characterized as being sensitive to the DNA-damaging agents. The assay, and subsequent analysis, has been used to identify six genes not formerly known to be involved in the damage response, whose deletion renders the yeast sensitive to UV or MMS treatment. The recently identified genes include three uncharacterized ORFs, as well as genes that encode protein products implicated in ubiquitination, gene silencing, and transport across the mitochondrial membrane. Epistatsis analysis of four of the genes was performed to determine the DNA damage repair pathways in which the protein products function.

Parallel identification of new genes in Saccharomyces cerevisiae

Short open reading frames (ORFs) occur frequently in primary genome sequence. Distinguishing bona fide small genes from the tens of thousands of short ORFs is one of the most challenging aspects of genome annotation. Direct experimental evidence is often required. Here we use a combination of expression profiling and mass spectrometry to verify the independent transcription of 138 and the translation of 50 previously nonannotated genes in the Saccharomyces cerevisiae genome. Through combined evidence, we propose the addition of 62 new genes to the genome and provide experimental support for the inclusion of 10 previously identified genes.

Functional profiling of the Saccharomyces cerevisiae genome

Determining the effect of gene deletion is a fundamental approach to understanding gene function. Conventional genetic screens exhibit biases, and genes contributing to a phenotype are often missed. We systematically constructed a nearly complete collection of gene-deletion mutants (96% of annotated open reading frames, or ORFs) of the yeast Saccharomyces cerevisiae. DNA sequences dubbed 'molecular bar codes' uniquely identify each strain, enabling their growth to be analysed in parallel and the fitness contribution of each gene to be quantitatively assessed by hybridization to high-density oligonucleotide arrays. We show that previously known and new genes are necessary for optimal growth under six well-studied conditions: high salt, sorbitol, galactose, pH 8, minimal medium and nystatin treatment. Less than 7% of genes that exhibit a significant increase in messenger RNA expression are also required for optimal growth in four of the tested conditions. Our results validate the yeast gene-deletion collection as a valuable resource for functional genomics.

The paradoxical population genetics of Plasmodium falciparum

Among the leading causes of death in African children is cerebral malaria caused by the parasitic protozoan Plasmodium falciparum. Endemic forms of this disease are thought to have originated in central Africa 5000-10000 years ago, coincident with the innovation of slash-and-burn agriculture and the diversification of the Anopheles gambiae complex of mosquito vectors. Population genetic studies of P. falciparum have yielded conflicting results. Some evidence suggests that today's population includes multiple ancient lineages pre-dating human speciation. Other evidence suggests that today's population derives from only one, or a small number, of these ancient lineages. Resolution of this issue is important for the evaluation of the long-term efficacy of drug and immunological control strategies.

Application of high-density array-based signature-tagged mutagenesis to discover novel Yersinia virulence-associated genes

Yersinia pestis, the causative agent of plague, and the enteropathogen Yersinia pseudotuberculosis have nearly identical nucleotide similarity yet cause markedly different diseases. To investigate this conundrum and to study Yersinia pathogenicity, we developed a high-density oligonucleotide array-based modification of signature-tagged mutagenesis (STM). Y. pseudotuberculosis YPIII mutants constructed with the tagged transposons were evaluated in the murine yersiniosis infection model. The DNA tags were amplified using biotinylated primers and hybridized to high-density oligonucleotide arrays containing DNA complementary to the tags. Comparison of the hybridization signals from input pools and output pools identified a mutant whose relative abundance was significantly reduced in the output pool. Sequence data from 31 transposon insertion regions was compared to the complete Y. pestis CO92 genome sequence. The 26 genes present in both species were found to be almost identical, but five Y. pseudotuberculosis genes identified through STM did not have counterparts in the Y. pestis genome and may contribute to the different tropisms in these closely related pathogens. Potential virulence genes identified include those involved in lipopolysaccharide biosynthesis, adhesion, phospholipase activity, iron assimilation, and gene regulation. The phospholipase A (PldA) mutant exhibited reduced phospholipase activity compared to the wild-type strain and in vivo attenuation of the mutant was confirmed. The combination of optimized double tag sequences and high-density array hybridization technology offers improved performance, efficiency, and reliability over classical STM and permits quantitative analysis of data.

Replication dynamics of the yeast genome

Oligonucleotide microarrays were used to map the detailed topography of chromosome replication in the budding yeast Saccharomyces cerevisiae. The times of replication of thousands of sites across the genome were determined by hybridizing replicated and unreplicated DNAs, isolated at different times in S phase, to the microarrays. Origin activations take place continuously throughout S phase but with most firings near mid-S phase. Rates of replication fork movement vary greatly from region to region in the genome. The two ends of each of the 16 chromosomes are highly correlated in their times of replication. This microarray approach is readily applicable to other organisms, including humans.